Transcriptome analysis reveals global gene expression changes of Chilo suppressalis in response to sublethal dose of chlorantraniliprole

Sublethal acetamiprid affects reproduction, development and disrupts gene expression in Binodoxys communis

At present, understanding of neonicotinoid toxicity in arthropods remains limited. We here evaluated the lethal and sublethal effects of acetamiprid in F0 and F1 generations of Binodoxys communis using a range of sublethal concentrations. The 10% lethal concentration (LC10) and half lethal concentration (LC25) of ACE had negative effects on the B. communis survival rate, adult longevity, parasitism rate, and emergence rate, and significantly prolonged the duration of the developmental cycle. ACE also had intergenerational effects, with some biological indices affected in the F1 generation after pesticide exposure. Transcriptomic analysis demonstrated that differentially expressed genes were enriched in specific pathways including the amino acid metabolism, carbohydrate metabolism, energy metabolism, exogenous metabolism, signal transduction, and glutathione metabolism pathways. These results indicated strong contact toxicity of ACE to B. communis, which may inhibit their biological control capacity. These results improve our understanding of the toxicological mechanisms of parasitic natural enemies in response to insecticide exposure.

Isolation, identification and transcriptome analysis of triadimefon-degrading strain Enterobacter hormaechei TY18

Triadimefon, a type of triazole systemic fungicide, has been extensively used to control various fungal diseases. However, triadimefon could lead to severe environmental pollution, and even threatens human health. To eliminate triadimefon residues, a triadimefon-degrading bacterial strain TY18 was isolated from a long-term polluted site and was identified as Enterobacter hormaechei. Strain TY18 could grow well in a carbon salt medium with triadimefon as the sole nitrogen source, and could efficiently degrade triadimefon. Under triadimefon stress, a total of 430 differentially expressed genes (DEGs), including 197 up-regulated and 233 down-regulated DEGs, were identified in strain TY18 using transcriptome sequencing (RNA-Seq). Functional classification and enrichment analysis revealed that these DEGs were mainly related to amino acid transport and metabolism, carbohydrate transport and metabolism, small molecule and pyrimidine metabolism. Interestingly, the DEGs encoding monooxygenase and hydrolase activity acting on carbon-nitrogen were highly up-regulated, might be mainly responsible for the metabolism in triadimefon. Our findings in this work suggest that strain E. hormaechei TY18 could efficiently degrade triadimefon for the first time. They provide a great potential to manage triadimefon biodegradation in the environment successfully.

Comparative Transcriptome Analysis of Galeruca daurica Reveals Cold Tolerance Mechanisms

Galeruca daurica (Joannis) is a pest species with serious outbreaks in the Inner Mongolian grasslands in recent years, and its larvae and eggs are extremely cold-tolerant. To gain a deeper understanding of the molecular mechanism of its cold-tolerant stress response, we performed de novo transcriptome assembly of G. daurica via RNA-Seq and compared the differentially expressed genes (DEGs) of first- and second-instar larvae grown and developed indoors and outdoors, respectively. The results show that cold tolerance in G. daurica is associated with changes in gene expression mainly involved in the glycolysis/gluconeogenesis pathway, the fatty acid biosynthesis pathway and the production of heat shock proteins (HSPs). Compared with the control group (indoor), the genes associated with gluconeogenesis, fatty acid biosynthesis and HSP production were up-regulated in the larvae grown and developed outdoors. While the changes in these genes were related to the physiological metabolism and growth of insects, it was hypothesized that the proteins encoded by these genes play an important role in cold tolerance in insects. In addition, we also investigated the expression of genes related to the metabolic pathway of HSPs, and the results show that the HSP-related genes were significantly up-regulated in the larvae of G. daurica grown and developed outdoors compared with the indoor control group. Finally, we chose to induce significant expression differences in the Hsp70 gene (Hsp70A1, Hsp70-2 and Hsp70-3) via RNAi to further illustrate the role of heat stress proteins in cold tolerance on G. daurica larvae. The results show that separate and mixed injections of dsHSP70A1, dsHsp70-2 and dsHsp70-3 significantly reduced expression levels of the target genes in G. daurica larvae. The super-cooling point (SCP) and the body fluid freezing point (FP) of the test larvae were determined after RNAi using the thermocouple method, and it was found that silencing the Hsp70 genes significantly increased the SCP and FP of G. daurica larvae, which validated the role of heat shock proteins in the cold resistance of G. daurica larvae. Our findings provide an important theoretical basis for further excavating the key genes and proteins in response to extremely cold environments and analyzing the molecular mechanism of cold adaptation in insects in harsh environments.

Discovery of 3-Formyl-N-(un)Substituted Benzylindole Pyrimidines as an Acaricidal Agent and Their Mechanism of Action

To discover the pronounced acaricide candidate, herein, a series of 3-formyl-N-(un)substituted benzylindole pyrimidines were prepared by structural modification of indoles at the N-1 and C-3 positions via the successive Vilsmeier-Haack-Arnold (VHA), aldol condensation, and cyclization reactions. The steric structures of nine compounds were undoubtedly confirmed by X-ray single-crystallography. Against Tetranychus cinnabarinus Boisduval, compounds V-15, V-31, V-34, V-42, V-44, and V-60 exhibited promising acaricidal activity with LC50 values of 0.299-0.481 mg/mL. In particular, compound V-34 displayed 4.2 times the acaricidal activity of its precursor 6-methylindole. Scanning electron microscopy (SEM) imaging revealed that the construction of the cuticle layer of V-34-treated T. cinnabarinus was seriously destroyed. Furthermore, RNA-Seq analysis indicated that compound V-34 could regulate the homeostasis metabolism of T. cinnabarinus through arachidonic acid and linoleic acid metabolism and lysosome pathways. These results suggested that compound V-34 can be further studied as a lead acaricidal agent.

Comparative phosphoproteomics analysis provides insights into the responses of Chilo suppressalis to sublethal chlorantraniliprole exposure

BACKGROUND

