Transcriptome analysis and identification of P450 genes relevant to imidacloprid detoxification in Bradysia odoriphaga

Aliens in the CYPome of the black fungus gnat, Bradysia coprophila

The diverse cytochrome P450 enzymes of insects play essential physiological roles and also play important roles in the metabolism of environmental chemicals such as insecticides. We manually curated the complement of P450 (CYP) genes, or CYPome, of the black fungus gnat, Bradysia (Sciara) coprophila (Diptera, Sciaroidea), a species with a variable number of chromosomes. This CYPome carries two types of "alien" P450 genes. The first type of alien P450s was found among the 163 CYP genes of the core genome (autosomes and X). They consist of 28 sequences resulting from horizontal gene transfer, with closest sequences not found in insects, but in other arthropods, often Collembola. These genes are not contaminants, because they are expressed genes with introns, found in synteny with regular dipteran genes, also found in B. odoriphaga and B. hygida. Two such "alien" genes are representatives of CYP clans not otherwise found in insects, a CYP53 sequence related to fungal CYP53 genes, and a CYP19-like sequence similar to some collembolan sequences but of unclear origin. The second type of alien P450s are represented by 99 sequences from germline-restricted chromosomes (GRC). While most are P450 pseudogenes, 33 are apparently intact, with half being more closely related to P450s from Cecidomyiidae than from Sciaridae, thus supporting the hypothesis of a cross-family hybridization origin of the GRC.

Nanoparticle-Shielded dsRNA Delivery for Enhancing RNAi Efficiency in Cotton Spotted Bollworm Earias vittella (Lepidoptera: Nolidae)

The spotted bollworm Earias vittella (Lepidoptera: Nolidae) is a polyphagous pest with enormous economic significance, primarily affecting cotton and okra. However, the lack of gene sequence information on this pest has a significant constraint on molecular investigations and the formulation of superior pest management strategies. An RNA-seq-based transcriptome study was conducted to alleviate such limitations, and de novo assembly was performed to obtain transcript sequences of this pest. Reference gene identification across E. vittella developmental stages and RNAi treatments were conducted using its sequence information, which resulted in identifying transcription elongation factor (TEF), V-type proton ATPase (V-ATPase), and Glyceraldehyde -3-phosphate dehydrogenase (GAPDH) as the most suitable reference genes for normalization in RT-qPCR-based gene expression studies. The present study also identified important developmental, RNAi pathway, and RNAi target genes and performed life-stage developmental expression analysis using RT-qPCR to select the optimal targets for RNAi. We found that naked dsRNA degradation in the E. vittella hemolymph is the primary reason for poor RNAi. A total of six genes including Juvenile hormone methyl transferase (JHAMT), Chitin synthase (CHS), Aminopeptidase (AMN), Cadherin (CAD), Alpha-amylase (AMY), and V-type proton ATPase (V-ATPase) were selected and knocked down significantly with three different nanoparticles encapsulated dsRNA conjugates, i.e., Chitosan-dsRNA, carbon quantum dots-dsRNA (CQD-dsRNA), and Lipofectamine-dsRNA conjugate. These results demonstrate that feeding nanoparticle-shielded dsRNA silences target genes and suggests that nanoparticle-based RNAi can efficiently manage this pest.

Biological and physiological responses of two Bradysia pests, Bradysia odoriphaga and Bradysia difformis, to Dinotefuran and Lufenuron

