Structure and Function of Cytochrome P450S in Insect Adaptation to Natural and Synthetic Toxins: Insights Gained from Molecular Modeling

Identification and Functional Characterization of Carboxylesterase Genes Involved in Spirodiclofen Resistance in Panonychus citri (McGregor)

Overexpression of carboxyl/cholinesterase (CCE) genes has been reported to be associated with many cases of pesticide resistance in arthropods. However, it has been rarely documented that CCE genes participate in spirodiclofen resistance in Panonychus citri. In previous research, we found that spirodiclofen resistance is related to increased P450 and CCE enzyme activities in P. citri. In this study, we identified two CCE genes, PcCCE3 and PcCCE5, which were significantly upregulated in spirodiclofen-resistant strain and after exposure to spirodiclofen. RNA interference of PcCCE3 and PcCCE5 increased the spirodiclofen susceptibility in P. citri. In vitro metabolism indicated that PcCCE3 and PcCCE5 could interact with spirodiclofen, but metabolites were detected only in the PcCCE3 treatment. Our results indicated that PcCCE3 participates in spirodiclofen resistance through direct metabolism, and PcCCE5 may be involved in the spirodiclofen resistance by passive binding and sequestration, which provides new insights into spirodiclofen resistance in P. citri.

Transcription Factors AhR and ARNT Regulate the Expression of CYP6SX1 and CYP3828A1 Involved in Insecticide Detoxification in Bradysia odoriphaga

Aryl hydrocarbon receptor (AhR) and aryl hydrocarbon receptor nuclear translocator (ARNT) mediate the responses of adaptive metabolism to various xenobiotics. Here, we found that BoAhR and BoARNT are highly expressed in the midgut of Bradysia odoriphaga larvae. The expression of BoAhR and BoARNT was significantly increased after exposure to imidacloprid and phoxim. The knockdown of BoAhR and BoARNT significantly decreased the expression of CYP6SX1 and CYP3828A1 as well as P450 enzyme activity and caused a significant increase in the sensitivity of larvae to imidacloprid and phoxim. Exposure to β-naphthoflavone (BNF) significantly increased the expression of BoAhR, BoARNT, CYP6SX1, and CYP3828A1 as well as P450 activity and decreased larval sensitivity to imidacloprid and phoxim. Furthermore, CYP6SX1 and CYP3828A1 were significantly induced by imidacloprid and phoxim, and the silencing of these two genes significantly reduced larval tolerance to imidacloprid and phoxim. Taken together, the BoAhR/BoARNT pathway plays key roles in larval tolerance to imidacloprid and phoxim by regulating the expression of CYP6SX1 and CYP3828A1.

BioProtIS: Streamlining protein-ligand interaction pipeline for analysis in genomic and transcriptomic exploration

The identification of protein-ligand interactions plays a pivotal role in elucidating biological processes and discovering potential bioproducts. Harnessing the capabilities of computational methods in drug discovery, we introduce an innovative Inverted Virtual Screening (IVS) pipeline. This pipeline Integrated molecular dynamics and docking analyses to ensure that protein structures are not only energetically favorable but also representative of stable conformations. The primary objective of this pipeline is to automate and streamline the analysis of protein-ligand interactions at both genomic and transcriptomic scales. In the contemporary post-genomic era, high-throughput computational screening for bioproducts, biological systems, and therapeutic drugs has become a cornerstone practice. This approach offers the promise of cost-effectiveness, time efficiency, and optimization of laboratory work. Nevertheless, a notable deficiency persists in the availability of efficient pipelines capable of automating the virtual screening process, seamlessly integrating input and output, and leveraging the full potential of open-source tools. To bridge this critical gap, we have developed a versatile pipeline known as BioProtIS. This tool seamlessly integrates a suite of state-of-the-art tools, including Modeller, AlphaFold, Gromacs, FPOCKET, and AutoDock Vina, thus facilitating the streamlined docking of ligands with an expansive repertoire of proteins sourced from genomes and transcriptomes, and substrates. To assess the pipeline's performance, we employed the transcriptomes of Cereus jamacaru (a cactus species) and Aspisoma lineatum (firefly), along with the genome of Homo sapiens. This integration not only improves the accuracy of ligand-protein interactions by minimizing replicability deviations but also optimizes the discovery process by enabling the simultaneous evaluation of multiple substrates. Furthermore, our pipeline accommodates distinct testing scenarios, such as blind docking or site-specific targeting, which are invaluable in applications ranging from drug repositioning to the exploration of new allosteric binding sites and toxicity assessments. BioProtIS has been designed with modularity at its core. This inherent flexibility empowers users to make custom modifications directly within the source code, tailoring the pipeline to their specific research needs. Moreover, it lays the foundation for seamless integration of diverse docking algorithms in future iterations, promising ongoing advancements in the field of computational biology. This pipeline is available for free distribution and can be download at: https://github.com/BBMDO/BioProtIS.

MicroRNA PC-5p-3991_515 mediates triflumezopyrim susceptibility in the small brown planthopper through regulating the post-transcriptional expression of P450 CYP417A2

BACKGROUND

Cytochrome P450 monooxygenases (P450s) are recognized as a major contributor to metabolic resistance in insects to most insecticides, through gene overexpressions and protein mutations. MicroRNA (miRNA), an important post-transcriptional regulator, has been reported to promote insecticide resistance by mediating the expression of detoxification enzyme genes.

RESULTS

In the present study, we reported that a novel microRNA PC-5p-3991_515 was involved in the post-transcriptional regulation of CYP417A2 and mediated the triflumezopyrim susceptibility in the small brown planthopper (SBPH), Laodelphax striatellus (Fallén). The tissue expression profiles showed that CYP417A2 was highly expressed in fat body. CYP417A2 was significantly up-regulated at 12, 36, 60, 84 and 108 h after the triflumezopyrim treatment. RNA interference (RNAi) against CYP417A2 significantly increased triflumezopyrim susceptibility in SBPH. According to the prediction by miRanda and TargetScan software, three miRNAs were indicated to bind to CYP417A2. However, when oversupply of agomir, only two miRNAs, PC-3p-625_4405 and PC-5p-3991_515, significantly increased the susceptibility to triflumezopyrim and decreased CYP417A2 levels. Furthermore, PC-5p-3991_515 was confirmed to be involved in the post-transcriptional regulation of CYP417A2 by dual luciferase reporter assay. Meanwhile, PC-5p-3991_515 was co-localized with CYP417A2 in the midgut in situ hybridization.

CONCLUSION

Our findings revealed that the novel microRNA, PC-5p-3991_515, post-transcriptionally regulated CYP417A2 expression, which then mediated the triflumezopyrim susceptibility in SBPH. © 2023 Society of Chemical Industry.

The Survival and Physiological Response of Calliptamus abbreviatus Ikonn (Orthoptera: Acrididae) to Flavonoids Rutin and Quercetin

Insect-resistant substances from plants are important natural resources that human beings can potentially develop and use to control pests. In this study, we explored the adverse effects of rutin and quercetin on grasshopper (Calliptamus abbreviatus), as well as the insect's physiological response to these substances in laboratory and field experiments. These two plant compounds exhibited toxic effects on C. abbreviatus, with quercetin showing a stronger toxicity, indicated by a lower survival, slower development, and higher induced gene expression and activities of UDP-glucuronosyltransferase, cytochrome P450s, superoxide dismutase, peroxidase and catalase, compared to rutin. These compounds, especially quercetin, have the potential to be developed as biopesticides to control grasshoppers.

Comprehensive analysis of the overexpressed cytochrome P450-based insecticide resistance mechanism in Spodoptera litura

Cytochrome P450s play critical roles in the metabolic resistance of insecticides in insects. Previous findings showed that enhanced P450 activity was an important mechanism mediating indoxacarb resistance, and multiple P450 genes were upregulated in indoxacarb resistant strains of Spodoptera litura. However, the functions of these P450 genes in insecticide resistance remain unknown. Here, the P450 inhibitor PBO effectively decreased the resistance of S. litura to indoxacarb. Ten upregulated P450 genes were characterized, all of which were overexpressed in response to indoxacarb induction. Knockdown of nine P450 genes decreased cell viability against indoxacarb, and further silencing of three genes (CYP339A1, CYP340G2, CYP321A19) in larvae enhanced the sensitivity to indoxacarb. Transgenic overexpression of these three genes increased resistance to indoxacarb in Drosophila melanogaster. Moreover, molecular modeling and docking predicted that these three P450 proteins could bind tightly to indoxacarb and N-decarbomethoxylated metabolite (DCJW). Interestingly, these three P450 genes may also mediate cross-resistance to chlorantraniliprole, λ-cyhalothrin and imidacloprid. Additionally, heterologous expression and metabolic assays confirmed that three recombinant P450s could effectively metabolize indoxacarb and DCJW. This study strongly demonstrates that multiple overexpressed mitochondrial and microsomal P450 genes were involved in insecticide resistance in S. litura.

