Characterization and expression of the cytochrome P450 gene family in diamondback moth, Plutella xylostella (L.)

Genome-Wide Comparative Analysis of the Cytochrome P450 Monooxygenase Family in 19 Aphid Species and Their Expression Analysis in 4 Cereal Crop Aphids

Cytochrome P450 monooxygenases (CYP450s) play a variety of physiological roles, including pesticide resistance, plant allelochemical detoxification, and hormone metabolism catalysis. However, limited information is available on the classification and expression profiles of the CYP450 gene family in aphid species. This is the first study to identify the cytochrome P450 gene family in 19 aphid species at the whole genome level. A total of 1100 CYP450 genes were identified in 19 aphid species. Three hundred CYP450 genes belonged to six cereal crop aphid species, which were further classified into four subfamilies according to the phylogenetic relationship. The conserved motifs, exon-intron structures, and genomic organization of the same subfamilies were similar. Predictions of subcellular localization revealed that the endoplasmic reticulum harbored the majority of CYP450 proteins. In Sitobion avenae and Rhopalosiphum maidis, the increase in the CYP450 gene was primarily caused by segmental duplication events. However, only tandem duplication occurred in the CYP450 gene family of Diuraphis noxia, Rhopalosiphum padi, Schizaphis graminum, and Sitobion miscanthi. Synteny analysis found three continuous colinear CYP450 gene pairs among six cereal crop aphid species. Furthermore, we obtained the expression profiles of four cereal crop aphids, including R. padi, D. noxia, S. graminum, and S. avenae. Differential expression analysis provided growth stage specificity genes, tissue specificity genes, organ specificity genes and some detoxification metabolic genes among these four cereal crop aphids. Meanwhile, their expression patterns were showed. The related functions and pathways of CYP450s were revealed by GO and KEGG enrichment analysis. Above all, we picked the differentially expressed CYP450 genes from all of the differentially expressed genes (DEGs). These differentially expressed CYP450 genes provided some new potential candidates for aphid control and management. This work establishes the foundation for further investigations into the regulatory functions of the CYP450 gene family in aphid species and beyond.

Molecular mechanisms of cytochrome P450-mediated detoxification of tetraniliprole, spinetoram, and emamectin benzoate in the fall armyworm, Spodoptera frugiperda (J.E. Smith)

The fall armyworm (FAW) Spodoptera frugiperda (J.E. Smith) is a highly damaging invasive omnivorous pest that has developed varying degrees of resistance to commonly used insecticides. To investigate the molecular mechanisms of tolerance to tetraniliprole, spinetoram, and emamectin benzoate, the enzyme activity, synergistic effect, and RNA interference were implemented in S. frugiperda. The functions of cytochrome P450 monooxygenase (P450) in the tolerance to tetraniliprole, spinetoram, and emamectin benzoate in S. frugiperda was determined by analysing changes in detoxification metabolic enzyme activity and the effects of enzyme inhibitors on susceptibility to the three insecticides. 102 P450 genes were screened via transcriptome and genome, of which 67 P450 genes were differentially expressed in response to tetraniliprole, spinetoram, and emamectin benzoate and validated by quantitative real-time PCR. The expression patterns of CYP9A75, CYP340AA4, CYP340AX8v2, CYP340L16, CYP341B15v2, and CYP341B17v2 were analysed in different tissues and at different developmental stages in S. frugiperda. Silencing CYP340L16 significantly increased the susceptibility of S. frugiperda to tetraniliprole, spinetoram, and emamectin benzoate. Furthermore, knockdown of CYP340AX8v2, CYP9A75, and CYP341B17v2 significantly increased the sensitivity of S. frugiperda to tetraniliprole. Knockdown of CYP340AX8v2 and CYP340AA4 significantly increased mortality of S. frugiperda to spinetoram. Knockdown of CYP9A75 and CYP341B15v2 significantly increased the susceptibility of S. frugiperda to emamectin benzoate. These results may help to elucidate the mechanisms of tolerance to tetraniliprole, spinetoram and emamectin benzoate in S. frugiperda.

Characterization of emamectin benzoate resistance in the diamondback moth, Plutella xylostella (Lepidoptera: Plutellidae)

BACKGROUND

Plutella xylostella (L.) is a destructive pest of cruciferous crops worldwide that has evolved resistance to many insecticides. Here we examined the mode of inheritance, cross-resistance profile, and potential mechanisms of emamectin benzoate resistance in a field-derived strain of P. xylostella from Japan.

RESULTS

A field-collected population of P. xylostella, was found to exhibit strong (> 150-fold) resistance to emamectin benzoate in insecticide bioassays when compared with a laboratory susceptible strain. Genetic analysis showed that resistance is inherited as an autosomal, recessive trait, and is conferred by a single or a few closely linked loci. The emamectin benzoate resistant strain also exhibited resistance to abamectin, lepimectin, chlorantraniliprole, lufenuron, spinetoram, indoxacarb, fipronil, dieldrin, endosulfan and lambda-cyhalothrin, demonstrating a remarkable multi-resistance profile. Insecticide bioassays employing inhibitors of detoxification enzymes revealed that piperonyl butoxide (PBO) increased the toxicity of emamectin benzoate in the resistant strain by ten-fold indicating the potential involvement of cytochrome P450 monooxygenases in avermectin resistance. Furthermore, cloning and sequencing of the primary receptor of avermectins, the GluCl channel, revealed the absence of target-site mutations in the resistant strain.

CONCLUSIONS

Our data on the mode of inheritance and mechanisms of resistance to emamectin benzoate in a P. xylostella strain from Japan provide a foundation for the development of regional resistance management strategies. However, the high levels of phenotypic resistance in this strain to a diverse range of other insecticide classes available for control illustrate the challenges associated with the sustainable control of this important pest. © 2023 Society of Chemical Industry.

Recent insights into pesticide resistance mechanisms in Plutella xylostella and possible management strategies

Insects are incredibly successful and diverse organisms, but they also pose a significant threat to agricultural crops, causing potential losses of up to US$470 billion. Among these pests, Plutella xylostella (Linnaeus), a devastating insect that attacks cruciferous vegetables, alone results in monetary losses of around US$4-5 billion worldwide. While insecticides have effectively protected plants under field conditions, their use comes with various environmental and mammalian hazards. Additionally, insects are developing resistance to commonly used insecticides, rendering management strategies less effective. Arthropods employ a range of behavioral and biochemical mechanisms to cope with harmful chemicals, which contribute to the development of resistance. Understanding these mechanisms is crucial for addressing the issue of resistance. It is imperative to integrate strategies that can delay the development of resistance and enhance the efficiency of insecticides. Therefore, we present an overview of insecticide resistance in insects, focusing on P. xylostella, to provide insights into the current resistance status of this pest and propose tactics that can improve the effectiveness of insecticides.

A Comparative Study of Transcriptional Regulation Mechanism of Cytochrome P450 CYP6B7 between Resistant and Susceptible Strains of Helicoverpa armigera

Cytochrome P450 CYP6B7 has previously been proved to be associated with fenvalerate-resistance in Helicoverpa armigera. Here, how CYP6B7 is regulated and involved in the resistance of H. armigera is studied. Seven base differences (M1-M7) were found in CYP6B7 promoter between a fenvalerate-resistant (HDTJFR) and a susceptible (HDTJ) strain of H. armigera. M1-M7 sites in HDTJFR were mutated into the corresponding base in HDTJ, and pGL3-CYP6B7 reporter genes with different mutation sites were constructed. Fenvalerate-induced activities of reporter genes mutated at M3, M4, and M7 sites were significantly reduced. Transcription factors Ubx and Br, whose binding sites contain M3 and M7, respectively, were overexpressed in HDTJFR. Knockdown of Ubx and Br results in significant expression inhibition of CYP6B7 and other resistance-related P450 genes, and increase of sensitivity of H. armigera to fenvalerate. These results indicate that Ubx and Br regulate the expression of CYP6B7 to mediate the fenvalerate-resistance in H. armigera.

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

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

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

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

Cytochrome P450 gene CYP6BQ8 mediates terpinen-4-ol susceptibility in the red flour beetle, Tribolium castaneum (Herbst) (Coleoptera: Tenebrionidae)

