Identification of transcriptome and fluralaner responsive genes in the common cutworm Spodoptera litura Fabricius, based on RNA-seq

Rapid Adaptation and Interspecific Introgression in the North American Crop Pest Helicoverpa zea

Insect crop pests threaten global food security. This threat is amplified through the spread of nonnative species and through adaptation of native pests to control measures. Adaptations such as pesticide resistance can result from selection on variation within a population, or through gene flow from another population. We investigate these processes in an economically important noctuid crop pest, Helicoverpa zea, which has evolved resistance to a wide range of pesticides. Its sister species Helicoverpa armigera, first detected as an invasive species in Brazil in 2013, introduced the pyrethroid-resistance gene CYP337B3 to South American H. zea via adaptive introgression. To understand whether this could contribute to pesticide resistance in North America, we sequenced 237 H. zea genomes across 10 sample sites. We report H. armigera introgression into the North American H. zea population. Two individuals sampled in Texas in 2019 carry H. armigera haplotypes in a 4 Mbp region containing CYP337B3. Next, we identify signatures of selection in the panmictic population of nonadmixed H. zea, identifying a selective sweep at a second cytochrome P450 gene: CYP333B3. We estimate that its derived allele conferred a ∼5% fitness advantage and show that this estimate explains independently observed rare nonsynonymous CYP333B3 mutations approaching fixation over a ∼20-year period. We also detect putative signatures of selection at a kinesin gene associated with Bt resistance. Overall, we document two mechanisms of rapid adaptation: the introduction of fitness-enhancing alleles through interspecific introgression, and selection on intraspecific variation.

Analysis of Genes Associated with Feeding Preference and Detoxification in Various Developmental Stages of Aglais urticae

Herbivorous insects and host plants have developed a close and complex relationship over a long period of co-evolution. Some plants provide nutrients for insects, but plants' secondary metabolites also influence their growth and development. Urtica cannabina roots and leaves are poisonous, yet Aglais urticae larvae feed on them, so we aimed to clarify the mechanism enabling this interaction. At present, studies on the detoxification mechanism of the A. urticae are rare. In our study, first, we used the A. urticae larval odor selection behavior bioassay and choice feeding preference assay to analyze the feeding preferences of A. urticae on its host plant, U. cannabina. Next, we used transcriptome sequencing to obtain the unigenes annotated and classified by various databases, such as KEGG and GO. In this study, we found that U. cannabina could attract A. urticae larvae to feed via scent, and the feeding preference assay confirmed that larvae preferred U. cannabina leaves over three other plants: Cirsium japonicum, Cannabis sativa, and Arctium lappa. The activity of detoxifying enzymes GST and CarE changed in larvae that had consumed U. cannabina. Furthermore, through transcriptomic sequencing analysis, 77,624 unigenes were assembled from raw reads. The numbers of differentially expressed genes were calculated using pairwise comparisons of all life stages; the expression of detoxification enzyme genes was substantially higher in larvae than in the pupal and adult stages. Finally, we identified and summarized 34 genes associated with detoxification enzymes, such as UDP-glucose 4-epimerase gene, 5 Glutathione S-transferase genes, 4 Carboxylesterase genes, 4 Cytochrome P450 genes, 10 ATP-binding cassette genes, 4 Superoxide dismutase, and Peroxidase. Moreover, we identified 28 genes associated with the development of A. urticae. The qRT-PCR results were nearly consistent with the transcriptomic data, showing an increased expression level of four genes in larvae. Taken together, this study examines the correlation between A. urticae and host plants U. cannabina, uncovering a pronounced preference for A. urticae larvae toward host plants. Consistent with RNA-seq, we investigated the mechanism of A. urticae's interaction with host plants and identified detoxification-related genes. The present study provides theoretical support for studying insect adaptation mechanisms and biological control.

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

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

Genomic signals of local adaptation across climatically heterogenous habitats in an invasive tropical fruit fly (Bactrocera tryoni)