The role of P450 enzymes in malaria and other vector-borne infectious diseases

Vector-borne infectious diseases are still an important global health problem. Malaria is the most important among them, mainly pediatric, life-threatening disease. Malaria and other vector-borne disorders caused by parasites, bacteria, and viruses have a strong impact on public health and significant economic costs. Most vector-borne diseases could be prevented by vector control, with attention to the ecological and biodiversity conservation aspects. Chemical control with pesticides and insecticides is widely used as a measure of prevention although increasing resistance to insecticides is a serious issue in vector control. Metabolic resistance is the most common mechanism and poses a big challenge. Insect enzyme systems, including monooxygenase CYP P450 enzymes, are employed by vectors mainly to metabolize insecticides thus causing resistance. The discovery and application of natural specific inhibitors/blockers of vector P450 enzymes as synergists for commonly used pesticides will contribute to the "greening" of insecticides. Besides vector CYPs, host CYP enzymes could also be exploited to fight against vector-borne diseases: using mostly their detoxifying properties and involvement in the immune response. Here, we review published research data on P450 enzymes from all players in vector-borne infections, that is, pathogens, vectors, and hosts, regarding the potential role of CYPs in disease. We discuss strategies on how to exploit cytochromes P450 in vector-borne disease control.

CYP4CL2 Confers Metabolic Resistance to Pyridaben in the Citrus Pest Mite Panonychus citri

The citrus red mite Panonychus citri has developed strong resistance to acaricides. Cytochrome P450 monooxygenases (P450s) can detoxify pesticides and are involved in pesticide resistance in many insects. Here, a pyridaben-resistant P. citri strain showed cross-resistance to cyenopyrafen, bifenazate, fenpyroximate, and tolfenpyrad. Piperonyl butoxide, a P450 inhibitor, significantly increased the toxicity of pyridaben to resistant (Pyr_Rs) and susceptible (Pyr_Control) P. citri strains. P450 activity was significantly higher in Pyr_Rs than in Pyr_Control. Analyses of RNA-Seq data identified a P450 gene (CYP4CL2) that is potentially involved in pyridaben resistance. Consistently, it was up-regulated in two field-derived resistant populations (CQ_WZ and CQ_TN). RNA interference-mediated knockdown of CYP4CL2 significantly decreased the pyridaben resistance in P. citri. Transgenic Drosophila melanogaster expressing CYP4CL2 showed increased pyridaben resistance. Molecular docking analysis showed that pyridaben could bind to several amino acids at substrate recognition sites in CYP4CL2. These findings shed light on P450-mediated pyridaben resistance in pest mites.

Evolution and functional role prediction of the CYP6DE and CYP6DJ subfamilies in Dendroctonus (Curculionidae: Scolytinae) bark beetles

Dendroctonus-bark beetles are natural components and key ecological agents of coniferous forests. They spend most of their lives under the bark, where they are exposed to highly toxic terpenes present in the oleoresin. Cytochrome P450 (CYP) is a multigene family involved in the detoxification of these compounds. It has been demonstrated that CYP6DE and CYP6DJ subfamilies hydroxylate monoterpenes, whose derivatives can act as pheromone synergist compounds or be pheromones themselves in these insects. Given the diversity and functional role of CYPs, we investigated whether these cytochromes have retained their function throughout the evolution of these insects. To test this hypothesis, we performed a Bayesian phylogenetic analysis to determine phylogenetic subgroups of cytochromes in these subfamilies. Subgroups were mapped and reconciled with the Dendroctonus phylogeny. Molecular docking analyses were performed with the cytochromes of each subgroup and enantiomers of α-pinene and β-pinene, (+)-3-carene, β-myrcene and R-(+)-limonene. In addition, functional divergence analysis was performed to identify critical amino acid sites that influence changes in catalytic site conformation and/or protein folding. Three and two phylogenetic subgroups were recovered for the CYP6DE and CYP6DJ subfamilies, respectively. Mapping and reconciliation analysis showed different gain and loss patterns for cytochromes of each subgroup. Functional predictions indicated that the cytochromes analyzed are able to hydroxylate all monoterpenes; however, they showed preferential affinities to different monoterpenes. Functional divergence analyses indicated that the CYP6DE subfamily has experimented type I and II divergence, whereas the CYP6DJ subfamily has evolved under strong functional constraints. Results suggest cytochromes of the CYP6DE subfamily evolve to reinforce their detoxifying capacity hydroxylating mainly α- and β-pinene to (+) and (-)-trans-verbenol, being the negative enantiomer used as a pheromone by several Dendroctonus species; whereas cytochromes of the CYP6DJ subfamily appear to retain their original function related to the detoxification of these compounds.

Sublethal effects of acetamiprid and afidopyropen on Harmonia axyridis: insights from transcriptomics analysis

Evaluating the sublethal effects of insecticide is crucial for protecting and utilizing natural enemies. In this study, we determined the sublethal effects of acetamiprid and afidopyropen on Harmonia axyridis (Pallas) and explored the potential molecular mechanisms underlying these effects through transcriptomics analysis. The results showed that sublethal concentrations of acetamiprid significantly reduced the adult fecundity and longevity of F₀H. axyridis and decreased the survival time and survival rate of the F₁ generation. Sublethal concentrations of afidopyropen prolonged the developmental time of 4th instar larvae in the F₀ generation. Additionally, acetamiprid and afidopyropen treatments significantly decreased the predation of H. axyridis. Furthermore, transcriptome sequencing analysis revealed that several P450 and UGT genes expressed differently when H. axyridis were exposed to sublethal concentrations of acetamiprid and afidopyropen, suggesting that the differential expression of detoxifying genes might be involved in the response and detoxification metabolism of acetamiprid and afidopyropen in H. axyridis. Our findings demonstrate that sublethal concentrations of acetamiprid adversely influences the development and predation of H. axyridis, while afidopyropen has limited effects on H. axyridis. These results are helpful for protecting and utilizing natural enemies and guiding the scientific use of pesticides in the field.

Transcription factor CncC potentially regulates cytochrome P450 CYP321A1-mediated flavone tolerance in Helicoverpa armigera

Insect P450s play crucial roles in metabolizing insecticides and toxic plant allelochemicals. In this study, our results demonstrate that Helicoverpa armigera can adapt to a lower concentration of flavone (a flavonoid phytochemical), and P450 activities and CYP321A1 transcript levels significantly increase after exposure to flavone. RNAi-mediated knockdown of CYP321A1 significantly reduced the tolerance of H. armigera larvae to flavone. In addition, the regulatory mechanisms driving CYP321A1 induction following exposure to flavone were investigated. Flavone exposure significantly increased H₂O₂ generation in the larval midgut. The mRNA levels of HaCncC and HaMaf-s significantly increased in the midgut of H. armigera after exposure to flavone. Knockdown of HaCncC significantly inhibited expression of flavone-induced CYP321A1 and resulted in a decrease in flavone induction of CYP321A1. HaCncC knockdown significantly reduced the tolerance of H. armigera larvae to flavone. Taken together, these results indicate that HaCncC regulates expression of the CYP321A1 gene responsible for flavone tolerance in H. armigera.

Potential key genes involved in metabolic resistance to malathion in the southern house mosquito, Culex quinquefasciatus, and functional validation of CYP325BC1 and CYP9M12 as candidate genes using RNA interference

BACKGROUND

Metabolic detoxification is one of the major mechanisms contributing to the development of resistance in mosquitoes, including the southern house mosquito, Culex quinquefasciatus. The three major detoxification supergene families, cytochrome P450s, glutathione S-transferases and general esterases, have been demonstrated to play an important role in metabolic resistance. In this study, we performed differential gene expression analysis based on high-throughput transcriptome sequencing on samples from four experimental groups to give insight into key genes involved in metabolic resistance to malathion in Cx. quinquefasciatus. We conducted a whole transcriptome analysis of field captured wild Cx. quinquefasciatus from Harris County (WI), Texas and a malathion susceptible laboratory-maintained Sebring colony (CO) to investigate metabolic insecticide resistance. Field captured mosquitoes were also phenotypically classified into the malathion resistant and malathion susceptible groups following a mortality response measure conducted using a Centers for Disease Control and Prevention (CDC) bottle assay. The live (MR) and dead (MS) specimens from the bottle assay, along with an unselected WI sample and a CO sample were processed for total RNA extraction and subjected to whole-transcriptome sequencing.

