Diversity and evolution of the P450 family in arthropods

Designing and applying a methodology to assess sperm cell viability and DNA damage in a model amphipod

Sperm quality is defined as the sperm cell ability to successfully fertilize eggs and allow normal embryo development⁠. Few studies explore sperm quality using aquatic invertebrates. Parhyale hawaiensis is a marine amphipod with a circumtropical distribution and considered a model for evolution, development, and ecotoxicological studies. We aimed to develop a methodology to collect sperm cells of P. hawaiensis and evaluate their viability and DNA damage (comet assay). We directly exposed the sperm cells to different mutagenic agents to optimize/develop the protocols. Then, as a proof of concept, we exposed the males to mutagenic compounds (EMS, benzo[a]pyrene (BaP), azo and anthraquinone dyes) at non-lethal concentrations verified by the proposed viability test and analyzed their sperm cells for DNA damage (comet assay). Organisms exposed to EMS presented a clear concentration response in the DNA damage response. We also showed that BaP was able to induce a statistically significant increase in DNA damage of the sperm cells. For the two dyes, although DNA damage increased, statistically differences were not observed. We believe we successfully developed a test to detect genotoxicity of chemicals in sperm cells using an invertebrate model. The protocol for sperm cell viability needs to be further explored with different chemicals to verify its utility as a toxicity endpoint. The developed genotoxicity test has the advantages to employ organisms that are easily cultivated in reduced space, use simple laboratory resources and reduced amount of material and reagents. Positive responses with this model could be used to disclose new germ cell mutagen candidates which could be further confirmed in vertebrates' systems.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Identification and functional validation of P450 genes associated with pyrethroid resistance in the malaria vector Anopheles sinensis (Diptera Culicidae)

Cytochrome P450 monooxygenases (P450s), a multifunctional protein superfamily, are one of three major classes of detoxification enzymes. However, the diversity and functions of P450 genes from pyrethroid-resistant populations of Anopheles sinensis have not been fully explored. In this study, P450 genes associated with pyrethroid resistance were systematically screened using RNA-seq in three field pyrethroid-resistant populations (AH-FR, CQ-FR, YN-FR) and one laboratory resistant strain (WX-LR) at developmental stages, tissues, and post blood-meal in comparison to the laboratory susceptible strain (WX-LS) in An. sinensis. Importantly, the expression of significantly upregulated P450s was verified using RT-qPCR, and the function of selected P450s in pyrethroid detoxification was determined with RNA interference using four laboratory pyrethroid-resistant strains (WX-LR, AH-LR, CQ-LR, YN-LR). Sixteen P450 genes were significantly upregulated in at least one field-resistant population, and 44 were significantly upregulated in different developmental stages, tissues or post blood-meal. A total of 19 P450s were selected to verify their association with pyrethroid resistance, and four of them (AsCYP6P3v1, AsCYP6P3v2, AsCYP9J10, and AsCYP9K1) demonstrated significant upregulation in laboratory pyrethroid-resistant strains using RT-qPCR. Knockdown of these four genes all significantly reduced pyrethroid resistance and increased the mortality by 57.19% (AsCYP6P3v1 and AsCYP6P3v2 knockdown group), 38.39% (AsCYP9K1 knockdown group) and 48.87% (AsCYP9J10 knockdown group) in An. sinensis by RNAi, which determined the pyrethroid detoxification function of these four genes. This study revealed the diversity of P450 genes and provided functional evidence for four P450s in pyrethroid detoxification in An. sinensis for the first time, which increases our understanding of the pyrethroid resistance mechanism, and is of potential value for pyrethroid resistance detection and surveillance.

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

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

Cytochrome P450 monooxygenase systems: Diversity and plasticity for adaptive stress response

Superfamily of cytochromes P450 (CYPs) is composed of heme-thiolate-containing monooxygenase enzymes, which play crucial roles in the biosynthesis, bioactivation, and detoxification of a variety of organic compounds, both endogenic and exogenic. Majority of CYP monooxygenase systems are multi-component and contain various redox partners, cofactors and auxiliary proteins, which contribute to their diversity in both prokaryotes and eukaryotes. Recent progress in bioinformatics and computational biology approaches make it possible to undertake whole-genome and phylogenetic analyses of CYPomes of a variety of organisms. Considerable variations in sequences within and between CYP families and high similarity in secondary and tertiary structures between all CYPs along with dramatic conformational changes in secondary structure elements of a substrate binding site during catalysis have been reported. This provides structural plasticity and substrate promiscuity, which underlie functional diversity of CYPs. Gene duplication and mutation events underlie CYP evolutionary diversity and emergence of novel selectable functions, which provide the involvement of CYPs in high adaptability to changing environmental conditions and dietary restrictions. In our review, we discuss the recent advancements and challenges in the elucidating the evolutionary origin and mechanisms underlying the CYP monooxygenase system diversity and plasticity. Our review is in the view of hypothesis that diversity of CYP monooxygenase systems is translated into the broad metabolic profiles, and this has been acquired during the long evolutionary time to provide structural plasticity leading to high adaptative capabilities to environmental stress conditions.

Comparative Genomics Uncovers the Evolutionary Dynamics of Detoxification and Insecticide Target Genes Across 11 Phlebotomine Sand Flies

Sand flies infect more than 1 million people annually with Leishmania parasites and other bacterial and viral pathogens. Progress in understanding sand fly adaptations to xenobiotics has been hampered by the limited availability of genomic resources. To address this gap, we sequenced, assembled, and annotated the transcriptomes of 11 phlebotomine sand fly species. Subsequently, we leveraged these genomic resources to generate novel evolutionary insights pertaining to their adaptations to xenobiotics, including those contributing to insecticide resistance. Specifically, we annotated over 2,700 sand fly detoxification genes and conducted large-scale phylogenetic comparisons to uncover the evolutionary dynamics of the five major detoxification gene families: cytochrome P450s (CYPs), glutathione-S-transferases (GSTs), UDP-glycosyltransferases (UGTs), carboxyl/cholinesterases (CCEs), and ATP-binding cassette (ABC) transporters. Using this comparative approach, we show that sand flies have evolved diverse CYP and GST gene repertoires, with notable lineage-specific expansions in gene groups evolutionarily related to known xenobiotic metabolizers. Furthermore, we show that sand flies have conserved orthologs of (i) CYP4G genes involved in cuticular hydrocarbon biosynthesis, (ii) ABCB genes involved in xenobiotic toxicity, and (iii) two primary insecticide targets, acetylcholinesterase-1 (Ace1) and voltage gated sodium channel (VGSC). The biological insights and genomic resources produced in this study provide a foundation for generating and testing hypotheses regarding the molecular mechanisms underlying sand fly adaptations to xenobiotics.

The cytochrome P450 subfamilies CYP392A and CYP392D are key players in acaricide metabolism in Tetranychusurticae

The spider mite Tetranychus urticae is a major agricultural pest with a global distribution, extremely diverse host range and a remarkable ability to develop resistance to a wide variety of acaricides. P450 mono-oxygenases have been frequently associated with resistance development in this species. In particular enzymes of the CYP392A-subfamily were shown to metabolize a number of key acaricides, including abamectin, amitraz, fenpyroximate and the active metabolite of pyflubumide. However, transcriptomic studies comparing highly resistant and susceptible populations have often revealed high expression of members of the CYP392D-subfamily, but these have been only poorly studied. Here, we conducted a meta-analysis of gene expression data of 20 populations and identified two key enzymes of this family, CYP392D2 and CYP392D8, whose expression is associated with resistance. We subsequently functionally expressed these enzymes, together with CYP392A11 and CYP392A16 as known metabolizers, and compared their potential to accept a wide diversity of acaricides as substrate. This study overall confirms previous discovered substrates for CYP392A11 and CYP392A16, but also reveals unreported metabolic activity towards new acaricides. These include carbaryl, chlorpyrifos and etoxazole for CYP392A16 and carbaryl, chlorpyrifos and NNI-0711-NH pyflubumide for CYP392A11. For the newly studied CYP392D-family, we show that CYP392D2 metabolizes pyridaben, fenpyroximate, etoxazole and chlorpyrifos, while CYP392D8 metabolizes carbaryl, fenazaquin and tebufenpyrad. Last, we observed that both CYP392A- and CYP392D-subfamily enzymes activate chlorpyrifos to its corresponding oxon. Our study indicates that there is both overlap and specificity in the activity of A- and D-subfamily enzymes against acaricides and model substrates. With the recent advent of highly efficient CRISPR/Cas9 gene editing protocols in T. urticae, the way is now paved to conduct further genetic experiments revealing and quantifying the role of these enzymes in the resistance phenotype in field populations.