Local adaptation plays a key role in the successful establishment of pest populations in new environments by enabling them to tolerate novel biotic and abiotic conditions experienced outside their native range. However, the genomic underpinnings of such adaptive responses remain unclear, especially for agriculturally important pests. We investigated population genomic signatures in the tropical/subtropical Queensland fruit fly, Bactrocera tryoni, which has an expanded range encompassing temperate and arid zones in Australia, and tropical zones in the Pacific Islands. Using reduced representation sequencing data from 28 populations, we detected allele frequency shifts associated with the native/invasive status of populations and identified environmental factors that have likely driven population differentiation. We also determined that precipitation, temperature, and geographic variables explain allelic shifts across the distribution range of B. tryoni. We found spatial heterogeneity in signatures of local adaptation across various climatic conditions in invaded areas. Specifically, disjunct invasive populations in the tropical Pacific Islands and arid zones of Australia were characterised by multiple significantly differentiated single nucleotide polymorphisms (SNPs), some of which were associated with genes with well-understood function in environmental stress (e.g., heat and desiccation) response. However, invasive populations in southeast Australian temperate zones showed higher gene flow with the native range and lacked a strong local adaptive signal. These results suggest that population connectivity with the native range has differentially affected local adaptive patterns in different invasive populations. Overall, our findings provide insights into the evolutionary underpinnings of invasion success of an important horticultural pest in climatically distinct environments.

Inducible Cytochrome P450s in the Fat Body and Malpighian Tubules of the Polyphagous Pests of Spodoptera litura Confer Xenobiotic Tolerance

Cytochrome P450 plays vital roles in detoxifying xenobiotics. In this study, SlCYP340A and SlCYP340L expression in the Spodoptera litura fat body and SlCYP332A1, SlCYP6AB12, SlCYP6AB58, SlCYP6AB59, and SlCYP6AN4 expression in the Malpighian tubules were significantly upregulated after cyantraniliprole exposure, and SlCYP6AB58 and SlCYP6AB59 expression levels were simultaneously increased in the Malpighian tubules after gossypol treatment. Drosophila ectopically expressing candidate P450 genes showed that SlCYP332A1, SlCYP6AB12, SlCYP6AB59, SlCYP6AN4, and SlCYP340A conferred cyantraniliprole tolerance. The overexpression of SlCYP6AB58 and SlCYP6AB59 in Drosophila increased the number of eggs laid under the gossypol treatment. Moreover, the knockdown of SlCYP332A1, SlCYP6AB12, SlCYP6AB59, SlCYP6AN4, and SlCYP340A increased S. litura mortality under the cyantraniliprole treatment. Homology modeling and molecular docking results suggested that candidate P450 has the potential to bind with cyantraniliprole. These results indicate that the CYP3 and CYP4 genes participate in cyantraniliprole detoxification and that SlCYP6AB59 may be simultaneously involved in the gossypol tolerance of S. litura.

Physiological Function of RDL1 and RDL2 Subunits of the Ionotropic GABA Receptor in the Spodoptera litura with the CRISPR/Cas9 System In Vivo

In insect ionotropic γ-aminobutyric acid receptor (iGABAR) subunits, only resistance to dieldrin (RDL) can be individually and functionally expressed in vitro. In lepidopteran, two to three RDL subtypes are identified; however, their physiological roles have not been distinguished in vivo. In this study, SlRdl1 and SlRdl2 of S. litura were individually knocked out using CRISPR/Cas9, respectively. The mortality and larval and pupal duration of KOSlRdl1 and KOSlRdl2 were increased. The flight time and distance were increased by 43.30%-80.66% and 58.96%-198.22%, respectively, in KOSlRdl1. The GABA-induced current was significantly decreased by 53.57%-74.28% and 46.91%-63.34% in the ventral nerve cord, and the GABA titer was significantly reduced by 17.65%-28.05% and 19.85%-42.46% in KOSlRdl1 and KOSlRdl2, respectively. In conclusion, SlRdl1 and SlRdl2 are necessary for the transmission of GABA-induced neural signals; however, only SlRdl1 could regulate the flight capability of S. litura. Our results provided a new avenue to study lepidopteran iGABARs.

Synergistic effects of chlorantraniliprole and camptothecin on physiological impairments, histopathological, biochemical changes, and genes responses in the larvae midgut of Spodoptera frugiperda