RESULTS

We demonstrated that the genes coding for detoxification enzymes, particularly cytochrome P450s, were highly up-regulated in the MR group compared to the MS group with similar up-regulation observed in the WI group compared to the CO group. A total of 1,438 genes were differentially expressed in comparison between MR and MS group, including 614 up-regulated genes and 824 down-regulated genes. Additionally, 1,871 genes were differentially expressed in comparison between WI and CO group, including 1,083 up-regulated genes and 788 down-regulated genes. Further analysis on differentially expressed genes from three major detoxification supergene families in both comparisons resulted in 16 detoxification genes as candidates potentially associated with metabolic resistance to malathion. Knockdown of CYP325BC1 and CYP9M12 using RNA interference on the laboratory-maintained Sebring strain significantly increased the mortality of Cx. quinquefasciatus after exposure to malathion.

CONCLUSION

We generated substantial transcriptomic evidence on metabolic detoxification of malathion in Cx. quinquefasciatus. We also validated the functional roles of two candidate P450 genes identified through DGE analysis. Our results are the first to demonstrate that knockdown of CYP325BC1 and CYP9M12 both significantly increased malathion susceptibility in Cx. quinquefasciatus, indicating involvement of these two genes in metabolic resistance to malathion.

Multiple-P450 Gene Co-Up-Regulation in the Development of Permethrin Resistance in the House Fly, Musca domestica

This paper reports a study conducted at the whole transcriptome level to characterize the P450 genes involved in the development of pyrethroid resistance, utilizing expression profile analyses of 86 cytochrome P450 genes in house fly strains with different levels of resistance to pyrethroids/permethrin. Interactions among the up-regulated P450 genes and possible regulatory factors in different autosomes were examined in house fly lines with different combinations of autosomes from a resistant house fly strain, ALHF. Eleven P450 genes that were significantly up-regulated, with levels > 2-fold those in the resistant ALHF house flies, were in CYP families 4 and 6 and located on autosomes 1, 3 and 5. The expression of these P450 genes was regulated by trans- and/or cis-acting factors, especially on autosomes 1 and 2. An in vivo functional study indicated that the up-regulated P450 genes also conferred permethrin resistance in Drosophila melanogaster transgenic lines. An in vitro functional study confirmed that the up-regulated P450 genes are able to metabolize not only cis- and trans-permethrin, but also two metabolites of permethrin, PBalc and PBald. In silico homology modeling and the molecular docking methodology further support the metabolic capacity of these P450s for permethrin and substrates. Taken together, the findings of this study highlight the important function of multi-up-regulated P450 genes in the development of insecticide resistance in house flies.

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 ² = 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.

Structural comparison of the cytochrome P450 enzymes CYP106A1 and CYP106A2 provides insight into their differences in steroid conversion

Understanding the structural basis of the selectivity of steroid hydroxylation requires detailed structural and functional investigations on various steroid hydroxylases with different selectivities, such as the bacterial cytochrome P450 enzymes. Here, the crystal structure of the cytochrome P450 CYP106A1 from Priestia megaterium was solved. CYP106A1 exhibits a rare additional structural motif of a cytochrome P450, a sixth β-sheet. The protein was found in different unusual conformations corresponding to both open and closed forms even when crystallized without any known substrate. The structural comparison of CYP106A1 with the previously investigated CYP106A2, including docking studies for both isoforms with the substrate cortisol, reveals a completely different orientation of the steroid molecule in the active sites. This distinction convincingly explains the experimentally observed differences in substrate conversion and product formation by the two enzymes.

Effects of phytocompound Precocene 1 on the expression and functionality of the P450 gene in λ-cyhalothrin-resistant Spodoptera litura (Fab.)

Spodoptera litura (Fabricius) is an agriculturally significant polyphagous insect pest that has evolved a high level of resistance to conventional insecticides. A dietary assay was used in this work to assess the resilience of field populations of S. litura to λ-cyhalothrin. Analysis of the function and expression of the cytochrome P450 gene was used to test the sensitivity of S. litura larvae to sub-lethal concentrations of the insecticidal plant chemical Precocene 1, both by itself and in combination with λ-cyhalothrin. The activity of esterase enzymes (α and β) was found to decrease 48 h post treatment with Precocene 1. The activity of GST enzyme and cytochrome P450 increased with Precocene 1 treatment post 48 h, however. Expression studies revealed the modulation by Precocene 1 of cytochrome P450 genes, CYP4M16, CYP4M15, CYP4S8V4, CYP4G31, and CYP4L10. While CYP4M16 expression was stimulated the most by the synergistic Precocene 1 + λ-cyhalothrin treatment, expression of CYP4G31 was the most down-regulated by Precocene 1 exposure. Hence, it is evident that λ-cyhalothrin-resistant pest populations are still sensitive to Precocene 1 at a sublethal concentration that is nevertheless capable of hindering their development. Precocene 1 can therefore be considered a potent candidate for the effective management of insecticide-resilient S. litura.

Knockdown of CYP6CR2 and CYP6DE5 reduces tolerance to host plant allelochemicals in the Chinese white pine beetle Dendroctonus armandi

Bark beetles rely on detoxifying enzymes to resist the defensive terpenoids of the host tree. Insect cytochrome P450 (CYPs) plays a key role in the detoxification of pesticides and plant allelochemicals. CYP6 family is unique to Insecta, and its biochemical function is basically related to the metabolism of exogenous substances. In this study, we sequenced and characterized the full-length cDNAs of two CYP6 genes from Chinese white pine beetle, Dendroctonus armandi. Spatiotemporal expression profiling revealed that the expression of CYP6CR2 and CYP6DE5 was higher in larval and adult stages of D. armandi than that in other developmental stages, and that two genes predominantly expressed in brain, midgut, fat body, Malpighian tubules or hemolymph. The expression of CYP6CR2 and CYP6DE5 was significantly induced after feeding on the phloem of Pinus armandii and exposure to six stimuli [(±)- α -pinene, (-)-α-pinene, (-)-β-pinene, (+)-3-carene, (±)-limonene and turpentine]. Importantly, silencing CYP6CR2 and CYP6DE5 separately could increase the sensitivity, led to a significant reduction of the activity of P450, resulting a significant increase in adult mortality after treatment with terpenoids. The comprehensive results of this study showed that in the process of host selection and colonization, the functions of CYPs were mainly to hydrolyze the chemical defense of the host and degrade odor molecules. These findings may help to develop new treatments to control this important pest.

CRISPR/Cas9-Mediated Knockout Reveals the Involvement of CYP304F1 in β-Cypermethrin and Chlorpyrifos Resistance in Spodoptera litura

Functions of insect CYP2 clan P450s in insecticide resistance are relatively less reported. In Spodoptera litura, a gene from the CYP2 clan (CYP304F1) was validated to be up-regulated significantly in a pyrethroid- and organophosphate-resistant population (QJ) than a susceptible population by RNA-Seq and qRT-PCR. Spatial-temporal expression indicated the high expression of CYP304F1 in the fourth, fifth, and sixth instar larvae and the metabolism-related tissue fat body and malpighian tubules. CYP304F1 was knocked out by CRISPR/Cas9, and a homozygous population (QJ-CYP304F1) with a G-base deletion at exon 2 was obtained after selection. Bioassay results showed that the LD₅₀ values to β-cypermethrin and chlorpyrifos in the QJ-CYP304F1 population decreased significantly, and the resistance ratio was both 1.81-fold in the QJ population compared with that in the QJ-CYP304F1 population. The toxicity of fenvalerate, cyhalothrin, or phoxim showed no significant change. These results suggested that CYP304F1 is involved in β-cypermethrin and chlorpyrifos resistance in S. litura.

Revealing the difference of α-amylase and CYP6AE76 gene between polyphagous Conogethes punctiferalis and oligophagous C. pinicolalis by multiple-omics and molecular biological technique

Background

