Integrated analysis of cytochrome P450 gene superfamily in the red flour beetle, Tribolium castaneum

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.

The Role of TcCYP6K1 and TcCYP9F2 Influences Trehalose Metabolism under High-CO2 Stress in Tribolium castaneum (Coleoptera)

Cytochrome P450 monooxygenases (CYP), crucial detoxification enzymes in insects, are involved in the metabolism of endogenous substances as well as the activation and degradation of exogenous compounds. In this study, T. castaneum was utilized to investigate the roles of TcCYP6K1 and TcCYP9F2 genes influencing in the trehalose metabolism pathway under high-CO2 stress. By predicting the functional sequences of TcCYP6K1 and TcCYP9F2 genes and analyzing their spatiotemporal expression patterns, it was discovered that both genes belong to the CYP3 group and exhibit high expression levels during the larval stage, decreasing during the pupal stage, while showing high expression in the fatty body, intestine, and malpighian tubules. Furthermore, following the knockdown of TcCYP6K1 and TcCYP9F2 genes in combination with treating larvae with 75% CO2, it was observed that larval mortality increased, and glycogen content significantly decreased, while trehalose content increased significantly. Additionally, membrane-bound trehalase enzyme activity declined, TPS gene expression was significantly upregulated, GS gene expression was significantly downregulated, and ATP content showed a marked decrease. In conclusion, CYP genes are critical responsive genes of T. castaneum to high CO2 levels, potentially impacting the insect's resistance to carbon dioxide through their involvement in the synthesis or breakdown of the carbohydrate metabolism pathway. These findings could serve as a theoretical basis for the utilization of novel pesticides in low-oxygen grain storage techniques and offer new insights for environmentally friendly pest control strategies in grain storage.

Population-specific thermal responses contribute to regional variability in arbovirus transmission with changing climates

Temperature is increasing globally, and vector-borne diseases are particularly responsive to such increases. While it is known that temperature influences mosquito life history traits, transmission models have not historically considered population-specific effects of temperature. We assessed the interaction between Culex pipiens population and temperature in New York State (NYS) and utilized novel empirical data to inform predictive models of West Nile virus (WNV) transmission. Genetically and regionally distinct populations from NYS were reared at various temperatures, and life history traits were monitored and used to inform trait-based models. Variation in Cx. pipiens life history traits and population-dependent thermal responses account for a predicted 2.9°C difference in peak transmission that is reflected in regional differences in WNV prevalence. We additionally identified genetic signatures that may contribute to distinct thermal responses. Together, these data demonstrate how population variation contributes to significant geographic variability in arbovirus transmission with changing climates.

Strong resistance to β-cyfluthrin in a strain of the beetle Alphitobius diaperinus: a de novo transcriptome analysis

The lesser mealworm, Alphitobius diaperinus, is an invasive tenebrionid beetle and a vector of pathogens. Due to the emergence of insecticide resistance and consequent outbreaks that generate significant phytosanitary and energy costs for poultry farmers, it has become a major insect pest worldwide. To better understand the molecular mechanisms behind this resistance, we studied a strain of A. diaperinus from a poultry house in Brittany that was found to be highly resistant to the β-cyfluthrin. The strain survived β-cyfluthrin exposures corresponding to more than 100 times the recommended dose. We used a comparative de novo RNA-Seq approach to explore genes expression in resistant versus sensitive strains. Our de novo transcriptomic analyses showed that responses to β-cyfluthrin likely involved a whole set of resistance mechanisms. Genes related to detoxification, metabolic resistance, cuticular hydrocarbon biosynthesis and proteolysis were found to be constitutively overexpressed in the resistant compared to the sensitive strain. Follow-up enzymatic assays confirmed that the resistant strain exhibited high basal activities for detoxification enzymes such as cytochrome P450 monooxygenase and glutathione-S-transferase. The in-depth analysis of differentially expressed genes suggests the involvement of complex regulation of signaling pathways. Detailed knowledge of these resistance mechanisms is essential for the establishment of effective pest control.

Comparative transcriptomics reveals the conservation and divergence of reproductive genes across three sympatric Tomicus bark beetles

Three tree-killing bark beetles belonging to the genus Tomicus, Tomicus yunnanensis, Tomicus brevipilosus and Tomicus minor (Coleoptera; Curculionidae, Scolytinae), are serious wood-borers with larvae feeding on the phloem tissues of Pinus yunnanensis. The three Tomicus beetles, in some cases, coexist in a same habitat, providing a best system for exploring the conservation and divergence of reproductive genes. Here, we applied comparative transcriptomics and molecular biology approaches to characterize reproductive-related genes in three sympatric Tomicus species. Illumina sequencing of female and male reproductive systems and residual bodies generated a large number of clean reads, representing 185,920,232 sequences in T. yunnanensis, 169,153,404 in T. brevipilosus and 178,493,176 in T. minor that were assembled into 32,802, 56,912 and 33,670 unigenes, respectively. The majority of the genes had detectable expression in reproductive tissues (FPKM >1), particularly those genes in T. brevipilosus accounting for 76.61 % of the total genes. From the transcriptomes, totally 838 genes encoding 463 detoxification enzymes, 339 chemosensory membrane proteins and 36 ionotropic glutamate receptors (iGluRs) were identified, including 622 reproductive tissue-expressed genes. Of these, members of carboxylesterases (COEs), ionotropic receptors (IRs), sensory neuron membrane proteins (SNMPs) and iGluRs were highly conserved in gene numbers and sequence identities across three Tomicus species. Further, expression profiling analyses revealed a number of genes expressed in reproductive tissues and the diverse expression characteristics in these beetles. The results provide evidence for the conservation and differences of reproductive genes among three sympatric closely related beetles, helping understand their different reproductive strategies and the maximization of the reproductive success.

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.

Overexpression of cytochrome P450 and esterase genes involved in permethrin resistance in larvae and adults of Culex quinquefasciatus

Mosquitoes are important vectors of several arthropod-borne diseases, which remain a priority for epidemiological research. Mosquito vector control strategies have traditionally relied on chemical insecticides such as synthetic pyrethroids. However, the indiscriminate use of pesticides has resulted in the development of resistance in many mosquito species. In insects, resistance evolves primarily through the overexpression of one or more gene products from the cytochrome P450, carboxylesterase, and glutathione superfamilies. The current study examined the expression of cytochrome P450 CYP6M2, CYP6AA7, CYP6Z2, CYP9J34, α-Esterase, Esterase B1, and neuroactin genes in larvae and adults of a permethrin-resistant (PerRes) and susceptible (Sus) Culex quinquefasciatus strains. The results showed that the CYP6AA7 gene was overexpressed (10-fold) in larvae and adults with PerRes (p < 0.01) followed by CYPJ34 (9.0-fold) and CYP6Z2 (5.0-fold) compared to the Sus, whereas fewer changes in CYP6M gene expression were observed in PerRes adults (p < 0.05), and no expression was found in larvae. The esterase gene was overexpressed in PerRes larvae (9.0-fold) followed by adults (2.5-fold) compared to the susceptible strain. Based on data, the present study suggests that cytochrome P450, CYP6AA7, CYP6Z2, CYP9J34, α-Esterase, Esterase B1, and neuroactin genes were involved in permethrin resistance in larval and adult Cx. quinquefasciatus.

