Comparative analysis of the UDP-glycosyltransferase multigene family in insects

Genome-wide phylogenetic analysis and expansion of gene families involved in detoxification in Smittia aterrima (Meigen)and Smittia pratorum (Goetghebuer) (Diptera, Chironomidae)

Smittia aterrima (Meigen, 1818) and Smittia pratorum (Goetghebuer, 1927) are important indicator insects for aquatic environments, showing extensive tolerance to the environment. However, the genome-wide phylogenetic relationships and characteristics of the detoxification mechanisms in S. aterrima and S. pratorum remain unclear. Based on the genomes of the two species obtained in our preliminary studies and nine genomes from the NCBI database, we found that chironomids diverged from other mosquitoes approximately 200 million years ago (MYA), and S. aterrima and S. pratorum diverged about 30 MYA according to phylogenetic analysis. Gene family evolution analysis showed significant expansion of 43 and 15 gene families in S. aterrima and S. pratorum, respectively, particularly those related to detoxification pathways. Positive selection analysis reveals that genes under positive selection are crucial for promoting environmental adaptation. Additionally, the detoxification-associated gene families including Cytochrome P450 (CYP), Glutathione S-transferases (GST), ATP-binding cassette (ABC), carboxylesterase (CCE), and UDP-glucuronosyltransferase (UGT) were annotated. Our analysis results show that these five detoxification gene families have significantly expanded in the chironomid genomes. This study highlights the genome evolution of chironomids and their responses to mechanisms of tolerance to environmental challenges.

The polyvalent sequestration ability of an economically important beetle

Many specialized herbivorous insects sequester single classes of toxic secondary metabolites from their host plants as protection against natural enemies. If and how herbivores can use multiple classes of plant toxins across the large chemical diversity of plants for self-protection is unknown. We show that the polyphagous adults of the beetle Diabrotica virgifera are capable of selectively accumulating benzoxazinoids, cucurbitacins, and glucosinolates but not cyanogenic glycosides. Female beetles transfer the sequestered defense metabolites into their eggs, protecting them against generalist predators. Eggs containing a mixture of toxins are better protected than eggs with individual toxins. This work shows how herbivores can exploit plant chemical diversity to their own benefit as a novel adaptive mechanism that contributes to the structuring of multitrophic interaction networks.

Knockout of two uridine diphosphate-glycosyltransferase genes increases the susceptibility of Spodoptera litura to Bacillus thuringiensis toxins

Uridine diphosphate-glycosyltransferases (UGTs) catalyze sugar conjugation of endogenous and exogenous molecules in insects. In this study, 45 putative UGT genes in 11 families were identified from the genome of S. litura. Exposure to Bt toxins in 5th-instar larvae of the WT strain led to a significant upregulation of midgut UGT40 expression, particularly of SlUGT40D20, SlUGT40D22, and SlUGT40F25. This upregulation was not observed following exposure to chemical pesticides. Knockout of the UGT genes SlUGT40D20 and SlUGT40D22 in S. litura (mutant strains SlUGT40D20-KO and SlUGT40D22-KO) via CRISPR/Cas9-mediated mutagenesis increased susceptibility of S. litura to Bacillus thuringiensis (Bt) insecticidal proteins. However, in comparison with the wild-type (WT) strain, the mutants did not change susceptibility to chemical pesticides. Observations of 5th-instar larval midgut by electron microscopy revealed severe damage to the midgut epithelium caused by Cry1Ac toxin at 10 μg/g in the SlUGT40D20-KO strain compared to the WT. SDS-PAGE and LC MS/MS analyses identified a specific protein band corresponding to putative proteoglycans in the peritrophic matrix of the WT strain, which was absent in the SlUGT40D20-KO strain. Our study suggests an inverse correlation between expression of some UGTs and the susceptibility of S. litura larvae to some Bt toxins.

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.

Trait responses, nonconsumptive effects, and the physiological basis of Helicoverpa armigera to bat predation risk

Predation reduces the population density of prey, affecting its fitness and population dynamics. Few studies have connected trait changes with fitness consequences in prey and the molecular basis and metabolic mechanisms of such changes in bat-insect systems. This study focuses on the responses of Helicoverpa armigera to different predation risks, focusing on echolocating bats and their calls. Substantial modifications were observed in the nocturnal and diurnal activities of H. armigera under predation risk, with enhanced evasion behaviors. Accelerated development and decreased fitness were observed under predation risks. Transcriptomic and metabolomic analyses indicated that exposure to bats induced the upregulation of amino acid metabolism- and antioxidant pathway-related genes, reflecting shifts in resource utilization in response to oxidative stress. Exposure to bat predation risks enhanced the activity of DNA damage repair pathways and suppressed energy metabolism, contributing to the observed trait changes and fitness decreases. The current results underscore the complex adaptive strategies that prey species evolve in response to predation risk, enhancing our understanding of the predator-prey dynamic and offering valuable insights for innovative and ecologically informed pest management strategies.

Physiological effects of field concentrations and sublethal concentrations of sulfoxaflor on Apis mellifera

BACKGROUND

