# Phenotypic insecticide resistance in arbovirus mosquito vectors in Catalonia and its capital Barcelona (Spain)

Krijn Paaijmans, Marco Brustollin, Carles Aranda, Roger Eritja, Sandra

Talavera, Nonito Pages, Silvie Huijben

► To cite this version:

Krijn Paaijmans, Marco Brustollin, Carles Aranda, Roger Eritja, Sandra Talavera, et al.. Phenotypic insecticide resistance in arbovirus mosquito vectors in Catalonia and its capital Barcelona (Spain). PLoS ONE, 2019, 14 (7), 13 p. 10.1371/journal.pone.0217860. [https://doi.org/10.1371/journal.pone.0217860](https://doi.org/10.1371/journal.pone.0217860). This is an open access article distributed under a Creative Commons Attribution 4.0 International License.

HAL Id: hal-02620538

https://hal.inrae.fr/hal-02620538

Submitted on 25 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
RESEARCH ARTICLE

# Phenotypic insecticide resistance in arbovirus mosquito vectors in Catalonia and its capital Barcelona (Spain)

Krijn Paaijmans<sup>1,2,3✉</sup>, Marco Brustollin<sup>4,5✉</sup>, Carles Aranda<sup>4,6</sup>, Roger Eritja<sup>6,7</sup>, Sandra Talavera<sup>4</sup>, Nonito Pagès<sup>4,8,9‡</sup>, Silvie Huijben<sup>1,2‡*</sup>

1 ISGlobal, Barcelona, Spain, 2 School of Life Sciences, Center for Evolution and Medicine, Arizona State University, Tempe, AZ, United States of America, 3 The Biodesign Center for Immunotherapy, Vaccines and Virotherapy, Arizona State University, Tempe, AZ, United States of America, 4 Centre de Recerca en Sanitat Animal (CReSA IRTA), Barcelona, Spain, 5 The Center for Infectious Disease Dynamics, and the Huck Institutes of The Life Sciences, Millennium Science Complex, Pennsylvania State University, University Park, PA, United States of America, 6 Servei de Control de Mosquits, Consell Comarcal del Baix Llobregat, Barcelona, Spain, 7 CREAF, Cerdanyola del Vallès, Spain, 8 CIRAD, UMR ASTRE, Petit Bourg, Guadeloupe, France, 9 ASTRE, CIRAD, INRA, Montpellier University, Montpellier, France

✉ These authors contributed equally to this work.

‡ These authors also contributed equally to this work.

* silvie.huijben@asu.edu

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

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

OPEN ACCESS

**Citation:** Paaijmans K, Brustollin M, Aranda C, Eritja R, Talavera S, Pagès N, et al. (2019) Phenotypic insecticide resistance in arbovirus mosquito vectors in Catalonia and its capital Barcelona (Spain). PLoS ONE 14(7): e0217860.
https://doi.org/10.1371/journal.pone.0217860

**Editor:** Vincent CORBEL, Institut de recherche pour le développement, FRANCE

**Received:** November 28, 2018

**Accepted:** May 20, 2019

**Published:** July 5, 2019

**Copyright:** © 2019 Paaijmans 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.

**Data Availability Statement:** All relevant data are within the manuscript and its Supporting Information files.

**Funding:** This work was financially supported by Programa Recercaixa 2011. Silvie Huijben was supported by a Branco Weiss Fellowship. 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.

## Abstract

A range of mosquito species that belong to the Culicidae family are responsible for the worldwide transmission of infectious arboviral diseases such as dengue fever, Zika, West Nile fever and Chikungunya fever. Spain is at risk of arbovirus outbreaks, as various arboviral diseases are frequently introduced and it has established competent vector populations. Autochthonous human cases of West Nile virus have been reported infrequently since 2004, and since October 2018 three autochthonous human case of dengue fever have been confirmed. In response to an outbreak of any arboviral disease, space spraying or fogging will be implemented to control adult mosquito populations. To ensure adulticiding is cost-effective, the insecticide susceptibility status of vectors throughout Catalonia, an autonomous region in north-eastern Spain, was assessed through standardized WHO tube and CDC bottle bioassays. All *Culex pipiens* populations tested were resistant to at least one of the pyrethroids tested, whereas *Aedes albopictus* populations were susceptible to all pyrethroids tested. More detailed studies on the *Cx. pipiens* populations from the Barcelona area (the capital and largest city of Catalonia) revealed resistance to all four classes of public health insecticides available (pyrethroids, carbamates, organophosphates and organochlorides). All *Ae. albopictus* populations were susceptible to those classes, except for one of the tests performed with pirimiphos-methyl (an organophosphate). Pyrethroids are currently the first line chemical class to be used in space spray operations in response to an outbreak of an arboviral disease. While pyrethroids can be effective in reducing *Ae. albopictus* populations, this class may not be effective to control *Cx. pipiens* populations.
Introduction

A range of mosquito species that belong to the Culicidae family are responsible for the world-
wide transmission of infectious arboviral diseases such as dengue fever, Zika, West Nile fever
and Chikungunya fever. On the European continent, autochthonous disease cases (i.e. infec-
tions acquired by mosquito transmission in a currently non-endemic area) have been observed
in several Mediterranean countries. Autochthonous transmission of dengue virus (DENV) has
occurred in France in 2010 [1], Croatia in 2010 [2] and on the Portuguese archipelago Madeira
in 2012 [3], with more than 2000 local cases reported in the latter. Autochthonous West Nile
virus (WNV) cases have been seen in France [4, 5] and Italy in 2011 [6] and the same countries
experienced outbreaks of Chikungunya virus (CHIKV), Italy in 2007 [7] and 2017 [8] and
France in 2010 [9] and 2017 [10].

In Spain, few clinical cases of WNV have been observed in 2004, 2010 and 2016 [11–13]
and serological evidence confirms its presence in Catalonia [14]. Mosquitoes and birds main-
tain the WNV enzootic transmission cycle [15] and periodic cases have been identified in
horses since 2010 [16]. Very recent (October 2018), three cases of autochthonous dengue have
been confirmed in Spain [17], and the virus has been detected in its local vector *Aedes albopic-
tus* [18]. This outbreak was arguably expected given that key ingredients for transmission were
present. First, pathogens such as DENV, Zika virus (ZIKV) and CHIKV are frequently
imported by travelers [19–21]. Second, arboviruses such as DENV, ZIKV, CHIKV and Yellow
Fever virus (YFV) can be transmitted by *Ae. albopictus* (the Asian tiger mosquito), which is
abundant across the Mediterranean coast [22]. In Spain, this species was first detected in the
Barcelona area (Sant Cugat del Vallès) in 2004 [23], but is slowly spreading across the Spanish
mainland [24]. *Ae. aegypti*, another -typically more competent- vector of the same arboviruses,
has recently been detected on the Canary Islands [25]. Other viruses such as WNV and Rift
Valley fever virus (RVFV) can be transmitted by *Ae. albopictus*, *Aedes caspius* and several
*Culex* species, including *Cx. pipiens* s.l. and *Cx. theileri* [26, 27]. And *Cx. pipiens* was also
found to be infected with Usutu virus (USUV) in Catalonia [28]. All these aforementioned
mosquito vectors are abundant in Spain [29, 30], although their role as vector has not been
clearly established for some viruses yet.

As there is no specific prophylactic treatment for most of these arboviruses, disease pre-
vention and control are currently limited to vector control interventions. Interventions in
Spain include (i) adulticiding (space spraying by ground, also known as fogging), (ii) larvi-
ciding, and (iii) larval source reduction via public cooperation and cleaning of landfills [31].
Larviciding is currently used in different mosquito habitats such as flooded areas, agricul-
tural areas and urban sites during mosquito season, usually between early spring and
autumn, with biological insecticides, mainly *Bacillus thuringiensis* var. *israelensis* (Bti) prod-
ucts. Occasionally ground fogging is implemented with different pyrethroids, depending on
mosquito density populations and risk of vector transmission, but adulticide use is infre-
quent. Main regions where mosquito activity is important are the Baix Llobregat in Barce-
lona, the Ebro delta in Tarragona and the Roses Bay in Girona. Also, many municipalities
such as Barcelona city conduct mosquito control operations. In emergency situations, such
as sudden disease outbreaks, the application of adulticides (fogging on the ground with
pyrethroids) is the planned intervention to rapidly target local vector populations [32]. To
ensure this intervention is cost-effective, an insecticide needs to be selected to which the
vectors are susceptible. As such, the level of insecticide susceptibility in a range of compe-
tent vectors needs to be monitored regularly, to allow for quick and effective decision mak-
ing. Here we present data on the insecticide susceptibility in various potential arbovirus
vectors (*Cx. pipiens*, *Ae. albopictus* and *Ae.caspius*) in Catalonia.
**Materials and methods**

This paper combines different insecticide resistance datasets that have been collected in different geographical regions by different research teams for different mosquito species over the past 5 years. The methodology and results sections are presented per mosquito species for clarity.

**Study area**

The potential arbovirus vectors *Cx. pipiens*, *Ae. albopictus* and *Ae.caspius* were collected between 2012 and 2017 in Catalonia, an autonomous region in north-eastern Spain along the Mediterranean coast (Fig 1, drawn using public Global Administrative Areas (GADM) data [33] in the base package of R [34]). The region is about 32,000 km² large, contains approximately 7.5 million inhabitants and has a coastal and inland Mediterranean climate. Subsequent tests focused on the Baix Llobregat region in 2015 to 2017, an area in close proximity to Barcelona, the capital and largest city of Catalonia. The Llobregat delta is one of Catalonia's most important wetland zones, and the park is mainly used for recreation. In addition, there are large agricultural areas and over 800,000 inhabitants living in several towns. The area receives a high number of international visitors every year with the international airport Barcelona-El Prat situated in the Baix Llobregat study area which handled over 44 million passengers in 2016 [35]. Of the 54 mosquito species found in Spain [30], including the invasive *Ae. albopictus* [23], 47 have been detected in Catalonia and 18 in the Baix Llobregat region [36]. Of these, *Cx. pipiens*, *Ae. albopictus* and *Ae. caspius* are the species that have most negative impact in the human population, not only for their ability to transmit human pathogens but also for their feeding habits. Aggressive species like *Ae. albopictus* and *Ae. caspius* can decrease human quality of life in affected communities by reducing the time people can spend outdoors [37] and/or increasing the need for preventive measures (i.e. repellents or nets).

No specific permissions were required for mosquito sampling in these locations, and sampling was carried out under the supervision of public health authorities of the Generalitat of Catalonia. Field studies did not involve endangered or protected species.

**Mosquito collections and rearing**

**Culex pipiens.** Immature stages of *Culex pipiens* were collected from various geographical locations in Catalonia (Fig 1) between 2012 and 2014 and reared to adults in the laboratory on a TetraMin Tropical fish flakes diet at 23°C, 80%RH. Emerging mosquitoes had *ad libitum* access to 10% sucrose solution. To identify the correct taxonomic biotype (*Cx. pipiens pipiens* biotype *pipiens*, *Cx. pipiens pipiens* biotype *molestus* and *hybrid* between *pipiens* and *molestus* biotype) 10 specimens for each locality were selected and analysed by PCR [38]. In 2015 and 2016, immature stages of *Cx. pipiens* were collected from Torrelles de Llobregat (Fig 1) and reared to adult mosquitoes on a diet of Frippak 3CD fish food (Baasrode, Belgium) at 25°C and >60% RH. Emerging mosquitoes had *ad libitum* access to 10% honey-water solution. In 2017, egg rafts of *Cx. pipiens* were collected at el Prat de Llobregat (Fig 1). Eggs were allowed to hatch and the larvae were reared to adult in the laboratory under the same conditions as in 2015 and 2016. Taxonomic biotype was not assessed for the 2015–2017 mosquito specimens.

**Aedes albopictus.** Eggs of *Ae. albopictus* were collected between 2012 and 2014 throughout Catalonia (Fig 1) using oviposition traps (ovitraps) consisting in black plastic 5L bucket half filled with water. A permeable wooden stick made of high-density fibreboard was added as oviposition medium. Eggs were stored under dry conditions and incubated with dechlorinated tap water until L1 larvae emerged. Larvae were maintained on a TetraMin Tropical fish flakes diet at 23°C, 80%RH. In 2016, adult *Ae. albopictus* were collected in Cornellà de
<figure><img src="image_4.png" /><figcaption>Fig 1. Map of Catalonia (Spain) highlighting the areas where *Aedes albopictus* (green symbols; 2012–2014, 2016), *Culex pipiens* (blue symbols; 2012–2014, 2015, 2016, 2017) and *Aedes caspius* (red symbol; 2017) were collected for the insecticide susceptibility assays. Diamonds indicate mosquitoes collected as eggs, triangles show sites with immature collections, and squares indicate adult collections.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0217860.g001

Llobregat (Fig 1) using BG-Sentinel traps (Biogents AG, Regensburg, Germany) equipped with CO₂ (dry ice) and a BG-Lure cartridge. Adult mosquitoes were provided 10% honey water *ad libitum* until the susceptibility test (6 days post-collection) at 25°C and >60% RH. For a second set of susceptibility tests in 2016, *Ae. albopictus* eggs were collected using ovitraps in a rural area in el Prat de Llobregat (Fig 1). Eggs were stored under dry conditions and incubated with dechlorinated tap water until L1 larvae emerged. Emerged larvae were reared to adult mosquitoes on a diet of Frippak 3CD fish food at 25°C and >60% RH.
**Aedes caspius.** A period of heavy rainfall allowed us to collect fourth instar larvae and pupae of the floodwater mosquito *Ae. caspius* at La Ricarda, a protected wetland in el Prat de Llobregat (Fig 1) in 2017, which were reared to adults on a diet of Frippak 3CD fish food at 25°C and >60% RH.

## Insecticide susceptibility assays

Insecticide susceptibility was tested using two different methods: CDC bottle bioassays and WHO tube tests.

**CDC bottle bioassays.** CDC bottle assays were performed following the guidelines of the Centers for Disease Control and Prevention [39]. Analytical standard lambda-cyhalothrin, delta-methrin, permethrin, pirimiphos-methyl, 4,4'-DDT and bendiocarb (Pestanal Analytical standard, Sigma-Aldrich) were dissolved in pure acetone to diagnostic concentrations (Table 1). All bottles (treatment and controls) were coated with insecticides on the same day of each test. To assess insecticide susceptibility levels, up to 25 two-to-six day old non-blood-fed female mosquitoes were introduced into each bottle, with four insecticide-coated and one control bottle in total (tests with *Ae. albopictus* from Figueras and Llorenç del Penedes contained half the numbers).

**WHO tube test.** WHO tube tests were performed according to the guidelines published by the WHO [40] using WHO-standard insecticide-treated papers (Table 1). For each insecticide tested, approximately 100 mosquitoes (4 replicates of 25) for the exposed group and 50 mosquitoes (2 replicates of 25) for the control group were kept in a holding tube lined with untreated paper. Next, they were exposed for 1 hour to insecticide-treated or control paper,

<table><caption>Table 1. Available information on the discriminating bioassay concentrations of insecticides for determining susceptibility of adult anopheline and culicine mosquitoes using WHO and CDC insecticide susceptibility tests for the insecticides used in our studies.</caption><thead><tr><th rowspan="2"></th><th colspan="3">Pyrethoids</th><th colspan="2">Carbamates</th><th rowspan="2">Organophosphates<br>Pirimiphos-methyl</th><th rowspan="2">Organochlorines<br>DDT</th></tr><tr><th>Deltamethrin</th><th>Lambda-cyhalothrin</th><th>Permethrin</th><th>Bendiocarb</th><th>Propoxur</th></tr></thead><tbody><tr><td colspan="8"><strong>WHO tube assay (%)</strong></td></tr><tr><td>Recommended dose <em>Anopheles</em> spp<sup>a</sup></td><td>0.05</td><td>0.05</td><td>0.75</td><td>0.1</td><td>0.1</td><td>0.25</td><td>4</td></tr><tr><td>Recommended dose <em>Aedes</em> spp<sup>a</sup></td><td>0.03<sup>b</sup></td><td>0.03</td><td>0.25</td><td>-</td><td>-</td><td>0.21<sup>c</sup></td><td>-</td></tr><tr><td><em>Cx. pipiens</em> tests 2015<sup>a</sup></td><td>0.05</td><td>-</td><td>-</td><td>-</td><td>0.1</td><td>0.25</td><td>4</td></tr><tr><td colspan="8"><strong>CDC bottle bio-assay (µg/bottle)</strong></td></tr><tr><td>Recommended dose <em>Anopheles</em> spp<sup>a,d</sup></td><td>12.5</td><td>12.5</td><td>21.5</td><td>12.5</td><td>-</td><td>20</td><td>100</td></tr><tr><td>Recommended dose <em>Aedes</em> spp<sup>a,d</sup></td><td>10</td><td>10</td><td>15</td><td>12.5</td><td>-</td><td>-</td><td>75</td></tr><tr><td><em>Ae. albopictus</em> tests 2012–2013<sup>e</sup></td><td>1000</td><td>1000</td><td>3125</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td><em>Ae. albopictus</em> tests 2016</td><td>10</td><td>-</td><td>-</td><td>10<sup>f</sup></td><td>-</td><td>20</td><td>75</td></tr><tr><td><em>Ae. caspius</em> tests 2017</td><td>10</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>75</td></tr><tr><td>Recommended dose <em>Culex</em> spp</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td><em>Cx. pipiens</em> tests 2012–2013<sup>e</sup></td><td>250</td><td>250</td><td>2500</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td><em>Cx. pipiens</em> tests 2016</td><td>10<sup>g</sup>/25<sup>h</sup></td><td>-</td><td>-</td><td>10<sup>f</sup></td><td>-</td><td>20<sup>g</sup>/40<sup>h</sup></td><td>75</td></tr><tr><td><em>Cx. pipiens</em> tests 2017</td><td>25<sup>h</sup></td><td>-</td><td>-</td><td>25<sup>h</sup></td><td>-</td><td>40<sup>h</sup></td><td>200<sup>h</sup></td></tr></tbody></table>

<footnote><sup>a</sup> WHO guidelines [40]</footnote>

<footnote><sup>b</sup> tentative [40],</footnote>

<footnote><sup>c</sup> determined for *Anopheles*, tentative for *Aedes* [40] (we assume a typo, as the value for *Anopheles* is 0.25),</footnote>

<footnote><sup>d</sup> CDC guidelines [39],</footnote>

<footnote><sup>e</sup> WHO guidelines for *Culex quinquefasciatus* [41], converted to µg/bottle</footnote>

<footnote><sup>f</sup> lower dose by error</footnote>

<footnote><sup>g</sup> value for *Aedes* used</footnote>

<footnote><sup>h</sup> twice the *Anopheles* dose, as recommended by Dr. William Brogdon (CDC, Atlanta, personal communication)</footnote>
respectively, using WHO-standard exposure tubes and procedures. After exposure, mosqui-
toes were transferred to paper holding cups and had *ad libitum* access to sugar water (10%
household sugar solution).

**Insecticide exposures.** *Cx pipiens* mosquitoes collected from various areas within Catalonia in 2012–2014 were tested using the CDC bottle assay with lambda-cyhalothrin, deltamethrin and permethrin. *Cx pipiens* mosquitoes from El Prat de Llobregat in 2015 were exposed to deltamethrin, propoxur, pirimiphos-methyl and 4,4'-DDT using WHO tube test. The collected *Cx pipiens* mosquitoes from 2016–2017 were exposed to deltamethrin, pirimiphos-methyl, 4,4'-DDT and bendiocarb in the CDC bioassay (Table 1).

*Ae. albopictus* mosquitoes collected in 2012/2014 were exposed to lambda-cyhalothrin, deltamethrin and permethrin, and in 2016, these mosquitoes were tested with deltamethrin, pirimiphos-methyl, 4,4'-DDT and bendiocarb. All tests on *Ae. albopictus* mosquitoes were performed with CDC bottle assays. *Ae. caspius* mosquitoes from 2017 were exposed to DDT and deltamethrin in the CDC bottle assay.

**Insecticide concentrations.** Of note is that the diagnostic dosages for insecticide resistance are not clear for each mosquito species or indeed some genera. As there were no established diagnostic doses for *Culex* species in general, but see [41] for guidance on individual mosquito species for some insecticides, two sets of susceptibility tests were conducted in 2016. In the first test, *Cx. pipiens* were exposed to diagnostic doses of pirimiphos-methyl, DDT and deltamethrin concentrations as determined for *Aedes* species. In the second test, mosquitoes from the same population were exposed to twice the diagnostic dose for *Anopheles* species (personal communication Dr. William Brogdon, CDC, Atlanta). For bendiocarb 10µg/bottle was accidentally used opposed to the recommended diagnostic dose of 12.5µg/bottle. The different recommended concentrations known for each mosquito genera and the insecticide concentrations used in this study are summarized in Table 1.

**Data analysis**

For the WHO tube bioassay mortality was assessed twenty-four hours post-exposure in each
replicate. In the CDC bottle bioassay, susceptibility was determined by quantifying mortality
at a diagnostic time in each bottle. The diagnostic time was 30 minutes for for lambda-cyhalo-
thrin, deltamethrin, permethrin, pirimiphos-methyl, and bendiocarb and 45 minutes for
DDT. A mosquito was identified as 'dead' if she lacked the ability to stand on her feet at a gen-
tle rotation of the bottle or tube. When control mortality was between 5 and 20%, mortality
was corrected using Abbott's formula [40], with negative values being rounded to zero (note
that control mortality is shown in S1 Table for *Cx. pipiens*, S2 Table for *Ae. albopictus* and S3
Table for *Ae. caspius*).

