Global transcriptome profiling and functional analysis reveal that tissue-specific constitutive overexpression of cytochrome P450s confers tolerance to imidacloprid in palm weevils in date palm fields

Deorphanizing an odorant receptor tuned to palm tree volatile esters in the Asian palm weevil sheds light on the mechanisms of palm tree selection

The Asian palm weevil, Rhynchophorus ferrugineus, is a tremendously important agricultural pest primarily adapted to palm trees and causes severe destruction, threatening sustainable palm cultivation worldwide. The host plant selection of this weevil is mainly attributed to the functional specialization of odorant receptors (ORs) that detect palm-derived volatiles. Yet, ligands are known for only two ORs of R. ferrugineus, and we still lack information on the mechanisms of palm tree detection. This study identified a highly expressed antennal R. ferrugineus OR, RferOR2, thanks to newly generated transcriptomic data. The phylogenetic analysis revealed that RferOR2 belongs to the major coleopteran OR group 2A and is closely related to a sister clade containing an R. ferrugineus OR (RferOR41) tuned to the non-host plant volatile and antagonist, α-pinene. Functional characterization of RferOR2 via heterologous expression in Drosophila olfactory neurons revealed that this receptor is tuned to several ecologically relevant palm-emitted odors, most notably ethyl and methyl ester compounds, but not to any of the pheromone compounds tested, including the R. ferrugineus aggregation pheromone. We did not evidence any differential expression of RferOR2 in the antennae of both sexes, suggesting males and females detect these compounds equally. Next, we used the newly identified RferOR2 ligands to demonstrate that including synthetic palm ester volatiles as single compounds and in combinations in pheromone-based mass trapping has a synergistic attractiveness effect to R. ferrugineus aggregation pheromone, resulting in significantly increased weevil catches. Our study identified a key OR from a palm weevil species tuned to several ecologically relevant palm volatiles and represents a significant step forward in understanding the chemosensory mechanisms of host detection in palm weevils. Our study also defines RferOR2 as an essential model for exploring the molecular basis of host detection in other palm weevil species. Finally, our work showed that insect OR deorphanization could aid in identifying novel behaviorally active volatiles that can interfere with weevil host-searching behavior in sustainable pest management applications.

Genome-wide identification, molecular evolution and gene expression of P450 gene family in Cyrtotrachelus buqueti

BACKGROUND

Insect Cytochrome P450 monooxygenase (CYPs or P450s) plays an important role in detoxifying insecticides, causing insect populations to develop resistance. However, the molecular functions of P450 gene family in Cyrtotrachelus buqueti genome are still lacking.

RESULTS

In this study, 71 CbuP450 genes have been identified. The amino acids length of CbuP450 proteins was between 183 aa ~ 1041 aa. They are proteins with transmembrane domains. The main component of their secondary structure is α-helix and random coils. Phylogenetic analysis showed that C. buqueti and Rhynchophorus ferrugineus were the most closely related. This gene family has 29 high-frequency codons, which tend to use A/T bases and A/T ending codons. Gene expression analysis showed that CbuP450_23 in the female adult may play an important role on high temperature resistance, and CbuP450_17 in the larval may play an important role on low temperature tolerance. CbuP450_10, CbuP450_17, CbuP450_23, CbuP450_10, CbuP450_16, CbuP450_20, CbuP450_23 and CbuP450_ 29 may be related to the regulation of bamboo fiber degradation genes in C. buqueti. Protein interaction analysis indicates that most CbuP450 proteins are mainly divided into three aspects: encoding the biosynthesis of ecdysteroids, participating in the decomposition of synthetic insecticides, metabolizing insect hormones, and participating in the detoxification of compounds.

CONCLUSIONS

We systematically analyzed the gene and protein characteristics, gene expression, and protein interactions of CbuP450 gene family, revealing the key genes involved in the stress response of CbuP450 gene family in the resistance of C. buqueti to high or low temperature stress, and identified the key CbuP450 proteins involved in important life activity metabolism. These results provided a reference for further research on the function of P450 gene family in C. buqueti.

Taxonomy, Biology, Symbionts, Omics, and Management of Rhynchophorus Palm Weevils (Coleoptera: Curculionidae: Dryophthorinae)

