Insecticide resistance in the bed bug comes with a cost.

Adaptation to new environmental stress is often associated with an alteration of one or more life history parameters. Insecticide resistant populations of insects often have reduced fitness relative to susceptible populations in insecticide free environments. Our previous work showed that three populations of bed bugs, Cimex lectularius L., evolved significantly increased levels of resistance to one product containing both β-cyfluthrin and imidacloprid insecticides with only one generation of selection, which gave us an opportunity to explore potential tradeoffs between life history parameters and resistance using susceptible and resistant strains of the same populations. Life history tables were compiled by collecting weekly data on mortality and fecundity of bugs from each strain and treatment throughout their lives. Selection led to a male-biased sex ratio, shortened oviposition period, and decreased life-time reproductive rate. Generation time was shortened by selection, a change that represents a benefit rather than a cost. Using these life history characteristics we calculated that there would be a 90% return to pre-selection levels of susceptibility within 2- 6.5 generations depending on strain. The significant fitness costs associated with resistance suggest that insecticide rotation or utilization of non-insecticidal control tactics could be part of an effective resistance management strategy.

Adaptation to a new environmental stress is often associated with an alteration of one or more life history parameters12. Ultimately, these tradeoffs may be the result of a physiological constraint, such as shunting resources into survival in the new environment; resources that could then not be used for egg production or rapid development. Genetic correlations between life history traits and genes that influence resistance may explain the intergenerational response. In the case of insecticide resistance, increased production of enzymes leading to insecticide detoxification or increased production of cuticular components that reduce penetration of the toxicant may have correlated effects on life history characters. Among various species, reduced longevity, delayed maturation and decreased egg production have been observed to accompany insecticide resistance345678. The net result of the adaptations to insecticide exposure is enhanced fitness in insecticide-treated habitats, but decreased fitness in insecticide-free areas.

Such tradeoffs open the possibility of insecticide resistance management by rotation between compounds with not only different modes of action but also different modes of detoxification (and hence different physiological costs) or to non-insecticidal control tactics. When there are negative life history tradeoffs associated with resistance in insecticide-free areas, populations should revert toward susceptiblity91011. If no costs existed in resistant individuals relative to susceptible, insecticide resistance management via rotation would be ineffective12.

In the past ten to fifteen years there has been a resurgence of pyrethroid resistant populations of Cimex lectularius L., the bed bug13141516. Pyrethroid resistance has been found in North America, Australia, Asia and Europe and is widespread throughout the United States and presumably elsewhere17181920. In the United States, resistance to pyrethroid-only insecticides has prompted a shift by urban pest management professionals to commercial insecticide products containing both a pyrethroid and a neonicotinoid2122. These two classes of insecticides act at different target sites on the insect neuron2324.

Our previous work investigating the evolutionary response of multiple populations of bed bugs to one of these combination products, Temprid SC® (β-cyfluthrin and imidacloprid), showed that resistance began to evolve in one generation in the laboratory22. This rapid evolution under laboratory conditions gave us an opportunity to explore the hypothesis that life history costs would be associated with decreased susceptibility to the combination insecticide in an environment no longer containing the insecticide.

Results

Mortality levels caused by the selected exposure times to label-rate Temprid SC in the second filial generation of the unselected line were very similar to the parental lines as would be expected in the absence of genetic drift (Table 1). Mortality in selected lines of LA1, CIN1 and NY1 were reduced by the selection that had occurred two generations earlier. This historical selection affected some, but not all, measured life history parameters (Figs. 1 and 2, see Supplemental Table S1 for analysis of variance results). Longevity, per cent hatch and per cent reaching the adult stage were not affected by the history of selection (Figs. 1, 2A,B). Unexpectedly, the sex ratio at maturity was biased towards females in unselected lines and biased towards males in selected lines (Fig. 2C). The generation time (G) was significantly less in all populations with a history of insecticide exposure compared to those not selected (Fig. 2D). Generation times were significantly reduced by selection by 5.1, 10.3 and 20.2% in LA1, CIN1, and NY1, respectively. Reproductive rates (R) were significantly lower in selected lines compared to their unselected counterparts (Fig. 2E). Reproductive rates were significantly reduced by 15.5, 42.5 and 35.3% in LA1, CIN1 and NY1, respectively. Additionally, the oviposition duration from the first egg laid to the last egg laid was significantly shorter for selected compared to unselected groups (Fig. 2F). Oviposition durations were reduced by 13.8, 29.4, and 46.9% in LA1, CIN1, and NY1, respectively.

Table 1

Percent mortality ± s.e.m. of all strains in the selected and unselected (parental and F2) groups and Z-statistics for the F2 generation.

Straina	Parental	F2
		Unselectedb	Selected	ZRep 1c	ZRep 2