Spodoptera frugiperda is an economically important agricultural pest and poses a serious threat to food security globally. Its management is gravely challenged by its high polyphagous nature, strong migratory ability, and massive fecundity. Chlorantraniliprole (CHL) is widely utilized in controlling S. frugiperda, its intensive application and over-reliance pose adverse health risks, development of resistance, toxicity to beneficial insects, natural enemies, and environmental contamination. To address S. frugiperda resistance to CHL and its inherent challenges, this study explores the synergistic effects of camptothecin (CPT) with CHL in its management. The binary mixed adversely induced the larvae weight and mortality when compared to single-treated. CHL + CPT (1:20 mg/L) had the highest larvae mortality of (73.80 %) with a high antagonistic factor (0.90), while (1:10 mg/L) with (66.10%) mortality exhibited a high synergistic factor (1.43). Further, CHL + CPT (1:10 mg/L) considerably altered the midgut epithelial cell, peritrophic membrane, microvilli, basement membrane, and regenerative cells. For biochemical analysis, CHL + CPT (1:10 mg/L) significantly decreased glutathione-S-transferase (1-chloro-2,4-dinitrobenzene CDNB) and cytochrome P450 (7-ethoxycoumarin O-deethylation) activities in the midgut in a dose and time dependent manner. Based on RNA-Seq analysis, a total of 4,373 differentially expressed genes (DEGs) were identified from the three treatments. CPT vs CK (Control) had 1694 (968 up-, 726 down-regulated), CHL vs CK with 1771 (978 up-, 793 down-regulated), and CHL + CPT vs CK had 908 (394 up-, 514 down-regulated) DEGs. The enrichment analysis disclosed significant pathways such as metabolism of xenobiotics by cytochrome P450, glutathione metabolism, TOLL and IMD (Immune Deficiency) signaling pathway, longevity regulating pathway. This study provides basis to expatiate on the molecular toxicological mechanism of CHL + CPT in management of fall armyworm.

Inducible Gut-Specific Carboxylesterase SlCOE030 in Polyphagous Pests of Spodoptera litura Conferring Tolerance between Nicotine and Cyantraniliprole

Insecticides tolerance in herbivorous arthropods is associated with preadaptation to host plant allelochemicals. However, how plant secondary metabolites activate detoxifying metabolic genes to develop tolerance remains unclear. Herein, the tolerance of Spodoptera litura larvae to cyantraniliprole was increased after nicotine exposure. An S. litura α esterase, SlCOE030, was predominantly expressed in the midgut and induced after exposure to cyantraniliprole, nicotine, and cyantraniliprole plus nicotine. Drosophila melanogaster with ectopically overexpressed SlCOE030 enhanced cyantraniliprole and nicotine tolerance by 4.91- and 2.12-fold, respectively. Compared to UAS-SlCOE030 and Esg-GAL4 lines, the Esg > SlCOE030 line laid more eggs after nicotine exposure. SlCOE030 knockdown decreased the sensitivity of nicotine-treated S. litura larvae to cyantraniliprole. Metabolism assays indicated that recombinant SlCOE030 protein metabolizes cyantraniliprole. Homology modeling and molecular docking analysis demonstrated that SlCOE030 exhibits effective affinities for cyantraniliprole and nicotine. Thus, insect CarEs may result in the development of cross-tolerance between synthetic insecticides and plant secondary metabolites.

De Novo Transcriptome Assembly and Analysis of Longevity Genes Using Subterranean Termite (Reticulitermes chinensis) Castes

The longevity phenomenon is entirely controlled by the insulin signaling pathway (IIS-pathway). Both vertebrates and invertebrates have IIS-pathways that are comparable to one another, though no one has previously described de novo transcriptome assembly of IIS-pathway-associated genes in termites. In this research, we analyzed the transcriptomes of both reproductive (primary kings “PK” and queens “PQ”, secondary worker reproductive kings “SWRK” and queens “SWRQ”) and non-reproductive (male “WM” and female “WF” workers) castes of the subterranean termite Reticulitermes chinensis. The goal was to identify the genes responsible for longevity in the reproductive and non-reproductive castes. Through transcriptome analysis, we annotated 103,589,264 sequence reads and 184,436 (7G) unigenes were assembled, GC performance was measured at 43.02%, and 64,046 sequences were reported as CDs sequences. Of which 35 IIS-pathway-associated genes were identified, among 35 genes, we focused on the phosphoinositide-dependent kinase-1 (Pdk1), protein kinase B2 (akt2-a), tuberous sclerosis-2 (Tsc2), mammalian target of rapamycin (mTOR), eukaryotic translation initiation factor 4E (EIF4E) and ribosomal protein S6 (RPS6) genes. Previously these genes (Pdk1, akt2-a, mTOR, EIF4E, and RPS6) were investigated in various organisms, that regulate physiological effects, growth factors, protein translation, cell survival, proliferation, protein synthesis, cell metabolism and survival, autophagy, fecundity rate, egg size, and follicle number, although the critical reason for longevity is still unclear in the termite castes. However, based on transcriptome profiling, the IIS-pathway-associated genes could prolong the reproductive caste lifespan and health span. Therefore, the transcriptomic shreds of evidence related to IIS-pathway genes provide new insights into the maintenance and relationships between biomolecular homeostasis and remarkable longevity. Finally, we propose a strategy for future research to decrypt the hidden costs associated with termite aging in reproductive and non-reproductive castes.