Palm weevils, Rhynchophorus spp., are destructive pests of native, ornamental, and agricultural palm species. Of the 10 recognized species, two of the most injurious species, Rhynchophorus ferrugineus and Rhynchophorus palmarum, both of which have spread beyond their native range, are the best studied. Due to its greater global spread and damage to edible date industries in the Middle East, R. ferrugineus has received more research interest. Integrated pest management programs utilize traps baited with aggregation pheromone, removal of infested palms, and insecticides. However, weevil control is costly, development of resistance to insecticides is problematic, and program efficacy can be impaired because early detection of infestations is difficult. The genome of R. ferrugineus has been sequenced, and omics research is providing insight into pheromone communication and changes in volatile and metabolism profiles of weevil-infested palms. We outline how such developments could lead to new control strategies and early detection tools.

Sustainable Pest Management in Date Palm Ecosystems: Unveiling the Ecological Dynamics of Red Palm Weevil (Coleoptera: Curculionidae) Infestations

The red palm weevil (RPW) poses a significant threat to date palm ecosystems, highlighting the necessity of sustainable pest management strategies that carefully consider the delicate ecological balance within these environments. This comprehensive review delves into innovative approaches to sustainable pest management, specifically focusing on date palm, and seeks to unravel the intricate ecological dynamics underlying RPW infestations. We thoroughly analyze biocontrol methods, eco-friendly chemical interventions, and integrated pest management (IPM) strategies, aiming to minimize the ecological impact while effectively addressing RPW infestations. By emphasizing the interplay of both living organisms (biotic) and environmental factors (abiotic) in shaping RPW dynamics, we advocate for a holistic and sustainable management approach that ensures the long-term resilience of date palm ecosystems. This review aims to contribute to an ecologically sound framework for pest management, promoting the sustainability and vitality of date palm ecosystems amidst the challenges posed by the RPW.

Identification and functional study of detoxification-related genes in response to tolfenpyrad stress in Glyphodes pyloalis Walker (Lepidoptera: Pyralidae)

Glyphodes pyloalis Walker (G. pyloalis) is a common destructive mulberry pest. Due to the long-term and frequent use of insecticides, it has developed tolerance to commonly used insecticides. Tolfenpyrad (TFP) is a novel pyrazole heterocyclic insecticide. In order to understand the TFP detoxification mechanism of G. pyloalis larvae, we first estimated the LC30 dose of TFP for 3rd instar G. pyloalis larvae. Next, we identified genes that were differentially expressed in 3rd instar G. pyloalis larvae treated with TFP compared to the control group by transcriptome sequencing. In total, 86,949,569 and 67,442,028 clean reads were obtained from TFP-treated and control G. pyloalis larvae, respectively. A total of 5588 differentially expressed genes (DEGs) were identified in TFP-treated and control G. pyloalis larvae, of which 3084 genes were upregulated and 2504 genes were downregulated. We analyzed the expression of 43 candidate detoxification enzyme genes associated with insecticide tolerance using qPCR. According to the spatiotemporal expression pattern of DEGs, we found that CYP6ABE1, CYP333A36 and GST-epsilon8 were highly expressed in the midgut, while CarEs14 was strongly expressed in haemolymph. Furthermore, we successfully knocked down these genes by RNA interference. After silencing CYP6ABE1 and CYP333A36, bioassay showed that the mortality rate of TFP-treated G. pyloalis larvae was significantly higher compared to the control group. This study provides a theoretical foundation for understanding the sensitivity of G. pyloalis to TFP and establish the basis for the effective and green management of this pest.

Genome-Wide Identification and Expression Analysis of the Cytochrome P450 Gene Family in Bemisia tabaci MED and Their Roles in the Insecticide Resistance

