The genomics of insecticide resistance

Towards holistic insect monitoring: species discovery, description, identification and traits for all insects

Holistic insect monitoring needs scalable techniques to overcome taxon biases, determine species abundances, and gather functional traits for all species. This requires that we address taxonomic impediments and the paucity of data on abundance, biomass and functional traits. We here outline how these data deficiencies could be addressed at scale. The workflow starts with large-scale barcoding (megabarcoding) of all specimens from mass samples obtained at biomonitoring sites. The barcodes are then used to group the specimens into molecular operational taxonomic units that are subsequently tested/validated as species with a second data source (e.g. morphology). New species are described using barcodes, images and short diagnoses, and abundance data are collected for both new and described species. The specimen images used for species discovery then become the raw material for training artificial intelligence identification algorithms and collecting trait data such as body size, biomass and feeding modes. Additional trait data can be obtained from vouchers by using genomic tools developed by molecular ecologists. Applying this pipeline to a few samples per site will lead to greatly improved insect monitoring regardless of whether the species composition of a sample is determined with images, metabarcoding or megabarcoding. This article is part of the theme issue 'Towards a toolkit for global insect biodiversity monitoring'.

Microbial Symbiont-Based Detoxification of Different Phytotoxins and Synthetic Toxic Chemicals in Insect Pests and Pollinators

Insects are the most diverse form of life, and as such, they interact closely with humans, impacting our health, economy, and agriculture. Beneficial insect species contribute to pollination, biological control of pests, decomposition, and nutrient cycling. Pest species can cause damage to agricultural crops and vector diseases to humans and livestock. Insects are often exposed to toxic xenobiotics in the environment, both naturally occurring toxins like plant secondary metabolites and synthetic chemicals like herbicides, fungicides, and insecticides. Because of this, insects have evolved several mechanisms of resistance to toxic xenobiotics, including sequestration, behavioral avoidance, and enzymatic degradation, and in many cases had developed symbiotic relationships with microbes that can aid in this detoxification. As research progresses, the important roles of these microbes in insect health and function have become more apparent. Bacterial symbionts that degrade plant phytotoxins allow host insects to feed on otherwise chemically defended plants. They can also confer pesticide resistance to their hosts, especially in frequently treated agricultural fields. It is important to study these interactions between insects and the toxic chemicals they are exposed to in order to further the understanding of pest insect resistance and to mitigate the negative effect of pesticides on nontarget insect species like Hymenopteran pollinators.

Identifying potential insecticide resistance markers through genomic-level comparison of Bemisia tabaci (Gennadius) lines

The invasive whitefly (Bemisia tabaci) MED is one of the most economically damaging plant pests. The extensive use of insecticide over decades has led to that the invasive B. tabaci MED has developed resistance to a wide range of insecticide classes, but little is known about the genetic background associated with resistance. To this end, we conducted a comparative genome-wide analysis of single-base nucleotide polymorphisms between MED whitefly lines collected from fields that were recently infested and an insecticide-susceptible MED whitefly line collected in 1976. First, low-coverage genome sequencings were conducted on DNA isolated from individual whiteflies. The sequencing results were evaluated using an available B. tabaci MED genome as a reference. Significant genetic differences were discovered between MED whitefly lines collected from fields that were recently infested and an insecticide-susceptible MED whitefly line based on the principal component analyses. Top GO categories and KEGG pathways that might be involved in insecticide resistance development were identified, and several of them have not been previously associated with resistance. Additionally, we identified several genetic loci with novel variations including Cytochrome P450 monooxygenases (P450s), UDP-glucuronosyltransferases (UGTs), Glutathione S-transferases (GSTs), esterase, carboxyl-esterases (COE), ABC transporters, fatty acyl-CoA reductase, voltage-gated sodium channels, GABA receptor, and cuticle proteins (CPs) that were previously reported to have close associations with pesticide resistance in well-studied insect groups that provide an essential resource for the design of insecticide resistance-linked loci arrays insecticide. Our results was obtained solely on resequencing genome data sets, more pesticide bio-assays combined with omics datasets should be further used to verify the markers identified here.

Transcriptomic analysis of the fat body of resistant and susceptible silkworm strains, Bombyx mori (Lepidoptera), after oral treatment with fenpropathrin

The widespread use of pyrethroid pesticides has brought serious economic losses in sericulture, but there is still no viable solution. The key to solving the problem is to improve silkworm resistance to pesticides, which depends on understanding the resistance mechanism of silkworms to pesticides. This study aimed to use transcriptomes to understand the underlying mechanism of silkworm resistance to fenpropathrin, which will provide a theoretical molecular reference for breeding pesticide-resistant silkworm varieties. In this study, the fat bodies of two strains with differential resistance after 12 h of fenpropathrin feeding were analyzed using RNA-Seq. After feeding fenpropathrin, 760 differentially expressed genes (DEGs) were obtained in the p50(r) strain and 671 DEGs in the 8y strain. The DEGs involved in resistance to fenpropathrin were further identified by comparing the two strains, including 207 upregulated DEGs in p50(r) and 175 downregulated DEGs in 8y. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that these fenpropathrin-related DEGs are mainly enriched in the metabolism and transporter pathways. Moreover, 28 DEGs involved in the metabolic pathway and 18 in the transporter pathway were identified. Furthermore, organic cation transporter protein 6 (BmOCT6), a transporter pathway member, was crucial in enhancing the tolerance of BmN cells to fenpropathrin. Finally, the knockdown of the expression of the homologs of BmOCT6 in Glyphodes pyloalis (G. pyloalis) significantly decreased the resistant level of larvae to fenpropathrin. The findings showed that the metabolism and transporter pathways are associated with resistance to fenpropathrin in silkworm, and OCT6 is an effective and potential target not only for silkworm breeding but also for pest biocontrol.