Resistance to the insecticide tested was identified as per the current WHO guidelines [40]. Mosquito populations were classified as susceptible if mortality ranges from 98 to 100% 24h after insecticide exposure. Mosquito mortality between 90 and 97% was defined as suspected resistance but would require follow-up testing to confirm resistance (if resistance is again below 98%). If mortality is less than 90%, resistance is confirmed and -if at least 100 mosqui-toes are tested- no additional bioassays are needed.

**Results**

Culex pipiens

Multiple biotypes have been detected in population from the following localities: Bellaterra
(biotypes *pipiens* and *molestus*); Bellvis (biotypes *pipiens* and *molestus*); Gava' (biotypes *pipiens*
and *hybrid*); Santa Coloma de Cervelló (biotypes *pipiens* and *hybrid*), Vic (biotypes *pipiens*,
molestus and *hybrid*). The population from Empuriabrava was characterized by the presence of a unique biotype: *Cx. pipiens pipiens* biotype *molestus*.

Resistance in *Cx. pipiens* to deltamethrin was observed in all six locations studied in Catalonia during 2012–2014 (Table 2). Resistance to lambda-cyhalothrin was seen in two of the six areas and suspected in the other three areas. Resistance to permethrin was detected in one of the six areas and suspected in two other areas. Mortality in response to deltamethrin exposure was particularly low in Bellvis (4%).