The whitefly, Bemisia tabaci MED (Hemiptera: Aleyrodidae), is an omnivorous agricultural pest, which causes huge economic losses to agriculture and is highly resistant to many pesticides. The overexpression of cytochrome P450 may play an important role in host adaptation and insecticide resistance in B. tabaci MED. Therefore, the present study systematically analyzed the cytochrome P450 gene family at the genome-wide level to understand its function in B. tabaci MED. Our analysis identified 58 cytochrome P450 genes in B. tabaci MED, among which 24 were novel. Phylogenetic analysis revealed broad functional and species-specific diversification in B. tabaci MED P450, suggesting the role of multiple P450 genes in detoxifying. Reverse transcription-real time quantitative PCR (RT-qPCR) showed that CYP4CS2, CYP4CS5, CYP4CS6, CYP4CS8, CYP6DW4, CYP6DW5, CYP6DW6, CYP6DZ8, and CYP6EN1 genes increased significantly after two days of exposure to imidacloprid. Interestingly, all nine genes belonged to the CYP4 and CYP6 families. A decrease in the expression of five genes (CYP6DW4, CYP6DW5, CYP6DW6, CYP6DZ8, and CYP4CS6) via RNA interference (RNAi) resulted in a significant increase in the mortalities of whiteflies when exposed to imidacloprid. These results indicate that the overexpression of the P450 genes may play an essential role in imidacloprid tolerance of B. tabaci MED. Thus, the present study provides basic information on P450 genes in B. tabaci MED, which will further help elucidate the insecticide resistance mechanism in the agricultural pest whitefly.

Omics in the Red Palm Weevil Rhynchophorus ferrugineus (Olivier) (Coleoptera: Curculionidae): A Bridge to the Pest

The red palm weevil (RPW), Rhynchophorus ferrugineus (Coleoptera: Curculionidae), is the most devastating pest of palm trees worldwide. Mitigation of the economic and biodiversity impact it causes is an international priority that could be greatly aided by a better understanding of its biology and genetics. Despite its relevance, the biology of the RPW remains poorly understood, and research on management strategies often focuses on outdated empirical methods that produce sub-optimal results. With the development of omics approaches in genetic research, new avenues for pest control are becoming increasingly feasible. For example, genetic engineering approaches become available once a species's target genes are well characterized in terms of their sequence, but also population variability, epistatic interactions, and more. In the last few years alone, there have been major advances in omics studies of the RPW. Multiple draft genomes are currently available, along with short and long-read transcriptomes, and metagenomes, which have facilitated the identification of genes of interest to the RPW scientific community. This review describes omics approaches previously applied to RPW research, highlights findings that could be impactful for pest management, and emphasizes future opportunities and challenges in this area of research.

Molecular Mechanisms Underlying Metabolic Resistance to Cyflumetofen and Bifenthrin in Tetranychus urticae Koch on Cowpea

Tetranychus urticae Koch (T. urticae) is one of the most tremendous herbivores due to its polyphagous characteristics, and is resistant to most acaricides. In this study, enzyme-linked immunosorbent assay (ELISA), transcriptome sequencing (RNA-seq) and quantitative real-time PCR (qRT-PCR) were carried out to analyze the mechanisms of T. urticae metabolic resistance to cyflumetofen and bifenthrin on cowpea. The enzyme activity of UDP-glucuronosyltransferases (UGTs) and carboxylesterases (CarEs) in the cyflumetofen-resistant (R_cfm) strain significantly decreased, while that of cytochrome P450 monooxygenases (P450s) significantly increased. Meanwhile, the activities of glutathione-S-transferases (GSTs), CarEs and P450s in the bifenthrin-resistant (R_bft) strain were significantly higher than those in the susceptible strain (Lab_SS). According to the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analyses, in the R_cfm mite strain, two carboxyl/cholinesterase (CCE) genes and two P450 genes were upregulated and one gene was downregulated, namely CYP392E7; in the R_bft mite strain, eleven CCE, nine UGT, two P450, four GST and three ABC genes were upregulated, while four CCE and three P450 genes were downregulated. Additionally, 94 differentially expressed genes (DEGs) were common to the two resistant groups. Specifically, TuCCE46 and TuCCE70 were upregulated in both resistant groups. Furthermore, the qRT-PCR validation data were consistent with those from the transcriptome sequencing analysis. Specifically, TuCCE46 (3.37-fold) was significantly upregulated in the R_cfm strain, while in the R_bft strain, TeturUGT22 (5.29-fold), teturUGT58p (1.74-fold), CYP392A11 (2.89-fold) and TuGSTd15 (5.12-fold) were significantly upregulated and TuCCE01 (0.13-fold) and CYP392A2p (0.07-fold) were significantly downregulated. Our study indicates that TuCCE46 might play the most important role in resistance to cyflumetofen, and TuCCE01, teturUGT58p, teturUGT22, CYP392A11, TuGSTd15, TuGSTm09 and TuABCG-13 were prominent in the resistance to bifenthrin. These findings provide further insight into the critical genes involved in the metabolic resistance of T. urticae to cyflumetofen and bifenthrin.