Bees (Apis mellifera), as important pollinators of agricultural crops, are at risk when pesticides are used. Sulfoxaflor is a new insecticide which acts on the nicotinic acetylcholine receptor (nAChR) in a similar way to neonicotinoids. The goal of this study is to evaluate the toxicity of sulfoxaflor and its effect on the A. mellifera exposure.

RESULTS

Initially, developmental indicators such as larval survival, pupation, and eclosion were inhibited by 5.0 mg/L (field concentration) sulfoxaflor. In the pupal stage, fat content was significantly increased, while the glycogen content decreased. In addition, A. mellifera heads were treated with 2.0 mg/L (sublethal concentration) of sulfoxaflor and analyzed by RNA sequencing. The transcriptome results indicated that 2.0 mg/L amounts of sulfoxaflor have adverse effects on the immune, digestive, and nervous systems. Sulfoxaflor down-regulated the expression of many genes involved in immunity, detoxification, the myosin cytoskeleton, sensory neurons, and odor-binding proteins.

CONCLUSION

Field concentration and sublethal concentration were used for the combined analysis of honeybees. The effect of sublethal concentration of sulfoxaflor on honeybees was studied for the first time from the perspective of transcriptome sequencing of honeybee head. A preliminary study was carried out on the stress of sulfoxaflor at sublethal concentration on honeybee workers, which has certain research significance and can provide theoretical basis for the use of sulfoxaflor in the field environment. © 2024 Society of Chemical Industry.

Viral- and fungal-mediated behavioral manipulation of hosts: summit disease

Summit disease, in which infected hosts seek heights (gravitropism), first noted in modern times by nineteenth-century naturalists, has been shown to be induced by disparate pathogens ranging from viruses to fungi. Infection results in dramatic changes in normal activity patterns, and such parasite manipulation of host behaviors suggests a strong selection for convergent outcomes albeit evolved via widely divergent mechanisms. The two best-studied examples involve a subset of viral and fungal pathogens of insects that induce "summiting" in infected hosts. Summiting presumably functions as a means for increasing the dispersal of the pathogen, thus significantly increasing fitness. Here, we review current advances in our understanding of viral- and fungal-induced summit disease and the host behavioral manipulation involved. Viral genes implicated in this process include a host hormone-targeting ecdysteroid UDP-glucosyltransferase (apparently essential for mediating summit disease induced by some viruses but not all) and a protein tyrosine phosphatase, with light dependance implicated. For summit disease-causing fungi, though much remains obscure, targeting of molting, circadian rhythms, sleep, and responses to light patterns appear involved. Targeting of host neuronal pathways by summit-inducing fungi also appears to involve the production of effector molecules and secondary metabolites that affect host muscular, immune, and/or neurological processes. It is hypothesized that host brain structures, particularly Mushroom Bodies (no relation to the fungus itself), important for olfactory association learning and control of locomotor activity, are critical targets for mediating summiting during infection. This phenomenon expands the diversity of microbial pathogen-interactions and host dynamics. KEY POINTS: • Summit disease or height seeking (gravitropism) results from viral and fungal pathogens manipulating insect host behaviors presumably to increase pathogen dispersal. • Insect baculoviruses and select fungal pathogens exhibit convergent evolution in host behavioral manipulation but use disparate molecular mechanisms. • Targets for affecting host behavior include manipulation of host hormones, feeding, locomotion, and immune, circadian, and neurological pathways.

A Chromosome-Scale and Annotated Reference Genome Assembly of Plecia longiforceps Duda, 1934 (Diptera: Bibionidae)

Urbanization is a leading factor effecting global biodiversity, driving rapid evolutionary processes in the local biota. Species that adapt and proliferate in city environments can become pests, with human activities facilitating their dispersal and excessive outbreaks. Here we present the first genome data of Plecia longiforceps, a lovebug pest in Eastern Asia with intensive aggregations recently occurring in the Seoul Metropolitan Area of Korea. PacBio HiFi and ONT Pore-C sequencing data were used to construct a highly continuous assembly with a total size of 707 Mb and 8 major pseudochromosomes, its integrity supported by the N50 length of 98.1 Mb and 96.8% BUSCO completeness. Structural and functional annotation using transcriptome data and ab initio predictions revealed a high proportion (69.3%) of repeat sequences, and synteny analysis with Bibio marci showed high levels of genomic collinearity. The genome will serve as an essential resource for both population genomics and molecular research on lovebug dispersal and outbreaks, and also implement studies on the eco-evolutionary processes of insects in urbanizing habitats.

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.

Knocking Out of UDP-Glycosyltransferase Gene UGT2B10 via CRISPR/Cas9 in Helicoverpa armigera Reveals Its Function in Detoxification of Insecticides

