How varying parameters impact insecticide resistance bioassay: An example on the worldwide invasive pest Drosophila suzukii

Monitoring pesticide resistance is essential for effective and sustainable agricultural practices. Bioassays are the basis for pesticide-resistance testing, but devising a reliable and reproducible method can be challenging because these tests are carried out on living organisms. Here, we investigated five critical parameters and how they affected the evaluation of resistance to the organophosphate phosmet or the pyrethroid lambda-cyhalothrin using a tarsal-contact protocol on Drosophila suzukii, a worldwide invasive pest. Three of the parameters were related to insect biology: (i) sex, (ii) age of the imago (adult stage) and (iii) genetic diversity of the tested population. The two remaining parameters were linked to the experimental setup: (iv) the number of individuals tested per dose and (v) the duration of exposure to the active ingredient. Results showed that response to insecticide differed depending on sex, males being twice as susceptible to phosmet as females. Age principally affected young females’ susceptibility to phosmet, because 0–24 hour-old flies were twice as susceptible as 24–48 hour-old and 72–96 hour-old females. Genetic diversity had no observable effect on resistance levels. The precision and accuracy of the median lethal dose (LD50) were greatly affected by the number of individuals tested per dose with a threshold effect. Finally, optimal duration of exposure to the active ingredient was 24 h, as we found an underestimation of mortality when assessed between 1 and 5 h after exposure to lambda-cyhalothrin. None of the main known point mutations on the para sodium channel gene associated with a knockdown effect were observed. Our study demonstrates the importance of calibrating the various parameters of a bioassay to develop a reliable method. It also provides a valuable and transferable protocol for monitoring D. suzukii resistance worldwide.


Introduction
Bioassays are considered as the gold standard for pesticide resistance testing. They can detect new resistances and can assess resistance in an integrative manner, whether the underlying mechanism is known or unknown [1]. With the growing ambition to employ sustainable practices of pest control, monitoring pesticide resistance is becoming the key to implementing efficient, integrated pest management strategies. Accurate and reliable phenotypic data, in addition to molecular data, are also essential to address numerous theoretical evolutionary questions linked to pesticide resistance (see [2] for examples with insecticides). As in most experiments involving living material, the output of pesticide resistance bioassays can be affected by many parameters, depending on the type of bioassay used and the study species.
When dealing with a new species/active substance pair, dose-response analyses are the most frequently used approach. The standard procedure consists in exposing insects to a range of doses of the active substance of interest and in counting the number of dead individuals for each dose. Using regression methods, it is then possible to determine various parameters of interest, and in particular the median lethal dose (LD 50 ), the dose that kills 50% of the exposed individuals. Results of bioassays may be affected by experimental conditions and design [3] and also by the biological characteristics of the insect species [4][5][6][7][8]. For each pest, it is necessary to develop, or at least to adapt, a standard experimental protocol to obtain reliable and repeatable results that will allow comparisons of pesticide susceptibility over space and time.
Because most insecticides targeting D. suzukii belong to the group of contact insecticides, we chose to use tarsal-contact bioassays. This type of bioassay is easy to set up and relatively inexpensive, facilitating its use worldwide. Organophosphates and pyrethroids have been the most effective and the most used insecticides in the field in France [42][43][44]47, 49]; we therefore focused our experiments on two active substances of these families commonly used in Europe against D. suzukii: the organophosphate phosmet and the pyrethroid lambda-cyhalothrin. The formulated products act via contact, ingestion and inhalation for phosmet and via contact and ingestion for lambda-cyhalothrin.
The main aim of this study was to illustrate and to investigate how various biological and technical parameters can affect the accuracy and reliability of a pesticide resistance bioassay.
The key biological parameters we explored were (i) sex, (ii) age of the imago (adult stage) and (iii) genetic diversity of the tested population. Two methodological parameters were also tested: (iv) number of individuals used to evaluate the level of insecticide resistance and (v) duration of exposure to the active substance. One interesting and valuable output of this work is a reliable, accurate and standardized tarsal-contact bioassay to test insecticide resistance.
Our work highlights the essential factors to control in the development of reliable bioassays on any other pest species.