Gut transcriptome analysis of P450 genes and cytochrome P450 reductase in three moth species feeding on gymnosperms or angiosperms

Cytochrome P450 enzymes (P450s, CYPs) are a superfamily of heme–thiolate proteins involved in the metabolism of endogenous and exogenous substances in insects. In this study, the identification of putative P450 proteins was done and the elimination of the repeated sequences resulted in 57 proteins from Gastropacha populifolia, 63 proteins from Dendrolimus punctatus, and 53 proteins from Dendrolimus tabulaeformis. The putative P450 proteins were aligned together with seven other insect species based on five conserved domains. A total of ten co-orthologous groups were identified. Interestingly, one co-orthologous gene, CYP4g15 in CYP4 clan, was identified and its 3D structure analysis showed that the highly conserved sites of the predicted motifs were close to the active sites of P450. Furthermore, this study revealed that insect CYP4g15 and two bacteria cytochrome P450 were monophyletic. This suggests that insects CYP4g15 are not only functionally conserved but also an ancient gene originating from different bacteria species.

Genome-wide and expression-profiling analyses of the cytochrome P450 genes in Tenebrionidea

Cytochrome P450 monooxygenases (CYPs) are present in almost all areas of the tree of life. As one of the largest and most diverse superfamilies of multifunctional enzymes, they play important roles in the metabolism of xenobiotics and biosynthesis of endogenous compounds, shaping the success of insects. In this study, the CYPome (an omics term for all the CYP genes in a genome) diversification was examined in the four Tenebrionidea species through genome-wide analysis. A total of 483 CYP genes were identified, of which 103, 157, 122, and 101 were respectively deciphered from the genomes of Tebebrio molitor, Asbolus verucosus, Hycleus cichorii and Hycleus phaleratus. These CYPs were classified into four major clans (mitochondrial, CYP2, CYP3, and CYP4), and clans CYP3 and CYP4 are most diverse. Phylogenetic analysis showed that most CYPs of these Tenebrionidea beetles from each clan had a very close 1:1 orthology to each other, suggesting that they originate closely and have evolutionally conserved function. Expression analysis at different developmental stages and in various tissues showed the life stage-, gut-, salivary gland-, fat body-, Malpighian tubule-, antennae-, ovary- and testis-specific expression patterns of T. molitor CYP genes, implying their various potential roles in development, detoxification, immune response, digestion, olfaction, and reproduction. Our studies provide a platform to understand the evolution of Tenebrionidea CYP gene superfamily, and a basis for further functional investigation of the T. molitor CYPs involved in various biological processes.

Transcriptome and population structure of glassy-winged sharpshooters (Homalodisca vitripennis) with varying insecticide resistance in southern California

Background

Homalodisca vitripennis Germar, the glassy-winged sharpshooter, is an invasive insect in California and a critical threat to agriculture through its transmission of the plant pathogen, Xylella fastidiosa. Quarantine, broad-spectrum insecticides, and biological control have been used for population management of H. vitripennis since its invasion and subsequent proliferation throughout California. Recently wide-spread neonicotinoid resistance has been detected in populations of H. vitripennis in the southern portions of California’s Central Valley. In order to better understand potential mechanisms of H. vitripennis neonicotinoid resistance, we performed RNA sequencing on wild-caught insecticide-resistant and relatively susceptible sharpshooters to profile their transcriptome and population structure.

Results