Identification and expression dynamics of CYPome across different developmental stages of Maconellicoccus hirsutus (Green)

Maconellicoccus hirsutus is a highly polyphagous insect pest, posing a substantial threat to various crop sp., especially in the tropical and sub-tropical regions of the world. While extensive physiological and biological studies have been conducted on this pest, the lack of genetic information has hindered our understanding of the molecular mechanisms underlying its growth, development, and xenobiotic metabolism. The Cytochrome P450 gene, a member of the CYP gene superfamily ubiquitous in living organisms is associated with growth, development, and the metabolism of both endogenous and exogenous substances, contributing to the insect's adaptability in diverse environments. To elucidate the specific role of the CYP450 gene family in M. hirsutus which has remained largely unexplored, a de novo transcriptome assembly of the pink mealybug was constructed. A total of 120 proteins were annotated as CYP450 genes through homology search of the predicted protein sequences across different databases. Phylogenetic studies resulted in categorizing 120 CYP450 genes into four CYP clans. A total of 22 CYP450 families and 30 subfamilies were categorized, with CYP6 forming the dominant family. The study also revealed five genes (Halloween genes) associated with the insect hormone biosynthesis pathway. Further, the expression of ten selected CYP450 genes was studied using qRT-PCR across crawler, nymph, and adult stages, and identified genes that were expressed at specific stages of the insects. Thus, the findings of this study reveal the expression dynamics and possible function of the CYP450 gene family in the growth, development, and adaptive strategies of M. hirsutus which can be further functionally validated.

Differentiation in detoxification gene complements, including neofunctionalization of duplicated cytochrome P450 genes, between lineages of cotton bollworm, Helicoverpa armigera

Here we investigate the evolutionary dynamics of five enzyme superfamilies (CYPs, GSTs, UGTs, CCEs and ABCs) involved in detoxification in Helicoverpa armigera. The reference assembly for an African isolate of the major lineages, H. a. armigera, has 373 genes in the five superfamilies. Most of its CYPs, GSTs, UGTs and CCEs and a few of its ABCs occur in blocks and most of the clustered genes are in subfamilies specifically implicated in detoxification. Most of the genes have orthologues in the reference genome for the Oceania lineage, H. a. conferta. However, clustered orthologues and subfamilies specifically implicated in detoxification show greater sequence divergence and less constraint on non-synonymous differences between the two assemblies than do other members of the five superfamilies. Two duplicated CYPs, which were found in the H. a. armigera but not H. a. conferta reference genome, were also missing in 16 Chinese populations spanning two different lineages of H. a. armigera. The enzyme produced by one of these duplicates has higher activity against esfenvalerate than a previously described chimeric CYP mutant conferring pyrethroid resistance. Various transposable elements were found in the introns of most detoxification genes, generating diverse gene structures. Extensive resequencing data for the Chinese H. a. armigera and H. a. conferta lineages also revealed complex copy number polymorphisms in 17 CCE001s in a cluster also implicated in pyrethroid metabolism, with substantial haplotype differences between all three lineages. Our results suggest that cotton bollworm has a versatile complement of detoxification genes which are evolving in diverse ways across its range.

Exploratory comparative transcriptomic analysis reveals potential gene targets associated with Cry1A.105 and Cry2Ab2 resistance in fall armyworm (Spodoptera frugiperda)

Genetically modified (GM) crops, expressing Bacillus thuringiensis (Bt) insecticidal toxins, have substantially transformed agriculture. Despite rapid adoption, their environmental and economic benefits face scrutiny due to unsustainable agricultural practices and the emergence of resistant pests like Spodoptera frugiperda, known as the fall armyworm (FAW). FAW's adaptation to Bt technology in corn and cotton compromises the long-term efficacy of Bt crops. To advance the understanding of the genetic foundations of resistance mechanisms, we conducted an exploratory comparative transcriptomic analysis of two divergent FAW populations. One population exhibited practical resistance to the Bt insecticidal proteins Cry1A.105 and Cry2Ab2, expressed in the genetically engineered MON-89Ø34 - 3 maize, while the other population remained susceptible to these proteins. Differential expression analysis supported that Cry1A.105 and Cry2Ab2 significantly affect the FAW physiology. A total of 247 and 254 differentially expressed genes were identified in the Cry-resistant and susceptible populations, respectively. By integrating our findings with established literature and databases, we underscored 53 gene targets potentially involved in FAW's resistance to Cry1A.105 and Cry2Ab2. In particular, we considered and discussed the potential roles of the differentially expressed genes encoding ABC transporters, G protein-coupled receptors, the P450 enzymatic system, and other Bt-related detoxification genes. Based on these findings, we emphasize the importance of exploratory transcriptomic analyses to uncover potential gene targets involved with Bt insecticidal proteins resistance, and to support the advantages of GM crops in the face of emerging challenges.

The moulting arthropod: a complete genetic toolkit review

Exoskeletons are a defining character of all arthropods that provide physical support for their segmented bodies and appendages as well as protection from the environment and predation. This ubiquitous yet evolutionarily variable feature has been instrumental in facilitating the adoption of a variety of lifestyles and the exploitation of ecological niches across all environments. Throughout the radiation that produced the more than one million described modern species, adaptability afforded by segmentation and exoskeletons has led to a diversity that is unrivalled amongst animals. However, because of the limited extensibility of exoskeleton chitin and cuticle components, they must be periodically shed and replaced with new larger ones, notably to accommodate the growing individuals encased within. Therefore, arthropods grow discontinuously by undergoing periodic moulting events, which follow a series of steps from the preparatory pre-moult phase to ecdysis itself and post-moult maturation of new exoskeletons. Each event represents a particularly vulnerable period in an arthropod's life cycle, so processes must be tightly regulated and meticulously executed to ensure successful transitions for normal growth and development. Decades of research in representative arthropods provide a foundation of understanding of the mechanisms involved. Building on this, studies continue to develop and test hypotheses on the presence and function of molecular components, including neuropeptides, hormones, and receptors, as well as the so-called early, late, and fate genes, across arthropod diversity. Here, we review the literature to develop a comprehensive overview of the status of accumulated knowledge of the genetic toolkit governing arthropod moulting. From biosynthesis and regulation of ecdysteroid and sesquiterpenoid hormones, to factors involved in hormonal stimulation responses and exoskeleton remodelling, we identify commonalities and differences, as well as highlighting major knowledge gaps, across arthropod groups. We examine the available evidence supporting current models of how components operate together to prepare for, execute, and recover from ecdysis, comparing reports from Chelicerata, Myriapoda, Crustacea, and Hexapoda. Evidence is generally highly taxonomically imbalanced, with most reports based on insect study systems. Biases are also evident in research on different moulting phases and processes, with the early triggers and late effectors generally being the least well explored. Our synthesis contrasts knowledge based on reported observations with reasonably plausible assumptions given current taxonomic sampling, and exposes weak assumptions or major gaps that need addressing. Encouragingly, advances in genomics are driving a diversification of tractable study systems by facilitating the cataloguing of putative genetic toolkits in previously under-explored taxa. Analysis of genome and transcriptome data supported by experimental investigations have validated the presence of an "ultra-conserved" core of arthropod genes involved in moulting processes. The molecular machinery has likely evolved with elaborations on this conserved pathway backbone, but more taxonomic exploration is needed to characterise lineage-specific changes and novelties. Furthermore, linking these to transformative innovations in moulting processes across Arthropoda remains hampered by knowledge gaps and hypotheses based on untested assumptions. Promisingly however, emerging from the synthesis is a framework that highlights research avenues from the underlying genetics to the dynamic molecular biology through to the complex physiology of moulting.

Toxicity of cockroach gel baits to the oothecal parasitoid Aprostocetus hagenowii (Hymenoptera: Eulophidae) and implications for cockroach integrated pest management

