Neofunctionalization of Duplicated P450 Genes Drives the Evolution of Insecticide Resistance in the Brown Planthopper

Multiple Mutations in the β1 Subunit of the Nicotinic Acetylcholine Receptor Confer Resistance to Neonicotinoids

The evolution of resistance to neonicotinoid insecticides threatens global agriculture. To elucidate its molecular basis, we employed Drosophila melanogaster as a model system to investigate resistance-associated mutations in the β1 subunit of nicotinic acetylcholine receptors (nAChRs). Using a CRISPR/Cas9-mediated allele replacement, we generated homozygous knock-in mutants (V62I, V101I, R81E, and A60T,R81E) without apparent fitness costs. Toxicity bioassays revealed that these mutations confer varying resistance levels, with the R81E mutation exhibiting over 225-fold resistance to thiamethoxam, clothianidin, and dinotefuran. A heteropentameric α1β1 nAChR model, generated using an AI-based protein-ligand prediction (Chai-1), showed that imidacloprid binds at the orthosteric site, where R81 forms a critical electrostatic interaction. Residues A60, V62, and V101, positioned further from the binding site, showed spatial distances correlated with their resistance ratios. These findings provide genetic and structural insights into neonicotinoid resistance mechanisms, offering a foundation for the design of next-generation insecticides and resistance management strategies.

microRNA Targeting Cytochrome P450 Is Involved in Chlorfenapyr Tolerance in the Silkworm, Bombyx mori (Lepidoptera: Bombycidae)

We first measured the content of chlorfenapyr and tralopyril in silkworm larvae using HPLC, revealing that chlorfenapyr can be biotransformed into tralopyril in silkworms. Then, a differential transcriptomic database of small RNA was constructed through Illumina RNA-Sequencing. qRT-PCR was conducted to determine the expression levels of Bmo-miR-6497-5p and the target CYP450 gene, and Bmo-miR-6497-5p was significantly upregulated in the L3 silkworm larvae 24, 48, and 72 h after they were treated with chlorfenapyr. Furthermore, the target P450 gene CYP337A2 was downregulated at these time points. Dual-luciferase validation revealed that the luciferase activity significantly decreased after Bmo-miR-6497-5p bound to CYP337A2. In addition, miRNA mimics/inhibitor injection and bioassays of chlorfenapyr and tralopyril revealed that the mortality of third silkworm larvae injected with the antagomir of Bmo-miR-6497-5p was significantly increased after exposure to a sublethal concentration of chlorfenapyr. These results imply that Bmo-miR-6497-5p targets CYP337A2, regulating its expression. Also, silkworms increase their tolerance to chlorfenapyr by upregulating Bmo-miR-6497-5p expression, thereby inhibiting the biotransformation of chlorfenapyr to toxic tralopyril catalyzed by CYP337A2. The present study reveals the function of microRNA in silkworm tolerance to chlorfenapyr and improves understanding regarding insecticide resistance in Lepidopteran insects.

Cythochrome P450-mediated dinotefuran resistance in onion thrips, Thrips tabaci

Onion thrips, Thrips tabaci, have developed resistance to many insecticides, and over the last decade, resistant populations have spread widely across Japan. The cytochrome P450 (CYP) family, a widely conserved detoxification enzyme that metabolizes xenobiotics such as insecticides and phytochemicals, is believed to play important roles in the development of resistance in T. tabaci. However, CYPs involved in insecticide resistance in T. tabaci remain unclear. To comprehensively identify CYPs in T. tabaci, the genome sequences of the thelytokous T. tabaci (ANO strain) were constructed, and 18,965 genes (protein coding) were predicted. We identified 127 CYP genes in the predicted gene set by manual curation, and 38 of these CYP genes belonged to the CYP3 clan, including genes from the CYP6 family, which is one of the most important CYP families involved in resistance to neonicotinoids in many insect pests. To identify the CYPs involved in resistance to dinotefuran, which is one of the neonicotinoids used to control T. tabaci, RNA sequencing of dinotefuran-resistant and dinotefuran-susceptible strains was performed. Results revealed that, TtCYP3652A1, which belongs to the thrips-specific CYP3652A subfamily in the CYP3 clan, was significantly upregulated in the resistant strain. In vitro CYP metabolism assays using insect cells were conducted for TtCYP3652A1 and five highly expressed CYP6 genes. Only TtCYP3652A1 significantly metabolized dinotefuran, which is considered to contribute to detoxification of dinotefuran. As no amino acid mutations were identified in the known target-site genes of neonicotinoids, TtCYP3652A1 was considered to be the main factor involved in the resistance to dinotefuran.

Functional analysis of AccCPR in Apis cerana cerana under pesticide and heavy metal stress

NADPH-cytochrome P450 reductase (CPR) plays important roles in the metabolism of both endogenous and exogenous compounds through cytochrome P450, and is also involved in the detoxification of insecticides mediated by cytochrome P450. However, the CPR from Apis cerana cerana has not been well characterized and its function is still undescribed. This study isolated the CPR gene from Apis cerana cerana and investigated its functional role in the resistance to pesticide and heavy metal stress. Bioinformatic analysis revealed significant homology between the gene and its counterparts in other species. Functional investigations demonstrated diverse expression and localization patterns of this gene, with AccCPR primarily expressed in muscular tissues and the gut, suggesting its potential roles in flight activities and intestinal barrier function of bees. Furthermore, the expression levels of this gene were significantly modulated under pesticide and heavy metal stress. Notably, the overexpression of AccCPR led to a marked alteration the tolerance to external stressors in E. coli. Additionally, the silencing of the AccCPR gene resulted in a significant decrease in antioxidant enzyme activity and the expression levels of genes associated with antioxidant functions. Consequently, the mortality rate of Apis cerana cerana under imidacloprid stress was significantly elevated. Taken together, our findings suggest that AccCPR may play a pivotal role in the resistance of Apis cerana cerana to abiotic stresses such as pesticides and heavy metals by regulating antioxidant pathways.

Chemosensory Proteins Protect Nilaparvata lugens from Imidacloprid by Sequestering the Insecticide and Facilitating Metabolic Detoxification

The involvement of chemosensory proteins (CSPs) in binding to insecticides has been implicated. However, our understanding of CSP-mediated insecticide resistance remains limited. Herein, 15 CSP genes were identified and characterized from Nilaparvata lugens. Expression analysis identified six CSPs with overexpression in the imidacloprid-resistant strain, whose involvement in imidacloprid resistance was validated by RNA interference. Among them, four CSPs were successfully expressed using a prokaryotic expression system, and their binding affinities to imidacloprid were confirmed through fluorescence competitive binding assays. Knockdown of them impaired the capacity of N. lugens to metabolize imidacloprid and inhibited the activity of metabolic detoxification pathways, while their overexpression in Escherichia coli enhanced bacterial metabolic efficiency toward imidacloprid. Furthermore, the transcriptional regulation of CSP2 and CSP15 was found to be mediated by AhR/ARNT and CncC/MafK. These findings suggest that the overexpression of CSPs in N. lugens promotes imidacloprid resistance by sequestering the insecticide and enhancing metabolic detoxification.

Differentially spliced mitochondrial CYP419A1 contributes to ethiprole resistance in Nilaparvata lugens

The brown planthopper Nilaparvata lugens is one of the most economically important pests of cultivated rice in Southeast Asia. Extensive use of insecticide treatments, such as imidacloprid, fipronil and ethiprole, has resulted in the emergence of multiple resistant strains of N. lugens. Previous investigation of the mechanisms of resistance to imidacloprid and ethiprole demonstrated that overexpression and qualitative changes in the cytochrome P450 gene CYP6ER1 lead to enhanced metabolic detoxification of these compounds. Here, we present the identification of a secondary mechanism enhancing ethiprole resistance mediated by differential splicing and overexpression of CYP419A1, a planthopper-specific, mitochondrial P450 gene. Although metabolic resistance to insecticides is usually mediated by overexpression of P450 genes belonging to either CYP 3 or 4 clades, we validate the protective effect of over-expression of CYP419A1, in vivo, using transgenic Drosophila melanogaster. Additionally, we report some unusual features of both the CYP419A1 gene locus and protein, which include, altered splicing associated with resistance, a non-canonical heme-binding motif and an extreme 5' end extension of the open reading frame. These results provide insight into the molecular mechanisms underpinning resistance to insecticides and have applied implications for the control of a highly damaging crop pest.

Cis-regulation of the CYP6CS1 gene and its role in mediating cross-resistance in a pymetrozine-resistant strain of Nilaparvata lugens

Pymetrozine is currently one of the primary insecticides used to control the brown planthopper, Nilaparvata lugens Stål (Hemiptera: Delphacidae), but the long-term effectiveness of this chemical is threatened by growing issues of resistance. Previous studies in a laboratory selected strain of N. lugens, Pym-R, have shown that resistance to pymetrozine can evolve without target-site mutations. A key candidate gene identified is the cytochrome P450 gene CYP6CS1, which is overexpressed in the resistant Pym-R strain compared to the laboratory susceptible strain, Pym-S. In this study, we provide a deeper characterization of the regulatory mechanism and phenotypic effects of CYP6CS1 by comparing the resistant and susceptible variants of this gene. Using artificial constructs in Luciferase activity assays, we elucidate the role of indels in the overexpression of CYP6CS1 in the resistant strain. Additionally, transgenic Drosophila experiments also revealed that the CYP6CS1 gene not only contributes to resistance against pymetrozine, but is able to confer moderate to low cross-resistance to several other pesticides. This research provides vital insights into the possible genetic mechanisms that may contribute to pymetrozine resistance in field populations. Future work will aim to examine the relevance of CYP6CS1 variation in the field with the aim of developing diagnostic markers of resistance.