Additional testing with *Cx. Pipiens* populations from the Barcelona area in 2015 showed this species was resistant to all four public health insecticides: deltamethrin (81% mortality), DDT (6%), pirimiphos-methyl (0%) and propoxur (0%) (Table 2). Subsequent tests (CDC bottle bioassay) confirmed resistance to the three chemical classes tested (pyrethroid deltamethrin: 96% [2016] and 96% [2017] mortality; organophospate pirimiphos-methyl: 0% [2016] and 33% [2017]; organochloride DDT: 78% [2016] and 68% [2017]). Deltamethrin and pirimiphos-methyl exposures were repeated in 2016 with a higher diagnostic dose (twice the dose for *Anopheles*, Table 1), which showed clear resistance to pirimiphos-methyl (0% mortality) and suspected resistance to deltamethrin (97.5%) (Table 2). Resistance to bendiocarb (46% mortality) was observed in 2017.

## Aedes albopictus

No resistance to pyrethroids was detected in the *Ae. albopictus* populations tested throughout Catalonia during 2012–2014 (Table 3). Additional tests in 2016 with mosquitoes from the Barcelona area showed that *Ae. albopictus* populations was fully susceptible to deltamethrin, bendiocarb and DDT, but showed levels of resistance to pirimiphos-methyl (86.52% mortality) in one out of the two tests performed.