Aprostocetus hagenowii (Ratzeburg) is a parasitoid wasp that parasitizes the oothecae of peridomestic pest cockroaches. A. hagenowii has been used in integrated pest management (IPM) programs for cockroach control but little is known about how this parasitoid responds to the insecticides commonly used for cockroach management. Five insecticidal gel bait products containing indoxacarb, clothianidin, fipronil, dinotefuran, or abamectin B1 were tested for their toxicity towards A. hagenowii and the American cockroach, Periplaneta americana (L.; Blattodea: Blattidae), a host of A. hagenowii and a common pest. All baits were tested as fresh and 1-d aged deposits. Indoxacarb was the only active ingredient that did not cause significant (P < 0.05) A. hagenowii mortality compared to the control in both the fresh and aged gel experiments (Median survival time [MST]s: 168 h fresh, 72 h aged). Clothianidin caused the lowest A. hagenowii MSTs across experiments (24 h, fresh and aged). All baits caused significant P. americana mortality as fresh and 1-d aged deposits (P < 0.05). Indoxacarb appears most compatible with A. hagenowii in cockroach IPM.

The natural anthraquinone dye emodin: Eco/genotoxicological characterization for aquatic organisms

Emodin is an anthraquinone secondary metabolite produced by several species of plants and fungi. Emodin is known for its pharmacological versatility, and, in the textile industry, for its good dyeing properties. However, its use in the textile industry can result in the formation and disposal of large volumes of wastewater. Emodin mutagenicity has been shown in bacteria and in human cells, but little is known about its possible toxic, genotoxic, or mutagenic effects in aquatic organisms. We have evaluated the eco/genotoxicity of emodin to aquatic organisms. Emodin was toxic to Daphnia similis (EC50 = 130 μg L-1) and zebrafish embryos (LC50 = 25 μg L-1). No toxicity was observed for Raphidocelis subcapitata, Ceriodaphnia dubia, or Parhyale hawaiensis. Additional biochemistry/molecular studies are needed to elucidate the toxic/mutagenic pathways of emodin in aquatic organisms. The PNEC value for emodin was 0.025 μg L-1. In addition to mutagenicity in the Salmonella/microsome assay, emodin was mutagenic in the micronucleus assay in the amphipod P. hawaiensis. Among the anthraquinone dyes tested to date, natural or synthetic, emodin was the most toxic to aquatic species.

Trial and error: New insights into recombinant expression of membrane-bound insect cytochromes P450 in Escherichia coli systems

Insect cytochromes P450 (CYP450s) are key enzymes responsible for a wide array of oxidative transformations of both endogenous and exogenous substrates. However, there is currently no a universal guideline established for heterologous expression of membrane-bound CYP450s, which hampers their downstream biochemical and structural studies. In this study, we conducted large-scale screening of protein overexpression in Escherichia coli using 71 insect CYP450 sequences and optimized the expression of a difficult-to-express CYP450 (CYP6HX3) using eight different optimizations, including selection of host strains and expression vectors, alternative of leader signal peptides, and N-terminal modifications. We confirmed that 1) Only insect CYP450s belonging to the CYP347 family could be expressed with N-terminal fusion of ompA2+ signal peptide in E. coli expression system. 2) E. coli Lemo 21 (DE3) effectively improved the expression of CYP6HX3 in the plasma membrane. 3) A brick-red appearance occurred frequently in the expressed thallus or membrane proteins, but this phenomenon could not necessarily indicate successful overexpression of target CYP450s. These findings provide new insights into the recombinant expression of insect CYP450s in E. coli systems and will facilitate the theoretical approaches for functional expression and production of eukaryotic CYP450s.

Role of CYP311A1 in wing development of Drosophila melanogaster

Lipid homeostasis is crucial for growth and development of organisms. Several cytochrome P450 monooxygenases (CYPs) are involved in lipid metabolism. The function of Cyp311a1 in the anterior midgut as a regulator of phosphatidylethanolamine (PE) metabolism in Drosophila melanogaster has been demonstrated, as depletion of Cyp311a1 caused larval growth arrest that was partially rescued by supplying PE. In this study, we investigated the role of CYP311A1 in wing morphogenesis in Drosophila. Using the GAL4-UAS system, Cyp311a1 was selectively knocked down in the wing disc. A deformed wing phenotype was observed in flies with reduced Cyp311a1 transcripts. BODIPY and oil red O staining revealed a reduction of neutral lipids in the wing disc after the depletion of Cyp311a1. In addition, we observed an enhanced sensitivity to Eosin Y penetration in the wings of Cyp311a1 knocked-down flies. Moreover, the reduction of CYP311A1 function in developing wings does not affect cell proliferation and apoptosis, but entails disordered Phalloidin or Cadherin distribution, suggesting an abnormal cell morphology and cell cortex structure in wing epithelial cells. Taken together, our results suggest that Cyp311a1 is needed for wing morphogenesis by participating in lipid assembly and cell homeostasis.

Genome sequencing, comparative analysis, and gene expression responses of cytochrome P450 genes in Oryzias curvinotus provide insights into environmental adaptation

The mangrove fish (Oryzias curvinotus) serves as a model for researching environmental adaptation and sexual development. To further such research, we sequenced and assembled a high-quality 842 Mb reference genome for O. curvinotus. Comparative genomic analysis revealed 891 expanded gene families, including significantly expanded cytochrome P450 (CYP) detoxification genes known to be involved in xenobiotic defense. We identified 69 O. curvinotus CYPs (OcuCYPs) across 18 families and 10 clans using multiple methods. Extensive RNA-seq and qPCR analysis demonstrated diverse spatiotemporal expression patterns of OcuCYPs by developmental stage, tissue type, sex, and pollutant exposure (17β-estradiol (E2) and testosterone (MT)). Many OcuCYPs exhibited sexual dimorphism in gonads, suggesting reproductive roles in steroidogenesis, while their responsiveness to model toxicants indicates their importance in environmental adaptation through enhanced detoxification. Pathway analysis highlighted expanded CYP genes in arachidonic acid metabolism, drug metabolism, and steroid hormone biosynthesis. This chromosome-level genomic resource provides crucial biological insights to elucidate the functional roles of expanded CYPs in environmental adaptation, sexual development, early life history, and conservation in the anthropogenically impacted mangrove habitats of O. curvinotus. It also enables future ecotoxicology research leveraging O. curvinotus as a pollution sentinel species.

CYP4G subfamily genes mediate larval integument development in Spodoptera frugiperda

Cytochrome P450 (CYP) 4G subfamily is closely related to the synthesis of cuticular hydrocarbons, leading to the enhanced desiccation and insecticide resistance of pests. However, functions of CYP4Gs in larval integument development remain unknown in Spodoptera frugiperda (J. E. Smith) (Lepidoptera: Noctuidae), which is a major transboundary migratory pest and become a common pest in China. On the basis of the genome and transcriptome datasets of S. frugiperda, CYP4G74, CYP4G75, CYP4G108, and CYP4G109 were identified, which contained the conserved domains of P450s and CYP4Gs. The spatial and temporal expression analysis showed that CYP4G74 and CYP4G75 were significantly highly expressed in adults and larval integuments, while CYP4G108 and CYP4G109 had low expressions in larval integuments. After silencing CYP4G74 and CYP4G75 by RNA interference, abnormal integument development occurred in larvae, some of which became smaller and dead, indicating important roles of CYP4G74 and CYP4G75 in the synthesis and development of integuments. The results clarify the functions of CYP4Gs in S. frugiperda and provide potential targets for the control of this pest.

Cyp6g2 is the major P450 epoxidase responsible for juvenile hormone biosynthesis in Drosophila melanogaster

BACKGROUND

Juvenile hormones (JH) play crucial role in regulating development and reproduction in insects. The most common form of JH is JH III, derived from MF through epoxidation by CYP15 enzymes. However, in the higher dipterans, such as the fruitfly, Drosophila melanogaster, a bis-epoxide form of JHB3, accounted most of the JH detected. Moreover, these higher dipterans have lost the CYP15 gene from their genomes. As a result, the identity of the P450 epoxidase in the JH biosynthesis pathway in higher dipterans remains unknown.

RESULTS

In this study, we show that Cyp6g2 serves as the major JH epoxidase responsible for the biosynthesis of JHB3 and JH III in D. melanogaster. The Cyp6g2 is predominantly expressed in the corpus allatum (CA), concurring with the expression pattern of jhamt, another well-studied gene that is crucial in the last steps of JH biosynthesis. Mutation in Cyp6g2 leads to severe disruptions in larval-pupal metamorphosis and exhibits reproductive deficiencies, exceeding those seen in jhamt mutants. Notably, Cyp6g2-/-::jhamt2 double mutants all died at the pupal stage but could be rescued through the topical application of JH analogs. JH titer analyses revealed that both Cyp6g2-/- mutant and jhamt2 mutant lacking JHB3 and JH III, while overexpression of Cyp6g2 or jhamt caused a significant increase in JHB3 and JH III titer.