Transcriptome Analysis for Identification of Genes Related to Growth and Development, Digestion and Detoxification, Olfaction in the Litchi Stink Bug Tessaratoma papillosa

Tessaratoma papillosa is a major pest of Litchi chinensis and Dimocarpus longan. Adult and nymph secretions are not only harmful to plants but also to humans. At present, there are not a lot of research on T. papillosa, especially omics research. We used high-throughput sequencing technology to sequence the T. papillosa transcriptome and obtained 67,597 unigenes homologous to Halyomorpha halys (88.03%). Subsequently, RNA-SEQ and comparative analyses were performed on the 14 different developmental stages and tissues. A total of 462 unigenes related to growth and development, 1,851 unigenes related to digestion and detoxification, and 70 unigenes related to olfaction were obtained. Moreover, expression analysis showed that the T. papillosa major life activities genes are uniformly expressed across all developmental states. However, the adult midgut gene expression patterns were utterly different from that of the nymphs. Similarly, female fat body genes exhibited distinct expression patterns compared to that of males and nymphs. Thus, different developmental stages and physiological functions affect gene expression patterns. We also found that most of the differential genes were associated with cellular maintenance. This study will help understand the growth and development of litchi stink bugs, their choice of host plants, food digestion and detoxification, and their reproductive behavior. In addition, this result can provide reference information for some target genes in the process of control of T. papillosa.

Selective Sweeps in a Nutshell: The Genomic Footprint of Rapid Insecticide Resistance Evolution in the Almond Agroecosystem

Among the most familiar forms of human-driven evolution on ecological time scales is the rapid acquisition of resistance to pesticides by insects. Since the widespread adoption of synthetic organic insecticides in the mid-twentieth century, over 500 arthropod species have evolved resistance to at least one insecticide. Efforts to determine the genetic bases of insecticide resistance have historically focused on individual loci, but the availability of genomic tools has facilitated the screening of genome-wide characteristics. We resequenced three contemporary populations of the navel orangeworm (Amyelois transitella), the principal pest of almond orchards in California, differing in bifenthrin resistance status to examine insecticide-induced changes in the population genomic landscape of this species. We detected an exceptionally large region with virtually no polymorphisms, extending to up to 1.3 Mb in the resistant population. This selective sweep includes genes associated with pyrethroid and DDT resistance, including a cytochrome P450 gene cluster and the gene encoding the voltage-gated sodium channel para. Moreover, the sequence along the sweep is nearly identical in the genome assembled from a population founded in 1966, suggesting that the foundation for insecticide resistance may date back a half-century, when California’s Central Valley experienced massive area-wide applications of DDT for pest control.

Comparative CYP-omic analysis between the DDT-susceptible and -resistant Drosophila melanogaster strains 91-C and 91-R

BACKGROUND

Cytochrome P450 monooxygenases (P450s) are involved in the biosynthesis of endogenous intracellular compounds and the metabolism of xenobiotics, including chemical insecticides. We investigated the structural and expression level variance across all P450 genes with respect to the evolution of insecticide resistance under multigenerational dichlorodiphenyltrichloroethane (DDT) selection.

RESULTS

RNA-sequencing (RNA-seq) and reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) indicated that the transcript levels of seven P450 genes were significantly up-regulated and three P450 genes were down-regulated in the DDT-resistant strain 91-R, as compared to the control strain 91-C. The overexpression of Cyp6g1 was associated with the presence of an Accord and an HMS-Beagle element insertion in the 5' upstream region in conjunction with copy number variation in the 91-R strain, but not in the 91-C strain. A total of 122 (50.2%) fixed nonsynonymous (amino acid-changing) mutations were found between 91-C and 91-R, and 20 (8.2%) resulted in amino acid changes within functional domains. Three P450 proteins were truncated as a result of premature stop codons and fixed between strains.

CONCLUSION

Our results demonstrate that a combination of changes in P450 protein-coding regions and transcript levels are possibly associated with DDT resistance, and thereby suggest that selection for variant function may occur within this gene family in response to chronic DDT exposure. © 2018 Society of Chemical Industry.

The unusual reproductive system of head and body lice (Pediculus humanus)