<table><caption>**Table 2.** Insecticide susceptibility of *Cx. pipiens* collected from various localities in Catalunya (Spain). Percentage indicates percent mortality (WHO tube tests: 24h following 1h exposure; CDC bottle assays: At discriminating exposure time of 30 min (45 min for DDT)); number between parentheses indicates the number of mosquitoes tested.</caption><thead><tr><th rowspan="3"></th><th colspan="6">Pyrethroids</th><th colspan="2">Carbamates</th><th colspan="2">Organochloride</th><th colspan="2">Organophosphate</th></tr><tr><th colspan="2">Permethrin</th><th colspan="2">Deltamethrin</th><th colspan="2">Lamda-cyhalothrin</th><th colspan="2">Bendiocarb</th><th colspan="2">Propoxur</th><th colspan="2">DDT</th><th colspan="2">Pirimiphos-methyl</th></tr><tr><th>CDC</th><th>WHO</th><th>CDC</th><th>WHO</th><th>CDC</th><th>WHO</th><th>CDC</th><th>WHO</th><th>WHO</th><th>CDC</th><th>WHO</th><th>WHO</th><th>CDC</th></tr></thead><tbody><tr><td colspan="14"><strong>2012–2014</strong></td></tr><tr><td>Bellaterra</td><td>97% (100)</td><td>-</td><td>80% (100)</td><td>100% (100)</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Bellvis</td><td>81% (100)</td><td>-</td><td>4% (100)</td><td>55% (100)</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Empuriabrava</td><td>95% (100)</td><td>-</td><td>62% (100)</td><td>97% (100)</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Gavà</td><td>99% (100)</td><td>-</td><td>35% (100)</td><td>85% (100)</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Santa Coloma de Cervelló</td><td>98% (100)</td><td>-</td><td>63% (100)</td><td>93% (100)</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Vic</td><td>99% (100)</td><td>-</td><td>87% (100)</td><td>95% (100)</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td colspan="14"><strong>2015</strong></td></tr><tr><td>Torrelles de Llobregat</td><td>-</td><td>81.2% (101)</td><td>-</td><td>-</td><td>-</td><td>-</td><td>0% (91)</td><td>6.3% (96)</td><td>-</td><td>0% (98)</td><td>-</td><td>-</td><td>-</td></tr><tr><td colspan="14"><strong>2016</strong></td></tr><tr><td>Torrelles de Llobregat</td><td>-</td><td>-</td><td>96.1% (85)<sup>1</sup><br>97.5% (91)<sup>2</sup></td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>78.2% (58)</td><td>-</td><td>0% (71)<sup>3</sup><br>0% (49)<sup>4</sup></td><td>-</td></tr><tr><td colspan="14"><strong>2017</strong></td></tr><tr><td>El Prat de Llobregat</td><td>-</td><td>-</td><td>96.2% (106)<sup>2</sup></td><td>-</td><td>46.2% (104)</td><td>-</td><td>-</td><td>-</td><td>68.3% (101)</td><td>-</td><td>33.3% (90)</td><td>-</td></tr></tbody></table>