Sublethal exposure to insecticides causes changes in insect behaviors and physiologies including feeding, mobility, communication, hormone homeostasis, development and fecundity, however, the underlying molecular mechanisms were largely unclear. Our previous studies revealed that sublethal chlorantraniliprole exposure disturbed the hormone homeostasis, reduced the weight and longevity and prolonged the developmental duration of Chilo suppressalis. In the present study, the potential phosphorylation modification regulation mechanisms in C. suppressalis in response to sublethal chlorantraniliprole exposure were explored using comparative and quantitative phosphoproteomics.

RESULTS

A total of 2640 phosphopeptides belonging to 1144 phosphoproteins were identified, among which 446 phosphopeptides derived from 303 unique phosphoproteins were differentially phosphorylated between the chlorantraniliprole-treated and control larvae. The phosphorylation levels of differentially phosphorylated phosphopeptides were further validated using parallel reaction monitoring (PRM). Functional classification and protein-protein interaction of the differentially phosphorylated proteins (DPPs) were analyzed. Generalized analysis of the DPPs and the differentially expressed genes (DEGs) identified in our previous study showed that sublethal chlorantraniliprole exposure significantly changed the transcription and phosphorylation levels of genes/proteins associated with carbohydrate and lipid metabolism, cytoskeleton, signal transduction, transcription, translation and post-translational modification, leading to the dysfunctions of energy metabolism, transcription regulation, protein synthesis and modification, and signal transduction in C. suppressalis. Further analysis of the phosphorylation motifs in DPPs revealed that the MAPKs, CDKs, CaMK II, PKA, PKC and CK II protein kinases might be directly responsible for the phosphoproteomics response of C. suppressalis to chlorantraniliprole treatment.

CONCLUSION

Our results provide abundant phosphorylation information for characterizing the protein modification in insects, and also provide valuable insights into the molecular mechanisms of insect post-translational modifications in response to sublethal insecticide exposure. © 2023 Society of Chemical Industry.

Sublethal effects of niclosamide on the aquatic snail Pomacea canaliculata

Pomacea canaliculata is a malignant invasive aquatic snail found worldwide, and niclosamide (NS) is one of the primary agents used for its control. NS applied to water will exist in non-lethal concentrations for some time due to degradation or water exchange, thus resulting in sublethal effects on environmental organisms. To identify sublethal effects of NS on Pomacea canaliculata, we studied the aspects of histopathology, oxygen-nitrogen ratio (R_O∶N), enzyme activity determination, and gene expression. After LC₃₀ NS treatment (0.310 g/L), many muscle fibers of the feet degenerated and some acinar vesicles of the hepatopancreas collapsed and dissolved. The oxygen-nitrogen ratio (R_O∶N) decreased significantly from 15.0494 to 11.5183, indicating that NS had changed the metabolic mode of Pomacea canaliculata and shifted it primarily to protein catabolism. Transcriptome analysis identified the sublethal effects of LC₃₀ NS on the snails at the transcriptional level. 386, 322, and 583 differentially expressed genes (DEGs) were identified in the hepatopancreas, gills, and feet, respectively. GO (Gene Ontology) functional analysis and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway annotations showed that DEGs in the hepatopancreas were mainly enriched for sugar metabolism, protein biosynthesis, immune response, and amino acid metabolism functional categories; DEGs in the gills were mainly enriched for ion transport and amino acid metabolism; DEGs in the feet were mainly enriched for transmembrane transport and inositol biosynthesis. In the future, we will perform functional validation of key genes to further explain the molecular mechanism of sublethal effects.

Risk assessment of predatory lady beetle Propylea japonica's multi-generational exposure to three non-insecticidal agrochemicals

The effects of non-insecticidal agrochemicals on pest natural predators remain largely unexplored except bees and silkworm. The herbicide quizalofop-p-ethyl (QpE), fungicide thiophanate-methyl (TM), and plant growth regulator mepiquat chloride (MC) have been extensively applied as non-insecticidal agrochemicals. Here, we systematically evaluated multiple effects of these 3 non-insecticidal agrochemicals on three generations of Propylea japonica, an important agroforestry predatory beetle, including the effects on its development, reproduction, enterobacteria, and transcriptomic response. The results showed that QpE exhibited a hormetic effect on P. japonica, thus significantly increasing the survival rate of generation 2 (F2) females, generation 3 (F3) females, and F3 males and body weight of F3 males. However, three successive generations exposed to TM and MC had no significant effect on longevity, body weight, survival rate, pre-oviposition period, and fecundity of P. japonica. Additionally, we investigated the effects of MC, TM, and QpE exposure on gene expression and gut bacterial community of F3 P. japonica. Under MC, TM, and QpE exposure, the overwhelming genes of P. japonica (99.90 %, 99.45 %, and 99.7 %) remained unaffected, respectively. Under TM and MC exposure, differentially expressed genes (DEGs) were not significantly enriched in any KEGG pathway, indicating TM and MC did not significantly affect functions of P. japonica, but under QpE exposure, the expression levels of drug metabolism-related genes were down-regulated. Although QpE treatment did not affect gut dominant bacterial community composition, it significantly increased relative abundances of detoxification metabolism-related bacteria such as Wolbachia, Pseudomonas and Burkholderia in P. japonica. However, TM and MC had no significant effect on the gut bacterial community composition and relative abundance in P. japonica. This study revealed for the first time the mechanism by which P. japonica might compensate for gene downregulation-induced detoxification metabolism decline through altering symbiotic bacteria under QpE exposure. Our findings provide reference for the rational application of non-insecticidal agrochemicals.

PI3K/Akt/CncC signaling pathway mediates the response to EPN-Bt infection in Holotrichia parallela larvae

BACKGROUND

Combining the entomopathogenic nematode (EPN), Heterorhabditis beicherriana LF strain, and Bacillus thuringiensis (Bt) HBF-18 strain is a practical strategy to manage the larvae of Holotrichia parallela Motschulsky (white grubs). However, the mechanisms underlying the larval defense response to this combined biocontrol strategy are unknown.

RESULTS