Function of Cytochrome P450s and Gut Microbiome in Biopesticide Adaptation of Grapholita molesta on Different Host Diets

Insects that feed on various host plants possess diverse xenobiotic adaptations; however, the underlying mechanisms are poorly understood. In the present study, we used Grapholita molesta, which shifts feeding sites from peach shoots to apple fruits, as a model to explore the effects of shifts in host plant diet on the profiles of cytochrome P450s and the gut bacteria microbiome, as well as their effects on biopesticide adaptation. We found that the sensitivity of the fruit-feeding G. molesta to emamectin benzoate biopesticide was significantly lower than that of the shoot-feeding larvae. We also found that the P450 enzyme activity and the expression of nine cytochrome P450s were enhanced in G. molesta fed on Fuji apples compared to those fed on peach shoots. The survival rates of G. molesta exposed to emamectin benzoate significantly decreased as each of three of four emamectin benzoate-inducted cytochrome P450 genes were silenced. Furthermore, we discovered the gut bacteria dynamics of G. molesta changed with the host shift and the structure of the gut bacteria microbiome was determined by the final diet ingested; additionally, the dysbiosis of the gut microbiota induced by antibiotics could significantly increase the sensitivity to emamectin benzoate. Taken together, our results suggest that the expression of P450s and the composition of the gut bacteria microbiome promote adaptation to emamectin benzoate in G. molesta, providing new insights into the molecular mechanisms underlying xenobiotic adaptation in this notorious pest.

Roles of insect odorant binding proteins in communication and xenobiotic adaptation

Odorant binding proteins (OBPs) are small water-soluble proteins mainly associated with olfaction, facilitating the transport of odorant molecules to their relevant receptors in the sensillum lymph. While traditionally considered essential for olfaction, recent research has revealed that OBPs are engaged in a diverse range of physiological functions in modulating chemical communication and defense. Over the past 10 years, emerging evidence suggests that OBPs play vital roles in purifying the perireceptor space from unwanted xenobiotics including plant volatiles and pesticides, potentially facilitating xenobiotic adaptation, such as host location, adaptation, and pesticide resistance. This multifunctionality can be attributed, in part, to their structural variability and effectiveness in transporting, sequestering, and concealing numerous hydrophobic molecules. Here, we firstly overviewed the classification and structural properties of OBPs in diverse insect orders. Subsequently, we discussed the myriad of functional roles of insect OBPs in communication and their adaptation to xenobiotics. By synthesizing the current knowledge in this field, our review paper contributes to a comprehensive understanding of the significance of insect OBPs in chemical ecology, xenobiotic adaptation, paving the way for future research in this fascinating area of study.

Exploiting fungi in bioremediation for cleaning-up emerging pollutants in aquatic ecosystems

Aquatic pollution negatively affects water bodies, marine ecosystems, public health, and economy. Restoration of contaminated habitats has attracted global interest since protecting the health of marine ecosystems is crucial. Bioremediation is a cost-effective and eco-friendly way of transforming hazardous, resistant contaminants into environmentally benign products using diverse biological treatments. Because of their robust morphology and broad metabolic capabilities, fungi play an important role in bioremediation. This review summarizes the features employed by aquatic fungi for detoxification and subsequent bioremediation of different toxic and recalcitrant compounds in aquatic ecosystems. It also details how mycoremediation may convert chemically-suspended matters, microbial, nutritional, and oxygen-depleting aquatic contaminants into ecologically less hazardous products using multiple modes of action. Mycoremediation can also be considered in future research studies on aquatic, including marine, ecosystems as a possible tool for sustainable management, providing a foundation for selecting and utilizing fungi either independently or in microbial consortia.

Sublethal effects of halofenozide on larval development and detoxification in Phaedon brassicae (Coleoptera: Chrysomelidae)

The brassica leaf beetle, Phaedon brassicae, is a serious defoliator of cruciferous crops. Halofenozide (Hal), an ecdysone agonist, is a new class of insect growth-regulating insecticide. Our preliminary experiment revealed the outstanding larval toxicity of Hal against P. brassicae. However, the metabolic degradation of this compound in insects remains unclear. In this study, oral administration of Hal at LC10 and LC25 caused severe separation of the cuticle and epidermis, leading to larval molting failure. Sublethal dose exposure also significantly reduced the larval respiration rate as well as their pupation rates and pupal weights. Conversely, the activities of the multifunctional oxidase, carboxylesterase (CarE), and glutathione S-transferase (GST) were significantly enhanced in Hal-treated larvae. Further analysis using RNA sequencing identified 64 differentially expressed detoxifying enzyme genes, including 31 P450s, 13 GSTs, and 20 CarEs. Among the 25 upregulated P450s, 22 genes were clustered into the CYP3 clan, and the other 3 genes belonged to the CYP4 clan. Meanwhile, 3 sigma class GSTs and 7 epsilon class GSTs were dramatically increased, accounting for the majority of the upregulated GSTs. Moreover, 16 of the 18 overexpressed CarEs were clustered into the coleopteran xenobiotic-metabolizing group. These results showed the augmented expression of detoxification genes in P. brassicae after exposed to sublethal dose of Hal, and helped to better understand the potential metabolic pathways that could contribute to the reduced sensitivity to Hal in this pest. Overall, a deep insight into the detoxification mechanisms would provide practical guidance for the field management of P. brassicae.

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 Analysis of the Cytochrome P450 Gene Family in Bemisia tabaci MED and Their Roles in the Insecticide Resistance

The whitefly, Bemisia tabaci MED (Hemiptera: Aleyrodidae), is an omnivorous agricultural pest, which causes huge economic losses to agriculture and is highly resistant to many pesticides. The overexpression of cytochrome P450 may play an important role in host adaptation and insecticide resistance in B. tabaci MED. Therefore, the present study systematically analyzed the cytochrome P450 gene family at the genome-wide level to understand its function in B. tabaci MED. Our analysis identified 58 cytochrome P450 genes in B. tabaci MED, among which 24 were novel. Phylogenetic analysis revealed broad functional and species-specific diversification in B. tabaci MED P450, suggesting the role of multiple P450 genes in detoxifying. Reverse transcription-real time quantitative PCR (RT-qPCR) showed that CYP4CS2, CYP4CS5, CYP4CS6, CYP4CS8, CYP6DW4, CYP6DW5, CYP6DW6, CYP6DZ8, and CYP6EN1 genes increased significantly after two days of exposure to imidacloprid. Interestingly, all nine genes belonged to the CYP4 and CYP6 families. A decrease in the expression of five genes (CYP6DW4, CYP6DW5, CYP6DW6, CYP6DZ8, and CYP4CS6) via RNA interference (RNAi) resulted in a significant increase in the mortalities of whiteflies when exposed to imidacloprid. These results indicate that the overexpression of the P450 genes may play an essential role in imidacloprid tolerance of B. tabaci MED. Thus, the present study provides basic information on P450 genes in B. tabaci MED, which will further help elucidate the insecticide resistance mechanism in the agricultural pest whitefly.

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.