<footnote n="1">¹10 µg/bottle;</footnote>

<footnote n="2">²25 µg/bottle;</footnote>

<footnote n="3">³20 µg/bottle;</footnote>

<footnote n="4">⁴40 µg/bottle</footnote>
<table><caption>**Table 3.** Insecticide susceptibility of *Ae. albopictus* collected from various localities in Catalunya (Spain). Percentage indicates percent mortality at discriminating exposure time of 30 min (45 min for DDT); number between parentheses indicates the number of mosquitoes tested.</caption><thead><tr><th rowspan="3"></th><th colspan="3">Pyrethroids</th><th>Carbamates</th><th>Organochloride</th><th>Organophosphate</th></tr><tr><th>Permethrin</th><th>Deltamethrin</th><th>Lamda-cyhalthrin</th><th>Bendiocarb</th><th>DDT</th><th>Pirimiphos-methyl</th></tr><tr><th>CDC</th><th>CDC</th><th>CDC</th><th>CDC</th><th>CDC</th><th>CDC</th></tr></thead><tbody><tr><td colspan="7"><strong>2012-2014</strong></td></tr><tr><td>Constantí</td><td>100% (100)</td><td>99% (100)</td><td>100% (100)</td><td>-</td><td>-</td><td>-</td></tr><tr><td>El Prat</td><td>100% (100)</td><td>99% (100)</td><td>100% (100)</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Figueras</td><td>100% (50)</td><td>100% (50)</td><td>100% (50)</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Lloret de Mar</td><td>100% (100)</td><td>100% (100)</td><td>100% (100)</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Llorenç del Penedes</td><td>100% (50)</td><td>100% (50)</td><td>100% (50)</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Mont-Roig del Camp</td><td>100% (100)</td><td>98% (100)</td><td>100% (100)</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Sant Julià de Ramis</td><td>100% (100)</td><td>100% (100)</td><td>100% (100)</td><td>-</td><td>-</td><td>-</td></tr><tr><td colspan="7"><strong>2016</strong></td></tr><tr><td>Cornella de Llobregat</td><td>-</td><td>100% (77)</td><td>-</td><td>100% (88)</td><td>100% (86)</td><td>86.5% (88)</td></tr><tr><td>El Prat de Llobregat</td><td>-</td><td>100% (82)</td><td>-</td><td>-</td><td>-</td><td>98.5% (86)</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0217860.t003

## Aedes caspius

*Ae. caspius* (collected from the Barcelona area in 2017) was fully susceptible to DDT and delta-methrin (both 100% mortality, other chemical classes not tested) (Table 4).

## Discussion

Here we provide evidence that several *Cx. pipiens* populations from Catalonia and Barcelona, which are potential vectors of arboviruses like WNV, USUV and RVFV [26, 27, 42] are resistant to some of the pyrethroid insecticides evaluated. We observed spatial variation in resistance phenotype, with high levels of resistance in Bellvis and to a lesser extent in Gava. One explanation for this observation could be the extensive use of pyrethroids in the agricultural activity in the past and at present. In contrast, no large and repetitive adult and larval control have been applied in the area with these products nor in the past nor at present. Since the 1980's, mosquito control in the area has been conducted with larvicides, formerly Temephos and since 1992 with Bti. No pyrethroids have been used as larvicides. Few adulticides have been conducted with Malathion, except for occasional adulticiding with pyrethroids. The only possible continuous application of adulticides using this family of insecticides originates from domestic insecticide devices that are largely used by public in general. Moreover, the *Cx. pipiens* populations from the Barcelona area are resistant to insecticides belonging to all four classes of public health insecticides available: Resistance is clearly demonstrated for DDT (organochloride), pirimiphos-methyl (organophosphate), propoxur (carbamate) and various pyrethroids. Again, one explanation is the widely use of DDT in agriculture and public health in the past in all Spain until 1977 when it was banned and the possible illegal use has been certainly almost inexistent. On the other hand, insecticide resistant genes may have migrated

<table><caption>**Table 4.** Insecticide susceptibility of *Ae. caspius* collected from Baix Llobregat (Barcelona, Spain) in 2017. Percentage indicates percent mortality at discriminating exposure time of 30 min (45 min for DDT); number between parentheses indicates the number of mosquitoes tested.</caption><thead><tr><th rowspan="2"></th><th>Pyrethroids</th><th>Organochloride</th></tr><tr><th>Deltamethrin</th><th>DDT</th></tr></thead><tbody><tr><td>La Ricarda</td><td>100% (85)</td><td>100% (43)</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0217860.t004
from other parts of Europe and/or the world, as extensive and long-distance migration of
genes has been observed in *Cx. pipiens* [43].

In contrast, *Ae. albopictus*, a potential vector of DENV, ZIKV, CHIKV and USUV [42] was
susceptible to all four classes of public health insecticides across different years and methodol-
ogies used. However, a single observation suggested some level of resistance to the organo-
phosphate pirimiphos-methyl, which needs to be confirmed. These findings are similar to
those observed in other European *Ae. albopictus* populations, such as in Switzerland [44].
However, pyrethroid resistance seems to be emerging in European *Ae. albopictus* populations,
as observed in Italy [45] and in Spain, with possible resistance to cypermethrin in Barcelona
[46], and the kdr mutation (V1016G) has been recently reported in Italy [47].

*Ae. albopictus* was first detected in Spain in 2004 close to the Llobregat study area [23]. Although it appears to be spreading across Spain, recent genetic analysis revealed that present populations of *Ae. albopictus* are likely linked to multiple independent introduction events into Spain [48]. The absence of observed resistance to all four classes of insecticides in the populations tested suggests that these incoming populations may not have experienced significant insecticide pressure before being introduced in Spain (either through vector control interventions, or via agricultural applications).

*Ae. caspius* was fully susceptible to deltamethrin, similar to observations in Greece [49], and to DDT. Its observed susceptibility is likely linked with its ecological niche: it is constrained to swampy natural areas, environments that are not human-made and where typically few pesticides are applied. This contrasts with *Cx. pipiens* habitats that are often man-made and have been shown to contain large insecticide quantities in the past [50].

