Transcriptome and Metabolome Analysis Reveals the Importance of Amino-Acid Metabolism in Spodoptera Frugiperda Exposed to Spinetoram

Simple Summary The spinetoram was high toxicity toward Spodoptera frugiperda larvae in, and was recommended to control Spodoptera frugiperda in China. In order to explore the effect of spinetoram on S. frugiperda, we studied the effect of spinetoram on amino acid metabolism of S. frugiperda by transcriptome analysis and LC-MS/MS analysis. The results showed this pesticide induced the change of arginine, glutamic acid, aspartic acid, and lysine cause in larvae of S. frugiperda, which was one of the factors of larval death. And, the down-regulation of phenylalanine may retard the tricarboxylic acid cycle to produce GTP. The decrease in citric acid is associated with an increase in leucine. However, enhancement of glucose metabolism and tryptophan provides the basis for the restoration of normal physiological functions of S. frugiperda larvae. Abstract Pests are inevitably exposed to sublethal and lethal doses in the agroecosystem following the application of pesticides indispensable to protect food sources. The effect of spinetoram on amino-acid metabolism of fall armyworm, Spodoptera frugiperda (J.E. Smith), was investigated, at the dose of LC10 and LC90, by transcriptome and LC-MS/MS analysis. Using statistics-based analysis of both POS and NEG mode, a total of 715,501 metabolites in S. frugiperda were significantly changed after spinetoram treatment. The enhancement of glucose metabolism provides energy support for detoxification in larvae. The decrease in valine and isoleucine is associated with an increase in leucine, without maintaining the conservation of citric acid in the larvae. The down-regulation of phenylalanine may retard the tricarboxylic acid cycle to produce GTP. The abundance of lysine was decreased in response to spinetoram exposure, which damages the nervous system of the larvae. The abundance of arginine increases and causes non-functional contraction of the insect’s muscles, causing the larva to expend excess energy. Tryptophan provides an important substrate for eliminating ROS. The changes in glutamic acid, aspartic acid, and lysine cause damage to the nerve centers of the larvae. The results of transcriptome and LC-MS/MS analysis revealed the effects of pesticide exposure on amino-acid metabolism of S. frugiperda successfully and provide a new overview of the response of insect physio-biochemistry against pesticides.