CONCLUSIONS

These findings collectively established that Cyp6g2 as the major JH epoxidase in the higher dipterans and laid the groundwork for the further understanding of JH biosynthesis. Moreover, these findings pave the way for developing specific Cyp6g2 inhibitors as insect growth regulators or insecticides.

Cytochrome P450 enzymes in the black-spotted frog (Pelophylax nigromaculatus): molecular characterization and upregulation of expression by sulfamethoxazole

Cytochrome P450 (CYP) enzymes are crucial for the detoxification of xenobiotics, cellular metabolism, and homeostasis. This study investigated the molecular characterization of CYP enzymes in the black-spotted frog, Pelophylax nigromaculatus, and examined the regulation of CYP expression in response to chronic exposure to the antibiotic sulfamethoxazole (SMX) at various environmental concentrations (0, 1, 10, and 100 μg/L). The full-length cDNA of Pn-CYP26B1 was identified. The sequence included open reading frames of 1,536 bp, encoding proteins comprising 511 amino acids. The signature motif, FxxGxxxCxG, was highly conserved when compared with a number of selected animal species. SMX significantly upregulated the expression of the protein CYP26B1 in frog livers at concentrations of 1 and 10 μg/L. SMX showed an affinity for CYP26B1 of -7.6 kcal/mol, indicating a potential mechanism for SMX detoxification or adaptation of the frog. These findings contributed to our understanding of the environmental impact of antibiotics on amphibian species and underscored the importance of CYP enzymes in maintaining biochemical homeostasis under exposure to xenobiotic stress.

UDP-glycosyltransferases act as key determinants of host plant range in generalist and specialist Spodoptera species

Phytophagous insects have evolved sophisticated detoxification systems to overcome the antiherbivore chemical defenses produced by many plants. However, how these biotransformation systems differ in generalist and specialist insect species and their role in determining insect host plant range remains an open question. Here, we show that UDP-glucosyltransferases (UGTs) play a key role in determining the host range of insect species within the Spodoptera genus. Comparative genomic analyses of Spodoptera species that differ in host plant breadth identified a relatively conserved number of UGT genes in generalist species but high levels of UGT gene pseudogenization in the specialist Spodoptera picta. CRISPR-Cas9 knockouts of the three main UGT gene clusters of Spodoptera frugiperda revealed that UGT33 genes play an important role in allowing this species to utilize the poaceous plants maize, wheat, and rice, while UGT40 genes facilitate utilization of cotton. Further functional analyses in vivo and in vitro identified the UGT SfUGT33F32 as the key mechanism that allows generalist S. frugiperda to detoxify the benzoxazinoid DIMBOA (2,4-dihydroxy-7-methoxy-2H-1,4-benzoxazin-3(4H)-one), a potent insecticidal phytotoxin produced by poaceous plants. However, while this detoxification capacity is conserved in several generalist Spodoptera species, Spodoptera picta, which specializes on Crinum plants, is unable to detoxify DIMBOA due to a nonfunctionalizing mutation in SpUGT33F34. Collectively, these findings provide insight into the role of insect UGTs in host plant adaptation, the mechanistic basis of evolutionary transitions between generalism and specialism and offer molecular targets for controlling a group of notorious insect pests.

Preliminary Study on the Pathogenic Mechanism of Jujube Flower Disease in Honeybees (Apis mellifera ligustica) Based on Midgut Transcriptomics

Honeybees are prone to poisoning, also known as jujube flower disease, after collecting nectar from jujube flowers, resulting in the tumultuous demise of foragers. The prevalence of jujube flower disease has become one of the main factors affecting the development of the jujube and beekeeping industries in Northern China. However, the pathogenic mechanisms underlying jujube flower disease in honeybees are poorly understood. Herein, we first conducted morphological observations of the midgut using HE-staining and found that jujube flower disease-affected honeybees displayed midgut damage with peritrophic membrane detachment. Jujube flower disease was found to increase the activity of chitinase and carboxylesterase (CarE) and decrease the activity of superoxide dismutase (SOD), catalase (CAT), glutathione S-transferase (GST), and the content of CYP450 in the honeybee midgut. Transcriptomic data identified 119 differentially expressed genes in the midgut of diseased and healthy honeybees, including CYP6a13, CYP6a17, CYP304a1, CYP6a14, AADC, and AGXT2, which are associated with oxidoreductase activity and vitamin binding. In summary, collecting jujube flower nectar could reduce antioxidant and detoxification capacities of the honeybee midgut and, in more severe cases, damage the intestinal structure, suggesting that intestinal damage might be the main cause of honeybee death due to jujube nectar. This study provides new insights into the pathogenesis of jujube flower disease in honeybees.

Overexpression of CYP6CX4 contributing to field-evolved resistance to flupyradifurone, one novel butenolide insecticide, in Bemisia tabaci from China

Bemisia tabaci is a formidable insect pest worldwide, and exhibits significant resistance to various insecticides. Flupyradifurone is one novel butenolide insecticide and has emerged as a new weapon against B. tabaci, but field-evolved resistance to this insecticide has become a widespread concern. To unravel the mechanisms of field-evolved flupyradifurone resistance, we conducted a comprehensive investigation into susceptibility of twenty-one field populations within the Beijing-Tianjin-Hebei Region of China. Alarmingly, thirteen of these populations displayed varying degrees of resistance, ranging from low to medium levels, and building upon our prior findings, we meticulously cloned and characterized the CYP6CX4 gene in B. tabaci. Our investigations unequivocally confirmed the association between CYP6CX4 overexpression and flupyradifurone resistance in three of the thirteen resistant strains via RNA interference. To further validate our findings, we introduced CYP6CX4 overexpression into a transgenic Drosophila melanogaster line, resulting in a significant development of resistance to flupyradifurone in D. melanogaster. Additionally, homology modeling and molecular docking analyses showed the stable binding of flupyradifurone to CYP6CX4, with binding free energy of -6.72 kcal mol-1. Collectively, our findings indicate that the induction of CYP6CX4 exerts one important role in detoxification of flupyradifurone, thereby promoting development of resistance in B. tabaci.

De novo assembly and annotation of Popillia japonica's genome with initial clues to its potential as an invasive pest

BACKGROUND

The spread of Popillia japonica in non-native areas (USA, Canada, the Azores islands, Italy and Switzerland) poses a significant threat to agriculture and horticulture, as well as to endemic floral biodiversity, entailing that appropriate control measures must be taken to reduce its density and limit its further spread. In this context, the availability of a high quality genomic sequence for the species is liable to foster basic research on the ecology and evolution of the species, as well as on possible biotechnologically-oriented and genetically-informed control measures.

RESULTS

The genomic sequence presented and described here is an improvement with respect to the available draft sequence in terms of completeness and contiguity, and includes structural and functional annotations. A comparative analysis of gene families of interest, related to the species ecology and potential for polyphagy and adaptability, revealed a contraction of gustatory receptor genes and a paralogous expansion of some subgroups/subfamilies of odorant receptors, ionotropic receptors and cytochrome P450s.

CONCLUSIONS

The new genomic sequence as well as the comparative analyses data may provide a clue to explain the staggering invasive potential of the species and may serve to identify targets for potential biotechnological applications aimed at its control.

Analysis of insecticide resistance and de novo transcriptome assembly of resistance associated genes in the European grapevine moth, Lobesia botrana (Lepidoptera: Tortricidae)

The European grapevine moth Lobesia botrana (Denis & Shiffermüller 1776) is an economically important pest of the vine-growing areas worldwide. Chemical insecticides have been used for its control; however, its resistance status is largely unknown in many regions. We monitored the susceptibility of several L. botrana populations from Greece and Turkey. In addition, based on RNAseq transcriptome analysis, we identified and phylogenetically classify the cytochrome P450 genes of L. botrana, as well as analysed target site sequences and looked for the presence of known resistance mutations. Resistance against chlorantraniliprole, alpha-cypermethrin, spinetoram, etofenprox, and acetamiprid was very low (below 2.5-fold in all cases, compared to a reference strain from Greece) in all populations from Greece that were included in the study. However, resistance against indoxacarb (4-30-fold), spinosad (5-59-fold), and deltamethrin (18-30 fold) was detected in the L. botrana populations from Turkey, compared to a reference population from Turkey. De novo transcriptome assembly and manual annotation, and subsequent PCR-based analysis of insecticide target sequences (i.e. voltage-gated sodium channel - VGSC: target of pyrethroids and oxadiazines; nicotinic acetylcholine receptor subunit a6 - nAChR_α6: target of spinosad; ryanodine receptor - RyR: target of diamides; glutamate-gated chloride channel - GluCl: target of avermectins and; acetylcholinesterase - AChE: target of organophosphates) showed the absence of known resistance mutations in all specimens from both countries. Finally, the L. botrana CYPome (116 genes) was manually analysed and phylogenetically characterised, to provide resources for future studies that will aim the analysis of metabolic resistance.