Insecticidal and P450 mediate metabolism of fluralaner against red imported fire ant, Solenopsis invicta (Hymenoptera: Formicidae)

The red imported fire ant (Solenopsis invicta), a worldwide invasive and polyphagous pest, and often nests in residential areas. Finding an alternative pesticide that is both effective on S. invicta and environmentally friendly is urgent and crucial. Fluralaner, a novel isoxazoline insecticide, has been proven to possess selective toxicity for insects versus mammals and has been safe for mammals and non-target organisms, suggesting its potential in pest management. However, little toxicity information is available for the controlment of S. invicta. In this article, we studied the toxicity of fluralaner against S. invicta systematically, and the roles of metabolism-related enzymes in the metabolism process of fluralaner. The toxicity results showed that the topical application and feeding application were all effective for S. invicta. Moreover, fluralaner can be transmitted among workers by contacting and feeding which leads to a toxic reaction among nestmates. By exploring the biochemistry change, we found cytochrome P450 monooxygenase (P450) may be involved in the detoxification of fluralaner as well as carboxylesterase (CarE), but not glutathione S-transferase (GST). Synergism assays gave solid evidence in which piperonyl butoxide, an activity inhibitor of P450, increased the toxicity of fluralaner to S. invicta. Importantly, with the RNAi treatment, four of S.invicta P450 genes were significantly inhibited and showed more sensitivity to fluralaner at LC₅₀ concentration. Our result indicated that fluralaner could be a potential alternative pesticide in S. invicta control. And CYP9AS16, CYP6AS161, CYP6SQ20, and CYP336A45 genes were closely associated with the metabolism process of fluralaner.

Transcriptomics and metagenomics of common cutworm (Spodoptera litura) and fall armyworm (Spodoptera frugiperda) demonstrate differences in detoxification and development

BACKGROUND

Spodoptera litura is an important polyphagous pest that causes significant damage to the agricultural sector. We performed RNA-seq of 15 S. litura individuals from larval (fifth and sixth instar larvae), chrysalis, and adult developmental stages. We also compared the S. litura transcriptome data with Spodoptera frugiperda across the same developmental stages, which was sequenced in our previous study.

RESULTS

A total of 101,885 differentially expressed transcripts (DETs) were identified in S. litura. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses indicated that S. litura may undergo active xenobiotic and detoxifying metabolism during its larval and adult stages, which may explain difficulties with current population control measures. We also found that DETs of single-copy orthologous genes between S. litura and S. frugiperda were involved in basic metabolism and development. However, energy and metabolic processes genes had a higher expression in S. litura, whereas nervous and olfactory function genes had a higher expression in S. frugiperda. Metagenomics analysis in larval S. litura and S. frugiperda revealed that microbiota participate in the detoxification and metabolism processes, but the relative abundance of detoxification-related microbiota was more abundant in S. frugiperda. Transcriptome results also confirmed the detoxification-related pathway of S. frugiperda was more abundant than in S. litura.

CONCLUSIONS

Significant changes at transcriptional level were identified during the different development stages of S. litura. Importantly, we also identified detoxification associated genes and gut microbiota between S. litura and S. frugiperda at different developmental stages, which will be valuable in revealing possible mechanisms of detoxification and development in these two lepidopterans.

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.

Spodoptera frugiperda Sf9 cells as a model system to investigate the role of detoxification gene expression in response to xenobiotics