Introduction
The fall armyworm, Spodoptera frugiperda (J.E. Smith) (Lepidoptera: Noctuidae), is one of the migratory polyphagous pests of the neotropics of the Americas [1]. Due to the stronger migration ability of the adults, they rapidly spread from America to over 100 countries worldwide including China [2,3]. It has been reported to attack 353 host plant species from 76 plant families, and the economic damage from S. frugiperda can reach 73% worldwide from maize alone [4]. In China alone, the estimated damage area loss due to S. frugiperda was more than 134.77 million ha in 2020 [5]. Because chemical insecticides have rapid effects, especially during pest outbreak periods, the farmers often use massive amounts of chemical insecticides to control pests. However, the S. frugiperda has evolved resistance to various insecticides, and it has been confirmed that the S. frugiperda has potential resistance genes to carbamate and organophosphorus insecticides in many provinces of China [6,7]. Spinetoram is a new member of the spinosyn family, and it acts on the insect nervous system at the nicotinic acetylcholine receptor (nAChR) and gammaaminobutyric acid receptor (GABA) [8]. Due to its high toxicity toward S. frugiperda larvae in laboratory bioassays and field applications, spinetoram was recommended to control S. frugiperda by the Ministry of Agriculture and Rural Affairs of the People's Republic of China in 2020 [9].
Insecticides are mainly used by foliar spray and seed coating, so they exert toxic effects via ingestion. When insecticides are influenced by environmental degradation or via heterogeneous spatial coverage on individual plants, the exposed individuals are influenced by the sublethal dose of insecticides [10]. So, the study of sublethal effects is recommended as an emerging tool to assess the environmental risk of insecticides on arthropods. The sublethal dose of insecticides can negatively affect the larval duration and larval weight in Spodoptera exigua [11], pupal weight, and pupation rate in Plutella xylostella [12], adult longevity and reproduction in Spodoptera litura [13]. Several scientists have recently reported that the sublethal effects of pesticide have negative impacts on the visual learning and immunocompetence of honeybees [14] and the embryo development of Eriopis connexa [15]. The lethal dose of imidacloprid have negative impacts on aversive learning and movement of honeybees and that chronic exposure effects were dose-dependent for the insecticide [16].
With the development of third-generation sequencing technology, it gives us a chance to comprehensively understand gene expression after pesticides are used to treat the insect. The transcriptome can intuitively observe the changes in insects' energy metabolism under pesticide stress. Some studies have shown the down-regulation of the transcriptional levels of trehalase and trehalose phosphate synthase (tps) were observed in Lymantria dispar and Apis mellifera after sublethal exposure to different insecticides, respectively [17,18]. In addition, detoxification enzymes are known to be critical for metabolizing insecticides in invertebrates, and they can rapidly increase their activity in response to chemical stress [19]. Wang [20] observed that the expression of CYP321B1 (Cytochrome P450) in the midgut of S. litura could be significantly up-regulated by the sublethal concentration of betachlorpyrifos and chlorpyrifos. Moreover, Glutathione S-transferase (GST) is involved in dispelling harmful oxygen-free radicals caused by insecticide exposure [21,22], and GST has been found to be up-regulated in the resistance to phoxim and fenpropathrin in silkworm, Bombyx mori [23].
Metabolomics is an emerging technology that provides the real-time detection of the adaptive response of the organisms to environmental changes, pathophysiological stimuli, and genetic modifications [24,25]. So, metabolomics is a promising potential approach that could aid in understanding the stress effects of toxic pesticides on the exposed individuals. Imidacloprid induced high energy consumption, excitotoxicity, and oxidative stress (OS) in Aphis gossypii when exposed to a sublethal and lethal doses, but Parkinson's disease (PD) was detected only in the lethal dose [26]. Thiacloprid could induce a reduction in serotonin activity in the honeybee brain, which was implicated in learning and behavior development [27]. Additionally, metabolomics can provide host insect metabolic pathways response during pathogenic microorganism infection. It was recently reported that glycolysis and glutaminolysis are not essential during BmNPV infection in silkworms [28]. In general, metabolomics provides a reasonable explanation for the biological functions of individuals caused by the insecticides and pathogenic microorganism stress.
In previous studies, we examined the influence of spinetoram on the detoxifying enzymes and biological characteristics in S. frugiperda [9,29]. To observe how metabolism and its underlying transcriptional regulation are triggered during spinetoram stress, our analysis will reveal the importance of amino acid metabolism based on a metabolome and transcriptome dataset of S. frugiperda larvae taken at different doses of spinetoram stress. Herein, we discovered the connections among important metabolites in physiological regulation response. It will be extremely useful to understand potential metabolic molecular mechanisms under insecticide stress on insects.

Insect Rearing
S. frugiperda larvae were obtained from a laboratory colony from the Institute of Plant Protection of the Chinese Academy of Agricultural Sciences (IPP, CAAS).

Acute Toxicity Assay
Diet-overlay bioassay was conducted according to the method used by . Spinetoram was diluted with DMSO before the experiment. Formulated insecticides were diluted to generate five serial concentrations with distilled water containing 0.1% Triton X-100. Control groups were treated with 0.1% Triton X-100 only. Liquid artificial diet was dispensed into each well of a 24-well plate. After the diet cooled and solidified, 50 µL of insecticide solution was applied to the surface of each well. After the diet was dried at room temperature, one 3rd instar larva was placed in each well, and 72 larvae were tested for each concentration. Mortality was estimated after 2 days. If larvae died or could not move coordinately, they were considered dead. The LC 10 , LC 50 , and LC 90 values for spinetoram were calculated using the Polo Plus software (v1.0, LeOra Software, Parma, MO, USA).

Sample Preparation
Third instar larvae of S. frugiperda were treated with LC 10 and LC 90 dose of spinetoram for 48 h with the same method as described in 2.3, before the sample collection, while a control group was fed under the same conditions with 0.1% Triton X-100 solution. Each treatment and control group comprised of ten replications that included four replications that were taken for transcriptome, and others were used for the non-targeted metabolomics. It was ensured that the content of total tissue was not less than 50 mg. All samples were stored at −80 • C for further analyses. Metabolomic and RNA-seq analysis were performed by Biomarker Technologies Co., Ltd. (Beijing, China).