Insights into the adaptive evolution of chromosome and essential traits through chromosome-level genome assembly of Gekko japonicus

Gekko japonicus possesses flexible climbing and detoxification abilities under insectivorous habits. Still, the evolutionary mechanisms behind these traits remain unclarified. This study presents a chromosome-level G. japonicus genome, revealing that its evolutionary breakpoint regions were enriched with specific repetitive elements and defense response genes. Gene families unique to G. japonicus and positively selected genes are mainly enriched in immune, sensory, and nervous pathways. Expansion of bitter taste receptor type 2 primarily in insectivorous species could be associated with toxin clearance. Detox cytochrome P450 in G. japonicus has undergone more birth and death processes than biosynthesis-type P450 genes. Proline, cysteine, glycine, and serine in corneous beta proteins of G. japonicus might influence flexibility and setae adhesiveness. Certain thermosensitive transient receptor potential channels under relaxed purifying selection or positive selection in G. japonicus might enhance adaptation to climate change. This genome assembly offers insights into the adaptive evolution of gekkotans.

Ontogeny of the Cytochrome P450 Superfamily in the Ornate Spiny Lobster (Panulirus ornatus)

Cytochrome P450s (CYP450s) are a versatile superfamily of enzymes known to undergo rapid evolution. They have important roles across growth and development pathways in crustaceans, although it is difficult to characterise orthologs between species due to their sequence diversity. Conserved CYP450s enzymes in crustaceans are those associated with ecdysteroidogenesis: synthesising and breaking down the active moult hormone, 20-hydroxyecdysone. The complex life cycle of the ornate spiny lobster, Panulirus ornatus, relies on moulting in order to grow and develop. Many of these diverse life stages have been analysed to establish a comprehensive transcriptomic database for this species. The transcripts putatively encoding for CYP450s were mapped using transcriptomic analysis and identified across growth and development stages. With the aid of phylogeny, 28 transcripts of 42 putative P. ornatus CYP450s were annotated, including the well conserved Halloween genes, which are involved in ecdysteroidogenesis. Expression patterns across the life stages determined that only a subset of the CYP450s can be detected in each life stage or tissue. Four Shed transcripts show overlapping expression between metamorphosis and adult tissues, suggesting pleotropic functions of the multiple Shed orthologs within P. ornatus.

A chromosome-level genome assembly of the Rhus gall aphid Schlechtendalia chinensis provides insight into the endogenization of Parvovirus-like DNA sequences

The Rhus gall aphid, Schlechtendalia chinensis, feeds on its primary host plant Rhus chinensis to induce galls, which have economic importance in medicines and the food industry. Rhus gall aphids have a unique life cycle and are economically beneficial but there is huge gap in genomic information about this group of aphids. Schlechtendalia chinensis induces rich-tannin galls on its host plant and is emerging as a model organism for both commercial applications and applied research in the context of gall production by insects. Here, we generated a high-quality chromosome-level assembly for the S. chinensis genome, enabling the comparison between S. chinensis and non-galling aphids. The final genome assembly is 344.59 Mb with 91.71% of the assembled sequences anchored into 13 chromosomes. We predicted 15,013 genes, of which 14,582 (97.13%) coding genes were annotated, and 99% of the predicted genes were anchored to the 13 chromosomes. This assembly reveals the endogenization of parvovirus-related DNA sequences (PRDs) in the S. chinensis genome, which could play a role in environmental adaptations. We demonstrated the characterization and classification of cytochrome P450s in the genome assembly, which are functionally crucial for sap-feeding insects and have roles in detoxification and insecticide resistance. This genome assembly also revealed the whole genome duplication events in S. chinensis, which can be considered in comparative evolutionary analysis. Our work represents a reference genome for gall-forming aphids that could be used for comparative genomic studies between galling and non-galling aphids and provides the first insight into the endogenization of PRDs in the genome of galling aphids. It also provides novel genetic information for future research on gall-formation and insect-plant interactions.

Contrasting roles of cytochrome P450s in amitraz and chlorfenapyr resistance in the crop pest Tetranychus urticae

The molecular mechanisms of amitraz and chlorfenapyr resistance remain only poorly understood for major agricultural pests and vectors of human diseases. This study focusses on a multi-resistant field strain of the crop pest Tetranychus urticae, which could be readily selected in the laboratory to high levels of amitraz and chlorfenapyr resistance. Toxicity experiments using tralopyril, the active toxophore of chlorfenapyr, suggested decreased activation as a likely mechanism underlying resistance. Starting from the same parental strain, transcriptome profiling revealed that a cluster of detoxifying genes was upregulated after amitraz selection, but unexpectedly downregulated after chlorfenapyr selection. Further functional validation associated the upregulation of CYP392A16 with amitraz metabolism and the downregulation of CYP392D8 with reduced activation of chlorfenapyr to tralopyril. Genetic mapping (QTL analysis by BSA) was conducted in an attempt to unravel the genetic mechanisms of expression variation and resistance. This revealed that chlorfenapyr resistance was associated with a single QTL, while 3 QTLs were uncovered for amitraz resistance. Together with the observed contrasting gene expression patterns, we argue that transcriptional regulators most likely underly the distinct expression profiles associated with resistance, but these await further functional validation.

Identification of inducible CYP3 and CYP4 genes associated with abamectin tolerance in the fat body and Malpighian tubules of Spodoptera litura

Abamectin, as a broad-spectrum bioinsecticide, has been widely used for the control of Lepidoptera insects, resulting in different levels of resistance to abamectin in Spodoptera litura. Cytochrome P450 monooxygenases (P450s) are known for their important roles in insecticide detoxification. In this study, the expression of SlCYP6B40, SlCYP4L12 and SlCYP9A32 in the fat body, and SlCYP4S9, SlCYP6AB12, SlCYP6AB58, SlCYP9A75a and SlCYP9A75b in Malpighian tubules was found to be significantly upregulated after abamectin exposure. SlCYP6AE44 and SlCYP6AN4 were simultaneously upregulated in these two tissues after abamectin exposure. Ectopically overexpressed SlCYP6AE44, SlCYP9A32 and SlCYP4S9 in transgenic Drosophila conferred tolerance to abamectin. In addition, homology modeling and molecular docking results suggested that SlCYP6AE44, SlCYP9A32 and SlCYP4S9 may be capable of binding with abamectin. These results demonstrate that upregulation of CYP3 and CYP4 genes may contribute to abamectin detoxification in S. litura and provide information for evidence-based insecticide resistance management strategies.

Insights into unique features of Drosophila CYP4G enzymes

The cytochrome P450 enzymes of the CYP4G subfamily are some of the most intriguing insect P450s in terms of structure and function. In Drosophila, CYP4G1 is highly expressed in the oenocytes and is the last enzyme in the biosynthesis of cuticular hydrocarbons, while CYP4G15 is expressed in the brain and is of unknown function. Both proteins have a CYP4G-specific and characteristic amino acid sequence insertion corresponding to a loop between the G and H helices whose function is unclear. Here we address these enigmatic structural and functional features of Drosophila CYP4Gs. First, we used reverse genetics to generate D. melanogaster strains in which all or part of the CYP4G-specific loop was removed from CYP4G1. We showed that the full loop was not needed for proper folding of the P450, but it is essential for function, and that just a short stretch of six amino acids is required for the enzyme's ability to make hydrocarbons. Second, we confirmed by immunocytochemistry that CYP4G15 is expressed in the brain and showed that it is specifically associated with the cortex glia cell subtype. We then expressed CYP4G15 ectopically in oenocytes, revealing that it can produce of a blend of hydrocarbons, albeit to quantitatively lower levels resulting in only a partial rescue of CYP4G1 knockdown flies. The CYP4G1 structural variants studied here should facilitate the biochemical characterization of CYP4G enzymes. Our results also raise the question of the putative role of hydrocarbons and their synthesis by cortex glial cells.

A life-history allele of large effect shortens developmental time in a wild insect population