It is important to keep in mind that obtaining meaningful results for some chemical classes, particular insecticides and mosquito species is difficult. There is no agreement yet on the discriminating concentrations for *Aedes* in the WHO tube bioassays [51] and data are missing for several chemical classes and insecticides. The situation for *Culex* vectors is worse: there are no guidelines on diagnostic concentrations apart from some scarce information on some species and few insecticides [41]. Regarding the CDC bottle bioassay [39], the diagnostic doses and times have been established for most insecticides, but not for pirimiphos-methyl (*Aedes*), and there is again no information for *Culex* species. As such, the observed resistance in *Cx. pipiens* may not be 'true' resistance, but simply an artefact of the doses and/or diagnostic time used. Preferably, dose-response curves are generated to identify the discriminating dose for each species and each insecticide, after which longitudinal susceptibility studies starting prior to insecticide use are implemented to observe changes over time. However, this is rarely practiced and indeed not the case for our study area.

In addition, it is key to consider that the outcomes presented pertain to standardized laboratory bio-assays, following existing guidelines and using suggested diagnostic dosages. Susceptibility in such conditions do not necessarily reflect susceptibility in programmatically applied interventions in field situations. Several other factors, not considered in these standard susceptibility tests, have been shown to impact insecticide toxicity, such as exposure temperature [52], variations in the larval environment [53], mosquito age [54], (multiple) blood feeding events [55] and infections with pathogens [56]. Moreover, actual doses and contact times with the insecticide following space spraying may vary in the field and differ from those tested in the lab. This will likely impact the efficacy of fogging campaigns and -as a consequence- disease transmission potential. Furthermore, mosquitoes that do survive an insecticide exposure may still die within a few days, reducing their transmission potential compared to unexposed mosquitoes, a phenomenon referred to as delayed mortality [57]. In other words, the classification 'resistant' does not necessarily imply an insecticide will not be effective in the field. This
ties in with the much larger debate on the impact (or not) of insecticide resistance on disease transmission [58–60].

## Conclusions

At the moment insecticides remain one of the corner stones of arboviral disease outbreak control and prevention. The results presented here indicate that *Ae. albopictus* populations across Catalonia can be controlled by space spraying with pyrethroids and maybe -at least in the Barcelona area- with organophosphates, the two chemical classes approved for space spraying by WHO [61]. For residual spray applications (of e.g. mosquito resting areas inside homes or cattle stables/sheds) the carbamate class can be used (again only confirmed for Barcelona). Having said that, the number of chemical classes approved for adulticiding by the European Union is limited, with pyrethroids being exclusively recommended in local and national guidelines [62]. Unfortunately, the situation for *Cx. pipiens* appears worse (but note the discussion on appropriate diagnostic doses/times) as resistance to all four classes of public health insecticides has been observed. As we rely heavily on a very limited choice of chemicals (pyrethroids are currently the first line chemical class for adulticiding in response to outbreaks of any arboviral disease), proper resistance management strategies should be in place [51] for any vector species that is targeted. This necessitates continuous entomological surveillance to monitor vector population dynamics and their insecticide susceptibility status [63].

## Supporting information

S1 Table. Control mortality during insecticide susceptibility testing of *Cx. pipiens* collected from various localities in Catalunya (Spain).

(DOCX)

S2 Table. Control mortality during insecticide susceptibility testing of *Ae. albopictus* collected from various localities in Catalunya (Spain).

(DOCX)

S3 Table. Control mortality during insecticide susceptibility testing of *Ae. caspius* collected from Baix Llobregat (Barcelona, Spain) in 2017.

(DOCX)

## Acknowledgments

We are grateful to the SCM staff for their assistance during mosquito collections. ISGlobal is a member of the CERCA Programme, Generalitat de Catalunya.

## Author Contributions

**Conceptualization:** Krijn Paaijmans, Marco Brustollin, Carles Aranda, Roger Eritja, Sandra Talavera, Nonito Pagès, Silvie Huijben.

**Data curation:** Krijn Paaijmans, Marco Brustollin, Silvie Huijben.

**Formal analysis:** Krijn Paaijmans, Marco Brustollin, Silvie Huijben.

**Investigation:** Krijn Paaijmans, Marco Brustollin, Silvie Huijben.

**Writing – original draft:** Krijn Paaijmans, Silvie Huijben.

**Writing – review & editing:** Marco Brustollin, Carles Aranda, Roger Eritja, Sandra Talavera, Nonito Pagès.
## References

1. La Ruche G, Souarès Y, Armengaud A, Peloux-Petiot F, Delaunay P, Desprès P, et al. First two autochthonous dengue virus infections in metropolitan France. Euro Surveill. 2010; 15(39):pii = 19676.

2. Gjenero-Margan I, Aleraj B, Krajcar D, Lesnikar V, Klobucar A, Pem-Novosel I, et al. Autochthonous dengue fever in Croatia, August-September 2010. Euro Surveill. 2011; 16(9):pii = 19805.

3. ECDC. Communicable Disease Threats Report: Week 1, 30 December-5 January 2013. European Centre for Disease Prevention and Control, 2013.

4. Durand JP, Simon F, Tolou H. West Nile virus: in France again, in humans and horses. Rev Prat. 2004; 54(7):703-10. PMID: 15253284

5. Bahuon C, Marcillaud-Pitel C, Bournez L, Leblond A, Beck C, Hars J, et al. West Nile virus epizootics in the Camargue (France) in 2015 and reinforcement of surveillance and control networks. Rev Sci Tech. 2016; 35(3):811-24. https://doi.org/10.20506/rst.35.3.2571 PMID: 28332648

6. Bagnarelli P, Marinelli K, Trotta D, Monachetti A, Tavio M, Del Gobbo R, et al. Human case of autochthonous West Nile virus lineage 2 infection in Italy, September 2011. Euro Surveill. 2011; 16(43):pii = 20002.

7. Angelini R, Finarelli A, Angelini P, Po C, Petropulacos K, Silvi G, et al. Chikungunya in north-eastern Italy: a summing up of the outbreak. Euro Surveill. 2007; 12(47):pii = 3313.

8. Venturi G, Di Luca M, Fortuna C, Remoli ME, Riccardo F, Severini F, et al. Detection of a chikungunya outbreak in Central Italy, August to September 2017. Euro Surveill. 2017; 22(39):pii = 17-00646.

9. Grandadam M, Caro V, Plumet S, Thiberge J-M, Souarès Y, Failloux A-B, et al. Chikungunya Virus, Southeastern France. Emerg Infect Dis. 2011; 17(5):910-3. https://doi.org/10.3201/eid1705.101873 PMID: 21529410

10. Calba C, Guerbois-Galla M, Franke F, Jeannin C, Auzet-Caillaud M, Grard G, et al. Preliminary report of an autochthonous chikungunya outbreak in France, July to September 2017. Euro Surveill. 2017; 22 (39):pii = 17-00647.

11. García-Bocanegra I, Jaén-Téllez JA, Napp S, Arenas-Montes A, Fernández-Morente M, Fernández-Molera V, et al. West Nile fever outbreak in horses and humans, Spain, 2010. Emerg Infect Dis. 2011; 17(12):2397. https://doi.org/10.3201/eid1712.110651 PMID: 22172565