The activities of some antioxidant enzymes (SOD, POD, CAT) and some detoxifying enzymes (AChE, P-450, CarE, GST) in grubs showed an activation-inhibition trend throughout the EPN-Bt exposure time course. Eight potentially key antioxidant and detoxifying enzyme genes in response to EPN-Bt infection were identified from the midgut of grubs through RNA sequencing. After silencing CAT, CarE18, and GSTs1, the enzyme activities were significantly decreased by 30.29%, 68.80%, and 34.63%, respectively. Meanwhile, the mortality of grubs was increased by 18.40%, 46.30%, and 42.59% after exposure to EPN-Bt for 1 day. Interestingly, the PI3K/Akt signaling pathway was significantly enriched in KEGG enrichment analysis, and the expression levels of phosphatidylinositol 3-kinase (PI3K), protein kinase B (Akt), cap 'n' collar isoform-C (CncC), kelch-like ECH-associated protein 1 (Keap1), and CarE18 were all up-regulated when exposed to EPN-Bt for 1 day. Furthermore, RNAi-mediated PI3K silencing showed a similar down-regulated trend between PI3K/Akt/CncC and CarE18. Moreover, silencing PI3K rendered grubs more susceptible to EPN-Bt and accelerated symbiotic bacteria multiplication in grubs.

CONCLUSION

These results suggest that the PI3K/Akt/CncC pathway mediates the expression of CarE18 and participates in the defense response of H. parallela larvae against EPN-Bt infection. Our data provide valuable insights into the design of appropriate management strategies for this well-known agricultural pest. © 2022 Society of Chemical Industry.

Integration of miRNA and mRNA expression profiles in Asian spongy moth Lymantria dispar in response to cyantraniliprole

The Asian spongy moth, Lymantria dispar, is a worldwide forest pest that damages >500 plant species. Nowadays, chemical control is the most widely used method because of its rapidity and effectiveness, but the insecticide resistance is a growing concern for spongy moth. As important post-transcriptional regulators of gene expression, whether microRNAs (miRNAs) are involved in insecticide tolerance is little known in spongy moth. Therefore, an integrated analysis of miRNA and mRNA was performed on L. dispar larvae treated with cyantraniliprole. Compared to the control group, a total of 432 differentially expressed genes (DEGs) and 23 differentially expressed miRNAs (DEMs) were identified in L. dispar larvae under cyantraniliprole exposure. Among them, twelve DEGs encoding detoxification enzymes/proteins were further analyzed. Twenty-one genes related to insecticide tolerance were predicted by 11 DEMs, of which 25 miRNA-mRNA interactions were identified. In the miRNA-mRNA network, novel-miR-4 and mmu-miR-3475-3p were involved in the response of L. dispar to cyantraniliprole stress by regulating five genes associated with detoxification, respectively. The P450 gene CYP4C1 (c34384.graph_c0) was the only DEG related to detoxification in the network, which was regulated by novel-miR-4. The expression levels of ten DEMs were confirmed by quantitative reverse transcription PCR (RT-qPCR) and the trends were consistent with miRNA-seq. This study identified some candidate miRNAs and mRNAs related to cyantraniliprole tolerance in L. dispar, which provides valuable transcriptomic information for revealing the molecular mechanisms of insect tolerance and developing novel insecticides.

Gene expression differentials driven by mass rearing and artificial selection in black soldier fly colonies

The micro-evolutionary forces that shape genetic diversity during domestication have been assessed in many plant and animal systems. However, the impact of these processes on gene expression, and consequent functional adaptation to artificial environments, remains under-investigated. In this study, whole-transcriptome dynamics associated with the early stages of domestication of the black soldier fly (BSF), Hermetia illucens, were assessed. Differential gene expression (DGE) was evaluated in relation to (i) generational time within the cultured environment (F2 vs. F3), and (ii) two selection strategies [no artificial selective pressure (NS); and selection for greater larval mass (SEL)]. RNA-seq was conducted on 5th instar BSF larvae (n = 36), representing equal proportions of the NS (F2 = 9; F3 = 9) and SEL (F2 = 9; F3 = 9) groups. A multidimensional scaling plot revealed greater gene expression variability within the NS and F2 subgroups, while the SEL group clustered separately with lower levels of variation. Comparisons between generations revealed 898 differentially expressed genes (DEGs; FDR-corrected p < 0.05), while between selection strategies, 213 DEGs were observed (FDR-corrected p < 0.05). Enrichment analyses revealed that metabolic, developmental, and defence response processes were over-expressed in the comparison between F2 and F3 larvae, while metabolic processes were the main differentiating factor between NS and SEL lines. This illustrates the functional adaptations that occur in BSF colonies across generations due to mass rearing; as well as highlighting genic dynamics associated with artificial selection for production traits that might inform future selective breeding strategies.

Transcriptomic Analysis Reveals the Detoxification Mechanism of Chilo suppressalis in Response to the Novel Pesticide Cyproflanilide

Chilo suppressalis is one of the most damaging rice pests in China’s rice-growing regions. Chemical pesticides are the primary method for pest control; the excessive use of insecticides has resulted in pesticide resistance. C. suppressalis is highly susceptible to cyproflanilide, a novel pesticide with high efficacy. However, the acute toxicity and detoxification mechanisms remain unclear. We carried out a bioassay experiment with C. suppressalis larvae and found that the LD10, LD30 and LD50 of cyproflanilide for 3rd instar larvae was 1.7 ng/per larvae, 6.62 ng/per larvae and 16.92 ng/per larvae, respectively. Moreover, our field trial results showed that cyproflanilide had a 91.24% control efficiency against C. suppressalis. We investigated the effect of cyproflanilide (LD30) treatment on the transcriptome profiles of C. suppressalis larvae and found that 483 genes were up-regulated and 305 genes were down-regulated in response to cyproflanilide exposure, with significantly higher CYP4G90 and CYP4AU10 expression in the treatment group. The RNA interference knockdown of CYP4G90 and CYP4AU10 increased mortality by 20% and 18%, respectively, compared to the control. Our results indicate that cyproflanilide has effective insecticidal toxicological activity, and that the CYP4G90 and CYP4AU10 genes are involved in detoxification metabolism. These findings provide an insight into the toxicological basis of cyproflanilide and the means to develop efficient resistance management tools for C. suppressalis.