Cytochrome P450 proteins (CYPs) in insects can encode various detoxification enzymes and catabolize heterologous substances, conferring tolerance to insecticides. This study describes the identification of a P450 gene (CYP6BQ8) from Tribolium castaneum (Herbst) and investigation of its spatiotemporal expression profile and potential role in the detoxification of terpinen-4-ol, a component of plant essential oils. The developmental expression profile showed that TcCYP6BQ8 expression was relatively higher in early- and late-larval stages of T. castaneum compared with other developmental stages. Tissue expression profiles showed that TcCYP6BQ8 was mainly expressed in the head and integument of both larvae and adults. The expression profiling of TcCYP6BQ8 in developmental stages and tissues is closely related to the detoxification of heterologous substances. TcCYP6BQ8 expression was significantly induced after exposure to terpinen-4-ol, and RNA interference against TcCYP6BQ8 increased terpinen-4-ol-induced larval mortality from 47.78 to 66.67%. This indicates that TcCYP6BQ8 may be involved in T. castaneum's metabolism of terpinen-4-ol. Correlation investigation between the CYP6BQ8 gene and terpinen-4-ol resistance in T. castaneum revealed that the TcCYP6BQ8 gene was one of the factors behind T. castaneum's resistance to terpinen-4-ol. This discovery may provide a new theoretical foundation for future regulation of T. castaneum.

Genome-Wide Identification and Expression Pattern of Cytochrome P450 Genes in the Social Aphid Pseudoregma bambucicola

Cytochrome P450 monooxygenases (P450s) have a variety of functions, including involvement in the metabolism of exogenous substances and the synthesis and degradation of endogenous substances, which are important for the growth and development of insects. Pseudoregma bambucicola is a social aphid that produces genetically identical but morphologically and behaviorally distinct first-instar soldiers and normal nymphs within colonies. In this study, we identified 43 P450 genes based on P. bambucicola genome data. Phylogenetic analysis showed that these genes were classified into 4 clans, 13 families, and 23 subfamilies. The CYP3 and CYP4 clans had a somewhat decreased number of genes. In addition, differential gene expression analysis based on transcriptome data showed that several P450 genes, including CYP18A1, CYP4G332, and CYP4G333, showed higher expression levels in soldiers compared to normal nymphs and adult aphids. These genes may be candidates for causing epidermal hardening and developmental arrest in soldiers. This study provides valuable data and lays the foundation for the study of functions of P450 genes in the social aphid P. bambucicola.

Insights into insecticide-resistance mechanisms in invasive species: Challenges and control strategies

Threatening the global community is a wide variety of potential threats, most notably invasive pest species. Invasive pest species are non-native organisms that humans have either accidentally or intentionally spread to new regions. One of the most effective and first lines of control strategies for controlling pests is the application of insecticides. These toxic chemicals are employed to get rid of pests, but they pose great risks to people, animals, and plants. Pesticides are heavily used in managing invasive pests in the current era. Due to the overuse of synthetic chemicals, numerous invasive species have already developed resistance. The resistance development is the main reason for the failure to manage the invasive species. Developing pesticide resistance management techniques necessitates a thorough understanding of the mechanisms through which insects acquire insecticide resistance. Insects use a variety of behavioral, biochemical, physiological, genetic, and metabolic methods to deal with toxic chemicals, which can lead to resistance through continuous overexpression of detoxifying enzymes. An overabundance of enzymes causes metabolic resistance, detoxifying pesticides and rendering them ineffective against pests. A key factor in the development of metabolic resistance is the amplification of certain metabolic enzymes, specifically esterases, Glutathione S-transferase, Cytochromes p450 monooxygenase, and hydrolyses. Additionally, insect guts offer unique habitats for microbial colonization, and gut bacteria may serve their hosts a variety of useful services. Most importantly, the detoxification of insecticides leads to resistance development. The complete knowledge of invasive pest species and their mechanisms of resistance development could be very helpful in coping with the challenges and effectively developing effective strategies for the control of invasive species. Integrated Pest Management is particularly effective at lowering the risk of chemical and environmental contaminants and the resulting health issues, and it may also offer the most effective ways to control insect pests.

Mobilome of the Rhus Gall Aphid Schlechtendalia chinensis Provides Insight into TE Insertion-Related Inactivation of Functional Genes

Transposable elements (TEs) comprise a considerable proportion of insect genomic DNA; how they contribute to genome structure and organization is still poorly understood. Here, we present an analysis of the TE repertoire in the chromosome-level genome assembly of Rhus gall aphid Schlechtendalia chinensis. The TE fractions are composed of at least 32 different superfamilies and many TEs from different families were transcriptionally active in the S. chinensis genome. Furthermore, different types of transposase-derived proteins were also found in the S. chinensis genome. We also provide insight into the TEs related insertional inactivation, and exogenization of TEs in functional genes. We considered that the presence of TE fragments in the introns of functional genes could impact the activity of functional genes, and a large number of TE fragments in introns could lead to the indirect inactivation of functional genes. The present study will be beneficial in understanding the role and impact of TEs in genomic evolution of their hosts.

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

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

Cytochrome P450 Genes Expressed in Phasmatodea Midguts

Cytochrome P450s (CYPs) are xenobiotic detoxification genes found in most eukaryotes, and linked in insects to the tolerance of plant secondary chemicals and insecticide resistance. The number and diversity of CYP clans, families, and subfamilies that an organism produces could correlate with its dietary breadth or specialization. This study examined the CYP diversity expressed in the midguts of six species of folivorous stick insects (Phasmatodea), to identify their CYP complement and see if any CYPs correlate with diet toxicity or specialization, and see what factors influenced their evolution in this insect order. CYP genes were mined from six published Phasmatodea transcriptomes and analyzed phylogenetically. The Phasmatodea CYP complement resembles that of other insects, though with relatively low numbers, and with significant expansions in the CYP clades 6J1, 6A13/14, 4C1, and 15A1. The CYP6 group is known to be the dominant CYP family in insects, but most insects have no more than one CYP15 gene, so the function of the multiple CYP15A1 genes in Phasmatodea is unknown, with neofunctionalization following gene duplication hypothesized. No correlation was found between CYPs and diet specialization or toxicity, with some CYP clades expanding within the Phasmatodea and others likely inherited from a common ancestor.

Genome-Wide Identification of P450 Genes in Chironomid Propsilocerus akamusi Reveals Candidate Genes Involved in Gut Microbiota-Mediated Detoxification of Chlorpyrifos

Chironomids commonly dominate macroinvertebrate assemblages in aquatic habitats and these non-biting midges are known for their ability to tolerate contaminants. Studies regarding the interplay between gut microbiota and host detoxification ability is currently a point of interest. Cytochrome P450s (P450s) are critical metabolic enzymes in which a subset is involved in xenobiotic detoxification. In this study, we first conducted an integrated global investigation of P450s based on the whole genomic sequence of Propsilocerus akamusi and retrieved a series of 64 P450 genes which were further classified into 4 clans and 25 families on the basis of phylogenetic relationships. With assistance of RNA-Seq and RT-qPCR validation, the expression profile of screened PaP450s in guts was compared between chlorpyrifos-challenged larvae with deficient gut microbiota (GD) and those with a conventional gut community (CV). An increasing prevalence of chlorpyrifos from sublethal to lethal dosages induced a greater mortality rate of individuals coupled with remarkable downregulation of 14 P450s in GD larval guts when compared to CV ones. Moreover, it turned out that the decreased level of PaCYP3998B1 and PaCYP3987D1 might imply impaired host endogenous detoxification capability potentiated by gut dysbiosis, reflected by a remarkably severe mortality in GD larvae treated with lethal chlorpyrifos. Collectively, our study unveiled candidate P450 genes that might be mediated by gut symbionts in chlorpyrifos-challenged P. akamusi larvae, possibly facilitating further understanding of the detoxified mechanism that chironomids might employ to alleviate poisonousness.

Functional analysis of a cytochrome P450 gene CYP9Z6 responding to terpinen-4-ol in the red flour beetle, Tribolium castaneum

