Comparative Proteomics and Metabonomics Analysis of Different Diapause Stages Revealed a New Regulation Mechanism of Diapause in Loxostege sticticalis (Lepidoptera: Pyralidae)

Histone acetylation is an important epigenetic mechanism that has been shown to play a role in diapause regulation. To explore the physiological and molecular mechanisms of histone deacetylase in the diapause process, LC-MS/MS analysis was used to perform TMT proteomic and metabolomic analysis on non-diapause (ND), pre-diapause (PreD), diapause (D), cold treatment (CT), and post-diapause (RD) stages of the meadow moth. A total of 5367 proteins were identified by proteomics, including 1179 differentially expressed proteins. We found 975 (602 up-regulated and 373 down-regulated), 997 (608 up-regulated and 389 down-regulated), 1119 (726 up-regulated and 393 down-regulated), 1179 (630 up-regulated and 549 down-regulated), 94 (51 up-regulated and 43 down-regulated), 111 (63 up-regulated and 48 down-regulated), 533 (243 up-regulated and 290 down-regulated), 58 (31 up-regulated and 27 down-regulated), and 516 (228 up-regulated and 288 down-regulated) proteins in ND and PreD, ND and D, ND and CT, ND and RD, PreD and D, PreD and CT, PreD and RD, D and CT, D and RD, and CT and RD stages, respectively. A total of 1255 differentially expressed metabolites were annotated by metabolomics. Through KEGG analysis and time series analysis of differentially expressed metabolites, we found that phospholipids were annotated in significantly different modules, demonstrating their important role in the diapause process of the meadow moth. Using phospholipids as an indicator for weighted gene co-expression network analysis, we analyzed the most relevant differentially expressed proteins in the module and found that ribosomal 40s and 60s subunits were the most relevant proteins for diapause. Because there have been studies that have shown that histone deacetylase is associated with the diapause of meadow moths, we believe that histone deacetylase regulates the 40s and 60s subunits of ribosomes, which in turn affects the diapause of meadow moths. This finding expands our understanding of the regulation of meadow moth diapause and provides new insights into its control mechanism.


Introduction
The meadow moth is a major migratory pest in China and even the whole northern hemisphere.Diapause is a defensive measure evolved by the meadow moth to adapt to the cold climate in northern China during its migration [1].It is of great significance to study the diapause mechanism of the meadow moth for scientific prevention and control of the occurrence and damage of meadow moth.Previous studies have shown that the entry of insects into diapause is not only the closure of gene expression [2], but also involves the complex regulation of metabolites and proteins in insect cells.The rapid and comprehensive progress of proteomics and metabonomics sequencing technology has been widely used to study various biological functions and processes of insects.In this study, we sequenced the non-diapause period, pre-diapause period, diapause period, cold treatment period, and post-diapause period of meadow moth by proteomics and metabonomics, in order to further understand the regulation mechanism of meadow moth diapause.
Diapause is an adaptive strategy employed by insects to enhance survival capabilities in response to adverse environmental conditions and seasonal changes.The meadow moth (Loxostege sticticalis) over-winters through diapause, with larvae entering diapause under short-day conditions in autumn.They survive the winter at temperatures as low as −20 • C and terminate diapause in the following spring, with the diapause period lasting up to six months.The meadow moth is a long-day developmental-type insect, where under long-day conditions (≥16 h), larvae almost completely avoid diapause; however, under short-day conditions of 15 h, some larvae enter diapause, and the number of diapausing larvae increases as the photoperiod shortens [3].Since the diapause stage of the meadow moth is spent encased in a cocoon and buried in the soil, it is widely believed that the photoperiod experienced during the late developmental stages is the primary light factor influencing diapause in this species [4].
Histone acetylation is a type of histone modification that can function independently of genomic changes at any stage of insect development, and may therefore play a role in regulating changes in gene expression necessary for diapause initiation, maintenance, and termination [5].In the cell, the helical chromatin structure is formed by packaging DNA into nucleosomes composed of histones H2a, H2b, H3, and H4.Histones H3 and H4 possess a tail that can be covalently modified to increase or decrease interactions between the DNA chain and nucleosomes, thereby altering transcriptional activity at specific DNA loci [6].Increasing the interaction between DNA and histones typically results in a compact chromatin structure, while reducing the interaction between histones and DNA creates regions of relaxed chromatin with helicase activity.Evidence that histone acetylation affects diapause transcription comes from studies on the southern ground cricket Allonemobius socius [7], the northern house mosquito Culex pipiens [8], and the flesh fly Sarcophaga bullata [9].Cui et al. conducted transcriptome analysis of diapause-regulating genes in the fall armyworm, finding that histone deacetylases play a significant role in regulating diapause [10].
Ribosomes, fundamental sites of protein synthesis in eukaryotes, comprise 60s and 40s subunits, which are assembled in the nucleus and independently transported out of the nucleus.The intricate biogenesis of ribosomes plays a critical role in regulating biological development at multiple levels, including translational regulation in multicellular organisms.Xu discovered that histone deacetylases in rice are involved in the deacetylation of ribosomal proteins, potentially affecting ribosomal stability and translational efficiency [11].Urwanisch found that histone deacetylases can down-regulate ribosomal proteins [12].
Previous studies have investigated the metabolomics of diapause in various insects: in diapausing wheat blossom midge (Sitodiplosis mosellana), significant increases were observed in trehalose and glycerol compared to non-diapausing individuals [13].In diapausing Asian tiger mosquito Aedes albopictus, 84.7% of diapause-specific metabolites were categorized as "lipids and lipid-like molecules" [14].These metabolites, showing significant differences between diapause and non-diapause stages in insects, likely play roles in energy provision and cold resistance during diapause.
This study utilized proteomics and metabolomics techniques to identify and investigate changes in protein and metabolite levels during different diapause stages of the oriental armyworm, providing foundational research data on proteins and metabolites associated with diapause.This aims to better our understanding of the molecular regulatory mechanisms underlying diapause formation in the oriental armyworm.