LA1	82.2 (±0.4)	83.3 (±6.7)	30.0 (±5.0)	7.20*	4.60*
CIN1	83.1 (±0.8)	81.7 (±0.0)	3.4 (±3.4)	7.20*	8.27*
NY1	86.4 (±1.7)	79.2 (±4.2)	36.7 (±8.4)	5.11*	4.38*

*p < 0.001

aMortality of strains was evaluated following the protocol of Gordon et al.22 by exposing groups of bugs to residual deposits of Temprid SC. Individuals from each strain were exposed for strain-specific exposure times (LA1 0.1 h, CIN1 1 h, NY1 19 h) calculated to kill 80 per cent of the population.

bAverage mortality of selection replicates 1 and 2 to strain specific ET80s.

cA test for a significant difference between selected and unselected proportions in the F2 generation for each replicate and strain [Analytical Software. Statistix 8.0 for Windows. Tallahassee, FL (2003)]. N = 60 for each replicate, strain and treatment.

Figure 1

Survival (l; left axis and curve) and fecundity (m; right axis and histogram) over time for unselected and selected strains. Survival analysis found no significant differences in longevity due to selection with insecticide or strain (see Supplemental Table S3 for full descriptive statistics from survival analysis). Adult survival reached 50% at 31.1 (±0.94), 31.6 (±2.60) and 34.0 (±3.30) weeks for the unselected LA1, CIN1 and NY1 strains, respectively, and at 32.5 (±1.58), 30.6 (±2.34) and 28.6 (±2.77) weeks for the paired selected strains. The average time from the egg to adult molt was 10.8 (±0.63), 10.0 (±0.70), 10.3 (±0.25) weeks for the unselected LA1, CIN1 and NY1 strains, respectively, and 10.5 (±0.87), 9.3 (±0.63), 10.3 (±0.85) weeks for the paired selected strains. The average number of eggs laid during an individual female’s lifetime was 165.8 (±17.7), 154.0 (±24.9), 126.8 (±11.6) for the unselected LA1, CIN1 and NY1 strains, respectively, and 134.2 (±11.5), 84.9 (±9.7), 65.7 (±13.9) for the paired selected strains.

Figure 2

Life history parameters investigated. Selection did not affect the per cent of eggs that hatched or the per cent of eggs that reached adulthood. However, selection significantly decreased the proportion of individuals that were female, generation time, reproductive rate and oviposition duration. All results from the analysis of variance are presented in Supplemental Table S1, and p-values depicted in the figure represents the effect of selection on the specific parameter investigated.

Discussion