Transcriptome Analysis
The total RNA was extracted from each sample with TRIzol method. The quality of the total RNA from individual tissue samples was evaluated with electrophoresis in 1% agarose gels, and the RNA was quantified spectrophotometrically with a NanoDrop 2000 spectrophotometer (Thermo Scientific, Waltham, MA, USA) and an Agilent Bioanalyzer 2100 (Agilent, Santa Clara, CA, USA).
Based on sequence similarity, the final assembled unigenes with length of at least 200 bp were searched using BLASTx against the following databases: Gene Ontology (GO), euKaryotic Orthologous Groups (KOG), Clusters of Orthologous Groups (COG), Swiss-Prot, Pfam, Kyoto Encyclopedia of Genes and Genomes (KEGG), evolutionary genealogy of genes: Non-supervised Orthologous Groups (eggNOG), and NCBI non-redundant nucleotide collection (nr/nt). E-value threshold of 10 −5 was taken as the criterion of significant similarity.

Metabolite Extraction
An amount of 50 mg of insect sample was transferred to an Eppendorf tube. Then, the following steps were taken: Take 50 mg of one sample and homogenize it with 500 uL of ice-cold methanol/water (70%, v/v). Homogenate the mixture at 30 Hz for 2 min. After homogenization, shake the mixture for 5 min and incubate it on ice for 15 min. Centrifuge the mixture at 12,000 rpm at 4 • C for 10 min and suck supernatant 400 µL into another centrifuge tube. Add 500 µL of ethyl acetate/methanol (V, 1:3) into the original centrifuge tube, oscillate the mixture for 5 min, and incubate it on ice for 15 min. Then, centrifuge it at 12,000 rpm at 4 • C for 10 min and take 400 µL of supernatant. Merge the two supernatants and concentrate it. Then, add 100 µL of 70% methanol water into the dried product and perform ultrasonic treatment for 3 min. Finally, centrifuge it at 12,000 rpm at 4 • C for 3 min and control 60 µL of supernatant for LC-MS/MS analysis.