Original Data Analysis, Protein and Metabolite Identification of Meadow Moth with Different Diapause Treatments
In order to explore the key proteins in the diapause process of meadow moths, TMT quantitative proteomics and LC-MS technology were used to study and compare the differences in protein abundance in five stages of diapause of meadow moth larvae.The molecular weight of most identified proteins was less than 101 kda (Supplementary Figure S1).Most peptides are 6-44 amino acids in length (Supplementary Figure S2).Of the 13,111 proteins identified, 5367 were qualified, and the details are in the attached File S3.
In order to accurately measure the metabolic changes in meadow moth larvae at different diapause stages, six biologically replicated meadow moth larvae at different diapause stages were determined by LC-MS.A total of 2020 metabolites were identified and determined.Its details are shown in the attached File S4.
Proteins were quantified based on data correlation acquisition (DDA) proteomics.A total of 26,502 peptides were detected (Figure 1A).The 26,502 peptides were annotated in five databases (EggNOG, GO, KEGG, Pfam, SubCell-Location) [15].In total, 5367 proteins were annotated in the SubCell-Location database.The minimum number of genes annotated in the GO database was 3664 (Figure 1B).
roles in energy provision and cold resistance during diapause.
This study utilized proteomics and metabolomics techniques to identify and inv gate changes in protein and metabolite levels during different diapause stages of the ental armyworm, providing foundational research data on proteins and metabolites a ciated with diapause.This aims to better our understanding of the molecular regula mechanisms underlying diapause formation in the oriental armyworm.

Original Data Analysis, Protein and Metabolite Identification of Meadow Moth with Different Diapause Treatments
In order to explore the key proteins in the diapause process of meadow moths, T quantitative proteomics and LC-MS technology were used to study and compare the ferences in protein abundance in five stages of diapause of meadow moth larvae.The lecular weight of most identified proteins was less than 101 kda (Supplementary Fi S1).Most peptides are 6-44 amino acids in length (Supplementary Figure S2).Of 13,111 proteins identified, 5367 were qualified, and the details are in the attached file In order to accurately measure the metabolic changes in meadow moth larvae at ferent diapause stages, six biologically replicated meadow moth larvae at different pause stages were determined by LC-MS.A total of 2020 metabolites were identified determined.Its details are shown in the attached File S4.
Proteins were quantified based on data correlation acquisition (DDA) proteomic total of 26,502 peptides were detected (Figure 1A).The 26,502 peptides were annotate five databases (EggNOG, GO, KEGG, Pfam, SubCell-Location) [15].In total, 5367 prot were annotated in the SubCell-Location database.The minimum number of genes a tated in the GO database was 3664 (Figure 1B).