Insight into weevil biology from a reference quality genome of the boll weevil, Anthonomus grandis grandis Boheman (Coleoptera: Curculionidae)

The boll weevil, Anthonomus grandis grandis Boheman, is one of the most historically impactful insects due to its near destruction of the US cotton industry in the early 20th century. Contemporary efforts to manage this insect primarily use pheromone baited traps for detection and organophosphate insecticides for control, but this strategy is not sustainable due to financial and environmental costs. We present a high-quality boll weevil genome assembly, consisting of 306 scaffolds with approximately 24,000 annotated genes, as a first step in the identification of gene targets for novel pest control. Gene content and transposable element distribution are similar to those found in other Curculionidae genomes; however, this is the most contiguous and only assembly reported to date for a member in the species-rich genus Anthonomus. Transcriptome profiles across larval, pupal, and adult life stages led to identification of several genes and gene families that could present targets for novel control strategies.

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

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

Identification and characterization of CYP307A1 as a molecular target for controlling the small hive beetle, Aethina tumida

BACKGROUND

The molting hormone 20-hydroxyecdysone (20E) plays a key role in insect development, metamorphosis, and reproduction. Previous studies have shown that ecdysteroid metabolism is regulated by a series of CYP genes in most of the insect species. However, the roles of these CYP genes in a Coleopteran beetle, Aethina tumida (small hive beetle, SHB) have not yet been explored.

RESULTS

In the current study, we identified seven CYP genes (six Halloween genes and one AtCYP18A1 gene) related to 20E metabolism. Reverse transcription quantitative polymerase chain reaction (RT-qPCR) showed that AtCYP307A1 and AtCYP307B1 were primarily expressed in the embryonic stage and in the cephalothorax of larvae. RNA interference (RNAi) screening revealed that suppression of AtCYP307A1 expression caused a lethal phenotype during the larval-pupal metamorphosis. Furthermore, Hematoxylin and Eosin staining of the integument showed that the RNAi of AtCYP307A1 inhibited the apolysis and degradation of the old cuticle. In addition, silencing of AtCYP307A1 resulted in significant down-regulation of 20E titers and the expression levels of 20E signaling pathway genes. Finally, the AtCYP307A1 RNAi phenotype was rescued by topical application of 20E.

CONCLUSION

Our studies suggest that AtCYP307A1 involved in 20E synthesis is indispensable during the larval-pupal metamorphosis of beetles, which could serve as a putative insecticide target for pest control. © 2022 Society of Chemical Industry.

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.

Comparative Transcriptome Analysis to Reveal Differentially Expressed Cytochrome P450 in Response to Imidacloprid in the Aphid Lion, Chrysoperla zastrowi sillemi (Esben-Petersen)

The aphid lion, Chrysoperla zastrowi sillemi (Neuroptera: Chrysopidae) is a highly effective beneficial predator of many agricultural pests and has developed resistance to several insecticides. Understanding the molecular mechanism of insecticide resistance in the predators is crucial for its effective application in IPM programs. Therefore, transcriptomes of imidacloprid-resistant and susceptible strains have been assessed using RNA-seq. Cytochrome P450 is one of the important gene families involved in xenobiotic metabolism. Hence, our study focused on the CYP gene family where mining, nomenclature, and phylogenetic analysis revealed a total of 95 unique CYP genes with considerable expansion in CYP3 and CYP4 clans. Further, differential gene expression (DGE) analysis revealed ten CYP genes from CYP3 and CYP4 clans to be differentially expressed, out of which nine genes (CYP4419A1, CYP4XK1, CYP4416A10, CYP4416A-fragment8, CYP6YL1, CYP6YH6, CYP9GK-fragment16, CYP9GN2, CYP9GK6) were downregulated and one (CYP9GK3) was upregulated in the resistant strain as compared to the susceptible strain. Expression validation by quantitative real-time PCR (qRT-PCR) is consistent with the DGE results. The expansion and differential expression of CYP genes may be an indicator of the capacity of the predator to detoxify a particular group of insecticides.

Identification and characterization of candidate detoxification genes in Pharsalia antennata Gahan (Coleoptera: Cerambycidae)

The wood-boring beetles, including the majority of Cerambycidae, have developed the ability to metabolize a variety of toxic compounds derived from host plants and the surrounding environment. However, detoxification mechanisms underlying the evolutionary adaptation of a cerambycid beetle Pharsalia antennata to hosts and habitats are largely unexplored. Here, we characterized three key gene families in relation to detoxification (cytochrome P450 monooxygenases: P450s, carboxylesterases: COEs and glutathione-S-transferases: GSTs), by combinations of transcriptomics, gene identification, phylogenetics and expression profiles. Illumina sequencing generated 668,701,566 filtered reads in 12 tissues of P. antennata, summing to 100.28 gigabases data. From the transcriptome, 215 genes encoding 106 P450s, 77 COEs and 32 GSTs were identified, of which 107 relatives were differentially expressed genes. Of the identified 215 genes, a number of relatives showed the orthology to those in Anoplophora glabripennis, revealing 1:1 relationships in 94 phylogenetic clades. In the trees, P. antennata detoxification genes mainly clustered into one or two subfamilies, including 64 P450s in the CYP3 clan, 33 COEs in clade A, and 20 GSTs in Delta and Epsilon subclasses. Combining transcriptomic data and PCR approaches, the numbers of detoxification genes expressed in abdomens, antennae and legs were 188, 148 and 141, respectively. Notably, some genes exhibited significantly sex-biased levels in antennae or legs of both sexes. The findings provide valuable reference resources for further exploring xenobiotics metabolism and odorant detection in P. antennata.

Chromosome-level genome assembly of Monochamus saltuarius reveals its adaptation and interaction mechanism with pine wood nematode