Bradysia odoriphaga and Bradysia difformis are destructive root maggots that cause severe losses to vegetables, flowers and edible fungi. Due to the long-term dependence on single pesticides, Bradysia resistance to insecticides has increased, and field control efficacy has decreased obviously. To screen alternative insecticides, and compare the insecticide susceptibility of these two species, we tested the toxicity of eight insecticides to B. odoriphaga and B. difformis, and measured the sublethal effects of Dinotefuran and Lufenuron on life-history parameters and detoxification enzyme activities. Bioassay results indicated that Dinotefuran and Lufenuron had relatively higher toxicity to B. odoriphaga and B. difformis compared to other neonicotinoid and insect growth regulator insecticides, respectively. Significant adverse impacts caused by sublethal concentrations (LC₂₀) of Dinotefuran and Lufenuron on the life-history parameters of F0 and F1 generations of B. odoriphaga and B. difformis were observed. These included reduced survival, prolonged larval development and reduced adult longevity and fecundity. B. odoriphaga had greater resistance and adaptation to insecticides than B. difformis, and an LC₂₀ concentration of Dinotefuran stimulated the reproduction of B. odoriphaga F1 generation and increased the life table parameters. Detoxifying enzymes (CarE and GSTs) and P450 activities fluctuated after a sublethal concentration (Dinotefuran and Lufenuron) treatment, and at the peak value of enzyme activities, the enhancement of detoxifying enzymes of B. odoriphaga was significantly higher than that of B. difformis. These results indicated that Dinotefuran and Lufenuron should be considered as alternatives to other insecticides for control of root maggots. B. odoriphaga exhibited stronger adaptation to insecticides than B. difformis. These data provide guidance for control of root maggots, and the basic information presented here can help reveal the differences in adaptive mechanisms between B. odoriphaga and B. difformis.

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.

Integration analysis of PacBio SMRT- and Illumina RNA-seq reveals P450 genes involved in thiamethoxam detoxification in Bradysia odoriphaga

The sciarid fly Bradysia odoriphaga is a serious pest of Chinese chive (Liliaceae). Neonicotinoid insecticides including thiamethoxam have been used for B. odoriphaga control. However, thiamethoxam resistance in B. odoriphaga has developed in recent years. To identify potential genes involved in detoxification metabolism of thiamethoxam in B. odoriphaga, a PacBio single-molecule real-time (SMRT) transcriptome sequencing and Illumina RNA-seq analysis on thiamethoxam treated B. odoriphaga were performed to explore differentially expressed genes in B. odoriphaga. After SMRT sequencing, analysis of Illumina RNA-Seq data showed a total of 172 differentially expressed genes (DEGs) after thiamethoxam treatment, among which eight upregulated DEGs were P450 genes that may be related to thiamethoxam metabolism. The qRT-PCR results of the eight up-regulated P450 unigenes after thiamethoxam treatment were consistent with RNA-Seq data. Furthermore, oral delivery mediated RNA interference of the eight upregulated P450 transcripts followed by insecticide bioassay was conducted, and three P450 unigenes were verified to be related to thiamethoxam detoxification in B. odoriphaga. This study provides new information about the P450 genes involved in thiamethoxam detoxification in B. odoriphaga.

The phylogenetic and evolutionary analyses of detoxification gene families in Aphidinae species

Detoxification enzymes play significant roles in the interactions between insects and host plants, wherein detoxification-related genes make great contributions. As herbivorous pests, aphids reproduce rapidly due to parthenogenesis. They are good biological materials for studying the mechanisms that allow insect adaptation to host plants. Insect detoxification gene families are associated with insect adaptation to host plants. The Aphidinae is the largest subfamily in the Aphididae with at least 2483 species in 256 genera in 2 tribes: the Macrosiphini (with 3/4 of the species) and the Aphidini. Most aphid pests on crops and ornamental plants are Aphidinae. Members of the Aphidinae occur in nearly every region of the world. The body shape and colour vary significantly. To research the role that detoxification gene families played in the process of aphid adaptation to host evolution, we analyzed the phylogeny and evolution of these detoxification gene families in Aphidinae. In general, the P450/GST/CCE gene families contract, whereas the ABC/UGT families are conserved in Aphidinae species compared to these families in other herbivorous insects. Genus-specific expansions of P450 CYP4, and GST Delta have occurred in the genus Acyrthosiphon. In addition, the evolutionary rates of five detoxification gene families in the evolution process of Aphidinae are different. The comparison of five detoxification gene families among nine Aphidinae species and the estimated relative evolutionary rates provided herein support an understanding of the interaction between and the co-evolution of Aphidinae and plants.

