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

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.

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.

Role of the epsilon glutathione S-transferases in xanthotoxin tolerance in Spodoptera litura

Spodoptera litura, a polyphagous lepidopteran pest, demonstrates a remarkable capacity to adapt to varying host plants by efficiently detoxifying phytochemicals. However, the underlying mechanism for this adaptation is not well understood. Herein, twenty eplison glutathione S-transferase genes (GSTes) were characterized and their roles in phytochemical tolerance were analyzed in S. litura. Most of the GSTe genes were mainly expressed in the larval midgut and fat body. Exposure to the phytochemicals, especially xanthotoxin, induced the expression of most GSTe genes. Molecular docking analysis revealed that xanthotoxin could form stable bonds with six xanthotoxin-responsive GSTes, with binding free energies ranging from -36.44 to -68.83 kcal mol-1. Knockdown of these six GSTe genes increased the larval susceptibility to xanthotoxin. Furthermore, xanthotoxin exposure significantly upregulated the expression of two transcription factor genes CncC and MafK. Silencing of either CncC or MafK reduced the expression of GSTe16, which exhibited the largest change in response to xanthotoxin. Additionally, analysis of the promoter sequence of GSTe16 revealed the presence of seven CncC/Maf binding sites. Luciferase reporter assays showed that CncC and MafK enhanced the expression of GSTe16, leading to the increased xanthotoxin tolerance in S. litura. These findings provide insight into the functions and transcriptional regulatory mechanisms of GSTes, thereby enhancing our understanding of the role of GSTs in the adaptation of lepidopteran pests to phytochemicals.

CYP321A Subfamily P450s Contribute to the Detoxification of Phytochemicals and Pyrethroids in Spodoptera litura

Although the induction of cytochrome P450 monooxygenases involved in insect detoxification has been well documented, the underlying regulatory mechanisms remain obscure. In Spodoptera litura, CYP321A subfamily members were effectively induced by exposure to flavone, xanthotoxin, curcumin, and λ-cyhalothrin, while knockdown of the CYP321A genes increased larval susceptibility to these xenobiotics. Homology modeling and molecular docking analyses showed that these four xenobiotics could stably bind to the CYP321A enzymes. Furthermore, two transcription factor genes, CncC and MafK, were significantly induced by the xenobiotics. Knockdown of CncC or MafK reduced the expression of four CYP321A genes and increased larval susceptibility to the xenobiotics. Dual-luciferase reporter assays showed that cotransfection of reporter plasmids carrying the CYP321A promoter with CncC and/or MafK-expressing constructs significantly magnified the promoter activity. These results indicate that the induction of CYP321A subfamily members conferring larval detoxification capability to xenobiotics is mediated by the activation of CncC and MafK.

Inducible Cytochrome P450s in the Fat Body and Malpighian Tubules of the Polyphagous Pests of Spodoptera litura Confer Xenobiotic Tolerance

Cytochrome P450 plays vital roles in detoxifying xenobiotics. In this study, SlCYP340A and SlCYP340L expression in the Spodoptera litura fat body and SlCYP332A1, SlCYP6AB12, SlCYP6AB58, SlCYP6AB59, and SlCYP6AN4 expression in the Malpighian tubules were significantly upregulated after cyantraniliprole exposure, and SlCYP6AB58 and SlCYP6AB59 expression levels were simultaneously increased in the Malpighian tubules after gossypol treatment. Drosophila ectopically expressing candidate P450 genes showed that SlCYP332A1, SlCYP6AB12, SlCYP6AB59, SlCYP6AN4, and SlCYP340A conferred cyantraniliprole tolerance. The overexpression of SlCYP6AB58 and SlCYP6AB59 in Drosophila increased the number of eggs laid under the gossypol treatment. Moreover, the knockdown of SlCYP332A1, SlCYP6AB12, SlCYP6AB59, SlCYP6AN4, and SlCYP340A increased S. litura mortality under the cyantraniliprole treatment. Homology modeling and molecular docking results suggested that candidate P450 has the potential to bind with cyantraniliprole. These results indicate that the CYP3 and CYP4 genes participate in cyantraniliprole detoxification and that SlCYP6AB59 may be simultaneously involved in the gossypol tolerance of S. litura.

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.

Contribution of UDP-glycosyltransferases to chlorpyrifos resistance in Nilaparvata lugens