Conogethes pinicolalis has been thought as a Pinaceae-feeding variant of the yellow peach moth, Conogethes punctiferalis. The divergence of C. pinicolalis from the fruit-feeding moth C. punctiferalis has been reported in terms of morphology, ecology, and genetics, however there is a lack of detailed molecular data. Therefore, in this study, we investigated the divergence of C. pinicolalis from C. punctiferalis from the aspects of transcriptomics, proteomics, metabolomics and bioinformatics.

Results

The expression of 74,611 mRNA in transcriptome, 142 proteins in proteome and 218 metabolites in metabolome presented significantly differences between the two species, while the KEGG results showed the data were mainly closely related to metabolism and redox. Moreover, based on integrating system-omics data, we found that the α-amylase and CYP6AE76 genes were mutated between the two species. Mutations in the α-amylase and CYP6AE76 genes may influence the efficiency of enzyme preference for a certain substrate, resulting in differences in metabolic or detoxifying ability in both species. The qPCR and enzyme activity test also confirmed the relevant gene expression.

Conclusions

These findings of two related species and integrated networks provide beneficial information for further exploring the divergence in specific genes, metabolism, and redox mechanism. Most importantly, it will give novel insight on species adaptation to various diets, such as from monophagous to polyphagous.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08753-9.

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.

Overexpression of Multiple Cytochrome P450 Genes Conferring Clothianidin Resistance in Bradysia odoriphaga

Cytochrome P450 monooxygenases (P450s) play important roles in the detoxification metabolism of xenobiotics and are involved in the resistance of insects to many insecticides. In this study, piperonyl butoxide (PBO), an inhibitor of P450 enzyme activity, significantly increased the toxicity of clothianidin in the clothianidin-resistant (CL-R) population of Bradysia odoriphaga. The enzyme activity of P450 in the CL-R population was significantly higher than that in the SS population. Furthermore, four P450 genes were found to be significantly overexpressed in the CL-R population. Tissue-specific expression analysis indicates that CYP9J57, CYP3828A1, CYP6SX1, and CYP6QE1 were most highly expressed in the midgut and/or Malpighian tubules. After exposure to LC₃₀ of clothianidin, the expression levels of the four P450 genes were significantly upregulated. The RNAi-mediated knockdown of CYP9J57, CYP3828A1, and CYP6QE1 significantly increased the susceptibility of B. odoriphaga to clothianidin. These results suggest that P450 genes are involved in clothianidin resistance in B. odoriphaga. This provides a better understanding of P450-mediated clothianidin resistance in B. odoriphaga and will contribute to the management of insect resistance to insecticides.

A chromosome-level genome assembly of the pollinating fig wasp Valisia javana

Fig wasp has always been thought the species-specific pollinator for their host fig (Moraceae, Ficus) and constitute a model system with its host to study co-evolution and co-speciation. The availability of a high-quality genome will help to further reveal the mechanisms underlying these characteristics. Here, we present a high-quality chromosome-level genome for Valisa javana developed by a combination of PacBio long-read and Illumina short-read. The assembled genome size is 296.34 Mb from 13 contigs with a contig N50 length of 26.76 kb. Comparative genomic analysis revealed expanded and positively selected genes related to biological features that aid fig wasps living in syconium of its highly specific host. Protein-coding genes associated with chemosensory, detoxification and venom genes were identified. Several differentially expressed genes in transcriptome data of V. javana between odor-stimulated samples and the controls have been identified in some olfactory signal transduction pathways, e.g. olfactory transduction, cAMP, cGMP-PKG, Calcim, Ras and Rap1. This study provides a valuable genomic resource for a fig wasp, and sheds insight into further revealing the mechanisms underlying their adaptive traits to their hosts in different places and co-speciation with their host.

Occurrence, Density, and Transcriptomic Response of the Leafhopper Erythroneura sudra (Hemiptera: Cicadellidae) When Confronted With Different Fruit Tree Species

The leafhopper, Erythroneura sudra (Distant) is becoming a dominant insect pest, and usually can cause significant damage to fruit production in northern China. We studied the occurrence and density of E. sudra on three fruit tree species and its transcriptomic responses when it was fed on leaves of these tree species. A higher density and survival rate of E. sudra were recorded when it fed on leaves of peach (Amygdalus persica L.) (Rosales: Rosaceae) and cherry (Cerasus pseudocerasus Lindl) (Rosales: Rosaceae) than on apple (Malus domestica Mill) (Rosales: Rosaceae). Also, feeding on M. domestica induced the largest variation in transcriptomic profiles in E. sudra. In total, 166 genes were differentially expressed (89 upregulated and 77 downregulated) in E. sudra when it fed on M. domestica, compared to when it fed on the other two tree species. The upregulated genes were mainly related to 'response to oxidative stress', 'stress-resistance', and 'xenobiotic metabolic process'. The downregulated genes were mainly related to 'structural constituent of cuticle', 'biosynthetic process', and 'development regulation'. These results suggested that M. domestica significantly changed the expression of many genes and consequently caused lower occurrence and density of E. sudra. Such information could enhance our understanding of the leafhopper-host plant relationship. Additionally, it can contribute to the improvement of current control strategies for this pest.

Contribution of multiple overexpressed carboxylesterase genes to indoxacarb resistance in Spodoptera litura

BACKGROUND

As an important family of detoxification enzymes, carboxylesterases (CarEs) have important roles in the development of insecticide resistance in almost all agricultural pests. Previous studies have suggested that enhancement of CarE activity is an important mechanism mediating indoxacarb resistance in Spodoptera litura, and several CarE genes have been found to be overexpressed in indoxacarb-resistant strains. However, the functions of these CarE genes in indoxacarb resistance needs to be further investigated.

RESULTS

The synergist triphenyl phosphate effectively reduced the resistance of S. litura to indoxacarb, suggesting an involvement of CarEs in indoxacarb resistance. Among seven identified S. litura CarE genes (hereafter SlituCOE), six were overexpressed in two indoxacarb-resistant strains, but there were no significant differences in gene copy number. Knockdown of SlituCOE009 and SlituCOE050 enhanced indoxacarb sensitivity in both susceptible and resistant strains, whereas knockdown of SlituCOE090, SlituCOE093 and SlituCOE074 enhanced indoxacarb sensitivity in only the resistant strain. Knockdown of the sixth gene, SlituCOE073, did not have any effect. Furthermore, simultaneous knockdown of the five SlituCOE genes had a greater effect on increasing indoxacarb sensitivity than silencing them individually. By contrast, overexpression of the five SlituCOE genes individually in Drosophila melanogaster significantly decreased the toxicity of indoxacarb to transgenic fruit flies. Furthermore, modeling and docking analysis indicated that the catalytic pockets of SlituCOE009 and SlituCOE074 were ideally shaped for indoxacarb and N-decarbomethoxylated metabolite (DCJW), but the binding affinity for DCJW was stronger than for indoxacarb.

CONCLUSION

This study reveals that multiple overexpressed CarE genes are involved in indoxacarb resistance in S. litura.

Pyrethroid resistance in the New World malaria vector Anopheles albimanus is mediated by cytochrome P450 CYP6P5

Pyrethroid resistance in the malaria vector Anopheles albimanus presents an obstacle to malaria elimination in the Americas. Here, An. albimanus CYP6P5 (the most overexpressed P450 in a Peruvian population) was functionally characterized. Recombinant CYP6P5 metabolized the type II pyrethroids, deltamethrin and α-cypermethrin with comparable affinities (K_M of 3.3 μM ± 0.4 and 3.6 μM ± 0.5, respectively), but exhibited a 2.7-fold higher catalytic rate for α-cypermethrin (k_cat of 6.02 min^-1 ± 0.2) versus deltamethrin (2.68 min^-1 ± 0.09). Time-course assays revealed progressive depletion of the above pyrethroids with production of four HPLC-detectable metabolites. Low depletion was obtained with type I pyrethroid, permethrin. Transgenic expression in Drosophila melanogaster demonstrated that overexpression of CYP6P5 alone conferred type II pyrethroid resistance, with only 16% and 55.3% mortalities in flies exposed to 0.25% α-cypermethrin and 0.15% deltamethrin, respectively. Synergist bioassays using P450 inhibitor piperonylbutoxide significantly recovered susceptibility (mortality = 73.6%, p < 0.001) in synergized flies exposed to 4% piperonylbutoxide, plus 0.25% α-cypermethrin, compared to non-synergized flies (mortality = 4.9%). Moderate resistance was also observed towards 4% DDT. These findings established the preeminent role of CYP6P5 in metabolic resistance in An. albimanus, highlighting challenges associated with deployment of insecticide-based control tools in the Americas.

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.