We identified 81 differentially expressed genes with higher expression in resistant individuals. The significant largest differentially expressed candidate gene linked to resistance status was a cytochrome P450 gene with similarity to CYP6A9. Furthermore, we observed an over-enrichment of GO terms representing functions supportive of roles in resistance mechanisms (cytochrome P450s, M13 peptidases, and cuticle structural proteins). Finally, we saw no evidence of broad-scale population structure, perhaps due to H. vitripennis' relatively recent introduction to California or due to the relatively small geographic scale investigated here.

Conclusions

In this work, we characterized the transcriptome of insecticide-resistant and susceptible H. vitripennis and identified candidate genes that may be involved in resistance mechanisms for this species. Future work should seek to build on the transcriptome profiling performed here to confirm the role of the identified genes, particularly the cytochrome P450, in resistance in H. vitripennis. We hope this work helps aid future population management strategies for this and other species with growing insecticide resistance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08939-1.

Comparative Transcriptome Analysis to Reveal Differentially Expressed Cytochrome P450 in Response to Imidacloprid in the Aphid Lion, Chrysoperla zastrowi sillemi (Esben-Petersen)

The aphid lion, Chrysoperla zastrowi sillemi (Neuroptera: Chrysopidae) is a highly effective beneficial predator of many agricultural pests and has developed resistance to several insecticides. Understanding the molecular mechanism of insecticide resistance in the predators is crucial for its effective application in IPM programs. Therefore, transcriptomes of imidacloprid-resistant and susceptible strains have been assessed using RNA-seq. Cytochrome P450 is one of the important gene families involved in xenobiotic metabolism. Hence, our study focused on the CYP gene family where mining, nomenclature, and phylogenetic analysis revealed a total of 95 unique CYP genes with considerable expansion in CYP3 and CYP4 clans. Further, differential gene expression (DGE) analysis revealed ten CYP genes from CYP3 and CYP4 clans to be differentially expressed, out of which nine genes (CYP4419A1, CYP4XK1, CYP4416A10, CYP4416A-fragment8, CYP6YL1, CYP6YH6, CYP9GK-fragment16, CYP9GN2, CYP9GK6) were downregulated and one (CYP9GK3) was upregulated in the resistant strain as compared to the susceptible strain. Expression validation by quantitative real-time PCR (qRT-PCR) is consistent with the DGE results. The expansion and differential expression of CYP genes may be an indicator of the capacity of the predator to detoxify a particular group of insecticides.

Repeated exposure of fluazinam fungicides affects gene expression profiles yet carries no costs on a nontarget pest

Fungicides are used to control pathogenic fungi of crop species, but they have also been shown to alter behavioral, life history and fitness related traits of nontarget insects. Here, we tested the fungicide effects on feeding behavior, survival and physiology of the nontarget pest insect, the Colorado potato beetle (CPB) (Leptinotarsa decemlineata). Feeding behavior was studied by a choice test of adult beetles, which were allowed to choose between a control and a fungicide (fluazinam) treated potato leaf. Larval survival was recorded after 24 and 72 h exposure to control and fungicide-treated leaves with 2 different concentrations. The adults did not show fungicide avoidance behavior. Similarly, survival of the larvae was not affected by the exposure to fungicides. Finally, to understand the effects of fungicides at the physiological level (gene expression), we tested whether the larval exposure to fungicide alter the expression of 5 metabolic pathway and stress associated genes. Highest concentration and 72-h exposure caused upregulation of 1 cytochrome P450 (CYP9Z14v2) and 1 insecticide resistance gene (Ldace1), whereas metabolic detoxification gene (Ugt1) was downregulated. At 24-h exposure, highest concentration caused downregulation of another common detoxification gene (Gs), while both exposure times to lowest concentration caused upregulation of the Hsp70 stress tolerance gene. Despite these overall effects, there was a considerable amount of variation among different families in the gene expression levels. Even though the behavioral effects of the fungicide treatments were minor, the expression level differences of the studied genes indicate changes on the metabolic detoxifications and stress-related pathways.

Haplotype-resolved genome assembly enables gene discovery in the red palm weevil Rhynchophorus ferrugineus