Developmental time is a key life-history trait with large effects on Darwinian fitness. In many insects, developmental time is currently under strong selection to minimize ecological mismatches in seasonal timing induced by climate change. The genetic basis of responses to such selection, however, is poorly understood. To address this problem, we set up a long-term evolve-and-resequence experiment in the beetle Tribolium castaneum and selected replicate, outbred populations for fast or slow embryonic development. The response to this selection was substantial and embryonic developmental timing of the selection lines started to diverge during dorsal closure. Pooled whole-genome resequencing, gene expression analysis and an RNAi screen pinpoint a 222 bp deletion containing binding sites for Broad and Tramtrack upstream of the ecdysone degrading enzyme Cyp18a1 as a main target of selection. Using CRISPR/Cas9 to reconstruct this allele in the homogenous genetic background of a laboratory strain, we unravel how this single deletion advances the embryonic ecdysone peak inducing dorsal closure and show that this allele accelerates larval development but causes a trade-off with fecundity. Our study uncovers a life-history allele of large effect and reveals the evolvability of developmental time in a natural insect population.

Functional orthologs of honeybee CYP6AQ1 in stingless bees degrade the butenolide insecticide flupyradifurone

Flupyradifurone (FPF), a novel butenolide insecticide binding to nicotinic acetylcholine receptors (nAChRs), has been shown to be less acutely toxic to western honey bees (Apis mellifera) than other insecticides such as neonicotinoids sharing the same target-site. A previous study revealed that this is due to enhanced oxidative metabolism of FPF, mediated by three cytochrome P450 monooxygenases (P450s), including CYP6AQ1. Therefore, we followed a toxicogenomics approach and investigated the potential role of functional CYP6AQ1 orthologs in FPF metabolism from eight different bee species, including stingless bees (Tribe: Meliponini). We conducted a phylogenetic analysis on four stingless bee species, including Frieseomelitta varia, Heterotrigona itama, Melipona quadrifasciata and Tetragonula carbonaria to identify CYP6AQ1-like functional orthologs. Three non-Meliponini, but tropical bee species, i.e., Ammobates syriacus, Euglossa dilemma and Megalopta genalis were analyzed as well. We identified candidate P450s in all (neo)tropical species with greater than 61% and 67% predicted protein sequence identities when compared to A. mellifera CYP6AQ1 and Bombus terrestris CYP6AQ26, respectively. Heterologous expression in High Five insect cells of these functional orthologs revealed a common coumarin substrate profile and a preference for the O-debenzylation of bulkier substrates. Competition assays using the fluorescent probe substrate 7-benzyloxymethoxy-4-trifluoromethylcoumarin (BOMFC) with these enzymes indicated inhibition of BOMFC metabolism by increasing concentrations of FPF. Furthermore, UPLC-MS/MS analysis revealed the capacity of all CYP6AQ1-like orthologs to metabolize FPF by hydroxylation in vitro at various levels, indicating a conserved FPF detoxification potential in different (neo)tropical bee species including Meliponini. This research, employing a toxicogenomics approach, provides important insights into the potential of stingless and other tropical bee species to detoxify FPF, and highlights the significance of investigating the detoxification mechanisms of insecticides in non-Apis bee species by molecular tools to inform risk assessment and conservation efforts.

Transcription factor CncC regulates the expression of antennal CYP6MU1 gene responsible for trans-2-hexen-1-al and nonanal recognition in Locusta migratoria

Cytochrome P450 monooxygenases (P450s) are a superfamily of multifunctional heme-containing proteins and could function as odorant-degrading enzymes (ODEs) in insect olfactory systems. In our previous study, we identified a P450 gene from the antennal transcriptome of Locusta migratoria, LmCYP6MU1, which could be induced by a variety of volatiles. However, the regulatory mechanisms of this gene in response to volatiles remain unknown. In current study, we investigated the tissues and development stages expression patterns of LmCYP6MU1 and determined its olfactory function in the recognition of the main host plant volatiles which induced LmCYP6MU1 expression. The results showed that LmCYP6MU1 was antenna-rich and highly expressed throughout the antennal developmental stages of locusts. LmCYP6MU1 played important roles in the recognition of trans-2-hexen-1-al and nonanal. Insect CncC regulates the expression of P450 genes. We tested whether LmCncC regulates LmCYP6MU1 expression. It was found that LmCncC knockdown in the antennae resulted in the downregulation of LmCYP6MU1 and repressed the volatiles-mediated induction of LmCYP6MU1. LmCncC knockdown reduced the electroantennogram (EAG) and behavioral responses of locusts to volatiles. These results suggested that LmCncC could regulate the basal and volatiles-mediated inducible expression of LmCYP6MU1 responsible for the recognition of trans-2-hexen-1-al and nonanal. These findings provide an original basis for understanding the regulation mechanisms of LmCncC on LmCYP6MU1 expression and help us better understand the LmCncC-mediated olfactory plasticity.

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 of five environmentally responsive gene families in a pine-feeding sawfly, Neodiprion lecontei (Hymenoptera: Diprionidae)

A central goal in evolutionary biology is to determine the predictability of adaptive genetic changes. Despite many documented cases of convergent evolution at individual loci, little is known about the repeatability of gene family expansions and contractions. To address this void, we examined gene family evolution in the redheaded pine sawfly Neodiprion lecontei, a noneusocial hymenopteran and exemplar of a pine-specialized lineage evolved from angiosperm-feeding ancestors. After assembling and annotating a draft genome, we manually annotated multiple gene families with chemosensory, detoxification, or immunity functions before characterizing their genomic distributions and molecular evolution. We find evidence of recent expansions of bitter gustatory receptor, clan 3 cytochrome P450, olfactory receptor, and antimicrobial peptide subfamilies, with strong evidence of positive selection among paralogs in a clade of gustatory receptors possibly involved in the detection of bitter compounds. In contrast, these gene families had little evidence of recent contraction via pseudogenization. Overall, our results are consistent with the hypothesis that in response to novel selection pressures, gene families that mediate ecological interactions may expand and contract predictably. Testing this hypothesis will require the comparative analysis of high-quality annotation data from phylogenetically and ecologically diverse insect species and functionally diverse gene families. To this end, increasing sampling in under-sampled hymenopteran lineages and environmentally responsive gene families and standardizing manual annotation methods should be prioritized.

Social isolation shortens lifespan through oxidative stress in ants

Social isolation negatively affects health, induces detrimental behaviors, and shortens lifespan in social species. Little is known about the mechanisms underpinning these effects because model species are typically short-lived and non-social. Using colonies of the carpenter ant Camponotus fellah, we show that social isolation induces hyperactivity, alters space-use, and reduces lifespan via changes in the expression of genes with key roles in oxidation-reduction and an associated accumulation of reactive oxygen species. These physiological effects are localized to the fat body and oenocytes, which perform liver-like functions in insects. We use pharmacological manipulations to demonstrate that the oxidation-reduction pathway causally underpins the detrimental effects of social isolation on behavior and lifespan. These findings have important implications for our understanding of how social isolation affects behavior and lifespan in general.

Effects of molting on the expression of ecdysteroid biosynthesis genes in the Y-organ of the blackback land crab, Gecarcinus lateralis

A pair of Y-organs (YOs) synthesize ecdysteroids that initiate and coordinate molting processes in decapod crustaceans. The YO converts cholesterol to secreted products through a biosynthetic pathway involving a Rieske oxygenase encoded by Neverland (Nvd) and cytochrome P450 monooxygenases encoded by Halloween genes Spook (Spo; Cyp307a1), Phantom (Phm; Cyp306a1), Disembodied (Dib; Cyp302a1), and Shadow (Sad; Cyp315a1). NAD kinase (NADK) and 5-aminolevulinic acid synthase (ALAS) support ecdysteroid synthesis in insects. A 20-hydroxylase, encoded by Shed in decapods and Shade in insects, converts ecdysone to the active hormone 20-hydroxyecdysone (20E). 20E is inactivated by cytochrome P450 26-hydroxylase (Cyp18a1). Contigs encoding these eight proteins were extracted from a Gecarcinus lateralis YO transcriptome and their expression was quantified by quantitative polymerase chain reaction. mRNA levels of Gl-Spo and Gl-Phm were four orders of magnitude higher in YO than those in nine other tissues, while mRNA levels of Gl-NADK and Gl-ALAS were similar in all ten tissues. In G. lateralis induced to molt by multiple leg autotomy, YO mRNA levels of Gl-Nvd, Gl-Spo, Gl-Phm, Gl-NADK, and Gl-ALAS were highest in intermolt and premolt stages and lower in postmolt. Gl-Dib mRNA level was not affected by molt stage. mRNA level of Gl-Sad, which converts 2-deoxyecdysone to ecdysone, was higher in mid- and late premolt stages, when YO ecdysteroidogenic capacity is greatest. Gl-Cyp18a1 mRNA level was highest in intermolt, decreased in premolt stages, and was lowest in postmolt. In animals induced to molt by eyestalk ablation, YO mRNA levels of all eight genes were not correlated with increased hemolymph 20E titers. These results suggest that YO ecdysteroidogenic genes are differentially regulated at transcriptional and translational levels.