LC-MS/MS Analysis
The derivatized samples were analyzed on a Waters UPLC Xevo G2-XS QT of chromatography system coupled to an ultra-high-performance liquid chromatography-mass spectroscopy (UHPLC-MS, ExionLC, AB SCIEX; Waters, Manchester, UK). Acquity UPLC HSS T3 ((1.8 µm, 2.1 × 100 mm, Waters (Waters Technologies)) was utilized to separate the derivatives in electrospray ionization (ESI) positive (POS) and negative (NEG) ion mode. In each data acquisition cycle, dual-channel data acquisition for low and high collision energy can be carried out simultaneously. The minimal impact energy is 2 V, the high impact energy is 10-40 V, and the scanning frequency is 0.

Processing and Statistical Analysis of Metabolomics Data
The original data file obtained by LC-MS analysis is firstly converted into mzML format by ProteoWizard software. Peak extraction, alignment, and retention time correction are performed by XCMS program. The "SVR" method was used to correct the peak area. Peaks with a deletion rate > 50% in each group of samples were filtered. After that, metabolic identification information was obtained by searching the laboratory's self-built database and integrating the public database and mtDNA. Finally, statistical analysis was carried out by the R program. Statistical analysis includes univariate analysis and multivariate analysis. Univariate statistical analysis includes Student's t-test and variance multiple analysis. Multivariate statistical analysis includes principal component analysis (PCA), partial least squares discriminant analysis (PLS-DA), and orthogonal partial least squares discriminant analysis (OPLS-DA).

Toxicity of Spinetoram to S. Frugiperda 3rd Instar Larvae
The toxicity of spinetoram to 3rd larvae of S. frugiperda was determined by the method of diet-overlay bioassay. The LC 10 value was estimated as 0.172 mg/L (95% limit: 0.125 and 0.214 mg/L, respectively), and the LC 90 value was estimated as 0.725 mg/L (95% limit: 0.648 and 0.836 mg/L, respectively). The LC 10 and LC 90 values were used as the stress doses on S. frugiperda for transcriptome and non-targeted metabolomics experiments, respectively.

Transcriptome Analysis of S. Frugiperda in Responses to Spinetoram Application
For transcriptome profiling, three treatments: Control (CK), Spin-10 (LC 10 ), and Spin-90 (LC 90 ), in four replicates at 48 h post S. frugiperda stress were used to construct 12 cDNA libraries. A total of 83.29 Gb clean data was obtained and the Q30 base percentage was 93.03% and above. Clean reads were aligned with a designated reference genome (assembly ZJU_Sfru_1.0) and the comparison efficiency ranged from 82.35 to 84.17% (Table 1). To evaluate biological replicate correlation in control and different doses, we analyzed Pearson's correlation coefficient analysis among all treatments ( Figure 1A). The abundance of down-regulated genes was higher than up-regulated genes as compared to control at LC 10 dose ( Figure 1B), among them, 1627 unigenes were up-regulated and the remaining 1831 unigenes were down-regulated. However, the abundance of genes was reversed compared to control at LC 90 dose ( Figure 1B,C), among them, 1804 unigenes were up-regulated and the remaining 1580 unigenes were down-regulated. The color difference indicates high (red) and low (green) expressions.   To explore the response of different processes in control and toxic treatments, we analyzed the KEGG enrichment analysis of common differentially expressed genes (DEGs). The sets of DEGs were assigned to significant KEGG pathways (p > 0.05), which in Control-LC 10 , Control-LC 90, and LC 10 -LC 90 at 48 h, were listed in Figure 1D. The KEGG pathways were analyzed for genes enriched in a higher number of metabolic processes in all treatments ( Figure 1F) and genes enriched in a higher level of cytochrome P450 in LC 10 vs. LC 90 . Moreover, we observed a range of important metabolic pathway genes and detoxification enzyme genes after spinetoram stress, for example, 5-hydroxytryptamine receptor, acetylcholinesterase, catalase, connectin, superoxide dismutase, cytochrome P450, glutamate synthase, etc. (Table 2).
Both positive and negative ion modes of ESI were used during LC-MS/MS analysis. A comprehensive and extensive analysis of metabolic profiles from LC 10 , LC 90, and control groups were performed by multiple orthogonal methods. Firstly, the principal component analysis (PCA) was performed to analyze the metabolic profiles. The PCA showed a no-complete separation among the LC 10 , LC 90 , and control groups at pos modes ( Figure 2A). The principal component analysis (PCA) model was performed to analyze the metabolism profiles. The PCA model showed a rough separation among the LC 10 , LC 90 , and control groups at the pro-model. The PCA models obtained from LC 10 -control and LC 90 -control showed a complete separation, indicating a significant difference in metabolic characteristics. However, no-complete separations were obtained from LC 10 -control and LC 90 -control at neg-modes ( Figure 2B). The PCA model showed a rough separation among the LC10, LC90, and control groups at the pro-model. The PCA models obtained from LC10-control and LC90-control showed a complete separation, indicating a significant difference in metabolic characteristics. However, no-complete separations were obtained from LC10-control and LC90-control at neg-modes ( Figure 2B). Using statistics based on metabolite identity determined by analysis of both POS and NEG mode, a total of 715, 501 metabolites in S. frugiperda were significantly changed after spinetoram infection, respectively. A total of 66 metabolites were significantly up-regulated, and 225 metabolites were significantly down-regulated at LC10-pos and LC90-pos ( Figure 3A). In the other model, a total of 114 and 161 metabolites were significantly upregulated at LC10-neg and LC90-neg while 94 and 132 metabolites were significantly downregulated at these time points ( Figure 3A). To further observe patterns of overall metabolite abundance, the heat map analysis was performed, which displayed a significant difference in the abundance of metabolites between control group, LC10, and LC90. In addition, we have further identified metabolites that associated with spinetoram treatment ( Figure 3B,C). After 48 h of stress, most metabolites were significantly down-regulated compared to the control samples via pos-mode ( Figure 3A,B,C). However, most metabolites were significantly up-regulated compared to the control samples via neg-mode (Figure 3A,D,E).   Using statistics based on metabolite identity determined by analysis of both POS and NEG mode, a total of 715,501 metabolites in S. frugiperda were significantly changed after spinetoram infection, respectively. A total of 66 metabolites were significantly upregulated, and 225 metabolites were significantly down-regulated at LC 10 -pos and LC 90 -pos ( Figure 3A). In the other model, a total of 114 and 161 metabolites were significantly upregulated at LC 10 -neg and LC 90 -neg while 94 and 132 metabolites were significantly downregulated at these time points ( Figure 3A). To further observe patterns of overall metabolite abundance, the heat map analysis was performed, which displayed a significant difference in the abundance of metabolites between control group, LC 10 , and LC 90 . In addition, we have further identified metabolites that associated with spinetoram treatment (Figure 3B,C). After 48 h of stress, most metabolites were significantly down-regulated compared to the control samples via pos-mode ( Figure 3A-C). However, most metabolites were significantly up-regulated compared to the control samples via neg-mode ( Figure 3A,D,E). Metabolic pathway disruption was further investigated in S. frugiperda exposed to LC10 and LC90 of spinetoram using MetaboAnalyst, based on metabolites (Figure 4). The most significantly modified pathways between LC10 and the control group were the biosynthesis of amino acids, carbon metabolism, ubiquinone and other terpenoid-quinone Metabolic pathway disruption was further investigated in S. frugiperda exposed to LC 10 and LC 90 of spinetoram using MetaboAnalyst, based on metabolites (Figure 4). The most significantly modified pathways between LC 10 and the control group were the biosynthesis of amino acids, carbon metabolism, ubiquinone and other terpenoid-quinone biosynthesis and glycine, serine, and threonine metabolism ( Figure 4A). The changes in cysteine and methionine metabolism, arachidonic acid metabolism, glycerophospholipid metabolism, linoleic acid metabolism, and terpenoid backbone biosynthesis were characterized in the LC 90 -control group at neg-mode ( Figure 4B). Furthermore, the most significant changes were observed in ABC transports, biosynthesis of amino acids, bile secretion, arginine and proline metabolism, phenylalanine metabolism, pantothenate, and CoA biosynthesis in pos-mode ( Figure 4C,D). In contrast with the observation that metabolites were mainly enriched in amino acid metabolism pathways, metabolites at neg-mode were enriched in amio sugar and nucleotide sugar metabolism between LC 10 and LC 90 groups ( Figure 4E,F). Moreover, the significant metabolic pathways were summarized to reveal the effects of the spinetoram in S. frugiperda, based on biomarkers (Table 3).

Discussion
Metabonomic has been playing an important role in understanding metabolic information [30]. In recent years, it has been widely used in the study of diverse areas, including diseases, plants, microbiology, nutrition, toxicology, and insects [31]. During our study, we focused on the toxicity of pesticides of fall armyworm larvae following spinetoram treatment and identified 715 and 501 metabolites with differential abundance at LC 10 and LC 90 in different models, respectively. Among all differential metabolites, including diverse amino acids, other compounds were confirmed and were found to be involved in carbon metabolism, serine and threonine metabolism, glycerophospholipid metabolism, linoleic acid metabolism, ABC transports, pantothenate, and CoA biosynthesis. Carbohydrate metabolisms were extremely varied in S. frugiperda exposed to a low dose. Glucose and lactic acid both play important roles in glycolysis and glucogenesis as part of the starch and sucrose metabolic pathways [32]. In our study, glucose was up-regulated in spinetoram studies relative to the control group, indicating an increase in energetics necessary for detoxification following pesticide exposure. Similar phenomena were observed in other insects in which supply of glucose was up-regulated after exposure to neonicotinoid to supply energy for detoxification of pesticides [27].
Amino acids are biologically important molecules that could also serve as an energy source [33]. In the treatment group, the metabolism of valine and isoleucine were down-regulated to reduce the accumulation of acetyl-CoA. The metabolism of leucine was up-regulated to make up for deficiency in the accumulation of acetyl-CoA and maintain the conservation of citric acid. Phenylalanine metabolism could produce tyrosine, and fumaric acid was the end product of tyrosine [34,35]. So, the down-regulation of phenylalanine may reduce the amount of fumaric acid, and then retard the tricarboxylic (TCA) cycle to produce GTP under spinetoram stress. In addition, phenylalanine are essential amino acids for mammals; they can be synthesized by the gut microbiota [36][37][38]. The LC 90 treatment of spinetoram may affect the activity of advocate microorganisms and reduce the efficiency of phenylalanine synthesis, which may be the reason for the decrease in the abundance of phenylalanine under LC 90 treatment. In animals, lysine can also induce apoptosis in neuronal cells [39]. The abundance of lysine was decreased after spinetoram exposure, which damages the nervous system of the larvae. The arginine could be converted into phosphoarginine by using arginine kinase (AK) [40]. The abundance concentration of arginine was increased under spinetoram stress, which may positively impact energy utilization and muscle contraction. The arginine was increased causing non-functional contraction of the insect's muscles, stimulating the larva to spend excess energy. Additionally, pipecolic acid is a major metabolite of lysine and may facilitate GABAerguc transmissions by stimulating GABA release, inhibiting GABA uptake, or enhancing GABA A receptor responses [41][42][43][44]. The reduction of lysine in the LC 90 group indicated pipecolic acid was influenced, which is essential for the transmission of neural substrates in S. frugiperda larvae. Amino acid variation was one of the important factors that causes the death of S. frugiperda larvae exposed to spinetoram.
Glutamine is utilized in the cell as an essential source of carbon for the synthesis of macromolecules and the production of energy [45]. Glutamine was decreased in the LC 90 group, indicating that S. frugiperda reduced production of energy. Glutamine can be hydrolyzed by glutaminase to glutamate and ammonia and is an important source of glutamate. However, the decrease in glutamine did not affect the decrease in glutamate abundance in S. frugiperda after exposure to spinetoram. On the contrary, the abundance of glutamate was increased in the LC 90 group when exposed to spinetoram. Further, exposure to spinetoram also led to increased glutamate levels, which is also associated with increased intestinal energy requirements [46]. Glutamate serves as one of the major excitatory neurotransmitters in the central nervous system, which induces excitotoxicity in individuals when exposed to neurotoxicants and forms a glutamate-gated chloride channel (GluCls) with ionic inhibitory glutamate receptor (IGuRs) [47,48]. The nAchRs, by binding with nicotine, can stimulate the release of glutamate [49]; this may be a cause of the abundance of glutamate being increased after exposure to spinetoram. Aspartic acid also performed as an excitatory neurotransmitter in the central nervous system [50], so the abundant concentration of aspartic acid could lead to excitotoxicity in the spinetoram treatment. Therefore, the glutamine, glutamate, and aspartic metabolites could be potentially used as biomarkers of excitotoxicity in S. frugiperda in response to the stress of spinetoram. Excitotoxicity could be one reason for the high mortality of S. frugiperda in the spinetoram exposure.
Reactive oxygen species (ROS) could be formed after the application of the insecticide, which resulted in DNA damage and cell apoptosis in exposed individuals [51]. When the organism is in a certain pathological state (e.g., nicotinic toxins), the oxidative capacity in the body exceeds the antioxidant capacity and causes oxidative stress [52,53]. Some researchers indicated that the excessive production of ATP and the glutamate excitotoxicity were induced by bioenergetic disruption, which stimulates the formation of ROS [54,55]. So, increased levels of glutamate are another evidence that the larvae were being affected by an oxidative stress response. The high energy consumption and excitotoxicity elevated the ROS in S. frugiperda when exposed to LC 10 and LC 90 concentrations of the spinetoram. Most organisms have evolved antioxidant defense mechanisms to maintain cellular redox homeostasis (e.g., antioxidant enzymes: catalase (CAT), superoxide dismutase (SOD), glutathione reductase (GR), glutathione S-transferase (GST), and glutathion (GSH)), which are significantly affected by pesticides exposure [56,57]. Meanwhile, ascorbate (ASA) is the primary cellular antioxidant and is altered the transport and metabolic mechanisms responsible for dehydroascorbic acid (DHA) recycling by oxidative stress [58,59]. The ASA-GSH cycle, in which ASA participates, is considered to be an essential system for scavenging reactive oxygen species under environmental stress [60]. The up-regulation of ASA in the LC 90 group further indicates the increased ability in scavenging of ROS from oxidative pressure in S. frugiperda. Vitamin E (alpha-, beta-, gamma-, and deltatocopherol) also serves as one of the major antioxidants in organisms, whereas antioxidant activity tocopherol was decreased by the addition of (−)-epicatechin (−)-epigallocatechin gallate [61]. We found that the down-regulation of delta-tocopherol was influenced by the up-regulation of epigallocatechin gallate in the LC 90 group. Arachidonic acid, per se, can help alleviate oxidative stress [62], and the up-regulation of arachidonic acid has a positive impact on scavenging ROS in spinetoram stress. Moreover, under the spinetoram, the increased arachidonic acid can induce inflammation, which may be a further reason for larvae death. Nicotinamide adenine dinucleotide (NAD) and its phosphorylated form, NADP+, are the major coenzymes of redox reactions in central metabolic pathways. Besides its crucial role in cellular metabolism, NAD also serves as the substrate of several families of regulatory proteins: protein deacetylases (sirtuins), ADP-ribosyltransferases (ARTs), and poly (ADP-ribose) polymerases (PARPs), which govern vital processes including gene expression, DNA repair, apoptosis, aging, cell cycle progression, and many others [63][64][65]. We found that the NAD and NADP+ were up-regulated at LC 90 group, compared to the control group. An increase in the levels of NAD and NADP+ confirms the induction of oxidative stress on exposure with these pesticides. In addition, it has been shown that NADH and NADPH are able to directly scavenge oxidants such as hypohalous acids produced during inflammation [66]; they play an active role in the larval removal of LC 90 spinetoram on themselves. However, they did not change at low dose concentration, and SOD may play a role in the larval antioxidant response ( Table 2).
Tryptophan is an essential amino acid for protein synthesis, and some important physiological metabolism, and participates in immune function [67]. 5-hydroxytryptamine (5-HT), which is synthesized by tryptophan, plays an important role in central neurotransmitters [68,69]. Degradation of tryptophan, resulting in lower blood levels and impaired cerebral production and release of serotonin, is enhanced by inter alia inflammation, pregnancy, and stress in all species investigated, including humans [70]. The spinetoraminduced ROS resulted in the decrease in tryptophan levels in larvae and the decrease in 5-HT levels. However, we found 5-HT were up-regulated, compared with the control group, in spinetoram stress and have a positive impact on scavenging ROS from oxidative pressure. We inferred that the decrease in tryptophan was attributed to increased NAD, NADP, and 5-HT synthesis in insects to counter oxidative stress in larvae. Immune responses and metabolic regulation are tightly coupled in all animals including insects, but the underlying mechanistic connections are far from clear [71,72]. 5-HT is mainly found in the mucous membranes of the digestive tract and is stored in small amounts in the central nervous system. It must be pointed out that only a minority of the metabolites in the hemolymph may have originated from hemocytes and that most metabolites that enter the circulatory system are produced by other tissues, mainly the fat body [28]. Pyrethroids inhibit KYNA synthesis and tryptophan metabolism, leading to de novo generation of nicotinamide adenine dinucleotide coenzyme (NAD+) [73], and may be one of the reasons for the increase in NAD under LC 90 treatment. Therefore, we should take into account the changes in metabolites affected by toxic substances or pathogenic microorganisms on different tissues or sites, which will be more conducive to our understanding of insect metabolic responses.
That nAChRs are widely expressed in the central nervous system and numerous functions of acetylcholine signals, have been elucidated, including the induction of dopamine (DA) [74,75]. Sporadic or idiopathic PD is a neurodegenerative disease, and the hallmark symptoms of PD are related to DA deficiency [76,77]. DA, one of the monoamines, is a neurotransmitter synthesized from tyrosine. The decrease in tyrosine was detected in the LC 90 group, indicating that the pathogenesis of the PD appeared in insects by exposure to a lethal dose of spinetoram, and it may further explain the high mortality in S. frugiperda after exposure to a lethal dose of spinetoram. Chlorfenapyr-treated zebrafish brain metabolome revealed a significant decrease in the metabolites such as tyrosine in high concentration [78] and is consistent with our experimental results. Additionally, sequential injections with paraquat resulted in a significant loss of dopaminergic neurons in the substantia nigra (SN) in a previous study [79]. Some scholars believed that the cause of PD is the loss of DA neurons and the subsequent imbalance of the DA/acetylcholine concentration ratio [80]. The pathogenesis of PD can depend on the self-regulation of the central nervous system to dynamically reshape the balance of DA/acetylcholine. So, we observed the decrease in acetylcholine in the LC 90 group compared with the LC 10 and control groups. We inferred that the appearance of this phenomenon also provides the conditions for the larvae to return to normal activities.

Conclusions
Physiological and biochemical changes in S. frugiperdai under the sublethal and lethal stress of spinetoram were analyzed based on transcriptome and LC-MS/MS analysis by identifying the regulated metabolites. Comprehensively, abnormal metabolism of amino acids causes energy consumption, excitotoxicity, ROS, and PD, which were the main cause of larval death under spinetoram exposure. The accumulations of lysine indicated the nonfunctional contraction of the insect's muscles and caused the larvae to spend excess energy. The accumulations of glutamate and asparagine caused the excitotoxicity in S. frugiperda in the LC 90 group. The increase in the biosynthesis of antioxidants demonstrated the response of S. frugiperda scavenged ROS under spinetoram exposure. The decrease in PD in the LC 90 group suggested that spinetoram exposure was an external cause for the occurrence of PD. In general, our findings provide new insights into the changes of amino acids in larvae under pesticide exposure.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to restriction e.g., privacy.