Insect reproduction is extremely variable, but the implications of alternative genetic systems are often overlooked in studies on the evolution of insecticide resistance. Both ecotypes of Pediculus humanus (Phthiraptera: Pediculidae), the human head and body lice, are human ectoparasites, the control of which is challenged by the recent spread of resistance alleles. The present study conclusively establishes for the first time that both head and body lice reproduce through paternal genome elimination (PGE), an unusual genetic system in which males transmit only their maternally derived chromosomes. Here, we investigate inheritance patterns of parental genomes using a genotyping approach across families of both ecotypes and show that heterozygous males exclusively or preferentially pass on one allele only, whereas females transmit both in a Mendelian fashion. We do however observe occasional transmission of paternal chromosomes through males, representing the first known case of PGE in which whole-genome meiotic drive is incomplete. Finally, we discuss the potential implications of this finding for the evolution of resistance and invite the development of new theoretical models of how this knowledge might contribute to increasing the success of pediculicide-based management schemes.

DIRProt: a computational approach for discriminating insecticide resistant proteins from non-resistant proteins

BACKGROUND

Insecticide resistance is a major challenge for the control program of insect pests in the fields of crop protection, human and animal health etc. Resistance to different insecticides is conferred by the proteins encoded from certain class of genes of the insects. To distinguish the insecticide resistant proteins from non-resistant proteins, no computational tool is available till date. Thus, development of such a computational tool will be helpful in predicting the insecticide resistant proteins, which can be targeted for developing appropriate insecticides.

RESULTS

Five different sets of feature viz., amino acid composition (AAC), di-peptide composition (DPC), pseudo amino acid composition (PAAC), composition-transition-distribution (CTD) and auto-correlation function (ACF) were used to map the protein sequences into numeric feature vectors. The encoded numeric vectors were then used as input in support vector machine (SVM) for classification of insecticide resistant and non-resistant proteins. Higher accuracies were obtained under RBF kernel than that of other kernels. Further, accuracies were observed to be higher for DPC feature set as compared to others. The proposed approach achieved an overall accuracy of >90% in discriminating resistant from non-resistant proteins. Further, the two classes of resistant proteins i.e., detoxification-based and target-based were discriminated from non-resistant proteins with >95% accuracy. Besides, >95% accuracy was also observed for discrimination of proteins involved in detoxification- and target-based resistance mechanisms. The proposed approach not only outperformed Blastp, PSI-Blast and Delta-Blast algorithms, but also achieved >92% accuracy while assessed using an independent dataset of 75 insecticide resistant proteins.

CONCLUSIONS

This paper presents the first computational approach for discriminating the insecticide resistant proteins from non-resistant proteins. Based on the proposed approach, an online prediction server DIRProt has also been developed for computational prediction of insecticide resistant proteins, which is accessible at http://cabgrid.res.in:8080/dirprot/ . The proposed approach is believed to supplement the efforts needed to develop dynamic insecticides in wet-lab by targeting the insecticide resistant proteins.

Quantitative trait locus mapping and functional genomics of an organophosphate resistance trait in the western corn rootworm, Diabrotica virgifera virgifera

The western corn rootworm, Diabrotica virgifera virgifera, is an insect pest of corn and population suppression with chemical insecticides is an important management tool. Traits conferring organophosphate insecticide resistance have increased in frequency amongst D. v. virgifera populations, resulting in the reduced efficacy in many corn-growing regions of the USA. We used comparative functional genomic and quantitative trait locus (QTL) mapping approaches to investigate the genetic basis of D. v. virgifera resistance to the organophosphate methyl-parathion. RNA from adult methyl-parathion resistant and susceptible adults was hybridized to 8331 microarray probes. The results predicted that 11 transcripts were significantly up-regulated in resistant phenotypes, with the most significant (fold increases ≥ 2.43) being an α-esterase-like transcript. Differential expression was validated only for the α-esterase (ST020027A20C03), with 11- to 13-fold greater expression in methyl-parathion resistant adults (P < 0.05). Progeny with a segregating methyl-parathion resistance trait were obtained from a reciprocal backcross design. QTL analyses of high-throughput single nucleotide polymorphism genotype data predicted involvement of a single genome interval. These data suggest that a specific carboyxesterase may function in field-evolved corn rootworm resistance to organophosphates, even though direct linkage between the QTL and this locus could not be established.

Identification of a major Quantitative Trait Locus determining resistance to the organophosphate temephos in the dengue vector mosquito Aedes aegypti

Organophosphate insecticides (OP) have extensively been used to control mosquitoes, such as the vector Aedes aegypti. Unfortunately, OP resistance has hampered control programs worldwide. We used Quantitative Trait Locus (QTL) mapping to evaluate temephos resistance in two F1 intercross populations derived from crosses between a resistant Ae. aegypti strain (RecR) and two susceptible strains (MoyoD and Red). A single major effect QTL was identified on chromosome 2 of both segregating populations, named rtt1 (resistance to temephos 1). Bioinformatics analyses identified a cluster of carboxylesterase genes (CCE) within the rtt1 interval. qRT-PCR demonstrated that different CCEs were up-regulated in F2 resistant individuals from both crosses. However, none exceeded the 2-fold expression. Primary mechanisms for temephos resistance may vary between Ae. aegypti populations, yet also appear to support previous findings suggesting that multiple linked esterase genes may contribute to temephos resistance in the RecR strain as well as other populations.

Draft genome of the most devastating insect pest of coffee worldwide: the coffee berry borer, Hypothenemus hampei