Comparative transcriptome analysis of Liriomyza trifolii (Burgess) and Liriomyza sativae (Blanchard) (Diptera: Agromyzidae) in response to rapid cold hardening

The ability of insets to react efficiently to fluctuation in temperature is crucial for them to survive in variable surroundings. Rapid cold hardening (RCH) is a process that increase cold tolerance in most insect species. The molecular mechanisms of RCH remain largely unknown, and whether it is associated with transcriptional changes is unclear. In this study, we compared the transcriptomes of Liriomyza trifolii and L. sativae exposed to RCH to investigate the transcript abundance due to RCH in both species. RNA-seq revealed 93,166 assembled unigenes, and 34,303 of these were annotated in the L. trifolii and L. sativae transcriptome libraries. After a 4-h treatment at 1°C (RCH) compared with control, 268 and 606 unigenes were differentially expressed in L. trifolii and L. sativae, respectively. When comparing pupae exposed to 2h cold shock directly with pupae went through 4h acclimation prior to 2h cold shock, 60 and 399 unigenes were differentially expressed in L trifolii and L sativae, respectively. Genes that were commonly expressed in both L. trifolii and L. sativae, included cytochrome P450, cuticular protein, glucose dehydrogenase, solute carrier family 22 and cationic amino acid transporter. Additionally, several pathways including galactose metabolism and peroxisome were significantly enriched during RCH. Our results show that the transcriptional response is correlated with RCH in the pupal stage of the two Liriomyza species, but more transcriptional changes were identified in L sativae than in L. trifolii.

Insecticide resistance monitoring in field populations of Chilo suppressalis Walker (Lepidoptera: Crambidae) from central China

Chilo suppressalis Walker (Lepidoptera: Crambidae) is a devastating rice crop pest in China. Chemical insecticides have been used to effectively managing C. suppressalis field populations in most of China’s agricultural regions. However, due to the intensive and extensive application of these insecticides, C. suppressalis has developed widespread resistance to many active ingredients. Thus, insecticide resistance development is a genuine concern for all crop specialists and growers. In this study, using the topical application method, we assessed the susceptibility of forty-six field populations of C. suppressalis to three insecticides in three Central Chinese provinces from 2010 to 2021. Our findings revealed that field populations of C. suppressalis built moderate to high levels of resistance to triazophos (Resistance Ratio (RR) = 41.9–250.0 folds), low to moderate levels of resistance to chlorpyrifos (RR = 9.5–95.2 folds), with the exception of the Zhijiang population in 2013 and the Xinyang population in 2015 at 4.8 folds and 3.4 folds resistance rates, respectively, despite showing susceptibility, and low and moderate levels of resistance to abamectin (RR = 4.1–53.5 folds). There were significant correlations between the activity of the detoxification enzymes (CarE) and the log LD₅₀ values of triazophos. These results should help effective insecticide resistance management strategies reduce the risk of rapid build-ups of resistance to insecticides and slow down the process of selection for insecticide resistance.

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.

Integrated miRNA and transcriptome profiling to explore the molecular mechanism of Spodoptera frugiperda larval midgut in response to azadirachtin exposure

As a destructive agricultural pest, Spodoptera frugiperda has spread worldwide in the past few years. Azadirachtin, an environmentally friendly and most promising compound, showed adverse effects, including mortality and growth inhibition, against S. frugiperda. While the effects of azadirachtin on the midgut of this pest remain to be determined. In this study, structural damage was observed in the larval midguts of S. frugiperda with azadirachtin exposure. RNA-seq on the larval midguts with different azadirachtin treatments was performed. Compared to the control group, a total of 3344 and 4759 differentially expressed genes (DEGs) were identified in the midguts with 0.1 and 0.5 μg/g azadirachtin exposure, respectively. Among them, the DEGs encoding detoxification enzymes/proteins, immune-related proteins, digestion and absorption-related proteins, and transcript factors were further analyzed. High-throughput sequencing was also used for the identification of differentially expressed microRNAs in different treatments. A total of 153 conserved miRNAs and 147 novel miRNAs were identified, of which 11 and 29 miRNAs were affected by 0.1 and 0.5 μg/g azadirachtin treatments, respectively. The integrated analysis found that 13 and 178 miRNA versus mRNA pairs were acquired in the samples with 0.1 and 0.5 μg/g azadirachtin treatments, respectively. The results of high-throughput sequencing were confirmed by real-time quantitative polymerase chain reaction (RT-qPCR). These results provide useful information for revealing the molecular mechanism of S. frugiperda larval midgut in response to azadirachtin.

Biochemical toxicity and transcriptome aberration induced by dinotefuran in Bombyx mori

Dinotefuran is a third-generation neonicotinoid pesticide and is increasingly used in agricultural production, which has adverse effects on nontarget organisms. However, the research on the impact of dinotefuran on nontarget organisms is still limited. Here the toxic effects of dinotefuran on an important economic species and a model lepidopteran insect, Bombyx mori, were investigated. Exposure to different doses of dinotefuran caused physiological disorders or death. Cytochrome P450, glutathione S-transferase, carboxylesterase, and UDP glycosyl-transferase activities were induced in the fat body at early stages after dinotefuran exposure. By contrast, only glutathione S-transferase activity was increased in the midgut. To overcome the lack of sensitivity of the biological assays at the individual organism level, RNA sequencing was performed to measure differential expressions of mRNA from silkworm larvae after dinotefuran exposure. Differential gene expression profiling revealed that various detoxification enzyme genes were significantly increased after dinotefuran exposure, which was consistent with the upregulation of the detoxifying enzyme. The global transcriptional pattern showed that the physiological responses induced by dinotefuran toxicity involved multiple cellular processes, including energy metabolism, oxidative stress, detoxification, and other fundamental physiological processes. Many metabolism processes, such as carbon metabolism, fatty acid biosynthesis, pyruvate metabolism, and the citrate cycle, were partially repressed in the midgut or fat body. Furthermore, dinotefuran significantly activated the MAPK/CREB, CncC/Keap1, PI3K/Akt, and Toll/IMD pathways. The links between physiological, biochemical toxicity and comparative transcriptomic analysis facilitated the systematic understanding of the integrated biological toxicity of dinotefuran. This study provides a holistic view of the toxicity and detoxification metabolism of dinotefuran in silkworm and other organisms.