The Role of Cytochrome P450s in Insect Toxicology and Resistance

Insect cytochrome P450 monooxygenases (P450s) perform a variety of important physiological functions, but it is their role in the detoxification of xenobiotics, such as natural and synthetic insecticides, that is the topic of this review. Recent advances in insect genomics and postgenomic functional approaches have provided an unprecedented opportunity to understand the evolution of insect P450s and their role in insect toxicology. These approaches have also been harnessed to provide new insights into the genomic alterations that lead to insecticide resistance, the mechanisms by which P450s are regulated, and the functional determinants of P450-mediated insecticide resistance. In parallel, an emerging body of work on the role of P450s in defining the sensitivity of beneficial insects to insecticides has been developed. The knowledge gained from these studies has applications for the management of P450-mediated resistance in insect pests and can be leveraged to safeguard the health of important beneficial insects.

Identification and characterization of CYPs induced in the Drosophila antenna by exposure to a plant odorant

Members of the cytochrome p450 (CYP) enzyme family are abundantly expressed in insect olfactory tissues, where they are thought to act as Odorant Degrading Enzymes (ODEs). However, their contribution to olfactory signaling in vivo is poorly understood. This is due in part to the challenge of identifying which of the dozens of antennal-expressed CYPs might inactivate a given odorant. Here, we tested a high-throughput deorphanization strategy in Drosophila to identify CYPs that are transcriptionally induced by exposure to odorants. We discovered three CYPs selectively upregulated by geranyl acetate using transcriptional profiling. Although these CYPs are broadly expressed in the antenna in non-neuronal cells, electrophysiological recordings from CYP mutants did not reveal any changes in olfactory neuron responses to this odorant. Neurons were desensitized by pre-exposing flies to the odorant, but this effect was similar in CYP mutants. Together, our data suggest that the induction of a CYP gene by an odorant does not necessarily indicate a role for that CYP in neuronal responses to that odorant. We go on to show that some CYPs have highly restricted expression patterns in the antenna, and suggest that such CYPs may be useful candidates for further studies on olfactory CYP function.

RNA-Seq analysis of the blue light-emitting Orfelia fultoni (Diptera: Keroplatidae) suggest photoecological adaptations at the molecular level

Bioluminescence in Diptera is found in the Keroplatidae family, within Arachnocampininae and Keroplatinae subfamilies, with reported occurrences in Oceania, Eurasia, and Americas. Larvae of Orfelia fultoni, which inhabit stream banks in the Appalachian Mountains, emit the bluest bioluminescence among insects, using it for prey attraction, similarly to Arachnocampa spp. Although bioluminescence has a similar prey attraction function, the systems of Arachonocampininae and Keroplatinae subfamilies are morphologically/biochemically distinct, indicating different evolutionary origins. To identify the possible coding genes associated with physiological control, ecological adaptations, and origin/evolution of bioluminescence in the Keroplatinae subfamily, we performed the RNA-Seq analysis of O. fultoni larvae during day and night and compared it with the transcriptomes of Arachnocampa luminosa, and reanalyzed the previously published proteomic data of O. fultoni against the RNA-Seq dataset. The abundance of chaperones/heat-shock and hexamerin gene products at night and in luciferase enriched fractions supports their possible association and participation in bioluminescence. The low diversity of copies/families of opsins indicate a simpler visual system in O. fultoni. Noteworthy, gene products associated with silk protein biosynthesis in Orfelia were more similar to Lepidoptera than to the Arachnocampa, indicating that, similarly to the bioluminescent systems, at some point, the biochemical apparatus for web construction may have evolved independently in Orfelia and Arachnocampa.

Genome mapping coupled with CRISPR gene editing reveals a P450 gene confers avermectin resistance in the beet armyworm

The evolution of insecticide resistance represents a global constraint to agricultural production. Because of the extreme genetic diversity found in insects and the large numbers of genes involved in insecticide detoxification, better tools are needed to quickly identify and validate the involvement of putative resistance genes for improved monitoring, management, and countering of field-evolved insecticide resistance. The avermectins, emamectin benzoate (EB) and abamectin are relatively new pesticides with reduced environmental risk that target a wide number of insect pests, including the beet armyworm, Spodoptera exigua, an important global pest of many crops. Unfortunately, field resistance to avermectins recently evolved in the beet armyworm, threatening the sustainable use of this class of insecticides. Here, we report a high-quality chromosome-level assembly of the beet armyworm genome and use bulked segregant analysis (BSA) to identify the locus of avermectin resistance, which mapped on 15-16 Mbp of chromosome 17. Knockout of the CYP9A186 gene that maps within this region by CRISPR/Cas9 gene editing fully restored EB susceptibility, implicating this gene in avermectin resistance. Heterologous expression and in vitro functional assays further confirm that a natural substitution (F116V) found in the substrate recognition site 1 (SRS1) of the CYP9A186 protein results in enhanced metabolism of EB and abamectin. Hence, the combined approach of coupling gene editing with BSA allows for the rapid identification of metabolic resistance genes responsible for insecticide resistance, which is critical for effective monitoring and adaptive management of insecticide resistance.

Genetic variation at the Cyp6m2 putative insecticide resistance locus in Anopheles gambiae and Anopheles coluzzii

Background

The emergence of insecticide resistance is a major threat to malaria control programmes in Africa, with many different factors contributing to insecticide resistance in its vectors, Anopheles mosquitoes. CYP6M2 has previously been recognized as an important candidate in cytochrome P450-mediated detoxification in Anopheles. As it has been implicated in resistance against pyrethroids, organochlorines and carbamates, its broad metabolic activity makes it a potential agent in insecticide cross-resistance. Currently, allelic variation within the Cyp6m2 gene remains unknown.

Methods

Here, Illumina whole-genome sequence data from Phase 2 of the Anopheles gambiae 1000 Genomes Project (Ag1000G) was used to examine genetic variation in the Cyp6m2 gene across 16 populations in 13 countries comprising Anopheles gambiae and Anopheles coluzzii mosquitoes. To identify whether these alleles show evidence of selection either through potentially modified enzymatic function or by being linked to variants that change the transcriptional profile of the gene, hierarchical clustering of haplotypes, linkage disequilibrium, median joining networks and extended haplotype homozygosity analyses were performed.

Results

Fifteen missense biallelic substitutions at high frequency (defined as > 5% frequency in one or more populations) are found, which fall into five distinct haplotype groups that carry the main high frequency variants: A13T, D65A, E328Q, Y347F, I359V and A468S. Despite consistent reports of Cyp6m2 upregulation and metabolic activity in insecticide resistant Anophelines, no evidence of directional selection is found occurring on these variants or on the haplotype clusters in which they are found.

Conclusion

These results imply that emerging resistance associated with Cyp6m2 is potentially driven by distant regulatory loci such as transcriptional factors rather than by its missense variants, or that other genes are playing a more significant role in conferring metabolic resistance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-03757-4.

A high-quality, chromosome-level genome assembly of the Black Soldier Fly (Hermetia illucens L.)

Hermetia illucens L. (Diptera: Stratiomyidae), the Black Soldier Fly (BSF) is an increasingly important species for bioconversion of organic material into animal feed. We generated a high-quality chromosome-scale genome assembly of the BSF using Pacific Bioscience, 10X Genomics linked read and high-throughput chromosome conformation capture sequencing technology. Scaffolding the final assembly with Hi-C data produced a highly contiguous 1.01 Gb genome with 99.75% of scaffolds assembled into pseudochromosomes representing seven chromosomes with 16.01 Mb contig and 180.46 Mb scaffold N50 values. The highly complete genome obtained a Benchmarking Universal Single-Copy Orthologs (BUSCO) completeness of 98.6%. We masked 67.32% of the genome as repetitive sequences and annotated a total of 16,478 protein-coding genes using the BRAKER2 pipeline. We analyzed an established lab population to investigate the genomic variation and architecture of the BSF revealing six autosomes and an X chromosome. Additionally, we estimated the inbreeding coefficient (1.9%) of the lab population by assessing runs of homozygosity. This provided evidence for inbreeding events including long runs of homozygosity on chromosome 5. The release of this novel chromosome-scale BSF genome assembly will provide an improved resource for further genomic studies, functional characterization of genes of interest and genetic modification of this economically important species.

Multiple cytochrome P450 genes: conferring high levels of permethrin resistance in mosquitoes, Culex quinquefasciatus