Monochamus saltuarius (Coleoptera: Cerambycidae) was reported as the vector beetle of the pine wood nematode (PWN, Bursaphelenchus xylophilus) in Japan and Europe. It was first reported to transmitted the PWN to native Pinus species in 2018 in Liaoning Province, China. However, the lack of genomic resources has limited the in-depth understanding of its interspecific relationship with PWN. Here, we obtained a chromosome-level reference genome of M. saltuarius combining Illumina, Nanopore and Hi-C sequencing technologies. We assembled the scaffolds into ten chromosomes (including an X chromosome) and obtained a 682.23 Mb chromosome-level genome with a N50 of 73.69 Mb. In total, 427.67 Mb (62.69 %) repeat sequences were identified and 14, 492 protein-coding genes were predicted, of which 93.06 % were annotated. We described the mth/mthl, P450, OBP and OR gene families associated with the vector beetle's development and resistance, as well as the host selection and adaptation, which serve as a valuable resource for understanding the host adaptation in insects during evolution. This high quality reference genome of M. saltuarius also provide new avenues for researching the mechanism of this synergistic damage between vector beetles and PWN.

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.

Three cytochrome P450 CYP4 family genes regulated by the CncC signaling pathway mediate phytochemical susceptibility in the red flour beetle, Tribolium castaneum

BACKGROUND

Insect cytochrome P450 monooxygenases (P450s) play a crucial role in phytochemical metabolism and tolerance. Three P450 genes (TcCYP4Q3, TcCYP4Q5, and TcCYP4Q7) are associated with the response of eugenol in Tribolium castaneum. However, the responding mechanisms of these P450 genes to eugenol remain unknown.

RESULTS

Here, spatiotemporal expression profiling revealed that TcCYP4Q3 and TcCYP4Q5 were most highly expressed in late adult, while TcCYP4Q7 was predominantly expressed in late larva; and all of these three P450 genes were mainly expressed in the fat body of larvae. Furthermore, the expressions of these three P450 genes were significantly up-regulated after exposure to eugenol, and depletion of them enhanced the susceptibility of beetles to eugenol. Interestingly, RNA interference (RNAi) against the CncC gene, a transcription factor of CncC signaling pathway associated with regulation of insect P450s in response to phytochemicals, reduced the transcripts of these three P450 genes following exposure to eugenol. Investigation of CncC signaling pathway showed that this pathway could be activated by eugenol.

CONCLUSION

Altogether, the results indicate that these three P450 genes are regulated by CncC signaling pathway to participate in the susceptibility of Tribolium castaneum to phytochemicals. These findings will aid implications for the development of novel therapeutics to control pest. © 2022 Society of Chemical Industry.

A chromosome-level genome assembly of the orange wheat blossom midge, Sitodiplosis mosellana Géhin (Diptera: Cecidomyiidae) provides insights into the evolution of a detoxification system

The orange wheat blossom midge Sitodiplosis mosellana Géhin (Diptera: Cecidomyiidae), an economically important pest, has caused serious yield losses in most wheat-growing areas worldwide in the past half-century. A high-quality chromosome-level genome for S. mosellana was assembled using PacBio long read, Illumina short read, and Hi-C sequencing technologies. The final genome assembly was 180.69 Mb, with contig and scaffold N50 sizes of 998.71 kb and 44.56 Mb, respectively. Hi-C scaffolding reliably anchored 4 pseudochromosomes, accounting for 99.67% of the assembled genome. In total, 12,269 protein-coding genes were predicted, of which 91% were functionally annotated. Phylogenetic analysis indicated that S. mosellana and its close relative, the swede midge Contarinia nasturtii, diverged about 32.7 MYA. The S. mosellana genome showed high chromosomal synteny with the genome of Drosophila melanogaster and Anopheles gambiae. The key gene families involved in the detoxification of plant secondary chemistry were analyzed. The high-quality S. mosellana genome data will provide an invaluable resource for research in a broad range of areas, including the biology, ecology, genetics, and evolution of midges, as well as insect–plant interactions and coevolution.

The phylogenetic and evolutionary analyses of detoxification gene families in Aphidinae species

Detoxification enzymes play significant roles in the interactions between insects and host plants, wherein detoxification-related genes make great contributions. As herbivorous pests, aphids reproduce rapidly due to parthenogenesis. They are good biological materials for studying the mechanisms that allow insect adaptation to host plants. Insect detoxification gene families are associated with insect adaptation to host plants. The Aphidinae is the largest subfamily in the Aphididae with at least 2483 species in 256 genera in 2 tribes: the Macrosiphini (with 3/4 of the species) and the Aphidini. Most aphid pests on crops and ornamental plants are Aphidinae. Members of the Aphidinae occur in nearly every region of the world. The body shape and colour vary significantly. To research the role that detoxification gene families played in the process of aphid adaptation to host evolution, we analyzed the phylogeny and evolution of these detoxification gene families in Aphidinae. In general, the P450/GST/CCE gene families contract, whereas the ABC/UGT families are conserved in Aphidinae species compared to these families in other herbivorous insects. Genus-specific expansions of P450 CYP4, and GST Delta have occurred in the genus Acyrthosiphon. In addition, the evolutionary rates of five detoxification gene families in the evolution process of Aphidinae are different. The comparison of five detoxification gene families among nine Aphidinae species and the estimated relative evolutionary rates provided herein support an understanding of the interaction between and the co-evolution of Aphidinae and plants.

The Genome of Rhyzopertha dominica (Fab.) (Coleoptera: Bostrichidae): Adaptation for Success

The lesser grain borer, Rhyzopertha dominica (F.) (Coleoptera: Bostrichidae), is a major global pest of cereal grains. Infestations are difficult to control as larvae feed inside grain kernels, and many populations are resistant to both contact insecticides and fumigants. We sequenced the genome of R. dominica to identify genes responsible for important biological functions and develop more targeted and efficacious management strategies. The genome was assembled from long read sequencing and long-range scaffolding technologies. The genome assembly is 479.1 Mb, close to the predicted genome size of 480.4 Mb by flow cytometry. This assembly is among the most contiguous beetle assemblies published to date, with 139 scaffolds, an N₅₀ of 53.6 Mb, and L₅₀ of 4, indicating chromosome-scale scaffolds. Predicted genes from biologically relevant groups were manually annotated using transcriptome data from adults and different larval tissues to guide annotation. The expansion of carbohydrase and serine peptidase genes suggest that they combine to enable efficient digestion of cereal proteins. A reduction in the copy number of several detoxification gene families relative to other coleopterans may reflect the low selective pressure on these genes in an insect that spends most of its life feeding internally. Chemoreceptor genes contain elevated numbers of pseudogenes for odorant receptors that also may be related to the recent ontogenetic shift of R. dominica to a diet consisting primarily of stored grains. Analysis of repetitive sequences will further define the evolution of bostrichid beetles compared to other species. The data overall contribute significantly to coleopteran genetic research.

The Birth-and-Death Evolution of Cytochrome P450 Genes in Bees

