Comparative analysis of the UDP-glycosyltransferase multigene family in insects

Characterization of the UDP-glycosyltransferase UGT33D1 in silkworm Bombyx mori

Uridine diphosphate (UDP)-glycosyltransferases (UGTs) are important metabolizing enzymes functioning by adding a sugar moiety to a small lipophilic substrate molecule and play critical roles in drug/toxin metabolism for all realms of life. In this study, the silkworm Bombyx mori UGT33D1 gene was characterized in detail. UGT33D1 was found localized in the endoplasmic reticulum (ER) compartment just like other animal UGTs and was mainly expressed in the silkworm midgut. We first reported that UGT33D1 was important to BmNPV infection, as silencing UGT33D1 inhibited the BmNPV infection in silkworm BmN cells, while overexpressing the gene promoted viral infection. The molecular pathways regulated by UGT33D1 were analysed via transcriptome sequencing upon UGT33D1 knockdown, highlighting the important role of the gene in maintaining a balanced oxidoreductive state of the organism. In addition, proteins that physically interact with UGT33D1 were identified through immunoprecipitation and mass spectrometry analysis, which includes tubulin, elongation factor, certain ribosomal proteins, histone proteins and zinc finger proteins that had been previously reported for human UGT-interacting proteins. This study provided preliminary but important functional information on UGT33D1 and is hoped to trigger deeper investigations into silkworm UGTs and their functional mechanisms.

Disarming the defenses: Insect detoxification of plant defense-related specialized metabolites

The ability of certain insects to feed on plants containing toxic specialized metabolites may be attributed to detoxification enzymes. Representatives of a few large families of detoxification enzymes are widespread in insect herbivores acting to functionalize toxins and conjugate them with polar substituents to decrease toxicity, increase water solubility and enhance excretion. Insects have also developed specific enzymes for coping with toxins that are activated upon plant damage. Another source of detoxification potential in insects lies in their microbiomes, which are being increasingly recognized for their role in processing plant toxins. The evolution of insect detoxification systems to resist toxic specialized metabolites in plants may in turn have selected for the great diversity of such metabolites found in nature.

A UDP-glycosyltransferase gene PcUGT202A9 was associated with abamectin resistance in Panonychus citri (McGregor)

Panonychus citri (McGregor) strains have developed a high level of resistance to abamectin, but the underlying molecular mechanism is unknown. Uridine diphosphate (UDP)-glycosyltransferases (UGTs) are critical for the removal of a variety of exogenous and endogenous substances. In this study, an enzyme activity assay revealed that UGTs potentially contribute to P. citri abamectin resistance. Spatiotemporal expression profiles showed that only PcUGT202A9 was significantly overexpressed in the abamectin-resistant strain (AbR) at all developmental stages. Moreover, UGT activity decreased significantly, whereas abamectin susceptibility increased significantly, in AbR after PcUGT202A9 was silenced. Three-dimensional modeling and molecular docking analyses revealed that PcUGT202A9 can bind stably to abamectin. Recombinant PcUGT202A9 activity was detected when α-naphthol was used, but the enzymatic activity was inhibited by abamectin (50 % inhibitory concentration: 803.3 ± 14.20 μmol/L). High-performance liquid chromatography and mass spectrometry analyses indicated that recombinant PcUGT202A9 can effectively degrade abamectin and catalyze the conjugation of UDP-glucose to abamectin. These results imply PcUGT202A9 contributes to P. citri abamectin resistance.

Duplications and Losses of the Detoxification Enzyme Glycosyltransferase 1 Are Related to Insect Adaptations to Plant Feeding

Insects have developed sophisticated detoxification systems to protect them from plant secondary metabolites while feeding on plants to obtain necessary nutrients. As an important enzyme in the system, glycosyltransferase 1 (GT1) conjugates toxic compounds to mitigate their harm to insects. However, the evolutionary link between GT1s and insect plant feeding remains elusive. In this study, we explored the evolution of GT1s across different insect orders and feeding niches using publicly available insect genomes. GT1 is widely present in insect species; however, its gene number differs among insect orders. Notably, plant-sap-feeding species have the highest GT1 gene numbers, whereas blood-feeding species display the lowest. GT1s appear to be associated with insect adaptations to different plant substrates in different orders, while the shift to non-plant feeding is related to several losses of GT1s. Most large gene numbers are likely the consequence of tandem duplications showing variations in collinearity among insect orders. These results reveal the potential relationships between the evolution of GT1s and insect adaptation to plant feeding, facilitating our understanding of the molecular mechanisms underlying insect-plant interactions.

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.

CRISPR/Cas9-Based Functional Characterization of SfUGT50A15 Reveals Its Roles in the Resistance of Spodoptera frugiperda to Chlorantraniliprole, Emamectin Benzoate, and Benzoxazinoids

UDP-glycosyltransferases (UGTs) are a diverse superfamily of enzymes. Insects utilize uridine diphosphate-glucose (UDP-glucose) as a glycosyl donor for glycosylation in vivo, involved in the glycosylation of lipophilic endosymbionts and xenobiotics, including phytotoxins. UGTs act as second-stage detoxification metabolizing enzymes, which are essential for the detoxification metabolism of insecticides and benzoxazine compounds. However, the UGT genes responsible for specific glycosylation functions in S. frugiperda are unclear at present. In this study, we utilized CRISPR/Cas9 to produce a SfUGT50A15-KO strain to explore its possible function in governing sensitivity to chemical insecticides or benzoxazinoids. The bioassay results suggested that the SfUGT50A15-KO strain was significantly more sensitive to chlorantraniliprole, emamectin benzoate, and benzoxazinoids than the wild-type strains. This finding suggests that the overexpression of the SfUGT50A15 gene may be linked to S. frugiperda resistance to pesticides (chlorantraniliprole and emamectin benzoate) as well as benzoxazinoids (BXDs).

Comparative genomics provides insights into molecular adaptation to hypermetamorphosis and cantharidin metabolism in blister beetles (Coleoptera: Meloidae)

Blister beetles (Coleoptera: Meloidae) are currently subdivided into three subfamilies: Eleticinae (a basal group), Nemognathinae, and Meloinae. These are all characterized by the endogenous production of the defensive terpene cantharidin (CA), whereas the two most derived subfamilies show a hypermetamorphic larval development. Here, we provide novel draft genome assemblies of five species sampled across the three blister beetle subfamilies (Iselma pallidipennis, Stenodera caucasica, Zonitis immaculata, Lydus trimaculatus, and Mylabris variabilis) and performed a comparative analysis with other available Meloidae genomes and the closely-related canthariphilous species (Pyrochroa serraticornis) to disclose adaptations at a molecular level. Our results highlighted the expansion and selection of genes potentially responsible for CA production and metabolism, as well as its mobilization and vesicular compartmentalization. Furthermore, we observed adaptive selection patterns and gain of genes devoted to epigenetic regulation, development, and morphogenesis, possibly related to hypermetamorphosis. We hypothesize that most genetic adaptations occurred to support both CA biosynthesis and hypermetamorphosis, two crucial aspects of Meloidae biology that likely contributed to their evolutionary success.

Uridine diphosphate glucosyltransferases are involved in spinosad resistance in western flower thrips Frankliniella occidentalis (Pergande)