The red palm weevil Rhynchophorus ferrugineus (Coleoptera: Curculionidae) is an economically-important invasive species that attacks multiple species of palm trees around the world. A better understanding of gene content and function in R. ferrugineus has the potential to inform pest control strategies and thereby mitigate economic and biodiversity losses caused by this species. Using 10x Genomics linked-read sequencing, we produced a haplotype-resolved diploid genome assembly for R. ferrugineus from a single heterozygous individual with modest sequencing coverage (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim$$\end{document}∼ 62x). Benchmarking against conserved single-copy Arthropod orthologs suggests both pseudo-haplotypes in our R. ferrugineus genome assembly are highly complete with respect to gene content, and do not suffer from haplotype-induced duplication artifacts present in a recently published hybrid assembly for this species. Annotation of the larger pseudo-haplotype in our assembly provides evidence for 23,413 protein-coding loci in R. ferrugineus, including over 13,000 predicted proteins annotated with Gene Ontology terms and over 6000 loci independently supported by high-quality Iso-Seq transcriptomic data. Our assembly also includes 95% of R. ferrugineus chemosensory, detoxification and neuropeptide-related transcripts identified previously using RNA-seq transcriptomic data, and provides a platform for the molecular analysis of these and other functionally-relevant genes that can help guide management of this widespread insect pest.

Antennal transcriptome sequencing and identification of candidate chemoreceptor proteins from an invasive pest, the American palm weevil, Rhynchophorus palmarum

For decades, the American palm weevil (APW), Rhynchophorus palmarum, has been a threat to coconut and oil palm production in the Americas. It has recently spread towards North America, endangering ornamental palms, and the expanding date palm production. Its behavior presents several parallelisms with a closely related species, R. ferrugineus, the red palm weevil (RPW), which is the biggest threat to palms in Asia and Europe. For both species, semiochemicals have been used for management. However, their control is far from complete. We generated an adult antennal transcriptome from APW and annotated chemosensory related gene families to obtain a better understanding of these species' olfaction mechanism. We identified unigenes encoding 37 odorant-binding proteins (OBPs), ten chemosensory proteins (CSPs), four sensory neuron membrane proteins (SNMPs), seven gustatory receptors (GRs), 63 odorant receptors (ORs), and 28 ionotropic receptors (IRs). Noticeably, we find out the R. ferrugineus pheromone-binding protein and pheromone receptor orthologs from R. palmarum. Candidate genes identified and annotated in this study allow us to compare these palm weevils' chemosensory gene sets. Most importantly, this study provides the foundation for functional studies that could materialize as novel pest management strategies.

Pheromone receptor of the globally invasive quarantine pest of the palm tree, the red palm weevil (Rhynchophorus ferrugineus)

Palm trees are of immense economic, sociocultural, touristic, and patrimonial significance all over the world, and date palm-related knowledge, traditions, and practices are now included in UNESCOs list of the Intangible Cultural Heritage of Humanity. Of all the pests that infest these trees, the red palm weevil (RPW), Rhynchophorus ferrugineus (Olivier), is its primary enemy. The RPW is a category-1 quarantine insect pest that causes enormous economic losses in palm tree cultivation worldwide. The RPW synchronizes mass gathering on the palm tree for feeding and mating, regulated by a male-produced pheromone composed of two methyl-branched compounds, (4RS, 5RS)-4-methylnonan-5-ol (ferrugineol) and 4(RS)-methylnonan-5-one (ferrugineone). Despite the importance of odorant detection in long-range orientation towards palm trees, palm colonization, and mating, the pheromone receptor has not been identified in this species. In this study, we report the identification and characterization of the first RPW pheromone receptor, RferOR1. Using gene silencing and functional expression in Drosophila olfactory receptor neurons, we demonstrate that RferOR1 is tuned to ferrugineol and ferrugineone and binds five other structurally related molecules. We reveal the lifetime expression of RferOR1, which correlates with adult mating success irrespective of age, a factor that could explain the wide distribution and spread of this pest. As palm weevils are challenging to control based on conventional methods, elucidation of the mechanisms of pheromone detection opens new routes for mating disruption and the early detection of this pest via the development of pheromone receptor-based biosensors.

Function Analysis of P450 and GST Genes to Imidacloprid in Aphis craccivora (Koch)