The birth-and-death model of multigene family evolution describes how gene families evolve and diversify through duplication and deletion. The cytochrome P450s are one of the most diverse and well-studied multigene families, involved in both physiological and xenobiotic functions. Extensive studies of insect P450 genes have demonstrated their role in insecticide resistance. Bees are thought to experience toxin exposure through their diet of nectar and pollen, as well as the resin-collecting behavior exhibited by some species. Here, we describe the repertoire of P450 genes in the orchid bee Euglossa dilemma. Male orchid bees form perfume bouquets used in courtship displays by collecting volatile compounds, resulting in exposure to compounds known to be toxic. In addition, we conducted phylogenetic and selection analyses across ten bee species encompassing three bee families. We find that social behavior and resin collection are not correlated with the repertoire of P450 present in a bee species. However, our analyses revealed that P450 clades can be classified as stable and unstable, and that genes involved in xenobiotic metabolism are more likely to belong to unstable clades. Furthermore, we find that unstable clades are under more dynamic evolutionary pressures and exhibit signals of adaptive evolution. This work highlights the complexity of multigene family evolution, revealing that multiple factors contribute to the diversification, stability, and dynamics of this gene family. Furthermore, we provide a resource for future detailed studies investigating the function of different P450s in economically important bee species.

Full-length transcriptome analysis of Spodoptera frugiperda larval brain reveals detoxification genes

Background

Spodoptera frugiperda (J. E. Smith), commonly known as fall armyworm (FAW), is one of the most destructive agricultural pests in the world and has posed a great threat to crops. The improper use of insecticides has led to rapid development of resistance. However, the genetic data available for uncovering the insecticide resistance mechanisms are scarce.

Methods

In this study, we used PacBio single-molecule real-time (SMRT) sequencing aimed at revealing the full-length transcriptome profiling of the FAW larval brain to obtain detoxification genes.

Results

A total of 18,642 high-quality transcripts were obtained with an average length of 2,371 bp, and 11,230 of which were successfully annotated in six public databases. Among these, 5,692 alternative splicing events were identified.

Cytochrome P450s Are Essential for Insecticide Tolerance in the Endoparasitoid Wasp Meteorus pulchricornis (Hymenoptera: Braconidae)

With the widespread application of insecticides, parasitoid wasps may also be under risk when exposed to insecticides directly at their free-living stages. The endoparasitoid wasp Meteorus pulchricornis is the predominant natural enemy of many lepidopteran pests, such as Spodoptera litura and Helicoverpa armigera. The cytochrome P450 monooxygenases constitute a ubiquitous and complex superfamily of hydrophobic, haem-containing enzymes. P450s are involved in the detoxification of many xenobiotics. However, their exact roles in the tolerance mechanism in parasitoids toward insecticides has received less attention. Here, 28 P450 genes in M. pulchricornis were identified from a previously constructed transcriptome dataset. These P450 genes belonged to CYP2, -3, and -4, and mitochondrial clans. Subsequently, eight candidate MpulCYPs were selected from four CYP clans to validate their expression patterns under phoxim, cypermethrin, and chlorfenapyr exposure by qRT-PCR. The results showed that all three insecticides had significant effects on the expression of MpulCYPs. To further study the function of P450s, CYP369B3 was silenced, and its expression levels of CYP369B3 were significantly decreased. Survival analysis indicated that after dsRNA injection, the mortality rate of wasps was significantly increased when M. pulchricornis females were exposed to insecticides compared to control groups. Our findings provide a theoretical base for elucidating the mechanism of insecticide tolerance and promote functional research on P450 genes in parasitoid wasps.

Transcriptome analysis unravels RNAi pathways genes and putative expansion of CYP450 gene family in cotton leafhopper Amrasca biguttula (Ishida)

Cotton Leafhopper, Amrasca biguttula is an important pest of cotton and okra in the Indian subcontinent. Presently limited genomic/transcriptomic information is available for this insect in any of open source databases. The present study reports the first assembled and annotated de novo transcriptome of cotton leafhopper. Out of 75,551 transcripts, 39,613 CDS (Coding Sequence) were predicted with 35,282 showing positive blast hits with NCBI nr database. The Gene ontology (GO) analysis annotated 7431 CDS  with KEGG pathway categorizing these CDS into 22 different functional groups. The majority of CDS were annotated in signal transduction and transport catabolism pathways. The sequence data was screened for RNAi pathway genes and presence of 37 transcripts associated with this process confirmed the existence of robust RNAi machinery. The role of core RNAi machinery genes (Dicer-2, Ago-2, Piwi and Staufen) has been validated through dsRNA feeding studies. The data resource has also been used to identify potential RNAi targets and genes associated with insecticide detoxification specifically CYP 450 family. The current study provides a useful sequence resource which can be used to initiate molecular studies in this insect with emphasis on insecticide resistance, RNAi and functional genomics.

Comparative analyses of transcriptional responses of Dectes texanus LeConte (Coleoptera: Cerambycidae) larvae fed on three different host plants and artificial diet

Dectes texanus is an important coleopteran pest of soybeans and cultivated sunflowers in the Midwestern United States that causes yield losses by girdling stems of their host plants. Although sunflower and giant ragweed are primary hosts of D. texanus, they began colonizing soybeans approximately 50 years ago and no reliable management method has been established to prevent or reduce losses by this pest. To identify genes putatively involved when feeding soybean, we compared gene expression of D. texanus third-instar larvae fed soybean to those fed sunflower, giant ragweed, or artificial diet. Dectes texanus larvae differentially expressed 514 unigenes when fed on soybean compared to those fed the other diet treatments. Enrichment analyses of gene ontology terms from up-regulated unigenes in soybean-fed larvae compared to those fed both primary hosts highlighted unigenes involved in oxidoreductase and polygalacturonase activities. Cytochrome P450s, carboxylesterases, major facilitator superfamily transporters, lipocalins, apolipoproteins, glycoside hydrolases 1 and 28, and lytic monooxygenases were among the most commonly up-regulated unigenes in soybean-fed larvae compared to those fed their primary hosts. These results suggest that D. texanus larvae differentially expressed unigenes involved in biotransformation of allelochemicals, digestion of plant cell walls and transport of small solutes and lipids when feeding in soybean.

The genome of the European estuarine calanoid copepod Eurytemora affinis: Potential use in molecular ecotoxicology

In this study, we sequenced and assembled the genome of a European estuarine calanoid copepod using Oxford Nanopore PromethION and Illumina HiSeq 2500 platforms. The length of the assembled genome was 776.1 Mb with N50 = 474.9 kb (BUSCO 85.9%), and the genome consisted of 2473 contigs. A total of 18,014 genes were annotated and orthologous gene clusters were analyzed in comparison to other copepods. In addition, genome-wide identification of cytochrome P450s, glutathione S-transferases, and ATP-binding cassette transporters in E. affinis was performed to determine gene repertoire of these detoxification-related gene families. Results revealed the presence of species-specific gene inventories, indicating that these gene families have evolved through species-specific gene loss/expansion processes, possibly due to adaptation to different environmental stressors. Our study provides a new inventory of the European estuarine calanoid copepod E. affinis genome with emphasis on phase I, II, and III detoxification systems.