DEP Analysis of Meadow Moth in Different Diapause Treatments
Principal component analysis (PCA) of the omics dataset showed that PC 1 (45.7 and PC 2 (24.00%) were the two main components (Figure 2A), and biological replica had good clustering.In these five different diapause stages, the DEPs whose prote expression was up-regulated or down-regulated ≥ 2 times were identified, and the pr ome was compared every two stages to identify DEPs.A total of 1179 non-redun DEPs were identified (Figure 2B).

DEP Analysis of Meadow Moth in Different Diapause Treatments
Principal component analysis (PCA) of the omics dataset showed that PC 1 (45.70%) and PC 2 (24.00%) were the two main components (Figure 2A), and biological replication had good clustering.In these five different diapause stages, the DEPs whose proteome expression was up-regulated or down-regulated ≥ 2 times were identified, and the proteome was compared every two stages to identify DEPs.A total of 1179 non-redundant DEPs were identified (Figure 2B).The heat map and cluster analysis of DEPs in different comparison groups are shown in Figure 3, indicating that there are significant differences in the proteins of meadow moth larvae at different diapause stages.The heat map and cluster analysis of DEPs in different comparison groups are shown in Figure 3, indicating that there are significant differences in the proteins of meadow moth larvae at different diapause stages.

Functional Classification of DEPs in Different Diapause Stages
We conducted a Gene Ontology (GO) functional analysis of differentially expressed proteins (DEPs) across various comparison groups to identify enriched biological processes, cellular components, and molecular functions (Figure 4).Across the groups ND vs. PreD, ND vs. D, ND vs. CT, ND vs. RD, and PreD vs. CT, the top three categories of proteins annotated in biological process, cellular component, and molecular function were consistent.Specifically, the biological processes included the macromolecule metabolic process, the organonitrogen compound metabolic process, and the cellular nitrogen compound metabolic process.The cellular components comprised integral components of membranes, intracellular organelles, and membrane-bounded organelles.In molecular function, the categories were cation binding, nucleoside phosphate binding, and nucleotide binding.
In the ND vs. PreD group, the most highly enriched biological processes were macromolecule metabolic processes (15.49%), followed by organonitrogen compound metabolic processes (14.95%) and cellular nitrogen compound metabolic processes (11.01%).For molecular functions, the category with the highest enrichment was cation binding proteins (21.26%).In the PreD vs. D group, the most highly enriched biological processes were the organonitrogen compound metabolic process (17.13%), the macromolecule metabolic process (14.89%), and the cellular nitrogen compound metabolic process (10.01%).For molecular functions, the most enriched categories were cation binding (28.88%), followed by nucleic acid binding (13.16%) and then anion binding (12.97%).
Among the five samples, except for 1255 co-expressed metabolites, the remaining metabolites were not expressed in all samples; 240 metabolites were expressed separately in the ND sample, while 18 metabolites in the RD sample were not expressed in other samples.The remaining samples were expressed in two or more samples (Figure 7A).KEGG enrichment analysis was performed on 1255 co-expressed differential metabolites selected, and the most significant 20 metabolic pathways were identified, including amino acid metabolism, lipid metabolism, nucleotide metabolism, etc. (Figure 7B).Among the five samples, except for 1255 co-expressed metabolites, the remaining metabolites were not expressed in all samples; 240 metabolites were expressed separately in the ND sample, while 18 metabolites in the RD sample were not expressed in other samples.The remaining samples were expressed in two or more samples (Figure 7A).KEGG enrichment analysis was performed on 1255 co-expressed differential metabolites selected, and the most significant 20 metabolic pathways were identified, including amino acid metabolism, lipid metabolism, nucleotide metabolism, etc. (Figure 7B).Among the five samples, except for 1255 co-expressed metabolites, the rema metabolites were not expressed in all samples; 240 metabolites were expressed sepa in the ND sample, while 18 metabolites in the RD sample were not expressed in samples.The remaining samples were expressed in two or more samples (Figur KEGG enrichment analysis was performed on 1255 co-expressed differential metab selected, and the most significant 20 metabolic pathways were identified, including a acid metabolism, lipid metabolism, nucleotide metabolism, etc. (Figure 7B).Time series analysis can analyze the dynamic expression and clustering of genes [16].We performed basic clustering analysis on differentially accumulated metabolites (DAMs) using STEM software (Version 1.3.13;https://www.cs.cmu.edu/~jernst/stem/accessed on 18 May 2024) and selected clusters with higher correlation.A total of 1255 non-redundant differential accumulation metabolites were divided into 30 files according to their expression trends.Among the 30 files (0-29), 5 files showed significant trends of 29, 3, 27, 28, 1.These five files were divided into four clusters based on similar accumulation trends (Figure 7C).
We performed KEGG metabolite annotation analysis on each cluster to identify enriched metabolites.Figure 8A-E show the first to fifth clusters, and we found that phospholipids have been annotated in 145, including substances such as phospholipids and sphingolipids.This indicates that phospholipids play an important role in the diapause process of meadow moths.Time series analysis can analyze the dynamic expression and clustering of genes [16].We performed basic clustering analysis on differentially accumulated metabolites (DAMs) using STEM software (Version 1.3.13;https://www.cs.cmu.edu/~jernst/stem/accessed on 18 May 2024) and selected clusters with higher correlation.A total of 1255 nonredundant differential accumulation metabolites were divided into 30 files according to their expression trends.Among the 30 files (0-29), 5 files showed significant trends of 29, 3, 27, 28, 1.These five files were divided into four clusters based on similar accumulation trends (Figure 7C).
We performed KEGG metabolite annotation analysis on each cluster to identify enriched metabolites.Figure 8A-E show the first to fifth clusters, and we found that phospholipids have been annotated in 145, including substances such as phospholipids and sphingolipids.This indicates that phospholipids play an important role in the diapause process of meadow moths.

WGCNA Analysis of Proteomics of Meadow Moth in Different Diapause Treatment Groups
After data preprocessing of differentially expressed proteins of meadow moth in different diapause treatment groups, 4897 non-redundant proteins were obtained, and 9 WGCNA modules were identified.The correlation between the light blue module and JH was the highest (r = 0.694, p = 4.1 × 10 −3 ).Therefore, 2073 proteins in the light blue module play an important role in changing the JH content of the meadow moth, thus participating in the diapause of the meadow moth (Figure 9A,B).

WGCNA Analysis of Proteomics of Meadow Moth in Different Diapause Treatment Groups
After data preprocessing of differentially expressed proteins of meadow moth in different diapause treatment groups, 4897 non-redundant proteins were obtained, and 9 WGCNA modules were identified.The correlation between the light blue module and JH was the highest (r = 0.694, p = 4.1 × 10 −3 ).Therefore, 2073 proteins in the light blue module play an important role in changing the JH content of the meadow moth, thus participating in the diapause of the meadow moth (Figure 9A,B).
The results of KEGG enrichment analysis and GO analysis of 2073 proteins showed that most proteins were related to metabolism.In the secondary metabolic pathway, RNA metabolism-related proteins, ribosomal complex-related proteins, energy metabolismrelated proteins, and oxidative phosphorylation-related proteins account for the majority (Figure 9C,D).
Chilo suppressalis of Lepidoptera Pyralidae was selected as the reference species of the turquoise module group of meadow moth, and the top 10% connectivity proteins were used as the hub proteins in the important module to construct the interaction network.The results showed that the centrality of three proteins in the interaction network was large.Therefore, the nodes of these three proteins play the most important role in ensuring the tight connection of the whole network.The compact centrality values of these three nodes are also high (Figure 10A and Table S3).However, the three proteins with the highest centrality were annotated.In addition to A0A3S2LT35, the other two proteins with serial 2024, 29, 3472 9 of 17 numbers xp9313648.1 and xp9307912.1 in NCBI were annotated with 60s ribosomal subunit and 40s ribosomal subunit, respectively.Studies have shown that histone deacetylase, which regulates lysine acetylation in rice, can regulate ribosomal subunits, thereby affecting plant photosynthesis [17].Therefore, we believe that 40s and 60s ribosomal subunits are regulated by histone acetylase and deacetylase to compress chromatin structure, which makes it difficult for RNA polymerase to enter, thus reducing protein-related expression to affect diapause.The results of KEGG enrichment analysis and GO analysis of 2073 proteins showed that most proteins were related to metabolism.In the secondary metabolic pathway, RNA metabolism-related proteins, ribosomal complex-related proteins, energy metabolism-related proteins, and oxidative phosphorylation-related proteins account for the majority (Figure 9C,D).
Chilo suppressalis of Lepidoptera Pyralidae was selected as the reference species of the turquoise module group of meadow moth, and the top 10% connectivity proteins were used as the hub proteins in the important module to construct the interaction network.The results showed that the centrality of three proteins in the interaction network was large.Therefore, the nodes of these three proteins play the most important role in ensuring the tight connection of the whole network.The compact centrality values of these three nodes are also high (Figure 10A and Table S5).However, the three proteins with the highest centrality were annotated.In addition to A0A3S2LT35, the other two proteins with serial numbers xp9313648.1 and xp9307912.1 in NCBI were annotated with 60s ribosomal subunit and 40s ribosomal subunit, respectively.Studies have shown that histone deacetylase, which regulates lysine acetylation in rice, can regulate ribosomal subunits, thereby affecting plant photosynthesis [17].Therefore, we believe that 40s and 60s ribosomal subunits are regulated by histone acetylase and deacetylase to compress chromatin structure, which makes it difficult for RNA polymerase to enter, thus reducing protein-related expression to affect diapause.  .Analysis of diapause-related protein interaction network of meadow moth.(A) Proteinprotein interaction network.Note: nodes represent proteins and edges represent interactions between two proteins.The size of the node is directly proportional to the connectivity (degree) of the node; that is, the more edges connected to the node, the greater the connectivity (degree) and the larger the node, indicating that the node gene is more important in the network.(B) Protein network center coefficient distribution map.

GSEA of Proteomics
Since there are only changes in light conditions and temperature during the whole diapause process of meadow moths, we analyzed the proteins related to the light transduction pathway (map04744, map04745), the circadian rhythm pathway (map04710, 04711, 04713) and the thermogenesis pathway (map04714) annotated for meadow moth by GSEA.GSEA can reveal whether some genes are enriched in the concentration of cir- Analysis of diapause-related protein interaction network of meadow moth.(A) Protein-protein interaction network.Note: nodes represent proteins and edges represent interactions between two proteins.The size of the node is directly proportional to the connectivity (degree) of the node; that is, the more edges connected to the node, the greater the connectivity (degree) and the larger the node, indicating that the node gene is more important in the network.(B) Protein network center coefficient distribution map.

GSEA of Proteomics
Since there are only changes in light conditions and temperature during the whole diapause process of meadow moths, we analyzed the proteins related to the light transduction pathway (map04744, map04745), the circadian rhythm pathway (map04710, 04711, 04713) and the thermogenesis pathway (map04714) annotated for meadow moth by GSEA.GSEA can reveal whether some genes are enriched in the concentration of circadian gene-related protein/thermogenic gene-related protein, and whether the genes in this gene concentration are randomly distributed or enriched at the top or bottom of the corresponding gene concentration.The results showed that except the PreD/CT and D/CT groups, the other groups were not significant (Figure 11A).In contrast, in the PreD/CT and D/RD groups, the thermogenic pathway-related protein (map04714) showed a significant effect (Figure 11B).From ND to PreD/D, the temperature did not change, but the proteins related to thermogenic pathway had a significant effect on it.In the case of no change in environmental factors, D and CT groups were at the stage of diapause maintenance, and the proteins related to the thermogenic pathway had a significant impact on the difference between D and CT.In this study, RSEM software with default parameters was used for analysis.The five stages of meadow moth larvae were divided into 10 groups, and GSEA was carried out for each group pairing.

Combined Proteomics and Metabonomics Analysis
In the D/CT group, because it is in the diapause maintenance stage and has a sign cant impact on the differences in proteins related to the light transduction pathway a thermogenic pathway, the D/CT group was analyzed by proteomics and metabonom The common annotation of proteomics and metabolomics detected 170 pathways (Fig 12A and Table S6).Two-way orthogonal partial least squares combined with discrimin analysis (O2PLS) was used to integrate the metabonomic and proteomic analysis (Fig 12B).It showed that the most closely related metabolites and proteins included ribosom glutathione metabolism, lysosome, nucleotide metabolism, lipid metabolism, tryptoph

Combined Proteomics and Metabonomics Analysis
In the D/CT group, because it is in the diapause maintenance stage and has a significant impact on the differences in proteins related to the light transduction pathway and thermogenic pathway, the D/CT group was analyzed by proteomics and metabo-nomics.The common annotation of proteomics and metabolomics detected 170 pathways (Figure 12A and Table S4).Two-way orthogonal partial least squares combined with discriminant analysis (O2PLS) was used to integrate the metabonomic and proteomic analysis (Figure 12B).It showed that the most closely related metabolites and proteins included ribosome, glutathione metabolism, lysosome, nucleotide metabolism, lipid metabolism, tryptophan metabolism, protein digestion and absorption (Figure 12C).

Discussion
As an important agricultural pest in China, studying its diapause is of great significance for scientific prevention and control of meadow moth outbreaks.During diapause, various changes occur in neuroendocrine, molecular, cellular, enzymatic, metabolic, endocrine, and behavioral aspects [18].Diapause typically occurs at specific developmental stages of insects, namely embryo, larva/nymph, pupa, and adult [19].It is divided into three phases: pre-diapause, diapause, and post-diapause.The pre-diapause phase involves the induction and preparation for diapause, encompassing the initiation, maintenance, and termination stages, while the post-diapause phase is characterized by a state of quiescence and rapid development [20,21].Typically, since diapause occurs at life cycle stages that resist unfavorable developmental environments, the intensity of diapause varies significantly among different insects; the duration of diapause can serve as a measure of its intensity.During diapause, most insects do not feed or, in some cases of larvae and adults, feed minimally.This indicates that insects must have sufficient food reserves before diapause to meet their metabolic needs during diapause, and still have enough reserves to complete development and resume activities post-diapause [21,22].
In this study, a total of 5367 DEP and 2020 DAM species were identified, which may include the proteins and metabolites most closely related to diapause of meadow moth.These proteins and metabolites include not only DEP and DAM in the prophase of diapause, diapause, and diapause termination, but also DEP and DAM, which play a key role in the diapause process of meadow moth.We further analyzed the differentially expressed proteins in proteomics and screened 1179 differentially accumulated proteins at

Discussion
As an important agricultural pest in China, studying its diapause is of great significance for scientific prevention and control of meadow moth outbreaks.During diapause, various changes occur in neuroendocrine, molecular, cellular, enzymatic, metabolic, endocrine, and behavioral aspects [18].Diapause typically occurs at specific developmental stages of insects, namely embryo, larva/nymph, pupa, and adult [19].It is divided into three phases: pre-diapause, diapause, and post-diapause.The pre-diapause phase involves the induction and preparation for diapause, encompassing the initiation, maintenance, and termination stages, while the post-diapause phase is characterized by a state of quiescence and rapid development [20,21].Typically, since diapause occurs at life cycle stages that resist unfavorable developmental environments, the intensity of diapause varies significantly among different insects; the duration of diapause can serve as a measure of its intensity.During diapause, most insects do not feed or, in some cases of larvae and adults, feed minimally.This indicates that insects must have sufficient food reserves before diapause to meet their metabolic needs during diapause, and still have enough reserves to complete development and resume activities post-diapause [21,22].
In this study, a total of 5367 DEP and 2020 DAM species were identified, which may include the proteins and metabolites most closely related to diapause of meadow moth.
These proteins and metabolites include not only DEP and DAM in the prophase of diapause, diapause, and diapause termination, but also DEP and DAM, which play a key role in the diapause process of meadow moth.We further analyzed the differentially expressed proteins in proteomics and screened 1179 differentially accumulated proteins at different diapause stages.Through comparative analysis, cluster analysis, and weighted gene coexpression network analysis, the important proteins related to carbohydrate metabolism, stress and defense response, lipid metabolism, and insect development were screened out.Based on the metabonomic time series and KEGG enrichment analysis of five groups of diapause larvae of meadow moth, 1255 different metabolites were identified, including amino acid metabolism, lipid metabolism, nucleotide metabolism, etc.
KEGG metabolite analysis shows that among 1255 metabolites, phospholipids are closely related to diapause.Most of the secondary metabolic pathways are important regulators in the diapause process of meadow moth.Phospholipids can promote lipid metabolism within insects and are integral components of cell membranes, capable of binding with membrane proteins to regulate cellular metabolism, energy conversion, cell repair, and signal transduction [23].In this study, certain phospholipids such as PI (Phosphatidyl inositol), PC (Phosphatidyl choline, also known as lecithin), PE (Phosphatidyl ethanolamine, also known as cephalin), PS (Phosphatidyl serine), and PA (Phosphatidic acid) were significantly higher in the diapause-affected Ostrinia furnacalis compared to the non-diapause group.Research on Bombyx mori has also shown that phospholipid levels positively correlate with cold resistance [24].Increased phospholipid content may influence lipid metabolism, transformation, and intercellular energy transfer during diapause in Ostrinia furnacalis, thereby enhancing cold tolerance.Lysophospholipids, derivatives of phospholipid hydrolysis, are more polar than phospholipids and can enhance membrane fluidity upon extensive contact with cell membranes [25].
To investigate whether circadian rhythm gene-related proteins and thermogenesis gene-related proteins play a role in protein variations during different diapause stages of the fall armyworm, we conducted Gene Set Enrichment Analysis (GSEA) on various treatment groups.The results showed that the role of circadian rhythm gene-related proteins in diapause and diapause release of meadow moth was not obvious.Heat-generating genes are the main reason for the differences in each stage.This may be because the regulation network of circadian rhythm signals is extremely complex, and the analysis result is not significant [26].
WGCNA analysis of differentially expressed proteins revealed that 2073 proteins were highly correlated with JH (r = 0.694, p = 4.1 × 10 −3 ).We believe that the 2073 proteins most related to JH in the light blue module are the key protein set for regulating diapause of meadow moth.We analyzed the protein interaction network of the light blue module and found that the 40s and 60s subunits of the ribosome played an important role in the whole network.Previous studies have shown that HDAC may play an important role in the diapause of meadow moth [10].At the same time, it was reported that lysine acetylation may increase the combination diversity of ribosomal activity, and then affect the histone modification in plant cells [11].Therefore, we speculate that the 40s and 60s subunits of ribosomes are regulated by histone deacetylase to compress the chromatin structure during diapause of meadow moth, which makes it difficult for RNA polymerase to enter, thus reducing protein-related expression.That is to say, there may be an interaction between HDAC and 40s and 60s subunits of ribosomes.These results expand the relevant regulatory network of meadow moth diapause and provide new insights.
However, the regulation of thermogenic pathways directly involved in energy metabolism is more direct than that of circadian rhythm genes.In the next step, we will further verify the functions of ribosomal subunits and HDAC.At the same time, we studied the influence mechanism of histone acetylase-related genes on diapause in meadow moth, and then determined the molecular mechanism of acetylation and endocrine regulation coupling.In addition, in order to mediate acetylation regulation, we will study the cellular signal transduction and co-regulation of acetylation in insect diapause.Acetylation may be a 2024, 29, 3472 13 of 17 molecular marker for diapause regulation or developmental decision in insects.It is of great significance to analyze the specific role of epigenetic code in diapause regulation, and then to understand the diapause mechanism of insects from a broader perspective.

Insect Collection
The caterpillars used in this experiment were collected from the agricultural fields of Dalad Banner and Otog Front Banner in Ordos City, Inner Mongolia Autonomous Region.To ensure that the experimental insects maintained consistent vitality, we reared the collected over-wintering generation of caterpillars in an insect-controlled climate chamber at 21 • C, with a humidity of 60-70% and an L:D ratio of 16:8 for three generations, after which they underwent diapause treatment.

Acquisition of Diapause and Diapause-Lifting Larvae
The diapause of the meadow moth is induced under the conditions of 21 • C and a 12:12 L:D cycle.In this laboratory, the diapause induction training of the meadow moth was conducted according to previous methods [4].The 1st to 4th instar larvae of the meadow moth were reared under 21 • C and a 16:8 L:D cycle, while the 5th instar larvae were induced under 21 • C and a 12:12 L:D cycle (ND).The feeding conditions for the control group and experimental group were identical, but the control group maintained a 16L:8D light condition after entering the 5th instar.After the control group's meadow moths burrowed into the soil to make cocoons and emerged, we collected the non-burrowed portion of the diapause-treated group as pre-diapause (PreD), while collecting the cocoons from the control group, cutting off one end of its cocoon, and observing its growth.After the control group pupated for over a week, the diapause-treated group remained in a larval state, and was named the diapause group (D).The meadow moths after two weeks of diapause were transferred to diapauserelease conditions (0-5 • C, 0L:24D) for 20 days, serving as the cold treatment group (CT).After being placed under normal cultivation conditions (21 ± 1 • C, 16L:8D) for 7-10 days, some of the moths in the bottle began to pupate, and those that remained in larval form were considered diapause-released larvae (RD).All tested meadow moths were placed in a rearing box measuring 28 × 18 × 13.5 cm, with 50 individuals per box.Five states of meadow moths (ND, PreD, D, CT, RD) were subjected to proteomic sequencing (Figure 13A).Twenty samples were collected from each treatment and replicated three times, with each sample collected near the end of the photoperiod.

Insect Collection
The caterpillars used in this experiment were collected from the agricultural fields of Dalad Banner and Otog Front Banner in Ordos City, Inner Mongolia Autonomous Region.To ensure that the experimental insects maintained consistent vitality, we reared the collected over-wintering generation of caterpillars in an insect-controlled climate chamber at 21 °C, with a humidity of 60-70% and an L:D ratio of 16:8 for three generations, after which they underwent diapause treatment.

Acquisition of Diapause and Diapause-Lifting Larvae
The diapause of the meadow moth is induced under the conditions of 21 °C and a 12:12 L:D cycle.In this laboratory, the diapause induction training of the meadow moth was conducted according to previous methods [4].The 1st to 4th instar larvae of the meadow moth were reared under 21 °C and a 16:8 L:D cycle, while the 5th instar larvae were induced under 21 °C and a 12:12 L:D cycle (ND).The feeding conditions for the control group and experimental group were identical, but the control group maintained a 16L:8D light condition after entering the 5th instar.After the control group's meadow moths burrowed into the soil to make cocoons and emerged, we collected the non-burrowed portion of the diapause-treated group as pre-diapause (PreD), while collecting the cocoons from the control group, cutting off one end of its cocoon, and observing its growth.After the control group pupated for over a week, the diapause-treated group remained in a larval state, and was named the diapause group (D).The meadow moths after two weeks of diapause were transferred to diapause-release conditions (0-5 °C, 0L:24D) for 20 days, serving as the cold treatment group (CT).After being placed under normal cultivation conditions (21 ± 1 °C, 16L:8D) for 7-10 days, some of the moths in the bottle began to pupate, and those that remained in larval form were considered diapause-released larvae (RD).All tested meadow moths were placed in a rearing box measuring 28 × 18 × 13.5 cm, with 50 individuals per box.Five states of meadow moths (ND, PreD, D, CT, RD) were subjected to proteomic sequencing (Figure 13A).Twenty samples were collected from each treatment and replicated three times, with each sample collected near the end of the photoperiod.

Sample Preparation for Proteomics Analysis
An appropriate amount of sample in the frozen state is taken out and transferred to the vibrating tube.An appropriate amount of protein lysate (8 M urea + 1% SDS, containing protease inhibitor) is added.A high-flux tissue grinder was used to vibrate the sample 3 times, for 40 s each time.Pyrolysis was performed on ice for 30 min, and vortex mixing was carried out for 5-10 s every 5 min.The supernatant was collected after centrifugation at 4 • C for 30 min at 12,000 G. Subsequently, the protein content was determined by the BCA method, and the operation was carried out in strict accordance with the instructions of the BCA reagent [26].
The final concentration of teab was made to be 100 mm by adding teab (triethylammonium bicarbonate buffer) to 100 µg of protein; TCEP (tris (2-carboxyethyl) phosphine) was added so that the final concentration of TCEP was 10 mm, and the reaction was carried out at 37 • C for 60 min; IAM (iodoacetamide) was added to make the final concentration of IAM 40 mm, and the mixture was kept in the dark for 40 min at room temperature.Each tube was added with pre-cooled acetone (acetone/sample v/v = 6:1) and was precipitated at −20 • C for 4 h.Centrifugation was performed at 10,000× g for 20 min, and the sediment was collected.The sample was fully dissolved with 100 µL of 100 mm teab.Trypsin was added according to an enzyme/protein ratio (m/m) of 1:50, and enzymolysis was performed overnight at 37 • C. The TMT reagent was taken out at 20 • C and allowed to return to room temperature; acetonitrile was added, vortexed, and centrifuged.Each tube was added with 100 µg of polypeptide into a TMT reagent and incubated at room temperature for 2 h.Hydroxylamine was added, reacted at room temperature for 30 min, mixed with the same amount of labeled product in a tube, and drained with a vacuum concentrator [27,28].

Sample Preparation for Metabonomic Analysis
The solid sample was taken into a 2 mL centrifuge tube, and a grinding ball with a diameter of 6 mm was added; 400 µL of extract (methanol/water = 4:1 (v:v)) containing 0.02 mg/mL of internal standard (l-2-chlorophenylalanine) was used for the extraction of metabolites.The sample solution was ground for 6 min (−10 • C, 50 Hz) with a frozen tissue grinder, and then subjected to low-temperature ultrasound for 30 min (5 • C, 40 kHz).The sample was then placed at −20 • C for 30 min and centrifuged for 15 min (4 • C, 13,000× g), and the supernatant was transferred to the injection vial with an internal cannula for analysis [29].

LC-MS/MS Analysis
Proteomics in the second dimension was analyzed by liquid chromatography tandem mass spectrometry (evosep one combined with obitrap exploris 480 mass spectrometer) (Shanghai Meiji Biomedical Technology Co., Ltd, Shanghai, China).The peptide was dissolved in a mass spectrometry loading buffer and then loaded on a C18 column (150 mm).The volume flow rate was 300 NL/min; 1% formic acid aqueous solution was in phase a, with 100% acetonitrile plus 0. 1% formic acid in phase B, 0~2 min, 5~5% B; 2-30 min, 5-38% B; 30-40 min, 38-90% B; 40-44 min, 90-90% B. The MS scanning range (M/z) was 350-1500, the acquisition mode was DDA, and the fragmentation mode was HCD.The resolution of primary mass spectrometry was 60,000, the resolution of secondary mass spectrometry was 15,000, and the dynamic exclusion time was 30 s. Turbotmt intelligent acquisition improves the resolution of reported ion isotopes.
The instrument platform for metabonomic LC-MS analysis is the uhplc-q exactive system (Shanghai Meiji Biomedical Technology Co., Ltd.) of semefier.Chromatographic conditions: a 3 µL sample was separated by an HSS T3 chromatographic column (100 mm × 2.1 mm I.D., 1.8 µm) and detected by mass spectrometry.Mass spectrum conditions: the sample mass spectrum signal was collected in the positive and negative ion scanning mode, and the mass scanning range was m/z: 70-1050.The ion spray voltage and positive ion voltage were 3500 V, the negative ion voltage was 2800 V, the sheath gas was 40 psi, the auxiliary heating gas was 10 psi, and the ion source heating temperature was 400 • C, with 20-40-60 V cyclic collision energy, an MS1 resolution of 70,000, and an MS2 resolution of 17,500.

Bioinformatics Analysis
The original files of mass spectrometry off the machine were analyzed by proteomediscoverertm software 2.4, and the database used for searching the database was the Meiji biological self-developed database.The false discovery rate (FDR) of peptide identification was set to FDR ≤ 0.01, the mass tolerance of peptides was set to 20 ppm, and the fragment mass tolerance was set to 0.02 da, containing at least one specific peptide [30].
The original data of LC-MS were imported into the metabonomics processing software progenesis Qi (Waters Corporation, Milford, MA, USA) for baseline filtering, peak identification, integration, retention time correction, and peak alignment.Finally, a data matrix of retention time, mass charge ratio, and peak intensity was obtained.At the same time, MS and MSMS mass spectrometry information was linked to the metabolism public databases HMDB (http://www.hmdb.ca/,accessed on 12 December 2023) and metlin (https://metlin.scripps.edu/,accessed on 17 December 2023) and meggith built its own database to match and obtain the metabolite information.We uploaded the databasesearched data matrix to the meggitt cloud platform (cloud.majorbio.com,accessed on 17 December 2023) for analysis.Firstly, the data matrix was preprocessed as follows: the data matrix used an 80% rule to remove missing values, that is, to retain at least 80% of non-zero variables in a group of samples, and then fill the gap value (the minimum value in the original matrix fills the gap value).In order to reduce the error caused by sample preparation and instrument instability, the total normalization method was used to normalize the response intensity of the sample mass spectrum peak to obtain the normalized data matrix.At the same time, the variables with a relative standard deviation (RSD) > 30% of QC samples were deleted, and log10 logarithm processing was performed to obtain the data matrix for subsequent analysis.Secondly, the ropls package (version1.6.2) in R language was used to perform principal component analysis (PCA) and orthogonal least partial binary discriminant analysis (OPLS-DA) on the preprocessed data matrix, and the stability of the model was evaluated by seven cycles of interactive validation.The selection of significantly different metabolites was determined based on the variable weight value (VIP) obtained by the opls-da model and the p value of Student's t test.The metabolites with vip > 1 and p < 0.05 were significantly different metabolites.
Differential metabolites were annotated by the KEGG database (https://www.kegg.jp/kegg/pathway.html,accessed on 17 December 2023) to obtain the pathways involved by differential metabolites.Python software (Version 3.10) package scipy.statswas used for pathway enrichment analysis, and a Fisher exact test was used to obtain the biological pathway that was most relevant to the experimental treatment.

Conclusions
In summary, through TMT proteomic and metabolomic analysis on non-diapause (ND), pre-diapause (PreD), diapause (D), cold treatment (CT), and post-diapause (RD) stages of the meadow moth, we found significant differences in protein and metabolic profiles between diapause and non diapause states.Then, we screened for differential proteins and metabolites among the different treatment groups of grassland moths mentioned above, and annotated a total of 1179 differential proteins and 1255 co expressed differential metabolites.We conducted KEGG analysis and time series analysis on the co expressed metabolites, and found that phospholipids were annotated in significantly different modules, indicating their important role in the diapause process of grassland moths.Using phospholipids as a weighted gene co expression network analysis indicator to analyze differentially expressed proteins, it was found that ribosomal 40s and 60s subunits were the most relevant proteins for diapause.Based on previous research in our laboratory, it has been found that the histone deacetylase of the meadow moth is related to its diapause.The results demonstrate that histone deacetylase regulates the 40s and 60s subunits of ribosomes, thereby affecting the diapause of the meadow moth.

Figure 1 .
Figure 1.Proteomic sample information of meadow moth in different experimental groups.(A formation on the number of peptides and proteins obtained by proteome.(B) The number of teins annotated in the five databases.

Figure 1 .
Figure 1.Proteomic sample information of meadow moth in different experimental groups.(A) Information on the number of peptides and proteins obtained by proteome.(B) The number of proteins annotated in the five databases.

Molecules 2024 , 20 Figure 2 .
Figure 2. Processing and proteomic sample information of meadow moth in different experimental groups.(A) PCA principal component scores of the proteome of meadow moth at different diapause stages.(B) The number of genes annotated in the five databases.Venn diagrams show the unique and common DEPs in 10 different control groups, including DEPs of up-regulated and down regulated genes (ND vs. PreD, ND vs. D, ND vs. CT, ND vs. RD, PreD vs. D, PreD vs. CT, PreD vs. RD, D vs. CT, D vs. RD, CT vs. RD).The heat map and cluster analysis of DEPs in different comparison groups are shown in

Figure 2 .
Figure 2. Processing and proteomic sample information of meadow moth in different experimental groups.(A) PCA principal component scores of the proteome of meadow moth at different diapause stages.(B) The number of genes annotated in the five databases.

Figure 2 .
Figure 2. Processing and proteomic sample information of meadow moth in different experimental groups.(A) PCA principal component scores of the proteome of meadow moth at different diapause stages.(B) The number of genes annotated in the five databases.Venn diagrams show the unique and common DEPs in 10 different control groups, including DEPs of up-regulated and down regulated genes (ND vs. PreD, ND vs. D, ND vs. CT, ND vs. RD, PreD vs. D, PreD vs. CT, PreD vs. RD, D vs. CT, D vs. RD, CT vs. RD).The heat map and cluster analysis of DEPs in different comparison groups are shown in Figure3, indicating that there are significant differences in the proteins of meadow moth larvae at different diapause stages.

Figure 3 .
Figure 3. Venn and thermogram analysis of differentially expressed proteins in meadow moth larvae at different diapause stages.(A) The Venn diagram of unique and common DEPs in up-regulated proteins in 10 different comparison groups, as well as the heat maps and cluster analysis of DEPs in different comparison groups.(B) The Venn diagram of unique and common DEPs in downregulated proteins in 10 different comparison groups, as well as the heat maps and cluster analysis

Figure 3 .
Figure 3. Venn and thermogram analysis of differentially expressed proteins in meadow moth larvae at different diapause stages.(A) The Venn diagram of unique and common DEPs in up-regulated proteins in 10 different comparison groups, as well as the heat maps and cluster analysis of DEPs in different comparison groups.(B) The Venn diagram of unique and common DEPs in down-regulated proteins in 10 different comparison groups, as well as the heat maps and cluster analysis of DEPs in different comparison groups.The blue part in the picture represents upregulated proteins, while the red part represents downregulated proteins.

Figure 4 .
Figure 4.The number of differentially expressed proteins of meadow moth larvae at different diapause stages is related to the molecular function and cellular components during diapause.A-J are the GO annotation diagram of ND vs. PreD, ND vs. D, ND vs. CT, ND vs. RD, PreD vs. D, PreD vs. CT, PreD vs. RD, D vs. CT, D vs. RD, and CT vs. RD.

Figure 4 .
Figure 4.The number of differentially expressed proteins of meadow moth larvae at different diapause stages is related to the molecular function and cellular components during diapause.(A-J) are the GO annotation diagram of ND vs. PreD, ND vs. D, ND vs. CT, ND vs. RD, PreD vs. D, PreD vs. CT, PreD vs. RD, D vs. CT, D vs. RD, and CT vs. RD.

Figure 5 .
Figure 5. Metabonomic sample information of meadow moth in different diapause treatment groups.(A) PCA of different samples.The scores of the first two PCs and the explained variance.(B) Classification diagram of KEGG compounds in the metabolome.(C) HMDB compound classification diagram.

Figure 5 .
Figure 5. Metabonomic sample information of meadow moth in different diapause treatment groups.(A) PCA of different samples.The scores of the first two PCs and the explained variance.(B) Classification diagram of KEGG compounds in the metabolome.(C) HMDB compound classification diagram.

Figure 6 .
Figure 6.Thermogram of metabolites of meadow moth at different diapause stages.Each sample has 6 biological replicates.

Figure 6 .
Figure 6.Thermogram of metabolites of meadow moth at different diapause stages.Each sample has 6 biological replicates.

Figure 6 .
Figure 6.Thermogram of metabolites of meadow moth at different diapause stages.Each s has 6 biological replicates.

Figure 7 .
Figure 7. Screening of differentially expressed metabolites in different diapause treatment groups.(A) Venn diagram of screening different metabolites among different experimental groups of meadow moth.(B) KEGG enrichment analysis of significant correlation groups, * p < 0.05, ** p < 0.01, and *** p < 0.001.(C).Time series analysis of meadow moth.

Figure 7 .
Figure 7. Screening of differentially expressed metabolites in different diapause treatment groups.(A) Venn diagram of screening different metabolites among different experimental groups of meadow moth.(B) KEGG enrichment analysis of significant correlation groups, p < 0.001 summarized with 3 asterisks [* * *].(C) Time series analysis of meadow moth.

Figure 8 .
Figure 8. KEGG compound analysis was performed on the five significant modules of time series analysis, with (A-E) representing time series analysis clusters 1 to 5.

Figure 8 .
Figure 8. KEGG compound analysis was performed on the five significant modules of time series analysis, with (A-E) representing time series analysis clusters 1 to 5.

Molecules 2024 , 20 Figure 9 .
Figure 9. Screening and annotation of JH-related proteins.(A,B) WGCNA analysis of five experimental treatment groups.Note: genes are divided into modules according to the expression trend.One branch represents one gene, and each color represents one module.Red indicates high abundance, blue indicates low relative expression, and gray indicates genes that are not classified into a specific module.(C) GO enrichment analysis of related proteins.(D) KEGG enrichment analysis of related proteins.

Figure 9 .
Figure 9. Screening and annotation of JH-related proteins.(A,B) WGCNA analysis of five experimental treatment groups.Note: genes are divided into modules according to the expression trend.One branch represents one gene, and each color represents one module.Red indicates high abundance, blue indicates low relative expression, and gray indicates genes that are not classified into a specific module.(C) GO enrichment analysis of related proteins.(D) KEGG enrichment analysis of related proteins.Molecules 2024, 29, x FOR PEER REVIEW 11 of 20

Figure 10
Figure10.Analysis of diapause-related protein interaction network of meadow moth.(A) Proteinprotein interaction network.Note: nodes represent proteins and edges represent interactions between two proteins.The size of the node is directly proportional to the connectivity (degree) of the node; that is, the more edges connected to the node, the greater the connectivity (degree) and the larger the node, indicating that the node gene is more important in the network.(B) Protein network center coefficient distribution map.

Figure 10 .
Figure 10.Analysis of diapause-related protein interaction network of meadow moth.(A) Protein-protein interaction network.Note: nodes represent proteins and edges represent interactions between two proteins.The size of the node is directly proportional to the connectivity (degree) of the node; that is, the more edges connected to the node, the greater the connectivity (degree) and the larger the node, indicating that the node gene is more important in the network.(B) Protein network center coefficient distribution map.

Figure 11 .
Figure 11.GSEA of proteomics of meadow moth larvae with different diapause treatments.Light transduction-related proteins and circadian rhythm proteins play a role in different stage meadow moth.(B) Thermogenic proteins play a role in different stages of meadow moth.The curve represents the dynamic enrichment score (ES) value, and the highest point represents the value of the gene set.

Figure 11 .
Figure 11.GSEA of proteomics of meadow moth larvae with different diapause treatments.(A) Light transduction-related proteins and circadian rhythm proteins play a role in different stages of meadow moth.(B) Thermogenic proteins play a role in different stages of meadow moth.The top curve represents the dynamic enrichment score (ES) value, and the highest point represents the ES value of the gene set.

Molecules 2024 , 20 Figure 12 .
Figure 12.Proteomics and metabolite group joint analysis between D and CT.(A) Proteomic and metabonomic annotations of the metabolic pathways.(B) Partial least squares discriminant analysis (PLS-DA) score maps of different strains of samples.(C) KEGG enriched the top 10 most closely related metabolites and proteins.

Figure 12 .
Figure 12.Proteomics and metabolite group joint analysis between D and CT.(A) Proteomic and metabonomic annotations of the metabolic pathways.(B) Partial least squares discriminant analysis (PLS-DA) score maps of different strains of samples.(C) KEGG enriched the top 10 most closely related metabolites and proteins.

Figure 13 .
Figure 13.Experimental design of proteomics and metabonomics analysis of meadow moth at different diapause stages using tandem mass spectrometry standard TMT and LC−MS methods.(A) Sampling diagram of different experimental groups of meadow moth.Note: as shown in the figure, the lighting conditions are L:D = 12:12, and ND is selected for the first day.The stage of low-

Figure 13 .
Figure 13.Experimental design of proteomics and metabonomics analysis of meadow moth at different diapause stages using tandem mass spectrometry standard TMT and LC−MS methods.(A) Sampling

:
Pathways co annotated by proteomics and metabonomics.Author Contributions: L.S. and L.X. conceived and designed the experiment.L.S. conducted the experiment.F.Y. provided an experimental analysis platform.J.F., Q.S. and H.L. provided the test insect sources.L.S., F.Y. and Q.Z.analyzed the data.L.S. wrote the manuscript.All authors have read and agreed to the published version of the manuscript.Funding: This study was supported by the National Key Research and Development Program of China (No. 2022YFD1400604) and the Inner Mongolia Science and Technology Plan Project (No. 2021GG0396).Institutional Review Board Statement: Not applicable.Informed Consent Statement: Not applicable.