Insecticides, especially pyrethroids, are the most important in the insect pest control and preventing insect vector-borne human diseases. However, insect pests, including mosquitoes, have developed resistance in the insecticides that used against them. Cytochrome P450s are associated with insecticide resistance through overexpression and detoxification mechanisms in insect species. In this study, we utilized a powerful tool, the RNAi technique, to determine the roles of key P450 genes overexpressed in permethrin resistant mosquitoes that confer insecticide resistance to unravel the molecular basis of resistance mechanisms in the mosquito Culex quinquefasciatus. The results showed that knockdown of 8 key P450 genes using RNAi techniques significantly decreased resistance to permethrin in resistant mosquitoes. In silico modeling and docking analysis further revealed the potential metabolic function of overexpressed P450 genes in the development of insecticide resistance in mosquitoes. These findings not only highlighted the functional importance of these P450 genes in insecticide resistance, but also revealed that overexpression of multiple P450 genes was responsible for the high levels of insecticide resistance in a mosquito population of Culex quinquefasciatus.

The genome of the stable fly, Stomoxys calcitrans, reveals potential mechanisms underlying reproduction, host interactions, and novel targets for pest control

Background

The stable fly, Stomoxys calcitrans, is a major blood-feeding pest of livestock that has near worldwide distribution, causing an annual cost of over $2 billion for control and product loss in the USA alone. Control of these flies has been limited to increased sanitary management practices and insecticide application for suppressing larval stages. Few genetic and molecular resources are available to help in developing novel methods for controlling stable flies.

Results

This study examines stable fly biology by utilizing a combination of high-quality genome sequencing and RNA-Seq analyses targeting multiple developmental stages and tissues. In conjunction, 1600 genes were manually curated to characterize genetic features related to stable fly reproduction, vector host interactions, host-microbe dynamics, and putative targets for control. Most notable was characterization of genes associated with reproduction and identification of expanded gene families with functional associations to vision, chemosensation, immunity, and metabolic detoxification pathways.

Conclusions

The combined sequencing, assembly, and curation of the male stable fly genome followed by RNA-Seq and downstream analyses provide insights necessary to understand the biology of this important pest. These resources and new data will provide the groundwork for expanding the tools available to control stable fly infestations. The close relationship of Stomoxys to other blood-feeding (horn flies and Glossina) and non-blood-feeding flies (house flies, medflies, Drosophila) will facilitate understanding of the evolutionary processes associated with development of blood feeding among the Cyclorrhapha.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-00975-9.

Primary Metabolism co-Opted for Defensive Chemical Production in the Carabid Beetle, Harpalus pensylvanicus

Of the approximately one million described insect species, ground beetles (Coleoptera: Carabidae) have long captivated the attention of evolutionary biologists due to the diversity of defensive compounds they synthesize. Produced using defensive glands in the abdomen, ground beetle chemicals represent over 250 compounds including predator-deterring formic acid, which has evolved as a defensive strategy at least three times across Insecta. Despite being a widespread method of defense, formic acid biosynthesis is poorly understood in insects. Previous studies have suggested that the folate cycle of one-carbon (C1) metabolism, a pathway involved in nucleotide biosynthesis, may play a key role in defensive-grade formic acid production in ants. Here, we report on the defensive gland transcriptome of the formic acid-producing ground beetle Harpalus pensylvanicus. The full suite of genes involved in the folate cycle of C1 metabolism are significantly differentially expressed in the defensive glands of H. pensylvanicus when compared to gene expression profiles in the rest of the body. We also find support for two additional pathways potentially involved in the biosynthesis of defensive-grade formic acid, the kynurenine pathway and the methionine salvage cycle. Additionally, we have found an array of differentially expressed genes in the secretory lobes involved in the biosynthesis and transport of cofactors necessary for formic acid biosynthesis, as well as genes presumably involved in the detoxification of secondary metabolites including formic acid. We also provide insight into the evolution of the predominant gene family involved in the folate cycle (MTHFD) and suggest that high expression of folate cycle genes rather than gene duplication and/or neofunctionalization may be more important for defensive-grade formic acid biosynthesis in H. pensylvanicus. This provides the first evidence in Coleoptera and one of a few examples in Insecta of a primary metabolic process being co-opted for defensive chemical biosynthesis. Our results shed light on potential mechanisms of formic acid biosynthesis in the defensive glands of a ground beetle and provide a foundation for further studies into the evolution of formic acid-based chemical defense strategies in insects.

Differential regulation of cytochrome P450 genes associated with biosynthesis and detoxification in bifenthrin-resistant populations of navel orangewom (Amyelois transitella)

Pyrethroid resistance was first reported in 2013 for the navel orangeworm, Amyelois transitella, but the genetic underpinnings of pyrethroid resistance are unknown. We investigated the role of cytochrome P450 monooxygenases (P450s) belonging to the CYP3 and CYP4 clans using colonies derived from individuals collected in 2016 from almond orchards in two counties. One colony (ALM) originated from an almond orchard in Madera County with no reported pyrethroid resistance and the second colony (R347) originated from the same Kern County orchard where pyrethroid resistance was first reported. We used high-throughput quantitative real-time PCR (qRT-PCR) analyses of 65 P450s in the CYP3 and CYP4 clans of A. transitella to identify P450s induced by bifenthrin and associated with pyrethroid resistance. Nine P450s were constitutively overexpressed in R347 compared to ALM, including CYP6AE54 (11.7-fold), belonging to a subfamily associated with metabolic pesticide detoxification in Lepidoptera and CYP4G89 (33-fold) belonging to a subfamily associated with cuticular hydrocarbon (CHC) synthesis and resistance via reduced pesticide penetrance. Cuticular hydrocarbons analysis revealed that R347 produced twice as many total CHCs in the egg and adult stages as ALM. Topical toxicity bioassays for R347 determined that egg mortality was reduced at low bifenthrin concentrations and larval mortality was reduced at high concentrations of bifenthrin compared to ALM. Our discovery of both changes in metabolism and production of CHCs for R347 have implications for the possible decreased efficacy of other classes of insecticide used to control this insect. The threat of widespread pyrethroid resistance combined with the potential for cross-resistance to develop through the mechanism of reduced penetrance warrants developing management strategies that facilitate insecticide passage across the cuticle.

Disruption of the cytochrome P450 CYP6BQ7 gene reduces tolerance to plant toxicants in the red flour beetle, Tribolium castaneum

In insects, the cytochrome P450 CYP6B family plays key roles in the detoxification of toxic plant substances. However, the function of CYP6 family genes in degrading plant toxicants in Tribolium castaneum, an extremely destructive global storage pest, have yet to be elucidated. In this study, a T. castaneum CYP gene, TcCYP6BQ7, was characterized. TcCYP6BQ7 expression was significantly induced after exposure to essential oil of the plant Artemisia vulgaris (EOAV). Spatiotemporal expression profiling revealed that TcCYP6BQ7 expression was higher in larval and adult stages of T. castaneum than in other developmental stages, and that TcCYP6BQ7 was predominantly expressed in the brain and hemolymph from the late larval stage. TcCYP6BQ7 silencing by RNA interference increased larvae mortality in response to EOAV from 49.67% to 71.67%, suggesting that this gene is associated with plant toxicant detoxification. Combined results from this study indicate that the CYP6 family gene TcCYP6BQ7 likely plays a pivotal role in influencing the susceptibility of T. castaneum to plant toxicants. These findings may have implications for the development of novel therapeutics to control this agriculturally important pest.

Roles of the variable P450 substrate recognition sites SRS1 and SRS6 in esfenvalerate metabolism by CYP6AE subfamily enzymes in Helicoverpa armigera

The cotton bollworm P450s of the clustered CYP6AE subfamily share high sequence identities but differ dramatically in their capacity to metabolize xenobiotics, especially esfenvalerate. Among them, CYP6AE17 has the highest sequence identity with CYP6AE18 but shows ~7-fold higher metabolic efficiency. CYP6AE11 is most active towards esfenvalerate but CYP6AE20 is inactive even though the enzymes share 54.8% sequence identity. Sequence analysis revealed the SRS1 (Substrate Recognition Site) and SRS6 between CYP6AE17 and CYP6AE18, and SRS1 between CYP6AE11 and CYP6AE20 are the most variable among all six SRSs. In order to identify the key factors that underlie the observed catalytic difference, we exchanged these SRS sequences between two pairs of P450s and studied the activity of the resulting hybrid mutants or chimeras. In vitro metabolism showed that the CYP6AE17/18 chimeras had 2- and 14-fold decreased activities and the CYP6AE18/17 chimeras had 6- and 10-fold increased activities to esfenvalerate. Meanwhile, after exchanging SRS1 with each other, the CYP6AE11/20 chimera folded incorrectly but the CYP6AE20/11 chimera gained moderate activity to esfenvalerate. Molecular modelling showed that amino acids variants within SRS1 or SRS6 change the shape and chemical environment of the active sites, which may affect the ligand-binding interactions. These results indicate that the protein structure variation resulting from the sequence diversity of SRSs promotes the evolution of insect chemical defense and contributes to the development of insect resistance to pesticides.

