Two point mutations (V419L and L925I) in the voltage-sensitive sodium channel alpha-subunit gene have been identified in deltamethrin-resistant bed bugs. A quantitative sequencing (QS) protocol was developed to establish a population-based genotyping method as a molecular resistance-monitoring tool based on the frequency of the two mutations. The nucleotide signal ratio at each mutation site was generated from sequencing chromatograms and plotted against the corresponding resistance allele frequency. Frequency prediction equations were generated from the plots by linear regression, and the signal ratios were shown to highly correlate with resistance allele frequencies (r2 > 0.9928). As determined by QS, neither mutation was found in a bed bug population collected in 1993. Populations collected in recent years (2007-2009), however, exhibited completely or nearly saturating L925I mutation frequencies and highly variable frequencies of the V419L mutation. In addition to QS, the filter contact vial bioassay (FCVB) method was established and used to determine the baseline susceptibility and resistance of bed bugs to deltamethrin and lambda-cyhalothrin. A pyrethroid-resistant strain showed >9,375- and 6,990-fold resistance to deltamethrin and lambda-cyhalothrin, respectively. Resistance allele frequencies in different bed bug populations predicted by QS correlated well with the FCVB results, confirming the roles of the two mutations in pyrethroid resistance. Taken together, employment of QS in conjunction with FCVB should greatly facilitate the detection and monitoring of pyrethroid-resistant bed bugs in the field. The advantages of FCVB as an on-site resistance-monitoring tool are discussed.

The adaptation of insect populations to insecticidal control is a continual threat to human health and sustainable agricultural practices, but many complex genomic mechanisms involved in this adaption remain poorly understood. This study applied a systems approach to investigate the interconnections between structural and functional variance in response to dichlorodiphenyltrichloroethane (DDT) within the Drosophila melanogaster strain 91-R. Directional selection in 6 selective sweeps coincided with constitutive gene expression differences in DDT resistant flies, including the most highly upregulated transcript, Unc-115 b, which plays a central role in axon guidance, and the most highly downregulated transcript, the angiopoietin-like CG31832, which is involved in directing vascular branching and dendrite outgrowth but likely may be under trans-regulatory control. Direct functions and protein–protein interactions mediated by differentially expressed transcripts control changes in cell migration, signal transduction, and gene regulatory cascades that impact the nervous system. Although changes to cellular stress response pathways involve 8 different cytochrome P450s, stress response, and apoptosis is controlled by a multifacetted regulatory mechanism. These data demonstrate that DDT selection in 91-R may have resulted in genome-wide adaptations that impacts genetic and signal transduction pathways that converge to modify stress response, cell survival, and neurological functions. This study implicates the involvement of a multigenic mechanism in the adaptation to a chemical insecticide, which impact interconnected regulatory cascades. We propose that DDT selection within 91-R might act systemically, wherein pathway interactions function to reinforce the epistatic effects of individual adaptive changes on an additive or nonadditive basis.

Adaptation of insect phenotypes for survival after exposure to xenobiotics can result from selection at multiple loci with additive genetic effects. To the authors’ knowledge, no selective sweep analysis has been performed to identify such loci in highly dichlorodiphenyltrichloroethane (DDT) resistant insects. Here we compared a highly DDT resistant phenotype in the Drosophila melanogaster (Drosophila) 91-R strain to the DDT susceptible 91-C strain, both of common origin. Whole genome re-sequencing data from pools of individuals was generated separately for 91-R and 91-C, and mapped to the reference Drosophila genome assembly (v. 5.72). Thirteen major and three minor effect chromosome intervals with reduced nucleotide diversity (π) were identified only in the 91-R population. Estimates of Tajima's D (D) showed corresponding evidence of directional selection in these same genome regions of 91-R, however, no similar reductions in π or D estimates were detected in 91-C. An overabundance of non-synonymous proteins coding to synonymous changes were identified in putative open reading frames associated with 91-R. Except for NinaC and Cyp4g1, none of the identified genes were the ‘usual suspects’ previously observed to be associated with DDT resistance. Additionally, up-regulated ATP-binding cassette transporters have been previously associated with DDT resistance; however, here we identified a structurally altered MDR49 candidate resistance gene. The remaining fourteen genes have not previously been shown to be associated with DDT resistance. These results suggest hitherto unknown mechanisms of DDT resistance, most of which have been overlooked in previous transcriptional studies, with some genes having orthologs in mammals.

The ATP-binding cassette (ABC) transporters represent a superfamily of proteins that have important physiological roles in both prokaryotes and eukaryotes. In insects, ABC transporters have previously been implicated in insecticide resistance. The 91-R strain of Drosophila melanogaster has been intensely selected with DDT over six decades. A recent selective sweeps analysis of 91-R implicated the potential role of MDR49, an ABC transporter, in DDT resistance, however, to date the details of how MDR49 may play a role in resistance have not been elucidated. In this study, we investigated the impact of structural changes and an alternative splicing event in MDR49 on DDT-resistance in 91-R, as compared to the DDT susceptible strain 91-C. We observed three amino acid differences in MDR49 when 91-R was compared with 91-C, and only one isoform (MDR49B) was implicated in DDT resistance. A transgenic Drosophila strain containing the 91-R-MDR49B isoform had a significantly higher LD50 value as compared to the 91-C-MDR49B isoform at the early time points (6 h to 12 h) during DDT exposure. Our data support the hypothesis that the MDR49B isoform, with three amino acid mutations, plays a role in the early aspects of DDT resistance in 91-R.

Dichloro-diphenyl-trichloroethane (DDT) resistance among arthropod species is a model for understanding the molecular adaptations in response to insecticide exposures. Previous studies reported that DDT resistance may involve one or multiple detoxification genes, such as cytochrome P450 monooxygenases (P450s), glutathione S-transferases (GSTs), esterases, and ATP binding cassette (ABC) transporters, or changes in the voltage-sensitive sodium channel. However, the possible involvement of microRNAs (miRNAs) in the post-transcriptional regulation of genes associated with DDT resistance in the Drosophila melanogaster strain 91-R remains poorly understood. In this study, the majority of the resulting miRNAs discovered in small RNA libraries from 91-R and the susceptible control strain, 91-C, ranged from 16–25 nt, and contained 163 precursors and 256 mature forms of previously-known miRNAs along with 17 putative novel miRNAs. Quantitative analyses predicted the differential expression of ten miRNAs between 91-R and 91-C, and, based on Gene Ontology and pathway analysis, these ten miRNAs putatively target transcripts encoding proteins involved in detoxification mechanisms. RT-qPCR validated an inverse correlation between levels of differentially-expressed miRNAs and their putatively targeted transcripts, which implies a role of these miRNAs in the differential regulation of detoxification pathways in 91-R compared to 91-C. This study provides evidence associating the differential expression of miRNAs in response to multigenerational DDT selection in Drosophila melanogaster and provides important clues for understanding the possible roles of miRNAs in mediating insecticide resistance traits.