Sample collection and preparation
The Ste-Foy population used in the experiments was originally established from approximately 20 D. suzukii females, collected from southeastern France in an urban area on raspberries (Sainte-Foy-les-Lyon, France) in May 2012 by Roland Allemand (Biometry and Evolutionary Biology Laboratory, CNRS-University of Lyon). The population was then mass-reared and maintained in standard drosophila vials (Ø 25 mm x h 95 mm) containing~10 ml of a standard food media composed of 10 g of agar diluted in 1 l of water, 60 g of glucose, 30 g of saccharose, 80 g of malted yeast, 20 g of yeast extract, 20 g of peptones, 0.5 g of magnesium sulfate, 0.5 g of calcium chloride and 1 g of nipagin, pre-diluted in 10 ml of ethanol. They were maintained in a climatic chamber at 23±1˚C, under a relative humidity of approximately 70% and a 16:8 h light-dark cycle. Temperature was strictly controlled because it can influence D. suzukii susceptibility to insecticides [54]. Then, 25-30 males and females per vial were left to mate and oviposit for four to seven days before being removed from the rearing vials. After two weeks, the newly emerged adults (imagoes) were isolated and transferred to a new vial to start a new generation. To maintain the Ste-Foy population, each generation consisted of approximately 20 vials and the emerged adults of all the vials were mixed together to maximize genetic variability before distributing them into new vials.
A low genetic diversity population, SF-IsoA, was generated from the Ste-Foy population by inbreeding the flies for three generations: each generation, a single brother and virgin sister pair were isolated and left to reproduce in a new vial. The SF-IsoA population was then maintained in the same conditions as the Ste-Foy population.