Effects of Amino Acid Mutation in Cytochrome P450 (CYP96A146) of Descurainia sophia on the Metabolism and Resistance to Tribenuron-Methyl

Cytochrome P450 monooxygenases (P450s) play important roles in herbicide resistance. In this study, there are four amino acid mutations (F39Y, H163Y, S203A, and V361E) between CYP96A146-S and CYP96A146-R, which were cloned, respectively, from susceptible (S) and tribenuron-methyl-resistant (TR) Descurainia sophia. The Arabidopsis expressing CYP96A146-S or CYP96A146-R showed resistance to tribenuron-methyl, carfentrazone-ethyl, and oxyfluorfen, while Arabidopsis transformed with CYP96A146-R or CYP96A146 with any two or three mutations of H163Y, S203A, or V361E exhibited significantly higher resistance to tribenuron-methyl than Arabidopsis expressing CYP96A146-S. The metabolic rates of tribenuron-methyl were significantly faster in Arabidopsis expressing CYP96A146-R than that with CYP96A146-S. The molecular dynamics simulation demonstrated that amino acid mutations did not affect the domain of the HEM ring, which could significantly enhance the volume of the catalytic pocket in P450 (CYP96A146), thereby increasing the collision rate between the catalytic pocket and tribenuron-methyl. Hence, the amino acid mutations may be one of the mechanisms underlying P450-mediated herbicide resistance.

A single mutation G454A in the P450 CYP9K1 drives pyrethroid resistance in the major malaria vector Anopheles funestus reducing bed net efficacy

Metabolic mechanisms conferring pyrethroid resistance in malaria vectors are jeopardizing the effectiveness of insecticide-based interventions, and identification of their markers is a key requirement for robust resistance management. Here, using a field-lab-field approach, we demonstrated that a single mutation G454A in the P450 CYP9K1 is driving pyrethroid resistance in the major malaria vector Anopheles funestus in East and Central Africa. Drastic reduction in CYP9K1 diversity was observed in Ugandan samples collected in 2014, with the selection of a predominant haplotype (G454A mutation at 90%), which was completely absent in the other African regions. However, 6 years later (2020) the Ugandan 454A-CYP9K1 haplotype was found predominant in Cameroon (84.6%), but absent in Malawi (Southern Africa) and Ghana (West Africa). Comparative in vitro heterologous expression and metabolism assays revealed that the mutant 454A-CYP9K1 (R) allele significantly metabolizes more type II pyrethroid (deltamethrin) compared with the wild G454-CYP9K1 (S) allele. Transgenic Drosophila melanogaster flies expressing 454A-CYP9K1 (R) allele exhibited significantly higher type I and II pyrethroids resistance compared to flies expressing the wild G454-CYP9K1 (S) allele. Furthermore, laboratory testing and field experimental hut trials in Cameroon demonstrated that mosquitoes harboring the resistant 454A-CYP9K1 allele significantly survived pyrethroids exposure (odds ratio = 567, P < 0.0001). This study highlights the rapid spread of pyrethroid-resistant CYP9K1 allele, under directional selection in East and Central Africa, contributing to reduced bed net efficacy. The newly designed DNA-based assay here will add to the toolbox of resistance monitoring and improving its management strategies.

F116I mutation in CYP9A25 associated with resistance to emamectin benzoate in Spodoptera litura

BACKGROUND

The F116V mutation in the substrate recognition site 1 (SRS1) of Spodoptera exigua CYP9A186 has been demonstrated to confer ~200-fold resistance to emamectin benzoate (EB). In this study, a novel mutation (F116I) in CYP9A25, orthologous to CYP9A186, was detected in a field population of Spodoptera litura (YJ22) collected from Yuanjiang, Yunnan province, China in 2022. The association of this mutation with EB resistance was investigated.

RESULTS

Two homozygous strains, YJ22-116F (wild-type at 116 position of CYP9A25) and YJ22-116I (mutant) were isolated from YJ22 through two rounds of crossing and DNA genotyping. Compared with YJ22-116F, the mutant strain YJ22-116I exhibited 31.8-fold resistance to EB. Resistance in YJ22-116I was shown to be incompletely dominant, and genetically linked with the F116I mutation. Further, heterologous expression and in vitro metabolism assays confirmed that the recombinant CYP9A25 protein with 116I mutation obtained metabolic capability against EB, whereas the wild-type CYP9A25 protein (with 116F) did not metabolize EB. Molecular modeling showed that the F116I mutation within SRS1 reduces the steric hindrance to substrate entry and improves ligand-binding interactions.

CONCLUSION

The causal association between the F116I mutation in CYP9A25 and medium-level EB resistance in S. litura has been verified. This finding is critical for the field monitoring of such mutations and thus for developing adaptive resistance management tactics for field populations of S. litura. © 2024 Society of Chemical Industry.

Two highly selected mutations in the tandemly duplicated CYP6P4a and CYP6P4b genes drive pyrethroid resistance in Anopheles funestus in West Africa

BACKGROUND

Gaining a comprehensive understanding of the genetic mechanisms underlying insecticide resistance in malaria vectors is crucial for optimising the effectiveness of insecticide-based vector control methods and developing diagnostic tools for resistance management. Considering the heterogeneity of metabolic resistance in major malaria vectors, the implementation of tailored resistance management strategies is essential for successful vector control. Here, we provide evidence demonstrating that two highly selected mutations in CYP6P4a and CYP6P4b are driving pyrethroid insecticide resistance in the major malaria vector Anopheles funestus, in West Africa.

RESULTS

Continent-wide polymorphism survey revealed escalated signatures of directional selection of both genes between 2014 and 2021. In vitro insecticide metabolism assays with recombinant enzymes from both genes showed that mutant alleles under selection exhibit higher metabolic efficiency than their wild-type counterparts. Using the GAL4-UAS expression system, transgenic Drosophila flies overexpressing mutant alleles exhibited increased resistance to pyrethroids. These findings were consistent with in silico predictions which highlighted changes in enzyme active site architecture that enhance the affinity of mutant alleles for type I and II pyrethroids. Furthermore, we designed two DNA-based assays for the detection of CYP6P4a-M220I and CYP6P4b-D284E mutations, showing their current confinement to West Africa. Genotype/phenotype correlation analyses revealed that these markers are strongly associated with resistance to types I and II pyrethroids and combine to drastically reduce killing effects of pyrethroid bed nets.

CONCLUSIONS

Overall, this study demonstrated that CYP6P4a and CYP6P4b contribute to pyrethroid resistance in An. funestus and provided two additional insecticide resistance molecular diagnostic tools that would contribute to monitoring and better management of resistance.

The Evolution and Mechanisms of Multiple-Insecticide Resistance in Rice Stem Borer, Chilo suppressalis Walker (Lepidoptera: Crambidae)

The emergence of insecticide resistance in the rice stem borer, Chilo suppressalis, is a growing threat to the sustainable control of this important insect crop pest. Thus, monitoring of C. suppressalis populations for insecticide resistance and characterization of the underlying genetic mechanisms is essential to inform rational control decisions and the development of resistance management strategies. Here, we monitored 126 C. suppressalis field populations from China for resistance evolution to four major insecticides: 53 for chlorantraniliprole, 50 for abamectin, 74 for triazophos, and 76 for spinetoram. Moderate to high levels of resistance were observed to all four insecticides. Investigation of the underlying resistance mechanisms revealed multiple mutations in the ryanodine receptor (RyR) and acetylcholinesterase 1 (AChE1), leading to target-site resistance to chlorantraniliprole and triazophos, respectively. In contrast, the absence of mutations in the glutamate-gated chloride channel (GluCl) and α6 nicotinic acetylcholine receptor (nAChR α6) subunit suggested that nontarget site mechanisms contribute to the multiple-insecticide resistance phenotypes observed in C. suppressalis. In this regard, we revealed overexpression of the uridine 5'-diphospho-glycosyltransferase UGT33AF1 and cytochrome P450 CYP6AB45 in C. suppressalis field populations. Functional characterization using transgenic Drosophila demonstrated that UGT33AF1 confers resistance against multiple insecticides in vivo, whereas CYP6AB45 does not appear to contribute to resistance. Collectively, our findings reveal the current status of resistance of C. suppressalis to insecticides in China and uncover a diverse profile of resistance mechanisms in this species. These findings provide a foundation for the development of sustainable strategies to effectively manage and control this pest.

Flavin-Dependent Monooxgenase Confers Resistance to Chlorantraniliprole and Spinetoram in the Rice Stem Borer Chilo suppressalis Walker (Lepidoptera: Crambidae)

Understanding the role of flavin-containing monooxygenases (FMOs) in the genetic mechanisms of insecticide resistance is essential for developing effective management strategies against the rice stem borer, Chilo suppressalis. In this study, we identified five FMO genes in C. suppressalis, examined their expression patterns, and revealed overexpression of FMO3B and FMO3C in field populations resistant to multiple insecticides, including chlorantraniliprole and spinetoram. Functional characterization using transgenic Drosophila indicated that FMO3B and FMO3C do not confer resistance to abamectin or methoxyfenozide but do mediate resistance to chlorantraniliprole and spinetoram. Knockdown of FMO3B and FMO3C increased sensitivity to these insecticides in C. suppressalis. Molecular docking studies indicated direct binding of chlorantraniliprole and spinetoram to these FMOs, underscoring their role in metabolic resistance. These findings indicate that FMOs are key enzymes in the metabolic resistance of C. suppressalis to chlorantraniliprole and spinetoram, enhancing our understanding of insecticide resistance and aiding the development of management strategies.