Spodoptera frugiperda (fall armyworm) is a highly destructive invasive pest that feeds on numerous crops including maize and rice. It has developed sophisticated mechanisms to detoxify xenobiotics such as secondary plant metabolites as well as manmade insecticides. The aim of the study was to explore the detoxification response to plant secondary metabolites and insecticides employing a S. frugiperda Sf9 cell model exposed to indole 3-carbinol (I3C) and methoprene. The cell Inhibitory Concentration 50 (IC₅₀) for these molecules was determined and IC₁₀, IC₂₀ and IC₃₀ doses were used to monitor the induction profiles of detoxification genes. Cytochrome P450 monooxygenases (P450s) of the CYP9A subfamily were the most inducible genes of the seven examined. Our results also showed the induction of the transcription factor Cap'n'collar isoform C (CncC). Transient transformation of Sf9 cells overexpressing CncC and its partner muscle aponeurosis fibromatosis (Maf) induces overexpression of CYP4M14, CYP4M15, CYP321A9 and GSTE1 while CYP9As were not induced. Next, we determined the capacity of recombinantly expressed CYP9A30, CYP9A31 and CYP9A32 to interact with methoprene and I3C. Fluorescence-based biochemical assays revealed an interaction of methoprene with functionally expressed CYP9A30, CYP9A31 and CYP9A32 whereas almost no interaction was detected for I3C, suggesting the ability of CYP9As to metabolize methoprene. Our results showed that Sf9 cells could be a useful model to decipher detoxification pathways of S. frugiperda.

Genomic and Transcriptomic Analysis Reveals Cuticular Protein Genes Responding to Different Insecticides in Fall Armyworm Spodoptera frugiperda

The fall armyworm (FAW), Spodoptera frugiperda, is a serious pest of crucial crops causing great threats to the food security of the world. It has evolved resistance to various insecticides, while the underlying molecular mechanisms remain largely unknown. Cuticular proteins (CPs), as primary components in cuticle, play an important role in insects' protection against environmental stresses. Few of them have been documented as participating in insecticide resistance in several insect species. In order to explore whether CP genes of the FAW exhibit a functional role in responding to insecticides stress, a total of 206 CPs, classified into eight families, were identified from the genome of the FAW through a homology-based approach coupled with manual efforts. The temporal expression profiles of all identified CP genes across developmental stages and their responses to 23 different insecticides were analyzed using the RNA-seq data. Expression profiling indicated that most of the CP genes displayed stage-specific expression patterns. It was found that the expression of 51 CP genes significantly changed after 48 h exposure to 17 different insecticides. The expression of eight CP genes responding to four insecticides were confirmed by RT-PCR analysis. The results showed that their overall expression profiles were consistent with RNA-seq analysis. The findings provide a basis for further functional investigation of CPs implied in insecticide stress in FAW.

Comparative analysis of the detoxification gene inventory of four major Spodoptera pest species in response to xenobiotics

The genus Spodoptera (Lepidoptera: Noctuidae) comprises some of the most polyphagous and destructive agricultural pests worldwide. The success of many species of this genus is due to their striking abilities to adapt to a broad range of host plants. Superfamilies of detoxification genes play a crucial role in the adaption to overcome plant defense mechanisms mediated by numerous secondary metabolites and toxins. Over the past decade, a substantial amount of expression data in Spodoptera larvae was produced for those genes in response to xenobiotics such as plant secondary metabolites, but also insecticide exposure. However, this information is scattered throughout the literature and in most cases does not allow to clearly identify candidate genes involved in host-plant adaptation and insecticide resistance. In the present review, we analyzed and compiled information on close to 600 pairs of inducers (xenobiotics) and induced genes from four main Spodoptera species: S. exigua, S. frugiperda, S. littoralis and S. litura. The cytochrome P450 monooxygenases (P450s; encoded by CYP genes) were the most upregulated detoxification genes across the literature for all four species. Most of the data was provided from studies on S. litura, followed by S. exigua, S. frugiperda and S. littoralis. We examined whether these detoxification genes were reported for larval survival under xenobiotic challenge in forward and reverse genetic studies. We further analyzed whether biochemical assays were carried out showing the ability of corresponding enzymes and transporters to breakdown and excrete xenobiotics, respectively. This revealed a clear disparity between species and the lack of genetic and biochemical information in S. frugiperda. Finally, we discussed the biological importance of detoxification genes for this genus and propose a workflow to study the involvement of these enzymes in an ecological and agricultural context.

Influence of three insecticides targeting GABA receptor on fall armyworm Spodoptera frugiperda: Analyses from individual, biochemical and molecular levels