Uridine diphosphate glucosyltransferases (UGTs) play crucial roles in the insect detoxification system and are associated with pesticide resistance. Our previous transcriptomic analysis of spinosad-susceptible (Ivf03) and resistant (NIL-R) Frankliniella occidentalis revealed numerous upregulated UGT genes in the NIL-R strain, suggesting their potential contribution to spinosad resistance. To investigate this hypothesis, here we conducted UGT activity assays and spinosad induction experiments, employing RNA interference (RNAi) techniques for gene function validation. We found significantly elevated UGT activity in the NIL-R strain compared to Ivf03, with 5-nitrouracil showing a substantial synergistic effect on the resistant strain. Eighteen UGT genes were identified in F. occidentalis, with gene expansion and duplication observed within families UGT466, 467, and 468. Ten out of the eighteen UGTs exhibited higher expression levels in NIL-R, specifically FoUGT466B1, FoUGT468A3, and FoUGT468A4 consistently being upregulated across nymphs, males, and females. RNAi-based functional validation targeting these three UGT genes led to increased susceptibility to spinosad in a life stage-, sex-, and dose-dependent manner. These results indicate that UGTs are indeed involved in spinosad resistance in F. occidentalis, and the effects are dependent on life stage, sex, and dose. Therefore, sustainable control for F. occidentalis resistance should always consider these differential responses.

Uridine Diphosphate-Glycosyltransferase RpUGT344D38 Contributes to λ-Cyhalothrin Resistance in Rhopalosiphum padi

Uridine diphosphate-glycosyltransferase (UGT) is a key phase II enzyme in the insect detoxification system. Pyrethroids are commonly used to control the destructive wheat aphid Rhopalosiphum padi. In this study, we found a highly expressed UGT gene, RpUGT344D38, in both λ-cyhalothrin (LCR)- and bifenthrin (BTR)-resistant strains of R. padi. After exposure to λ-cyhalothrin and bifenthrin, the expression levels of RpUGT344D38 were significantly increased in the resistant strains. Knockdown of RpUGT344D38 did not affect the resistance of BTR, but it did significantly increase the susceptibility of LCR aphids to λ-cyhalothrin. Molecular docking analysis demonstrated that RpUGT344D38 had a stable binding interaction with both bifenthrin and λ-cyhalothrin. The recombinant RpUGT344D38 was able to metabolize 50% of λ-cyhalothrin. This study provides a comprehensive analysis of the role of RpUGT344D38 in the resistance of R. padi to bifenthrin and λ-cyhalothrin, contributing to a better understanding of aphid resistance to pyrethroids.

Characterization and functional analysis of UDP-glycosyltransferases reveal their contribution to phytochemical flavone tolerance in Spodoptera litura

Efficient detoxification is the key factor for phytophagous insect to adapt to phytochemicals. However, the role of uridine diphosphate (UDP)-glycosyltransferases (UGTs) in insect anti-defense to phytochemical flavone is largely unknown. In this study, 52 UGT genes were identified in Spodoptera litura and they presented evident gene duplication. UGT played a crucial part in larval tolerance to flavone because the enzyme activity and transcriptional level of 77 % UGT members were remarkably upregulated by flavone administration and suppression of UGT enzyme activity and gene expressions significantly increased larval susceptibility to flavone. Bacteria coexpressing UGTs had high survival rates under flavone treatment and flavone was dramatically metabolized by UGT recombinant cells, which indicated the involvement of UGTs in flavone detoxification. What's more, ecdysone pathway was activated by flavone. Topical application of 20-hydroxyecdysone highly upregulated UGT enzyme activity and more than half of UGT expressions. The effects were opposite when ecdysone receptor (EcR) and ultraspiracle (USP)-mediated ecdysone signaling pathway was inhibited. Furtherly, promoter reporter assays of 5 UGT genes showed that their transcription activities were notably increased by cotransfection with EcR and USP. In consequence, this study suggested that UGTs were involved in flavone detoxification and their transcriptional expressions were regulated by ecdysone pathway.

Overexpression and nonsynonymous mutations of UDP-glycosyltransferases are potentially associated with pyrethroid resistance in Anopheles funestus

UDP-glycosyltransferases (UGTs) enzymes are pivotal in insecticide resistance by transforming hydrophobic substrates into more hydrophilic forms for efficient cell elimination. This study provides the first comprehensive investigation of Anopheles funestus UGT genes, their evolution, and their association with pyrethroid resistance. We employed a genome-wide association study using pooled sequencing (GWAS-PoolSeq) and transcriptomics on pyrethroid-resistant An. funestus, along with deep-targeted sequencing of UGTs in 80 mosquitoes Africa-wide. UGT310B2 was consistently overexpressed Africa-wide and significant gene-wise Fst differentiation was observed between resistant and susceptible populations: UGT301C2 and UGT302A3 in Malawi, and UGT306C2 in Uganda. Additionally, nonsynonymous mutations in UGT genes were identified. Gene-wise Tajima's D density curves provide insights into population structures within populations across these countries, supporting previous observations. These findings have important implications for current An. funestus control strategies facilitating the prediction of cross-resistance to other UGT-metabolised polar insecticides, thereby guiding more effective and targeted insecticide resistance management efforts.

Per- and polyfluoroalkyl substance (PFAS) mixtures induce gut microbiota dysbiosis and metabolic disruption in silkworm (Bombyx mori L.)

Mixed legacy and emerging per- and polyfluoroalkyl substances (PFASs) are commonly found in soil and dust; however, the potential toxicity of PFAS mixtures (mPFASs) in insects is unknown. Using 16S rRNA gene sequencing and transcriptome sequencing (RNA-Seq), we evaluated the adverse effects of mPFASs on silkworms, a typical lepidopteran insect. After exposure to mPFASs, the silkworm midgut was enriched with high levels of PFASs, which induced histopathological changes. The composition of the midgut microbiota was significantly affected by mPFAS exposure, and functional predictions revealed significant disruption of some metabolic pathways. RNA-seq analysis revealed that mPFASs significantly changed the transcription profiles. Functional enrichment analysis of the differentially expressed genes also revealed that biological processes related to metabolic pathways and the digestive system were significantly affected, similar to the results of the gut microbiota analysis, suggesting that mPFAS exposure had an adverse effect on the metabolic function of silkworms and may further affect their normal growth. Finally, the significant correlation between abundance changes in the gut microbiota and metabolism/digestion-related genes further highlighted the role of the gut microbiota in mPFAS-related processes affecting the metabolic functions of silkworms. To our knowledge, this study is the first to evaluate the toxic effects of mPFASs in insects and provide basic data for further PFAS toxicity investigations in insects and comprehensive ecological risk assessments of mPFASs.

Phylogenomics of the Ecdysteroid Kinase-like (EcKL) Gene Family in Insects Highlights Roles in Both Steroid Hormone Metabolism and Detoxification

The evolutionary dynamics of large gene families can offer important insights into the functions of their individual members. While the ecdysteroid kinase-like (EcKL) gene family has previously been linked to the metabolism of both steroid molting hormones and xenobiotic toxins, the functions of nearly all EcKL genes are unknown, and there is little information on their evolution across all insects. Here, we perform comprehensive phylogenetic analyses on a manually annotated set of EcKL genes from 140 insect genomes, revealing the gene family is comprised of at least 13 subfamilies that differ in retention and stability. Our results show the only two genes known to encode ecdysteroid kinases belong to different subfamilies and therefore ecdysteroid metabolism functions must be spread throughout the EcKL family. We provide comparative phylogenomic evidence that EcKLs are involved in detoxification across insects, with positive associations between family size and dietary chemical complexity, and we also find similar evidence for the cytochrome P450 and glutathione S-transferase gene families. Unexpectedly, we find that the size of the clade containing a known ecdysteroid kinase is positively associated with host plant taxonomic diversity in Lepidoptera, possibly suggesting multiple functional shifts between hormone and xenobiotic metabolism. Our evolutionary analyses provide hypotheses of function and a robust framework for future experimental studies of the EcKL gene family. They also open promising new avenues for exploring the genomic basis of dietary adaptation in insects, including the classically studied coevolution of butterflies with their host plants.

RpUGT344J7 is involved in the reproduction switch of Rhopalosiphum padi with holocyclic life cycle