The coffee berry borer, Hypothenemus hampei, is the most economically important insect pest of coffee worldwide. We present an analysis of the draft genome of the coffee berry borer, the third genome for a Coleopteran species. The genome size is ca. 163 Mb with 19,222 predicted protein-coding genes. Analysis was focused on genes involved in primary digestion as well as gene families involved in detoxification of plant defense molecules and insecticides, such as carboxylesterases, cytochrome P450, gluthathione S-transferases, ATP-binding cassette transporters, and a gene that confers resistance to the insecticide dieldrin. A broad range of enzymes capable of degrading complex polysaccharides were identified. We also evaluated the pathogen defense system and found homologs to antimicrobial genes reported in the Drosophila genome. Ten cases of horizontal gene transfer were identified with evidence for expression, integration into the H. hampei genome, and phylogenetic evidence that the sequences are more closely related to bacterial rather than eukaryotic genes. The draft genome analysis broadly expands our knowledge on the biology of a devastating tropical insect pest and suggests new pest management strategies.

A glycine insertion in the estrogen-related receptor (ERR) is associated with enhanced expression of three cytochrome P450 genes in transgenic Drosophila melanogaster

Insecticide-resistant Drosophila melanogaster strains represent a resource for the discovery of the underlying molecular mechanisms of cytochrome P450 constitutive over-expression, even if some of these P450s are not directly involved in the resistance phenotype. For example, in select 4,4'-dichlorodiphenyltrichloroethane (DDT) resistant strains the glucocorticoid receptor-like (GR-like) potential transcription factor binding motifs (TFBMs) have previously been shown to be associated with constitutively differentially-expressed cytochrome P450s, Cyp12d1, Cyp6g2 and Cyp9c1. However, insects are not known to have glucocorticoids. The only ortholog to the mammalian glucocorticoid receptor (GR) in D. melanogaster is an estrogen-related receptor (ERR) gene, which has two predicted alternative splice isoforms (ERRa and ERRb). Sequencing of ERRa and ERRb in select DDT susceptible and resistant D. melanogaster strains has revealed a glycine (G) codon insertion which was only observed in the ligand binding domain of ERR from the resistant strains tested (ERR-G). Transgenic flies, expressing the ERRa-G allele, constitutively over-expressed Cyp12d1, Cyp6g2 and Cyp9c1. Only Cyp12d1 and Cyp6g2 were over-expressed in the ERRb-G transgenic flies. Phylogenetic studies show that the G-insertion appeared to be located in a less conserved domain in ERR and this insertion is found in multiple species across the Sophophora subgenera.

Organophosphate and pyrethroid hydrolase activities of mutant Esterases from the cotton bollworm Helicoverpa armigera

Two mutations have been found in five closely related insect esterases (from four higher Diptera and a hymenopteran) which each confer organophosphate (OP) hydrolase activity on the enzyme and OP resistance on the insect. One mutation converts a Glycine to an Aspartate, and the other converts a Tryptophan to a Leucine in the enzymes’ active site. One of the dipteran enzymes with the Leucine mutation also shows enhanced activity against pyrethroids. Introduction of the two mutations in vitro into eight esterases from six other widely separated insect groups has also been reported to increase substantially the OP hydrolase activity of most of them. These data suggest that the two mutations could contribute to OP, and possibly pyrethroid, resistance in a variety of insects. We therefore introduced them in vitro into eight Helicoverpa armigera esterases from a clade that has already been implicated in OP and pyrethroid resistance. We found that they do not generally enhance either OP or pyrethroid hydrolysis in these esterases but the Aspartate mutation did increase OP hydrolysis in one enzyme by about 14 fold and the Leucine mutation caused a 4–6 fold increase in activity (more in one case) of another three against some of the most insecticidal isomers of fenvalerate and cypermethrin. The Aspartate enzyme and one of the Leucine enzymes occur in regions of the H. armigera esterase isozyme profile that have been previously implicated in OP and pyrethroid resistance, respectively.

Susceptibility of Rhipicephalus(Boophilus) microplus (Acari: Ixodidae) to pyrethroids and their associations in Pernambuco, Brazil

The synthetic pyrethroids and their associations have been widely used for controlling Rhipicephalus (Boophilus)microplus. The frequent use of acaricides has been inducing the development of resistance in the tick populations. The aim of this study was to assess the susceptibility of R. (B.)microplus populations to pyrethroids and their associations in the region of Garanhuns, Pernambuco, Brazil. In addition, the level of information among farm owners regarding tick control measures was investigated. Ticks were collected directly from naturally infested dairy cattle in the region and were exposed to pyrethroids and their associations. At the same time, an epidemiological questionnaire was applied with the aim of investigating the level of information among the farmers. The results reported here indicate thatR. (B.) micropluspopulations in the dairy region of Garanhuns show resistance to pyrethroids and their associations, except when the product is associated with piperonyl butoxide. Regarding the results from the epidemiological survey, it was seen that there is a considerable lack of information among the farmers in relation to ixodid control measures. The level of ticks resistance to acaricides varied widely across the region studied. No alternative control programs have been implemented among these farms, thus demonstrating that there is a need for more information relating to the biology and control of R.(B.) microplus.

Transcriptomic and proteomic responses of sweetpotato whitefly, Bemisia tabaci, to thiamethoxam

BACKGROUND