(+)-Sesamin-oxidising CYP92B14 shapes specialised lignan metabolism in sesame

Sesamum spp. (sesame) are known to accumulate a variety of lignans in a lineage-specific manner. In cultivated sesame (Sesamum indicum), (+)-sesamin, (+)-sesamolin and (+)-sesaminol triglucoside are the three major lignans found richly in the seeds. A recent study demonstrated that SiCYP92B14 is a pivotal enzyme that allocates the substrate (+)-sesamin to two products, (+)-sesamolin and (+)-sesaminol, through multiple reaction schemes including oxidative rearrangement of α-oxy-substituted aryl groups (ORA). In contrast, it remains unclear whether (+)-sesamin in wild sesame undergoes oxidation reactions as in S. indicum and how, if at all, the ratio of the co-products is tailored at the molecular level. Here, we functionally characterised SrCYP92B14 as a SiCYP92B14 orthologue from a wild sesame, Sesamum radiatum, in which we revealed accumulation of the (+)-sesaminol derivatives (+)-sesangolin and its novel structural isomer (+)-7´-episesantalin. Intriguingly, SrCYP92B14 predominantly produced (+)-sesaminol either through ORA or direct oxidation on the aromatic ring, while a relatively low but detectable level of (+)-sesamolin was produced. Amino acid substitution analysis suggested that residues in the putative distal helix and the neighbouring heme propionate of CYP92B14 affect the ratios of its co-products. These data collectively show that the bimodal oxidation mechanism of (+)-sesamin might be widespread across Sesamum spp., and that CYP92B14 is likely to be a key enzyme in shaping the ratio of (+)-sesaminol- and (+)-sesamolin-derived lignans from the biochemical and evolutionary perspectives.

Post-transcriptional modulation of cytochrome P450s, Cyp6g1 and Cyp6g2, by miR-310s cluster is associated with DDT-resistant Drosophila melanogaster strain 91-R

The role of miRNAs in mediating insecticide resistance remains largely unknown, even for the model species Drosophila melanogaster. Building on prior research, this study used microinjection of synthetic miR-310s mimics into DDT-resistant 91-R flies and observed both a significant transcriptional repression of computationally-predicted endogenous target P450 detoxification genes, Cyp6g1 and Cyp6g2, and also a concomitant increase in DDT susceptibility. Additionally, co-transfection of D. melanogaster S2 cells with dual luciferase reporter constructs validated predictions that miR-310s bind to target binding sites in the 3ʹ untranslated regions (3ʹ-UTR) of both Cyp6g1 and Cyp6g2 in vitro. Findings in the current study provide empirical evidence for a link between reduced miRNA expression and an insecticidal resistance phenotype through reduced targeted post-transcriptional suppression of transcripts encoding proteins involved in xenobiotic detoxification. These insights are important for understanding the breadth of adaptive molecular changes that have contributed to the evolution of DDT resistance in D. melanogaster.

Chemosensation and Evolution of Drosophila Host Plant Selection

The ability to respond to chemosensory cues is critical for survival of most organisms. Among insects, Drosophila melanogaster has the best characterized olfactory system, and the availability of genome sequences of 30 Drosophila species provides an ideal scenario for studies on evolution of chemosensation. Gene duplications of chemoreceptor genes allow for functional diversification of the rapidly evolving chemoreceptor repertoire. Although some species of the genus Drosophila are generalists for host plant selection, rapid evolution of olfactory receptors, gustatory receptors, odorant-binding proteins, and cytochrome P450s has enabled diverse host specializations of different members of the genus. Here, I review diversification of the chemoreceptor repertoire among members of the genus Drosophila along with co-evolution of detoxification mechanisms that may have enabled occupation of diverse host plant ecological niches.

Expression and functional analysis of P450 gene induced tolerance/resistance to lambda-cyhalothrin in quercetin fed larvae of beet armyworm Spodoptera exigua (Hübner)

Beet armyworm, Spodoptera exigua (Hübner) is an agronomical important and most devastating polyphagous pest that damages a variety of crops around the globe including China. Quercetin is one of the abundant dietary flavonoids and the important defense allelochemicals in plants. Therefore, the changes in insect detoxification enzymes activities in response to plants allelochemicals may result increased the sensitivity to insecticides. In this study, we examined the induced effect of quercetin on larval tolerance to lambda-cyhalothrin in S. exigua. Application of cytochrome P450 inhibitor piperonyl butoxide (PBO) significantly synergized the lambda-cyhalothrin toxicity in quercetin-fed S. exigua larvae. Moreover, larval weight significantly reduced in quercetin, lambda-cyhalothrin, and quercetin + lambda-cyhalothrin treatment. Furthermore, our results showed that the P450 detoxification enzyme effectively increased in all treatments as compared to the control. Quantitative Real-time PCR analysis revealed that expression level of CYP6AE10 significantly upregulated in larvae treated with quercetin, lambda-cyhalothrin and quercetin + lambda-cyhalothrin in the midgut and fat body respectively. In addition, RNAi mediated knockdown of CYP6AE10 in S. exigua larvae significantly decreased the transcription level of target cytochrome P450 gene followed by the exposure with quercetin, lambda-cyhalothrin, and quercetin + lambda-cyhalothrin. Similarly, the knockdown of CYP6AE10 by the injection of dsRNA led to increased mortality after the treatment with respective chemicals. Overall, these data showed that P450s might possibly play an important role in the metabolic adaptation of S. exigua larvae to its host plant defense allelochemicals as well as insecticides. In conclusion, S. exigua can take benefit from its host plant's secondary metabolites to elaborate its defense against synthetic insecticides.

Functional analysis of four upregulated carboxylesterase genes associated with fenpropathrin resistance in Tetranychus cinnabarinus (Boisduval)

BACKGROUND

Carboxylesterases (CarEs) are important in pesticide resistance. Four overexpressed CarE genes with inducible character were screened out in fenpropathrin-resistant Tetranychus cinnabarinus, but their functional roles remained to be further analyzed by RNAi and protein expression.

RESULTS

Feeding a single double-stranded (ds)RNA of each of four genes led to gene-specific downregulation of mRNA, decreased esterase activity and diminished resistance in T. cinnabarinus. More interestingly, feeding four dsRNAs simultaneously led to a more significant decrease in enzymatic activity and fold resistance than feeding a single dsRNA individually, suggesting that these CarE genes were involved in fenpropathrin-resistance and had cooperative roles. The gene CarE6 was regarded as the primary and representative candidate to be functionally expressed, because silencing of CarE6 led to the most significant decrease in resistance level. The activity of CarE6 protein was competitively inhibited by fenpropathrin. It could effectively decompose 41.7 ± 0.09% of fenpropathrin within 3 h, proving that CarE6 protein was capable of metabolizing fenpropathrin effectively in T. cinnabarinus.

CONCLUSION

The results confirm that four CarE genes are cooperatively involved in fenpropathrin resistance and the metabolic enzymes encoded by these overexpressed genes do indeed metabolize acaricide in resistant T. cinnabarinus in the evolution of acaricide resistance. © 2018 Society of Chemical Industry.

Molecular recognition of wood polyphenols by phase II detoxification enzymes of the white rot Trametes versicolor

Wood decay fungi have complex detoxification systems that enable them to cope with secondary metabolites produced by plants. Although the number of genes encoding for glutathione transferases is especially expanded in lignolytic fungi, little is known about their target molecules. In this study, by combining biochemical, enzymatic and structural approaches, interactions between polyphenols and six glutathione transferases from the white-rot fungus Trametes versicolor have been demonstrated. Two isoforms, named TvGSTO3S and TvGSTO6S have been deeply studied at the structural level. Each isoform shows two distinct ligand-binding sites, a narrow L-site at the dimer interface and a peculiar deep hydrophobic H-site. In TvGSTO3S, the latter appears optimized for aromatic ligand binding such as hydroxybenzophenones. Affinity crystallography revealed that this H-site retains the flavonoid dihydrowogonin from a partially purified wild-cherry extract. Besides, TvGSTO6S binds two molecules of the flavonoid naringenin in the L-site. These data suggest that TvGSTO isoforms could interact with plant polyphenols released during wood degradation.