Many aphid species exhibit both cyclical parthenogenesis (CP) and the obligate parthenogenesis (OP) life history, which are genetically determined. In CP aphid lineages, the parthenogenetic individuals can switch from asexual to sexual reproduction quickly in response to environmental factors such as changes in photoperiod and temperature. However, the OP aphid lineages do not undergo sexual reproduction under any conditions. So far, mechanisms underlying the reproduction switch in CP aphids have not been fully elucidated. Rhopalosiphum padi, a serious worldwide insect pest of wheat, has both CP and OP lineages. Uridine diphosphate-glycosyltransferases (UGTs) are enzymes that participate in the metabolic detoxification of xenobiotics. Here, we identified 43 RpUGT genes from R. padi genome and transcriptome sequences, and found that: (1) the UGT content of the CP lineage was significantly higher than that in the OP lineage at the key time points when CP lineage mainly produce virginoparae, gynoparae, and males under inducing condition, while there were no significant difference under normal conditions; (2) RpUGT344J7 gene was highly expressed during the time points when CP lineages produce gynopara and males; (3) the critical time points for CP lineages to produce virginoparaee, gynoparae, and males were affected when the CP lineages were injected with dsRpUGT344J7; (4) the knockdown of RpUGT344J7 caused a significant reduction in the total number of virginoparae, gynoparae, and males in the offspring under inducing condition. The findings contribute to our understanding of the molecular mechanisms underlying the quick shift from asexual to sexual reproduction in aphid species.

Whitefly-Plant Interactions: An Integrated Molecular Perspective

The rapid advances in available transcriptomic and genomic data and our understanding of the physiology and biochemistry of whitefly-plant interactions have allowed us to gain new and significant insights into the biology of whiteflies and their successful adaptation to host plants. In this review, we provide a comprehensive overview of the mechanisms that whiteflies have evolved to overcome the challenges of feeding on phloem sap. We also highlight the evolution and functions of gene families involved in host perception, evaluation, and manipulation; primary metabolism; and metabolite detoxification. We discuss the emerging themes in plant immunity to whiteflies, focusing on whitefly effectors and their sites of action in plant defense-signaling pathways. We conclude with a discussion of advances in the genetic manipulation of whiteflies and the potential that they hold for exploring the interactions between whiteflies and their host plants, as well as the development of novel strategies for the genetic control of whiteflies.

Over-expression of UDP-glycosyltransferase UGT353G2 confers resistance to neonicotinoids in whitefly (Bemisia tabaci)

The whitefly, Bemisia tabaci, comes up high metabolic resistance to most neonicotinoids in long-term evolution, which is the key problem of pest control. UGT glycosyltransferase, as a secondary detoxification enzyme, plays an indispensable role in detoxification metabolism. In this study, UGT inhibitors, 5-nitrouracil and sulfinpyrazone, dramatically augmented the toxic damage of neonicotinoids to B. tabaci. A UGT named UGT353G2 was identified in whitefly, which was notably up-regulated in resistant strain (3.92 folds), and could be induced by most neonicotinoids. Additionally, the using of RNA interference (RNAi) suppresses UGT353G2 substantially increased sensitivity to neonicotinoids in resistant strain. Our results support that UGT353G2 may be involved in the neonicotinoids resistance of whitefly. These findings will help further verify the functional role of UGTs in neonicotinoid resistance.

Transposon accumulation at xenobiotic gene family loci in aphids

The evolution of resistance is a major challenge for the sustainable control of pests and pathogens. Thus, a deeper understanding of the evolutionary and genomic mechanisms underpinning resistance evolution is required to safeguard health and food production. Several studies have implicated transposable elements (TEs) in xenobiotic-resistance evolution in insects. However, analyses are generally restricted to one insect species and/or one or a few xenobiotic gene families (XGFs). We examine evidence for TE accumulation at XGFs by performing a comparative genomic analysis across 20 aphid genomes, considering major subsets of XGFs involved in metabolic resistance to insecticides: cytochrome P450s, glutathione S-transferases, esterases, UDP-glucuronosyltransferases, and ABC transporters. We find that TEs are significantly enriched at XGFs compared with other genes. XGFs show similar levels of TE enrichment to those of housekeeping genes. But unlike housekeeping genes, XGFs are not constitutively expressed in germline cells, supporting the selective enrichment of TEs at XGFs rather than enrichment owing to chromatin availability. Hotspots of extreme TE enrichment occur around certain XGFs. We find, in aphids of agricultural importance, particular enrichment of TEs around cytochrome P450 genes with known functions in the detoxification of synthetic insecticides. Our results provide evidence supporting a general role for TEs as a source of genomic variation at host XGFs and highlight the existence of considerable variability in TE content across XGFs and host species. These findings show the need for detailed functional verification analyses to clarify the significance of individual TE insertions and elucidate underlying mechanisms at TE-XGF hotspots.

UDP-glucuronosyltransferase is involved in susceptibility of Chironomus kiiensis Tokunaga, 1936 (Diptera: Chironomidae) to insecticides

UDP-glucuronosyltransferases (UGTs) could transform various exogenous and endogenous compounds, which help detoxification of pesticides in insects. To investigate the role of UGTs in the detoxification metabolism of insecticides in Chironomus kiiensis, CkUGT302M1, CkUGT302N1, CkUGT308N1 and CkUGT36J1 genes were identified with 1449-1599 bp encoding 482-532 amino acids. Four UGT genes shared 40.86∼53.36% identity with other homologous insect species, and expressed in all developmental stages, notably in the larval and adult stages. Expression of CkUGTs was higher in the gastric caecum, midgut and head. Moreover, CkUGTs expression and activity were significantly increased in C. kiiensis larvae in exposure to sublethal concentrations of carbaryl, deltamethrin and phoxim, respectively. To further explore the functions of UGT genes, the CkUGT308N1 was effectively silenced in 4th instar C. kiiensis larvae by RNA interference, which resulted in the mortality of dsCkUGT308N1 treated larvae increased by 71.43%, 111.11% and 62.50% under sublethal doses of carbaryl, deltamethrin and phoxim at the 24-h time point, respectively. The study revealed that the CkUGT308N1 gene in C. kiiensis could contribute to the metabolism of pesticides and provide a scientific basis for evaluating the water pollution of pesticides.

Chlorpyrifos and chlorfenapyr resistance in Spodoptera frugiperda (Lepidoptera: Noctuidae) relies on UDP-glucuronosyltransferases

Fall armyworm, Spodoptera frugiperda (J. E. Smith), has become an important agricultural pest worldwide. S. frugiperda is mainly controlled by the chemical insecticides, whereas the frequent application of insecticides would result in the resistance development. Insect uridine diphosphate-glucuronosyltransferases (UGTs), as phase II metabolism enzymes, play vital roles in the breakdown of endobiotic and xenobiotics. In this study, 42 UGT genes were identified by RNA-seq, including 29 UGT genes were elevated compared to the susceptible population, and the transcript levels of 3 UGTs (UGT40F20, UGT40R18, and UGT40D17) were increased by more than 2.0-fold in the field populations. Expression pattern analysis revealed that S. frugiperda UGT40F20, UGT40R18, and UGT40D17 were increased by 6.34-, 4.26-, and 8.28-fold, compared the susceptible populations, respectively. The expression of UGT40D17, UGT40F20, and UGT40R18 was affected after exposure to phenobarbital, chlorpyrifos, chlorfenapyr, sulfinpyrazone, and 5-nitrouracil. The induced expression of UGT genes may have improved UGT enzymatic activity, while the inhibition of UGTs genes expression may decreased UGT enzymatic activity. Sulfinpyrazone, and 5-nitrouracil, significantly increased the toxicity of chlorpyrifos and chlorfenapyr, as well as phenobarbital significantly reduced the toxicity of chlorpyrifos and chlorfenapyr against the susceptible populations and field populations of S. frugiperda. The suppression of UGTs (UGT40D17, UGT40F20, and UGT40R18) significantly increased the insensitivity of the field populations to chlorpyrifos and chlorfenapyr. These findings strongly supported our viewpoint that UGTs may play a critical role in insecticide detoxification. This study provides a scientific basis for the management of S. frugiperda.