Functionally characterized arthropod pest and pollinator cytochrome P450s associated with xenobiotic metabolism

The cytochrome P450 family (P450s) of arthropods includes diverse enzymes involved in endogenous essential physiological functions and in the oxidative metabolism of xenobiotics, insecticides and plant allelochemicals. P450s can also establish insecticide selectivity in bees and pollinators. Several arthropod P450s, distributed in different phylogenetic groups, have been associated with xenobiotic metabolism, and some of them have been functionally characterized, using different in vitro and in vivo systems. The purpose of this review is to summarize scientific publications on arthropod P450s from major insect and mite agricultural pests, pollinators and Papilio sp, which have been functionally characterized and shown to metabolize xenobiotics and/or their role (direct or indirect) in pesticide toxicity or resistance has been functionally validated. The phylogenetic relationships among these P450s, the functional systems employed for their characterization and their xenobiotic catalytic properties are presented, in a systematic approach, including critical aspects and limitations. The potential of the primary P450-based metabolic pathway of target and non-target organisms for the development of highly selective insecticides and resistance-breaking formulations may help to improve the efficiency and sustainability of pest control.

Selection and Validation of Reference Genes for RT-qPCR Normalization in Bradysia odoriphaga (Diptera: Sciaridae) Under Insecticides Stress

Bradysia odoriphaga (Diptera: Sciaridae) is the most serious root maggot pest which causes substantial damage to the Chinese chive. Organophosphate (OP) and neonicotinoid insecticides are widely used chemical pesticides and play important roles in controlling B. odoriphaga. However, a strong selection pressure following repeated pesticide applications has led to the development of resistant populations of this insect. To understand the insecticide resistance mechanism in B. odoriphaga, gene expression analysis might be required. Appropriate reference gene selection is a critical prerequisite for gene expression studies, as the expression stability of reference genes can be affected by experimental conditions, resulting in biased or erroneous results. The present study shows the expression profile of nine commonly used reference genes [elongation factor 1α (EF-1α), actin2 (ACT), elongation factor 2α (EF-2α), glucose-6-phosphate dehydrogenase (G6PDH), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ribosomal protein L10 (RPL10), ribosomal protein S3 (RPS3), ubiquitin-conjugating enzyme (UBC), and α-tubulin (TUB)] was systematically analyzed under insecticide stress. Moreover, we also evaluated their expression stability in other experimental conditions, including developmental stages, sexes, and tissues. Five programs (NormFinder, geNorm, BestKeeper, RefFinder, and ΔCt) were used to validate the suitability of candidate reference genes. The results revealed that the most appropriate sets of reference genes were RPL10 and ACT across phoxim; ACT and TUB across chlorpyrifos and chlorfluazuron; EF1α and TUB across imidacloprid; EF1α and EF2α across developmental stages; RPL10 and TUB across larvae; EF1α and ACT across tissues, and ACT and G6PDH across sex. These results will facilitate the standardization of RT-qPCR and contribute to further research on B. odoriphaga gene function under insecticides stress.

Integration of transcriptome and proteome reveals molecular mechanisms underlying stress responses of the cutworm, Spodoptera litura, exposed to different levels of lead (Pb)

