Search

Abstract Background The western corn rootworm, Diabrotica virgifera virgifera, is a pervasive pest of maize in North America and Europe, which has adapted to current pest management strategies. In advance of an assembled and annotated D. v. virgifera genome, we developed transcriptomic resources to use in identifying candidate genes likely to be involved in the evolution of resistance, starting with members of the ATP-binding cassette (ABC) transporter family. Results In this study, 65 putative D. v. virgifera ABC (DvvABC) transporters were identified within a combined transcriptome assembly generated from embryonic, larval, adult male, and adult female RNA-sequence libraries. Phylogenetic analysis placed the deduced amino-acid sequences of the DvvABC transporters into eight subfamilies (A to H). To supplement our sequence data with functional analysis, we identified orthologs of Tribolium castaneum ABC genes which had previously been shown to exhibit overt RNA interference (RNAi) phenotypes. We identified eight such D. v. virgifera genes, and found that they were functionally similar to their T. castaneum counterparts. Interestingly, depletion of DvvABCB_39715 and DvvABCG_3712 transcripts in adult females produced detrimental reproductive and developmental phenotypes, demonstrating the potential of these genes as targets for RNAi-mediated insect control tactics. Conclusions By combining sequence data from four libraries covering three distinct life stages, we have produced a relatively comprehensive de novo transcriptome assembly for D. v. virgifera. Moreover, we have identified 65 members of the ABC transporter family and provided the first insights into the developmental and physiological roles of ABC transporters in this pest species.

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 Drosophila melanogaster 91-R and 91-C strains are of common origin, however, 91-R has been intensely selected for dichlorodiphenyltrichloroethane (DDT) resistance over six decades while 91-C has been maintained as the non-selected control strain. These fly strains represent a unique genetic resource to understand the accumulation and fixation of mutations under laboratory conditions over decades of pesticide selection. Considerable research has been done to investigate the differential expression of genes associated with the highly DDT resistant strain 91-R, however, with the advent of whole genome sequencing we can now begin to develop an in depth understanding of the genomic changes associated with this intense decades-long xenobiotic selection pressure. Here we present the first whole genome sequencing analysis of the 91-R and 91-C fly strains to identify genome-wide structural changes within the open reading frames. Between-strain changes in allele frequencies revealed a higher percent of new alleles going to fixation for the 91-R strain, as compared to 91-C (P