The role of insect UDP-glycosyltransferases (UGTs) in the detoxification of insecticides has rarely been reported. A UGT gene UGT2B10 was previously found overexpressed in a fenvalerate-resistant strain of Helicoverpa armigera. Herein, UGT2B10 was cloned, and its involvement in insecticide detoxification was investigated. UGT2B10 was highly expressed in the larvae, mainly in the fat body and midgut. Treatment with UGT inhibitors 5-nitrouracil and sulfinpyrazone significantly enhanced the fenvalerate toxicity. Knocking down UGT2B10 by RNAi significantly increased the larvae mortality by 17.89%. UGT2B10 was further knocked out by CRISPR/Cas9, and a homozygous strain (HD-dUGT2B10) with a C-base deletion at exon 2 was obtained. The sensitivity of HD-dUGT2B10 to fenvalerate, deltamethrin, cyantraniliprole, acetamiprid, and lufenuron increased significantly, with sensitivity index increased 2.523-, 2.544-, 2.250-, 2.473-, and 3.556-fold, respectively. These results suggested that UGT2B10 was involved in the detoxification of H. armigera to insecticides mentioned above, shedding light upon further understanding of the detoxification mechanisms of insecticides by insect UGTs.

Insights into the Detoxification of Spruce Monoterpenes by the Eurasian Spruce Bark Beetle

Plant defence mechanisms, including physical barriers like toughened bark and chemical defences like allelochemicals, are essential for protecting them against pests. Trees allocate non-structural carbohydrates (NSCs) to produce secondary metabolites like monoterpenes, which increase during biotic stress to fend off pests like the Eurasian spruce bark beetle, ESBB (Ips typographus). Despite these defences, the ESBB infests Norway spruce, causing significant ecological damage by exploiting weakened trees and using pheromones for aggregation. However, the mechanism of sensing and resistance towards host allelochemicals in ESBB is poorly understood. We hypothesised that the exposure of ESBB to spruce allelochemicals, especially monoterpenes, leads to an upsurge in the important detoxification genes like P450s, GSTs, UGTs, and transporters, and at the same time, genes responsible for development must be compromised. The current study demonstrates that exposure to monoterpenes like R-limonene and sabiene effectively elevated detoxification enzyme activities. The differential gene expression (DGE) analysis revealed 294 differentially expressed (DE) detoxification genes in response to R-limonene and 426 DE detoxification genes in response to sabiene treatments, with 209 common genes between the treatments. Amongst these, genes from the cytochrome P450 family 4 and 6 genes (CP4 and CP6), esterases, glutathione S-transferases family 1 (GSTT1), UDP-glucuronosyltransferase 2B genes (UDB), and glucose synthesis-related dehydrogenases were highly upregulated. We further validated 19 genes using RT-qPCR. Additionally, we observed similar high expression levels of detoxification genes across different monoterpene treatments, including myrcene and α-pinene, suggesting a conserved detoxification mechanism in ESBB, which demands further investigation. These findings highlight the potential for molecular target-based beetle management strategies targeting these key detoxification genes.

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.

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

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

Uridine diphosphate (UDP)-glycosyltransferases (UGTs) are associated with insecticide resistance in the major malaria vectors Anopheles gambiae s.l. and Anopheles funestus

Malaria remains one of the highest causes of morbidity and mortality, with 249 million cases and over 608,000 deaths in 2022. Insecticides, which target the Anopheles mosquito vector, are the primary method to control malaria. The widespread nature of resistance to the most important insecticide class, the pyrethroids, threatens the control of this disease. To reverse the stall in malaria control there is urgent need for new vector control tools, which necessitates understanding the molecular basis of pyrethroid resistance. In this study we utilised multi-omics data to identify uridine-diphosphate (UDP)-glycosyltransferases (UGTs) potentially involved in resistance across multiple Anopheles species. Phylogenetic analysis identifies sequence similarities between Anopheline UGTs and those involved in agricultural pesticide resistance to pyrethroids, pyrroles and spinosyns. Expression of five UGTs was characterised in An. gambiae and An. coluzzii to determine constitutive over-expression, induction, and tissue specificity. Furthermore, a UGT inhibitor, sulfinpyrazone, restored susceptibility to pyrethroids and DDT in An. gambiae, An. coluzzii, An. arabiensis and An. funestus, the major African malaria vectors. Taken together, this study provides clear association of UGTs with pyrethroid resistance as well as highlighting the potential use of sulfinpyrazone as a novel synergist for vector control.

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).

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.

Brain Gene Expression of Foraging Behavior and Social Environment in Ceratina calcarata

Rudimentary social systems have the potential to both advance our understanding of how complex sociality may have evolved and our understanding of how changes in social environment may influence gene expression and cooperation. Recently, studies of primitively social Hymenoptera have greatly expanded empirical evidence for the role of social environment in shaping behavior and gene expression. Here, we compare brain gene expression profiles of foragers across social contexts in the small carpenter bee, Ceratina calcarata. We conducted experimental manipulations of field colonies to examine gene expression profiles among social contexts including foraging mothers, regular daughters, and worker-like dwarf eldest daughters in the presence and absence of mother. Our analysis found significant differences in gene expression associated with female age, reproductive status, and social environment, including circadian clock gene dyw, hexamerin, and genes involved in the regulation of juvenile hormone and chemical communication. We also found that candidate genes differentially expressed in our study were also associated with division of labor, including foraging, in other primitively and advanced eusocial insects. Our results offer evidence for the role of the regulation of key developmental hormones and circadian rhythms in producing cooperative behavior in rudimentary insect societies.

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.