Exploring the role of detoxification genes in the resistance of Bursaphelenchus xylophilus to different exogenous nematicidal substances using transcriptomic analyses

Bursaphelenchus xylophilus (Pine wood nematode, PWN) has become a worldwide forest disease due to its rapid infection ability, high lethality and difficulty in control. The main means of countering B. xylophilus is currently chemical control, but nematicides can present problems such as environmental pollution and drug resistance. The development of novel environmentally-friendly nematicides has thus become a focus of recent research. In this study, BxUGT3 and BxUGT34, which might be related to detoxification, were investigated by comparing transcriptomic and WGCNA approaches. Three other genes with a similar expression pattern, BxUGT13, BxUGT14, and BxUGT16, were found by gene family analysis. Further bioassays and qPCR assays confirmed that these five genes showed significant changes in transcript levels upon exposure to α-pinene and carvone, demonstrating that they respond to exogenous nematicidal substances. Finally, RNAi and bioassays showed that B. xylophilus with silenced BxUGT16 had increased mortality in the face of α-pinene and carvone stress, suggesting that BxUGT16 plays an important role in detoxification. Taken together, this study used novel molecular research methods, explored the detoxification mechanism of B. xylophilus at a transcriptomic level, and revealed a molecular target for the development of novel biopesticides.

UGT2B13 and UGT2C1 are involved in lambda-cyhalothrin resistance in Rhopalosiphum padi

Uridine diphosphate-glucuronosyltransferases (UGTs) are major multifunctional detoxification phase II enzymes involved in the metabolic detoxification of xenobiotics. However, their roles in insecticides resistance are still unclear. In this study, we identified two UGTs genes (UGT2B13 and UGT2C1) in Rhopalosiphum padi, a serious insect pest of wheat worldwide. Bioassays results showed that the resistance ratio of R. padi resistance strain (LC-R) to lambda-cyhalothrin (LC) was 2963.8 fold. The roles of UGT2B13 and UGT2C1 in lambda-cyhalothrin resistance were evaluated. Results indicated that the UGTs contents were significantly increased in the LC resistant strain of R. padi. UGT2B13 and UGT2C1 were significantly overexpressed in the LC-R strain. Transcription levels of UGT2B13 and UGT2C1 were relatively higher in the gut of LC-R strain. RNA interference (RNAi) of UGT2B13 or UGT2C1 significantly decreased the UGTs contents of the LC-R aphids and increased mortality of R. padi exposure to the LC50 concentration of LC. This study provides a new view that UGTs are involved in LC resistance of R. padi. The findings will promote further work to detailed the functions of UGTs in the metabolism resistance of insects to insecticides.

Overexpression of the UDP-glycosyltransferase UGT34A23 confers resistance to the diamide insecticide chlorantraniliprole in the tomato leafminer, Tuta absoluta

The tomato leafminer, Tuta absoluta, is an invasive crop pest that has evolved resistance to many of the insecticides used for its control. To facilitate the investigation of the underpinning mechanisms of resistance in this species we generated a contiguous genome assembly using long-read sequencing data. We leveraged this genomic resource to investigate the genetic basis of resistance to the diamide insecticide chlorantraniliprole in Spanish strains of T. absoluta that exhibit high levels of resistance to this insecticide. Transcriptomic analyses revealed that, in these strains, resistance is not associated with previously reported target-site mutations in the diamide target-site, the ryanodine receptor, but rather is associated with the marked overexpression (20- to >100-fold) of a gene encoding a UDP-glycosyltransferase (UGT). Functional expression of this UGT, UGT34A23, via ectopic expression in Drosophila melanogaster demonstrated that it confers strong and significant resistance in vivo. The genomic resources generated in this study provide a powerful resource for further research on T. absoluta. Our findings on the mechanisms underpinning resistance to chlorantraniliprole will inform the development of sustainable management strategies for this important pest.

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.

Adaptive Changes in Detoxification Metabolism and Transmembrane Transport of Bombyx mori Malpighian Tubules to Artificial Diet

The high adaptability of insects to food sources has contributed to their ranking among the most abundant and diverse species on Earth. However, the molecular mechanisms underlying the rapid adaptation of insects to different foods remain unclear. We explored the changes in gene expression and metabolic composition of the Malpighian tubules as an important metabolic excretion and detoxification organ in silkworms (Bombyx mori) fed mulberry leaf and artificial diets. A total of 2436 differentially expressed genes (DEGs) and 245 differential metabolites were identified between groups, with the majority of DEGs associated with metabolic detoxification, transmembrane transport, and mitochondrial function. Detoxification enzymes, such as cytochrome P450 (CYP), glutathione-S-transferase (GST), and UDP-glycosyltransferase, and ABC and SLC transporters of endogenous and exogenous solutes were more abundant in the artificial diet group. Enzyme activity assays confirmed increased CYP and GST activity in the Malpighian tubules of the artificial diet-fed group. Metabolome analysis showed increased contents of secondary metabolites, terpenoids, flavonoids, alkaloids, organic acids, lipids, and food additives in the artificial diet group. Our findings highlight the important role of the Malpighian tubules in adaptation to different foods and provide guidance for further optimization of artificial diets to improve silkworm breeding.

Expression profile of the entire detoxification gene inventory of the western honeybee, Apis mellifera across life stages

The western honeybee, Apis mellifera, is a managed pollinator of many crops and potentially exposed to a wide range of foreign compounds, including pesticides throughout its life cycle. Honeybees as well as other insects recruit molecular defense mechanisms to facilitate the detoxification of xenobiotic compounds. The inventory of detoxification genes (DETOXome) is comprised of five protein superfamilies: cytochrome P450 monooxygenases (P450), carboxylesterases, glutathione S-transferases (GST), UDP-glycosyl transferases (UGT) and ATP-binding cassette (ABC) transporters. Here we characterized the gene expression profile of the entire honeybee DETOXome by analyzing 47 transcriptomes across the honeybee life cycle, including different larval instars, pupae, and adults. All life stages were well separated by principal component analysis, and K-means clustering revealed distinct temporal patterns of gene expression. Indeed, >50% of the honeybee detoxification gene inventory is found in one cluster and follows strikingly similar expression profiles, i.e., increased expression during larval development, followed by a sharp decline after pupation and a steep increase again in adults. This cluster includes 29 P450 genes dominated by CYP3 and CYP4 clan members, 15 ABC transporter genes mostly belonging to the ABCC subfamily and 13 carboxylesterase genes including almost all members involved in dietary/detox and hormone/semiochemical processing. RT-qPCR analysis of selected detoxification genes from all families revealed high expression levels in various tissues, especially Malpighian tubules, fatbody and midgut, supporting the view that these tissues are essential for metabolic clearance of environmental toxins and pollutants in honeybees. Our study is meant to spark further research on the molecular basis of detoxification in this critical pollinator to better understand and evaluate negative impacts from potentially toxic substances. Additionally, the entire gene set of 47 transcriptomes collected and analyzed provides a valuable resource for future honeybee research across different disciplines.

Intraspecific diversity in the mechanisms underlying abamectin resistance in a cosmopolitan pest

Pesticide resistance relies on a myriad of mechanisms, ranging from single mutations to a complex and polygenic architecture, and it involves mechanisms such as target-site insensitivity, metabolic detoxification, or a combination of these, with either additive or synergistic effects. Several resistance mechanisms against abamectin, a macrocyclic lactone widely used in crop protection, have been reported in the cosmopolitan pest Tetranychus urticae. However, it has been shown that a single mechanism cannot account for the high levels of abamectin resistance found across different mite populations. Here, we used experimental evolution combined with bulked segregant analyses to map quantitative trait loci (QTL) associated with abamectin resistance in two genetically unrelated populations of T. urticae. In these two independent QTL mapping experiments, three and four QTLs were identified, of which three were shared between experiments. Shared QTLs contained genes encoding subunits of the glutamate-gated chloride channel (GluCl) and harboured previously reported mutations, including G314D in GluCl1 and G326E in GluCl3, but also novel resistance candidate loci, including DNA helicases and chemosensory receptors. Surprisingly, the fourth QTL, present only in only one of the experiments and thus unique for one resistant parental line, revealed a non-functional variant of GluCl2, suggesting gene knock-out as resistance mechanism. Our study uncovers the complex basis of abamectin resistance, and it highlights the intraspecific diversity of genetic mechanisms underlying resistance in a cosmopolitan pest.

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.