Heavy metals are major environmental pollutants that affect organisms across different trophic levels. Herbivorous insects play an important role in the bioaccumulation, and eventually, biomagnification of these metals. Although effects of heavy metal stress on insects have been well-studied, the molecular mechanisms underlying their effects remain poorly understood. Here, we used the RNA-Seq profiling and isobaric tags for relative and absolute quantitation (iTRAQ) approaches to unravel these mechanisms in the polyphagous pest Spodoptera litura exposed to lead (Pb) at two different concentrations (12.5 and 100 mg Pb/kg; PbL and PbH, respectively). Altogether, 1392 and 1630 differentially expressed genes (DEGs) and 58, 114 differentially expressed proteins (DEPs) were identified in larvae exposed to PbL and PbH, respectively. After exposed to PbL, the main up-regulated genes clusters and proteins in S. litura larvae were associated with their metabolic processes, including carbohydrate, protein, and lipid metabolism, but the levels of cytochrome P450 associated with the pathway of xenobiotic biodegradation and metabolism were found to be decreased. In contrast, the main up-regulated genes clusters and proteins in larvae exposed to PbH were enriched in the metabolism of xenobiotic by cytochrome P450, drug metabolism-cytochrome P450, and other drug metabolism enzymes, while the down-regulated genes and proteins were found to be closely related to the lipid (lipase) and protein (serine protease, trypsin) metabolism and growth processes (cuticular protein). These findings indicate that S. litura larvae exposed to PbL could enhance food digestion and absorption to prioritize for growth rather than detoxification, whereas S. litura larvae exposed to PbH reduced food digestion and absorption and channelized the limited energy for detoxification rather than growth. These contrasting results explain the dose-dependent effects of heavy metal stress on insect life-history traits, wherein low levels of heavy metal stress induce stimulation, while high levels of heavy metal stress cause inhibition at the transcriptome and proteome levels.

Identification of functional cytochrome P450 and ferredoxin from Streptomyces sp. EAS-AB2608 by transcriptional analysis and their heterologous expression

Bioconversion using microorganisms and their enzymes is an important tool in many industrial fields. The discovery of useful new microbial enzymes contributes to the development of industries utilizing bioprocesses. Streptomyces sp. EAS-AB2608, isolated from a soil sample collected in Japan, can convert the tetrahydrobenzotriazole CPD-1 (a selective positive allosteric modulator of metabotropic glutamate receptor 5) to its hydroxylated form at the C4-(R) position. The current study was performed to identify the genes encoding the enzymes involved in CPD-1 bioconversion and to verify their function. To identify gene products responsible for the conversion of CPD-1, we used RNA sequencing to analyze EAS-AB2608; from its 8333 coding sequences, we selected two genes, one encoding cytochrome P450 (easab2608_00800) and the other encoding ferredoxin (easab2608_00799), as encoding desirable gene products involved in the bioconversion of CPD-1. The validity of this selection was tested by using a heterologous expression approach. A bioconversion assay using genetically engineered Streptomyces avermitilis SUKA24 ∆saverm3882 ∆saverm7246 co-expressing the two selected genes (strain ES_SUKA_63) confirmed that these gene products had hydroxylation activity with respect to CPD-1, indicating that they are responsible for the conversion of CPD-1. Strain ES_SUKA_63 also showed oxidative activity toward other compounds and therefore might be useful not only for bioconversion of CPD-1 but also as a tool for synthesis of drug metabolites and in optimization studies of various pharmaceutical lead compounds. We expect that this approach will be useful for bridging the gap between the latest enzyme optimization technologies and conventional enzyme screening using microorganisms. KEY POINTS: • Genes easab2608_00800 (cyp) and easab2608_00799 (fdx) were selected by RNA-Seq. • Selection validity was evaluated by an engineered S. avermitilis expression system. • Strain ES_SUKA_63 showed oxidative activity toward CPD-1 and other compounds.

Molecular characterization and expression profiles of nicotinic acetylcholine receptors in Bradysia odoriphaga