Aphis craccivora (Koch) is an economically important pest that affects legumes in worldwide. Chemical control is still the primary efficient method for A. craccivora management. However, the mechanism underlying insecticide resistance in A. craccivora has not been elucidated. A previous study observed that piperonyl butoxide (PBO) and diethyl maleate (DEM) significantly synergized imidacloprid in A. craccivora field populations, indicating that cytochrome P450 (P450) and glutathione S-transferase (GST) genes may play pivotal roles in imidacloprid resistance. In this study, 38 P450 genes and 10 GST genes were identified in A. craccivora through transcriptomic analysis. The expression levels of these P450 and GST genes were measured in susceptible (SUS) strains of A. craccivora under imidacloprid treatment with LC₁₅, LC₅₀, and LC₈₅ doses. The expression levels of CYP18A1, CYP6CY21, CYP6DA1, CYP6DA2, CYP4CJ1, CYP4CJ2, and CYP380C6 were up-regulated in the three treatments. Most of these genes belong to CYP3 and CYP4 Clans. In addition, the expression levels of all P450 and GST genes in A. craccivora were also measured in the Juye (JY) and Linqing (LQ) field populations. The expression levels of CYP6DA2, CYP4CJ1, and CYP380C6 were up-regulated in the SUS strain after imidacloprid treatment at three doses, and these genes were overexpressed in the JY population. Furthermore, the sensitivity of A. craccivora to imidacloprid was significantly increased after knockdown of CYP380C6 and CYP6DA2 through RNA interference. These results may help to elucidate the mechanisms underlying of imidacloprid resistance in A. craccivora.

Susceptibility of Four Species of Aphids in Wheat to Seven Insecticides and Its Relationship to Detoxifying Enzymes

Sitobion avenae (Fabricius), Rhopalosiphum padi (Linnaeus), Schizaphis graminum (Rondani), and Metopolophium dirhodum (Walker) (Hemiptera: Aphididae) are important pests of wheat and other cereals worldwide. In this study, the susceptibilities of four wheat aphid species to seven insecticides were assessed. Furthermore, the activities of carboxylesterase (CarE), glutathione S-transferase (GSTs), and cytochrome P450 monooxygenase (P450s) were determined in imidacloprid treated and untreated aphids. The results showed that the susceptibilities of four wheat aphid species to tested insecticides are different and M. dirhodum has shown higher tolerance to most insecticides. Relatively higher CarE and GST activities were observed in M. dirhodum, and P450s activities increased significantly in response to imidacloprid treatment. Moreover, susceptibility to imidacloprid were increased by the oxidase inhibitor piperonyl butoxide in M. dirhodum (20-fold). The results we have obtained imply that P450s may play an important role in imidacloprid metabolic process in M. dirhodum. We suggest that a highly species-specific approach is essential for managing M. dirhodum.

High-throughput profiling of diapause regulated genes from Trichogramma dendrolimi (Hymenoptera: Trichogrammatidae)

Background

The parasitoid wasp, Trichogramma dendrolimi, can enter diapause at the prepupal stage. Thus, diapause is an efficient preservation method during the mass production of T. dendrolimi. Previous studies on diapause have mainly focused on ecological characteristics, so the molecular basis of diapause in T. dendrolimi is unknown. We compared transcriptomes of diapause and non-diapause T. dendrolimi to identify key genes and pathways involved in diapause development.

Results

Transcriptome sequencing was performed on diapause prepupae, pupae after diapause, non-diapause prepupae, and pupae. Analysis yielded a total of 87,022 transcripts with an average length of 1604 bp. By removing redundant sequences and those without significant BLAST hits, a non-redundant dataset was generated, containing 7593 sequences with an average length of 3351 bp. Among them, 5702 genes were differentially expressed. The result of Gene Ontology (GO) enrichment analysis revealed that regulation of transcription, DNA-templated, oxidation-reduction process, and signal transduction were significantly affected. Ten genes were selected for validation using quantitative real-time PCR (qPCR). The changes showed the same trend as between the qPCR and RNA-Seq results. Several genes were identified as involved in diapause, including ribosomal proteins, zinc finger proteins, homeobox proteins, forkhead box proteins, UDP-glucuronosyltransferase, Glutathione-S-transferase, p53, and DNA damage-regulated gene 1 (pdrg1). Genes related to lipid metabolism were also included.

