# Insecticide-Driven Patterns of Genetic Variation in the Dengue Vector Aedes aegypti in Martinique Island

Christophe Paupy, Sébastien Marcombe, Margot Paris, Charline Bringuier, André Yebakima, Fabrice Chandre, Jean-Philippe David, Vincent Corbel, Laurence Despres

► To cite this version:

Christophe Paupy, Sébastien Marcombe, Margot Paris, Charline Bringuier, André Yebakima, et al. Insecticide-Driven Patterns of Genetic Variation in the Dengue Vector Aedes aegypti in Martinique Island. PLoS ONE, 2013, 8 (10), pp.e77857. 10.1371/journal.pone.0077857. hal-03002444

HAL Id: hal-03002444

https://hal.umontpellier.fr/hal-03002444

Submitted on 26 May 2021

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
# Insecticide-Driven Patterns of Genetic Variation in the Dengue Vector *Aedes aegypti* in Martinique Island

Sébastien Marcombe¹*, Margot Paris²*anc*, Christophe Paupy¹,³, Charline Bringuier², André Yebakima⁴, Fabrice Chandre¹, Jean-Philippe David², Vincent Corbel¹*ab, Laurence Despres²*

1 UMR MIVEGEC (IRD 224-CNRS 5290-UM1-UM2), Institut de recherche pour le développement (IRD), Montpellier, France, 2 Laboratoire d'Ecologie Alpine (LECA), UMR 5553 CNRS- Grenoble Alpes Université, Grenoble, France, 3 Centre International de Recherches Médicales de Franceville (CIRMF), Franceville, Gabon, 4 Centre de Démoustication, Conseil Général de la Martinique, Fort de France, Martinique, France

## Abstract

Effective vector control is currently challenged worldwide by the evolution of resistance to all classes of chemical insecticides in mosquitoes. In Martinique, populations of the dengue vector *Aedes aegypti* have been intensively treated with temephos and deltamethrin insecticides over the last fifty years, resulting in heterogeneous levels of resistance across the island. Resistance spreading depends on standing genetic variation, selection intensity and gene flow among populations. To determine gene flow intensity, we first investigated neutral patterns of genetic variability in sixteen populations representative of the many environments found in Martinique and experiencing various levels of insecticide pressure, using 6 microsatellites. Allelic richness was lower in populations resistant to deltamethrin, and consanguinity was higher in populations resistant to temephos, consistent with a negative effect of insecticide pressure on neutral genetic diversity. The global genetic differentiation was low, suggesting high gene flow among populations, but significant structure was found, with a pattern of isolation-by-distance at the global scale. Then, we investigated adaptive patterns of divergence in six out of the 16 populations using 319 single nucleotide polymorphisms (SNPs). Five SNP outliers displaying levels of genetic differentiation out of neutral expectations were detected, including the *kdr*-V1016I mutation in the voltage-gated sodium channel gene. Association tests revealed a total of seven SNPs associated with deltamethrin resistance. Six other SNPs were associated with temephos resistance, including two non-synonymous substitutions in an alkaline phosphatase and in a sulfotransferase respectively. Altogether, both neutral and adaptive patterns of genetic variation in mosquito populations appear to be largely driven by insecticide pressure in Martinique.

**Citation:** Marcombe S, Paris M, Paupy C, Bringuier C, Yebakima A, et al. (2013) Insecticide-Driven Patterns of Genetic Variation in the Dengue Vector *Aedes aegypti* in Martinique Island. PLoS ONE 8(10): e77857. doi:10.1371/journal.pone.0077857

**Editor:** John Vontas, University of Crete, Greece

**Received** June 25, 2013; **Accepted** September 4, 2013; **Published** October 18, 2013

**Copyright:** © 2013 Marcombe 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:** Funding came from the Conseil Général de la Martinique University J. Fourier. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

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

## Introduction

During their life-span, mosquitoes face multiple selective pressures, including climatic constraints linked to an obligate aquatic larval stage, together with intense insecticide pressure. Insecticide pressure has led to the spreading of resistance alleles in all target species, and the managing of insecticide resistance in vectors of pathogens represents a major challenge for humanity. Originating from Africa, *Aedes aegypti* (Diptera, Culicidae) has spread throughout the tropics in historical times [1]. This mosquito species is the primary vector for dengue fever, yellow fever, and chikungunya virus in urban

and sub-urban tropical areas across the world. In Martinique Island, a French overseas department located in the eastern Caribbean in the Lesser Antilles, the dengue virus (*Flaviviridae*) is transmitted by *Ae. aegypti* following an endemo-epidemic pattern and is now a major public health issue [2]. In the absence of vaccine or specific treatment against the disease, vector control by the use of chemical or biological agents against *Ae. aegypti* is currently the main method for reducing dengue transmission, in addition to environmental management and educational programs

Over the last fifty years, larval and adult *Ae. aegypti* populations from Martinique have been intensively targeted

<footnote>* E-mail: laurence.despres@ujf-grenoble.fr</footnote>

<footnote>‡ These authors contributed equally to this work.</footnote>

<footnote>a Current address: Plant Ecological Genetics, Institute of Integrative Biology, ETH Zurich, Zurich, Switzerland</footnote>

<footnote>b Current address: Department of Entomology, Faculty of Agriculture, Kasetsart University, Bangkok, Thailand</footnote>
with insecticide applications. Organochlorines such as DDT were first used since the 1950s and were replaced in the early 1990s by pyrethroid-based formulations (i.e. deltamethrin) for adult control by aerial spraying applications [3]. In addition to adult control, the organophosphate temephos (Abate®) was used for decades for larval control and there was a switch in 2009 to the bio-insecticide *Bacillus thuringiensis var. israelensis* (Bti). Because of the intense pesticide selection pressure in Martinique, vector control is now facing operational challenges with the emergence and development of insecticide resistance in *Ae. aegypti* populations [4] [5]. Resistance to organophosphates and pyrethroids has been reported since the 1980s and the 1990s respectively. A previous study has shown that deltamethrin and temephos resistances were widespread in Martinique. Molecular screening for common insecticide target-site mutations revealed the presence of the "knock-down resistance" V1016I *kdr* mutation at high frequency (>87%), and real time quantitative RT-qPCR showed the potential involvement of several candidate detoxification genes in insecticide resistance [3]. This study also showed that the insecticide resistance in *Ae. aegypti* populations was not homogeneously distributed across Martinique and that there were potential links existing between resistance distribution and vector control operations but also other factors such as climate and urbanization. Indeed, these factors are likely affecting the dynamics and the selection of insecticide resistance mechanisms in multiple ways [6]. The speed of resistance spreading depends on standing genetic variation in natural populations, selection intensity, local adaptation and gene flow among populations. The different climate and geological areas of Martinique with the northern part, mountainous, volcanic and rainy and the southern part slightly rugged, sunny and dry, together with the constant insecticide selection pressures and the well-developed road network that may favour passive gene flows, could lead to specific patterns of genetic diversity and population structure. The analyse of both neutral and adaptive genetic variation using genome-wide distributed molecular markers represents a powerful tool to better understand the genetic structure of these populations as well as the molecular mechanisms underlying insecticide resistance in Martinique Island.

First, using extensive field sampling (16 populations and 31 sub-populations) in different ecotypes of the island and seven microsatellite loci, we investigated neutral patterns of genetic variability and population structure at macro and microgeographic scales across *Ae. aegypti* populations from Martinique. Then, we investigated the association between insecticide resistance (temephos and deltamethrin) and 319 single nucleotide polymorphisms (SNPs) in 6 populations, in order to identify candidate genes for temephos and deltamethrin resistances in *Ae. aegypti*. We found that both neutral (microsatellites) and adaptive (SNPs) genetic variation is affected by insecticide pressure in *Ae. aegypti* populations in Martinique.

# Materials and Methods

## Sampling localities and bioassays