Bradysia odoriphaga is a destructive insect pest, damaging more than 30 crop species. Nicotinic acetylcholine receptors (nAChRs) mediating fast excitatory transmission in the central nervous system in insects are the molecular targets of some economically important insecticides including imidacloprid, which has been widely used to control B. odoriphaga in China since 2013. However, the clear characterization about nAChRs in B. odoriphaga is still unknown. Hence, our objective is to identify and characterize the nAChR gene family in B. odoriphaga based on the transcriptome database and sequence, phylogenetic and expression profiles analysis. In this study, we cloned seven nAChR subunit genes from B. odoriphaga, including Boα1, Boα2, Boα3, Boα7, Boα8, Boβ1 and Boβ3. Sequence analysis revealed that the seven nAChR subunits of B. odoriphaga shared the typical structural features with Drosophila melanogaster nAChR α1 subunit, including an extracellular N-terminal domain containing six functional loops (loop A-F), a signature Cys-loop with two disulfide bond-forming cysteines separated by 13 amino acid residues, and four typical transmembrane helices (TM1-TM4) in the C-terminal region. Phylogenetic analysis suggested that seven nAChR subunit genes in B. odoriphaga are evolutionarily conserved among four model insects, including D. melanogaster, Bombyx mori, Apis mellifera and Tribolium castaneum. Meanwhile, nAChR α4, α5, α6 and β2 subunit genes may potentially exist in B. odoriphaga, which need further study. Furthermore, quantitative real-time PCR analysis revealed the specific expression pattern of nAChR subunits in three body parts including head, thorax and abdomen, and developmental expression pattern of nAChR subunits throughout the B. odoriphaga life cycle. These results provided necessary information for further investigating the diverse functions of nAChRs in B. odoriphaga.

Survey of the Bradysia odoriphaga Transcriptome Using PacBio Single-Molecule Long-Read Sequencing

The damage caused by Bradysia odoriphaga is the main factor threatening the production of vegetables in the Liliaceae family. However, few genetic studies of B. odoriphaga have been conducted because of a lack of genomic resources. Many long-read sequencing technologies have been developed in the last decade; therefore, in this study, the transcriptome including all development stages of B. odoriphaga was sequenced for the first time by Pacific single-molecule long-read sequencing. Here, 39,129 isoforms were generated, and 35,645 were found to have annotation results when checked against sequences available in different databases. Overall, 18,473 isoforms were distributed in 25 various Clusters of Orthologous Groups, and 11,880 isoforms were categorized into 60 functional groups that belonged to the three main Gene Ontology classifications. Moreover, 30,610 isoforms were assigned into 44 functional categories belonging to six main Kyoto Encyclopedia of Genes and Genomes functional categories. Coding DNA sequence (CDS) prediction showed that 36,419 out of 39,129 isoforms were predicted to have CDS, and 4319 simple sequence repeats were detected in total. Finally, 266 insecticide resistance and metabolism-related isoforms were identified as candidate genes for further investigation of insecticide resistance and metabolism in B. odoriphaga.

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

Background

Cytochrome P450-dependent monooxygenases (P450s), constituting one of the largest and oldest gene superfamilies found in many organisms from bacteria to humans, play a vital role in the detoxification and inactivation of endogenous toxic compounds. The use of various insecticides has increased over the last two decades, and insects have developed resistance to most of these compounds through the detoxifying function of P450s. In this study, we focused on the red palm weevil (RPW), Rhynchophorus ferrugineus, the most devastating pest of palm trees worldwide, and demonstrated through functional analysis that upregulation of P450 gene expression has evolved as an adaptation to insecticide stress arising from exposure to the neonicotinoid-class systematic insecticide imidacloprid.

Results

Based on the RPW global transcriptome analysis, we identified 101 putative P450 genes, including 77 likely encoding protein coding genes with ubiquitous expression. A phylogenetic analysis revealed extensive functional and species-specific diversification of RPW P450s, indicating that multiple CYPs actively participated in the detoxification process. We identified highly conserved paralogs of insect P450s that likely play a role in the development of resistance to imidacloprid: Drosophila Cyp6g1 (CYP6345J1) and Bemisia tabaci CYP4C64 (CYP4LE1). We performed a toxicity bioassay and evaluated the induction of P450s, followed by the identification of overexpressed P450s, including CYP9Z82, CYP6fra5, CYP6NR1, CYP6345J1 and CYP4BD4, which confer cross-resistance to imidacloprid. In addition, under imidacloprid insecticide stress in a date palm field, we observed increased expression of various P450 genes, with CYP9Z82, CYP4BD4, CYP6NR1 and CYP6345J1 being the most upregulated detoxification genes in RPWs. Expression profiling and cluster analysis revealed P450 genes with multiple patterns of induction and differential expression. Furthermore, we used RNA interference to knock down the overexpressed P450s, after which a toxicity bioassay and quantitative expression analysis revealed likely candidates involved in metabolic resistance against imidacloprid in RPW. Ingestion of double-stranded RNA (dsRNA) successfully knocked down the expression of CYP9Z82, CYP6NR1 and CYP345J1 and demonstrated that silencing of CYP345J1 and CYP6NR1 significantly decreased the survival rate of adult RPWs treated with imidacloprid, indicating that overexpression of these two P450s may play an important role in developing tolerance to imidacloprid in a date palm field.