Full length sequencing reveals novel transcripts of detoxification genes along with related alternative splicing events and lncRNAs in Phyllotreta striolata

The striped flea beetle, Phyllotreta striolata (Fabricius), damages crops in the Brassicaceae. The genetic data for this pest are insufficient to reveal its insecticide resistance mechanisms or to develop molecular markers for resistance monitoring. We used PacBio Iso-Seq technology to sequence the full-length transcriptome of P. striolata. After isoform sequence clustering and removal of redundant transcripts, a total of 41,293 transcripts were obtained, and 35,640 of these were annotated in the database of gene products. Structure analysis uncovered 4,307 alternative splicing events, and 3,836 sequences were recognized as lncRNAs. Transcripts with the complete coding region of important detoxification enzymes were further classified. There were 57 transcripts of P450s distributed in CYP2, CYP3, CYP4, and Mito CYP clades, 29 transcripts of ESTs from 4 functional groups, 17 transcripts of GSTs classified into 5 families, 51 transcripts of ABCs distributed in 6 families, and 19 transcripts of UGTs. Twenty-five lncRNAs were predicted to be regulators of these detoxification genes. Full-length transcriptome sequencing is an efficient method for molecular study of P. striolata and it is also useful for gene function analysis.

A chromosome-level assembly of the harlequin ladybird Harmonia axyridis as a genomic resource to study beetle and invasion biology

The harlequin ladybird, Harmonia axyridis (Pallas), is a well-known model organism for genetic studies and is also a well-studied natural enemy used for pest control. It became an invasive species after being introduced to North America and Europe as a pest control agent. Though two genome assemblies for this insect have been previously reported, a high-quality genome assembly at the chromosome level is still not available. Here, we obtained a new chromosome-level genome assembly of H. axyridis by combining various sequencing technologies, namely Illumina short reads, PacBio long reads, 10X Genomics and Hi-C. The chromosome-level genome assembly is 423 Mb with a scaffold N50 of 45.92 Mb. Using Hi-C data 1,897 scaffolds were anchored to eight chromosomes. A total of 730,068 repeat sequences were identified, making up 51.2% of the assembled genome. After masking these repeat sequences, we annotated 22,810 protein-encoding genes. The X chromosome and Y-linked scaffolds were also identified by resequencing male and female genomes and calculating the male to female coverage ratios. Two gene families associated with environmental adaptation, odorant receptor and cytochrome P450, were analysed and showed no obvious expansion in H. axyridis. We successfully constructed a putative biosynthesis pathway of harmonine, a defence compound in the haemolymph of H. axyridis, which is a key factor for H. axyridis strong immunity. The chromosome-level genome assembly of H. axyridis is a helpful resource for studies of beetle biology and invasive biology.

The role of cytochrome P450-mediated detoxification in insect adaptation to xenobiotics

Insect cytochrome P450 monooxygenases (P450s) are well known to be involved in metabolic detoxification of xenobiotics, such as phytochemicals, insecticides and environmental pollutants. Enhanced metabolic detoxification is closely associated with the constitutive overexpression and induction of P450s. In general, multiple insect P450s are co-responsible for xenobiotic detoxification. Considering the capacity of P450s to respond to a wide range of xenobiotics, synergistic interactions between natural and synthetic xenobiotics and P450-mediated cross-tolerance/resistance are ubiquitous. Recent studies have indicated that both transcription factors and signaling pathways are involved in the regulation of P450 genes in xenobiotic responses. This article reviews our current understanding of P450-mediated detoxification in insect adaptation to xenobiotics and highlights recent progress in the molecular basis of P450 regulation.

Transcriptomic analysis of s-methoprene resistance in the lesser grain borer, Rhyzopertha dominica, and evaluation of piperonyl butoxide as a resistance breaker

BACKGROUND

The lesser grain borer, Rhyzopertha dominica is a serious pest of stored grains. Fumigation and contact insecticides play a major role in managing this pest globally. While insects are developing genetic resistance to chemicals, hormonal analogues such as s-methoprene play a key role in reducing general pest pressure as well as managing pest populations that are resistant to fumigants and neurotoxic contact insecticides. However, resistance to s-methoprene has been reported in R. dominica with some reports showing a remarkable high resistance, questioning the use of this compound and other related analogues in grain protection. The current study attempts to identify possible molecular mechanisms that contribute in resistance to s-methoprene in R. dominica.

RESULTS

Transcriptome analysis of resistant and susceptible strains of this pest species identified a set of differentially expressed genes related to cytochrome P450s, indicating their potential role in resistance to s-methoprene. Laboratory bioassays were performed with s-methoprene treated wheat grains in presence and absence of piperonyl butoxide (PBO), a cytochrome P450 inhibitor. The results indicate that PBO, when applied alone, at least at the concentration tested here, had no effect on R. dominica adult emergence, but has a clear synergistic effect to s-methoprene. The number of produced progeny decreased in presence of the inhibitor, especially in the resistant strain. In addition, we also identified CYP complement (CYPome) of R. dominica, annotated and analysed phylogenetically, to understand the evolutionary relationships with other species.

CONCLUSIONS

The information generated in current study suggest that PBO can effectively be used to break resistance to s-methoprene in R. dominica.

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

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

Off-target effects of RNAi correlate with the mismatch rate between dsRNA and non-target mRNA

RNAi is a potent technique for the knockdown of target genes. However, its potential off-target effects limit the widespread applications in both reverse genetic analysis and genetic manipulation. Previous efforts have uncovered rules underlying specificity of siRNA-based silencing, which has broad applications in humans, but the basis for specificity of dsRNAs, which are better suited for use as insecticides, is poorly understood. Here, we investigated the rules governing dsRNA specificity. Mutational analyses showed that dsRNAs with >80% sequence identity with target genes triggered RNAi efficiently. dsRNAs with ≥16 bp segments of perfectly matched sequence or >26 bp segments of almost perfectly matched sequence with one or two mismatches scarcely distributed (single mismatches inserted between ≥5 bp matching segments or mismatched couplets inserted between ≥8 bp matching segments) also able to trigger RNAi. Using these parameters to predict off-target risk, dsRNAs can be designed to optimize specificity and efficiency, paving the way to the widespread, rational application of RNAi in pest control.

Genomic and transcriptomic analysis unveils population evolution and development of pesticide resistance in fall armyworm Spodoptera frugiperda