The fall armyworm (FAW) Spodoptera frugiperda (Lepidoptera: Noctuidae) is a severe agricultural pest, which has invaded into China in 2019 and caused heavy damage to maize. The γ-aminobutyric acid receptor (GABAR)-targeted insecticides including broflanilide, fluralaner and fipronil exhibit high toxicity towards lepidopteran pests. However, whether they could be used for control of FAW and their possible mode of action in FAW remain unclear. In this study, broflanilide, fluralaner and fipronil exhibited high oral toxicity in FAW larvae with median lethal dose (LD₅₀) values of 0.677, 0.711, and 23.577 mg kg^-1 (active ingredient/ artificial food), respectively. In the electrophysiological assay, fluralaner and fipronil could strongly inhibit GABA-induced currents of homomeric FAW resistance to dieldrin 1 (RDL1) receptor with median inhibitory concentration (IC₅₀) values of 5.018 nM (95% confidence interval (CI) 2.864-8.789) and 8.595 nM (95% CI 5.105-14.47), respectively, whereas broflanilide could not. In addition, the cytochrome P450 (P450), glutathione-S-transferase (GST) and carboxylesterase (CarE) activities were positively response to broflanilide, P450 and GST to fluralaner, and GST and CarE to fipronil, respectively, compared with those of control. In conclusion, we firstly reported a notable insecticidal activity of three representative GABAR-targeted insecticides to FAW in vivo, and in vitro using electrophysiological assay. The GST is the primary detoxification enzyme for three tested insecticides. Our results would guide the rotational use of GABAR-targeted insecticides in field.

Identification and characterization of amphibian SLC26A5 using RNA-Seq

Background

Prestin (SLC26A5) is responsible for acute sensitivity and frequency selectivity in the vertebrate auditory system. Limited knowledge of prestin is from experiments using site-directed mutagenesis or domain-swapping techniques after the amino acid residues were identified by comparing the sequence of prestin to those of its paralogs and orthologs. Frog prestin is the only representative in amphibian lineage and the studies of it were quite rare with only one species identified.

Results

Here we report a new coding sequence of SLC26A5 for a frog species, Rana catesbeiana (the American bullfrog). In our study, the SLC26A5 gene of Rana has been mapped, sequenced and cloned successively using RNA-Seq. We measured the nonlinear capacitance (NLC) of prestin both in the hair cells of Rana’s inner ear and HEK293T cells transfected with this new coding gene. HEK293T cells expressing Rana prestin showed electrophysiological features similar to that of hair cells from its inner ear. Comparative studies of zebrafish, chick, Rana and an ancient frog species showed that chick and zebrafish prestin lacked NLC. Ancient frog’s prestin was functionally different from Rana.

Conclusions

We mapped and sequenced the SLC26A5 of the Rana catesbeiana from its inner ear cDNA using RNA-Seq. The Rana SLC26A5 cDNA was 2292 bp long, encoding a polypeptide of 763 amino acid residues, with 40% identity to mammals. This new coding gene could encode a functionally active protein conferring NLC to both frog HCs and the mammalian cell line. While comparing to its orthologs, the amphibian prestin has been evolutionarily changing its function and becomes more advanced than avian and teleost prestin.

Dynamic changes of transcriptome of fifth-instar spodoptera litura larvae in response to insecticide

Spodoptera litura is a major insect with a cosmopolitan distribution and strong resistance to multiple insecticides. Determining the molecular basis and key candidate genes of the insecticide resistance of S. litura may help in managing this insect. In this study, fifth-instar S. litura larvae were subjected to transcriptome analysis at 6, 12, 24, 48, and 72 h after feeding on an LC20 dose of avermectin. The result showed that genes responding to avermectin changed dynamically with different gene counts and resistance mechanisms at the fifth instar based on a metabolic pathway map. These responses included degrading the insecticide by a series of P450 and glutathione-S-transferase enzymes starting at the 12 h time point, with subsequent increases in the number of genes involved and shifts to TOLL and immune deficiency (IMD) pathways at 48 h after feeding the insecticide. Weighted correlation network analysis (WGCNA) determined a co-expression module related to the avermectin response at 12 and 24 h (r = 0.403, p = 0.0371; r = 0.436, p = 0.023), in which a hub gene (LOC111358940) related to metalloproteinase activity was identified. In addition, Analysis of the genes in the co-expression module further revealed that eight genes encoding UDP-glucuronosyltransferases were directly associated with insecticide response in S. litura. These results provide better understanding of the avermectin response mechanism of S. litura and may be useful in developing improved control strategies for this species.