Two laboratory strains from French Polynesia (Bora-Bora) and Benin (SBE) and sixteen field-caught populations were used in the study. The two laboratory strains are susceptible to all insecticides and have been used as reference strains for resistance assays. *Ae. aegypti* was sampled from individual houses as larvae or pupae in sixteen localities of Martinique from January to March 2009 (Figure 1, Table S1). Each population was constituted with about one thousand of larvae originated from 10 to 20 larval breeding sites (domestic breeding habitats) located in a 250 meters radius. Larvae and pupae were reared until the adult stage (F0 generation) and then morphologically identified. The adults were kept in a cage for copulation and females were blood-fed on rabbit (protocol CEEA-LR-13002 agreed by The French National Committee of ethics for animal experimentation) to obtain a F1 progeny. When a convenient number of F1 eggs were obtained, F0 were stored at -20°C until molecular analysis. Larvae and adults obtained from the F1 progeny were used for bioassays. The chosen populations covered most of the island ecotypes (coastal, mountainous, rural and urban) and two of them were sampled in two small islands (Ilet Long and Ilet Anonyme) located one kilometer away from the Atlantic coast (Figure 1, Table S1). This sampling also covered a large range of insecticide pressures, from no treatment up to 81 interventions between 2006 and 2009 (Table S1). No specific permits were required for mosquito collection. We confirm that most of the sampling locations were not privately-owned. We asked the owners to sample in their back yards if needed. We also confirm that the locations were not protected in any way and that the field and laboratory studies did not involve endangered or protected species.

As described in [3], larval and adult bioassays were performed following WHO protocols [7] to measure the resistant ratios ($RR_{50}$ and $RR_{95}$) of twelve field-caught populations, the SBE strain (from Benin) and using the laboratory Bora-Bora strain as the susceptible reference strain. Larval and adult bioassays were carried out using temephos (Organophosphate) and deltamethrin (Pyrethroid), respectively.

## Microsatellite genotyping and population genetic analyses

Microsatellite genotyping and population genetic analyses were implemented at the macro-geographic level using 16 F0 samples of 30 females and 10 males. In 3 localities, Fort-de France, St-Anne and Vauclin, individuals from 10, 4 and 4 distinct larval breeding sites respectively, were used for a micro-geographic level analysis. DNA was extracted from adult mosquitoes of both sexes using 2% CTAB as described in [8]. DNA pellets were resuspended in 20 µL sterile water and stored at -20 °C until analysis. Individual genotypes were scored at a set of seven microsatellite markers 34/72, AEDGA, AED19, AEDC, A10, 38/38 and H08. The microsatellites and PCR conditions have previously been detailed [9,10]. DNA was amplified in a GeneAmp 9600 thermal cycler (PerkinElmer, USA). PCR products were analyzed using ABI Prism™ 3100
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Location of *Ae. aegypti* populations sampling sites in Martinique, and the proportion of assignment to each of the four genetic clusters identified by STRUCTURE.</figcaption></figure>

doi: 10.1371/journal.pone.0077857.g001

(Applied Biosystems), and genotyped with GeneMapper (Applied Biosystems) (for details see 9). The microsatellite 38/38 was not polymorphic in our populations and was not kept for further analyses.

Genetic diversity by locus and sample and overall was characterized by estimating the number of alleles and allelic richness (Rs, [11]) using the software FSTAT2.9.3.2 [12]. Linkage disequilibrium, deviations from Hardy-Weinberg equilibrium (HWE) and genetic differentiation indices were assessed using GENEPOP4.7 software [13]. FIS and FST estimates were calculated according to [14] and tested for statistical significance with exact tests available in GENEPOP4.7. The overall significance of multiple tests was estimated by Fisher's combined probability tests. Nominal significance levels for multiple testing were corrected using the sequential Bonferroni procedure [15]. Genetic isolation by geographic distance was assessed with Mantel's test in GENEPOP4.7, using the correlation between genetic and geographic distances by the regression of pairwise FST / (1 - FST) on the natural logarithm (Ln) of straight-line geographic distance [16]. Genetic structure was also studied at the macro and micro-geographic levels by examining the partition of neutral genetic variation using an AMOVA (Analysis of Molecular Variance, Arlequin 3.1 [17]). A bayesian clustering method implemented in STRUCTURE 2.3 [18] was used to assign individuals from all 16 populations to a pre-determined

number of clusters (K=1-16) based on multi-locus microsatellite data. For each run, a burn-in period of 10,000 steps was followed by 30,000 iterations under the admixture model and the assumption of correlated allele frequencies among populations. To assist the clustering, sampling locations were used as prior information (LOCPRIOR model). For each K of 1-16, 20 runs were performed. Estimated log probabilities (Ln P(D)) were averaged across runs and compared to determine the posterior probability of each K. K with the highest likelihood was retained. Pearson's correlation between genetic diversity estimates (allelic richness and FIS) and environmental variables (altitude, human density) and resistance ratio were tested using the "stat" package in the R software version 2.10 [19].

## SNP genotyping