12. Kaptoul D, Viladrich PF, Domingo C, Niubo J, Martinez-Yelamos S, De Ory F, et al. West Nile virus in Spain: report of the first diagnosed case (in Spain) in a human with aseptic meningitis. Scand J Infect Dis. 2007; 39(1):70-1. https://doi.org/10.1080/00365540600740553 PMID: 17366016

13. Centro de Coordinación de Alertas y Emergencias Sanitarias. Evaluación de riesgo del virus del Nilo Occidental. Caso humano y focos en equinos de la CA de Andalucía 2016 [cited 2018 18 March]. Available from: https://www.msssi.gob.es/profesionales/saludPublica/ccayes/alertasActual/docs/ERR_VNO_2016_31_08_DEF.pdf.

14. Bofill D, Domingo C, Cardeñosa N, Zaragoza J, de Ory F, Minguell S, et al. Human West Nile Virus Infection, Catalonia, Spain. Emerg Infect Dis. 2006; 12(7):1163-4. https://doi.org/10.3201/eid1207.060164 PMID: 16845777

15. Ferraguti M, JM LAP, Soriguer R, Llorente F, Jimenez-Clavero MA, Figuerola J. West Nile virus-neutralizing antibodies in wild birds from southern Spain. Epidemiol Infect. 2016; 144(9):1907-11. https://doi.org/10.1017/S0950268816000133 PMID: 26846720

16. Junta de Andalucia. Fiebre del Nilo Occidental (West Nile virus) 2017 [cited 2018 18 March]. Available from: http://www.juntadeandalucia.es/organismos/agriculturapescaydesarrollorural/areas/ganaderia/sanidad-animal/paginas/fiebre-nilo-occidental.html.

17. ECDC. Local transmission of dengue fever in France and Spain-2018 (Rapid risk assessment, 22 October 2018). 2018.

18. Aranda C, Martínez MJ, Montalvo T, Eritja R, Navero-Castillejos J, Herreros E, et al. Arbovirus surveillance: first dengue virus detection in local Aedes albopictus mosquitoes in Europe, Catalonia, Spain, 2015. 2018; 23(47):1700837.

19. Requena-Méndez A, García C, Aldasoro E, Vicente J, Martínez M, Pérez-Molina J, et al. Cases of chikungunya virus infection in travellers returning to Spain from Haiti or Dominican Republic, April-June 2014 Euro Surveill. 2014; 19(28):pii = 20853.

20. Neumayr A, Muñoz J, Schunk M, Bottieau E, Cramer J, Calleri G, et al. Sentinel surveillance of imported dengue via travellers to Europe 2012 to 2014: TropNet data from the DengueTools Research Initiative. Euro Surveill. 2017; 22(1):pii = 30433.

21. Sanchez BS, Tato LM, Martin SG, Perez E, Grasa C, Valderrama S, et al. Imported malaria in children in Madrid, Spain, 2007-2013. Enferm Infecc Microbiol Clin. 2017; 23(16):30395-0.
22. Collantes F, Delacour S, Alarcón-Elbal PM, Ruiz-Arrondo I, Delgado JA, Torrell-Sorio A, et al. Review of ten-years presence of Aedes albopictus in Spain 2004–2014: known distribution and public health concerns. Parasit Vectors. 2015; 8(1):655.

23. Aranda C, Eritja R, Roiz D. First record and establishment of the mosquito *Aedes albopictus* in Spain. Med Vet Entomol. 2006; 20:150–2. https://doi.org/10.1111/j.1365-2915.2006.00605.x PMID: 16608499

24. Alarcón-Elbal PM, Delacour-Estrella S, Ruiz-Arrondo I, Collantes F, Delgado JA, Morales-Bueno J. Updated distribution of Aedes albopictus (Diptera: Culicidae) in Spain: new findings in the mainland Spanish Levante, 2013. Mem Inst Oswaldo Cruz. 2014; 109(6):782–6. https://doi.org/10.1590/0074-0276140214 PMID: 25317706

25. Gobierno de Canaries. El Sistema de Vigilancia Entomológica de Canarias detecta la presencia de mosquito Aedes aegypti en Fuerteventura 2017 [cited 2018 18 March]. Available from: http://www.gobcan.es/noticias/lanzadera/90127/sistema-vigilancia-entomologica-canarias-detecta-presencia-mosquito-aedes-aegypti-fuerteventura.

26. Brustolin M, Talavera S, NuÑEz A, Santamaría C, Rivas R, Pujol N, et al. Rift Valley fever virus and European mosquitoes: vector competence of Culex pipiens and Stegomyia albopicta (= Aedes albopictus). Med Vet Entomol. 2017; 31(4):365–72. https://doi.org/10.1111/mve.12254 PMID: 28782121

27. Brustolin M, Talavera S, Santamaría C, Rivas R, Pujol N, Aranda C, et al. *Culex pipiens* and *Stegomyia albopicta* (= *Aedes albopictus*) populations as vectors for lineage 1 and 2 West Nile virus in Europe. Med Vet Entomol. 2016; 30(2):166–73. https://doi.org/10.1111/mve.12164 PMID: 26890285

28. Busquets N, Alba A, Allepuz A, Aranda C, Ignacio Nunez J. Usutu virus sequences in *Culex pipiens* (Diptera: Culicidae), Spain. Emerg Infect Dis. 2008; 14(5):861–3. https://doi.org/10.3201/eid1405.071577 PMID: 18439389

29. Bueno-Marí R, Bernués-Bañeres A, Jiménez-Peydró R. Updated checklist and distribution maps of mosquitoes (Diptera: Culicidae) of Spain. Eur Mosq Bull. 2012; 30:91–126.

30. Eritja R, Aranda C, Padrós J, Goula M, Lucientes J, Escosa R, et al. An annotated checklist and bibliography of the mosquitoes of Spain (Diptera: Culicidae). Eur Mosq Bull. 2000; 8:10–8.

31. Abramides GC, Roiz D, Guitart R, Quintana S, Giménez N. Control of the Asian tiger mosquito (Aedes albopictus) in a firmly established area in Spain: risk factors and people's involvement. T Roy Soc Trop Med H. 2013;107.

32. Catalunya AdSPd. Protocol per a la vigilància i el control de les arbovirosis transmeses per mosquits a Catalunya (actualitzat 28/06/2017). Barcelona: Generalitat de Catalunya. Departament de Salut, 2017.