STAT5B, Akt and p38 Signaling Activate FTZ-F1 to Regulate the Xenobiotic Tolerance-Related Gene SlCyp9a75b in Spodoptera litura

Cytochrome P450 monooxygenases in insects have been verified to implicated in insecticide and phytochemical detoxification metabolism. However, the regulation of P450s, which are modulated by signal-regulated transcription factors (TFs), is less well studied in insects. Here, we found that the Malpighian tubule specific P450 gene SlCYP9A75b in Spodoptera litura is induced by xenobiotics. The transgenic Drosophila bioassay and RNAi results indicated that this P450 gene contributes to α-cypermethrin, cyantraniliprole, and nicotine tolerance. In addition, functional analysis revealed that the MAPKs p38, PI3K/Akt, and JAK-STAT activate the transcription factor fushi tarazu factor 1 (FTZ-F1) to regulate CYP9A75b expression. These findings provide mechanistic insights into the contributions of CYP9A genes to xenobiotic detoxification and support the possible involvement of different signaling pathways and TFs in tolerance to xenobiotics in insects.

Genome structural variants shape adaptive success of an invasive urban malaria vector Anopheles stephensi

Global changes are associated with the emergence of several invasive species. However, the genomic determinants of the adaptive success of an invasive species in a new environment remain poorly understood. Genomic structural variants (SVs), consisting of copy number variants, play an important role in adaptation. SVs often cause large adaptive shifts in ecologically important traits, which makes SVs compelling candidates for driving rapid adaptations to environmental changes, which is critical to invasive success. To address this problem, we investigated the role SVs play in the adaptive success of Anopheles stephensi , a primary vector of urban malaria in South Asia and an invasive malaria vector in several South Asian islands and Africa. We collected whole genome sequencing data from 115 mosquitoes from invasive island populations and four locations from mainland India, an ancestral range for the species. We identified 2,988 duplication copy number variants and 16,038 deletions in these strains, with ∼50% overlapping genes. SVs are enriched in genomic regions with signatures of selective sweeps in the mainland and invasive island populations, implying a putative adaptive role of SVs. Nearly all high-frequency SVs, including the candidate adaptive variants, in the invasive island populations are present on the mainland, suggesting a major contribution of existing variation to the success of the island populations. Among the candidate adaptive SVs, three duplications involving toxin-resistance genes evolved, likely due to the widespread application of insecticides in India since the 1950s. We also identify two SVs associated with the adaptation of An. stephensi larvae to brackish water in the island and two coastal mainland populations, where the mutations likely originated. Our results suggest that existing SVs play a vital role in the evolutionary success of An. stephensi in new environmental conditions.

Cytochrome P450 Enzyme Design by Constraining the Catalytic Pocket in a Diffusion Model

Although cytochrome P450 enzymes are the most versatile biocatalysts in nature, there is insufficient comprehension of the molecular mechanism underlying their functional innovation process. Here, by combining ancestral sequence reconstruction, reverse mutation assay, and progressive forward accumulation, we identified 5 founder residues in the catalytic pocket of flavone 6-hydroxylase (F6H) and proposed a "3-point fixation" model to elucidate the functional innovation mechanisms of P450s in nature. According to this design principle of catalytic pocket, we further developed a de novo diffusion model (P450Diffusion) to generate artificial P450s. Ultimately, among the 17 non-natural P450s we generated, 10 designs exhibited significant F6H activity and 6 exhibited a 1.3- to 3.5-fold increase in catalytic capacity compared to the natural CYP706X1. This work not only explores the design principle of catalytic pockets of P450s, but also provides an insight into the artificial design of P450 enzymes with desired functions.

The G932C mutation of chitin synthase 1 gene (CHS1) mediates buprofezin resistance as confirmed by CRISPR/Cas9-mediated knock-in approach in the brown planthopper, Nilaparvata lugens

The brown planthopper (Nilaparvata lugens) is a major destructive rice pest in Asia. High levels of insecticide resistance have been frequently reported, and the G932C mutation in the chitin synthase 1 (CHS1) gene has been found to mediate buprofezin resistance. However, there has been no direct evidence to confirm the functional significance of the single G932C substitution mutation leading to buprofezin resistance in N. lugens. Here, we successfully constructed a knock-in homozygous strain (Nl-G932C) of N. lugens using CRISPR/Cas9 coupled with homology-directed repair (HDR). Compared with the background strain susceptible to buprofezin (Nl-SS), the knock-in strain (Nl-G932C) showed a 94.9-fold resistance to buprofezin. Furthermore, resistant strains (Nl-932C) isolated from the field exhibited a 2078.8-fold resistance to buprofezin, indicating that there are other mechanisms contributing to buprofezin resistance in the field. Inheritance analysis showed that the resistance trait is incomplete dominance. In addition, the Nl-G932C strain had a relative fitness of 0.33 with a substantially decreased survival rate, emergence rate, and fecundity. This study provided in vivo functional evidence for the causality of G932C substitution mutation of CHS1 with buprofezin resistance and valuable information for facilitating the development of resistance management strategies in N. lugens. This is the first example of using CRISPR/Cas9 gene-editing technology in a hemipteran insect to directly confirm the role of a candidate target site mutation in insecticide resistance.

Insecticide resistance reduces the profitability of insect-resistant rice cultivars

INTRODUCTION

Preventing crop yield loss caused by pests is critical for global agricultural production. Agricultural pest control has largely relied on chemical pesticides. The interaction between insecticide resistance and the adaptation of herbivorous pests to host plants may represent an emerging threat to future food security.

OBJECTIVES

This study aims to unveil genetic evidence for the reduction in the profitability of resistant cultivars derived from insecticide resistance in target pest insects.

METHODS

An experimental evolution system encompassing resistant rice and its major monophagous pest, the brown planthopper Nilaparvata lugens, was constructed. Whole genome resequencing and selective sweep analysis were utilized to identify the candidate gene loci related to the adaptation. RNA interference and induced expression assay were conducted to validate the function of the candidate loci.

RESULTS

We found that the imidacloprid-resistant population of N. lugens rapidly adapted to resistant rice IR36. Gene loci related to imidacloprid resistance may contribute to this phenomenon. Multiple alleles in the nicotinic acetylcholine receptor (nAChR)-7-like and P450 CYP4C61 were significantly correlated with changes in virulence to IR36 rice and insecticide resistance of N. lugens. One avirulent/susceptible genotype and two virulent/resistant genotypes could be inferred from the corresponding alleles. Importantly, we found that the virulent/resistant genotypes already exist in the wild in China, exhibiting increasing frequencies along with insecticide usage. We validated the relevance of these genotypes and the virulence to three more resistant rice cultivars. Knockdown of the above two genes in N. lugens significantly decreased both the resistance to imidacloprid and the virulence towards resistant rice.

CONCLUSION

Our findings provide direct genetic evidence to the eco-evolutionary consequence of insecticide resistance, and suggest an urgent need for the implementation of predictably sustainable pest management.

Genome assembly and transcriptomic analysis to elucidate the ability of Nasonovia ribisnigri to break host plant resistance

Aphid genomic resources enable the study of complex life history traits and provide information on vector biology, host adaption and speciation. The currant-lettuce aphid (Nasonovia ribisnigri (Hemiptera: Aphididae) (Mosley)) is a cosmopolitan pest of outdoor lettuce (Lactuca sativa (Asterales: Asteraceae) (Linnaeus)). Until recently, the use of resistant cultivars was an effective method for managing N. ribisnigri. A resistant cultivar containing a single gene (Nr-locus), introduced in the 1980s, conferred complete resistance to feeding. Overreliance of this Nr-locus in lettuce resulted in N. ribisnigri's ability to break resistance mechanism, with first reports during 2003. Our work attempts to understand which candidate gene(s) are associated with this resistance-breaking mechanism. We present two de novo draft assembles for N. ribisnigri genomes, corresponding to both avirulent (Nr-locus susceptible) and virulent (Nr-locus resistant) biotypes. Changes in gene expression of the two N. ribisnigri biotypes were investigated using transcriptomic analyses of RNA-sequencing (RNA-seq) data to understand the potential mechanisms of resistance to the Nr-locus in lettuce. The draft genome assemblies were 94.2% and 91.4% complete for the avirulent and virulent biotypes, respectively. Out of the 18,872 differentially expressed genes, a single gene/locus was identified in N. ribisnigri that was shared between two resistant-breaking biotypes. This locus was further explored and validated in Real-Time Quantitative Reverse Transcription PCR (qRT-PCR) experiments and has predicted localisations in both the cytoplasm and nucleus. This is the first study to provide evidence that a single gene/locus is likely responsible for the ability of N. ribisnigri to overcome the Nr-locus resistance in the lettuce host.

The molecular determinants of pesticide sensitivity in bee pollinators

Bees carry out vital ecosystem services by pollinating both wild and economically important crop plants. However, while performing this function, bee pollinators may encounter potentially harmful xenobiotics in the environment such as pesticides (fungicides, herbicides and insecticides). Understanding the key factors that influence the toxicological outcomes of bee exposure to these chemicals, in isolation or combination, is essential to safeguard their health and the ecosystem services they provide. In this regard, recent work using toxicogenomic and phylogenetic approaches has begun to identify, at the molecular level, key determinants of pesticide sensitivity in bee pollinators. These include detoxification systems that convert pesticides to less toxic forms and key residues in insecticide target-sites that underlie species-specific insecticide selectivity. Here we review this emerging body of research and summarise the state of knowledge of the molecular determinants of pesticide sensitivity in bee pollinators. We identify gaps in our knowledge for future research and examine how an understanding of the genetic basis of bee sensitivity to pesticides can be leveraged to, a) predict and avoid negative bee-pesticide interactions and facilitate the future development of pest-selective bee-safe insecticides, and b) inform traditional effect assessment approaches in bee pesticide risk assessment and address issues of ecotoxicological concern.