The sweetpotato whitefly, Bemisia tabaci (Hemiptera: Aleyrodidae), is one of the most widely distributed agricultural pests. Although it has developed resistance to many registered insecticides including the neonicotinoid insecticide thiamethoxam, the mechanisms that regulate the resistance are poorly understood. To understand the molecular basis of thiamethoxam resistance, "omics" analyses were carried out to examine differences between resistant and susceptible B. tabaci at both transcriptional and translational levels.

RESULTS

A total of 1,338 mRNAs and 52 proteins were differentially expressed between resistant and susceptible B. tabaci. Among them, 11 transcripts had concurrent transcription and translation profiles. KEGG analysis mapped 318 and 35 differentially expressed genes and proteins, respectively, to 160 and 59 pathways (p<0.05). Thiamethoxam treatment activated metabolic pathways (e.g., drug metabolism), in which 118 transcripts were putatively linked to insecticide resistance, including up-regulated glutathione-S-transferase, UDP glucuronosyltransferase, glucosyl/glucuronosyl transferase, and cytochrome P450. Gene Ontology analysis placed these genes and proteins into protein complex, metabolic process, cellular process, signaling, and response to stimulus categories. Quantitative real-time PCR analysis validated "omics" response, and suggested a highly overexpressed P450, CYP6CX1, as a candidate molecular basis for the mechanistic study of thiamethoxam resistance in whiteflies. Finally, enzymatic activity assays showed elevated detoxification activities in the resistant B. tabaci.

CONCLUSIONS

This study demonstrates the applicability of high-throughput omics tools for identifying molecular candidates related to thiamethoxam resistance in an agricultural important insect pest. In addition, transcriptomic and proteomic analyses provide a solid foundation for future functional investigations into the complex molecular mechanisms governing the neonicotinoid resistance in whiteflies.

Nuclear genetic diversity in human lice (Pediculus humanus) reveals continental differences and high inbreeding among worldwide populations

Understanding the evolution of parasites is important to both basic and applied evolutionary biology. Knowledge of the genetic structure of parasite populations is critical for our ability to predict how an infection can spread through a host population and for the design of effective control methods. However, very little is known about the genetic structure of most human parasites, including the human louse (Pediculus humanus). This species is composed of two ecotypes: the head louse (Pediculus humanus capitis De Geer), and the clothing (body) louse (Pediculus humanus humanus Linnaeus). Hundreds of millions of head louse infestations affect children every year, and this number is on the rise, in part because of increased resistance to insecticides. Clothing lice affect mostly homeless and refugee-camp populations and although they are less prevalent than head lice, the medical consequences are more severe because they vector deadly bacterial pathogens. In this study we present the first assessment of the genetic structure of human louse populations by analyzing the nuclear genetic variation at 15 newly developed microsatellite loci in 93 human lice from 11 sites in four world regions. Both ecotypes showed heterozygote deficits relative to Hardy-Weinberg equilibrium and high inbreeding values, an expected pattern given their parasitic life history. Bayesian clustering analyses assigned lice to four distinct genetic clusters that were geographically structured. The low levels of gene flow among louse populations suggested that the evolution of insecticide resistance in lice would most likely be affected by local selection pressures, underscoring the importance of tailoring control strategies to population-specific genetic makeup and evolutionary history. Our panel of microsatellite markers provides powerful data to investigate not only ecological and evolutionary processes in lice, but also those in their human hosts because of the long-term coevolutionary association between lice and humans.

Use of mutagenesis, genetic mapping and next generation transcriptomics to investigate insecticide resistance mechanisms

Insecticide resistance is a worldwide problem with major impact on agriculture and human health. Understanding the underlying molecular mechanisms is crucial for the management of the phenomenon; however, this information often comes late with respect to the implementation of efficient counter-measures, particularly in the case of metabolism-based resistance mechanisms. We employed a genome-wide insertional mutagenesis screen to Drosophila melanogaster, using a Minos-based construct, and retrieved a line (MiT[w⁻]3R2) resistant to the neonicotinoid insecticide Imidacloprid. Biochemical and bioassay data indicated that resistance was due to increased P450 detoxification. Deep sequencing transcriptomic analysis revealed substantial over- and under-representation of 357 transcripts in the resistant line, including statistically significant changes in mixed function oxidases, peptidases and cuticular proteins. Three P450 genes (Cyp4p2, Cyp6a2 and Cyp6g1) located on the 2R chromosome, are highly up-regulated in mutant flies compared to susceptible Drosophila. One of them (Cyp6g1) has been already described as a major factor for Imidacloprid resistance, which validated the approach. Elevated expression of the Cyp4p2 was not previously documented in Drosophila lines resistant to neonicotinoids. In silico analysis using the Drosophila reference genome failed to detect transcription binding factors or microRNAs associated with the over-expressed Cyp genes. The resistant line did not contain a Minos insertion in its chromosomes, suggesting a hit-and-run event, i.e. an insertion of the transposable element, followed by an excision which caused the mutation. Genetic mapping placed the resistance locus to the right arm of the second chromosome, within a ∼1 Mb region, where the highly up-regulated Cyp6g1 gene is located. The nature of the unknown mutation that causes resistance is discussed on the basis of these results.

Insecticide resistance mechanisms in the green peach aphid Myzus persicae (Hemiptera: Aphididae) I: A transcriptomic survey