33. University of California Berkeley. Global Administrative Areas (retreived 07/02/2018; http://www.gadm.org). 2018.

34. R Development Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing. Retrieved from: http://www.R-project.org. 2018.

35. Observatori del Turisme a Barcelona: ciutat i regió. Barcelona tourism activity report 2016. Barcelona: Observatori del Turisme a Barcelona: ciutat i regió, 2017.

36. Padrós J. Culícids del Baix Llobregat (Diptera:Culicidae) faunística i biologia. Barcelona: University of Barcelona; 1997.

37. Romi R. *Aedes albopictus* in Italia: un problema sanitario sottovalutato. Ann Ist Super Sanità. 2001; 37 (2):241-7. PMID: 11758282

38. Bahnck CM, Fonseca DM. Rapid assay to identify the two genetic forms of *Culex* (*Culex*) *pipiens* L. (Diptera: Culicidae) and hybrid populations. Am J Trop Med Hyg. 2006; 75(2):251–5. PMID: 16896127

39. Brogdon W, Chan A. Guideline for evaluating insecticide resistance in vectors using the CDC bottle bioassay. Atlanta: Centers for Disease Control and Prevention, 2010.

40. WHO. Test procedures for insecticide resistance monitoring in malaria vector mosquitoes. Second edition. Geneva: World Health Organization; 2016.

41. WHO. Vector resistance to pesticides: fifteenth report of the WHO expert committee on vector biology and control. Geneva: World Health Organization; 1992.

42. Calzolari M, Gaibani P, Bellini R, Defilippo F, Pierro A, Albieri A, et al. Mosquito, Bird and Human Surveillance of West Nile and Usutu Viruses in Emilia-Romagna Region (Italy) in 2010. PLoS ONE. 2012; 7 (5):e38058. https://doi.org/10.1371/journal.pone.0038058 PMID: 22666446

43. Chevillon C, Raymond M, Guillemaud T, Lenormand T, Pasteur N. Population genetics of insecticide resistance in the mosquito *Culex pipiens*. Biol J Linn Soc. 1999; 68(1):147–57.

44. Suter T, Crespo MM, de Oliveira MF, de Oliveira TSA, de Melo-Santos MAV, de Oliveira CMF, et al. Insecticide susceptibility of *Aedes albopictus* and *Ae. aegypti* from Brazil and the Swiss-Italian border region. Parasit Vectors. 2017; 10(1):017–2364.
45. Pichler V, Bellini R, Veronesi R, Arnoldi D, Rizzoli A, Lia RP, et al. First evidence of resistance to pyre-throid insecticides in Italian *Aedes albopictus* populations 26 years after invasion. Pest Manag Sci. 2018; 74(6):1319-27. https://doi.org/10.1002/ps.4840 PMID: 29278457

46. Bengoa M, Eritja R, Delacour S, Miranda MA, Sureda A, Lucientes J. First data on resistance to pyre-throids in wild populations of *Aedes albopictus* from Spain. J Am Mosq Control Assoc. 2017; 33(3):246-9. https://doi.org/10.2987/17-6636R.1 PMID: 28854116

47. Pichler V, Malandruccolo C, Serini P, Bellini R, Severini F, Toma L, et al. Phenotypic and genotypic pyrethroid resistance of *Aedes albopictus*, with focus on the 2017 chikungunya outbreak in Italy. Pest Manag Sci. in press.

48. Brustolin M. Autochthonous and invasive mosquitoes of Catalonia as vectors of zoonotic arboviruses [PhD-thesis]. Bellaterra: Universitat Autònoma de Barcelona; 2016.

49. Fotakis EA, Chaskopoulou A, Grigoraki L, Tsiamantas A, Kounadi S, Georgiou L, et al. Analysis of pop-ulation structure and insecticide resistance in mosquitoes of the genus Culex, Anopheles and Aedes from different environments of Greece with a history of mosquito borne disease transmission. Acta Trop. 2017; 174:29-37. https://doi.org/10.1016/j.actatropica.2017.06.005 PMID: 28606820

50. Eritja R. Anàlisi integrada sobre dues formes ecològiques de Culex (Culex) pipiens Linné 1758 (Diptera: Culicidae) al Baix Llobregat. Barcelona: University of Barcelona; 1998.

51. WHO. Global plan for insecticide resistance management in malaria vectors. Geneva: World Health Organization; 2012.

52. Glunt KD, Blanford JI, Paaijmans KP. Chemicals, climate, and control: Increasing the effectiveness of malaria vector control tools by considering relevant temperatures. PLoS Pathog. 2013; 9(10): e1003602. https://doi.org/10.1371/journal.ppat.1003602 PMID: 24098111

53. Owusu HF, Chitnis N, Müller P. Insecticide susceptibility of Anopheles mosquitoes changes in response to variations in the larval environment. Sci Rep. 2017; 7(1):3667. https://doi.org/10.1038/s41598-017-03918-z PMID: 28623302

54. Rajatileka S, Burhani J, Ranson H. Mosquito age and susceptibility to insecticides. Trans R Soc Trop Med Hyg. 2011; 105(5):247-53. https://doi.org/10.1016/j.trstmh.2011.01.009 PMID: 21353689

55. Oliver S, Brooke B. The effect of multiple blood-feeding on the longevity and insecticide resistant pheno-type in the major malaria vector Anopheles arabiensis (Diptera: Culicidae). Parasit Vectors. 2014; 7(1):390.

56. Alout H, Yameogo B, Djogbénou LS, Chandre F, Dabiré RK, Corbel V, et al. Interplay between *Plasmo-dium* infection and resistance to insecticides in vector mosquitoes. J Infect Dis. 2014; 210(9):1464-70. https://doi.org/10.1093/infdis/jiu276 PMID: 24829465

57. Viana M, Hughes A, Matthiopoulos J, Ranson H, Ferguson HM. Delayed mortality effects cut the malaria transmission potential of insecticide-resistant mosquitoes. Proc Natl Acad Sci USA. 2016; 113(32):8975-80. https://doi.org/10.1073/pnas.1603431113 PMID: 27402740

58. Rivero A, Vézilier J, Weill M, Read AF, Gandon S. Insecticide Control of Vector-Borne Diseases: When Is Insecticide Resistance a Problem? PLoS Pathog. 2010; 6(8):e1001000. https://doi.org/10.1371/journal.ppat.1001000 PMID: 20700451

59. Thomas MB, Read AF. The threat (or not) of insecticide resistance for malaria control. Proc Natl Acad Sci USA. 2016; 113(32):8900-2. https://doi.org/10.1073/pnas.1609889113 PMID: 27482081

60. Huijben S, Paaijmans KP. Putting evolution in elimination: Winning our ongoing battle with evolving malaria mosquitoes and parasites. Evolutionary Applications. 2018; 11:415-30. https://doi.org/10.1111/eva.12530 PMID: 29636796

61. WHO. Space spray application of insecticides for vector and public health pest control: A practitioner's guide. Geneva: World Health Organization 2003.

62. European Centre for Disease Prevention and Control. Vector control with a focus on *Aedes aegypti* and *Aedes albopictus* mosquitoes: literature review and analysis of information. Stockholm: ECDC, 2017.

63. WHO. Global vector control response 2017-2030. Geneva: World Health Organization, 2017.