New insights into chlorantraniliprole metabolic resistance mechanisms mediated by the striped rice borer cytochrome P450 monooxygenases: A case study of metabolic differences

The anthranilic diamide insecticide chlorantraniliprole has been extensively applied to control Lepidoptera pests. However, its overuse leads to the development of resistance and accumulation of residue in the environment. Four P450s (CYP6CV5, CYP9A68, CYP321F3, and CYP324A12) were first found to be constitutively overexpressed in an SSB CAP-resistant strain. It is imperative to further elucidate the molecular mechanisms underlying P450s-mediated CAP resistance for mitigating its environmental contamination. Here, we heterologously expressed these four P450s in insect cells and evaluated their abilities to metabolize CAP. Western blotting and reduced CO difference spectrum tests showed that these four P450 proteins had been successfully expressed in Sf9 cells, which are indicative of active functional enzymes. The recombinant proteins CYP6CV5, CYP9A68, CYP321F3, and CYP324A12 exhibited a preference for metabolizing the fluorescent P450 model probe substrates EC, BFC, EFC, and EC with enzyme activities of 0.54, 0.67, 0.57, and 0.46 pmol/min/pmol P450, respectively. In vitro metabolism revealed distinct CAP metabolic rates (0.97, 0.86, 0.75, and 0.55 pmol/min/pmol P450) and efficiencies (0.45, 0.37, 0.30, and 0.17) of the four recombinant P450 enzymes, thereby elucidating different protein catalytic activities. Furthermore, molecular model docking confirmed metabolic differences and efficiencies of these P450s and unveiled the hydroxylation reaction in generating N-demethylation and methylphenyl hydroxylation during CAP metabolism. Our findings not only first provide new insights into the mechanisms of P450s-mediated metabolic resistance to CAP at the protein level in SSB but also demonstrate significant differences in the capacities of multiple P450s for insecticide degradation and facilitate the evaluation and mitigation of toxic risks associated with CAP application in the environment.

Two Point Mutations in CYP4CE1 Promoter Contributed to the Differential Regulation of CYP4CE1 Expression by FoxO between Susceptible and Nitenpyram-Resistant Nilaparvata lugens

Four P450s were reported to be important for imidacloprid resistance in Nilaparvata lugens, a major insect pest on rice, which was confirmed in this study in an imidacloprid-resistant strain (ImiR). Here we found that only two (CYP4CE1 and CYP6ER1) from these four P450 genes were overexpressed in a nitenpyram-resistant strain (NitR) when compared to a susceptible strain (SUS). CYP4CE1 RNAi reduced nitenpyram and imidacloprid resistance in NitR and ImiR strains, with a greater reduction in nitenpyram resistance. The transcription factor FoxO mediated nitenpyram resistance in NitR and ImiR strains, but it was not differentially expressed among strains. The potential reason for the differential regulation of FoxO on CYP4CE1 expression was mainly from sequence differences in the CYP4CE1 promoter between susceptible and resistant insects. In six FoxO response elements predicted in the CYP4CE1 promoter, the single-nucleotide polymorphisms were frequently detected in over 50% of NitR and ImiR individuals. The luciferase reporter assays showed that two mutations, -650T/G and -2205T/A in two response elements at the positions of -648 and -2200 bp, mainly contributed to the enhanced regulation on CYP4CE1 expression by FoxO in resistant insects. The frequency was over 69% for both -650T/G and -2205T/A detected in NitR and ImiR individuals but less than 20% in SUS insects. In conclusion, CYP4CE1 overexpression importantly contributed to nitenpyram resistance in N. lugens, and two mutations in the CYP4CE1 promoter of resistant insects led to an enhanced regulation on CYP4CE1 expression by FoxO.

Up-regulation of CYP6G4 mediated by a CncC/maf binding-site-containing insertion confers resistance to multiple classes of insecticides in the house fly Musca domestica

Populations of many insect species have evolved a variety of resistance mechanisms in response to insecticide selection. Current knowledge about mutations responsible for insecticide resistance is largely achieved from studies on target-site resistance, while much less is known about metabolic resistance. Although it is well known that P450 monooxygenases are one of the major players involved in insecticide metabolism and resistance, understanding mutation(s) responsible for CYP-mediated resistance has been a big challenge. In this study, we used the house fly to pursue a better understanding of P450 mediated insecticide resistance at the molecular level. Metabolism studies illustrated that CYP6G4 had a broad-spectrum metabolic activity in metabolizing insecticides. Population genotyping revealed that the CYP6G4v1 allele harboring a DNA insertion (MdIS1) had been selected in many house fly populations on different continents. Dual luciferase reporter assays identified that the MdIS1 contained a CncC/Maf binding site, and electrophoretic mobility shift assay confirmed that transcription factor CncC was involved in the MdIS1-mediated regulation. This study highlights the common involvement of the CncC pathway in adaptive evolution, and provides an interesting case supportive of parallel evolution in P450-mediated insecticide resistance in insects.

Molecular Evolutionary Mechanisms of CYP6ER1vA-Type Variant Associated with Resistance to Neonicotinoid Insecticides in Field Populations of Nilaparvata lugens

The evolution of insecticide resistance has threatened the control of Nilaparvata lugens. Research on mechanisms behind neonicotinoid resistance in N. lugens remains incomplete. This study examined P450-mediated resistance to neonicotinoids in a resistant N. lugens strain (XA-2017-3G). The overexpression of CYP6ER1 in the XA-2017-3G strain plays a role in neonicotinoid resistance, as confirmed by RNA interference. Phenotypic analyses of CYP6ER1-mediated resistance in strains, including laboratory-susceptible, field-collected, and imidacloprid-laboratory further-selected strains, revealed that the vA-type/vL-type genotype exhibited greater resistance to neonicotinoids compared to the vA-type/vA-type genotype. The mRNA expression levels of CYP6ER1vA-type were closely correlated with the levels of neonicotinoid resistance in N. lugens strains, in which CYP6ER1vA-type overexpression is in part attributed to increased copy numbers of CYP6ER1. CYP6ER1vA-type-mediated neonicotinoid resistance was further confirmed by a CYP6ER1vA-type transgenic Drosophila melanogaster line. Taken together, our findings strongly suggest that the overexpression of CYP6ER1vA-type, which can be partially attributed to copy number variations, plays a crucial role in N. lugens resistance to neonicotinoids.

Overexpression of the F116V allele of CYP9A186 in transgenic Helicoverpa armigera confers high-level resistance to emamectin benzoate

Insect cytochrome P450s play important roles in the detoxification of xenobiotics and the metabolic resistance to insecticides. However, the approach for in vivo validation of the contribution of specific candidate P450s to resistance is still limited in most non-model insect species. Previous studies with heterologous expression and in vitro functional assays have confirmed that a natural substitution (F116V) in the substrate recognition site 1 (SRS1) of the CYP9A186 of Spodoptera exigua is a gain-of-function mutation, which results in detoxification capability of and thus high-level resistance to both emamectin benzoate (EB) and abamectin. In this study, we established an effective piggyBac-based transformation system in the serious agricultural pest Helicoverpa armigera and overexpressed in vivo a resistance P450 allele, CYP9A186-F116V, from another lepidopteran pest Spodoptera exigua. Bioassays showed that transgenic H. armigera larvae expressing CYP9A186-F116V obtained 358-fold and 38.6-fold resistance to EB and abamectin, respectively. In contrast, a transgenic line of Drosophila melanogaster overexpressing this P450 variant only confers ∼20-fold resistance to the two insecticides. This bias towards the resistance level revealed that closely related species might provide a more appropriate cellular environment for gene expression and subsequent toxicokinetics of insecticides. These results not only present an alternative method for in vivo functional characterization of P450s in H. armigera and other phylogenetically close species but also provide a valuable genetic engineering toolkit for the genetic manipulation of H. armigera.

Reverse genetic study reveals the molecular targets of chordotonal organ TRPV channel modulators

Insecticides have been widely used for the control of insect pests that have a significant impact on agriculture and human health. A better understanding of insecticide targets is needed for effective insecticide design and resistance management. Pymetrozine, afidopyropen and flonicamid are reported to target on proteins that located on insect chordotonal organs, resulting in the disruption of insect coordination and the inhibition of feeding. In this study, we systematically examined the susceptibility of six Drosophila melanogaster mutants (five transient receptor potential channels and one mechanoreceptor) to three commercially used insecticides, in order to identify the receptor subunits critical to the insect's response to insecticides. Our results showed that iav1, nan36aand wtrw1 mutants exhibited significantly reduced susceptibility to pymetrozine and afidopyropen, but not to flonicamid. The number of eggs produced by the three mutant females were significantly less than that of the w1118 strain. Meanwhile, the longevity of all male mutants and females of nan36a and wtrw1 mutants was significantly shorter than that of the w1118 strain as the control. However, we observed no gravitaxis defects in wtrw1 mutants and the anti-gravitaxis of wtrw1 mutants was abolished by pymetrozine. Behavioral assays using thermogenetic tools further confirmed the bioassay results and supported the idea that Nan as a TRPV subfamily member located in Drosophila chordotonal neurons, acting as a target of pymetrozine, which interferes with Drosophila and causes motor deficits with gravitaxis defects. Taken together, this study elucidates the interactions of pymetrozine and afidopyropen with TRPV channels, Nan and Iav, and TRPA channel, Wtrw. Our research provides another evidence that pymetrozine and afidopyropen might target on nan, iav and wtrw channels and provides insights into the development of sustainable pest management strategies.