Tribolium castaneum is an agricultural and stored pest found throughout the world. The cytochrome P450 genes of T. castaneum can encode various detoxification enzymes and catabolize heterologous substances, conferring tolerance to insecticides. Herein, we describe the identification of a P450 gene (CYP9Z6) from T. castaneum and investigated its expression profile and potential role in the detoxification of terpinen-4-ol. TcCYP9Z6 expression was significantly induced after exposure to terpinen-4-ol, and RNA-mediated silencing of TcCYP9Z6 increased terpinen-4-ol-induced larval mortality from 47.75% to 63.92%, showing that TcCYP9Z6 is closely related to the detoxification of terpinen-4-ol. The developmental expression profile revealed that TcCYP9Z6 was mainly expressed in late adults and late larvae. Tissue expression profiling revealed that the highest TcCYP9Z6 expression occurred in the head, in both the adult and the larval tissues, followed by the gut in larvae and the antennae in adults. These developmental stages and tissues with high TcCYP9Z6 expression are closely related to the detoxification of heterologous substances. These results indicated that TcCYP9Z6 may play a pivotal role in the detoxification of terpinen-4-ol, which provides support for using TcCYP9Z6 a potential gene for the RNAi-mediated prevention and control of T. castaneum.

Dynamics of λ-cyhalothrin disappearance and expression of selected P450 genes in bees depending on the ambient temperature

Temperature has a significant influence on the action of pyrethroids, and their effect increases with decreasing ambient temperature. Using gas chromatography, we assessed the degradation rate of λ-cyhalothrin, active ingredients (AI) of Karate Zeon 050 CS from pyrethroid group, in bees incubated for 48 h under different temperature conditions. With RT-qPCR method, we studied expression levels of selected cytochrome P450 genes after exposure to the plant protection product (PPP). The half-life of λ-cyhalothrin decreased from 43.32 to 17.33 h in the temperature range of 21–31°C. In animals incubated at 16°C, the AI half-life was even shorter and amounted to 10.19 h. The increase in temperature increased the expression of Cyp9Q1, Cyp9Q2, and Cyp9Q3 in the group of control bees. We showed a two-fold statistically significant increase in gene expression after treatment with PPP bees. The obtained results indicate that honey bees are characterized by susceptibility to pyrethroids that vary depending on the ambient temperature. This may be due to the different expressions of genes responsible for the detoxification of these PPPs at different temperatures.

A haploid diamondback moth (Plutella xylostella L.) genome assembly resolves 31 chromosomes and identifies a diamide resistance mutation

The diamondback moth, Plutella xylostella (L.), is a highly mobile brassica crop pest with worldwide distribution and can rapidly evolve resistance to insecticides, including group 28 diamides. Reference genomes assembled using Illumina sequencing technology have provided valuable resources to advance our knowledge regarding the biology, origin and movement of diamondback moth, and more recently with its sister species, Plutella australiana. Here we apply a trio binning approach to sequence and annotate a chromosome level reference genome of P. xylostella using PacBio Sequel and Dovetail Hi-C sequencing technology and identify a point mutation that causes resistance to commercial diamides. A P. xylostella population collected from brassica crops in the Lockyer Valley, Australia (LV-R), was reselected for chlorantraniliprole resistance then a single male was crossed to a P. australiana female and a hybrid pupa sequenced. A chromosome level 328 Mb P. xylostella genome was assembled with 98.1% assigned to 30 autosomes and the Z chromosome. The genome was highly complete with 98.4% of BUSCO Insecta genes identified and RNAseq informed protein prediction annotated 19,002 coding genes. The LV-R strain survived recommended field application doses of chlorantraniliprole, flubendiamide and cyclaniliprole. Some hybrids also survived these doses, indicating significant departure from recessivity, which has not been previously documented for diamides. Diamide chemicals modulate insect Ryanodine Receptors (RyR), disrupting calcium homeostasis, and we identified an amino acid substitution (I4790K) recently reported to cause diamide resistance in a strain from Japan. This chromosome level assembly provides a new resource for insect comparative genomics and highlights the emergence of diamide resistance in Australia. Resistance management plans need to account for the fact that resistance is not completely recessive.

Resistance in the Genus Spodoptera: Key Insect Detoxification Genes

The genus Spodoptera (Lepidoptera: Noctuidae) includes species that are among the most important crop pests in the world. These polyphagous species are able to feed on many plants, including corn, rice and cotton. In addition to their ability to adapt to toxic compounds produced by plants, they have developed resistance to the chemical insecticides used for their control. One of the main mechanisms developed by insects to become resistant involves detoxification enzymes. In this review, we illustrate some examples of the role of major families of detoxification enzymes such as cytochromes P450, carboxyl/cholinesterases, glutathione S-transferases (GST) and transporters such as ATP-binding cassette (ABC) transporters in insecticide resistance. We compare available data for four species, Spodoptera exigua, S. frugiperda, S. littoralis and S. litura. Molecular mechanisms underlying the involvement of these genes in resistance will be described, including the duplication of the CYP9A cluster, over-expression of GST epsilon or point mutations in acetylcholinesterase and ABCC2. This review is not intended to be exhaustive but to highlight the key roles of certain genes.

Interaction between Metarhizium anisopliae and Its Host, the Subterranean Termite Coptotermes curvignathus during the Infection Process

Metarhizium anisopliae (Metchnikoff) Sorokin, a pathogenic fungus to insects, infects the subterranean termite, Coptotermes curvignathus Holmgren, a devastating pest of plantation trees in the tropics. Electron microscopy and proteomics were used to investigate the infection and developmental process of M. anisopliae in C. curvignathus. Fungal infection was initiated by germ tube penetration through the host's cuticle as observed at 6 h post-inoculation (PI), after which it elongated into the host's integumental tissue. The colonization process continued as seen from dissemination of blastospores in the hemocoel at 96 h PI. At this time point, the emergent mycelia had mummified the host and forty-eight hours later, new conidia were dispersed on the termites' body surface. Meanwhile, hyphal bodies were observed in abundance in the intercellular space in the host's body. The proteomes of the pathogen and host were isolated separately using inoculated termite samples withdrawn at each PI-time point and analyzed in two-dimensional electrophoresis (2-DE) gels. Proteins expressed in termites showed evidence of being related to cell regulation and the immune response, while those expressed in M. anisopliae, to transportation and fungal virulence. This study provides new information on the interaction between termites and its entomopathogen, with potential utilization for developing future biopesticide to control the termite population.

Cytochrome P450 monooxygenase genes in the wild silkworm, Bombyx mandarina

Wild (Bombyx mandarina) and domestic silkworms (B. mori) are good models for investigating insect domestication, as 5000 years of artificial breeding and selection have resulted in significant differences between B. mandarina and B. mori. In this study, we improved the genome assemblies to the chromosome level and updated the protein-coding gene annotations for B. mandarina. Based on this updated genome, we identified 68 cytochrome P450 genes in B. mandarina. The cytochrome P450 repository in B. mandarina is smaller than in B. mori. Certain currently unknown key genes, rather than gene number, are critical for insecticide resistance in B. mandarina, which shows greater resistance to insecticides than B. mori. Based on the physical maps of B. mandarina, we located 66 cytochrome P450s on 18 different chromosomes, and 27 of the cytochrome P450 genes were concentrated into seven clusters. KEGG enrichment analysis of the P450 genes revealed the involvement of cytochrome P450 genes in hormone biosynthesis. Analyses of the silk gland transcriptome identified candidate cytochrome P450 genes (CYP306A) involved in ecdysteroidogenesis and insecticide metabolism in B. mandarina.

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.

Sequencing of Transcriptome and Small RNA Revealed the Xenobiotic Metabolism-Related Genes and Potential Regulatory miRNA in Asian Tramp Snail