Six out of the 16 populations (GRMN, ILA, ILO, LAM, RSAL and SAN) representative of the range of insecticide resistance levels observed in Martinique were analysed for 319 SNPs. 318 of these SNPs were developed in the *Ae. aegypti* laboratory strain Bora-Bora using transcriptomic data in [20]. Briefly, short cDNA tags were obtained from pools of larvae and sequenced using a Genome Analyzer I (Illumina Corporation, USA). The reads were filtered, mapped to Ae. aegypti genome assembly AaegL 1.1 (see VectorBase, http://aegypti vectorbase.org/) and 1520 SNPs with a coverage > 10x and a minimal allelic frequency of 5% were detected. The 318 SNPs used in this study were successfully transferred and genotyped in the natural *Ae. aegypti* populations of the Martinique Island [20]. They were distributed in 226 out of the 4,758 supercontigs in the reference genome. The knock down *kdr*-V1016I mutation in the voltage-gated sodium channel gene was also included in the SNP panel, as it was the only one found in a previous scan in Martinique populations among the four I1011M, I1011V, V1016I and V1016G *kdr* mutations [3] [21]. DNA was extracted from whole larvae or adults (30-40 individuals per population) from the F1 generation (also used for bioassays), using the DNeasy tissue Kit (Qiagen) according to the manufacturer's instructions. DNA concentration was determined for each sample with a NanoDrop™ ND-1000 (NanoDrop Technologies), normalized to 50ng/µl and genotyped using the DNA Illumina GoldenGate® array genotyping platform at the Genotyping Center GENTYANE in Clermont-Ferrand, France (http://www.ibisa.net/plateformes/GENTYANE,101.html). The fluorescent signal intensities of the specific alleles were determined by the BeadExpress Reader (Illumina Inc), and quantified and matched to genotype classes using GenomeStudio v3.1.14 (Illumina Inc) with manual adjustment when necessary. To ensure a good reliability of the genotyping data we used the GenCall50 (GC50) which indicates the reliability of each genotype call, the GenTrain which measures the cluster quality of the SNPs, and the call rate (CR) which represents the fraction of the samples that has been successfully genotyped for a given SNP. We kept SNPs with GC50 higher than 0.35, GenTrain higher than 0.50 and CR higher than 0.50.
## Detection of signature of positive selection (Fst outliers)

To detect SNPs displaying higher than expected differentiation under neutrality between the 6 populations (outliers) we used BayeScan v3 (freely available at http://cmpg.unibe.ch/software/bayescan/). This Bayesian method is based on a logistic regression model in which each logit value of genetic differentiation Fst between a given locus *i* in population *j* is decomposed as a linear combination of $\alpha_i$ and $\beta_j$, corresponding to a locus effect and to a population effect, respectively. To test for the significance of $\alpha_i$ for each locus, the posterior probability that $\alpha_i$ is different from zero is estimated using a Reversible Jump MCMC algorithm [22]. Priors for $\alpha_i$ and $\beta_j$ were Gaussian distributions with means 0 and -2 respectively and standard errors 1 and 1.8 as in [23]. The estimation of model parameters was automatically tuned on the basis of short pilot runs (10 pilot runs, length 2000). Preliminary tests indicated that a burn-in of 10,000 iterations was enough to attain convergence. We used a sample size of 10,000 and a thinning interval of 50, as suggested by [22], resulting in a total length of the chain of 500,000 iterations. Four independent runs were performed to check the consistency between detected outliers. The loci were ranked according to their estimated posterior probability value and all loci with a Bayes factor >3 were retained as outliers.

## Association with insecticide resistance

$RR_{95}$ and $RR_{50}$ were highly correlated ($r^2 = 0.92$ and 0.90, both *P* < 0.0001, for temephos and deltamethrin respectively) and we used $RR_{95}$ in association analysis in order to maximize the range of RR variation between the populations. Association analyses between $RR_{95}$ and SNPs were performed using the least squares fixed effects linear model implemented in the software TASSEL v3 [24]. This model tests for association between genotypes at each segregating marker and phenotype values, while accounting for population structure by using the 3 first principal components of a PCA derived from genetic markers as covariates. *P* values were corrected for multiple comparisons by permutation tests with 1000 permutations and 0.05 confidence level.

## Markers location in the reference genome and candidate genes identification

All SNPs detected either as *Fst* outliers, or as being significantly associated with a resistance phenotype, were blasted on the reference genome of *Ae. aegypti* (http://aegypti vectorbase.org/) to locate them inside or at proximity of potential candidate genes. The length of linkage disequilibrium can strongly vary from a genomic region to another and from a population to another, so that it is very difficult to assess an overall disequilibrium length in natural populations. In a closed population (a laboratory *Ae. aegypti* strain) under a standard pure-drift demographic model, overall linkage disequilibrium was estimated to be 20 Kb for neutral markers [25]. In a genomic region under selection, the amount of linkage disequilibrium can be considerably greater than that predicted for a neutral region. Furthermore, in natural population, gene flow can produce significant levels of

disequilibrium. Therefore we have arbitrarily chosen a conservative threshold of 500 Kb to search for candidate genes potentially under selection in the vicinity of outlier markers, and all the genes present up to 500 Kb each direction from the detected SNP were considered as potentially involved in the resistance phenotype and discussed.

## Results

### Bioassays

Bioassays showed that the twelve tested populations from Martinique were resistant to temephos and deltamethrin in comparison with the laboratory susceptible Bora-Bora strain and the SBE strain from Benin (Table 1). However, the resistance observed against temephos and deltamethrin was heterogeneous in the island. For temephos, $RR_{50}$ ranged from 13 fold (SJOS) to 36-fold (GMRN) and $RR_{95}$ ranged from 26-fold (SJOS) to 153-fold (VCLN). The populations from the two islets showed a lower resistance level to temephos with $RR_{50}$ of 5-fold for both islets, and $RR_{95}$ of 17 and 21-fold for ILA and ILO, respectively. Mortality after 24 h deltamethrin exposure ranged from 19% for AJPB to 90% for the RSAL population compared to the 100% mortality in the susceptible strains Bora-Bora and SBE. For each insecticide, resistance ratios were not significantly explained by the number of insecticide treatments, or by other environmental variables such as the type of locality (rural, urban), human density or altitude (Table S1).

### Microsatellite polymorphism analysis

**Genetic variability and Hardy-Weinberg expectations.** A total of 31 populations were successfully genotyped with the 7 microsatellites markers that were all polymorphic except 38/38 that was not kept for further analyses. Genetic diversity was relatively homogeneous among the different populations tested at both macro and micro-geographic levels (Table 2 and Table S2). Average allelic diversity was low (2.33 alleles/locus, range 1.73-2.38) as compared to the diversity found in African populations [10], but comparable to the diversity found in south American populations [26]. At the macro-geographic scale, when the populations were considered together, important heterozygote deficit were detected (34/72 and H08, and all loci; Table 2), suggesting that Martinique is not a single panmictic unit. However, within localities, departures from Hardy-Weinberg expectations and significant heterozygote deficits were detected only 3 times: in GRMN at locus 34/72, and in AJB and RS when all loci were considered together (Table 2). Overall, these results indicate that each sampled site can be considered as a panmictic population. At the micro-geographic level and considering the microsatellite markers individually, significant deviations were detected only in Fort-de-France samples (Table S2). Genetic diversity indices (*Nall*, *Rs* and $F_{IS}$) were not correlated (all *P* > 0.05) with the environmental variables measured (altitude, human density). However, allelic richness was negatively correlated with the level of resistance to deltamethrin (*t* = -2.66, *df* = 10, *P* = 0.02, *r* = -0.64) underlying a significant decrease of diversity in the most resistant populations. Furthermore, $F_{IS}$ was significantly correlated with the level of resistance to temephos (*t* = 2.43, *df* = 10).
<table><caption>**Table 1.** Resistance status of *Aedes aegypti* populations of Martinique to temephos (larvae) and deltamethrin (adults).</caption><thead><tr><th rowspan="2">Insecticide</th><th colspan="3">Temephos</th><th colspan="3">Deltamethrin</th><th rowspan="2">Mortality (24h)</th><th rowspan="2">n</th></tr><tr><th>RR<sub>50</sub> (ci RR<sub>50</sub>)</th><th>RR<sub>95</sub> (ci RR<sub>95</sub>)</th><th>n</th><th>RR<sub>50</sub> (ci RR<sub>50</sub>)</th><th>RR<sub>95</sub> (ci RR<sub>95</sub>)</th><th></th></tr></thead><tbody><tr><td>Bora<sup>1</sup></td><td>□</td><td>□</td><td>995</td><td>□</td><td>□</td><td>100%</td><td>204</td></tr><tr><td>SBE</td><td>0.7 (0.5-0.8)</td><td>0.83 (0.8-0.9)</td><td>502</td><td>1 (0.93-1)</td><td>0.9 (0.9-1)</td><td>100%</td><td>101</td></tr><tr><td>SPIER</td><td><b>13</b> (12-15)</td><td><b>54</b> (42-75)</td><td>498</td><td><b>4</b> (3.9-4.2)</td><td><b>5.8</b> (5.3-6.7)</td><td>70%</td><td>97</td></tr><tr><td>RSAL</td><td><b>29</b> (26-33)</td><td><b>82</b> (65-110)</td><td>501</td><td><b>3.7</b> (3.6-3.8)</td><td><b>4.7</b> (4.3-4.8)</td><td>90%</td><td>99</td></tr><tr><td>VCLN</td><td><b>27</b> (25-32)</td><td><b>153</b> (117-211)</td><td>699</td><td><b>5.7</b> (5.1-6.8)</td><td><b>7.1</b> (6-10.5)</td><td>20%</td><td>99</td></tr><tr><td>DIAM</td><td><b>23</b> (21-26)</td><td><b>90</b> (73-117)</td><td>702</td><td><b>5.4</b> (5-6.1)</td><td><b>6.2</b> (5.4-8.5)</td><td>34%</td><td>99</td></tr><tr><td>LAM</td><td><b>32</b> (28-36)</td><td><b>139</b> (113-176)</td><td>701</td><td><b>5.6</b> (6.8)</td><td><b>8</b> (6.5-12.8)</td><td>55%</td><td>94</td></tr><tr><td>GRMN</td><td><b>35</b> (30-45)</td><td><b>149</b> (102-217)</td><td>697</td><td><b>5.3</b> (4.9-6.1)</td><td><b>6.9</b> (5.9-9.6)</td><td>64%</td><td>101</td></tr><tr><td>SAN</td><td><b>19</b> (18-22)</td><td><b>56</b> (45-73)</td><td>699</td><td><b>6.7</b> (5.7-8.6)</td><td><b>12.2</b> (9-21.5)</td><td>42%</td><td>98</td></tr><tr><td>SJOS</td><td><b>13</b> (11-16)</td><td><b>27</b> (16-42)</td><td>696</td><td><b>5.8</b> (5.2-7)</td><td><b>8.5</b> (6.8-13.4)</td><td>29%</td><td>95</td></tr><tr><td>FDF</td><td><b>15</b> (14-16)</td><td><b>36</b> (30-45)</td><td>702</td><td><b>6.4</b> (5.5-8.6)</td><td><b>9.4</b> (7-18.2)</td><td>23%</td><td>101</td></tr><tr><td>AJPB</td><td><b>15</b> (13-17)</td><td><b>54</b> (43-74)</td><td>600</td><td><b>5.9</b> (5.2-6.9)</td><td><b>10.3</b> (8.2-15.8)</td><td>19%</td><td>103</td></tr><tr><td>Bora<sup>2</sup></td><td>□</td><td>□</td><td>602</td><td>□</td><td>□</td><td>100%</td><td>100</td></tr><tr><td>ILA</td><td><b>5</b> (3-7)</td><td><b>17</b> (7-39)</td><td>900</td><td><b>6.7</b> (6.6-6.8)</td><td><b>7.9</b> (-7.5-9)</td><td>27%</td><td>100</td></tr><tr><td>IOLO</td><td><b>5</b> (4-5)</td><td><b>21</b> (17-29)</td><td>900</td><td><b>5.1</b> (4.9-5.4)</td><td><b>7.9</b> (7.8-8.4)</td><td>38%</td><td>101</td></tr></tbody></table>

1. Bioassays were performed in two series corresponding to two generations of the laboratory susceptible Bora-Bora strain (Bora and Bora2). RR: resistant ratio; ci: confidence interval; LC₅₀ and LC₉₅ in mg/L were 0.003 and 0.06 respectively for Bora and 0.005 and 0.01mg/L for Bora2 with temephos. KD₅₀ and KD₉₅ were 14 and 21 min respectively (11 and 18min for Bora2) with deltamethrin.

doi: 10.1371/journal.pone.0077857.t001

= 10, P = 0.035, r = 0.609) showing an increased consanguinity in populations more resistant to temephos.

## Genetic differentiation and population structure

The global genetic differentiation index *F*<sub>ST</sub> was low (0.079) but the populations were significantly structured both at the macro-geographic and at the micro-geographic scales (Table 3). Pairwise *F*<sub>ST</sub> values were significant for 69% of the population pairs tested. Genetic and geographic distances were correlated at the Martinique scale (Mantel test; *P* = 0.009) and at the local scale in Vauclin (*P* = 0.042), but not in Fort-de-France and Sainte-Anne (*P* > 0.5). Four genetic clusters were identified in Martinique (Figure 1), one being mostly represented in the two islets and their closest eastern Vauclin population on the main island, and another corresponding to populations around the gulf of Fort-de-France (urban area). All the populations in Martinique were highly admixed, suggesting important gene flow across populations. The most homogenous populations were the two islets assigned >80% to one cluster and Fond-St-Denis assigned >90% to another cluster.

## SNP diversity

Average allelic diversity per locus per population was 1.76 (range 1.6-1.9) and the 6 populations were at the HW equilibrium except for a slight deficit in heterozygotes in Rivière Salée. The global *F*<sub>ST</sub> = 0.031 was low, suggesting large gene flow among these populations; the structure analysis based on SNPs identified the same 4 main clusters as identified based on microsatellites (data not shown). Pairwise *F*<sub>ST</sub> were generally smaller for SNPs than for microsatellites (Table S3).

## Detection of *F*<sub>ST</sub> outliers and association tests

Five out of the 319 SNPs were identified as *F*<sub>ST</sub> outliers across the 6 populations tested in Martinique (Table 4), including the *kdr*-V1016I mutation in the voltage-gated sodium channel. However, the *kdr* mutation was not detected as significantly associated with deltamethrin resistance. The resistance allele is found at high frequency in all populations on the main island (85% to 92%), confirming previous report [3], which may explain the failure of the GLM regression to detect significant association between deltamethrin resistance and the frequency of the resistant allele. From the five *F*<sub>ST</sub> outliers, two were significantly associated with deltamethrin resistance, and one with temephos resistance. Five other markers were associated with deltamethrin resistance and four other were associated with temephos resistance (Table 4). Taken separately, each marker explained from 9 to 17% of the resistance variability observed between populations.

## Markers location in the reference genome and candidate genes identification

Most of the SNPs identified were either inside a gene annotated in Vectorbase, or in the flanking 3' extremity of a gene, which is not surprising because SNPs have been designed based on transcriptome sequences anchored on the 3' end of transcripts [27]. The seven SNPs associated with deltamethrin resistance corresponded to six genes encoding respectively a leucin-rich membrane protein, two proteins associated with ribosomes, and three hypothetical proteins of unknown function. All SNPs were synonymous. Among the five SNPs associated with temephos resistance, two were non-synonymous mutations in an alkaline phosphatase and in a sulfotransferase, two were synonymous mutations in an
<table><caption>**Table 2.** Genetic variability and goodness of fit to Hardy-Weinberg expectation in *Ae. aegypti* populations at the macro-geographic level in Martinique.</caption><thead><tr><th></th><th>1 n=30</th><th>2n=30</th><th>3 n=30</th><th>4 n=30</th><th>5n=30</th><th>6 n=30</th><th>7 n=30</th><th>13 n=30</th><th>18 n=30</th><th>19 n=30</th><th>20 n=30</th><th>21 n=30</th><th>22 n=30</th><th>26 n=30</th><th>30 n=30</th><th>31 n=30</th><th>all n=30</th></tr></thead><tbody><tr><td><strong>3472</strong></td><td><em>N</em><sub>all</sub> 2</td><td>4</td><td>2</td><td>3</td><td>5</td><td>3</td><td>3</td><td>2</td><td>3</td><td>3</td><td>3</td><td>3</td><td>2</td><td>2</td><td>2</td><td>2</td><td>5</td></tr><tr><td></td><td><em>R</em><sub>s</sub> 1.33</td><td>3.29</td><td>1.36</td><td>2.96</td><td>3.26</td><td>2.36</td><td>2.76</td><td>2.00</td><td>2.72</td><td>2.71</td><td>2.78</td><td>2.71</td><td>1.78</td><td>2.00</td><td>2.00</td><td>2.00</td><td>2.66</td></tr><tr><td></td><td><em>F</em><sub>IS</sub> 0.000</td><td>0.285</td><td>0.000</td><td>0.233</td><td>0.264</td><td><strong>0.774</strong></td><td>-0.147</td><td>-0.364</td><td>0.137</td><td>-0.114</td><td>0.099</td><td>-0.130</td><td>-0.042</td><td>-0.515</td><td>-0.862</td><td>-0.262</td><td><strong>-0.009</strong></td></tr><tr><td><strong>Aedga</strong></td><td><em>N</em><sub>all</sub> 1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>2</td><td>1</td><td>1</td><td>1</td><td>1</td><td>4</td></tr><tr><td></td><td><em>R</em><sub>s</sub> 1.00</td><td>1.00</td><td>1.00</td><td>1.00</td><td>1.00</td><td>1.00</td><td>1.00</td><td>1.00</td><td>1.00</td><td>1.00</td><td>1.00</td><td>1.61</td><td>1.00</td><td>1.00</td><td>1.00</td><td>1.00</td><td>1.07</td></tr><tr><td></td><td><em>F</em><sub>IS</sub> NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>-0.020</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>0.328</td></tr><tr><td><strong>A10</strong></td><td><em>N</em><sub>all</sub> 3</td><td>4</td><td>3</td><td>3</td><td>2</td><td>4</td><td>4</td><td>4</td><td>3</td><td>3</td><td>2</td><td>4</td><td>3</td><td>3</td><td>4</td><td>4</td><td>6</td></tr><tr><td></td><td><em>R</em><sub>s</sub> 3.00</td><td>3.55</td><td>3.00</td><td>2.56</td><td>2.00</td><td>3.62</td><td>2.89</td><td>2.67</td><td>2.56</td><td>2.56</td><td>2.00</td><td>3.15</td><td>2.71</td><td>2.81</td><td>3.62</td><td>3.10</td><td>3.23</td></tr><tr><td></td><td><em>F</em><sub>IS</sub> -0.397</td><td>0.004</td><td>-0.103</td><td>-0.324</td><td>-0.063</td><td>0.138</td><td>0.092</td><td>0.077</td><td>-0.202</td><td>-0.217</td><td>0.392</td><td>0.020</td><td>-0.037</td><td>0.101</td><td>0.288</td><td>-0.021</td><td>-0.035</td></tr><tr><td><strong>Aed19</strong></td><td><em>N</em><sub>all</sub> 2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td></td><td><em>R</em><sub>s</sub> 2.00</td><td>2.00</td><td>2.00</td><td>2.00</td><td>2.00</td><td>2.00</td><td>2.00</td><td>2.00</td><td>2.00</td><td>2.00</td><td>2.00</td><td>2.00</td><td>2.00</td><td>2.00</td><td>2.00</td><td>2.00</td><td>2.00</td></tr><tr><td></td><td><em>F</em><sub>IS</sub> 0.366</td><td>0.116</td><td>-0.166</td><td>-0.241</td><td>0.102</td><td>-0.102</td><td>-0.080</td><td>-0.244</td><td>0.067</td><td>-0.120</td><td>-0.017</td><td>0.191</td><td>0.198</td><td>0.231</td><td>-0.200</td><td>0.087</td><td>-0.010</td></tr><tr><td><strong>Aedc</strong></td><td><em>N</em><sub>all</sub> 3</td><td>2</td><td>1</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>3</td><td>3</td><td>3</td><td>2</td><td>4</td><td>4</td><td>1</td><td>6</td></tr><tr><td></td><td><em>R</em><sub>s</sub> 2.59</td><td>1.89</td><td>1.00</td><td>1.93</td><td>1.99</td><td>1.93</td><td>1.97</td><td>1.38</td><td>1.97</td><td>2.04</td><td>2.51</td><td>1.89</td><td>1.76</td><td>3.06</td><td>2.03</td><td>1.00</td><td>2.16</td></tr><tr><td></td><td><em>F</em><sub>IS</sub> -0.118</td><td>-0.077</td><td>NA</td><td>-0.098</td><td>-0.191</td><td>-0.098</td><td>-0.137</td><td>0.000</td><td>0.151</td><td>-0.040</td><td>0.380</td><td>0.663</td><td>-0.040</td><td>-0.008</td><td>0.325</td><td>NA</td><td>-0.012</td></tr><tr><td><strong>H08</strong></td><td><em>N</em><sub>all</sub> 3</td><td>3</td><td>2</td><td>2</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>2</td><td>3</td><td>3</td><td>3</td><td>2</td><td>3</td><td>3</td><td>4</td></tr><tr><td></td><td><em>R</em><sub>s</sub> 2.45</td><td>2.57</td><td>2.00</td><td>2.00</td><td>2.36</td><td>3.00</td><td>2.92</td><td>2.63</td><td>2.98</td><td>2.00</td><td>3.00</td><td>2.73</td><td>2.82</td><td>2.00</td><td>2.84</td><td>2.37</td><td>2.85</td></tr><tr><td></td><td><em>F</em><sub>IS</sub> 0.649</td><td>0.004</td><td>-0.125</td><td>-0.250</td><td>-0.288</td><td>-0.251</td><td>-0.039</td><td>-0.192</td><td>0.345</td><td>-0.012</td><td>0.197</td><td>0.345</td><td>0.217</td><td>0.386</td><td>0.146</td><td>-0.157</td><td><strong>0.001</strong></td></tr><tr><td><strong>Mean across all loci</strong></td><td><em>N</em><sub>all</sub> 2.33</td><td>2.67</td><td>1.83</td><td>2.17</td><td>2.50</td><td>2.50</td><td>2.50</td><td>2.33</td><td>2.33</td><td>2.33</td><td>2.33</td><td>2.83</td><td>2.17</td><td>2.33</td><td>2.67</td><td>2.17</td><td>4.50</td></tr><tr><td></td><td><em>R</em><sub>s</sub> 2.06</td><td>2.38</td><td>1.73</td><td>2.08</td><td>2.10</td><td>2.32</td><td>2.26</td><td>1.95</td><td>2.21</td><td>2.05</td><td>2.22</td><td>2.35</td><td>2.01</td><td>2.15</td><td>2.25</td><td>1.91</td><td>2.33</td></tr><tr><td></td><td><em>F</em><sub>IS</sub> <strong>0.004</strong></td><td>0.069</td><td>-0.124</td><td>-0.134</td><td>-0.040</td><td>0.098</td><td>-0.041</td><td>-0.169</td><td>0.106</td><td>-0.111</td><td><strong>0.190</strong></td><td>0.152</td><td>0.087</td><td>0.029</td><td>-0.102</td><td>-0.088</td><td><strong>-0.014</strong></td></tr></tbody></table>

1 Ajoupa-Bouillon; 2: Saint-Pierre; 3: Fond-Saint-Denis; 4: Sainte-Marie; 5: Trinité; 6: Gros-Morne; 7: Saint-Joseph; 13: Fort-de-France; 18: Lamentin; 19: Saint-Esprit; 20: Rivière Salée; 21: Diamant; 22: Sainte-Anne 1; 26: Vauclin 1; 30: Ilet Anonyme; 31: Ilet Long. n: sample size. *N*<sub>all</sub>: number of scored alleles; *R*<sub>s</sub>: Allele richness; *F*<sub>IS</sub>: Inbreeding coefficient; In bold: significant deficit in heterozygotes (*P* < 0.05) after Bonferroni correction.
doi: 10.1371/journal.pone.0077857.t002

aldehyde dehydrogenase and an unknown protein, and the last one was in a non-coding region. Genes located on the same supercontig and within ± 500 Kb from a marker associated with resistance are presented in Table S4. Among these, genes encoding enzymes potentially involved in detoxification processes such as transferases, phosphatases and dehydrogenases were identified. Genes encoding transcription factors and miRNA potentially regulating gene expression were also identified.

Discussion

Population structure at local and global scales

Our study showed a low genetic diversity in *Ae. aegypti* populations from Martinique compared with levels observed in other parts of the world, especially in native African populations [9,10], but was consistent with a previous study based on isoenzymes in this Island [28]. The loss of genetic diversity observed here can be explained by the founder effect that occurred during the introduction of the mosquito in the Caribbean. Furthermore, this low genetic diversity may have been maintained since its introduction in the 17th century by several factors such as recurrent demographic crashes linked to the specific environmental variations in Martinique (i.e. rainfall and droughts) and more recently by anthropogenic actions such as vector control. Because of the recurrent

<table><caption>**Table 3.** Genetic differentiation estimated with 6 microsatellite markers at two geographic levels.</caption><thead><tr><th></th><th>Number of populations</th><th><em>F</em><sub>ST</sub><sup>a</sup></th><th><em>P</em></th><th>% pairwise differentiation<sup>b</sup></th></tr></thead><tbody><tr><td>Macro-geographic level</td><td>16</td><td>0.0719</td><td><strong>&lt; 10<sup>-4</sup></strong></td><td>69</td></tr><tr><td>Micro-geographic level</td><td></td><td></td><td></td><td></td></tr><tr><td>Fort de France</td><td>10</td><td>0.0839</td><td><strong>&lt; 10<sup>-4</sup></strong></td><td>69</td></tr><tr><td>Sainte-Anne</td><td>4</td><td>0.0741</td><td><strong>&lt; 10<sup>-4</sup></strong></td><td>100</td></tr><tr><td>Vauclin</td><td>4</td><td>0.0532</td><td><strong>&lt; 10<sup>-4</sup></strong></td><td>83</td></tr></tbody></table>

a. First indices estimated by population groups and *P*, probability associated with Fisher's exact test

b. Average pairwise population percentages showing a significant genetic differentiation after Bonferroni correction using a threshold of *P*<0.05
doi: 10.1371/journal.pone.0077857.t003

Dengue epidemics occurring in Martinique, vector control (i.e. source reduction and insecticide applications) is intense and may induce local bottleneck effects and a reduction of the genetic diversity in *Ae. aegypti* populations. In support to this hypothesis, we found significant correlations between neutral genetic diversity estimates and insecticide resistance ratio: allelic richness was negatively correlated to deltamethrin
<table><caption>**Table 4.** Summary of the SNPs detected as outliers in BayeScan analysis and/or as associated with a resistance phenotype.</caption><thead><tr><th>SNP name<br>(supercontig,gene)</th><th>Bayescan outlier</th><th>Associated<br>resistance<br>phenotype</th><th>perm_pR2</th><th>Supercontig<br>number</th><th>gene</th><th>function</th><th>5'-3' position</th><th>alleles</th><th>aa change</th><th>Supercontig size</th></tr></thead><tbody><tr><td>1.4_AAELO00245</td><td>yes</td><td>deltamethrin</td><td>0.006 0.093 1.4</td><td></td><td>AAELO000245</td><td>conserved hypothetical protein</td><td>5094449</td><td>[T/C]</td><td></td><td>5177111</td></tr><tr><td>1.250_AAELO07363</td><td>yes</td><td>deltamethrin</td><td>0.001 0.143 1.250</td><td></td><td>AAELO007363 (3' UTR)</td><td>leucine-rich transmembrane protein</td><td>337073</td><td>[A/C]</td><td></td><td>1628417</td></tr><tr><td>1.636_NA0299</td><td>yes</td><td>temephos</td><td>0.001 0.170 1.636</td><td></td><td></td><td></td><td>381176</td><td>[T/C]</td><td></td><td>653825</td></tr><tr><td>1.200_NA0113</td><td>yes</td><td></td><td>1.200</td><td></td><td></td><td></td><td>1695291</td><td>[T/G]</td><td></td><td>1852562</td></tr><tr><td>1.186_kdr-V1016I</td><td>yes</td><td></td><td>1.186</td><td></td><td></td><td></td><td>115680</td><td>[A/G]</td><td>Ile-Val</td><td>1957664</td></tr><tr><td>1.1168_AAELO14562_A</td><td>deltamethrin</td><td>0.001 0.135 1.1168</td><td></td><td></td><td>AAELO006019</td><td>voltage-gated sodium channel</td><td>115680</td><td>[A/G]</td><td></td><td>191961</td></tr><tr><td>1.1168_AAELO14562_B</td><td>deltamethrin</td><td>0.006 0.092 1.1168</td><td></td><td></td><td>AAELO14562</td><td>60S ribosomal protein L12</td><td>182837</td><td>[T/C]</td><td></td><td>191961</td></tr><tr><td>1.541_AAELO11089</td><td>deltamethrin</td><td>0.003 0.105 1.541</td><td></td><td></td><td>AAELO11089</td><td>Ribonucleoprotein</td><td>183099</td><td>[A/G]</td><td></td><td>819817</td></tr><tr><td>1.68_AAELO02794</td><td>deltamethrin</td><td>0.024 0.080 1.68</td><td></td><td></td><td>AAELO02794 (3' UTR)</td><td>hypothetical protein</td><td>1416497</td><td>[A/C]</td><td></td><td>2950385</td></tr><tr><td>1.70_AAELO02875</td><td>deltamethrin</td><td>0.029 0.101 1.70</td><td></td><td></td><td>AAELO002875</td><td>hypothetical protein</td><td>79260</td><td>[T/G]</td><td></td><td>2929944</td></tr><tr><td>1.83_AAELO03298</td><td>temephos</td><td>0.016 0.115 1.83</td><td></td><td></td><td>AAELO003298</td><td>alkaline phosphatase</td><td>254086</td><td>[T/A]</td><td>Thr--&gt;Ser</td><td>2974912</td></tr><tr><td>1.414_AAELO09633</td><td>temephos</td><td>0.001 0.154 1.414</td><td></td><td></td><td>AAELO009633</td><td>conserved hypothetical protein</td><td>1016060</td><td>[T/C]</td><td></td><td>1096716</td></tr><tr><td>1.355_AAELO08898</td><td>temephos</td><td>0.03 0.092 1.355</td><td></td><td></td><td>AAELO008898</td><td>sulfotransferase</td><td>665884</td><td>[A/C]</td><td>Asp--&gt;Lys</td><td>1256691</td></tr><tr><td>1.1002_AAELO14080</td><td>temephos</td><td>0.002 0.131 1.1002</td><td></td><td></td><td>AAELO14080</td><td>Aldehyde deshydrogenase</td><td>176807</td><td>[A/G]</td><td></td><td>296325</td></tr></tbody><tfoot><tr><td colspan="11">a. Fst indices estimated by population groups and P, probability associated with Fisher's exact test</td></tr><tr><td colspan="11">b. Average pairwise population percentages showing a significant genetic differentiation after Bonferroni correction using a threshold of P&lt;0.05</td></tr><tr><td colspan="11">doi: 10.1371/journal.pone.0077857.t004</td></tr></tfoot></table>
resistance, suggesting that population reduction following insecticide treatments leads to allelic loss. Furthermore, F<sub>IS</sub> was positively correlated with temephos resistance, meaning that more resistant populations are more consanguineous, a result expected if only individuals carrying resistance alleles and their progeny can survive.

Despite an overall low genetic differentiation among sites, most pairwise *F*<sub>ST</sub> were significant and there was a significant pattern of isolation-by-distance (IBD) at the whole Martinique scale. This suggests that gene flow mostly occurs at a local scale with high admixture occurring among neighbour sites. Furthermore, four genetic clusters were identified, with marked differences in the proportion of assignation of individuals to each cluster from a site to another. Northern populations were mostly assigned to one cluster (red on Figure 1) and Southern populations were more admixed with a predominance of another cluster (green on Figure 1). The genetic drift observed in Fond-Saint-Denis (mostly red cluster) and in the Islets (mostly blue cluster) could be explained by their geographic isolation, with restricted access in the mountainous area and only accessible by boat, respectively. At a more local scale, the IBD pattern was marginally significant in only one out of three sites analysed (Vauclin, *P* = 0.042), supporting the effect of genetic drift acting on mosquito populations. Genetic drift at local scale could be due to the recurrent insecticide applications leading to local bottlenecks. Furthermore, flight activity of *Ae. aegypti* is limited to few hundred meters [29] which is likely to increase the neutral genetic differentiation between the populations through genetic drift. No significant correlation was found between diversity estimates and environmental variables such as altitude or human density.

## Outlier detection and association with insecticide resistance

The genome scan analysis was performed on 6 populations representative of the diverse mosquito habitats found across Martinique Island, and of the most extreme insecticide resistance ratio values observed. BayeScan detected five outliers out of 319 SNPs, all with positive *αi* values, indicative of directional selection. These included the knockdown resistance *kdr* mutation (V1016I), which is a well-characterized mechanism of resistance to pyrethroid insecticides in insects (including *Ae. aegypti*) and is caused by a point mutation of the *para*-type sodium channel [30], and three other genomic regions that were associated with resistance to either temephos or deltamethrin. Therefore, only one out of five BayeScan outliers was not related to insecticide pressure. This outlier was not found significantly associated to altitude or to human density, but it might be linked to other environmental or ecological pressures that were not measured in the present study, including climate and/or agricultural practices. Despite a highly diversified landscape in Martinique, and a complex life cycle involving an aquatic larval stage and an aerial adult stage, the mosquito genome appears to be more constrained by insecticide pressure than by any other environmental pressure. Although identified as an outlier by the BayeScan approach, the *kdr* mutation V1016I was not significantly associated with deltamethrin resistance; this could be due to

the fact that this mutation was present at high frequency in all the main island populations (>85%, as also found in [3]), decreasing the power of the association study. Furthermore, other *kdr* mutations such as the *kdr* mutation F1534C found at a relative high frequency in Grand Cayman Island and shown to be linked with DDT and permethrin resistance [31] might also be involved in pyrethroid resistance in Martinique, and their occurrence should be searched for in Martinique. Association study revealed that nine additional SNPs not detected as outliers by BayeScan were significantly correlated with resistance to either insecticide. BayeScan is a conservative approach and is sensitive to the level of genetic differentiation across populations: if genetic differentiation is too small, the Reversible Jump MCMC algorithm used to estimate posterior probabilities can fail to reach convergence, and if it is too high, the power to detect outlier loci is low. In the association study, population structure is taken into account as a fixed effect in the general linear model fitting phenotypic to genotypic value at each locus under investigation. As it is a correlative approach, it is not sensitive to population processes such as drift, migration, and demographic fluctuations. On another hand, this correlative approach might detect false positives. We independently detected two markers significantly associated with deltamethrin resistance located at two different positions in the same gene (AAEL014562) suggesting that this association was not discovered by chance. A further way to test for the role of outliers in insecticide resistance would be to compare their relative frequencies in deads and alives after bioassays.

## Candidate genes for resistance to deltamethrin and temephos

Among the SNPs significantly associated with resistance, only two were non-synonymous. Both were associated with resistance to temephos, and both involved genes encoding for enzymes. One induced a change from a threonine to a serine in an alkaline phosphatase while the other induced a change from a lysine to an asparagine in a sulfotransferase (Table S4). Alkaline phosphatases (ALP) catalyse the transfer of a phosphate group to water (hydrolysis) or alcohol (transphosphorylation) [32]. ALP-based assays have been used to detect pesticides in the environment, including organophosphates [33] and might thus be involved in their detoxification. Sulfotransferases are known to be involved in phase II detoxification processes by adding a sulphate group to various substrates to facilitate their excretion [34,35]. These two genes (and point mutations) were never described previously as being linked to temephos resistance, but they represent valuable candidates for further functional validation. All the other SNPs significantly associated with insecticide resistance were synonymous, and therefore not associated to protein variations potentially contributing to resistance, but rather in the regulation of gene expression (transcription factors, miRNA) or being in linkage disequilibrium with a functional mutation (hitch-hiking with the adaptive mutation). Because we genotyped individuals from natural populations, the length of linkage disequilibrium is likely to be low, and we restricted our search for candidates to 500 Kb upstream or
downstream from each SNP (Table S4). Among candidates, we found several genes potentially involved in detoxification processes such as phosphatases, dehydrogenases and transferases. In mammals, alcohol and aldehyde dehydrogenases (ADH and ALDH) can contribute to the metabolism of pyrethroids through further processing of the metabolites generated by esterases, phenoxybenzyl alcool and phenoxybenzyl aldehyde, to phenoxibenzoic acid [36,37]. However, no classical detoxification enzymes found over-transcribed in insecticide-resistant populations from Martinique were identified [3,4]. Indeed, recent data suggest that cytochrome P450 monooxygenases play a crucial role in pyrethroid detoxification in mosquitoes [38] while carboxylesterases are known to confer resistance to organophosphates [39]. It is possible that mutations causing the over-transcription of these detoxification genes in Martinique are located too far away from our markers, or that regulators of detoxification enzymes are located far away from the gene they regulate (trans-regulators); indeed, we identified potential regulators of transcription situated close to SNPs associated with resistance (Table S4). Although this genome scan provides valuable data to identify insecticide resistance genes in mosquitoes, the main limits are the partial annotation of the *Ae. aegypti* genome, with many genes annotated as 'hypothetical proteins' of unknown function, and the poor assembly of the genome with 4,758 supercontigs. As a result, we restricted our search to genes located 500 Kb in either direction in the same supercontig as the SNP associated to resistance, a relatively short distance, that was even reduced for markers located closer than 500 Kb to supercontig boundaries (see Table S4).

Finally, there was no SNP or genomic region associated with both resistance to deltamethrin and to temephos, which suggests that different resistance mechanisms are evolving in natural populations against these two commonly used insecticides, temephos being used against larvae and deltamethrin against adults. These two insecticides have different targets, the acetylcholinesterase and the sodium

channel respectively, and the finding of different genes involved in the two types of resistance is not surprising. Functional validation of the precise role in each type of insecticide resistance for each candidate gene detected in this study is the next step to be undertaken.

## Supporting Information

<caption>**Table S1. GPS coordinates and characteristics of the 16 populations sampled in Martinique.**
(DOCX)</caption>

<caption>**Table S2. Genetic variability and goodness of fit to Hardy-Weinberg expectation in *Ae. aegypti* populations at the micro-geographic level in Martinique.**
(DOC)</caption>

<caption>**Table S3. Comparison of pairwise Fst obtained with microsatellites and SNPs for six populations of Martinique.**
(DOCX)</caption>

<caption>**Table S4. Genes up to 500 kb each direction from the SNP.**
(XLSX)</caption>

## Acknowledgements

We thank Serge Sélior for his help for mosquito collections, Erica Bianco for laboratory help, and Charles Poncet and Lydia Jaffrelo for SNPs genotyping.

## Author Contributions

Conceived and designed the experiments: MP JPD VC LD. Performed the experiments: SM MP. Analyzed the data: SM MP CP CP LD. Contributed reagents/materials/analysis tools: CP AY FC LD. Wrote the manuscript: SM MP LD.

## References

1. Soper FL (1977) Ventures in World Health. D.C.: Washington.

2. Luis San Martin J, Brathwaite O, Zambrano B, Orlando Solorzano J, Bouckenooghe A et al. (2010) The Epidemiology of Dengue in the Americas Over the Last Three Decades: A Worrisome Reality. Am J Trop Med Hyg 82: 128-135. doi:10.4269/ajtmh.2010.09-0346. PubMed: 20065008.

3. Marcombe S, Mathieu RB, Pocquet N, Riaz M-A, Poupardin R et al. (2012) Insecticide Resistance in the Dengue Vector Aedes aegypti from Martinique: Distribution, Mechanisms and Relations with Environmental Factors. PLOS ONE 7: e30989. PubMed: 22363529.

4. Marcombe S, Carron A, Darriet F, Etienne M, Agnew P et al. (2009) Reduced Efficacy of Pyrethroid Space Sprays for Dengue Control in an Area of Martinique with Pyrethroid Resistance. Am J Trop Med Hyg 80: 745-751. PubMed: 19407118.

5. Marcombe S, Darriet F, Tolosa M, Agnew P, Duchon S et al. (2011) Pyrethroid Resistance Reduces the Efficacy of Space Sprays for Dengue Control on the Island of Martinique (Caribbean). PLOS Neglected Trop Dis 5: e1202. PubMed: 21713017.

6. Nkya TE, Akhouayri I, Kisinza W, David J-P (2013) Impact of environment on mosquito response to pyrethroid insecticides: Facts, evidences and prospects. Insect Biochem Mol Biol 43: 407-416. doi: 10.1016/j.ibmb.2012.10.006. PubMed: 23123179.

7. WHO (2005) Guidelines for laboratory and field testing of mosquito larvicides. Geneva: World Health Organization.

8. Morlais I, Ponçon N, Simard F, Cohuet A, Fontenille D (2004) Intraspecific nucleotide variation in Anopheles gambiae: New insights into the biology of malaria vectors. Am J Trop Med Hyg 71: 795-802. PubMed: 15642974.

9. Paupy C, Brengues C, Kamgang B, Hervé J-P, Fontenille D et al. (2008) Gene flow between domestic and sylvan populations of Aedes aegypti (Diptera: Culicidae) in North Cameroon. J Med Entomol 45: 391-400. doi:10.1603/0022-2585(2008)45[391:GFBDAS]2.0.CO;2. PubMed: 18533431.

10. Paupy C, Brengues C, Ndiath O, Toty C, Hervé J-P et al. (2010) Morphological and genetic variability within Aedes aegypti in Niakhar, Senegal. Infect Genet Evol 10: 473-480. doi:10.1016/j.meegid. 2010.03.001. PubMed: 20223297.

11. el Mousadik A, Petit RJ (1996) Chloroplast DNA phylogeography of the argan tree of Morocco. Mol Ecol 5: 547-555. doi:10.1046/j.1365-294X. 1996.00123.x. PubMed: 8794562.

12. Goudet J (1995) FSTAT (Version 1.2): A computer program to calculate F-statistics. J Hered 86: 485-486.

13. Raymond M, Roussel F (1995) GENEPOP (Version 1.2) - Population genetics software for exact tests and oecumenicism. J Hered 86: 248-249.

14. Weir BS, Cockerham CC (1984) Estimating F-statistics for the analysis of population structure. Evolution 38: 1358-1370. doi:10.2307/2408641.
15. Holm S (1979) A simple sequentially rejective multiple test procedure. Scand J Statist 6: 65-70.

16. Rousset F (1997) Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance. Genetics 145: 1219-1228. PubMed: 9093870.

17. Excoffier L, Laval G, Schneider S (2005) Arlequin ver. 3.0: An integrated software package for population genetics data analysis. Evol Bioinform Online 1: 47-50.

18. Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155: 945-959. PubMed: 10835412.

19. Team RDC (2008) R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing.

20. Paris M, Marcombe S, Coissac E, Corbel V, David J-P et al. (2013) Investigating the genetics of Bti resistance using mRNA tag sequencing: application on laboratory strains and natural populations of the dengue vector Aedes aegypti. Evol Applications. doi:10.1111/eva. 12082.

21. Brengues C, Hawkes NJ, Chandre F, McCarroll L, Duchon S et al. (2003) Pyrethroid and DDT cross-resistance in Aedes aegypti is correlated with novel mutations in the voltage-gated sodium channel gene. Med Vet Entomol 17: 87-94. doi:10.1046/j. 1365-2915.2003.00412.x. PubMed: 12680930.

22. Foll M, Gaggiotti O (2008) A Genome-Scan Method to Identify Selected Loci Appropriate for Both Dominant and Codominant Markers: A Bayesian Perspective. Genetics 180: 977-993. doi:10.1534/genetics. 108.092221. PubMed: 18780740.

23. Beaumont MA, Balding DJ (2004) Identifying adaptive genetic divergence among populations from genome scans. Mol Ecol 13: 969-980. doi:10.1111/j.1365-294X.2004.02125.x. PubMed: 15012769.

24. Bradbury PJ, Zhang Z, Kroon DE, Casstevens TM, Ramdoss Y et al. (2007) TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23: 2633-2635. doi:10.1093/ bioinformatics/btm308. PubMed: 17586829.

25. Paris M, Despres L (2012) Identifying insecticide resistance genes in mosquito by combining AFLP genome scans and 454 pyrosequencing. Mol Ecol 21: 1672-1686. doi:10.1111/j.1365-294X.2012.05499.x. PubMed: 22348648.

26. Paupy C, Le Goff G, Brengues C, Guerra M, Revollo J et al. (2012) Genetic structure and phylogeography of Aedes aegypti, the dengue and yellow-fever mosquito vector in Bolivia. Infect Genet Evol 12: 1260-1269. doi:10.1016/j.meegid.2012.04.012. PubMed: 22522103.

27. Paris M, Melodelima C, Coissac E, Tetreau G, Reynaud S et al. (2012) Transcription profiling of resistance to Bti toxins in the mosquito Aedes

aegypti using next-generation sequencing. J Invert Pathol 109: 201-208. doi:10.1016/j.jip.2011.11.004. PubMed: 22115744.

28. Yebakima A, Charles C, Mousson L, Vazeille M, Yp-Tcha MM, et al. (2004) Genetic heterogeneity of the dengue vector Aedes aegypti in Martinique. Tropical; Medicine & International. Health 9: 582-587.

29. Honorio NA, Silva WD, Leite PJ, Goncalves JM, Lounibos LP et al. (2003) Dispersal of Aedes aegypti and Aedes albopictus (Diptera : Culicidae) in an urban endemic dengue area in the State of Rio de Janeiro, Brazil. Mem Inst Oswaldo Cruz 98: 191-198. doi:10.1590/ S0074-02762003000200005.

30. Saavedra-Rodriguez K, Urdaneta-Marquez L, Rajatileka S, Moulton M, Flores AE et al. (2007) A mutation in the voltage-gated sodium channel gene associated with pyrethroid resistance in Latin American Aedes aegypti. Insect Mol Biol 16: 785-798. doi:10.1111/j. 1365-2583.2007.00774.x. PubMed: 18093007.

31. Harris AF, Rajatileka S, Ranson H (2010) Pyrethroid Resistance in Aedes aegypti from Grand Cayman. American. J Trop Med Hyg 83: 277-284. doi:10.4269/ajtmh.2010.09-0623.

32. Chen QX, Zheng WZ, Lin JY, Shi Y, Xie WZ et al. (2000) Effect of metal ions on the activity of green crab (Scylla Serrata) alkaline phosphatase. Int J Biochem Cell Biol 32: 879-885. doi:10.1016/ S1357-2725(00)00026-1. PubMed: 10940645.

33. Van Dyk JS, Pletschke B (2011) Review on the use of enzymes for the detection of organochlorine, organophosphate and carbamate pesticides in the environment. Chemosphere 82: 291-307. doi:10.1016/ j.chemosphere.2010.10.033. PubMed: 21055790.

34. Gamage N, Barnett A, Hempel N, Duggleby RG, Windmill KF et al. (2006) Human sulfotransferases and their role in chemical metabolism. Toxicol Sci 90: 5-22. PubMed: 16322073.

35. Wang LQ, James MO (2006) Inhibition of sulfotransferases by xenobiotics. Curr Drug Metab 7: 83-104. doi: 10.2174/138920006774832596. PubMed: 16454694.

36. Choi J, Rose RL, Hodgson E (2002) In vitro human metabolism of permethrin: the role of human alcohol and aldehyde dehydrogenases. Pestic Biochem Physiol 74: 117-128. doi:10.1016/ S0048-3575(02)00154-2.

37. Hodgson E (2003) In vitro human phase I metabolism of xenobiotics I: Pesticides and related compounds used in agriculture and public health, May 2003. J Biochem Mol Toxicol 17: 201-206. doi:10.1002/jbt. 10080. PubMed: 12898643.

38. David JP, Ismail HM, Chandor-Proust A, Paine MJI (2013) Role of cytochrome P450s in insecticide resistance: impact on the control of mosquito-borne diseases and use of insecticides on Earth. Philos Trans R Soc Lond B-Biol Sci 368: 20120429. PubMed: 23297352.

39. Hemingway J (2000) The molecular basis of two contrasting metabolic mechanisms of insecticide resistance. Insect Biochem Mol Biol 30: 1009-1015. doi:10.1016/S0965-1748(00)00079-5. PubMed: 10989287.