Tarsal-contact bioassays
All the experiments were conducted using a common framework for the bioassays. Based on this framework, the five parameters of the protocol were modified to test for their effect on the assessed pesticide resistance. The tarsal-contact bioassay consists in exposing fly distal part of the legs (tarsi) to insecticide for a pre-determined time before assessing mortality. We used 20 ml scintillation glass vials (Ø 28 mm x h 61 mm), which have previously been used on D. suzukii with good results for resistance monitoring [48,52]. A volume of 500 μl of the insecticide solution at the desired concentration (dissolved in acetone) was deposited on the walls of the vials. The insecticide solution was uniformly distributed with acetone evaporation by rolling the vials at room temperature for 1 h, on a hot-dog roller. Seven doses of the active substance were tested, including a control with acetone only. The range of concentrations differed for the two insecticides and depended on what concentrations allowed a dose-response curve ranging from 0 to 100% mortality. For the organophosphate phosmet (95% of purity), an acetylcholinesterase inhibitor and obtained from Gowan, the chosen range of concentrations was 0, 4.7, 9.4, 18.9, 37.7, 75.5 and 150.9 mg/l. Lambda-cyhalothrin (92.8% of purity) is a pyrethroid that targets the sodium channel and was furnished by Syngenta France SAS. The selected concentration range was adjusted according to the bioassay and consisted in 7 doses chosen from the following concentrations: 0, 0.01, 0.05, 0.1, 0.25, 0.5, 1, 3 and 5 mg/l. Because contact time lasted up to 24 h, 100 μl of an agarose-sucrose mixture (5% sucrose, 8 g/l agarose) was provided in each vial to prevent death by starvation. For a bioassay, one to four vials for the same dose were used, depending on the number of flies tested. The number of D. suzukii flies of a uniform age class (the standard being 24 to 48 hours old in this study) from a single population ranged from 3 to 32 by vial. Males and females were not separated prior to the bioassay to avoid anesthesia (necessary to determine the sex of the flies). The vials were plugged with thin netting held in place by perforated caps to prevent flies from escaping and still allow good ventilation. The vials were maintained in a climatic chamber during the bioassay at 20 ±1˚C, under a relative humidity of approximately 75% and a 16:8 h light-dark cycle. After a certain duration of exposure (the standard being 24 h), the vials were briefly shaken and the numbers of dead, moribund and live flies were assessed. Individuals that could not remain on their legs and those that showed unusual behavior (i.e. uncertain or irregular flight, twitching legs and/or uneven movements) were considered as moribund. Immobile adult flies were considered dead. The sex of the flies for the different categories was determined during this mortality assessment by checking for the presence of black spots on the wings which characterizes males. The susceptibility of a population to the chemical was assessed by calculating the LD 50 value (see Statistical analysis section below).
This protocol was used for conducting five experiments with some variation regarding the D. suzukii individuals or the conditions of the bioassay (see Table 1). For each experiment, several bioassays were performed and each bioassay was done on a different date. Within a bioassay, multiple vials per dose could be used in order to have a limited number of flies per vial (up to 32). For Experiment 1 (24 bioassays), the test aimed to compare the susceptibility to phosmet according to sex at 24 h of exposure on flies from Ste-foy population aged of 24 to 48 h. In Experiment 2 (26 bioassays), we explored the influence of age class on susceptibility of Ste-Foy population to 24h of exposition to phosmet by testing males and females of three different age classes: 0 to 24 h, 24 to 48 h and 72 to 96 h. Experiment 3 (5 bioassays) consisted in observing the impact of the genetic diversity of a population on its resistance to phosmet after 24h of exposure. To do so, the reference population Ste-Foy was compared to the inbred SF-IsoA population and flies were 24 to 48 h old. The respective levels of genetic diversity of the two populations were assessed using molecular tools (see below). Experiment 4 (24 bioassays) investigated the impact of the number of flies tested per bioassay (mean number of flies per dose per bioassay ranging from 6 to 37 for the females and from 9 to 35 for the males, see Table 2) on the accuracy of the LD 50 estimates after 24 h of exposure of 24 to 48 h old flies from Ste-Foy to phosmet. Finally, Experiment 5 (3 bioassays) explored the effect of the duration exposure to pyrethroid insecticide (lambda-cyhalothrin) on the evaluation of susceptibility of Ste-Foy population. Unlike organophosphate insecticides, pyrethroid can potentially induce a knockdown effect (due to several mutations in the sodium channel molecular target), which may cause an erroneous mortality assessment depending on the time of observation after initial exposure to the chemical. That is why we conducted the bioassays with lambda-cyhalothrin in this Experiment 5. Mortality was assessed repeatedly on the same vials at 10 different times: 1, 2, 3, 4, 5, 20, 21, 22, 23 and 24 h after initial exposure. Assessing the neutral and adaptive genetic diversity of the populations DNA extraction. Adults (aged from 24 to 72 h old) were picked from the rearing stock in order to assess the genetic diversity of Ste-Foy and SF-IsoA populations. Flies were killed directly in 70% ethanol and the total DNA was extracted from the whole body within a week. DNA from 44 individual D. suzukii samples (30 females, 14 males) from each population (Ste-Foy and SF-IsoA) was extracted following a modified method based on Walsh et al. (1991) [55]. A volume of 100 μl of a 10% Chelex 100 solution (Bio-Rad) and 3% of 10 mg/ml proteinase K (Eurobio) was added to each sample. Then, each sample was crushed using 2 mm steel beads on a 1600 MiniG tissue homogenizer (Spex1 SamplePrep) at 1500 strokes/min for 5 sec, in a 96-well format PCR plate. Tissues were digested for 14 h at 56˚C with a Mastercycler thermocycler (Eppendorf) with a final temperature step of 30 min at 98˚C, the supernatant was used as DNA template for PCR reaction.
Microsatellite genotyping. The genetic diversity of the two studied populations was assessed using microsatellite markers. We used a slightly modified nested PCR approach described by Schuelke [56] coupled with a selection of 13 microsatellite markers developed by Fraimout et al. [57] (S1 Table). A forward specific primer was conjugated with a 5'-GTTGTAAAACGACGGCCAGT-3' M13-tail at its 5' end, the same labeled universal M13-tail was used for fluorescence detection. We used a 1:10 ratio with 1 unit tailed forward specific primer to 10 units labeled M13 tails and 10 specific reverse primers. PCR amplifications were carried out on a Mastercycler thermocycler (Eppendorf) in a 12 μl reaction volume containing 1X GoTaq1 Flexi Buffer, 1.5 mM MgCl 2 , 0.1 mg/ml bovine serum albumin (BSA), 200 μM of each dNTP, 0.4 μM of the labelled M13 tail, 0.4 μM of the specific reverse primer, 0.04 μM of M13-tailed specific forward primer, 1 unit of GoTaq1 Flexi DNA Polymerase (Promega) and 2 μl of DNA template. The PCR conditions were: 3 min at 95˚C followed by 30 cycles at 95˚C for 30 sec, 57˚C for 45 sec, 72˚C for 45 sec and 10 cycles at 95˚C for 30 sec, 54˚C for 30 sec, 72˚C for 45 sec with a final extension step at 72˚C for 20 min.
Checking for the kdr mutation. Knockdown resistance to pyrethroids is known to be induced by some mutations in its molecular target, the voltage-gated sodium channel. To assess the presence of the main kdr mutations that lead to the L1014F substitution [58], we adapted a PCR-RFLP protocol from Franck et al. [59]. Twenty adults (aged from 24 to 72 h old) were randomly picked from the rearing stock of Ste-Foy population as well as 20 flies from the same population that survived an exposition of 24 h at 0.25 mg/l of lambda-cyhalothrin (from Experiment 5) and their DNA was extracted as described above. A 371 bp PCR fragment was amplified on a Mastercycler thermocycler (Eppendorf) in a 12 μl reaction volume containing 1X GoTaq1 Flexi Buffer, 1.5 mM MgCl 2 , 0.1 mg/ml BSA, 200 μM of each dNTP, 0.4 μM of each primer (CKDR1 and Cgd2; see Table 3), 1 unit of GoTaq1 Flexi DNA Polymerase (Promega) and 2 μl of DNA template. The PCR conditions were: 3 min at 95˚C followed by 35 cycles at 95˚C for 30 sec, 56˚C for 45 sec, 72˚C for 45 sec with a final extension step at 72˚C for 20 min. The PCR products (5 μl) were digested at 37˚C for 16 h with 2 units of MluCI endonuclease and 1X of NEB Buffer (New England Biolabs) in 10 μl of reaction volume. The digested products (4 μl) were separated on a 2% agarose gel at 100 V for 30 min using the RunOneTM Electrophoresis System (Embi Technology). The size of the fragments was estimated by comparison with a 100 bp DNA ladder (Promega). The L1014F mutation reveals a restriction site recognized by the MluCI endonuclease, which cuts the 371 bp PCR fragment into two fragments of respectively 123 bp and 248 bp. Susceptible genotypes remain undigested.
Primers Ds-SKdr-F and Ds-Seq-R were designed specifically for this study. Primers CKDR1 and Cgd2, designed in our laboratory and used in previous studies [60,61], were used as is given their perfect homologies.
Voltage-gated sodium channel gene sequencing. To verify the absence of the L1014F substitution and other secondary substitutions previously described in other insect species (i.e.  Table 1), 1 unit of GoTaq1 Flexi DNA Polymerase (Promega) and 2 μl of DNA template. The PCR conditions were: 3 min at 95˚C followed by 35 cycles at 95˚C for 30 sec, 54˚C for 45 sec, 72˚C for 2 min with a final extension step at 72˚C for 20 min. The PCR products (20 μl) were sequenced (Eurofins Genomics) using primers Ds-SKdr-F, Ds-Seq-R and CKDR1 (see Table 3). DNA sequences were manually aligned and analyzed using Bioedit software [66].