The Asian tramp snail, Bradybaena similaris (Ferusssac), is an invasive land snail species and has been a rising agricultural pest in south of China. As a pest, it also plays a role in transmission of Angiostrongylus cantonensis. However, present studies on this species are rare and the molecular information is limited. For this purpose, we sequenced the transcriptome and small RNA of B. similaris collected from citrus orchards. In total, 89,747 unigenes with an N50 size of 1287 bp and an average length of 817 bp were generated from ∼8.9 Gb transcriptome and 31 Mb clean reads were generated from ∼36 Mb small RNA library. To demonstrate the usefulness of these two datasets, we analyzed a series of genes associated with xenobiotic metabolism and core RNAi machinery. Analysis of the transcripts resulted in annotation of 126 putative genes encoding cytochrome P450 monooxygenases (CYP, 45), carboxyl/cholinesterases (CCE, 13), glutathione-S-transferases (GST, 24), and ATP-binding cassette transporters (ABC, 44). Analysis of the small RNA detected 42 miRNAs. In addition, four genes involved in small RNA pathways (miRNA, piRNA, and siRNA) were identified, and a total of 430 genes that can be targeted by miRNAs were predicted. Moreover, we found that a few miRNAs could target certain genes involved in xenobiotic metabolism. Therefore, we believe that these two datasets and the characterization of the identified/predicted genes will facilitate the molecular study of this species as well as other land snails with agricultural importance.

Xenobiotic transcription factors CncC and maf regulate expression of CYP321A16 and CYP332A1 that mediate chlorpyrifos resistance in Spodoptera exigua

Insect cytochrome P450 s (P450 s) are associated with the metabolic detoxification of toxic xenobiotics and their constitutive upregulation is often associated with resistance to natural and synthetic toxicants. The P450 s CYP321A16 and CYP332A1 are constitutively overexpressed in an insecticide-resistant strain of beet armyworm, Spodoptera exigua. However, the function and upstream regulation of these two P450 s remain unknown. Here, we investigated the function of CYP321A16 and CYP332A1 in resistance to the organophosphate insecticide, chlorpyrifos and their regulation by the transcription factors CncC and Maf. Transgenic strains of Drosophila melanogaster expressing CYP321A16 or CYP332A1 showed higher levels of tolerance to chlorpyrifos than the control flies with the same genetic background. Furthermore, recombinant CYP321A16 and CYP332A1 proteins metabolized chlorpyrifos. Analysis of the putative promoter sequences of the genes coding for CYP321A16 and CYP332A1 revealed conserved CncC/Maf binding sites. Transfection of luciferase reporter plasmids containing the promoter of CYP450 gene together with CncC and Maf expression plasmids significantly enhanced the activity of the reporter. Promoter truncation identified a site in the promoter of CYP321A16 that is critical for the CncC/Maf binding. These data demonstrate that resistance to chlorpyrifos in S. exigua is conferred by the combined action of CYP321A16 and CYP332A1 and uncovered their regulation by the transcription factors CncC and Maf.

The genome of pest Rhynchophorus ferrugineus reveals gene families important at the plant-beetle interface

The red palm weevil, Rhynchophorus ferrugineus, infests palm plantations, leading to large financial losses and soil erosion. Pest-host interactions are poorly understood in R. ferrugineus, but the analysis of genetic diversity and pest origins will help advance efforts to eradicate this pest. We sequenced the genome of R. ferrugineus using a combination of paired-end Illumina sequencing (150 bp), Oxford Nanopore long reads, 10X Genomics and synteny analysis to produce an assembly with a scaffold N50 of ~60 Mb. Structural variations showed duplication of detoxifying and insecticide resistance genes (e.g., glutathione S-transferase, P450, Rdl). Furthermore, the evolution of gene families identified those under positive selection including one glycosyl hydrolase (GH16) gene family, which appears to result from horizontal gene transfer. This genome will be a valuable resource to understand insect evolution and behavior and to allow the genetic modification of key genes that will help control this pest.

Genome-wide characterization and transcriptomic analyses of neuropeptides and their receptors in an endoparasitoid wasp, Pteromalus puparum

In insects, neuropeptides constitute a group of signaling molecules that act in regulation of multiple physiological and behavioral processes by binding to their corresponding receptors. On the basis of the bioinformatic approaches, we screened the genomic and transcriptomic data of the parasitoid wasp, Pteromalus puparum, and annotated 36 neuropeptide precursor genes and 33 neuropeptide receptor genes. Compared to the number of precursor genes in Bombyx mori (Lepidoptera), Chilo suppressalis (Lepidoptera), Drosophila melanogaster (Diptera), Nilaparvata lugens (Hemiptera), Apis mellifera (Hymenoptera), and Tribolium castaneum (Coleoptera), P. puparum (Hymenoptera) has the lowest number of neuropeptide precursor genes. This lower number may relate to its parasitic life cycle. Transcriptomic data of embryos, larvae, pupae, adults, venom glands, salivary glands, ovaries, and the remaining carcass revealed stage-, sex-, and tissue-specific expression patterns of the neuropeptides, and their receptors. These data provided basic information about the identity and expression profiles of neuropeptides and their receptors that are required to functionally address their biological significance in an endoparasitoid wasp.

Substrate-specificity of cytochrome P450-mediated detoxification as an evolutionary strategy for specialization on furanocoumarin-containing hostplants: CYP6AE89 in parsnip webworms

The parsnip webworm, Depressaria pastinacella, is restricted to two hostplant genera containing six structurally diverse furanocoumarins. Of these, imperatorin is detoxified by a specialized cytochrome P450, CYP6AB3. A previous whole-larva transcriptome analysis confirmed the presence of nine transcripts that belong to the CYP6AE subfamily. Here, by examining midgut-specific gene expression patterns we determined that CYP6AE89 transcripts were highly expressed and furanocoumarin-inducible. Computer docking and energy-minimization of a CYP6AE89 model with all six furanocoumarins showed that 5-methoxylated bergapten and 8-methoxylated xanthotoxin had the smallest distances from the heme to the proton-donor residue in the catalytic I-helix, and that the 5,8-dimethoxylated isopimpinellin and bergapten had the smallest energy-minimized distance from the heme oxygen to the furan ring double bond. To evaluate this prediction, we expressed the CYP6AE89 protein in an Escherichia coli system, and used it to detect high catalytic activity against the two mono-methoxylated linear furanocoumarins - bergapten and xanthotoxin - and weak activity against isopimpinellin. Thus, CYP6AE89, like CYP6AB3, is probably specialized for detoxifying only a subset of hostplant furanocoumarins. A maximum-likelihood tree built with six representative lepidopterans with manually annotated cytochrome P450s shows that CYP6AE89 may have evolved much faster than the other CYP6AE proteins, possibly indicative of host selection pressure.

Comparative transcriptome profiling reveals candidate genes related to insecticide resistance of Glyphodes pyloalis

Glyphodes pyloalis Walker (Lepidoptera: Pyralididae) is a common pest in sericulture and has developed resistance to different insecticides. However, the mechanisms involved in insecticide resistance of G. pyloalis are poorly understood. Here, we present the first whole-transcriptome analysis of differential expression genes in insecticide-resistant and susceptible G. pyloalis. Clustering and enrichment analysis of DEGs revealed several biological pathways and enriched Gene Ontology terms were related to detoxification or insecticide resistance. Genes involved in insecticide metabolic processes, including cytochrome P450, glutathione S-transferases and carboxylesterase, were identified in the larval midgut of G. pyloalis. Among them, CYP324A19, CYP304F17, CYP6AW1, CYP6AB10, GSTs5, and AChE-like were significantly increased after propoxur treatment, while CYP324A19, CCE001c, and AChE-like were significantly induced by phoxim, suggesting that these genes were involved in insecticide metabolism. Furthermore, the sequence variation analysis identified 21 single nucleotide polymorphisms within CYP9A20, CYP6AB47, and CYP6AW1. Our findings reveal many candidate genes related to insecticide resistance of G. pyloalis. These results provide novel insights into insecticide resistance and facilitate the development of insecticides with greater specificity to G. pyloalis.

Identification of Differentially Expressed Proteins in Sugarcane in Response to Infection by Xanthomonas albilineans Using iTRAQ Quantitative Proteomics