Metal-Organic Framework-Based Insecticide and dsRNA Codelivery System for Insecticide Resistance Management

Targeted silencing of resistance-associated genes by specific double-stranded RNA (dsRNA) is an attractive strategy for overcoming insecticide resistance in insect pests. However, silencing target genes of insect pests by feeding on dsRNA transported via plants remains challenging. Herein, a codelivery system of insecticide and dsRNA is designed by encapsulating imidacloprid and dsNlCYP6ER1 within zeolitic imidazolate framework-8 (ZIF-8) nanoparticles to improve the susceptibility of Nilaparvata lugens (Stål) to imidacloprid. With an average particle size of 195 nm and a positive surface charge, the derived imidacloprid/dsNlCYP6ER1@ZIF-8 demonstrates good monodispersity. Survival curve results showed that the survival rates of N. lugens treated with imidacloprid and imidacloprid@ZIF-8 were 82 and 62%, respectively, whereas, in the imidacloprid/dsNlCYP6ER1@ZIF-8 treatment group, the survival rate of N. lugens is only 8%. Pot experiments demonstrate that the survival rate in the imidacloprid/dsNlCYP6ER1@ZIF-8 treatment group was much lower than that in the imidacloprid treatment group, decreasing from 54 to 24%. The identification of NlCYP6ER1 expression and the fluorescence tracking of ZIF-8 demonstrate that ZIF-8 can codeliver dsRNA and insecticide to insects via rice. Safety evaluation results showed that the dsNlCYP6ER1@ZIF-8 nanoparticle had desirable biocompatibility and biosafety to silkworm. This dsRNA and insecticide codelivery system may be extended to additional insecticides with potential resistance problems in the future, greatly enhancing the development of pest resistance management.

Divergent amplifications of CYP9A cytochrome P450 genes provide two noctuid pests with differential protection against xenobiotics

Here, we provide mechanistic support for the involvement of the CYP9A subfamily of cytochrome P450 monooxygenases in the detoxification of host plant defense compounds and chemical insecticides in Spodoptera exigua and Spodoptera frugiperda. Our comparative genomics shows that a large cluster of CYP9A genes occurs in the two species but with significant differences in its contents, including several species-specific duplicates and substantial sequence divergence, both between orthologs and between duplicates. Bioassays of CRISPR-Cas9 knockouts of the clusters show that, collectively, the CYP9As can detoxify two furanocoumarin plant defense compounds (imperatorin and xanthotoxin) and insecticides representing three different chemotypes (pyrethroids, avermectins, and oxadiazines). However, in vitro metabolic assays of heterologously expressed products of individual genes show several differences between the species in the particular CYP9As with activities against these compounds. We also find that the clusters show tight genetic linkage with high levels of pyrethroid resistance in field strains of the two species. We propose that their divergent amplifications of the CYP9A subfamily have not only contributed to the development of the broad host ranges of these species over long evolutionary timeframes but also supplied them with diverse genetic options for evolving resistance to chemical insecticides in the very recent past.

Structural modeling of cytochrome P450 51 from a deep-sea fish points to a novel structural feature in other CYP51s

Cytochromes P450 (CYP), enzymes involved in the metabolism of endogenous and xenobiotic substrates, provide an excellent model system to study how membrane proteins with unique functions have catalytically adapted through evolution. Molecular adaptation of deep-sea proteins to high hydrostatic pressure remains poorly understood. Herein, we have characterized recombinant cytochrome P450 sterol 14α-demethylase (CYP51), an essential enzyme of cholesterol biosynthesis, from an abyssal fish species, Coryphaenoides armatus. C. armatus CYP51 was heterologously expressed in Escherichia coli following N-terminal truncation and purified to homogeneity. Recombinant C. armatus CYP51 bound its sterol substrate lanosterol giving a Type I binding spectra (KD 15 μM) and catalyzed lanosterol 14α-demethylation turnover at 5.8 nmol/min/nmol P450. C. armatus CYP51 also bound the azole antifungals ketoconazole (KD 0.12 μM) and propiconazole (KD 0.54 μM) as determined by Type II absorbance spectra. Comparison of C. armatus CYP51 primary sequence and modeled structures with other CYP51s identified amino acid substitutions that may confer an ability to function under pressures of the deep sea and revealed heretofore undescribed internal cavities in human and other non-deep sea CYP51s. The functional significance of these cavities is not known. PROLOGUE: This paper is dedicated in memory of Michael Waterman and Tsuneo Omura, who as good friends and colleagues enriched our lives. They continue to inspire us.

CRISPR/Cas9-mediated knockout of NlCYP6CS1 gene reveals its role in detoxification of insecticides in Nilaparvata lugens (Hemiptera: Delphacidae)

BACKGROUND

The brown planthopper (Nilaparvata lugens) is one of the major rice insect pests in Asia. Recently, high levels of insecticide resistance have been frequently reported and cytochrome P450 monooxygenase (P450)-mediated metabolic detoxification is a common resistance mechanism in N. lugens. However, there has been no persuasive genetic method to prove the role of P450s in insecticide resistance in N. lugens.

RESULTS

Here, CRISPR/Cas9 system was used to disrupt the P450 gene NlCYP6CS1 to elucidate its role in insecticide resistance in field populations of N. lugens. We successfully constructed a homozygous strain (Nl6CS1-KO) with a 5-bp deletion and 1-bp insertion mutation of NlCYP6CS1. Compared with a background resistant strain (Nl-R), the susceptibility of knockout strain Nl6CS1-KO to imidacloprid, nitenpyram, thiamethoxam, dinotefuran, and pymetrozine was increased by 2.3-, 3.4-, 7.0-, 4.2- and 3.9-fold, respectively, but not significantly changed to triflumezopyrim, chlorpyrifos and buprofezin. Life table analysis demonstrated that the Nl6CS1-KO strain resembled the Nl-R strain in terms of egg and nymph developmental duration and adult lifespan, but differed from the Nl-R strain in the survival rate of eggs and nymphs, reproduction, and body weight.

CONCLUSIONS

Our study demonstrates the effect of functional deletion of NlCYP6CS1 on multiple insecticide resistance in N. lugens. For the first time, we applied CRISPR/Cas9 system to reveal the mechanism of insecticide resistance in N. lugens, which may shed light on similar studies in other hemipteran insects. © 2023 Society of Chemical Industry.

Development of an efficient insecticide substrate and inhibitor screening system of insect P450s using fission yeast

Metabolic resistance is one of the most frequent mechanisms of insecticide resistance, characterized by an increased expression of several important enzymes and transporters, especially cytochrome P450s (CYPs). Due to the large number of P450s in pests, determining the precise relationship between these enzymes and the insecticide substrates is a challenge. Herein, we developed a luminescence-based screening system for efficient identification of insecticide substrates and insect P450 inhibitors. We recombinantly expressed Bemisia tabaci CYP6CM1vQ (Bt CYP6CM1vQ) in the fission yeast Schizosaccharomyces pombe and subsequently permeabilized the yeast cells to convert them into "enzyme bags". We exploited these enzyme bags to screen the activity of twelve luciferin substrates and identified Luciferin-FEE as the optimal competing probe that was further used to characterize the metabolism of eight candidate commercial insecticides. Among them, Bt CYP6CM1vQ exhibited notable activity against pymetrozine and imidacloprid. Their binding modes were predicted by homology modeling and molecular docking, revealing the mechanisms of the metabolism. We also tested the inhibitory effect of eight known P450 inhibitors using our system and identified letrozole and 1-benzylimidazole as showing significant activity against Bt CYP6CM1vQ, with IC₅₀ values of 23.74 μM and 1.30 μM, respectively. Their potential to be developed as an insecticide synergist was further proven by an in vitro toxicity assay using imidacloprid-resistant Bemisia tabaci. Overall, our luciferin-based enzyme bag method is capable of providing a robust and efficient screening of insect P450 substrates and, more importantly, inhibitors to overcome the resistance.

A single point mutation in the Bemisia tabaci cytochrome-P450 CYP6CM1 causes enhanced resistance to neonicotinoids

The tobacco whitefly, Bemisia tabaci, is a polyphagous crop pest which causes high levels of economic damage across the globe. Insecticides are often required for the effective control of this species, among which the neonicotinoid class have been particularly widely used. Deciphering the mechanisms responsible for resistance to these chemicals is therefore critical to maintain control of B. tabaci and limit the damage it causes. An important mechanism of resistance to neonicotinoids in B. tabaci is the overexpression of the cytochrome P450 gene CYP6CM1 which leads to the enhanced detoxification of several neonicotinoids. In this study we show that qualitative changes in this P450 dramatically alter its metabolic capacity to detoxify neonicotinoids. CYP6CM1 was significantly over-expressed in two strains of B. tabaci which displayed differing levels of resistance to the neonicotinoids imidacloprid and thiamethoxam. Sequencing of the CYP6CM1 coding sequence from these strains revealed four different alleles encoding isoforms carrying several amino acid changes. Expression of these alleles in vitro and in vivo provided compelling evidence that a mutation (A387G), present in two of the CYP6CM1 alleles, results in enhanced resistance to several neonicotinoids. These data demonstrate the importance of both qualitative and quantitative changes in genes encoding detoxification enzymes in the evolution of insecticide resistance and have applied implications for resistance monitoring programs.