Statistical analysis
Moribund and dead individuals were combined and considered as dead. Prior to data analyses, we kept only bioassays with a mortality rate in the control less than 15%. Survival data were analyzed in a non-linear regression framework using the 'drc' package [67] in R [68]. For all experiments, we assumed a binomial distribution of errors. For Experiments 1, 2 and 3, survival data were fitted to the three-parameter log-normal model, which is equivalent to the classic probit model [69] with an additional parameter that takes into account the 'natural' mortality rate observed in the control of each categories. Due to the very low mortality rate for the different categories of Experiment 4 and Experiment 5 (<10%) as well as the sometimes limited number of individuals tested, survival data were fitted to the two-parameter log-normal model equivalent of the classic probit model [69]. For Experiments 1, 2 3 and 5, the results of the bioassays were homogeneous and were therefore pooled together. The results of the bioassays included in Experiment 4 were not pooled because the purpose of this experiment was to assess the effect of the mean number of individuals in the bioassay on the LD 50 estimates. Models were fitted for the different categories separately within each experiment. The 'drc' package allows the estimation of LD 50 and the associated 95% CI and standard errors. Difference in LD 50 between males and females (Experiment 1) was tested directly using the 'compParm' function implemented in the 'drc' package. For the other experiments, the datasets were split into male and female subsets prior to the analysis. Pairwise comparisons of LD 50 between the different age classes, the populations with different levels of genetic diversity and the different times of exposure were performed for each sex separately for Experiments 2, 3 and 5, respectively, using the 'compParm' function. LD 50 estimated on each bioassay of Experiment 4 separately were compared with the LD 50 estimated on the bioassays pooled by sex based on the overlapping of the 95% CI and the respective values. Additionally, for Experiment 3, two genetic diversity indices were computed for the Ste-Foy and the SF-IsoA populations. The number of alleles (N a ) and gene diversity (H e ) [70] were estimated for both populations and for all loci with GENEPOP V4.2.2 [71]. The dataset and the code used for the different analyses, as well as for the production of the figures, have been deposited in an online repository (doi: 10.5281/zenodo.2842939).