Sugarcane can suffer severe yield losses when affected by leaf scald, a disease caused by Xanthomonas albilineans. This bacterial pathogen colonizes the vascular system of sugarcane, which can result in reduced plant growth and plant death. In order to better understand the molecular mechanisms involved in the resistance of sugarcane to leaf scald, a comparative proteomic study was performed with two sugarcane cultivars inoculated with X. albilineans: one resistant (LCP 85-384) and one susceptible (ROC20) to leaf scald. The iTRAQ (isobaric tags for relative and absolute quantification) approach at 0 and 48 h post-inoculation (hpi) was used to identify and annotate differentially expressed proteins (DEPs). A total of 4295 proteins were associated with 1099 gene ontology (GO) terms by GO analysis. Among those, 285 were DEPs during X. albilineans infection in cultivars LCP 85-384 and ROC20. One hundred seventy-two DEPs were identified in resistant cultivar LCP 85-384, and 113 of these proteins were upregulated and 59 were downregulated. One hundred ninety-two DEPs were found in susceptible cultivar ROC20 and half of these (92) were upregulated, whereas the other half corresponded to downregulated proteins. The significantly upregulated DEPs in LCP 85-384 were involved in metabolic pathways, the biosynthesis of secondary metabolites, and the phenylpropanoid biosynthesis pathway. Additionally, the expression of seven candidate genes related to photosynthesis and glycolytic pathways, plant innate immune system, glycosylation process, plant cytochrome P450, and non-specific lipid transfer protein was verified based on transcription levels in sugarcane during infection by X. albilineans. Our findings shed new light on the differential expression of proteins in sugarcane cultivars in response to infection by X. albilineans. The identification of these genes provides important information for sugarcane variety improvement programs using molecular breeding strategies.

De novo transcriptome analysis and identification of genes associated with immunity, detoxification and energy metabolism from the fat body of the tephritid gall fly, Procecidochares utilis

The fat body, a multifunctional organ analogous to the liver and fat tissue of vertebrates, plays an important role in insect life cycles. The fat body is involved in protein storage, energy metabolism, elimination of xenobiotics, and production of immunity regulator-like proteins. However, the molecular mechanism of the fat body's physiological functions in the tephritid stem gall-forming fly, Procecidochares utilis, are still unknown. In this study, we performed transcriptome analysis of the fat body of P. utilis using Illumina sequencing technology. In total, 3.71 G of clean reads were obtained and assembled into 30,559 unigenes, with an average length of 539 bp. Among those unigenes, 21,439 (70.16%) were annotated based on sequence similarity to proteins in NCBI's non-redundant protein sequence database (Nr). Sequences were also compared to NCBI's non-redundant nucleotide sequence database (Nt), a manually curated and reviewed protein sequence database (SwissProt), and KEGG and gene ontology annotations were applied to better understand the functions of these unigenes. A comparative analysis was performed to identify unigenes related to detoxification, immunity and energy metabolism. Many unigenes involved in detoxification were identified, including 50 unigenes of putative cytochrome P450s (P450s), 18 of glutathione S-transferases (GSTs), 35 of carboxylesterases (CarEs) and 26 of ATP-binding cassette (ABC) transporters. Many unigenes related to immunity were identified, including 17 putative serpin genes, five peptidoglycan recognition proteins (PGRPs) and four lysozyme genes. In addition, unigenes potentially involved in energy metabolism, including 18 lipase genes, five fatty acid synthase (FAS) genes and six elongases of very long chain fatty acid (ELOVL) genes, were identified. This transcriptome improves our genetic understanding of P. utilis and the identification of a numerous transcripts in the fat body of P. utilis offer a series of valuable molecular resources for future studies on the functions of these genes.

The Roles of Four Novel P450 Genes in Pesticides Resistance in Apis cerana cerana Fabricius: Expression Levels and Detoxification Efficiency

Cytochrome P450 monooxygenases (P450s) are widely distributed multifunctional enzymes that play crucial roles in insecticide detoxification or activation. In this study, to ascertain the molecular mechanisms of P450s in the detoxification of Chinese honeybees, Apis cerana cerana Fabricius (A. c. cerana), we isolated and characterized four new P450 genes (Acc301A1, Acc303A1, Acc306A1, and Acc315A1). The open reading frames of the four genes are 1263 to 1608 bp in length and encode four predicted polypeptides of 499 to 517 amino acids in length. Real-time quantitative PCR (RT-qPCR) results showed that expression of all four genes was observed in all developmental stages. In addition, Western blot assays further indicated the RT-qPCR results that showed that the four genes were induced by pesticide (thiamethoxam, deltamethrin, dichlorovos, and paraquat) treatments. Furthermore, we also used double-stranded RNA-mediated RNA interference to investigate the functions of Acc301A1, Acc303A1,and Acc306A1 in the antioxidant defense of honeybees. RNA interference targeting Acc301A1, Acc303A1, and Acc306A1 significantly increased the mortality rate of A. c. cerana upon pesticide treatment. These results provide important evidence about the role of the four P450 genes involved in detoxification.

The genome of the marine monogonont rotifer Brachionus plicatilis: Genome-wide expression profiles of 28 cytochrome P450 genes in response to chlorpyrifos and 2-ethyl-phenanthrene

Brachionus spp. (Rotifera: Monogononta) are globally distributed in aquatic environments and play important roles in the aquatic ecosystem. The marine monogonont rotifer Brachionus plicatilis is considered a suitable model organism for ecology, evolution, and ecotoxicology. In this study, we assembled and characterized the B. plicatilis genome. The total length of the assembled genome was 106.9 Mb and the number of final scaffolds was 716 with an N50 value of 1.15 Mb and a GC content of 26.75%. A total of 20,154 genes were annotated after manual curation. To demonstrate the use of whole genome data, we targeted one of the main detoxifying enzyme of phase I detoxification system and identified in a total of 28 cytochrome P450 s (CYPs). Based on the phylogenetic analysis using the maximum likelihood, 28 B. plicatilis-CYPs were apparently separated into five different clans, namely, 2, 3, 4, mitochondrial (MT), and 46 clans. To better understand the CYPs-mediated xenobiotic detoxification, we measured the mRNA expression levels of 28 B. plicatilis CYPs in response to chlorpyrifos and 2-ethyl-phenanthrene. Most B. plicatilis CYPs were significantly modulated (P < 0.05) in response to chlorpyrifos and 2-ethyl-phenanthrene. In addition, xenobiotic-sensing nuclear receptor (XNR) response element sequences were identified in the 5 kb upstream of promoter regions of 28 CYPs from the genome of B. plicatilis, indicating that these XNR can be associated with detoxification of xenobiotics. Overall, the assembled B. plicatilis genome presented here will be a useful resource for a better understanding the molecular ecotoxicology in the view of molecular mechanisms underlying toxicological responses, particularly on xenobiotic detoxification in this species.

Differential Expression of Cytochrome P450 CYP6 Genes in the Brown Marmorated Stink Bug, Halyomorpha halys (Hemiptera: Pentatomidae)

Cytochrome (CYP) P450s are a superfamily of enzymes that detoxify xenobiotics and regulate numerous physiological processes in insects. The genomes of phytophagous insects usually contain large numbers of P450s, especially within the CYP3 clan. Within this clan, CYP6 subfamily members help detoxify plant host secondary metabolites. In this study, we analyzed three CYP6 genes in the highly polyphagous invasive pest, Halyomorpha halys (Stål), commonly known as brown marmorated stink bug. We characterized and validated the expression of HhCYP6BQ27, HhCYP6BK13, and HhCYP6BK24 among sexes, tissues (gut, fat body, and Malpighian tubules) and hosts (apple, corn, soybean). Sequence characterization by amino acid alignments confirmed the presence of conserved motifs typical of the P450 superfamily. No significant differences existed in gene expression among sexes or when fed different hosts, suggesting that these transcripts might have broad substrate specificities. However, significant differences in gene expression were observed among the tissues studied and were gene-dependent. Collectively, the results show that H. halys differentially expressed CYP6 genes among tissues, which may be related to important and specific physiological functions. This study has increased our understanding of H. halys biology that can be useful for functional studies and can potentially be exploited in developing sustainable pest management strategies.

Boosting the biosynthesis of betulinic acid and related triterpenoids in Yarrowia lipolytica via multimodular metabolic engineering