Metabolic functional redundancy of the CYP9A subfamily members leads to P450-mediated lambda-cyhalothrin resistance in Cydia pomonella

BACKGROUND

The evolution of insect resistance to pesticides poses a continuing threat to sustainable pest management. While much is known about the molecular mechanisms that confer resistance in model insects and few agricultural pests, far less is known about fruit pests. Field-evolved resistance to synthetic insecticides such as lambda-cyhalothrin has been widely documented in Cydia pomonella, a major invasive pest of pome fruit worldwide, and the increased production of cytochrome P450 monooxygenases (P450s) has been linked to resistance in field-evolved resistant populations. However, the underlying molecular mechanisms of P450-mediated insecticide resistance remain largely unknown.

RESULTS

Here we found that functional redundancy and preference of metabolism by P450s genes in the CYP9A subfamily confer resistance to lambda-cyhalothrin in Cydia pomonella. A total of four CYP9A genes, including CYP9A61, CYP9A120, CYP9A121, and CYP9A122, were identified from Cydia pomonella. Among these, CYP9A120, CYP9A121, and CYP9A122 were predominantly expressed in the midgut of larvae. The expression levels of these P450 genes were significantly induced by a lethal dose that would kill 10% (LD₁₀ ) of lambda-cyhalothrin and were overexpressed in a field-evolved lambda-cyhalothrin resistant population. Knockdown of CYP9A120 and CYP9A121 by RNA-mediated interference (RNAi) increased the susceptibility of larvae to lambda-cyhalothrin. In vitro assays demonstrated that recombinant P450s expressed in Sf9 cells can metabolize lambda-cyhalothrin, but with functional redundancy and divergence through regioselectivity of metabolism. CYP9A121 preferred to convert lambda-cyhalothrin to 2'-hydroxy-lambda-cyhalothrin, whereas CYP9A122 only generated 4'-hydroxy metabolite of lambda-cyhalothrin. Although possesses a relatively low metabolic capability, CYP9A120 balanced catalytic competence to generate both 2'- and 4'-metabolites.

CONCLUSION

Collectively, these results reveal that metabolic functional redundancy of three members of the CYP9A subfamily leads to P450-mediated lambda-cyhalothrin resistance in Cydia pomonella, thus representing a potential adaptive evolutionary strategy during its worldwide expansion. © 2022 Society of Chemical Industry.

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.

Genetics of tolerance in honeybees to the neonicotinoid clothianidin

The effects of neonicotinoid insecticides (NNIs) on honeybee health are intensely debated, with numerous studies showing negative effects of exposure, while others report no such effects. We carried out experiments to study the genetic and molecular basis of NNI tolerance in honeybees, which may underlie the discrepancies observed in the literature. We discovered that worker survival post-exposure to an acute oral dose of clothianidin is heritable (H 2 = 37.8%). Tolerance to clothianidin was not associated with differences in the expression of detoxification enzymes in our experiments. Instead, mutations in the primary neonicotinoid detoxification genes CYP9Q1 and CYP9Q3 were strongly associated with worker survival post-clothianidin exposure. In some instances, the strong association between CYP9Q haplotypes and worker survival was associated with the protein's predicted binding affinity for clothianidin. Our findings have implications regarding future toxicological studies utilizing honeybees as a model pollinator.

Importance of CYP6ER1 Was Different among Neonicotinoids in Their Susceptibility in Nilaparvata lugens

CYP6ER1 overexpression is a prevalent mechanism for neonicotinoid resistance in Nilaparvata lugens. Except for imidacloprid, the metabolism of other neonicotinoids by CYP6ER1 lacked direct evidence. In this study, a CYP6ER1 knockout strain (CYP6ER1^-/-) was constructed using the CRISPR/Cas9 strategy. The CYP6ER1^-/- strain showed much higher susceptibility to imidacloprid and thiacloprid with an SI (sensitivity index, LC₅₀ of WT/LC₅₀ of CYP6ER1^-/-) of over 100, which was 10-30 for four neonicotinoids (acetamiprid, nitenpyram, clothianidin, and dinotefuran) and less than 5 for flupyradifurone and sulfoxaflor. Recombinant CYP6ER1 showed the highest activity to metabolize imidacloprid and thiacloprid and moderate activity for the other four neonicotinoids. Main metabolite identification and oxidation site prediction revealed that CYP6ER1 activities were insecticide structure-dependent. The most potential oxidation site of imidacloprid and thiacloprid was located in the five-membered heterocycle with hydroxylation activity. For the other four neonicotinoids, the potential site was within the ring opening of a five-membered heterocycle, indicating N-desmethyl activity.

Alternative Splicing and Expression Reduction of P450 Genes Mediating the Oxidation of Chlorpyrifos Revealed a Novel Resistance Mechanism in Nilaparvata lugens

Cytochrome P450 enzymes metabolize various xenobiotics in insects. Compared to numerous P450s associated with insecticide detoxification and resistance, fewer have been identified to bioactivate proinsecticides in insects. Here we reported that two P450s, CYP4C62 and CYP6BD12, in Nilaparvata lugens could bioactivate chlorpyrifos, an organophosphorus insecticide, into its active ingredient chlorpyrifos-oxon in vivo and in vitro. RNAi knockdown of these two genes significantly reduced the sensitivity to chlorpyrifos and the formation of chlorpyrifos-oxon in N. lugens. Chlorpyrifos-oxon was generated when chlorpyrifos was incubated with the crude P450 enzyme prepared from N. lugens or recombinant CYP4C62 and CYP6BD12 enzymes. The expression reduction of CYP4C62 and CYP6BD12 and alternative splicing in CYP4C62 reduced the oxidation of chlorpyrifos into chlorpyrifos-oxon, which contributed importantly to chlorpyrifos resistance in N. lugens. This study revealed a novel mechanism of insecticide resistance due to the bioactivation reduction, which would be common for all currently used proinsecticides.

Microbiome variation correlates with the insecticide susceptibility in different geographic strains of a significant agricultural pest, Nilaparvata lugens

Microbiome-mediated insecticide resistance is an emerging phenomenon found in insect pests. However, microbiome composition can vary by host genotype and environmental factors, but how these variations may be associated with insecticide resistance phenotype remains unclear. In this study, we compared different field and laboratory strains of the brown planthopper Nilaparvata lugens in their microbiome composition, transcriptome, and insecticide resistance profiles to identify possible patterns of correlation. Our analysis reveals that the abundances of core bacterial symbionts are significantly correlated with the expression of several host detoxifying genes (especially NlCYP6ER1, a key gene previously shown involved in insecticides resistance). The expression levels of these detoxifying genes correlated with N. lugens insecticide susceptibility. Furthermore, we have identified several environmental abiotic factors, including temperature, precipitation, latitude, and longitude, as potential predictors of symbiont abundances associated with expression of key detoxifying genes, and correlated with insecticide susceptibility levels of N. lugens. These findings provide new insights into how microbiome-environment-host interactions may influence insecticide susceptibility, which will be helpful in guiding targeted microbial-based strategies for insecticide resistance management in the field.

The overexpression of cytochrome P450 genes confers buprofezin resistance in the brown planthopper, Nilaparvata lugens (Stål)

BACKGROUND

Buprofezin, an insect growth regulator, has been widely used to control brown planthopper (BPH), Nilaparvata lugens, one of the most destructive pests of rice crops in Asia. The intensive use of this compound has resulted in very high levels of resistance to buprofezin in the field, however, the underpinning mechanisms of resistance have not been fully resolved.

RESULTS

Insecticide bioassays using the P450 inhibitor piperonyl butoxide significantly synergized the toxicity of buprofezin in two resistant strains of BPH (BPR and YC2017) compared to a susceptible strain (Sus), suggesting P450s play a role in resistance to this compound. Whole transcriptome profiling identified 1110 genes that were upregulated in the BPR strain compared to the Sus strain, including 13 cytochrome P450 genes, eight esterases and one glutathione S-transferase. Subsequently, qPCR validation revealed that four of the P450 genes, CYP6ER1vA, CYP6CW1, CYP4C77, and CYP439A1 were significantly overexpressed in both the BRP and YC2017 strains compared with the Sus strain. Further functional analyses showed that only suppression of CYP6ER1vA, CYP6CW1, and CYP439A1 gene expression by RNA interference significantly increased the toxicity of buprofezin against BPH. However, only transgenic Drosophila melanogaster expressing CYP6ER1vA and CYP439A1 exhibited significant resistance to buprofezin. Finally, the BPR strain was found to exhibit modest but significant levels of resistance to acetamiprid, dinotefuran and pymetrozine.

CONCLUSIONS

Our findings provide strong evidence that the overexpression of CYP6ER1vA and CYP439A1 contribute to buprofezin resistance in BPH, and that resistance to this compound is associated with low-level resistance to acetamiprid, dinotefuran and pymetrozine. These results advance understanding of the molecular basis of BPH resistance to buprofezin and will inform the development of management strategies for the control of this highly damaging pest. © 2022 Society of Chemical Industry.

Pesticide contamination in agro-ecosystems: toxicity, impacts, and bio-based management strategies