Conclusion

Our study provides useful background information on imidacloprid-specific induction and overexpression of P450s, which may enable the development of diagnostic tools/markers for monitoring the spread of insecticide resistant RPWs. The observed trend of increasing tolerance to imidacloprid in the date palm field therefore indicated that strategies for resistance management are urgently needed.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5837-4) contains supplementary material, which is available to authorized users.

Identification of a novel cytochrome P450 CYP3356A1 linked with insecticide detoxification in Bradysia odoriphaga

BACKGROUND

Cytochrome P450 monooxygenases play an important role in the metabolic detoxification of insecticides in insect pests. However, little is known about the role of a specific P450 gene and its responses to insecticide exposure in Bradysia odoriphaga, a major pest in Chinese chive production.

RESULTS

In this study, a novel P450 gene, CYP3356A1, was cloned from Bradysia odoriphaga. The full-length cDNA sequence of CYP3356A1 is 2153 bp and its open reading frame (ORF) encodes 508 amino acids. Quantitative real time PCR(qRT-PCR) analyses in different tissues showed that CYP3356A1 expression was the highest in the Malpighian tubule. Moreover, among the different developmental stages of the insect, the highest expression of CYP3356A1 was found in fourth-instar larvae. Expression of CYP3356A1 was upregulated by treatment with imidacloprid, thiamethoxam, and β-cypermethrin at median lethal concentrations (LC₅₀ ). RNA interference (RNAi)-mediated silencing of CYP3356A1 significantly increased mortality by 36.90%, 25.17%, and 36.73 when fourth-instar B. odoriphaga larvae were exposed to imidacloprid, thiamethoxam, and β-cypermethrin, respectively, at the LC₅₀ dose.

CONCLUSION

These results demonstrate that CYP3356A1 is related to the detoxification of imidacloprid, thiamethoxam, and β-cypermethrin in B. odoriphaga. Moreover, the study also increased our understanding of the molecular mechanisms of insecticide detoxification in this pest insect. © 2018 Society of Chemical Industry.

Identification and functional analysis of a cytochrome P450 gene involved in imidacloprid resistance in Bradysia odoriphaga Yang et Zhang

Insect cytochrome P450 monooxygenases played an important role in detoxifying insecticides which potentially contributed to the metabolic resistance to insecticides. Bradysia odoriphaga, as a major pest of Chinese chive, was reported to be highly tolerant to neonicotinoid insecticides imidacloprid. In this study, a novel P450 gene, CYP6FV12, was cloned from B. odoriphaga. The full-length cDNA sequence of CYP6FV12 is 2520 bp long and its open reading frame (ORF) encodes 519 amino acids. Quantitative real-time PCR showed that the highest expression of CYP6FV12 was observed in fourth-instar larvae, which is 154.32-fold higher than that of eggs. Highest expression of CYP6FV12 was observed in the midgut, followed by fat body, which was 13.67 and 5.42-fold higher than that in cuticle, respectively. The expression of CYP6FV12 was significantly up-regulated in B. odoriphaga larvae after exposed to imidacloprid at the concentrations of 10, 30, 50, and 70 mg/L. Moreover, RNAi mediated silencing of CYP6FV12 increased mortality by 28.62% when the fourth-instar larvae were treated with imidacloprid. This is the first systematic study on isolated P450s gene involved in imidacloprid resistance in B. odoriphaga and increased our understanding of the molecular mechanisms of insecticide detoxification in this pest insect.