BACKGROUND

Betulinic acid is a pentacyclic lupane-type triterpenoid and a potential antiviral and antitumor drug, but the amount of betulinic acid in plants is low and cannot meet the demand for this compound. Yarrowia lipolytica, as an oleaginous yeast, is a promising microbial cell factory for the production of highly hydrophobic compounds due to the ability of this organism to accumulate large amounts of lipids that can store hydrophobic products and supply sufficient precursors for terpene synthesis. However, engineering for the heterologous production of betulinic acid and related triterpenoids has not developed as systematically as that for the production of other terpenoids, thus the production of betulinic acid in microbes remains unsatisfactory.

RESULTS

In this study, we applied a multimodular strategy to systematically improve the biosynthesis of betulinic acid and related triterpenoids in Y. lipolytica by engineering four functional modules, namely, the heterogenous CYP/CPR, MVA, acetyl-CoA generation, and redox cofactor supply modules. First, by screening 25 combinations of cytochrome P450 monooxygenases (CYPs) and NADPH-cytochrome P450 reductases (CPRs), each of which originated from 5 different sources, we selected two optimal betulinic acid-producing strains. Then, ERG1, ERG9, and HMG1 in the MVA module were overexpressed in the two strains, which dramatically increased betulinic acid production and resulted in a strain (YLJCC56) that exhibited the highest betulinic acid yield of 51.87 ± 2.77 mg/L. Then, we engineered the redox cofactor supply module by introducing NADPH- or NADH-generating enzymes and the acetyl-CoA generation module by directly overexpressing acetyl-CoA synthases or reinforcing the β-oxidation pathway, which further increased the total triterpenoid yield (the sum of the betulin, betulinic acid, betulinic aldehyde yields). Finally, we engineered these modules in combination, and the total triterpenoid yield reached 204.89 ± 11.56 mg/L (composed of 65.44% betulin, 23.71% betulinic acid and 10.85% betulinic aldehyde) in shake flask cultures.

CONCLUSIONS

Here, we systematically engineered Y. lipolytica and achieved, to the best of our knowledge, the highest betulinic acid and total triterpenoid yields reported in microbes. Our study provides a suitable reference for studies on heterologous exploitation of P450 enzymes and manipulation of triterpenoid production in Y. lipolytica.

The expression of Spodoptera exigua P450 and UGT genes: tissue specificity and response to insecticides

Cytochrome P450 and UDP-glucosyltransferase (UGT) as phase I and phase II metabolism enzymes, respectively, play vital roles in the breakdown of endobiotics and xenobiotics. Insects can increase the expression of detoxification enzymes to cope with the stress from xenobiotics including insecticides. However, the molecular mechanisms for insecticide detoxification in Spodoptera exigua remain elusive, and the genes conferring insecticide metabolisms in this species are less well reported. In this study, 68 P450 and 32 UGT genes were identified. Phylogenetic analysis showed gene expansions in CYP3 and CYP4 clans of P450 genes and UGT33 family of this pest. P450 and UGT genes exhibited specific tissue expression patterns. Insecticide treatments in fat body cells of S. exigua revealed that the expression levels of P450 and UGT genes were significantly influenced by challenges of abamectin, lambda-cyhalothrin, chlorantraniliprole, metaflumizone and indoxacarb. Multiple genes for detoxification were affected in expression levels after insecticide exposures. The results demonstrated that lambda-cyhalothrin, chlorantraniliprole, metaflumizone and indoxacarb induced similar responses in the expression of P450 and UGT genes in fat body cells; eight P450 genes and four UGT genes were co-up-regulated significantly, and no or only a few CYP/UGT genes were down-regulated significantly by these four insecticides. However, abamectin triggered a distinct response for P450 and UGT gene expression; more P450 and UGT genes were down-regulated by abamectin than by the other four compounds. In conclusion, P450 and UGT genes from S. exigua were identified, and different responses to abamectin suggest a different mechanism for insecticide detoxification.

Diversity and evolution of cytochromes P450 in stramenopiles

MAIN CONCLUSION

Comparative genomic analysis of cytochromes P450 revealed high diversification and dynamic changes in stramenopiles, associated with transcriptional responsiveness to various environmental stimuli. Comparative genomic and molecular evolution approaches were used to characterize cytochromes P450 (P450) diversity in stramenopiles. Phylogenetic analysis pointed to a high diversity of P450 in stramenopiles and identified three major clans. The CYP51 and CYP97 clans were present in brown algae, diatoms and Nannochloropsis gaditana, whereas the CYP5014 clan mainly includes oomycetes. Gene gain and loss patterns revealed that six CYP families-CYP51, CYP97, CYP5160, CYP5021, CYP5022, and CYP5165-predated the split of brown algae and diatoms. After they diverged, diatoms gained more CYP families, especially in the cold-adapted species Fragilariopsis cylindrus, in which eight new CYP families were found. Selection analysis revealed that the expanded CYP51 family in the brown alga Cladosiphon okamuranus exhibited a more relaxed selection constraint compared with those of other brown algae and diatoms. Our RNA-seq data further evidenced that most of P450s in Saccharina japonica are highly expressed in large sporophytes, which could potentially promote the large kelp formation in this developmental stage. A survey of Ectocarpus siliculosus and diatom transcriptomes showed that many P450s are responsive to stress, nutrient limitation or light quality, suggesting pivotal roles in detoxification or metabolic processes under adverse environmental conditions. The information provided in this study will be helpful in designing functional experiments and interpreting P450 roles in this particular lineage.

Structure and expression of sulfatase and sulfatase modifying factor genes in the diamondback moth, Plutella xylostella

The diamondback moth, Plutella xylostella (L.), uses sulfatases (SULF) to counteract the glucosinolate-myrosinase defensive system that cruciferous plants have evolved to deter insect feeding. Sulfatase activity is regulated by post-translational modification of a cysteine residue by sulfatase modifying factor 1 (SUMF1). We identified 12 SULF genes (PxylSulfs) and two SUMF1 genes (PxylSumf1s) in the P. xylostella genome. Phylogenetic analysis of SULFs and SUMFs from P. xylostella, Bombyx mori, Manduca sexta, Heliconius melpomene, Danaus plexippus, Drosophila melanogaster, Tetranychus urticae and Homo sapiens showed that the SULFs were clustered into five groups, and the SUMFs could be divided into two groups. Profiling of the expression of PxylSulfs and PxylSumfs by RNA-seq and by quantitative real-time polymerase chain reaction showed that two glucosinolate sulfatase genes (GSS), PxylSulf2 and PxylSulf3, were primarily expressed in the midgut of 3rd- and 4th-instar larvae. Moreover, expression of sulfatases PxylSulf2, PxylSulf3 and PxylSulf4 were correlated with expression of the sulfatases modifying factor PxylSumf1a. The findings from this study provide new insights into the structure and expression of SUMF1 and PxylSulf genes that are considered to be key factors for the evolutionary success of P. xylostella as a specialist herbivore of cruciferous plants.

Mechanisms of resistance to thiamethoxam and dinotefuran compared to imidacloprid in the brown planthopper: Roles of cytochrome P450 monooxygenase and a P450 gene CYP6ER1

The brown planthopper (BPH, Nilaparvata lugens) has developed high resistance to the first-generation neonicotinoids (imidacloprid). With commercialization and widespread field use of the second-(thiamethoxam) and third-(dinotefuran) generation neonicotinoids, resistance to these insecticides is also reported. We investigated the cytochrome P450 monooxygenase-mediated detoxification in thiamethoxam- and dinotefuran- resistant in comparison to imidacloprid-resistant strains of BPH. In the three moderately resistant BPH strains selected separately with the three insecticides from a same susceptible strain, P450 activities were significantly enhanced over the susceptible control. Seven of 26 tested P450 genes were up-regulated and CYP6ER1 was a strongly over-expressed gene in all the three resistant strains. Knockdown of CYP6ER1 in the susceptible insects reduced P450 activity, retarded nymph growth and significantly increased sensitivity to each one of the three neonicotinoids. Taken together, we show that enhanced P450 activity and over-expression of CYP6ER1 gene are involved in BPH resistance to thiamethoxam and dinotefuran as to imidacloprid. These findings are of significance in management thiamethoxam and dinotefuran resistance in the BPH, especially in the management of potential cross-resistance to the three generations of neonicotinoids.