SUPPLEMENTARY INFORMATION

The online version of this article (10.1007/s13205-021-02651-9) contains supplementary material, which is available to authorized users.

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.

Characterization of Rhopalosiphum padi takeout-like genes and their role in insecticide susceptibility

Due to the extensive use of chemical insecticides, the field populations of Rhopalosiphum padi, a serious wheat pest worldwide, have developed resistance to insecticides. Therefore, deep understanding of the mechanisms of the aphid's physiological response to insecticides would be of importance for the management of insecticide resistance in pests. Takeout belongs to a protein superfamily found exclusively in insects. Previous research showed that the takeout gene had various functions in insect physiology and behavior. However, few studies have explored the functions of takeout in insect insecticide susceptibility. The susceptibility of R. padi to imidacloprid and beta-cypermethrin was tested. Thirteen takeout-like genes were identified based on the genome database of R. padi. The number of exons was variable in these takeout-like genes, and nine highly conserved amino acids (two Cysteine, two Proline, four Glycine and one Aspartic acid) were identified. Expression levels of takeout-like-2, takeout-like-3, takeout-like-5, takeout-like-8, takeout-like-10 and takeout-like-11 were significantly increased after imidacloprid treatment; seven genes (takeout-like-1, takeout-like-2, takeout-like-5, takeout-like-6, takeout-like-7, takeout-like-8 and takeout-like-11) tended to be upregulated after beta-cypermethrin treatment. RNA interference results showed that the mortalities of R. padi injected with dsTOL-2, dsTOL-5, dsTOL-8, dsTOL-10 and dsTOL-11 were significantly increased after exposure to imidacloprid in comparison with control (injection of dsGFP). Under two sublethal concentrations of beta-cypermethrin, the silencing of takeout-like-2, takeout-like-5 and takeout-like-11 significantly increased the mortalities of R. padi. These results provide evidence for the involvement of takeout-like genes in insecticide susceptibility of R. padi, which improves our understanding the determinant of insecticide susceptibility.

Transcriptome analysis of Spodoptera litura reveals the molecular mechanism to pyrethroids resistance

Spodoptera litura is a destructive agricultural pest and has evolved resistance to multiple insecticides, especially pyrethroids. At present, the resistance mechanism to pyrethroids remains unclear. Four field-collected populations, namely CZ, LF, NJ and JD, were identified to have high resistance to pyrethroids comparing to pyrethroid-susceptible population (GX), with resistant ratio ranging from 11.5- to 9123.5-fold. To characterize pyrethroid resistance mechanism, the transcriptomes between two pyrethroid-resistant (LF and NJ) and a pyrethroid-susceptible (GX) populations were compared by RNA-sequencing. Results showed that multiple differentially expressed genes were enriched in metabolism-related GO terms and KEGG pathways. 35 up-regulated metabolism-related unigenes were selected to verify by qRT-PCR and 15 unigenes, including 4 cytochrome P450s (P450s), 5 glutathione S-transferase (GSTs), 1 UDP-glycosyltransferase (UGT), 4 carboxylesterases (COEs) and 1 and ATP-binding cassette transporters (ABC), were all up-regulated in the four pyrethroid-resistant populations. The expression levels of CYP3 and GST3, which were annotated as CYP6A13 and GSTE1, respectively, showed positive correlation with their pyrethroid resistance levels among the four pyrethroid-resistant populations. While the expression levels of CYP5, CYP12, COE4 and ABC5 showed good correlation with their pyrethroid resistance levels in at least three populations. UGT5 had the highest expression level among the tested UGT genes in the four pyrethroid-resistant populations. RNAi mediated silencing of CYP6 increased the cumulative mortality treated by beta cypermethrin and cyhalothrin significantly, while silencing of GST3 increased the cumulative mortality treated by fenvalerate significantly. CYP3, CYP5, CYP6, CYP12, GST3, COE4, UGT5 and ABC5 play important roles in pyrethroid resistance among the four pyrethroid-resistant populations. Our work provides a valuable clue for further study of pyrethroid resistance mechanisms in S. litura.