BACKGROUND

Insecticide resistance is one of the best examples of rapid micro-evolution found in nature. Since the development of the first synthetic insecticide in 1939, humans have invested considerable effort to stay ahead of resistance phenotypes that repeatedly develop in insects. Aphids are a group of insects that have become global pests in agriculture and frequently exhibit insecticide resistance. The green peach aphid, Myzus persicae, has developed resistance to at least seventy different synthetic compounds, and different insecticide resistance mechanisms have been reported worldwide.

METHODOLOGY/PRINCIPAL FINDINGS

To further characterize this resistance, we analyzed genome-wide transcriptional responses in three genotypes of M. persicae, each exhibiting different resistance mechanisms, in response to an anti-cholinesterase insecticide. The sensitive genotype (exhibiting no resistance mechanism) responded to the insecticide by up-regulating 183 genes primarily ones related to energy metabolism, detoxifying enzymes, proteins of extracellular transport, peptidases and cuticular proteins. The second genotype (resistant through a kdr sodium channel mutation), up-regulated 17 genes coding for detoxifying enzymes, peptidase and cuticular proteins. Finally, a multiply resistant genotype (carrying kdr and a modified acetylcholinesterase), up-regulated only 7 genes, appears not to require induced insecticide detoxification, and instead down-regulated many genes.

CONCLUSIONS/SIGNIFICANCE

This study suggests strongly that insecticide resistance in M. persicae is more complex that has been described, with the participation of a broad array of resistance mechanisms. The sensitive genotype exhibited the highest transcriptional plasticity, accounting for the wide range of potential adaptations to insecticides that this species can evolve. In contrast, the multiply resistant genotype exhibited a low transcriptional plasticity, even for the expression of genes encoding enzymes involved in insecticide detoxification. Our results emphasize the value of microarray studies to search for regulated genes in insects, but also highlights the many ways those different genotypes can assemble resistant phenotypes depending on the environmental pressure.

Laboratory and field comparisons of insecticides to reduce infestation of Drosophila suzukii in berry crops

BACKGROUND

The spotted wing Drosophila, Drosophila suzukii (Matsumura, 1931) (Diptera: Drosophilidae), is an invasive pest of small-fruit crops. Unlike most other Drosophila, this insect is able to oviposit into and damage ripe and ripening fruit, making it unmarketable. Because this is a new pest in the United States, it is necessary to identify registered insecticides to manage this insect effectively in conventional and organic production systems.

RESULTS

The present laboratory bioassays and field trials identified a number of insecticides representing various modes of action that are effective in controlling D. suzukii. Products that performed well in the laboratory bioassay also performed well in the field, indicating that screening of new chemistries in the laboratory is a worthy exercise. Field application of pyrethoids, organophosphates or spinosyns provided 5-14 days of residual control of D. suzukii. The efficacy of the neonicotinoids as adulticides was not satisfactory compared with the other contact-mode-of-action chemistries. Based on the zero tolerance by the small-fruit industry and the individual effects mentioned above, neonicotinoids are not currently recommended for D. suzukii management.

CONCLUSIONS

There are effective insecticides registered for controlling D. suzukii infestations in susceptible small-fruit crops.

In silico characterization and comparative genomic analysis of the Culex quinquefasciatus glutathione S-transferase (GST) supergene family

The glutathione S-transferases (GSTs) are phase II class of detoxification enzymes that are involved both directly and indirectly in insecticide resistance mechanisms. The Culex quinquefasciatus GST superfamily was analyzed by utilizing the public domain Culex genome sequence. In total, 35 cytosolic (seven classes) and 5 microsomal putatively active GSTs were retrieved, classified, and annotated. The study revealed the presence of three unclassified GSTs. Of 35 cytosolic GSTs, 65% contributed by insect specific Delta-Epsilon classes. Gene cluster analysis revealed that most of the genes of Delta, Epsilon, and Theta classes were organized into gene clusters. The gene organization analysis revealed the dominance of phase "0" introns in the Culex GST family. The studies on intron loss and gain events revealed that the Delta GSTs have experienced a higher number of loss and gains during their evolution. A positive correlation was observed between the phylogenetic relationship of members of the GST superfamily and their corresponding exon-intron organization. In addition, the genes within the gene clusters revealed the monophyletic phylogenetic relationship implying the importance of gene duplication events in the gene families' evolution. Finally, the comparative genomic analysis has shown a complex evolutionary scenario associated with the GST supergene family evolution in insects.

Selective sweeps at the organophosphorus insecticide resistance locus, Rop-1, have affected variation across and beyond the α-esterase gene cluster in the Australian sheep blowfly, Lucilia cuprina