A detoxification pathway initiated by a nuclear receptor TcHR96h in Tetranychus cinnabarinus (Boisduval)

Understanding the mechanism of detoxification initiation in arthropods after pesticide exposure is crucial. Although the identity of transcription factors that induce and regulate the expression of detoxification genes in response to pesticides is beginning to emerge, whether transcription factors directly interact with xenobiotics is unclear. The findings of this study revealed that a nuclear hormone receptor, Tetranychus cinnabarinus hormone receptor (HR) TcHR96h, regulates the overexpression of the detoxification gene TcGSTm02, which is involved in cyflumetofen resistance. The nuclear translocation of TcHR96h increased after cyflumetofen exposure, suggesting direct binding with cyflumetofen. The direct binding of TcHR96h and cyflumetofen was supported by several independent proteomic assays that quantify interactions with small molecules. Together, this study proposes a model for the initiation of xenobiotic detoxification in a polyphagous agricultural pest. These insights not only provide a better understanding of the mechanisms of xenobiotic detoxification and metabolism in arthropods, but also are crucial in understanding adaptation in polyphagous herbivores.

Identification of metabolizing enzyme genes associated with xenobiotics and odorants in the predatory stink bug Arma custos based on transcriptome analysis

The predatory stink bug, Arma custos, is a highly effective beneficial predator of crop pests. The lack of gene information related to xenobiotic detoxification and odorant degrading enzymes in the predator stink bugs to date has limited our ability for more in-depth studies of biological control. Hence, we conducted de novo assembly of the A. custos transcriptome from guts, antennae, and other tiussue samples of 5th instar larvae using Illumina sequencing technology. A total of 91, 50 and 23 genes of cytochrome P450 monooxygenases (CYPs), carboxyl/choline esterases (CCEs) and glutathione S-transferases (GSTs) genes were identified, respectively. Gene expansions of CYP3 and CYP4 clans and the hormone and pheromone processing CCE class were found in A. custos. Analysis of tissue-specific expression patterns showed that 37 CYPs, 14 CCEs and 8 GSTs were enriched in guts, and 6 CYPs, 5 CCEs and 2 GSTs were up-regulated in antennae, suggesting their potential roles on xenobiotics detoxification and ordorant degradation. Gene information data presented here could be useful for a deeper understanding of the ecology, physiology and behavior of this beneficial species and could be helpful to improve their bio-control efficiency.

A review of the molecular mechanisms of acaricide resistance in mites and ticks

The Arachnida subclass of Acari comprises many harmful pests that threaten agriculture as well as animal health, including herbivorous spider mites, the bee parasite Varroa, the poultry mite Dermanyssus and several species of ticks. Especially in agriculture, acaricides are often used intensively to minimize the damage they inflict, promoting the development of resistance. Beneficial predatory mites used in biological control are also subjected to acaricide selection in the field. The development and use of new genetic and genomic tools such as genome and transcriptome sequencing, bulked segregant analysis (QTL mapping), and reverse genetics via RNAi or CRISPR/Cas9, have greatly increased our understanding of the molecular genetic mechanisms of resistance in Acari, especially in the spider mite Tetranychus urticae which emerged as a model species. These new techniques allowed to uncover and validate new resistance mutations in a larger range of species. In addition, they provided an impetus to start elucidating more challenging questions on mechanisms of gene regulation of detoxification associated with resistance.

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

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

Identification and functional study of detoxification-related genes in response to tolfenpyrad stress in Glyphodes pyloalis Walker (Lepidoptera: Pyralidae)

Glyphodes pyloalis Walker (G. pyloalis) is a common destructive mulberry pest. Due to the long-term and frequent use of insecticides, it has developed tolerance to commonly used insecticides. Tolfenpyrad (TFP) is a novel pyrazole heterocyclic insecticide. In order to understand the TFP detoxification mechanism of G. pyloalis larvae, we first estimated the LC30 dose of TFP for 3rd instar G. pyloalis larvae. Next, we identified genes that were differentially expressed in 3rd instar G. pyloalis larvae treated with TFP compared to the control group by transcriptome sequencing. In total, 86,949,569 and 67,442,028 clean reads were obtained from TFP-treated and control G. pyloalis larvae, respectively. A total of 5588 differentially expressed genes (DEGs) were identified in TFP-treated and control G. pyloalis larvae, of which 3084 genes were upregulated and 2504 genes were downregulated. We analyzed the expression of 43 candidate detoxification enzyme genes associated with insecticide tolerance using qPCR. According to the spatiotemporal expression pattern of DEGs, we found that CYP6ABE1, CYP333A36 and GST-epsilon8 were highly expressed in the midgut, while CarEs14 was strongly expressed in haemolymph. Furthermore, we successfully knocked down these genes by RNA interference. After silencing CYP6ABE1 and CYP333A36, bioassay showed that the mortality rate of TFP-treated G. pyloalis larvae was significantly higher compared to the control group. This study provides a theoretical foundation for understanding the sensitivity of G. pyloalis to TFP and establish the basis for the effective and green management of this pest.

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

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

Comparative evolutionary analyses of eight whitefly Bemisia tabaci sensu lato genomes: cryptic species, agricultural pests and plant-virus vectors

BACKGROUND

The group of > 40 cryptic whitefly species called Bemisia tabaci sensu lato are amongst the world's worst agricultural pests and plant-virus vectors. Outbreaks of B. tabaci s.l. and the associated plant-virus diseases continue to contribute to global food insecurity and social instability, particularly in sub-Saharan Africa and Asia. Published B. tabaci s.l. genomes have limited use for studying African cassava B. tabaci SSA1 species, due to the high genetic divergences between them. Genomic annotations presented here were performed using the 'Ensembl gene annotation system', to ensure that comparative analyses and conclusions reflect biological differences, as opposed to arising from different methodologies underpinning transcript model identification.

RESULTS

We present here six new B. tabaci s.l. genomes from Africa and Asia, and two re-annotated previously published genomes, to provide evolutionary insights into these globally distributed pests. Genome sizes ranged between 616-658 Mb and exhibited some of the highest coverage of transposable elements reported within Arthropoda. Many fewer total protein coding genes (PCG) were recovered compared to the previously published B. tabaci s.l. genomes and structural annotations generated via the uniform methodology strongly supported a repertoire of between 12.8-13.2 × 10³ PCG. An integrative systematics approach incorporating phylogenomic analysis of nuclear and mitochondrial markers supported a monophyletic Aleyrodidae and the basal positioning of B. tabaci Uganda-1 to the sub-Saharan group of species. Reciprocal cross-mating data and the co-cladogenesis pattern of the primary obligate endosymbiont 'Candidatus Portiera aleyrodidarum' from 11 Bemisia genomes further supported the phylogenetic reconstruction to show that African cassava B. tabaci populations consist of just three biological species. We include comparative analyses of gene families related to detoxification, sugar metabolism, vector competency and evaluate the presence and function of horizontally transferred genes, essential for understanding the evolution and unique biology of constituent B. tabaci. s.l species.

CONCLUSIONS

These genomic resources have provided new and critical insights into the genetics underlying B. tabaci s.l. biology. They also provide a rich foundation for post-genomic research, including the selection of candidate gene-targets for innovative whitefly and virus-control strategies.

High-quality chromosome-level scaffolds of the plant bug Pachypeltis micranthus provide insights into the availability of Mikania micrantha control

BACKGROUND

The plant bug, Pachypeltis micranthus Mu et Liu (Hemiptera: Miridae), is an effective potential biological control agent for Mikania micrantha H.B.K. (Asteraceae; one of the most notorious invasive weeds worldwide). However, limited knowledge about this species hindered its practical application and research. Accordingly, sequencing the genome of this mirid bug holds great significance in controlling M. micrantha.

RESULTS