De novo transcriptomic analysis of the alimentary tract of the tephritid gall fly, Procecidochares utilis

The tephritid gall fly, Procecidochares utilis, is an important obligate parasitic insect of the malignant weed Eupatorium adenophorum which biosynthesizes toxic secondary metabolites. Insect alimentary tracts secrete several enzymes that are used for detoxification, including cytochrome P450s, glutathione S-transferases, and carboxylesterases. To explore the adaptation of P. utilis to its toxic host plant, E. adenophorum at molecular level, we sequenced the transcriptome of the alimentary tract of P. utilis using Illumina sequencing. Sequencing and de novo assembly yielded 62,443 high-quality contigs with an average length of 604 bp that were further assembled into 45,985 unigenes with an average length of 674 bp and an N50 of 983 bp. Among the unigenes, 30,430 (66.17%) were annotated by alignment against the NCBI non-redundant protein (Nr) database, while 16,700 (36.32%), 16,267 (35.37%), and 11,530 (25.07%) were assigned functions using the Clusters of Orthologous Groups (COG), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Ontology (GO) databases, respectively. Using the comprehensive transcriptome data set, we manually identified several important gene families likely to be involved in the detoxification of toxic compounds including 21 unigenes within the glutathione S-transferase (GST) family, 22 unigenes within the cytochrome P450 (P450) family, and 16 unigenes within the carboxylesterase (CarE) family. Quantitative PCR was used to verify eight, six, and two genes of GSTs, P450s, and CarEs, respectively, in different P. utilis tissues and at different developmental stages. The detoxification enzyme genes were mainly expressed in the foregut and midgut. Moreover, the unigenes were higher expressed in the larvae, pupae, and 3-day adults, while they were expressed at lower levels in eggs. These transcriptomic data provide a valuable molecular resource for better understanding the function of the P. utilis alimentary canal. These identified genes could be pinpoints to address the molecular mechanisms of P. utilis interacting with toxic plant host.

Genome-wide and expression-profiling analyses suggest the main cytochrome P450 genes related to pyrethroid resistance in the malaria vector, Anopheles sinensis (Diptera Culicidae)

BACKGROUND

Anopheles sinensis is one of the major malaria vectors. However, pyrethroid resistance in An. sinensis is threatening malaria control. Cytochrome P450-mediated detoxification is an important pyrethroid resistance mechanism that has been unexplored in An. sinensis. In this study, we performed a comprehensive analysis of the An. sinensis P450 gene superfamily with special attention to their role in pyrethroid resistance using bioinformatics and molecular approaches.

RESULTS

Our data revealed the presence of 112 individual P450 genes in An. sinensis, which were classified into four major clans (mitochondrial, CYP2, CYP3 and CYP4), 18 families and 50 subfamilies. Sixty-seven genes formed nine gene clusters, and genes within the same cluster and the same gene family had a similar gene structure. Phylogenetic analysis showed that most of An. sinensis P450s (82/112) had very close 1: 1 orthology with Anopheles gambiae P450s. Five genes (AsCYP6Z2, AsCYP6P3v1, AsCYP6P3v2, AsCYP9J5 and AsCYP306A1) were significantly upregulated in three pyrethroid-resistant populations in both RNA-seq and RT-qPCR analyses, suggesting that they could be the most important P450 genes involved in pyrethroid resistance in An. sinensis.

CONCLUSION

Our study provides insight on the diversity of An. sinensis P450 superfamily and basis for further elucidating pyrethroid resistance mechanism in this mosquito species. © 2018 Society of Chemical Industry.

Refining a steroidogenic model: an analysis of RNA-seq datasets from insect prothoracic glands

BACKGROUND

The prothoracic gland (PG), the principal steroidogenic organ of insects, has been proposed as a model for steroid hormone biosynthesis and regulation.

RESULTS

To validate the robustness of the model, we present an analysis of accumulated transcriptomic data from PGs of two model species, Drosophila melanogaster and Bombyx mori. We identify that the common core components of the model in both species are encoded by nine genes. Five of these are Halloween genes whose expression differs substantially between the PGs of these species.

CONCLUSIONS

We conclude that the PGs can be a model for steroid hormone synthesis and regulation within the context of mitochondrial cholesterol transport and steroid biosynthesis but beyond these core mechanisms, gene expression in insect PGs is too diverse to fit in a context-specific model and should be analysed within a species-specific framework.

Overexpression of cytochrome P450 CYP6BG1 may contribute to chlorantraniliprole resistance in Plutella xylostella (L.)

BACKGROUND

The diamondback moth Plutella xylostella (L.) is the most widely distributed pest of cruciferous crops and has developed resistance to most commonly used insecticides, including chlorantraniliprole. Resistance to chlorantraniliprole is likely caused by mutations of the target, the ryanodine receptor, and/or mediated by an increase in detoxification enzyme activities. Although target-site resistance is documented in detail, resistance mediated by increased metabolism has rarely been reported.

RESULTS

The activity of cytochrome P450 was significantly higher in two resistant P. xylostella populations than in a susceptible one. Among ten detected cytochrome P450 genes, CYP6BG1 was significantly overexpressed (over 80-fold) in a field-resistant population compared with expression in a susceptible one. Knockdown of CYP6BG1 by RNA interference dramatically reduced the 7-ethoxycoumarin-O-deethylase (7-ECOD) activity of P450 by 45.5% and increased the toxicity of chlorantraniliprole toward P. xylostella by 26.8% at 48 h postinjection of double-stranded RNA. By contrast, overexpression of CYP6BG1 in a transgenic Drosophila melanogaster line significantly decreased the toxicity of the insecticide to the transgenic flies.

CONCLUSIONS

Overexpression of CYP6BG1 may contribute to chlorantraniliprole resistance in P. xylostella. Our findings will provide new insights into the mechanisms of resistance to diamide insecticides in other insect pests. © 2017 Society of Chemical Industry.

Characterization of UDP-glucuronosyltransferase genes and their possible roles in multi-insecticide resistance in Plutella xylostella (L.)

BACKGROUND

Uridine diphosphate-glucuronosyltransferases (UGTs), as multifunctional detoxification enzymes, play important roles in the biotransformation of various compounds. However, their roles in insecticide resistance are still unclear. This study presents a genome-wide identification of the UGTs in diamondback moth, Plutella xylostella (L.), a notorious insect pest of cruciferous crops worldwide. The possible roles of these UGTs in insecticide resistance were evaluated.

RESULTS

A total of 21 putative UGTs in P. xylostella were identified. Quantitative real-time polymerase chain reaction (PCR)-based analyses showed that all the UGT genes were expressed in all tested developmental stages and tissues. Bioassay results indicated that a field-collected population (BL) was resistant to 9 of 10 commonly used insecticides, and 10 of 21 UGT mRNAs were upregulated in the BL population. Exposure to the LC₅₀ of each insecticide affected the expression of most UGT genes. Among these, the expression levels of UGT40V1, UGT45B1 and UGT33AA4 were induced by more than five insecticides, whereas indoxacarb and metaflumizone significantly repressed the expression of most UGT genes.

CONCLUSION

UGTs may play important roles in the metabolism of commonly used insecticides in P. xylostella. These findings provide valuable information for further research on the physiological and toxicological functions of specific UGT genes in P. xylostella. © 2017 Society of Chemical Industry.

Phylogenetic and functional characterization of ten P450 genes from the CYP6AE subfamily of Helicoverpa armigera involved in xenobiotic metabolism