Identification and Validation of ATP-Binding Cassette Transporters Involved in the Detoxification of Abamectin in Rice Stem Borer, Chilo suppressalis

Chilo suppressalis has developed high levels of resistance to abamectin in many areas of China, while the underline resistance mechanisms are largely unclear. ATP-binding cassette (ABC) transporters function in transporting a large diversity of substrates including insecticides and play important roles in the detoxification metabolism of insects. In this study, synergism bioassay revealed that the ABC transporters were involved in the detoxification of C. suppressalis to abamectin. Six ABC transporter genes were upregulated in C. suppressalis after abamectin exposure, among which five genes CsABCC8, CsABCE1, CsABCF1, CsABCF2, and CsABCH1 were induced in the detoxification-related tissues. In addition, the five ABC transporters were recombinantly expressed in Sf9 cells, and the cytotoxicity assay showed that the viabilities of cells expressing CsABCC8 or CsABCH1 were significantly increased when compared with the viabilities of cells expressing EGFP after abamectin, chlorantraniliprole, cyantraniliprole, fipronil, and chlorpyrifos treatment, respectively. Overexpression of CsABCE1 significantly increased the viabilities of cells to abamectin, chlorantraniliprole, deltamethrin, and indoxacarb exposure, respectively. These results suggested that CsABCC8, CsABCE1, and CsABCH1 might participate in the detoxification and transport of abamectin and several other classes of insecticides in C. suppressalis. Our study provides valuable insights into the transport-related detoxification mechanisms in C. suppressalis and other insects.

Identification and characterization of glutathione S-transferases and their potential roles in detoxification of abamectin in the rice stem borer, Chilo suppressalis

The glutathione S-transferases (GSTs) are a kind of metabolic enzymes and participate in the detoxification metabolism of xenobiotics in various organisms. In insects, GSTs play important roles in the development of insecticide resistance and antioxidant protection. The rice stem borer Chilo suppressalis is one of the most damaging pests in rice and has developed high levels of resistance to abamectin in many areas of China, whereas the potential resistance mechanisms of C suppressalis to abamectin are still unclear. In the present study, a total of 23 CsGSTs genes were identified from the C. suppressalis transcriptome and genome, including 21 cytosolic and two microsomal CsGSTs. The cytosolic CsGSTs were further classified into seven categories based on phylogenetic analysis, and their sequence characteristics and genome structures were also analyzed. Synergism study revealed that the susceptibility of C. suppressalis to abamectin was increased significantly when the CsGSTs were inhibited by diethyl maleate (DEM). Sixteen CsGSTs genes were up-regulated in C. suppressalis larvae after treatment with abamectin, among which four CsGSTs genes including CsGSTe2, CsGSTe4, CsGSTo4 and CsGSTu1 were significantly induced in the midgut and fat body tissues. These results indicated that CsGSTs were associated with the detoxification of C. suppressalis to abamectin, and CsGSTe2, CsGSTe4, CsGSTo4 and CsGSTu1 might play important roles in the insecticide detoxification or antioxidant protection in C. suppressalis. Our present study provides valuable information on C. suppressalis GSTs, and are helpful in understanding the contributions of GSTs in abamectin detoxification in C. suppressalis and other insects.

Application of transcriptomic analysis to unveil the toxicity mechanisms of fall armyworm response after exposure to sublethal chlorantraniliprole

The anthranilic diamide insecticide chlorantraniliprole is highly effective against Lepidoptera pests, but the underlying mechanisms of toxic effects of chlorantraniliprole exposures for adapting to the chemical environment are poorly known in fall armyworm (FAW), Spodoptera frugiperda (J.E.Smith). FAW being one of the most pests of maize in Latin America, suddenly appeared in China in 2019 and spread rapidly. In this study, using bioassay and transcriptomic and biochemical analyses, we comprehensively investigated gene expression changes of third instar larvae in response to different sublethal concentrations (LC10 and LC30) of chlorantraniliprole in this insect. Exposure to LC10 chlorantraniliprole (0.73 mg/L) causes 1266 differentially expressed genes (DEGs), of which 578 are up-regulated and 688 down-regulated. Exposure to LC30 (2.49 mg/L) causes differential expression of 3637 DEGs (1545 up-, 2092 down-regulated). Interestingly, the LC30 treatment led to a significant increase in the number of DEGs compared to that of the LC10, indicating a concentration effect manner. Moreover, enrichment analysis identified important DEGs belonging to specific categories, such as amino acid, carbohydrate, lipid, energy, xenobiotics metabolisms, signal transduction, and posttranslational modification pathways, and enzymes activities in enriched pathways were significantly altered at the LC10 and LC30, which matched transcriptome analysis to mediate toxic mechanisms. The DEGs encoding detoxification-related genes were identified and validated by quantitative real-time PCR (qRT-PCR), which correlated with the RNA-sequencing (RNA-seq) data. To our knowledge, these findings provide the first toxicity mechanisms for a better understanding of chlorantraniliprole action and detoxification in FAW and other insect pests at molecular level.

Effects of camptothecin on histological structures and gene expression profiles of fat bodies in Spodoptera frugiperda

Spodoptera frugiperda is a serious threat to global food production. Our previous study demonstrated that Camptothecin (CPT), a bioactive secondary metabolite from Camptotheca acuminata (Decne: Nyssaceae), exhibits adverse impact on the larval midgut of S. frugiperda and inhibits insect growth. However, effects of CPT on fat bodies of S. frugiperda larvae have not been examined yet. In the present study, we found that histological structures of fat bodies of S. frugiperda larvae were damaged in insects treated with CPT. Comparative transcriptomic analyses among different fat body samples from controls and insects treated with 1.0 and 5.0 μg/g CPT were performed. A total of 4212 and 5044 differentially expressed genes (DEGs) were identified in the samples treated with 1.0 and 5.0 μg/g CPT, respectively. Our data indicated that the pathways of detoxification, immune response, fatty acids, chitin, and hormone biosynthesis in fat bodies were affected by CPT treatments based on DEGs. These results provided a comprehensive view of the damage and gene expression changes in fat bodies of S. frugiperda after CPT exposure, which shall be useful to reveal the mechanism of CPT toxicity against S. frugiperda in future.