A major theoretical consequence of selection at a locus is the genetic hitchhiking of linked sites (selective sweep). The extent of hitchhiking around a gene is related to the strength of selection and the rate of recombination, with its impact diminishing with distance from the selected site. At the Rop-1 locus of the sheep blowfly, Lucilia cuprina, polymorphisms at two different sites within the LcαE7 gene encode forms of the protein that confer organophosphorus insecticide resistance. To assess the impact of selection at these two sites on variation around LcαE7, we sequenced regions within six other genes along chromosome IV across isogenic (IV) strains of L. cuprina. High levels of linkage disequilibrium, characterized by low haplotype number (K) and diversity (H), and significant R(2) values were observed for two genes, LcαE1 and LcαE10, both members of the same α-esterase gene cluster as LcαE7. A significant R(2) value was also observed for a gene predicted to be the next closest to LcαE7, AL03, but not for any of the other genes, LcRpL13a, Lcdsx, or LcAce. Skews in the site frequency spectra toward high-frequency variants were significant for LcαE1 (Fay and Wu's H = -2.91), LcαE10 (H = -1.85), and Lcdsx (H = -2.00). Since the selective sweeps, two forms of likely returning variation were observed, including variation in microsatellites in an intron of LcαE10 and a recombination event between LcαE7 and LcαE10. These data suggest that two incomplete soft sweeps have occurred at LcαE7 that have significantly affected variation across, and beyond, the α-esterase gene cluster of L. cuprina. The speed and impact of these selective sweeps on surrounding genomic variation and the ability of L. cuprina to respond to future environmental challenges are discussed.

Quantitative trait loci mapping of genome regions controlling permethrin resistance in the mosquito Aedes aegypti

The mosquito Aedes aegypti is the principal vector of dengue and yellow fever flaviviruses. Permethrin is an insecticide used to suppress Ae. aegypti adult populations but metabolic and target site resistance to pyrethroids has evolved in many locations worldwide. Quantitative trait loci (QTL) controlling permethrin survival in Ae. aegypti were mapped in an F(3) advanced intercross line. Parents came from a collection of mosquitoes from Isla Mujeres, México, that had been selected for permethrin resistance for two generations and a reference permethrin-susceptible strain originally from New Orleans. Following a 1-hr permethrin exposure, 439 F(3) adult mosquitoes were phenotyped as knockdown resistant, knocked down/recovered, or dead. For QTL mapping, single nucleotide polymorphisms (SNPs) were identified at 22 loci with potential antixenobiotic activity including genes encoding cytochrome P450s (CYP), esterases (EST), or glutathione transferases (GST) and at 12 previously mapped loci. Seven antixenobiotic genes mapped to chromosome I, six to chromosome II, and nine to chromosome III. Two QTL of major effect were detected on chromosome III. One corresponds with a SNP previously associated with permethrin resistance in the para sodium channel gene and the second with the CCEunk7o esterase marker. Additional QTL but of relatively minor effect were also found. These included two sex-linked QTL on chromosome I affecting knockdown and recovery and a QTL affecting survival and recovery. On chromosome II, one QTL affecting survival and a second affecting recovery were detected. The patterns confirm that mutations in the para gene cause target-site insensitivity and are the major source of permethrin resistance but that other genes dispersed throughout the genome contribute to recovery and survival of mosquitoes following permethrin exposure.

Mity model: Tetranychus urticae, a candidate for chelicerate model organism

Chelicerates (scorpions, horseshoe crabs, spiders, mites and ticks) are the second largest group of arthropods and are of immense importance for fundamental and applied science. They occupy a basal phylogenetic position within the phylum Arthropoda, and are of crucial significance for understanding the evolution of various arthropod lineages. Chelicerates are vectors of human diseases, such as ticks, and major agricultural pests, such as spider mites, thus this group is also of importance for both medicine and agriculture. The developmental genetics of chelicerates is poorly understood and a challenge for the future progress for many aspects of chelicerate biology is the development of a model organism for this group. Toward this end, we are developing a chelicerate genetic model: the two-spotted spider mite Tetranychus urticae. T. urticae has the smallest genome of any arthropod determined so far (75 Mbp, 60% of the size of the Drosophila genome), undergoes rapid development and is easy to maintain in the laboratory. These features make T. urticae a promising reference organism for the economically important, poorly studied and species-rich chelicerate lineage.

Contributions of genomics to life-history theory

Life-history theory seeks to understand the factors that produce variation in life histories that are found both among and within species. At the organismal level there is a well developed mathematical framework, and an important focus of the current research is determining the biological underpinnings of this framework, with particular attention to the causal mechanisms that underlie trade-offs. Genomic approaches are proving useful in addressing this issue.

Temporal synergism by microencapsulation of piperonyl butoxide and alpha-cypermethrin overcomes insecticide resistance in crop pests

A microencapsulated formulation that gives a burst release of piperonyl butoxide (PBO) several hours before a burst release of a conventional pyrethroid can effectively overcome metabolic resistance in Bemisia tabaci Gennadius, Helicoverpa armigera (Hübner), Aphis gossypii Glover and Myzus persicae Sulzer. This increase in efficacy against resistant pests was reflected in a field trial against B. tabaci on cotton, eliminating the need for two treatments. The ratio between the active insecticide and the synergist was found to be crucial in reducing resistance factors.

Sharing of transcription factors after gene duplication in the yeast Saccharomyces cerevisiae

In a set of 190 duplicate gene pairs in yeast Saccharomyces cerevisiae, the sharing of transcription factors tended to decrease with increased divergence in coding sequence, at both synonymous and nonsynonymous sites. Our results showed a significantly higher sharing of transcription factors by duplicated gene pairs falling within duplicated genomic blocks than in other duplicated gene pairs; and genes in duplicated blocks also showed significantly greater conservation at the coding sequence level. In spite of the overall trends, there were certain gene pairs, both in duplicated blocks and in other genomic regions, which were highly divergent in coding sequence and yet had identical patterns of transcription factor binding. These results suggest that functional differentiation of genes after duplication is a multi-dimensional process, with different duplicate pairs differentiating in different ways.

Seasonal fluctuation in susceptibility to insecticides within natural populations of Drosophila melanogaster. II. Features of genetic variation in susceptibility to organophosphate insecticides within natural populations of D. melanogaster

To elucidate genetic variation in susceptibility to organophosphate insecticides within natural populations of Drosophila melanogaster, we conducted an analysis of variance for mortality data sets of isofemale lines (10-286 lines) used in the previous studies. Susceptibility of isofemale lines to the three organophosphate insecticides was continuously distributed within each natural population, ranging from susceptible to resistant. Analysis of variance showed highly significant variation among isofemale lines in susceptibility to each insecticide for each natural population. Significant genetic variances in susceptibility to the three chemicals were estimated for the Katsunuma population; 0.0529-0.2722 for malathion, 0.0492-0.1603 for prothiophos, and 0.0469-0.1696 for fenitrothion. Contrary to the consistent seasonal tendency towards an increase in mean susceptibility in the fall, reported in the previous study, genetic variances in susceptibility to the three organophosphates did not change significantly in 1997 but tended to increase by 2- to 5-times in 1998. We tested whether both the observed situations, maintenance and increase in genetic variance in organophosphate resistance, can be generated under circumstances in which the levels of resistance to the three organophosphates tended to decrease, by conducting a simulation analysis, based on the hypothesis that resistant genotypes have lower fitnesses than susceptible ones under the density-independent condition. The simulation analysis generally explained the pattern in the mean susceptibility and genetic variances in susceptibility to the three organophosphates, observed in the Katsunuma population of D. melanogaster. It was suggested that the differences in the frequencies of resistance genes in the summer population could affect the patterns in genetic variance in organophosphate resistance in the fall population.

Genome-wide transcription profile of field- and laboratory-selected dichlorodiphenyltrichloroethane (DDT)-resistant Drosophila

Genome-wide microarray analysis (Affymetrix array) was used (i) to determine whether only one gene, the cytochrome P450 enzyme Cyp6g1, is differentially transcribed in dichlorodiphenyltrichloroethane (DDT)-resistant vs. -susceptible Drosophila; and (ii) to profile common genes differentially transcribed across a DDT-resistant field isolate [Rst(2)DDT(Wisconsin)] and a laboratory DDT-selected population [Rst(2)DDT(91-R)]. Statistical analysis (ANOVA model) identified 158 probe sets that were differentially transcribed among Rst(2)DDT(91-R), Rst(2)DDT(Wisconsin), and the DDT-susceptible genotype Canton-S (P < 0.01). The cytochrome P450 Cyp6a2 and the diazepam-binding inhibitor gene (Dbi) were over transcribed in the two DDT-resistant genotypes when compared to the wild-type Drosophila, and this difference was significant at the most stringent statistical level, a Bonferroni correction. The list of potential candidates differentially transcribed also includes 63 probe sets for which molecular function ontology annotation of the probe sets did not exist. A total of four genes (Cyp6a2, Dbi, Uhg1, and CG11176) were significantly different (P < 5.6 e(-06)) between Rst(2)DDT(91-R) and Canton-S. Additionally, two probe sets encoding Cyp12d1 and Dbi were significantly different between Rst(2)DDT(Wisconsin) and Canton-S after a Bonferroni correction. Fifty-two probe sets, including those associated with pesticide detoxification, ion transport, signal transduction, RNA transcription, and lipid metabolism, were commonly expressed in both resistant lines but were differentially transcribed in Canton-S. Our results suggest that more than Cyp6g1 is overtranscribed in field and laboratory DDT-resistant genotypes, and the number of commonalities suggests that similar resistance mechanisms may exist between laboratory- and field-selected DDT-resistant fly lines.

Microarray analysis of cytochrome P450 mediated insecticide resistance in Drosophila

Insecticide resistance in laboratory selected Drosophila strains has been associated with upregulation of a range of different cytochrome P450s, however in recent field isolates of D. melanogaster resistance to DDT and other compounds is conferred by one P450 gene, Cyp6g1. Using microarray analysis of all Drosophila P450 genes, here we show that different P450 genes such as Cyp12d1 and Cyp6a8 can also be selected using DDT in the laboratory. We also show, however, that a homolog of Cyp6g1 is over-expressed in a field resistant strain of D. simulans. In order to determine why Cyp6g1 is so widely selected in the field we examine the pattern of cross-resistance of both resistant strains and transgenic flies over-expressing Cyp6g1 alone. We show that all three DDT selected P450s can confer resistance to the neonicotinoid imidacloprid but that Cyp6a8 confers no cross-resistance to malathion. Transgenic flies over-expressing Cyp6g1 also show cross-resistance to other neonicotinoids such as acetamiprid and nitenpyram. We suggest that the broad level of cross-resistance shown by Cyp6g1 may have facilitated its selection as a resistance gene in natural Drosophila populations.