Effect of biological parameters on LD 50 estimates: Sex, individual fly age and population genetic diversity
The influence of sex on insecticide susceptibility. A total of 2745 females and 2341 males were tested in 24 bioassays. The LD 50 values for phosmet were estimated at 39.6 (95% CI: 37.2-42.0) and 19.8 (95% CI: 18.6-21.1) mg/l for females and males, respectively (Fig 1). The two-fold higher resistance to phosmet of females compared with males was highly significant (relative female:male potency = 1.99, t-value = 11.3, p-value < 0.001).
The  The influence of genetic diversity on insecticide susceptibility. Microsatellite genotyping confirmed that the genetic diversity of the SF-IsoA population was reduced compared with the Ste-Foy population at 13 loci (Fig 3 and S2 Table). The mean number of alleles (per microsatellite marker) and the mean gene diversity was higher in Ste-Foy (N a = 3.31; H e = 0.595) than in SF-IsoA (N a = 1.46; H e = 0.179), confirming that Ste-Foy and SF-IsoA populations are appropriate biological materials to test the influence of genetic diversity on resistance.

PLOS ONE
How varying parameters impact insecticide resistance bioassay

Effect of the experimental setup on LD 50 values: Number of tested individuals and exposure duration
The influence of the number of tested individuals on insecticide susceptibility. We analyzed the variation in the estimated LD 50 values according to the number of tested individuals per bioassay per dose of phosmet (from a minimum mean number per dose of 9.4 and 10 flies for males and females, respectively, to a maximum mean number per dose of 34.9 and 37.1 for males and females, respectively). Our results clearly highlight the importance of the number of tested flies on the accuracy of the estimation of LD 50 (Fig 5). For both females and males, an insufficient number of tested individuals led to imprecise and haphazard evaluations of LD 50 (Fig 5). The standard errors associated with the LD 50 values estimated with a limited number of individuals tended to be higher when fewer individuals were included in the experiment. In contrast, when the mean number of individuals tested per dose exceeded a threshold of approximatively 30, the LD 50 assessed for the various bioassays were similar and close to the value estimated with every tested individual pooled together as one unique total bioassay according to sex. Males with 29, 30 and 35 individuals per dose led to a similar estimation of LD 50 values with overlapping 95% confidence intervals: 21. The influence of the duration of insecticide exposure on D. suzukii susceptibility. Exposure duration to lambda-cyhalothrin had a significant effect on the LD 50 estimates in the study population, both for females and males. The three bioassays were very similar: only 2 out of 30 LD50 comparisons between bioassays at each time of exposure were significantly different. For this reason, subsequent analyses were performed on pooled data from the bioassays. During the first phase of the experiment (from 1 to 5 h duration), we observed a decrease in the LD 50 estimate. In the second phase, the LD 50 estimate was stable (from 20 to 24 h duration) (Fig 6A). The duration of exposure before scoring had a significant effect on the LD 50 estimate (e.g. relative potency 5h:24h = 1.55, t-value = 2.53, p-value = 0.011) (Fig 6B). There was a clear significant break between 5 h and 20 h of exposure time. In particular, the estimation of LD 50 for 1, 2, 3 or 4 h exposure was less precise than after longer exposure times, with larger standard errors and 95% confidence intervals. The LD 50 estimates for 20 to 24 h of exposure were similar (relative potency close to 1) and not significantly different (Fig 6B).
Detection of putative mutations in the voltage-gated sodium channel. The genotyping of 20 random individuals from Ste-Foy and 20 individuals from Ste-Foy that survived for 24 h at 0.25 mg/l of lambda-cyhalothrin did not reveal the presence of mutant L1014F kdr alleles. Therefore, the results observed in our lambda-cyhalothrin assays were not biased by this mutation. The partial sequencing of the gene encoding for the voltage-gated sodium channel of two flies from the Ste-Foy population that survived at 0.25 mg/l of lambda-cyhalothrin and exhibited a kdr-like phenotype did not reveal the presence of any of the 12 main resistance mutations described in the literature (GenBank accession nos. MK645039-MK645040).