Volatile β-Ocimene Can Regulate Developmental Performance of Peach Aphid Myzus persicae Through Activation of Defense Responses in Chinese Cabbage Brassica pekinensis

In nature, plants have evolved sophisticated defense mechanisms against the attack of pathogens and insect herbivores. Plant volatile-mediated plant-to-plant communication has been assessed in multitrophic systems in different plant species and different pest species. β-ocimene is recognized as an herbivore-induced plant volatile that play an important role in the chemical communication between plants and pests. However, it is still unclear whether β-ocimene can active the defense mechanism of Chinese cabbage Brassica pekinensis against the peach aphid Myzus persicae. In this study, we found that treatment of Chinese cabbage with β-ocimene inhibited the growth of M. persicae in terms of weight gain and reproduction. Moreover, β-ocimene treatment negatively influenced the feeding behavior of M. persicae by shortening the total feeding period and phloem ingestion and increasing the frequency of stylet puncture. When given a choice, winged aphids preferred to settle on healthy Chinese cabbage compared with β-ocimene-treated plants. In addition, performance of the parasitoid Aphidius gifuensis in terms of Y-tube olfaction and landings was better on β-ocimene-treated Chinese cabbage than on healthy plants. Furthermore, β-ocimene significantly increased the expression levels of salicylic acid and jasmonic acid marker genes and the accumulation of glucosinolates. Surprisingly, the transcriptional levels of detoxifying enzymes (CYP6CY3, CYP4, and GST) in aphids reared on β-ocimene-treated Chinese cabbage were significantly higher than those of aphids reared on healthy plants. In summary, our results indicated that β-ocimene can activate the defense response of Chinese cabbage against M. persicae, and that M. persicae can also adjust its detoxifying enzymes machinery to counter the host plant defense reaction.

Prokaryotic functional expression and activity comparison of three CYP9A genes from the polyphagous pest Helicoverpa armigera

Cytochrome P450s (CYPs or P450s) have been long recognized as very important enzymes in the metabolism of xenobiotic and endogenous compounds, but only a few CYPs have been functionally characterized in insects. The effort in functional characterization of insect P450s is heavily hindered by technical difficulties in preparing active, individual P450 enzymes directly from the target insect. In this paper, we describe the functional expression of two additional pyrethroid resistance-associated CYP9A genes (CYP9A12 and CYP9A17) from the polyphagous pest Helicoverpa armigera in the facile Escherichia coli. The functionality of E. coli produced CYP9A12, CYP9A14, and CYP9A17 was investigated and activities of these CYP9As were compared against three probe substrates after reconstitution with NADPH-dependent cytochrome P450 reductase. The results showed that active forms of CYP9A12 and CYP9A17 were expressed in E. coli with a content of about 1.0-1.5 nmol mg-1 protein in membrane preparations. In vitro assays showed that CYP9A14 was capable of catalyzing O-dealkylation of methoxyresorufin (MROD), ethoxyresorufin (EROD), and benzyloxyresorufin (BROD), while CYP9A12 and CYP9A17 exhibited only MROD and EROD activities. Kinetic studies demonstrated that CYP9A14 had the greatest k cat/K m value for MROD, and CYP9A17 for EROD, while the lowest k cat/K m values for both MROD and EROD were observed for CYP9A12. The distinct biochemical traits suggest that the three paralogous CYP9As may play different roles in xenobiotic metabolism in this important pest.

Assessment of biochemical mechanisms of tolerance to chlorpyrifos in ancient and contemporary Daphnia pulicaria genotypes

The evolution of tolerance to environmental contaminants in non-target taxa has been largely studied by comparing extant populations experiencing contrasting exposure. Previous research has demonstrated that "resurrected" genotypes from a population of Daphnia pulicaria express temporal variation in sensitivity to the insecticide chlorpyrifos. Ancient genotypes (1301-1646AD.) were on average more sensitive to this chemical compared to the contemporary genotypes (1967-1977AD.). To determine the physiological mechanisms of tolerance, a series of biochemical assays was performed on three ancient and three contemporary genotypes; these six genotypes exhibited the most sensitive and most tolerant phenotypes within the population, respectively. Metabolic tolerance mechanisms were evaluated using acute toxicity testing, while target-site tolerance was assessed via in vitro acetylcholinesterase (AChE) assays. Acute toxicity tests were conducted using i) the toxic metabolite chlorpyrifos-oxon (CPF-oxon) and ii) CPF-oxon co-applied with piperonyl butoxide (PBO), a known Phase-I metabolic inhibitor. Both series of toxicity tests reduced the mean variation in sensitivity between tolerant and sensitive genotypes. Exposure to CPF-O reduced the disparity from a 4.7-fold to 1.6-fold difference in sensitivity. The addition of PBO further reduced the variation to a 1.2-fold difference in sensitivity. In vitro acetylcholinesterase assays yielded no significant differences in constitutive activity or target-site sensitivity. These findings suggest that pathways involving Phase-I detoxification and/or bioactivation of chlorpyrifos play a significant role in dictating the microevolutionary trajectories of tolerance in this population.

Transcriptional profile and differential fitness in a specialist milkweed insect across host plants varying in toxicity

Interactions between plants and herbivorous insects have been models for theories of specialization and co-evolution for over a century. Phytochemicals govern many aspects of these interactions and have fostered the evolution of adaptations by insects to tolerate or even specialize on plant defensive chemistry. While genomic approaches are providing new insights into the genes and mechanisms insect specialists employ to tolerate plant secondary metabolites, open questions remain about the evolution and conservation of insect counterdefences, how insects respond to the diversity defences mounted by their host plants, and the costs and benefits of resistance and tolerance to plant defences in natural ecological communities. Using a milkweed-specialist aphid (Aphis nerii) model, we test the effects of host plant species with increased toxicity, likely driven primarily by increased secondary metabolites, on aphid life history traits and whole-body gene expression. We show that more toxic plant species have a negative effect on aphid development and lifetime fecundity. When feeding on more toxic host plants with higher levels of secondary metabolites, aphids regulate a narrow, targeted set of genes, including those involved in canonical detoxification processes (e.g., cytochrome P450s, hydrolases, UDP-glucuronosyltransferases and ABC transporters). These results indicate that A. nerii marshal a variety of metabolic detoxification mechanisms to circumvent milkweed toxicity and facilitate host plant specialization, yet, despite these detoxification mechanisms, aphids experience reduced fitness when feeding on more toxic host plants. Disentangling how specialist insects respond to challenging host plants is a pivotal step in understanding the evolution of specialized diet breadths.

A cytochrome P450 gene plays a role in the recognition of sex pheromones in the tobacco cutworm, Spodoptera litura

Cytochrome P450 (P450 or CYP) genes are involved in fundamental physiological functions, and might be also associated with the olfactory recognition of sex pheromones in beetles and moths. A P450 gene, Spodoptera litura CYP4L4 (SlituCYP4L4), was cloned for the first time from the antennae of S. litura. SlituCYP4L4 was almost exclusively expressed in the adult stage and predominantly expressed in the adult antennae. In situ hybridization showed that SlituCYP4L4 localized mainly at the base of the long sensilla trichoidea, which responds to sex pheromone components. Pretreatment with an S. litura sex pheromone significantly reduced the expression levels of SlituCYP4L4, consistent with other genes involved in sex pheromone recognition. The expression level of SlituCYP4L4 was different in moths collected with different ratios of sex pheromone lures and collected in different geographical locations. After gene knockdown of SlituCYP4L4 in the antennae, the electroantennogram (EAG) responses of male and female moths to (9Z,11E)-tetradecadienyl acetate or (9Z,12E)-tetradecadienyl acetate were significantly decreased. In contrast, EAG responses to plant volatiles and sex pheromones of other moth species were not significantly influenced in these moths. SlituCYP4L4 was also expressed in the gustatory tissues and sensilla, which suggests that SlituCYP4L4 may have other functions in the chemosensory system. Our results have shown for the first time the function of a CYP gene with appendage-specific expression in insect sex pheromone recognition, especially in adult moths.