The fall armyworm (FAW), Spodoptera frugiperda, is a destructive pest native to America and has recently become an invasive insect pest in China. Because of its rapid spread and great risks in China, understanding of FAW genetic background and pesticide resistance is urgent and essential to develop effective management strategies. Here, we assembled a chromosome-level genome of a male FAW (SFynMstLFR) and compared re-sequencing results of the populations from America, Africa, and China. Strain identification of 163 individuals collected from America, Africa and China showed that both C and R strains were found in the American populations, while only C strain was found in the Chinese and African populations. Moreover, population genomics analysis showed that populations from Africa and China have close relationship with significantly genetic differentiation from American populations. Taken together, FAWs invaded into China were most likely originated from Africa. Comparative genomics analysis displayed that the cytochrome p450 gene family is extremely expanded to 425 members in FAW, of which 283 genes are specific to FAW. Treatments of Chinese populations with twenty-three pesticides showed the variant patterns of transcriptome profiles, and several detoxification genes such as AOX, UGT and GST specially responded to the pesticides. These findings will be useful in developing effective strategies for management of FAW in China and other invaded areas.

Genome-enabled insights into the biology of thrips as crop pests

BACKGROUND

The western flower thrips, Frankliniella occidentalis (Pergande), is a globally invasive pest and plant virus vector on a wide array of food, fiber, and ornamental crops. The underlying genetic mechanisms of the processes governing thrips pest and vector biology, feeding behaviors, ecology, and insecticide resistance are largely unknown. To address this gap, we present the F. occidentalis draft genome assembly and official gene set.

RESULTS

We report on the first genome sequence for any member of the insect order Thysanoptera. Benchmarking Universal Single-Copy Ortholog (BUSCO) assessments of the genome assembly (size = 415.8 Mb, scaffold N50 = 948.9 kb) revealed a relatively complete and well-annotated assembly in comparison to other insect genomes. The genome is unusually GC-rich (50%) compared to other insect genomes to date. The official gene set (OGS v1.0) contains 16,859 genes, of which ~ 10% were manually verified and corrected by our consortium. We focused on manual annotation, phylogenetic, and expression evidence analyses for gene sets centered on primary themes in the life histories and activities of plant-colonizing insects. Highlights include the following: (1) divergent clades and large expansions in genes associated with environmental sensing (chemosensory receptors) and detoxification (CYP4, CYP6, and CCE enzymes) of substances encountered in agricultural environments; (2) a comprehensive set of salivary gland genes supported by enriched expression; (3) apparent absence of members of the IMD innate immune defense pathway; and (4) developmental- and sex-specific expression analyses of genes associated with progression from larvae to adulthood through neometaboly, a distinct form of maturation differing from either incomplete or complete metamorphosis in the Insecta.

CONCLUSIONS

Analysis of the F. occidentalis genome offers insights into the polyphagous behavior of this insect pest that finds, colonizes, and survives on a widely diverse array of plants. The genomic resources presented here enable a more complete analysis of insect evolution and biology, providing a missing taxon for contemporary insect genomics-based analyses. Our study also offers a genomic benchmark for molecular and evolutionary investigations of other Thysanoptera species.

Deciphering the molecular mechanisms of mother-to-egg immune protection in the mealworm beetle Tenebrio molitor

In a number of species, individuals exposed to pathogens can mount an immune response and transmit this immunological experience to their offspring, thereby protecting them against persistent threats. Such vertical transfer of immunity, named trans-generational immune priming (TGIP), has been described in both vertebrates and invertebrates. Although increasingly studied during the last decade, the mechanisms underlying TGIP in invertebrates are still elusive, especially those protecting the earliest offspring life stage, i.e. the embryo developing in the egg. In the present study, we combined different proteomic and transcriptomic approaches to determine whether mothers transfer a "signal" (such as fragments of infecting bacteria), mRNA and/or protein/peptide effectors to protect their eggs against two natural bacterial pathogens, namely the Gram-positive Bacillus thuringiensis and the Gram-negative Serratia entomophila. By taking the mealworm beetle Tenebrio molitor as a biological model, our results suggest that eggs are mainly protected by an active direct transfer of a restricted number of immune proteins and of antimicrobial peptides. In contrast, the present data do not support the involvement of mRNA transfer while the transmission of a "signal", if it happens, is marginal and only occurs within 24h after maternal exposure to bacteria. This work exemplifies how combining global approaches helps to disentangle the different scenarios of a complex trait, providing a comprehensive characterization of TGIP mechanisms in T. molitor. It also paves the way for future alike studies focusing on TGIP in a wide range of invertebrates and vertebrates to identify additional candidates that could be specific to TGIP and to investigate whether the TGIP mechanisms found herein are specific or common to all insect species.

Genome-wide and expression-profiling analyses of the cytochrome P450 genes in Bactrocera dorsalis (Hendel) and screening of candidate P450 genes associated with malathion resistance

BACKGROUND

Bactrocera dorsalis (Hendel) is a notorious agricultural pest worldwide, and its resistance to insecticides is a major obstacle in successful control. Cytochrome P450s (P450s) are major metabolic enzymes associated with insecticide resistance. The genome of B. dorsalis was sequenced recently, allowing an integrated genome-wide analysis of P450 genes (P450s) and the analysis of correlations between these genes and insecticide resistance in this pest.

RESULTS

Totally, 101 P450s were identified in the B. dorsalis genome and classified into four clans, 25 families and 57 subfamilies. Quantitative reverse transcription polymerase chain reaction results showed that most of these genes were highly expressed in adults (46) and in metabolic tissues, including the fatbody (63), midgut (61) and Malphagian tubules (66). In a malathion-resistant strain, 13 and 9 genes were significantly upregulated and downregulated, respectively, compared with a susceptible strain, and these genes were screened as candidate genes associated with malathion resistance.

CONCLUSION

This study provides useful information for understanding the evolution and potential functions of P450s in B. dorsalis, and the results lay the foundation for further studies on the correlations between P450s and malathion resistance in B. dorsalis. © 2020 Society of Chemical Industry.

Genome mining and functional analysis of cytochrome P450 genes involved in insecticide resistance in Leucinodes orbonalis (Lepidoptera: Crambidae)

Genome-wide analysis of cytochrome P450 monooxygenase (CYP) genes from the advanced genome project of the Leucinodes orbonalis and the expression analysis provided significant information about the metabolism-mediated insecticide resistance. A total of 72 putative CYP genes were identified from the genome and transcriptome of L. orbonalis. The genes were classified under 30 families and 46 subfamilies based on the standard nomenclature. In the present study, a novel CYP gene, CYP324F1, was identified and it has not been reported from any other living system so far. Biochemical assays showed enhanced titers (5.81-18.5-fold) of O-demethylase of CYP in five field-collected populations. We selected 34 homologous CYP gene sequences, seemed to be involved in insecticide resistance for primer design and quantitative real-time PCR studies. Among the many overexpressed genes (>10 fold), the expression levels of CYP324F1 and CYP306A1 were prominent across all the field populations as compared with the susceptible iso-female line. Oral delivery of ds-CYP324F1 and ds-CYP306A1 directed against CYP324F1 and CYP306A1 to the larvae of one of the insecticide resistance populations caused reduced expression of these two transcripts in a dose-dependent manner (53.4%-85.0%). It appears that the increased titer of O-demethylase is the result of increased transcription level of CYP genes in resistant populations. The data provide insight for identifying the novel resistance management strategies against L. orbonalis.