As a multigene superfamily of Phase II detoxification enzymes, uridine diphosphate (UDP)-glycosyltransferases (UGTs) play important roles in the metabolism of xenobiotics including insecticides. In this study, 5-nitrouracil, an inhibitor of UGT enzyme activity, effectively increased the toxicity of chlorpyrifos to the chlorpyrifos-resistant strain of Nilaparvata lugens, one of the most resistant rice pests. The enzyme content of UGT in the resistant strain was significantly higher than that in the susceptible strain. Among 20 identified UGT genes, UGT386H2, UGT386J2, UGT386N2 and UGT386P1 were found significantly overexpressed in the resistant strain and can be effectively induced by chlorpyrifos. These four UGT genes were most highly expressed in the midgut and/or fat body, two main insect detoxification tissues. Amino acid sequence alignments revealed that these four UGTs contained a variable N-terminal substrate-binding domain and a conserved C-terminal sugar donor-binding domain. Furthermore, homology modeling and molecular docking analyses showed that these UGTs could stably bind to chlorpyrifos and chlorpyrifos oxon, with the binding free energies from -19.4 to -110.62 kcal mol^-1. Knockdown of UGT386H2 or UGT386P1 by RNA interference dramatically increased the susceptibility of the resistant strain to chlorpyrifos. These findings suggest that overexpression of these two UGT genes contributes to chlorpyrifos resistance in N. lugens.

Differential regulation of the Tor gene homolog drives the red/green pigmentation phenotype in the aphid Myzuspersicae

In some aphid species, intraspecific variation in body colour is caused by differential carotenoid content: whilst green aphids contain only yellow carotenoids (β-, γ-, and β,γ-carotenes), red aphids additionally possess red carotenoids (torulene and 3,4-didehydrolycopene). Unusually, within animals who typically obtain carotenoids from their diet, ancestral horizontal gene transfer of carotenoid biosynthetic genes from fungi (followed by gene duplication), have imbued aphids with the intrinsic gene repertoire necessary to biosynthesise carotenoids. In the pea aphid, Acyrthosiphon pisum a lycopene (phytoene) desaturase gene (Tor) underpins the red/green phenotype, with this locus present in heterozygous form in red individuals but absent in green aphids, resulting in them being unable to convert lycopene into the red compounds 3,4-didehydrolycopene and torulene. The green peach aphid, Myzus persicae, separated from the pea aphid for ≈45MY also exists as distinct colour variable morphs, with both red and green individuals present. Here, we examined genomic data for both red and green morphs of M. persicae and identified an enlarged (compared to A. pisum) repertoire of 16 carotenoid biosynthetic genes (11 carotenoid desaturases and five carotenoid cyclase/synthase genes). From these, we identify the homolog of A. pisum Tor (here called carotene desaturase 2 or CDE-2) and show through 3D modelling that this homolog can accommodate the torulene precursor lycopene and, through RNA knockdown feeding experiments, demonstrate that disabling CDE-2 expression in red M. persicae clones results in green-coloured offspring. Unlike in A. pisum, we show that functional CDE-2 is present in the genomes of both red and green aphids. However, expression differences between the two colour morphs (350-700 fold CDE-2 overexpression in red clones), potentially driven by variants identified in upstream putative regulatory elements, underpin this phenotype. Thus, whilst aphids have a common origin of their carotenoid biosynthetic pathway, two aphid species separated for over 40MY have evolved very different drivers of intraspecific colour variation.

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.

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.

Pyrethrin-resembling pyrethroids are metabolized more readily than heavily modified ones by CYP9As from Helicoverpa armigera

The cotton bollworm, Helicoverpa armigera, is a polyphagous pest threatening many economically important crops worldwide. Until recently, synthetic pyrethroids remain in wide use for controlling pest insects including the cotton bollworm. Understanding the metabolic mechanism of pyrethroids in a given pest can provide significant implication for a smart choice of insecticides, and such information is useful for the development of novel selective and safe insecticides. In this study, we used complexes of recombinant H. armigera cytochrome P450 CYP9A and NADPH-dependent cytochrome P450 reductase to investigate the capacity of three CYP9A paralogs in the transformation of seven structurally different pyrethroids by metabolism assays. The results showed that the three paralogous CYP9As were able to metabolize multiple pyrethroids. Interestingly, all the three CYP9As transformed pyrethrin-resembling pyrethroids (e.g. bioallethrin) more efficiently than the heavily modified ones (e.g. bifenthrin). These findings suggest that herbivorous insects can cope with synthetic insecticides using their physiological systems that initially evolved to survive exposure to the defensive chemicals in their host plants, adding support to the pre-adaptation hypothesis.