The cotton bollworm, Helicoverpa armigera, is a generalist herbivore widely distributed over the world and is a major lepidopteran pest on cotton. Studies, especially from Asia, show that it relies on cytochrome P450 monooxygenases with broad substrate specificities to protect itself from pesticides. The number of P450s may have expanded in the processes of coping with the wide diversity of phytochemicals that the insect encounters among its numerous host plants. In order to examine the metabolic capabilities of these P450s, we focused here on all ten P450s of the Helicoverpa armigera CYP6AE subfamily, which can be easily induced by plant toxins and pyrethroids. These P450s, along with cytochrome P450 reductase, were heterologously expressed in insect cells and compared functionally. In vitro metabolism showed that all CYP6AE subfamily members can convert esfenvalerate to 4'-hydroxyesfenvalerate efficiently except CYP6AE20. In contrast, none of the recombinant CYP6AE enzymes could metabolise gossypol under our experimental conditions. Epoxidation capabilities were observed in the CYP6AE subfamily, aldrin can be converted to dieldrin at rates up to 0.45 ± 0.04 pmol/min/pmol P450. Seven P450s in this subfamily can metabolise imidacloprid, but with lower efficiency than Bemisia tabaci CYP6CM1vQ. CYP6AE20 had virtually no metabolic competence to these four compounds but could metabolise several model fluorogenic substrates. These results showed the broad substrate spectrum of H. armigera CYP6AE P450s and suggest a limited role of gossypol upon the evolution of H. armigera CYP6AE genes.

Cytochrome P450 diversification and hostplant utilization patterns in specialist and generalist moths: Birth, death and adaptation

Across insect genomes, the size of the cytochrome P450 monooxygenase (CYP) gene superfamily varies widely. CYPome size variation has been attributed to reciprocal adaptive radiations in insect detoxification genes in response to plant biosynthetic gene radiations driven by co-evolution between herbivores and their chemically defended hostplants. Alternatively, variation in CYPome size may be due to random "birth-and-death" processes, whereby exponential increase via gene duplications is limited by random decay via gene death or transition via divergence. We examined CYPome diversification in the genomes of seven Lepidoptera species varying in host breadth from monophagous (Bombyx mori) to highly polyphagous (Amyelois transitella). CYPome size largely reflects the size of Clan 3, the clan associated with xenobiotic detoxification, and to some extent phylogenetic age. Consistently across genomes, families CYP6, CYP9 and CYP321 are most diverse and CYP6AB, CYP6AE, CYP6B, CYP9A and CYP9G are most diverse among subfamilies. Higher gene number in subfamilies is due to duplications occurring primarily after speciation and specialization ("P450 blooms"), and the genes are arranged in clusters, indicative of active duplicating loci. In the parsnip webworm, Depressaria pastinacella, gene expression levels in large subfamilies are high relative to smaller subfamilies. Functional and phylogenetic data suggest a correlation between highly dynamic loci (reflective of extensive gene duplication, functionalization and in some cases loss) and the ability of enzymes encoded by these genes to metabolize hostplant defences, consistent with an adaptive, nonrandom process driven by ecological interactions.

Two genomes of highly polyphagous lepidopteran pests (Spodoptera frugiperda, Noctuidae) with different host-plant ranges

Emergence of polyphagous herbivorous insects entails significant adaptation to recognize, detoxify and digest a variety of host-plants. Despite of its biological and practical importance - since insects eat 20% of crops - no exhaustive analysis of gene repertoires required for adaptations in generalist insect herbivores has previously been performed. The noctuid moth Spodoptera frugiperda ranks as one of the world’s worst agricultural pests. This insect is polyphagous while the majority of other lepidopteran herbivores are specialist. It consists of two morphologically indistinguishable strains (“C” and “R”) that have different host plant ranges. To describe the evolutionary mechanisms that both enable the emergence of polyphagous herbivory and lead to the shift in the host preference, we analyzed whole genome sequences from laboratory and natural populations of both strains. We observed huge expansions of genes associated with chemosensation and detoxification compared with specialist Lepidoptera. These expansions are largely due to tandem duplication, a possible adaptation mechanism enabling polyphagy. Individuals from natural C and R populations show significant genomic differentiation. We found signatures of positive selection in genes involved in chemoreception, detoxification and digestion, and copy number variation in the two latter gene families, suggesting an adaptive role for structural variation.

Segmental duplications: evolution and impact among the current Lepidoptera genomes

Background

Structural variation among genomes is now viewed to be as important as single nucleoid polymorphisms in influencing the phenotype and evolution of a species. Segmental duplication (SD) is defined as segments of DNA with homologous sequence.

Results

Here, we performed a systematic analysis of segmental duplications (SDs) among five lepidopteran reference genomes (Plutella xylostella, Danaus plexippus, Bombyx mori, Manduca sexta and Heliconius melpomene) to understand their potential impact on the evolution of these species. We find that the SDs content differed substantially among species, ranging from 1.2% of the genome in B. mori to 15.2% in H. melpomene. Most SDs formed very high identity (similarity higher than 90%) blocks but had very few large blocks. Comparative analysis showed that most of the SDs arose after the divergence of each linage and we found that P. xylostella and H. melpomene showed more duplications than other species, suggesting they might be able to tolerate extensive levels of variation in their genomes. Conserved ancestral and species specific SD events were assessed, revealing multiple examples of the gain, loss or maintenance of SDs over time. SDs content analysis showed that most of the genes embedded in SDs regions belonged to species-specific SDs (“Unique” SDs). Functional analysis of these genes suggested their potential roles in the lineage-specific evolution. SDs and flanking regions often contained transposable elements (TEs) and this association suggested some involvement in SDs formation. Further studies on comparison of gene expression level between SDs and non-SDs showed that the expression level of genes embedded in SDs was significantly lower, suggesting that structure changes in the genomes are involved in gene expression differences in species.

Conclusions

The results showed that most of the SDs were “unique SDs”, which originated after species formation. Functional analysis suggested that SDs might play different roles in different species. Our results provide a valuable resource beyond the genetic mutation to explore the genome structure for future Lepidoptera research.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1007-y) contains supplementary material, which is available to authorized users.

A reference gene set for sex pheromone biosynthesis and degradation genes from the diamondback moth, Plutella xylostella, based on genome and transcriptome digital gene expression analyses

Background

Female moths synthesize species-specific sex pheromone components and release them to attract male moths, which depend on precise sex pheromone chemosensory system to locate females. Two types of genes involved in the sex pheromone biosynthesis and degradation pathways play essential roles in this important moth behavior. To understand the function of genes in the sex pheromone pathway, this study investigated the genome-wide and digital gene expression of sex pheromone biosynthesis and degradation genes in various adult tissues in the diamondback moth (DBM), Plutella xylostella, which is a notorious vegetable pest worldwide.

Results

A massive transcriptome data (at least 39.04 Gb) was generated by sequencing 6 adult tissues including male antennae, female antennae, heads, legs, abdomen and female pheromone glands from DBM by using Illumina 4000 next-generation sequencing and mapping to a published DBM genome. Bioinformatics analysis yielded a total of 89,332 unigenes among which 87 transcripts were putatively related to seven gene families in the sex pheromone biosynthesis pathway. Among these, seven [two desaturases (DES), three fatty acyl-CoA reductases (FAR) one acetyltransferase (ACT) and one alcohol dehydrogenase (AD)] were mainly expressed in the pheromone glands with likely function in the three essential sex pheromone biosynthesis steps: desaturation, reduction, and esterification. We also identified 210 odorant-degradation related genes (including sex pheromone-degradation related genes) from seven major enzyme groups. Among these genes, 100 genes are new identified and two aldehyde oxidases (AOXs), one aldehyde dehydrogenase (ALDH), five carboxyl/cholinesterases (CCEs), five UDP-glycosyltransferases (UGTs), eight cytochrome P450 (CYP) and three glutathione S-transferases (GSTs) displayed more robust expression in the antennae, and thus are proposed to participate in the degradation of sex pheromone components and plant volatiles.

Conclusions

To date, this is the most comprehensive gene data set of sex pheromone biosynthesis and degradation enzyme related genes in DBM created by genome- and transcriptome-wide identification, characterization and expression profiling. Our findings provide a basis to better understand the function of genes with tissue enriched expression. The results also provide information on the genes involved in sex pheromone biosynthesis and degradation, and may be useful to identify potential gene targets for pest control strategies by disrupting the insect-insect communication using pheromone-based behavioral antagonists.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3592-y) contains supplementary material, which is available to authorized users.