Discussion
Our study explored a selection of critical parameters to control in pest bioassays. From an applied point of view, this study defined a useful reproducible protocol for a reliable tarsalcontact bioassay to evaluate insecticide resistance in D. suzukii.
We found that D. suzukii does not respond in the same way to insecticide according to sex and age. In most studies on insecticide resistance of D. suzukii, tests are done on females only or on a mix of males and females without taking sex into account. A few experiments have reported that males have greater insecticide susceptibility than females [44,48,52]. This phenomenon has also been observed in various insect species and specifically in Diptera for various insecticides, including pyrethroids and organophosphates [72][73][74][75]. We observed the same phenomenon in our experiment with males being twice as susceptible as females. The origin of this difference has already been discussed in other studies, being sometimes attributed to the difference in size between sexes due to sexual dimorphism [76][77][78], although detoxification enzymes also seem to be involved in these differences [79,80]. Our results also showed an influence of age on susceptibility, but only for females. Susceptibility to pesticide can vary according to the development stage of the insect [81]. Smirle et al. [48] reported that malathion (an organophosphate) induces more toxicity on 5-8-day-old D. suzukii adults than on 2-dayold individuals. Most experiments use flies at least 5 days old, which is time-consuming in terms of sample preparation. No experiment has indicated insecticide susceptibility in D. suzukii flies younger than 24 h. In our experiment, newly hatched females (0-24 h) showed significantly higher susceptibility to insecticides than older age classes. This period corresponds to the sexual maturation stage in this synovigenic species [82,83]. During the ovarian maturation period, metabolic detoxification pathways and energy investments are different compared with the rest of the adult life [84,85]. Testing the effects of mating status on insecticide resistance might be an interesting avenue of research. Mating occurs mainly after the female refractory period of 1 to 2 days after emergence [82,86] and generates strong changes in female physiology and behavior in Drosophila [87,88]. On the contrary, young D. melanogaster flies have been described as more resistant to several abiotic stresses [89]. Our finding is of interest for establishing a reliable bioassay method, but may also have an applied interest in targeting treatments on the early life stage of young females.
Inbred populations with low genetic diversity may display lower resistance to pathogenic organisms [90] and may be more sensitive to environmental stress [91]. By creating a low genetic diversity population SF-IsoA, we were able to test the hypothesis that reduced genetic diversity may (i) reduce the variability of insecticide response estimates (i.e. LD 50 values) within the population and increase repeatability, and (ii) genetically fix the strain and obtain a stable reference population through generations. We therefore expected a different insecticide response in the SF-IsoA compared with the Ste-Foy population. However, no differences in susceptibility were detected in our comparison between original and inbred populations. Our finding is thus likely an illustration of the fact that the Ste-Foy population is initially devoid of resistance alleles and therefore allelic paucity does not change this status.
Regarding the effect of the experimental setup on insecticide susceptibility, the number of tested flies per bioassay and per dose and the duration of exposure to the tested insecticide are essential parameters to control. As suggested in previous studies, a range from 10 to 25 insects per dose is necessary [92,93]. Our results showed that the number of tested insects is the most important factor among those explored in this work. According to our data, a reliable mean number of tested insects per dose is a minimum of 29 males and 34 females. We therefore recommend about 30 individuals per dose for the tests to be accurate and reproducible. For most of the existing protocols on D. suzukii, the duration of exposure to the chemical is usually about 24 h [44, 48, 50, 51] and can reach 48 or 72 h [4,43,45,47]. A study on the mortality induced by ingestion of spinosad on this species emphasized the importance of time elapsed (5 days) before assessing mortality, suggesting that exposure duration is crucial in different types of bioassays [53]. Bioassays can be extremely time consuming, making it tempting to minimize exposure duration to reduce the total testing time. However, for pyrethroids, several species have shown a knockdown phenotype that can induce a transient state of apparent mortality in the insects a few minutes after exposure to the insecticide. In a study of mosquitoes exposed to lambda-cyhalothrin, the time to 50% knockdown (KT 50 ) ranged from 14 (Mansonia africana) to 128 min (Anopheles gambiae) with KT 95 values of 42 and 277 min, respectively [94]. In another study, the KT 50 for D. melanogaster flies ranged from 11 to 40 min, depending on the pyrethroid and the strain [95]. Therefore, assessing the mortality too early after exposure may distort results, while a late, albeit reliable, assessment [96] may result in lost time. We observed a few individuals with a kdr-like phenotype in our experiment. They showed temporary paralysis followed by complete recovery of their locomotive capacities. Additional research did not reveal any of the main resistance mutations or the kdr mutation at the expected locus (1014) of the para gene. Similar results have been reported from a study on D. suzukii adults exposed to spinosad [52]-to which resistance is emerging in the US [37]. This study's results showed highly variable LD 50 values, calculated from the large number of moribund flies due to an insufficient exposure time (6 h). To avoid potential misinterpretation due to a knockdown effect, a rapid detoxification mechanism (potentially present in susceptible insects) or a slow mode of action of an active substance, we recommend assessing insect mortality at least 24 h after insecticide exposure (or more depending on the mode of action of the insecticide). Because the phenotypic effect of the insecticide is estimated by visual assessment, it is of utmost importance to establish precise scoring criteria that can be standardized across experiments and laboratory workers, in particular to differentiate moribund individuals from the dead and live individuals.
As illustrated by the low to moderate levels of resistance of D. suzukii populations that have been described in the US [37], monitoring insecticide resistance of this pest is important to be able to react quickly and adapt control management. We hereby propose a validated and operational method for a reliable worldwide monitoring of D. suzukii resistance to insecticides, and highlighted several parameters that are essential to control for in the design of resistance assessment bioassays on other species as well.