CncC/Maf-mediated xenobiotic response pathway in insects

Insects have evolved resistance to almost all insecticides developed for their control. Multiple mechanisms of resistance, including enhanced metabolism and excretion of insecticides, target-site insensitivity, reduced penetration of insecticides, and avoidance behavior, have been reported. The genes coding for proteins involved in resistance have been identified in numerous insects. The enzymes and transporters required for all three phases of insecticide metabolism and excretion including cytochrome P450 monooxygenases, glutathione S-transferases, UDP-glucuronosyltransferases, carboxylesterases, and ATP-binding cassette transmembrane transporters have been identified. Recent research in multiple insect species identified CNC-bZIP transcription factor superfamily members as regulators of genes coding for enzymes and transporters involved in insecticide metabolic resistance. The information on the pathway including reactive oxygen species, cap "n" collar isoform-C, and its heterodimer partner, muscle aponeurosis fibromatosis transcription factors involved in overexpression of enzymes and transporters involved insecticide resistance will be summarized.

Tribolium castaneum: A Model for Investigating the Mode of Action of Insecticides and Mechanisms of Resistance

The red flour beetle, Tribolium castaneum, is a worldwide insect pest of stored products, particularly food grains, and a powerful model organism for developmental, physiological and applied entomological research on coleopteran species. Among coleopterans, T. castaneum has the most fully sequenced and annotated genome and consequently provides the most advanced genetic model of a coleopteran pest. The beetle is also easy to culture and has a short generation time. Research on this beetle is further assisted by the availability of expressed sequence tags and transcriptomic data. Most importantly, it exhibits a very robust response to systemic RNA interference (RNAi), and a database of RNAi phenotypes (iBeetle) is available. Finally, classical transposonbased techniques together with CRISPR/Cas-mediated gene knockout and genome editing allow the creation of transgenic lines. As T. castaneum develops resistance rapidly to many classes of insecticides including organophosphates, methyl carbamates, pyrethroids, neonicotinoids and insect growth regulators such as chitin synthesis inhibitors, it is further a suitable test system for studying resistance mechanisms. In this review, we will summarize recent advances in research focusing on the mode of action of insecticides and mechanisms of resistance identified using T. castaneum as a pest model.

Identification and characterization of detoxification genes in two cerambycid beetles, Rhaphuma horsfieldi and Xylotrechus quadripes (Coleoptera: Cerambycidae: Clytini)

The longhorned beetles, Rhaphuma horsfieldi and Xylotrechus quadripes, are two polyphagous insects with larvae feeding on different host plants. In this study, we identified and characterized three gene superfamilies of cytochrome P450s (CYPs), carboxylesterases (COEs) and glutathione-S-transferases (GSTs) involved in the detoxification of endobiotics (e.g., hormones and steroids) and xenobiotics (e.g., insecticides, sex pheromones and plant allelochemicals) through a combination approach of bioinformatics, phylogenetics, expression profiles and genomics. Transcriptome analyses led to the identification of 281 transcripts encoding 135 P450s, 108 COEs and 38 GSTs from the two beetles, coupled with comparative studies of detoxification genes among coleopteran species, suggesting a correlation between host range and the sizes of P450 or COE gene repertoires. The P450s of two beetles were phylogenetically classified into four clades, representing the majority of genes in the CYP3 clan. The COEs from R. horsfieldi and X. quadripes were separately grouped into 11 and 10 clades, and the GST superfamily was assigned into six clades. Expression profiles revealed that the detoxification genes were broadly expressed in various tissues as an implication of functional diversities. Ultimately and more importantly, five alternative splicing events in the Epsilon GSTs, including RhorGSTe7.1/GSTe7.2 and XquaGSTe3.1/GST3.2, were acquired in Coleoptera, in which these genes and their orthologs shared highly conserved gene structure. Our current study has complemented the resources for the detoxification genes in the family Cerambycidae, and allows for functional experiments to identify candidate molecular targets involved in pest resistance to insecticides like organophosphates, organochlorines and pyrethroids.

Physiological resistance alters behavioral response of Tetranychus urticae to acaricides

Multiple acaricide resistance in Tetranychus urticae continues to threaten crop production globally, justifying the need to adequately study resistance for sustainable pest management. Most studies on acaricide resistance have focused on the acute contact toxicity of acaricides with little or no information on the behavioral responses elicited after acaricide exposure. Furthermore, the impact of physiological resistance on these behavioral responses remains unknown in most pest species, including T. urticae. We tested the effect of acaricide resistance on contact toxicity, irritancy and repellency of mitochondrial electron transport inhibitor of complex I (MET-I) and mite growth inhibitor (MGI) acaricides on multiple T. urticae strains. We also tested whether acaricides with similar physiological target site/mode of action also elicit similar behavioral effects on T. urticae strains. MET-I acaricides (fenazaquin, fenpyroximate, and pyrabiden) and MGIs (clofentezine, hexythiazox and etoxazole) elicited a dose-dependent irritant and repellent effect on T. urticae. Selection of strains for physiological resistance to these acaricides affected the behavioral response of T. urticae, especially in MET-I resistant strains, that showed reduced irritancy and repellency to MET-I acaricides. Behavioral response also affected the oviposition of T. urticae, where strains generally showed preferential oviposition away from the acaricides. The outcome of this study highlights negative consequences of acaricide resistance that can potentially affect T. urticae management.

Transcriptome-Based Identification and Molecular Evolution of the Cytochrome P450 Genes and Expression Profiling under Dimethoate Treatment in Amur Stickleback (Pungitius sinensis)

Cytochrome P450s (CYPs) are a family of membrane-bound mono-oxygenase proteins, which are involved in cell metabolism and detoxification of various xenobiotic substances. In this study, we identified 58 putative CYP genes in Amur stickleback (Pungitius sinensis) based on the transcriptome sequencing. Conserved motif distribution suggested their functional relevance within each group. Some present recombination events have accelerated the evolution of this gene family. Moreover, a few positive selection sites were identified, which may have accelerated the functional divergence of this family of proteins. Expression patterns of these CYP genes were investigated and indicated that most were affected by dimethoate treatment, suggesting that CYPs were involved in the detoxication of dimethoate. This study will provide a foundation for the further functional investigation of CYP genes in fishes.