Continuous rise in application of pesticides in the agro-ecosystems in order to ensure food supply to the ever-growing population is of greater concern to the human health and the environment. Once entered into the agro-ecosystem, the fate and transport of pesticides is determined largely by the nature of pesticides and the soil attributes, in addition to the soil-inhabiting microbes, fauna, and flora. Changes in the soil microbiological actions, soil properties, and enzymatic activities resulting from pesticide applications are the important factors substantially affecting the soil productivity. Disturbances in the microbial community composition may lead to the considerable perturbations in cycling of major nutrients, metals, and subsequent uptake by plants. Indiscriminate applications are linked with the accumulation of pesticides in plant-based foods, feeds, and animal products. Furthermore, rapid increase in the application of pesticides having long half-life has also been reported to contaminate the nearby aquatic environments and accumulation in the plants, animals, and microbes surviving there. To circumvent the negative consequences of pesticide application, multitude of techniques falling in physical, chemical, and biological categories are presented by different investigators. In the present study, important findings pertaining to the pesticide contamination in cultivated agricultural soils; toxicity on soil microbes, plants, invertebrates, and vertebrates; effects on soil characteristics; and alleviation of toxicity by bio-based management approaches have been thoroughly reviewed. With the help of bibliometric analysis, thematic evolution and research trends on the bioremediation of pesticides in the agro-ecosystems have also been highlighted.

Functional Validation of the Roles of Cytochrome P450s in Tolerance to Thiamethoxam and Imidacloprid in a Field Population of Aphis gossypii

Field populations of Aphis gossypii (SDR) have evolved high resistance to neonicotinoids, including thiamethoxam and imidacloprid. Synergism bioassays and transcriptomic comparison of the SDR and susceptible (SS) strains revealed that the cytochrome P450s may contribute to the neonicotinoid resistance evolution. The transcripts of some P450s were constitutively overexpressed in the SDR strain, and many genes showed expression plasticity under insecticide exposure. Drosophila that ectopically expressed CYPC6Y9, CYP4CK1, CYP6DB1, and CYP6CZ1 showed greater resistance (>8.0-fold) to thiamethoxam, and Drosophila that expressed CYPC6Y9, CYP6CY22, CYP6CY18, and CYP6D subfamily genes showed greater resistance (>5-fold) to imidacloprid. Five P450 genes that caused thiamethoxam resistance also conferred resistance to α-cypermethrin. Furthermore, the knockdown of CYP4CK1, CYP6CY9, CYP6CY18, CYPC6Y22, CYP6CZ1, and CYP6DB1 dramatically increased the sensitivity of the SDR strain to thiamethoxam or imidacloprid. These results indicate the involvement of multiple P450 genes, rather than one key gene, in neonicotinoid resistance in field populations.

Functional Inquiry into ATP-Binding Cassette Transporter Genes Contributing to Spirotetramat Resistance in Aphis gossypii Glover

ATP-binding cassette (ABC) transporters regulate the efflux of a broad spectrum of substrates to extracellular transporting, which play an important role in the detoxification process in arthropods. Here, we described a comprehensive approach to explore the involvement of ABC transporters in spirotetramat resistance in cotton aphids. In this study, synergism bioassays showed 17.05% and 35.42% increases in the toxicity to spirotetramat with the ABC inhibitor verapamil in adult and 3rd instar nymph aphids of the SR strain, respectively. In a competitive assay based on the microinjection of a fluorescent ABC transporter substrate, verapamil (a general ABC inhibitor) and spirotetramat significantly inhibited the elimination of Texas Red. Based on transcriptome data of midguts of spirotetramat-susceptible (SS) and -resistant (SR) strains, the expression levels of ABCB4, ABCB5, ABCF2, MRP11, and MRP12 were significantly upregulated in the SR strain midgut compared to that of the SS strain. Gene functional analysis based on ectopic expression and RNA interference (RNAi) proved that ABCB4, ABCB5, ABCF2, MRP11, and MRP12 were involved in the tolerance of cotton aphids to spirotetramat. Moreover, the upregulated ABCF2, ABCB4, and ABCB5 in the midgut of the SR strain contributed more to the resistance of spirotetramat in in vitro functional analysis. In summary, these results demonstrate that candidate ABC transporter genes in the midgut tissue were involved in spirotetramat resistance, which will help reveal the relationship between ABC transporters and the development of spirotetramat resistance in field populations.

Sequential phase I metabolism of pyrethroids by duplicated CYP6P9 variants results in the loss of the terminal benzene moiety and determines resistance in the malaria mosquito Anopheles funestus

Pyrethroid resistance in Anopheles funestus is threatening the eradication of malaria. One of the major drivers of pyrethroid resistance in An. funestus are cytochrome P450 monooxygenases CYP6P9a and CYP6P9b, which are found upregulated in resistant An. funestus populations from Sub-Saharan Africa and are known to metabolise pyrethroids. Here, we have functionally expressed CYP6P9a and CYP6P9b variants and investigated their interactions with azole-fungicides and pyrethroids. Some azole fungicides such as prochloraz inhibited CYP6P9a and CYP6P9b at nanomolar concentrations, whereas pyrethroids were weak inhibitors (>100 μM). Amino acid sequence comparisons suggested that a valine to isoleucine substitution at position 310 in the active site cavity of CYP6P9a and CYP6P9b, respectively, might affect substrate binding and metabolism. We therefore swapped the residues by site directed mutagenesis to produce CYP6P9a^I310V and CYP6P9b^V310I. CYP6P9b^V310I produced stronger metabolic activity towards coumarin substrates and pyrethroids, particularly permethrin. The V310I mutation was previously also detected in a pyrethroid resistant field population of An. funestus in Benin. Additionally, we found the first metabolite of permethrin and deltamethrin after hydroxylation, 4'OH permethrin and 4'OH deltamethrin, were also suitable substrates for CYP6P9-variants, and were depleted by both enzymes to a higher extent than as their respective parent compounds (approximately 20% more active). Further, we found that both metabolites were toxic against An. funestus FANG (pyrethroid susceptible) but not towards FUMOZ-R (pyrethroid resistant) mosquitoes, the latter suggesting detoxification by overexpressed CYP6P9a and CYP6P9b. We confirmed by mass-spectrometric analysis that CYP6P9a and CYP6P9b are capable of cleaving phenoxybenzyl-ethers in type I pyrethroid permethrin and type II pyrethroid deltamethrin and that both enzymes preferentially metabolise trans-permethrin. This provides new insight into the metabolism of pyrethroids and a greater understanding of the molecular mechanisms of pyrethroid resistance in An. funestus.

The Cross-Resistance Pattern and the Metabolic Resistance Mechanism of Acetamiprid in the Brown Planthopper, Nilaparvata lugens (Stål)

Acetamiprid is widely used in paddy fields for controlling Nilaparvata lugens (Stål). However, the risk of resistance development, the cross-resistance pattern and the resistance mechanism of acetamiprid in this pest remain unclear. In this study, an acetamiprid-resistant strain (AC-R) was originated from a field strain (UNSEL) through successive selection with acetamiprid for 30 generations, which reached 60.0-fold resistance when compared with a laboratory susceptible strain (AC-S). The AC-R strain (G₃₀) exhibited cross-resistance to thiamethoxam (25.6-fold), nitenpyram (21.4-fold), imidacloprid (14.6-fold), cycloxaprid (11.8-fold), dinotefuran (8.7-fold), sulfoxaflor (7.6-fold) and isoprocarb (8.22-fold), while there was no cross-resistance to etofenprox, buprofezin and chlorpyrifos. Acetamiprid was synergized by the inhibitor piperonyl butoxide (2.2-fold) and the activity of cytochrome P450 monooxygenase was significantly higher in the AC-R strain compared with the AC-S strain, suggesting the critical role of P450. The gene expression results showed that the P450 gene CYP6ER1 was significantly overexpressed in AC-R compared with the AC-S and UNSEL strains. In addition, the RNA interference (RNAi) of CYP6ER1 significantly increased the susceptibility of AC-R to acetamiprid. Molecular docking predicted that acetamiprid and CYP6ER1 had close binding sites, and the nitrogen atoms had hydrogen bond interactions with CYP6ER1. These results demonstrated that the overexpression of CYP6ER1 contributed to acetamiprid resistance in N. lugens.

Heritable Epigenomic Modifications Influence Stress Resilience and Rapid Adaptations in the Brown Planthopper (Nilaparvata lugens)

DNA methylation in insects is integral to cellular differentiation, development, gene regulation, genome integrity, and phenotypic plasticity. However, its evolutionary potential and involvement in facilitating rapid adaptations in insects are enigmatic. Moreover, our understanding of these mechanisms is limited to a few insect species, of which none are pests of crops. Hence, we studied methylation patterns in the brown planthopper (BPH), a major rice pest, under pesticide and nutritional stress, across its life stages. Moreover, as the inheritance of epigenetic changes is fundamentally essential for acclimation, adaptability, and evolution, we determined the heritability and persistence of stress-induced methylation marks in BPH across generations. Our results revealed that DNA methylation pattern(s) in BPH varies/vary with environmental cues and is/are insect life-stage specific. Further, our findings provide novel insights into the heritability of stress-induced methylation marks in BPH. However, it was observed that, though heritable, these marks eventually fade in the absence of the stressors, thereby suggesting the existence of fitness cost(s) associated with the maintenance of the stressed epigenotype. Furthermore, we demonstrate how 5-azacytidine-mediated disruption of BPH methylome influences expression levels of stress-responsive genes and, thereby, highlight demethylation/methylation as a phenomenon underlying stress resilience of BPH.

Effects of nicotinic acetylcholine receptor subunit deletion mutants on insecticide susceptibility and fitness in Drosophila melanogaster

BACKGROUND

Nicotinic acetylcholine receptors (nAChRs) are major excitatory neurotransmitter receptors in insects and also the target site for many insecticides. Unfortunately, the effectiveness of these insecticides is diminishing as a consequence of the evolution of insecticide resistance. Further exploration of insecticide targets is important to sustainable pest management.