Effects of azadirachtin on detoxification-related gene expression in the fat bodies of the fall armyworm, Spodoptera frugiperda

The fall armyworm, Spodoptera frugiperda, has become a worldwide pest and threatens world food production. A previous study indicated that azadirachtin, the most effective botanical insecticide for S. frugiperda, inhibits larval growth of the insect. The effect of azadirachtin on the tissues of the larvae, however, remains to be determined. In this study, the effects of azadirachtin on the structure of fat bodies were analyzed. Comparative transcriptomic analysis was conducted between controls and samples treated with 0.1 μg/g azadirachtin for 7 days to explore potential relevant mechanisms. The expression of 5356 genes was significantly affected after azadirachtin treatment, with 3020 up-regulated and 2336 down-regulated. Among them, 137 encode detoxification enzymes, including 53 P450s, 20 GSTs, 27 CarEs, 16 UGTs, and 12 ABC transporters. Our results indicated that azadirachtin could destroy fat body structure and change the mRNA levels of detoxification-related genes. The up-regulated genes encoding detoxification enzymes might be related to detoxifying azadirachtin. Our results elucidate a preliminary mechanism of azadirachtin detoxification in the fat bodies of S. frugiperda larvae.

Contribution of UDP-glycosyltransferases to chlorpyrifos resistance in Nilaparvata lugens

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

Identification and characterization of UDP-glycosyltransferase genes and the potential role in response to insecticides exposure in Bactrocera dorsalis

BACKGROUND

The oriental fruit fly, Bactrocera dorsalis (Hendel) is a worldwide pest damaging a wide range of hosts. Due to the long-term indiscriminate use of insecticides, B. dorsalis has developed serious resistance to several insecticides. UDP-glycosyltransferases (UGTs) are secondary metabolic enzymes involved in biotransformation and play an important role in the metabolism of plant secondary metabolites and synthetic insecticides in insects. Thus, we suspect that UGTs in B. dorsalis play an important role in insecticide tolerance.

RESULTS

In this study, 31 UGT genes were identified in the genome of B. dorsalis, belonging to 13 subfamilies. Real-time quantitative polymerase chain reaction (RT-qPCR) results revealed that 12 UGT genes were highly expressed in the antennae, midgut, Malpighian tubule and fat body. The mRNA expressions of 17 UGT genes were up-regulated upon exposure to λ-cyhalothrin, imidacloprid, abamectin and chlorpyrifos. Knockdown of the selected five UGT genes (BdUGT301D2, BdUGT35F2, BdUGT36K2, BdUGT49D2, BdUGT50B5) by RNA interference increased the mortality of B. dorsalis from 9.29% to 27.22% upon exposure to four insecticides.

CONCLUSION

The abundance of UGTs in B. dorsalis is similar to other insect species, and 12 out of 31 UGTs were specifically expressed in metabolic tissues, suggesting a key role in detoxification. Down-regulation of five selected UGT genes increased the susceptibility of B. dorsalis to various insecticides, indicating that UGTs may play an important role in tolerance of B. dorsalis to multiple insecticides. © 2022 Society of Chemical Industry.

Silkworm arylsulfatase in the midgut content is expressed in the silk gland and fed via smearing on the food from the spinneret

We found the activity of arylsulfatase in the midgut contents of the silkworm, Bombyx mori. We identified a 60-kDa protein that comigrates with the activity on a column chromatography following ammonium sulfate precipitation. Based on its partial amino acid sequence, we searched for its coding gene using Basic Local Alignment Search Tool (BLAST) and identified KWMTBOMO05106. Transcriptional data suggest a specific expression of the gene in middle silk glands. The majority (80%) of arylsulfatase activity was found in the silk glands, concurring the specific transcription in the silk gland. Observing the feeding behaviour of the silkworm, we found that silkworms smear a mucus secretes from the spinneret on the food pellet as they feed on. Arylsulfatase activity was also detected in the food pellet bitten by the silkworm as well as in the gut content. Furthermore, arylsulfatase activity was not detected either in the food pellet and in the gut content when silkworms had obstructed the spinneret. These results suggest that arylsulfatase is secreted from the silk glands and may contribute to digestive function.

Improved genome assembly provides new insights into the environmental adaptation of the American cockroach, Periplaneta americana

The synanthropic pest and a model organism for entomological research, American cockroach, Periplaneta americana (Linnaeus), can survive in unfavorable environments for humans. To investigate the genetic mechanisms of success in environmental adaptation of P. americana, we de novo reassembled its whole genome based on next-generation sequencing and PacBio sequencing. The final genome reassembly consisted of approximately 3.34 Gb with scaffold N50 of 465.51 Kb. The completeness (95.4%) of the complete genome was evaluated with single-copy orthologous genes using BUSCO. We identified 18,618 protein-coding genes, 16,443 (88.32%) of which were well supported by public protein databases. We identified 482.04 Mb (approximately 14.45%) repeat elements, 1,385,093 perfect microsatellites simple sequence repeats in P. americana genome, which was higher than other four Blattaria insects. Comparative genomics analysis revealed obvious expansion in the gene families associated with chemoreception (olfactory receptors, gustatory receptors, ionotropic glutamate receptors, chemosensory protein, and sensory neuron membrane protein), which provided the necessary information for functional characterization of the chemosensory receptors of P. americana, with potential for new or refined applications of semiochemicals-based control of this pest insect. Similarly, gene families (cytochrome P450s, carboxyl/choline esterases, and UDP-glycosyl-transferases) encoding receptors for bitter or toxic substances and detoxification enzymes were obviously expanded in P. americana, enabling its ability to detect and detoxify many toxins. Enrichment analysis of positively selected genes in P. americana revealed items associated with metabolic process and catalytic activity, which possibly contributed to the pesticide resistance of P. americana. We also analyzed the homologs to antimicrobial peptide genes reported in the Drosophila genome, and identified two attacins and seven defensins in P. americana. Our data and findings will substantially facilitate molecular studies in P. americana, including elucidation of detoxification mechanisms of xenobiotic, as well as development of new pest management strategies for the control of pests like P. americana.

Transcriptomic Identification and Expression Profile Analysis of Odorant-Degrading Enzymes from the Asian Corn Borer Moth, Ostrinia furnacalis