Conclusions

We generated a large amount of transcriptome data from T. dendrolimi, providing a resource for future gene function research. The diapause-related genes identified help reveal the molecular mechanisms of diapause, in T. dendrolimi, and other insect species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07285-4.

Diverse patterns of constitutive and inducible overexpression of detoxifying enzyme genes among resistant Aphis glycines populations

Understanding the mechanisms of pyrethroid resistance is essential to the effective management of pesticide resistance in Aphis glycines Matsumura (Hemiptera: Aphididae). We mined putative detoxifying enzyme genes in the draft genome sequence of A. glycines for cytochrome oxidase P450 (CYP), glutathione-S-transferase (GST) and esterases (E4 and carboxylesterases-CES). Aphids from clonal populations resistant to pyrethroids from three sites in Minnesota, USA, were screened against a diagnostic LC99 concentration of either λ-cyhalothrin or bifenthrin and detoxifying enzyme genes expression in survivors was analyzed by qPCR. Their expression profiles were compared relative to a susceptible clonal population. We found 61 CYP (40 full-length), seven GST (all full-length), seven E4 (five full-length) and three CES (two full-length) genes, including 24 possible pseudogenes. The detoxifying enzymes had different expression profiles across resistant aphid populations, possibly reflecting differences in the genetic background and pyrethroid selection pressures as the number of constitutively overexpressed detoxifying enzyme genes was correlated with the level of resistance. Our findings will strengthen the understanding of the pyrethroid resistance mechanisms in A. glycines.

Susceptibility testing of Anopheles malaria vectors with the neonicotinoid insecticide clothianidin; results from 16 African countries, in preparation for indoor residual spraying with new insecticide formulations

BACKGROUND

In 2017, more than 5 million house structures were sprayed through the U.S. President's Malaria Initiative, protecting more than 21 million people in sub-Saharan Africa. New IRS formulations, SumiShield™ 50WG and Fludora Fusion™ WP-SB, became World Health Organization (WHO) prequalified vector control products in 2017 and 2018, respectively. Both formulations contain the neonicotinoid active ingredient, clothianidin. The target site of neonicotinoids represents a novel mode of action for vector control, meaning that cross-resistance through existing mechanisms is less likely. In preparation for rollout of clothianidin formulations as part of national IRS rotation strategies, baseline susceptibility testing was conducted in 16 countries in sub-Saharan Africa.

METHODS

While work coordinated by the WHO is ongoing to develop a suitable bottle bioassay procedure, there was no published guidance regarding clothianidin susceptibility procedures or diagnostic concentrations. Therefore, a protocol was developed for impregnating filter papers with 2% w/v SumiShield™ 50WG dissolved in distilled water. Susceptibility tests were conducted using insectary-reared reference Anopheles and wild collected malaria vector species. All tests were conducted within 24 h of treating papers, with mortality recorded daily for 7 days, due to the slow-acting nature of clothianidin against mosquitoes. Anopheles gambiae sensu lato (s.l.) adults from wild collected larvae were tested in 14 countries, with wild collected F0 Anopheles funestus s.l. tested in Mozambique and Zambia.

RESULTS

One-hundred percent mortality was reached with all susceptible insectary strains and with wild An. gambiae s.l. from all sites in 11 countries. However, tests in at least one location from 5 countries produced mortality below 98%. While this could potentially be a sign of clothianidin resistance, it is more likely that the diagnostic dose or protocol requires further optimization. Repeat testing in 3 sites in Ghana and Zambia, where possible resistance was detected, subsequently produced 100% mortality. Results showed susceptibility to clothianidin in 38 of the 43 sites in sub-Saharan Africa, including malaria vectors with multiple resistance mechanisms to pyrethroids, carbamates and organophosphates.

CONCLUSIONS

This study provides an interim diagnostic dose of 2% w/v clothianidin on filter papers which can be utilized by National Malaria Control Programmes and research organizations until the WHO concludes multi-centre studies and provides further guidance.