Transcriptome analysis reveals gene expression differences in Liriomyza trifolii exposed to combined heat and abamectin exposure

Liriomyza trifolii is an invasive pest of horticultural and vegetable crops that possesses robust competitive advantages that enable it to replace closely-related species. High temperatures often occur concomitantly with insecticide usage during L. trifolii outbreaks. In this study, we compared the transcriptomes of L. trifolii exposed to high temperature (40 °C T40), insecticide (LC₅₀ of technical grade abamectin, I50) and combined high temperature and abamectin exposure (IT5040, I50 followed by T40; and TI4050, T40 followed by I50). RNA-seq generated and revealed 44,633 unigenes with annotation data; these were compared with COG and KEGG databases for functional classification and enrichment analysis. Compared with the I50 treatment, COG classification indicated that ‘post-translational modification, protein turnover, chaperones’ was enriched in the IT5040 treatment. In the TI4050 treatment, ‘carbohydrate transport and metabolism’ was the most abundant group. The most enriched KEGG pathways in the TI4050 and IT5040 treatments were ‘longevity regulating pathway - multiple species’ and ‘protein processing in endoplasmic reticulum’, respectively. Subsequent annotation and enrichment analyses indicated that stress-related genes such as CYP450s and HSPs were differentially expressed in the I50 vs. TI4050 or I50 vs. IT5040 treatment groups. Three commercial insecticide formulations were also used to further verify the expression of selected differentially-expressed genes. This study will be conductive to consider the temperature effect on insecticide tolerance in L. trifolii, and provides a framework for improving the application efficiency of insecticides in hot weather, which will ultimately reduce the overuse of pesticides.

Flonicamid and knockdown of inward rectifier potassium channel gene CsKir2B adversely affect the feeding and development of Chilo suppressalis

BACKGROUND

The selective insecticide flonicamid shows highly insecticidal activities against piercing-sucking insects and has been widely used for the control of Hemipteran insect pests, whereas its effects on Lepidopteran insect pests remain largely unknown. Recently, inward rectifier potassium (Kir) channel has been verified to be a target of flonicamid, however, functional characterization of Lepidopteran Kir genes is still lacking.

RESULTS

Flonicamid shows no insecticidal toxicity against Chilo suppressalis larvae. However, the feeding and growth of larvae were reversibly inhibited by flonicamid (50-1200 mg L^-1 ). Flonicamid treatment also remarkably reduced and delayed the pupation and eclosion of Chilo suppressalis. Additionally, five distinct Kir channel genes (CsKir1, CsKir2A, CsKir2B, CsKir3A and CsKir3B) were cloned from Chilo suppressalis. Expression profiles analysis revealed that CsKir2A was predominately expressed in the hindgut of larvae, whereas CsKir2B had high expressions in the Malpighian tubules and hindgut. RNA interference (RNAi)-mediated knockdown of CsKir2B significantly reduced the growth and increased the mortalities of larvae, whereas silencing of CsKir2A had no obvious effects on Chilo suppressalis.

CONCLUSION

Flonicamid exhibits adverse effects on the growth and development of Chilo suppressalis. CsKir2B might be involved in the feeding behavior of Chilo suppressalis. These results provide valuable information on the effects of flonicamid on non-target insects as well as the function of insect Kir channels, and are helpful in developing new insecticide targeting insect Kir channels. © 2020 Society of Chemical Industry.

Identification of azadirachtin responsive genes in Spodoptera frugiperda larvae based on RNA-seq

The fall armyworm Spodoptera frugiperda (Lepidoptera: Noctuidae) is a polyphagous pest with 353 plant species as its hosts, including maize, sorghum, cotton, and rice. Azadirachtin is one of the most effective botanical insecticides. The effect of azadirachtin against S. frugiperda remains to be determined. Here we report strong growth inhibition of azadirachtin on S. frugiperda larvae under either 1.0 or 5.0 μg/g azadirachtin. To explore the relevant mechanisms, the larvae fed with normal artificial diet and with 1.0 μg/g azadirachtin exposure for 3 days were collected as samples for RNA-Seq. RNA-Seq on S. frugiperda larvae under different treatments identified a total of 24,153 unigenes, including 3494 novel genes, were identified. Among them, 1282 genes were affected by 1.0 μg/g azadirachtin exposure, with 672 up-regulated and 610 down-regulated. The impacted genes include 61 coding for detoxification enzymes (31 P450s, 7 GSTs, 11 CarEs, 7 UGTs and 5 ABC transporters), 31 for cuticle proteins, and several proteins involved in insect chitin and hormone biosynthesis. Our results indicated that azadirachtin could regulate the growth of S. frugiperda by affecting insect chitin and hormone biosynthesis pathway. The enhanced expression of detoxification enzymes might be related to detoxifying azadirachtin. These findings provided a foundation for further delineating the molecular mechanism of growth regulation induced by azadirachtin in S. frugiperda larvae.

Transcriptome sequencing reveals Cnaphalocrocis medinalis against baculovirus infection by oxidative stress

Cnaphalocrocis medinalis granulovirus (CnmeGV) is a potential microbial agent against the rice leaffolder. Innate immunity is essential for insects to survive pathogenic infection. Therefore, to clarify the immune response of Cnaphalocrocis medinalis to the viral colonization, the gene expression profile of C. medinalis infected with CnmeGV was constructed by RNA-seq. A total of 8,503 differentially expressed genes (DEGs) were found including 5,304 up-regulated and 3,199 down-regulated unigenes. Gene enrichment analysis indicated that these DEGs were mainly linked to protein synthesis and metabolic process as well as ribosome and virus-infection pathways. Specifically, a significantly up-regulated PiggyBac-like transposon gene was identified suggested that the enhancement of transposon activity is related to host immunity. Further, the DEGs encoding oxidative stress related genes were identified and validated by RT-qPCR. Overall, 9 antioxidant enzyme genes and 4 antioxidant protein genes were up-regulated, and the extensive glutathione S-transferase genes were down-regulated. Our results provide a basis for understanding the molecular mechanisms of baculovirus action and oxidative stress response in C. medinalis and other insects.