RESULTS

In order to validate the role of nAChR subunits in insecticide susceptibility and test whether the subunit's absence imposes the fitness cost on insects, we determined the susceptibility of eight nAChR subunit deletion mutants of Drosophila melanogaster to nine insecticides. These findings highlighted the specific resistance of the Dα6 deletion mutant to spinosyns. Although triflumezopyrim, dinotefuran and imidacloprid are competitive modulators of nAChRs, differences in susceptibility of the insect with different deletion mutants suggested that the target sites of these three insecticides do not overlap completely. Mutants showed decreased susceptibility to insecticides, accompanied by a reduction in fitness. The number of eggs produced by Dα1^attP , Dα2^attP , Dβ2^attP and Dβ3^attP females was significantly lesser than that of the vas-Cas9 strain as the control. In addition, adults of Dα2^attP , Dα3^attP and Dα7^attP strains showed lower climbing performance. Meanwhile, males of Dα3^attP , Dα5^attP , Dβ2^attP and Dβ3^attP , and females of Dβ2^attP showed significantly shorter longevity than those of the vas-Cas9 strain.

CONCLUSION

This study provides new insights into the interactions of different insecticides with different nAChRs subunit in D. melanogaster as a research model, it could help better understand such interaction in agricultural pests whose genetic manipulations for toxicological research are often challenging. © 2022 Society of Chemical Industry.

Single amino acid variations drive functional divergence of cytochrome P450s in Helicoverpa species

Divergence of gene function is a hallmark of evolution, but assessing such divergence in one species or between species requires information on functional alterations of the alleles and homologs. Here, we explore the functional divergence of two paralogs, CYP6AE19 and CYP6AE20, from Helicoverpa armigera, and two close orthologs, CYP6B8 and CYP6B7, from two related species (Helicoverpa zea and H. armigera); although there is high sequence identity within each pair of enzymes, the latter P450 of each pair has lost metabolic competence towards the plant allelochemical xanthotoxin. Multiple chimeric and single/double site mutants were created by exchanging the diverse substrate recognition sites (SRSs) and amino acids within each pair of P450s. Heterologous expression in Sf9 cells and in vitro metabolism studies showed that the exchange of SRS4 swapped the activity of CYP6AE19 and CYP6AE20, and subsequent site-directed mutagenesis demonstrated that the CYP6AE20 V318M substitution causes a gain-of-function towards xanthotoxin. Meanwhile, a single amino acid substitution (L489P) in SRS6 was found to swap activity between the CYP6B orthologs. Sequence alignments of CYP6AE paralogs and all reported insect xanthotoxin-metabolizing P450s suggest M318 and P489 are essential for the catalytic activities of CYP6AE paralogs and CYP6B orthologs, respectively, but P450s in different subfamilies may have different mechanisms towards the same substrate. Our findings demonstrate that a single amino acid substitution can suffice to alter substrate metabolism and this functional divergence resulting from natural mutations will help to further our understanding of the process of natural selection of P450 genes and their role in insect-host plant interactions.

Cross effects of heat stress and three insecticides on the survival of the codling moth Cydia pomonella (L.): Investigating the molecular and biochemical mechanisms

As temperature is expected to strongly increase in the future, understanding temperature-mediated toxicity of insecticides is determinant to assess pest management efficiency in a warming world. Investigating molecular and biochemical mechanisms associated with cross mechanisms of temperature and insecticides on pests' tolerance would also be useful in this context. This study aimed to investigate cross effects between temperature and insecticides on the survival of a major pest, the codling moth Cydia pomonella, and their underlying mechanisms. The effect of three insecticidal active ingredients, i.e. chlorantraniliprole, emamectin and spinosad, was assessed at different temperatures on: (i) C. pomonella larval survival; (ii) detoxification enzymes activities (cytochrome P450 multi-function oxygenases, carboxylesterases and glutathione S-transferases) and (iii) genes expression of some detoxification enzymes, heat shock proteins and receptors targeted by the insecticides. We observed a decreased efficiency of emamectin and spinosad at high temperature to control the codling moth while no influence of temperature on chlorantraniliprole efficacy was observed. Detoxification enzymes activities were improved by heat stress alone but not by double stress (temperature + insecticides). Moreover, two detoxification genes (Cyp9A61 and Gst1) were over-expressed by a single stress but not by two stresses while Hsp70 and Cyp6B2 genes may be involved in tolerance to two stresses in C. pomonella. These results confirmed the cross effects of temperature and insecticides on C. pomonella for emamectin and spinosad and provided clues to understand how temperature affects the susceptibility of C. pomonella to insecticides. They illustrate however the complexity of molecular and biochemical responses of individuals facing multiple stresses.

Adipokinetic hormone signaling regulates cytochrome P450-mediated chlorantraniliprole sensitivity in Spodoptera frugiperda (Lepidoptera: Noctuidae)

BACKGROUND

Fall armyworm (FAW, Spodoptera frugiperda) is one of the most destructive and invasive pests worldwide and causes significant economic losses. Intensive and frequent use of insecticides has led to the development of resistance in FAW. Adipokinetic hormone (AKH) have been proven to be involved in insecticide resistance in insects. However, the molecular mechanism underlying chlorantraniliprole resistance mediated by AKH signaling in FAW remains unclear.

RESULTS

The expression of SpfAKHR was highest in male adults and lowest in 1st instar larvae. SpfAKH was expressed the highest in eggs and the lowest in 6th instar larvae. AKH signaling was involved in the sensitivity of FAW to chlorantraniliprole through a toxicological bioassay, and the combination of chlorantraniliprole and bithionol (an inhibitor of key enzymes in the AKH pathway) significantly increased the mortality of FAW. Chlorantraniliprole significantly induced the expression of ten P450s, SpfAKH and SpfAKHR in FAW. RNA interference against SpfAKHR significantly decreased the P450 content, downregulated the expression of three P450 genes (SpfCYP6B50, SpfCYP321A9 and SpfCYP9A58) and inhibited the resistance of FAW to chlorantraniliprole. The topical application of AKH peptide significantly increased the P450 content, upregulated the expression of five P450 genes (SpfCYP321A9, SpfCY321A8, SpfCYP321A10, SpfCYP321A7 and SpfCYP6AB12), and enhanced the survival of FAW against chlorantraniliprole.

CONCLUSIONS

AKH plays an important role in enhancing chlorantraniliprole resistance in FAW by exerting a positive influence on P450 gene expression and P450 content. These results provide valuable insights into insecticide resistance regulation and FAW control strategies. © 2022 Society of Chemical Industry.

P450 gene duplication and divergence led to the evolution of dual novel functions and insecticide cross-resistance in the brown planthopper Nilaparvata lugens

The sustainable control of many highly damaging insect crop pests and disease vectors is threatened by the evolution of insecticide resistance. As a consequence, strategies have been developed that aim to prevent or delay resistance development by rotating or mixing insecticides with different modes of action (MoA). However, these approaches can be compromised by the emergence of mechanisms that confer cross-resistance to insecticides with different MoA. Despite the applied importance of cross-resistance, its evolutionary underpinnings remain poorly understood. Here we reveal how a single gene evolved the capacity to detoxify two structurally unrelated insecticides with different MoA. Using transgenic approaches we demonstrate that a specific variant of the cytochrome P450 CYP6ER1, previously shown to confer resistance to the neonicotinoid imidacloprid in the brown planthopper, N. lugens, also confers cross-resistance to the phenylpyrazole ethiprole. CYP6ER1 is duplicated in resistant strains, and we show that while the acquisition of mutations in two encoded substrate recognition sites (SRS) of one of the parologs led to resistance to imidacloprid, a different set of mutations, outside of known SRS, are primarily responsible for resistance to ethiprole. Epistatic interactions between these mutations and their genetic background suggest that the evolution of dual resistance from the same gene copy involved functional trade-offs in respect to CYP6ER1 catalytic activity for ethiprole versus imidacloprid. Surprisingly, the mutations leading to ethiprole and imidacloprid resistance do not confer the ability to detoxify the insecticide fipronil, another phenylpyrazole with close structural similarity to ethiprole. Taken together, these findings reveal how gene duplication and divergence can lead to the evolution of multiple novel functions from a single gene. From an applied perspective they also demonstrate how cross-resistance to structurally unrelated insecticides can evolve, and illustrate the difficulty in predicting cross-resistance profiles mediated by metabolic mechanisms.

The evolution of mechanisms that confer resistance to insecticides of different mode of action (MoA) threatens to undermine the primary strategy used to manage insecticide resistance based on the rotation of insecticides of different MoA. Despite its importance, our understanding of how cross-resistance evolves remains surprisingly poor. Here, we uncover the mechanisms by which a single gene, CYP6ER1, encoding a cytochrome P450 enzyme, evolved resistance to two insecticides of different MoA, namely the neonicotinoid imidacloprid and the phenylpyrazole ethiprole. We have previously shown that CYP6ER1 is duplicated in imidacloprid resistant strains of the brown planthopper, N. lugens. In the current study we demonstrate that two different sets of mutations, occurring in the same copy of CYP6ER1, were key to conferring the de novo ability to detoxify imidacloprid and ethiprole. We uncover complex interactions between these mutations and their genetic background that provide strong evidence of functional trade-offs associated with evolving dual resistance. Our findings provide fundamental and applied insights into the versatility of gene duplication in providing opportunities for functional innovation during the evolution of adaptive traits and reveal the mechanisms by which dual novel functions can arise from single genes.