The Asian corn borer moth Ostrinia furnacalis is an important lepidopteran pest of maize in Asia. Odorant-degrading enzymes (ODEs), including carboxylesterases (CCEs), glutathione S-transferases (GSTs), cytochrome P450s (CYPs), UDP-glycosyltransferases (UGTs), and aldehyde oxidases (AOXs), are responsible for rapid inactivation of odorant signals in the insect antennae. In this study, we performed a transcriptome assembly for the antennae of O. furnacalis to identify putative ODE genes. Transcriptome sequencing revealed 35,056 unigenes, and 21,012 (59.94%) of these were annotated by searching against the reference sequences in the NCBI non-redundant (NR) protein database. For functional classification, these unigenes were subjected to Gene Ontology (GO), Eukaryotic Orthologous Groups (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations. We identified 79 genes encoding putative ODEs: 19 CCEs, 17 GSTs, 24 CYPs, 13 UGTs, and 6 AOXs. BLASTX best hit results indicated that these genes shared quite high amino acid identities with their respective orthologs from other lepidopteran species. Reverse transcription-quantitative PCR showed that OfurCCE2, OfurCCE5, and OfurCCE18 were enriched in male antennae, while OfurCCE7 and OfurCCE10 were enriched in female antennae. OfurCCE14 and OfurCCE15 were expressed at near-equal amounts in the antennae of both sexes. Our findings establish a solid foundation for future studies aimed at understanding the olfactory functions of these genes in O. furnacalis.

Transcriptome Analysis to Identify Responsive Genes under Sublethal Concentration of Bifenazate in the Diamondback Moth, Plutella xylostella (Linnaeus, 1758) (Lepidoptera: Plutellidae)

Bifenazate is a novel acaricide that has been widely used to control spider mites. Interestingly, we found bifenazate had a biological activity against the diamondback moth (Plutella xylostella), one of the most economically important pests on crucifer crops around the world. However, the molecular mechanisms underlying the response of P. xylostella to bifenazate treatment are not clear. In this study, we first estimated the LC₃₀ dose of bifenazate for third-instar P. xylostella larvae. Then, in order to identify genes that respond to the treatment of this insecticide, the comparative transcriptome profiles were used to analyze the gene expression changes in P. xylostella larvae after exposure to LC₃₀ of bifenazate. In total, 757 differentially expressed genes (DEGs) between bifenazate-treated and control P. xylostella larvae were identified, in which 526 and 231 genes were up-regulated and down-regulated, respectively. The further Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that the xenobiotics metabolisms pathway was significantly enriched, with ten detoxifying enzyme genes (four P450s, five glutathione S-transferases (GSTs), and one UDP-Glucuronosyltransferase (UGT)) were up-regulated, and their expression patterns were validated by qRT-PCR as well. Interestingly, the present results showed that 17 cuticular protein (CP) genes were also remarkably up-regulated, including 15 CPR family genes. Additionally, the oxidative phosphorylation pathway was found to be activated with eight mitochondrial genes up-regulated in bifenazate-treated larvae. In contrast, we found some genes that were involved in tyrosine metabolism and purine pathways were down-regulated, indicating these two pathways of bifenazate-exposed larvae were significantly inhibited. In conclusion, the present study would help us to better understand the molecular mechanisms of sublethal doses of bifenazate detoxification and action in P. xylostella.

Genome analysis of Phrixothrix hirtus (Phengodidae) railroad worm shows the expansion of odorant-binding gene families and positive selection on morphogenesis and sex determination genes

Among bioluminescent beetles of the Elateroidea superfamily, Phengodidae is the third largest family, with 244 bioluminescent species distributed only in the Americas, but is still the least studied from the phylogenetic and evolutionary points of view. The railroad worm Phrixothrix hirtus is an essential biological model and symbolic species due to its bicolor bioluminescence, being the only organism that produces true red light among bioluminescent terrestrial species. Here, we performed partial genome assembly of P. hirtus, combining short and long reads generated with Illumina sequencing, providing the first source of genomic information and a framework for comparative analyses of the bioluminescent system in Elateroidea. This is the largest genome described in the Elateroidea superfamily, with an estimated size of ∼3.4 Gb, displaying 32 % GC content, and 67 % transposable elements. Comparative genomic analyses showed a positive selection of genes and gene family expansion events of growths and morphogenesis gene products, which could be associated with the atypical anatomical development and morphogenesis found in paedomorphic females and underdeveloped males. We also observed gene family expansion among distinct odorant-binding receptors, which could be associated with the pheromone communication system typical of these beetles, and retrotransposable elements. Common genes putatively regulating bioluminescence production and control, including two luciferase genes corresponding to lateral lanterns green-emitting and head lanterns red-emitting luciferases with 7 exons and 6 introns, and genes potentially involved in luciferin biosynthesis were found, indicating that there are no clear differences about the presence or absence of gene families associated with bioluminescence in Elateroidea.

Detoxication and bioconversion of aflatoxin B1 by yellow mealworms (Tenebrio molitor): A sustainable approach for valuable larval protein production from contaminated grain

Yellow mealworm (Tenebrio molitor) is a supplementary protein source for food and feed and represents a promising solution to manage grain contaminated with Aflatoxin B₁ (AFB₁). In this study, AFB₁ present in different concentrations in wheat bran was treated and removed via bioconversion by yellow mealworm of different instars, with emphasis on the bioconversion performance and metabolism of AFB₁. Upon application of wheat bran spiked with 100 μg/kg AFB₁ to 5th-6th instar yellow mealworms, the conversion rate of AFB₁ was up to 87.85 %. Low level of AFB₁ (< 2 μg/kg) was accumulated in the larval bodies, and the survival rate, development and nutrition contents of yellow mealworm were not significantly affected. It was revealed that 1 kg of wheat bran contaminated with AFB₁ increased the weight of yellow mealworms from 138 g to 469 g, containing approximately 103 g of protein. The bioconversion of AFB₁ by yellow mealworms led to generation of 13 metabolites in the frass and 3 metabolites in the larvae. AFB₁ was detoxicated and removed via phase I metabolism comprising reduction, dehydrogenation, hydration, demethylation, hydroxylation, decarbonylation and ketoreduction, followed by phase II metabolism involving conjugation of amino acid, glucoside and glutathione (GSH). The toxicity of AFB₁ metabolites was deemed lower than that of AFB₁ according to their structures. This study provides a sustainable approach and theoretical foundation on using yellow mealworms for cleaner grain contamination management and valuable larval protein production via bioconversion of food and feed contaminated by AFB₁.

RNAi-mediated knockdown of gut receptor-like genes prohibitin and α-amylase altered the susceptibility of Galleria mellonella to Cry1AcF toxin

Background

Due to the prolonged usage of Bt-based biopesticides and Bt-transgenic crops worldwide, insects are continually developing resistance against Cry toxins. This resistance may occur if any mechanistic step in the insecticidal process is disrupted possibly because of the alteration in Cry-receptor binding affinity due to mutation in receptor genes. Compared to other lepidopteran insects, Cry receptor-related research has made asymmetric progress in the model insect Galleria mellonella.

Results

Present study describes the molecular characterization and functional analysis of five Cry toxin receptor-related genes (prohibitin, GLTP, α-amylase, ADAM and UDP-GT) and a gut repair gene (arylphorin) from the gut tissues of G. mellonella. Protein–protein docking analysis revealed that Cry1AcF putatively binds with all the five candidate proteins, suggesting their receptor-like function. These receptor-like genes were significantly overexpressed in the gut tissues of fourth-instar G. mellonella larvae upon early exposure to a sub-lethal dose of Cry1AcF toxin. However, targeted knockdown (by using bacterially-expressed dsRNAs) of these genes led to variable effect on insect susceptibility to Cry1AcF toxin. Insects pre-treated with prohibitin and α-amylase dsRNA exhibited significant reduction in Cry1AcF-induced mortality, suggesting their probable role as Cry receptor. By contrast, insects pre-treated with GLTP, ADAM and UDP-GT dsRNA exhibited no significant decline in mortality. This maybe explained by the possibility of RNAi feedback regulation (as few of the receptors belong to multigene family) or redundant role of GLTP, ADAM and UDP-GT in Cry intoxication process.

Conclusion

Since the laboratory culture of G. mellonella develop Bt resistance quite rapidly, findings of the current investigation may provide some useful information for future Cry receptor-related research in the model insect.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08843-8.

Alternative transcript splicing regulates UDP-glucosyltransferase-catalyzed detoxification of DIMBOA in the fall armyworm (Spodoptera frugiperda)

Herbivorous insects often possess the ability to detoxify chemical defenses from their host plants. The fall armyworm (Spodoptera frugiperda), which feeds principally on maize, detoxifies the maize benzoxazinoid 2,4-dihydroxy-7-methoxy-1,4-benzoxazin-3-one (DIMBOA) by stereoselective re-glucosylation using a UDP-glucosyltransferase, SfUGT33F28. SfUGT33F28 activity is induced by feeding on a DIMBOA-containing diet, but how this induction is regulated is unknown. In the present work, we describe the alternative splicing of the SfUGT33F28 transcript. Variant transcripts are differentially expressed in response to DIMBOA, and this transcriptional response is mediated by an insect aryl hydrocarbon receptor. These variants have large deletions leading to the production of truncated proteins that have no intrinsic UGT activity with DIMBOA but interact with the full-length enzyme to raise or lower its activity. Therefore, the formation of SfUGT33F28 splice variants induces DIMBOA-conjugating UGT activity when DIMBOA is present in the insect diet and represses activity in the absence of this plant defense compound.

Host and gut microbiome modulate the antiparasitic activity of nectar metabolites in a bumblebee pollinator

Antimicrobial nectar secondary metabolites can support pollinator health by preventing or reducing parasite infections. To better understand the outcome of nectar metabolite-parasite interactions in pollinators, we determined whether the antiparasitic activity was altered through chemical modification by the host or resident microbiome during gut passage. We investigated this interaction with linden (Tilia spp.) and strawberry tree (Arbutus unedo) nectar compounds. Unedone from A. unedo nectar inhibited the common bumblebee gut parasite Crithidia bombi in vitro and in Bombus terrestris gynes. A compound in Tilia nectar, 1-[4-(1-hydroxy-1-methylethyl)-1,3-cyclohexadiene-1-carboxylate]-6-O-β-d-glucopyranosyl-β-d-glucopyranose (tiliaside), showed no inhibition in vitro at naturally occurring concentrations but reduced C. bombi infections of B. terrestris workers. Independent of microbiome status, tiliaside was deglycosylated during gut passage, thereby increasing its antiparasitic activity in the hindgut, the site of C. bombi infections. Conversely, unedone was first glycosylated in the midgut without influence of the microbiome to unedone-8-O-β-d-glucoside, rendering it inactive against C. bombi, but subsequently deglycosylated by the microbiome in the hindgut, restoring its activity. We therefore show that conversion of nectar metabolites by either the host or the microbiome modulates antiparasitic activity of nectar metabolites. This article is part of the theme issue 'Natural processes influencing pollinator health: from chemistry to landscapes'.

Parallel meta-transcriptome analysis reveals degradation of plant secondary metabolites by beetles and their gut symbionts

Switching to a new host plant is a driving force for divergence and speciation in herbivorous insects. This process of incorporating a novel host plant into the diet may require a number of adaptations in the insect herbivores that allow them to consume host plant tissue that may contain toxic secondary chemicals. As a result, herbivorous insects are predicted to have evolved efficient ways to detoxify major plant defenses and increase fitness by either relying on their own genomes or by recruiting other organisms such as microbial gut symbionts. In the present study we used parallel meta-transcriptomic analyses of Altica flea beetles and their gut symbionts to explore the contributions of beetle detoxification mechanisms versus detoxification by their gut consortium. We compared the gut meta-transcriptomes of two sympatric Altica species that feed exclusively on different host plant species as well as their F₁ hybrids that were fed one of the two host plant species. These comparisons revealed that gene expression patterns of Altica are dependent on both beetle species identity and diet. The community structure of gut symbionts was also dependent on the identity of the beetle species, and the gene expression patterns of the gut symbionts were significantly correlated with beetle species and plant diet. Some of the enriched genes identified in the beetles and gut symbionts are involved in the degradation of secondary metabolites produced by plants, suggesting that Altica flea beetles may use their gut microbiota to help them feed on and adapt to their host plants.

Identification and Characterization of UDP-Glycosyltransferase Genes in a Cerambycid Beetle, Pharsalia antennata Gahan, 1894 (Coleoptera: Cerambycidae)

The cerambycid beetle, Pharsalia antennata Gahan, 1894 (Coleoptera: Cerambycidae), is a wood-boring pest that spends most of its life cycle in the trunks or under the bark of trees. These distinctive biological characteristics make it likely that this beetle will encounter a number of plant defensive compounds, coupled with a broad range of host plants, possibly resulting in the overexpression or expansion of uridine diphosphate (UDP)-glycosyltransferase (UGT) genes. Here, we identified and characterized the UGT gene family in P. antennata through transcriptome data, sequence and phylogenetic analyses, and PCR and homology modeling approaches. In total, 59 transcripts encoding UGTs were identified, 34 of which harbored full-length sequences and shared high conservation with the UGTs of Anoplophora glabripennis. Of the 34 PantUGTs, only 31.78% amino acid identity was observed on average, but catalytic and sugar binding residues were highly conserved. Phylogenetic analyses revealed four Cerambycidae-specific clades, including 30 members from P. antennata. Combining the transcriptome and PCR data showed that PantUGTs had a wide tissue expression, and the majority of the genes were presented mainly in antennae or abdomens, suggesting their putative roles in olfaction and detoxification. This study provides, for the first time, information on the molecular and genetic basis of P. antennata, greatly enhancing our knowledge of the detoxification-related UGT gene family.

The Genomics and Population Genomics of the Light Brown Apple Moth, Epiphyas postvittana, an Invasive Tortricid Pest of Horticulture

Simple Summary

In this study, we produced a genomic resource for the light brown apple moth, Epiphyas postvittana, to understand the biological basis of adaptation to a high number of hosts (polyphagy) and the invasive nature of this and other lepidopteran pests. The light brown apple moth is an invasive pest of horticultural plants, with over 500 recorded plant hosts. With origins in Australia, the pest has subsequently spread to New Zealand, Hawaii, California and Europe, causing significant economic losses for fruit producers. Comparative genomic analyses with other lepidopteran genomes indicate that a high proportion of the genome is made up of repetitive sequences, with the majority of the known elements being DNA transposable elements and retrotransposons. Twenty gene families show significant expansions, including some likely to have a role in its pest status. Finally, population genomics, investigated by a RAD-tag approach, indicated likely patterns of invasion and admixture, with Californian moths most probably being derived from Australia.

Abstract

The light brown apple moth, Epiphyas postvittana is an invasive, polyphagous pest of horticultural systems around the world. With origins in Australia, the pest has subsequently spread to New Zealand, Hawaii, California and Europe, where it has been found on over 500 plants, including many horticultural crops. We have produced a genomic resource, to understand the biological basis of the polyphagous and invasive nature of this and other lepidopteran pests. The assembled genome sequence encompassed 598 Mb and has an N50 of 301.17 kb, with a BUSCO completion rate of 97.9%. Epiphyas postvittana has 34% of its assembled genome represented as repetitive sequences, with the majority of the known elements made up of longer DNA transposable elements (14.07 Mb) and retrotransposons (LINE 17.83 Mb). Of the 31,389 predicted genes, 28,714 (91.5%) were assigned to 11,438 orthogroups across the Lepidoptera, of which 945 were specific to E. postvittana. Twenty gene families showed significant expansions in E. postvittana, including some likely to have a role in its pest status, such as cytochrome p450s, glutathione-S-transferases and UDP-glucuronosyltransferases. Finally, using a RAD-tag approach, we investigated the population genomics of this pest, looking at its likely patterns of invasion.

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.

Untargeted metabolomics-based response analysis of temperature and insecticide exposure in Aedes aegypti

In this study, we utilized an untargeted NMR metabolomics approach to identify the vector response in terms of metabolic profiling after temperature and insecticide exposure in comparison with the control. Clearly, temperature and insecticide exposure cause changes in the underlying metabolism, and the NMR metabolomic profile enables a direct examination of the immediate response of the vector to cope up with these changes. The present study was designed in four parts: A-Aedes aegypti were exposed to 40 °C for one-hour, DDT-4%, malathion-5%, and deltamethrin-0.05% separately and, part B-D; one-hour exposure at 35 °C and 40 °C temperatures followed by one-hour exposure to insecticide. The resultant metabolite profiles were compared with the control. In response to temperature and insecticide exposure, several metabolites and altered pathways were identified. Citrate, maltose, lipids, Nicotinate, Choline, Pyruvate and β-hydroxybutyrate were found as important components of major biological pathways such as tri-carboxylic acid cycle, branched amino acid degradation, glycolysis/gluconeogenesis, amino acid metabolism, lipid and carbohydrate metabolism, nucleotide PRPP pathway, and phospholipid metabolism. Furthermore, the results also suggest that the changes imposed by exposure to temperature and insecticides individually, are reversed with combined exposure, thus negating the impact of each other and posing a threat to the control of Aedes-borne diseases such as dengue, chikungunya, Zika and yellow fever.

Multi-insecticide resistant malaria vectors in the field remain susceptible to malathion, despite the presence of Ace1 point mutations

Insecticide resistance in Anopheles mosquitoes is seriously threatening the success of insecticide-based malaria vector control. Surveillance of insecticide resistance in mosquito populations and identifying the underlying mechanisms enables optimisation of vector control strategies. Here, we investigated the molecular mechanisms of insecticide resistance in three Anopheles coluzzii field populations from southern Côte d’Ivoire, including Agboville, Dabou and Tiassalé. All three populations were resistant to bendiocarb, deltamethrin and DDT, but not or only very weakly resistant to malathion. The absence of malathion resistance is an unexpected result because we found the acetylcholinesterase mutation Ace1-G280S at high frequencies, which would typically confer cross-resistance to carbamates and organophosphates, including malathion. Notably, Tiassalé was the most susceptible population to malathion while being the most resistant one to the pyrethroid deltamethrin. The resistance ratio to deltamethrin between Tiassalé and the laboratory reference colony was 1,800 fold. By sequencing the transcriptome of individual mosquitoes, we found numerous cytochrome P450-dependent monooxygenases – including CYP6M2, CYP6P2, CYP6P3, CYP6P4 and CYP6P5 – overexpressed in all three field populations. This could be an indication for negative cross-resistance caused by overexpression of pyrethroid-detoxifying cytochrome P450s that may activate pro-insecticides, thereby increasing malathion susceptibility. In addition to the P450s, we found several overexpressed carboxylesterases, glutathione S-transferases and other candidates putatively involved in insecticide resistance.

Insecticide-based mosquito control has saved millions of lives from malaria and other vector-borne diseases. However, the emergence and increase of insecticide resistant Anopheles populations seriously threaten to derail malaria control programmes. Surveillance of insecticide resistance and understanding the underlying molecular mechanisms are key for choosing effective vector control strategies. Here, we characterised the degree and mechanisms of resistance in three malaria vector populations from Côte d’Ivoire. Our key finding was that these multi-insecticide resistant malaria vectors largely remained susceptible to malathion, despite the presence of a mutation in the target enzyme of this organophosphate insecticide that would typically confer resistance. Intriguingly, we found overexpression of metabolic P450 enzymes that are known to detoxify insecticides and activate pro-insecticides such as malathion. It is highly probable that, here, we observed P450-mediated negative cross-resistance for the first time in Anopheles field populations. Negative cross-resistance merits further investigation as advantage could be taken of this phenomenon in the fight against multi-resistant malaria vectors.

Transcriptomic analysis of salivary gland and proteomic analysis of oral secretion in Helicoverpa armigera under cotton plant leaves, gossypol, and tannin stresses

Gossypol and tannin are involved in important chemical defense processes in cotton plants. In this study, we used transcriptomics and proteomics to explore the changes in salivary gland functional genes and oral secretion (OS) proteins after feeding with artificial diet (containing gossypols and tannins) and cotton plant leaves. We found that dietary cotton plant leaves, gossypols and tannins exerted adverse impacts on the genes that regulated the functions of peptidase, GTPase, glycosyl hydrolases in the salivary glands of the Helicoverpa armigera (H. armigera). However, GST, UGT, hydrolases, and lipase genes were up-regulated to participate in the detoxification and digestive of H. armigera. The oral secretory proteins of H. armigera were significantly inhibited under the stress of gossypol and tannin, such as enzyme activity, but some proteins (such as PZC71358.1) were up-regulated and involved in immune and digestive functions. The combined analysis of transcriptomics and metabolomics showed a weak correlation, and the genes and proteins involved were mainly in digestive enzyme activities. Our work clarifies the deleterious physiological impacts of gossypols and tannins on H. armigera and reveals the mechanism by which H. armigera effectively mitigate the phytotoxic effects through detoxification and immune systems.

Knockdown of Uridine Diphosphate Glucosyltransferase 86Dg Enhances Susceptibility of Tribolium castaneum (Coleoptera: Tenebrionidae) to Artemisia vulgaris (Asterales: Asteraceae) Essential Oil

Uridine diphosphate glucosyltransferases (UGTs), which are phase II detoxification enzymes, are found in various organisms. These enzymes play an important role in the detoxification mechanisms of plant allelopathy and in insects. Artemisia vulgaris L. (Asterales: Asteraceae: Artemisia) essential oil has strong contact toxicity to Tribolium castaneum Herbst (Coleoptera: Tenebrionidae) larvae. However, the effect of A. vulgaris essential oil on UGTs is unclear. In this study, A. vulgaris essential oil was shown to significantly induce the expression of the TcUgt86Dg transcript. Furthermore, treatment of TcUgt86Dg-silenced individuals with A. vulgaris essential oil resulted in higher mortality than for the control individuals, indicating that TcUgt86Dg is involved in detoxification of A. vulgaris essential oil in T. castaneum. The developmental expression profile showed that the expression of TcUgt86Dg in late adults was higher than in other developmental stages. Furthermore, the expression profile in adult tissues revealed higher expression of TcUgt86Dg in the head, antenna, fat body, and accessory gland than in other tissues. These data show that TcUgt86Dg may be involved in the metabolism of exogenous toxins by T. castaneum; thus, our results have elucidated one possible mechanism of resistance to A. vulgaris essential oil and provide a theoretical basis for a control scheme for T. castaneum.

Integration of transcriptome and proteome reveals molecular mechanisms underlying stress responses of the cutworm, Spodoptera litura, exposed to different levels of lead (Pb)

Heavy metals are major environmental pollutants that affect organisms across different trophic levels. Herbivorous insects play an important role in the bioaccumulation, and eventually, biomagnification of these metals. Although effects of heavy metal stress on insects have been well-studied, the molecular mechanisms underlying their effects remain poorly understood. Here, we used the RNA-Seq profiling and isobaric tags for relative and absolute quantitation (iTRAQ) approaches to unravel these mechanisms in the polyphagous pest Spodoptera litura exposed to lead (Pb) at two different concentrations (12.5 and 100 mg Pb/kg; PbL and PbH, respectively). Altogether, 1392 and 1630 differentially expressed genes (DEGs) and 58, 114 differentially expressed proteins (DEPs) were identified in larvae exposed to PbL and PbH, respectively. After exposed to PbL, the main up-regulated genes clusters and proteins in S. litura larvae were associated with their metabolic processes, including carbohydrate, protein, and lipid metabolism, but the levels of cytochrome P450 associated with the pathway of xenobiotic biodegradation and metabolism were found to be decreased. In contrast, the main up-regulated genes clusters and proteins in larvae exposed to PbH were enriched in the metabolism of xenobiotic by cytochrome P450, drug metabolism-cytochrome P450, and other drug metabolism enzymes, while the down-regulated genes and proteins were found to be closely related to the lipid (lipase) and protein (serine protease, trypsin) metabolism and growth processes (cuticular protein). These findings indicate that S. litura larvae exposed to PbL could enhance food digestion and absorption to prioritize for growth rather than detoxification, whereas S. litura larvae exposed to PbH reduced food digestion and absorption and channelized the limited energy for detoxification rather than growth. These contrasting results explain the dose-dependent effects of heavy metal stress on insect life-history traits, wherein low levels of heavy metal stress induce stimulation, while high levels of heavy metal stress cause inhibition at the transcriptome and proteome levels.