Induction of ER stress, antioxidant and detoxification response by sublethal doses of chlorantraniliprole in the silk gland of silkworm, Bombyx mori

Sublethal doses of chlorantraniliprole (CAP) disrupt spinning disorder in the silkworm Bombyx mori (B. mori) and cause reduced cocoon production. In the present study, we investigated the effects of trace amounts of CAP on morphology and gene expression of the B. mori silk gland, found the posterior silk gland cells were possessed of disintegrated Endoplasmic reticulum (ER), unevenly distributed chromatin after exposure to CAP (0.01 mg/L). Gene expression analysis revealed that IRE1 and ATF6 ER stress-signaling pathways were inhibited, the PERK/CncC pathway was activated. Digital gene expression (DGE) analysis showed that detoxification-related genes, antioxidant genes and genes involved in ER protein processing pathway were expressed differentially in CAP-treated silkworm larvae. Notably, the transcript levels of the detoxification-related genes (CYP4M5, CYP6AB4, GSTD3 and GSTS1) and the antioxidant genes (CAT, TPX and SOD) were significantly increased, and the expression of ER protein processing-related genes (Sec61β, Sec61γ, Sec23α and ERGIC-53) was significantly decreased after CAP exposure. The results showed that sublethal doses of CAP exposure caused ER stress, oxidative damage to the silk gland and the perturbation of protein processing in ER, thereby probably leading to abnormal growth of the silk gland and triggering the spinning failure in silkworm.

Identification and transcriptional response of ATP-binding cassette transporters to chlorantraniliprole in the rice striped stem borer, Chilo suppressalis

BACKGROUND

As the largest transporter gene family in metazoans, ATP-binding cassette (ABC) transporters regulate the efflux of a broad spectrum of substrates from the cytoplasm to the outside of the cell. In arthropods, ABC transporters are involved in phase III of the detoxification process, and play important roles in the metabolism and transport of insecticides.

RESULTS

We identified 54 ABC transporters from the genome and transcriptome of Chilo suppressalis, one of the most damaging pests of rice in China. The identified ABC transporters were classified into eight subfamilies (ABCA to ABCH) based on NCBI BLAST and phylogenetic analysis. Synergism studies showed that treatment with verapamil, a potent inhibitor of ABC transporters, resulted in significantly increased toxicity of chlorantraniliprole against C. suppressalis larvae. Among the 21 tested ABC genes, three ABC transporter genes including CsABCC8, CsABCG1C and CsABCH1 were significantly upregulated after chlorantraniliprole treatment.

CONCLUSION

ABC transporters play important roles in the detoxification and transport of chlorantraniliprole in C. suppressalis. The results from our study provide valuable information on C. suppressalis ABC transporters, and are helpful in understanding the roles of ABC transporters in chlorantraniliprole resistance mechanisms in C. suppressalis and other insect pests. © 2020 Society of Chemical Industry.

Comparative transcriptome analysis of three invasive leafminer flies provides insights into interspecific competition

Liriomyza spp. (Diptera: Agromyzidae) represent a group of economically-significant highly polyphagous pests of plants grown in field and greenhouse conditions. Liriomyza spp. share similar biological and morphological characteristics, and complex interspecific interactions have been documented among these species in various geographical regions. Where the displacement of one of these species by the other has been studied, no unique mechanisms have been identified as causing it. The impact of competitive factors (such as, insecticide tolerance, thermotolerance, and adaptability to cropping systems) may be unique to specific geographic regions of Liriomyza spp., but more research is needed to confirm these hypotheses. In this study, RNA-seq was used to determine the transcriptomes of three closely-related leafminers, e.g. L. sativae, L. trifolii, and L. huidobrensis. Over 20 Gb of clean reads were generated and assembled into unique transcriptomes, and 38,747 unigenes were annotated in different databases. In pairwise comparisons, L. trifolii and L. sativae had more up-regulated genes than L. huidobrensis. With respect to common differentially-expressed genes (Co-DEGs), the three leafminers exhibited distinct groups of highly-expressed gene clusters. When genes related to competitive factors were compared, expression patterns in L. trifolii and L. sativae were more closely related to each other than to L. huidobrensis. The data suggest that DEGs involved in competitive factors may play a key role in competition and displacement of leafminers. The divergent genes identified in this study will be valuable in revealing possible mechanisms of invasion, displacement and interspecific competition in Liriomyza spp.

Metabolic and transcriptome responses of RNAi-mediated AMPKα knockdown in Tribolium castaneum

Background

The AMP-activated protein kinase (AMPK) is an intracellular fuel sensor for lipid and glucose metabolism. In addition to the short-term regulation of metabolic enzymes by phosphorylation, AMPK may also exert long-term effects on the transcription of downstream genes through the regulation of transcription factors and coactivators. In this study, RNA interference (RNAi) was conducted to investigate the effects of knockdown of TcAMPKα on lipid and carbohydrate metabolism in the red flour beetle, Tribolium castaneum, and the transcriptome profiles of dsTcAMPKα-injected and dsEGFP-injected beetles under normal conditions were compared by RNA-sequencing.

Results

RNAi-mediated suppression of TcAMPKα increased whole-body triglyceride (TG) level and the ratio between glucose and trehalose, as was confirmed by in vivo treatment with the AMPK-activating compound, 5-Aminoimidazole-4-carboxamide1-β-D-ribofuranoside (AICAR). A total of 1184 differentially expressed genes (DEGs) were identified between dsTcAMPKα-injected and dsEGFP-injected beetles. These include genes involved in lipid and carbohydrate metabolism as well as insulin/insulin-like growth factor signaling (IIS). Real-time quantitative polymerase chain reaction analysis confirmed the differential expression of selected genes. Interestingly, metabolism-related transcription factors such as sterol regulatory element-binding protein 1 (SREBP1) and carbohydrate response element-binding protein (ChREBP) were also significantly upregulated in dsTcAMPKα-injected beetles.