Three different strains of bed bugs with different initial levels of susceptibility to Temprid SC incurred significant life history costs after selection with this pyrethroid/neonicotinoid combination product22. The unselected strains varied in their levels of resistance to Temprid SC as is reflected in the exposure times that were required to lead to ca. 80% mortality (i.e., LA1bed bugs with different insecticide susceptibility profiles found that insecticide resistant strains had fitness costs compared to more susceptible strains25. Further, and more convincing, support for the hypothesis of a cost associated with resistance comes from the observation that parental selection for resistance resulted in decreased Rs in the F2 of all three populations. Similar tradeoffs between insecticide resistance and life history parameters have been recorded in other insect pests34567826.

In this system, one possible mechanism of the observed costs could be a tradeoff between production of detoxifying enzymes and allocation of resources for fecundity. A previously published study investigating the molecular mechanisms of resistance in the CIN1 selected strain found that four cytochrome P450s and one carboxylesterase was significantly over expressed compared to the CIN1 unselected line27. Increased detoxification is likely at least one mechanism of resistance in the LA1 and NY1 selected strains as well. Research investigating mechanisms of resistance in many populations of C. lectularius found that the P450 class of enzymes frequently confers a level of resistance to pyrethroid insecticides162728293031. Given this information, resources may be shunted from fecundity to the production of detoxifying enzymes. This implies a significant metabolic cost associated with the production of these enzymes. Further experiments involving the use of RNA interference could be used to elucidate the molecular mechanism of the observed costs and their correlation with detoxification-mediated insecticide resistance.

The R ratios of selected and unselected lines are 0.70, 0.42 and 0.57 for LA1, CIN1 and NY1, respectively. Numbers greater than 1 suggest an advantage, whereas numbers less than one suggest a disadvantage. Thus, a reversion to susceptibility should occur rapidly in insecticide-free environments. We calculated the time course to 50 and 90% recovery of pre-insecticide selection levels of susceptibility using a model proposed earlier32.

denotes the time required to reach a significant degree of susceptibility, denotes generation time, denotes the proportion of the population required to be susceptible, denotes the initial proportion susceptible, denotes the fitness of the susceptible and denotes the fitness of the resistant. The Rs for each treatment and strain were used for the fitness variables and . A return to 50% of pre-selection levels of susceptibility would take 3.03, 4.88 and 0.71 generations for LA1, CIN1 and NY1, respectively. Additionally, the same model predicted that the number of generations required for a strain to return to 90% level of pre-selection susceptibility was 6.51, 5.95 and 2.06 for LA1, CIN1 and NY1, respectively. Given that Gs were less than 20 weeks, reversion to 90% pre-selection levels of resistance should occur within ca. 1 to 2.5 years in the absence of selection.

Currently, the combination pyrethroid/neonicotinoid products are some of the most effective choices for control in the field2122. In theory, rotation to products utilizing alternative modes of action could reverse resistance31011 assuming that alleles for susceptibility still exist in the population. However, such reversion is complicated by public intolerance for resident populations (which contrasts with agricultural resistance management11), limited gene flow that would help re-establish susceptibility33, limited choices of alternative chemical classes21 and independence of actions of pest managers when a more coordinated approach might be required. Despite the complications, integration of, or rotation among, both chemical and non-chemical approaches (heat, vacuuming, encasement of beds, etc.) will be necessary for sustained management of populations of bed bugs.

Material and Methods

Insects

Three strains of bed bugs were used for this study. The LA1 strain was collected from Los Angeles in 2007 and was susceptible to pyrethroids13. The strain CIN1 was originally collected from Cincinnati, OH in 2005 and was resistant to pyrethroids13. Subsequently, its originally high level of pyrethroid resistance has declined but not to the level of the susceptible colony LA127. The NY1 strain was collected from New York City, NY in 2007 and was resistant to pyrethroids17. However, a reversion toward susceptibility has also been recorded for this strain, though not to the degree of CIN127. For each strain, two samples of bugs were selected overtime, and two separate lineages of selected and unselected strains were initiated for each strain (LA1, CIN1 and NY1; Fig. 3) by exposing these strains to residual deposits of the pyrethroid/neonicotinoid combination product Temprid SC for a time calculated to kill 80% of the population (ET80) at the label rate22. All subsequent evaluations of susceptibility were performed using the same exposure times and bioassay. Bugs from an untreated F1 generation were used to establish an F2 generation (Fig. 3). Insects were housed in incubators away from any insecticide exposure at 26.7° C, 65 ± 5% RH, and a photoperiod of 14:10 (L:D) h. All bed bugs were fed weekly on defibrinated rabbit blood (Quad Five, Ryegate, MT) that was warmed to 39 °C with a circulating water bath using an artificial membrane feeder34.

Figure 3

Experimental design for selection experiment. Bugs were exposed to label rate Temprid SC® for a time calculated to kill 80% of the respective strain of bed bugs22. Two parental lines of selected and unselected bugs were created for each strain. Using the subsequent F2 generation from these lines that receive no further exposure to insecticide, two samples were taken from each replicate to generate data. Replicates were not synchronous and separated by weeks.

Life history variables

Life history data were collected for each strain and treatment. A group of 20 eggs was gathered within a 24 hour period from three to 12 females (Supplemental Table S2) and maintained within a small petri dish (5.1 cm diameter) lined with black filter paper. These eggs were allowed to hatch and the resulting individuals were fed weekly for the duration of their lives. Adult offspring were allowed to mate ad libitum. Longevity and fecundity were recorded weekly in order to calculate per cent hatching, per cent reaching adulthood, proportion of females, G, R and oviposition duration. In addition, this information was used to generate weekly survival (l; female specific survival could not be determined initially due to a lack of information about sex until eclosion to the adult stage) and fecundity (m; in this case oviposition rates). These recorded values were used to calculate net R as the number of female offspring per individual per generation. Because sex determination was only made during the adult stage, estimates of l, m, and R were based on an assumption of a 1:1 sex ratio of eggs laid by the F2 generation. The F2 generation for each replicate was followed for up to 70 weeks when the last individual died. Within strains, samples of replicates 1 and 2 were separated in time by 2 to 4 weeks.

Data analysis

A nested analysis of variance was used to investigate the effects of strain, treatment and replicate(strain) on the hatch rate of eggs laid by F1 mothers, the per cent of those eggs that reach adulthood, the observed sex ratios after reaching maturity, G, R and oviposition duration [Systat software. SYSTAT 13. San Jose, CA (2008)]. Net R was calculated as the sum of the weekly l*m for each replicate with the assumption of a 1:1 sex ratio of eggs laid by F2 females. Generation time was calculated by taking the summation of l*m*x divided by the summation of l*m35. A Z-test was performed to investigate the contrast in susceptibility to Temprid SC at the ET80 of unselected and selected bugs in the F2 generation [Analytical Software. Statistix 8.0 for Windows. Tallahassee, FL (2003)]. Survival analysis was performed to investigate the effects of selection and strain on longevity [JMP® version 11.2.0. SAS Institute Inc., Cary, NC (2013)].

Additional Information

How to cite this article: Gordon, J. R. et al. Insecticide resistance in the bed bug comes with a cost. Sci. Rep. 5, 10807; doi: 10.1038/srep10807 (2015).