Here, 712.72 Mb high-quality chromosome-level scaffolds of P. micranthus were generated, of which 707.51 Mb (99.27%) of assembled sequences were anchored onto 15 chromosome-level scaffolds with contig N50 of 16.84 Mb. The P. micranthus genome had the highest GC content (42.43%) and the second highest proportion of repetitive sequences (375.82 Mb, 52.73%) than the three other mirid bugs (i.e., Apolygus lucorum, Cyrtorhinus lividipennis, and Nesidiocoris tenuis). Phylogenetic analysis showed that P. micranthus clustered with other mirid bugs and diverged from the common ancestor approximately 200 million years ago. Gene family expansion and/or contraction were analyzed, and significantly expanded gene families associated with P. micranthus feeding and adaptation to M. micrantha were manually identified. Compared with the whole body, transcriptome analysis of the salivary gland revealed that most of the upregulated genes were significantly associated with metabolism pathways and peptidase activity, particularly among cysteine peptidase, serine peptidase, and polygalacturonase; this could be one of the reasons for precisely and highly efficient feeding by the oligophagous bug P. micranthus on M. micrantha.

CONCLUSION

Collectively, this work provides a crucial chromosome-level scaffolds resource to study the evolutionary adaptation between mirid bug and their host. It is also helpful in searching for novel environment-friendly biological strategies to control M. micrantha.

Genome-Wide Identification of Detoxification Genes in Wild Silkworm Antheraea pernyi and Transcriptional Response to Coumaphos

For a half-century, the commercial wild silkworm, Antheraea pernyi, has been protected by coumaphos, which is an internal organophosphorus insecticide used to kill the potential parasitic fly larvae inside. Knowledge about the detoxification genes of A. pernyi as well as the detoxification mechanism for this species remains severely limited. In this study, we identified 281 detoxification genes (32 GSTs, 48 ABCs, 104 CYPs, and 97 COEs) in the genome of this insect, which are unevenly distributed over 46 chromosomes. When compared to the domesticated silkworm, Bombyx mori, a lepidopteran model species, A. pernyi has a similar number of ABCs, but a greater number of GSTs, CYPs, and COEs. By transcriptome-based expression analysis, we found that coumaphos at a safe concentration level significantly changed the pathways related to ATPase complex function and the transporter complex in A. pernyi. KEGG functional enrichment analysis indicated that protein processing in the endoplasmic reticulum was the most affected pathway after coumaphos treatment. Finally, we identified four significantly up-regulated detoxification genes (ABCB1, ABCB3, ABCG11, and ae43) and one significantly down-regulated detoxification gene (CYP6AE9) in response to coumaphos treatment, suggesting that these five genes may contribute to detoxification of coumaphos in A. pernyi. Our study provides the first set of detoxification genes for wild silkworms from Saturniidae and highlights the importance of detoxification gene repertoire in insect pesticide tolerance.

Knockdown of CYP9A9 increases the susceptibility to lufenuron, methoxyfenozide and a mixture of both in Spodoptera exigua

Lufenuron (LUF) and Methoxyfenozide (MET) as Insect Growth Regulators (IGRs) contribute to the current control of the catastrophic crop pest, Spodoptera exigua (Lepidoptera, Noctuidae). Yet S. exigua has evolved resistance to LUF and MET, which is possibly mediated by cytochrome P450 monooxygenases (P450s), particularly from the CYP3 clade family, as it plays a key role in the detoxification of insecticides. However, a mixture of LUF and MET (MML) (optimal ratio: 6:4) remains highly insecticidal. Here, we analysed the response of S. exigua to sublethal concentrations of LUF, MET, and MML via transcriptomics. Twelve differentially expressed genes (DEGs) encoding CYP3 clade members were observed in transcriptomes and CYP9A9 was significantly upregulated after treatment with LUF, MET, and MML. Further, CYP9A9 was most highly expressed in the midgut of L4 S. exigua larvae. RNAi-mediated knockdown of CYP9A9 reduced the activity of CYP450 and increased the susceptibility of S. exigua larvae to LUF, MET, and MML. Thus, CYP9A9 plays a key role in the detoxification of LUF, MET, and MML in S. exigua. These findings provide new insights into insecticidal actions of IGRs, which can be applied to the establishment of novel pest management strategies.

Increased metabolism in combination with the novel cytochrome b target-site mutation L258F confers cross-resistance between the Qo inhibitors acequinocyl and bifenazate in Tetranychus urticae

Acequinocyl and bifenazate are potent acaricides acting at the Q_o site of complex III of the electron transport chain, but frequent applications of these acaricides have led to the development of resistance in spider mites. Target-site resistance caused by mutations in the conserved cd1- and ef-helices of the Q_o pocket of cytochrome b has been elucidated as the main resistance mechanism. We therefore monitored Q_o pocket mutations in European field populations of Tetranychus urticae and uncovered a new mutation, L258F. The role of this mutation was validated by revealing patterns of maternal inheritance and by the independently replicated introgression in an unrelated susceptible genetic background. However, the parental strain exhibited higher resistance levels than conferred by the mutation alone in isogenic lines, especially for acequinocyl, implying the involvement of strong additional resistance mechanisms. This was confirmed by revealing a polygenic inheritance pattern with classical genetic crosses and via synergism experiments. Therefore, a genome-wide expression analysis was conducted that identified a number of highly overexpressed detoxification genes, including many P450s. Functional expression revealed that the P450 CYP392A11 can metabolize bifenazate by hydroxylation of the ring structure. In conclusion, the novel cytochrome b target-site mutation L258F was uncovered in a recently collected field strain and its role in acequinocyl and bifenazate resistance was validated. However, the high level of resistance in this strain is most likely caused by a combination of target-site resistance and P450-based increased detoxification, potentially acting in synergism.

Biochemical characterisation of Cytochrome P450 oxidoreductase from the cattle tick, Rhipicephalus microplus, highlights potential new acaricide target

Management of the cattle tick, Rhipicephalus microplus, presents a challenge because some populations of this cosmopolitan and economically important ectoparasite are resistant to multiple classes of acaricides. Cytochrome P450 oxidoreductase (CPR) is part of the cytochrome P450 (CYP450) monooxygenases that are involved in metabolic resistance by their ability to detoxify acaricides. Inhibiting CPR, the sole redox partner that transfers electrons to CYP450s, could overcome this type of metabolic resistance. This report represents the biochemical characterisation of a CPR from ticks. Recombinant CPR of R. microplus (RmCPR), minus its N-terminal transmembrane domain, was produced in a bacterial expression system and subjected to biochemical analyses. RmCPR displayed a characteristic dual flavin oxidoreductase spectrum. Incubation with nicotinamide adenine dinucleotide phosphate (NADPH) lead to an increase in absorbance between 500 and 600 nm with a corresponding appearance of a peak absorbance at 340-350 nm indicating functional transfer of electrons between NADPH and the bound flavin cofactors. Using the pseudoredox partner, kinetic parameters for both cytochrome c and NADPH binding were calculated as 26.6 ± 11.4 µM and 7.03 ± 1.8 µM, respectively. The turnover, K_cat, for RmCPR for cytochrome c was calculated as 0.08 s^-1 which is significantly lower than the CPR homologues of other species. IC₅₀ (Half maximal Inhibitory Concentration) values obtained for the adenosine analogues 2', 5' ADP, 2'- AMP, NADP⁺and the reductase inhibitor diphenyliodonium were: 140, 82.2, 24.5, and 75.3 µM, respectively. Biochemically, RmCPR resembles CPRs of hematophagous arthropods more so than mammalian CPRs. These findings highlight the potential of RmCPR as a target for the rational design of safer and potent acaricides against R. microplus.

A conserved hymenopteran-specific family of cytochrome P450s protects bee pollinators from toxic nectar alkaloids

Many plants produce chemical defense compounds as protection against antagonistic herbivores. However, how beneficial insects such as pollinators deal with the presence of these potentially toxic chemicals in nectar and pollen is poorly understood. Here, we characterize a conserved mechanism of plant secondary metabolite detoxification in the Hymenoptera, an order that contains numerous highly beneficial insects. Using phylogenetic and functional approaches, we show that the CYP336 family of cytochrome P450 enzymes detoxifies alkaloids, a group of potent natural insecticides, in honeybees and other hymenopteran species that diverged over 281 million years. We linked this function to an aspartic acid residue within the main access channel of CYP336 enzymes that is highly conserved within this P450 family. Together, these results provide detailed insights into the evolution of P450s as a key component of detoxification systems in hymenopteran species and reveal the molecular basis of adaptations arising from interactions between plants and beneficial insects.