Conclusions

AMPK plays a critical role in the regulation of beetle metabolism. The findings of DEGs involved in lipid and carbohydrate metabolism provide valuable insight into the role of AMPK signaling in the transcriptional regulation of insect metabolism.

Physiological adaptations to sugar-mimic alkaloids: Insights from Bombyx mori for long-term adaption and short-term response

Insects evolved adaptive plasticity to minimize the effects of the chemical defenses of their host plants. Nevertheless, the expressional response and adaptation of phytophagous specialists for long-term adaption and short-term response to host phytochemicals remains largely unexplored. The mulberry (Morus alba)-silkworm (Bombyx mori) interaction is an old and well-known model of plant-insect interaction. In this study, we examined the long-term adaption and short-term response of the mulberry-specialist silkworm to two sugar-mimic alkaloids in mulberry: the commonly encountered 1-deoxynojirimycin (1-DNJ) and occasionally encountered 1,4-dideoxy-1,4-imino-D-arabinitol (D-AB1), respectively. Global transcriptional patterns revealed that the physiological responses induced by the selective expression of genes involved in manifold cellular processes, including detoxification networks, canonical digestion processes, target enzymes, and other fundamental physiological processes, were crucial for regulating metabolic homeostasis. Comparative network analysis of the effects of exposure to D-AB1 and 1-DNJ supported the contention that B. mori produced similar and specific trajectories of changed gene expression in response to different sugar-mimic alkaloids. D-AB1 elicited a substantial proportion of downregulated genes relating to carbohydrate metabolism, catabolic process, lipid metabolism, and glycan biosynthesis and metabolism. This study dramatically expands our knowledge of the physiological adaptations to dietary sugar-mimic alkaloid intake and uncovered both metabolic evolutionarily responses and unique adaptive mechanisms previously unknown in insects.

Transcriptome analysis reveals gene expression changes of the fat body of silkworm (Bombyx mori L.) in response to selenium treatment

A comparative transcriptome analysis was conducted to investigate the gene expression changes in the fat body of silkworm after treatment with different concentrations (50 μM and 200 μM) of selenium (Se). 912 differential expression genes (DEGs) (371 up-regulated and 541 down-regulated) and 1420 DEGs (1078 up-regulated and 342 down-regulated) were identified in silkworm fat body treated with 50 μM and 200 μM of Se, respectively. In case of 50 μM group, DEGs were mainly enriched in the peroxisome pathway and fatty acid metabolism pathway, and later were associated with antioxidant defense and nutrition regulation. After 200 μM Se-treatment, DEGs were mainly located in the glycerolipid metabolism and arachidonic acid metabolism pathways, which further encoded detoxification related genes. Furthermore, 32 candidate DEGs from these pathways had been selected to confirm the RNA-seq data. Among these DEGs, 14 genes were up-regulated in the 50 μM Se-treated group (only three genes in the 200 μM Se-treated group) which were involved in lipid metabolism and antioxidant defense, and 13 up-regulated genes (only two genes were up-regulated in the 50 μM Se-treated group) were involved in detoxification of the 200 μM Se-treated group. These changes showed that lower concentration of Se could regulate the nutrition and promote antioxidation pathways; whereas, high levels of Se promoted the detoxification of silkworm. These findings can be helpful to understand the possible mechanisms of Se action and detoxification in silkworm and other insects.

Molecular Insights into the Insensitivity of Lepidopteran Pests to Cycloxaprid

Cycloxaprid (CYC) is effective in the control of hemipteran pests, but its bioactivity against lepidopteran pests is still unclear. Here, the bioactivity of CYC against lepidopteran pests was found to be much worse than that against hemipteran insects. To reveal the mechanism, the transcriptomes of CYC-treated and untreated Ostrinia furnacalis larvae were compared. Among the top 20 differentially expressed genes, 11 encode proteins involved in cuticle formation, while only one encodes a detoxifying enzyme. Thus, the cuticle appears to be important for the insensitivity of O. furnacalis to CYC. A pretreatment of O. furnacalis larvae with methoprene enhanced the bioactivity of CYC by 1.12-fold. Moreover, mixtures of CYC with graphene oxide increased the bioactivity of CYC by 1.88-fold. Because lepidopteran and hemipteran insects often harm crops at the same time, the work can help make full use of CYC and reduce the environmental impacts of using multiple pesticides.

Sublethal effects of chlorantraniliprole on molting hormone levels and mRNA expressions of three Halloween genes in the rice stem borer, Chilo suppressalis

While sublethal effects of insecticide on insect development have been widely studied, the underlying mechanisms remain elusive. Our previous studies revealed that sublethal concentrations of chlorantraniliprole significantly increased the juvenile hormone levels and resulted in both prolonged developmental time and reduced fecundity in Chilo suppressalis. In the present study, we evaluated the sublethal effects of chlorantraniliprole on molting hormone (MH) levels and mRNA expressions of three Halloween genes including CsCYP307A1, CsCYP306A1 and CsCYP314A1 in C. suppressalis. The results showed that the MH levels in different developmental stages of C. suppressalis were decreased after exposure to LC₁₀ and LC₃₀ of chlorantraniliprole. However, analysis of temporal expression profiles revealed that the mRNA levels of three Halloween genes were not closely correlated with the ecdysteroid titers in C. suppressalis. Notably, the transcript levels of CsCYP307A1, CsCYP306A1 and CsCYP314A1 were induced after treatment with sublethal concentrations of chlorantraniliprole in specific developmental stages. These results indicated that chlorantraniliprole had adverse effects on insect MH biosynthesis, and in addition to the involvement in MH biosynthesis, CsCYP307A1, CsCYP306A1 and CsCYP314A1 may also play important roles in the detoxification metabolism of chlorantraniliprole in C. suppressalis.