Transcriptome and metabolome comprehensive analysis reveal the molecular basis of slow-action and non-repellency of cycloxaprid against an eusocial pest, Solenopsis invicta

The eusocial pest, red imported fire ant (Solenopsis invicta), is a highly invasive species that poses significant threats to public safety, agriculture, and the ecological environment. Cycloxaprid, a newly identified effective, slow-acting, and non-repellent insecticide against S. invicta, allows contaminated individuals to transfer the insecticide among nestmates through body contact. However, the molecular-level changes occurring in S. invicta post cycloxaprid exposure and any molecular alterations contributing to the slow demise or decreased sensitivity remain unclear. In this study, transcriptomic and metabolomic techniques were used to investigate the molecular mechanisms of S. invicta exposed to cycloxaprid. Differential analysis results revealed 275, 323, and 536 differentially expressed genes at 12, 24, and 48 h, respectively. Genes involved in lipid and energy metabolism, DNA integration, and hormone synthesis were largely upregulated at 12 h, suggesting S. invicta might actively resist cycloxaprid impacts, and predominantly downregulated at 48 h, indicating further functional impairment and impending death. Also, we observed an imbalance in olfactory perception pathways at 12 h, which may indicate a disruption in the olfactory system of S. invicta. Metabolomic results showed that the regulation of most differential metabolites (DMs) was consistent with the expression changes of their related DEGs at different time points. Our study provides insights into the mechanism underlying slow-acting and non-repellent properties of cycloxaprid against S. invicta.


Introduction
The red imported fire ant, Solenopsis invicta, is a highly destructive eusocial pest that has invaded more than 20 countries and regions worldwide, including the United States, Australia, and China (Ascunce et al., 2011;Wetterer, 2013).This pest poses significant threats to public safety, agriculture, and the ecological environment (Gutrich et al., 2007;Vinson, 2013;Wang et al., 2020).Medically, S. invicta stings can cause intense pain, blistering, pustules, allergic reactions, and, in extreme cases, anaphylactic shock or fatal outcomes (Haddad Junior and Larsson, 2015;Bernaba et al., 2019).Also, S. invicta can cause damage to the roots and seeds of various crops, leading to reduced yields and economic losses in agricultural systems (Chan and Guénard, 2020).Additionally, the invasion of S. invicta significantly deceased the abundance and diversity of native ants and other arthropods (Morrison, 2002).Given these issues, S. invicta resulted in massive global economic costs for damage and management (Angulo et al., 2022).For example, S. invicta was estimated to lead to annual costs of USD 6 billion in the United States and AUD 1.65 billion in Australia (Pereira et al., 2002;Wylie and Janssen-May, 2017).
Eusocial insects exhibit complex behaviors and biological features (e.g., highly organized social structures, hierarchical communication systems, task specialization, and unique foodsharing behaviors) that differentiate them from solitary insects (Robinson, 1992;Gordon, 2019).These characteristics not only ensure the high adaptability and flexibility of eusocial pests to invade diverse habitats but also pose challenges to their control.Traditional insecticides can rapidly eliminate individuals of eusocial pests outside the nest, but they cannot easily reach the reproducers and larvae concealed deep within the nest, resulting in inadequate control at a colony level (Drees and Gold, 2003;Silverman and Brightwell, 2008;Rust and Su, 2012).Therefore, insecticides to control eusocial pests such as S. invicta should be slow-acting and non-repellent to ensure the horizontal transfer of insecticides among nestmates through physical contact and eventually eradicate the whole colony.Currently, insecticides such as imidacloprid, fipronil, rotenone, and spinosad have been widely applied for S. invicta management (Chen and Oi, 2020).
Cycloxaprid is a novel oxabridged cis-nitromethylene neonicotinoid insecticide characterized by a nitro (-NO 2 ) group in a cis-configuration (Shao et al., 2013a).This insecticide appears to activate a different site compared to traditional neonicotinoid insecticide on the nicotinic acetylcholine receptors (nAChRs) and exhibits high insecticidal activity (Shao et al., 2013b;Zhang et al., 2019).For now, it has been applied to control various pests, including A. gossypii, Bemisia tabaci, and Sitobion avenae (Cui et al., 2016;Wang et al., 2018a;Cui et al., 2020).For example, Cui et al. (2012) reported that cycloxaprid exhibited intense contact and root-systemic activity against S. avenae, negatively affecting its feeding behavior and growth rate.Also, exposure to sublethal levels of cycloxaprid significantly shortened the pre-oviposition period and life expectancy of Aphis gossypii and inhibited the activity of various detoxification enzymes (Cui et al., 2018).
Our previous studies showed that cycloxaprid is a slow-acting and non-repellent insecticide against S. invicta workers at certain concentrations (Zhang et al., 2022).For example, surfaces treated with cycloxaprid (0.3 μg/cm 2 ) led to high mortality in S. invicta workers.However, cycloxaprid caused a much slower speed of kill than other fire ant control agents such as bifenthrin.In addition, cycloxaprid did not significantly influence the foraging, digging, and nesting behaviors of S. invicta workers.Interestingly, both corpses and live ants previously exposed to cycloxaprid caused horizontal toxicity against untreated ants.However, the molecular-level changes within S. invicta post cycloxaprid exposure remain unclear.In this study, we analyzed changes in gene expression and metabolites in S. invicta workers after exposure to cycloxaprid for 12, 24, and 48 h.By examining the differences in genes and metabolites, we aimed to reveal the mechanism of action of cycloxaprid on S. invicta and discuss the potential reasons underlying its slow-acting and non-repellent properties.
2 Materials and methods

Ant collection and maintenance
Six S. invicta colonies were collected from Tianlu Lake Forest Park (23 °12′ N, 113 °8′ E), Guangzhou, China, on 3 November 2021.Mounds were more than 5 m apart at the collection site.Nest materials with immatures (eggs, larvae, and pupae) and adults (workers and queens) of S. invicta were rapidly shoveled into a plastic container (length: 50 cm, width: 37 cm, height: 26 cm; Kangyi Co., Ltd., Hubei, China).The walls of the containers were previously smeared with talcum powder, and therefore S. invicta cannot climb and escape (Ning et al., 2019).The collected S. invicta colonies were brought to the laboratory and reared under relatively constant conditions (a photoperiod of 14 h light: 10 h dark, and room temperature of 21 °C-23 °C).Ten days before the experiment, S. invicta colonies were extracted from the soil using the waterdropping method provided by Chen (2007).The extracted ants were transferred to another plastic container (the inner walls were smeared with the talcum powder as mentioned above).Frozen crickets, water, and 20% honey water were supplied.

Cycloxaprid treatment
Cycloxaprid was obtained from the School of Pharmacy, East China University of Science and Technology, Shanghai, China.Plastic boxes (length: 44 cm, width: 29 cm, height: 19 cm; Kangyi Co., Ltd., Hubei, China) with talcum powder-coated inner walls were used as experimental arenas.Cycloxaprid was dissolved in acetone (analytical reagent, Fuyu Chemical Co., Ltd., Tianjin, China) to prepare for the aimed concentrations.The required amount of solution was then evenly smeared onto the bottom of the boxes to reach the cycloxaprid concentration of 0.3 μg/cm 2 after complete evaporation of the acetone.The same amount of acetone (without cycloxaprid) was added to the control boxes.For each colony, 10 g of S. invicta workers were weighed and transferred to a cycloxapridtreated or untreated box, and 10-mL tubes containing water or 20% honey water were supplied.Here, we only tested worker ants because the immatures and reproductive castes (i.e., queens) may have different molecular responses that may interfere with the results of workers (Kohlmeier et al., 2017;Negroni et al., 2021).Each box was placed in a larger plastic container (length: 50 cm, width: 37 cm, height: 26 cm), where moistened paper towels were placed on the bottom to maintain humidity (Wen et al., 2022).The bioassays were maintained at room temperature (21 °C-23 °C).At 12, 24, and 48 h, 500 mg of live S. invicta workers were randomly collected from each experimental box and immediately transferred to liquid nitrogen, then stored in a −80 °C freezer.

Transcriptomic analysis
The transcriptome of S. invicta workers collected from cycloxaprid-treated or untreated boxes was analyzed using RNAseq.Solenopsis invicta colonies 1-3 were used in this test, with each colony as an independent biological replicate.In the sample preparation process, 20 mg S. invicta workers were pulverized into a fine powder under liquid nitrogen.Subsequently, the total RNA was isolated using TRIzol reagent (Invitrogen, USA), adhering strictly to the manufacturer's guidelines.The RNA purification, library construction, and sequencing were carried out utilizing the Illumina sequencing platform, with all services provided by Metware Biotechnology Co., Ltd.(Wuhan, China).We employed Fastp (v0.19.3) to filter out low-quality reads and adapter sequences (Chen et al., 2018).We also utilized HISAT (v2.1.0)to construct the index and align clean reads with the reference genome (https://www.ncbi.nlm.nih.gov/assembly/GCF_016802725.1/#/st)(Kim et al., 2015).StringTie (v1.3.4d) was adopted to predict new genes (Pertea et al., 2016), and the transcripts were annotated using Diamond (v2.1.8)(Buchfink et al., 2015).The task of counting aligned reads and calculating the Fragments Per Kilobase of transcript per Million (FPKM) of each gene according to the gene length was performed using FeatureCounts (v1.6.2) (Liao et al., 2014).The differential expression between cycloxapridtreated and untreated ants was analyzed by utilizing DESeq2 (v1.22.1) (Love et al., 2014) and EdgeR (Robinson et al., 2010), with p values corrected through the Benjamini-Hochberg method (Haynes, 2013).To identify differentially expressed genes (DEGs), a threshold of |log2foldchange| (log2FC) > 1 and a p-value < 0.05 were applied as the markers of significant differential expression.
The functional analyses with gene annotation and enrichment analyses were conducted to investigate the functions of DEGs at different time points.DEGs were annotated using the Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Non-Redundant Protein Sequence Database (NR), Kentucky Online Gateway (KOG), Swiss-Prot, and Pfam databases.We executed enrichment analysis using the hypergeometric test.In the case of KEGG, the test was carried out at the pathway level, whereas for GO, the testing was based on the specific GO term.

Quantitative real-time PCR
To ascertain the precision and dependability of the transcriptome, we conducted an analysis using qRT-PCR on the total RNA of samples previously treated or untreated with cycloxaprid.Solenopsis invicta colonies 1-3 were tested, with each colony as an independent biological replicate.Eight selected genes (AASS, ACO, ORAI, SCD, CYP450-1, CYP450-2, CYP450-3, GST) were chosen based on their significant differential expression (see results) as well as their important roles in key biological processes (i.e., fatty acid synthesis and metabolism [AASS, ACO, SCD], energy conversion [ORAI], and detoxification [CYP450, GST]).We employed Ribosomal protein S3 (RPS3; XM_011172385) as the internal reference gene for normalization (Chen et al., 2022).Specific gene primer pairs were crafted utilizing the NCBI Primer-BLAST tool (http://www.ncbi.nlm.nih.gov/tools/primerblast/).The primers used in this study are detailed in Supplementary Table S1.Each sample's total RNA (2 µg) was transcribed in reverse to cDNA utilizing HiScript III RT SuperMix for qPCR (+gDNA wiper) (Vazyme, Nanjing, China), adhering strictly to the manufacturer's guidelines.Subsequently, the qRT-PCR process was executed employing the Option Real-Time PCR System CFX96 (Bio-Rad, CA).The procedure conditions comprised an initial pre-denaturation phase at 95 °C lasting for 2 min, followed by 39 cycles of denaturation at 95 °C for 5 s and annealing at 60 °C for 30 s.Following this, a melting curve was created by gradually escalating the temperature from 65 °C to 95 °C.Each qPCR operation was carried out in triplicate as three technical replicates to assess the variability introduced by the experimental procedure.The 2 −ΔΔCT method (Livak and Schmittgen, 2001) was employed to determine the relative expression levels of the target genes.Statistical scrutiny was undertaken using a paired t-test, with a p-value < 0.05 indicating the statistical significance.

Metabolomic analysis
The metabolome of S. invicta workers collected from cycloxaprid-treated or untreated boxes was analyzed using Widely-targeted metabolomic analysis (Chen et al., 2013).Six S. invicta colonies were tested, each as an independent biological replicate.Samples kept at −80 °C were defrosted on ice and then pulverized using liquid nitrogen.A solution of 400 μL, comprising methanol and water in a ratio of 7:3, was blended with an internal standard and added to the 20 mg of the crushed sample.This composite was then agitated at 1,500 rpm for 5 min.Following a 15-min icing period, the samples were centrifuged at 12,000 rpm for 10 min at 4 °C.A portion of the supernatant (300 μL) was set aside and cooled at −20 °C for 30 min, then centrifuged at 12,000 rpm for 3 min at 4 °C.For Liquid chromatography-mass spectrometry (LC-MS) analysis, a 200 μL aliquot of the supernatant was processed using the Ultra-Performance Liquid Chromatography and Tandem Mass Spectrometry (UPLC-MS/MS) systems provided by Wuhan Metware Biotechnology Co., Ltd.(China).
To identify Differential metabolites (DMs), a threshold of |log2foldchange| (log2FC) > 1 and a VIP >1 were applied as the markers of significant difference.VIP values were derived from Partial least squares discriminant analysis (PLS-DA) outcomes using the MetaboAnalystR package in R. Before PLS-DA, the data underwent log-transformation and mean-centering.A permutation test involving 200 permutations was implemented to prevent overfitting.The R function 'plsr' was employed to perform unsupervised PLS-DA, following a unit variance-scaling of the data.Annotation of DMs was accomplished through the KEGG Compound database, with mapping to the KEGG Pathway database.Statistically significant enrichment pathways were identified using the p-value from a hypergeometric test for a specified list of metabolites.

Correlation analysis of transcriptomic and metabolomic data
A correlation analysis was carried out between transcriptome and metabolome to gain deeper insight into the regulatory network of S. invicta workers responding to cycloxaprid.The transcriptomic and metabolomic analysis data from samples of colonies 1-3 were used for correlation analysis.The Pearson Correlation Coefficient (PCC) for differentially expressed genes and metabolites was computed employing the 'cor' function within the R package.To showcase the fold changes in genes and metabolites across different groups, a nine-quadrant plot was generated with a PCC >0.6 and p-value <0.05, using the ggplot2 and getopt packages in R. For generating a correlation network diagram, genes and metabolites with a PCC >0.85 and a p-value <0.05 were chosen, using Cytoscape (v3.9.1) as per Shannon et al. (2003).

Transcriptomic analysis
Each sample yielded 42,596,782 to 56,609,672 clean reads, with Q30 percentages ranging from 90.38% to 92.81% and GC content ranging from 44.29% to 48.15% (Supplementary Table S2), reflecting the high quality of the transcriptome sequencing data.Following the alignment of clean reads to the reference genome of S. invicta, we identified 15,673 to 15,967 expressed genes for each sample.A total of 275, 323, and 536 DEGs at 12, 24, and 48 h were identified, respectively, using a threshold of |log2FC| > 1 and p-value < 0.05 (Figure 1A).Venn diagram (Figure 1B) revealed 23 shared DEGs between 12 and 24 h, 24 shared DEGs between 12 and 48 h, 56 shared DEGs between 24 and 48 h, and 9 shared DEGs across all three time points.
Functional analyses provided insights into significantly enriched pathways (Figure 2; Figure 3) and the involvement of upregulated and downregulated DEGs in these pathways (Supplementary Figure S1, S2).At 12 h, the significantly enriched pathways primarily consisted of upregulated DEGs associated with fatty acid metabolism, fatty acid biosynthesis, AMPK signaling, insulin signaling, acyltransferase activity, DNA integration, and insect hormone biosynthesis.At 24 h, the significantly enriched pathways also primarily comprised upregulated DEGs, which were involved in protein digestion and absorption, pancreatic secretion, proteolysis, extracellular region, and amino sugar and nucleotide sugar metabolism.At 48 h, the significantly enriched pathways included both upregulated and downregulated DEGs.Specifically, DEGs involved in fatty acid metabolism, fatty acid biosynthesis, AMPK signaling pathway, insulin signaling pathway, iron ion binding, and heme binding were downregulated.In contrast, some DEGs involved with transmembrane transport, extracellular region, and cAMP signaling were upregulated.In addition, 22 DEGs involved in sensory perception of smell, odorant binding, and olfactory receptor activity displayed an inconsistent distribution at 12 h, with 13 upregulated and 9 downregulated.

Quantitative real-time PCR
Accordance was observed between the qRT-PCR expression patterns and the FPKM values valued from RNA-seq data (Figure 4).At specific time points, statistical analysis revealed significant differences (p-value < 0.05) in the expression levels of selected genes between cycloxaprid-treated and untreated ants.

Metabolomic analysis
UPLC-MS/MS detected 24 classes and 1,135 metabolites (Figure 5A).The PLS-DA analysis results demonstrated a progressive separation between the experimental and control groups across the three time points (12, 24, and 48 h).The separation gradually became more evident as the exposure time increased, and by 48 h, the clustering was distinctly clear (Figure 5B).A total of 24, 23, and 68 DMs at 12, 24, and 48 h were identified, respectively, using a threshold of | log2FC| > 2 and VIP >0.05 (Figure 5C).Venn diagram (Figure 5D) revealed 5 shared DEGs between 12 and 24 h, 2 shared DEGs between 12 and 48 h, 4 shared DEGs between 24 and 48 h, and 3 shared DMs across all three time points.
Utilizing the KEGG database, 346 out of 1,135 metabolites were annotated.Additionally, KEGG enrichment analysis revealed that DMs were enriched in 11, 13, and 15 pathways at 12, 24, and 48 h, respectively (Figure 6).Pathways, including metabolic pathways, neuroactive ligand-receptor interaction, biosynthesis of cofactors, and tyrosine metabolism, were significantly enriched at all three time points.Differential abundance score plots illustrated the changing trends of enriched pathways (Supplementary Figure S3-S5).At 12 h, DMs in the fructose and mannose metabolism and antifolate resistance were upregulated, whereas DMs in the caffeine metabolism were downregulated.At 24 h, DMs in the pertussis and Fc epsilon RI signaling pathway were upregulated, whereas DMs in the caffeine metabolism were downregulated.At 48 h, DMs in the pertussis and riboflavin metabolism were upregulated, whereas DMs in the sphingolipid metabolism and D-Glutamine and D-glutamate metabolism were downregulated.

Correlation analysis of the transcriptomic and metabolomic data
Correlation analysis revealed a significant association between 226 DEGs and 54 DMs (|R| > 0.8, p < 0.05).The nine-quadrant diagram (Figure 7) showed that the regulatory relationships between DEGs and DMs changed over time post cycloxaprid exposure.There was a notably higher number of points in Quadrants 1 and 3 at 48 h (Figure 7C) compared with 12 and 24 h (Figures 7A,B), indicating a significant shift in the gene expression patterns that caused an increased number of metabolites influenced by both upregulated and downregulated DEGs.Meanwhile, the number of points in Quadrant 7 remained relatively stable across the time points, suggesting a persistent association between downregulated DEGs and decreased DMs throughout the experiment.In addition, we generated a series of correlation network diagrams (|R| > 0.85, p < 0.05) to provide a more extensive perspective on the interaction between DEGs and DMs in S. invicta post cycloxaprid exposure (Figure 8), which highlights the intricate relationships between the transcriptomic and metabolomic changes in response to cycloxaprid exposure.

Discussion
A few recent studies focused on the molecular effect of insecticides on S. invica.For example, Siddiqui et al. (2022a) utilized transcriptomic analysis to investigate the detoxification mechanisms of S. invicta in response to beta-cypermethrin and fipronil, revealing significant alterations in metabolic pathways, insulin signaling pathways, and fatty acid metabolism.This study also detected changes in the expression of detoxification genes within the families of acetylcholinesterase (AChE), carboxylesterase (CarE), glutathione S-transferase (GST), and cytochrome P450 (CYP450).Similarly, Siddiqui et al. (2022b) reported that indoxacarb treatment also affected these pathways and genes in S. invicta.In addition, Xiong et al. (2022) utilized qRT-PCR and RNA interference (RNAi) technologies to investigate the function of specific CYP450 family genes of S. invicta treated with fluralaner, and they found that several genes, such as CYP9AS16, CYP6AS161, CYP6SQ20, and CYP336A45, were closely related to the detoxification process of fluralaner.These findings align with the present study, providing a consistent view of the insecticide response at the molecular level in S. invicta.
Cycloxaprid has recently garnered substantial attention as a promising alternative to traditional insecticides, owing to its potential to counter resistance issues for agricultural pests (Wang et al., 2018a;Cui et al., 2020).For example, Zhang et al. (2018) found a slower rate of resistance development to cycloxaprid in Laodelphax striatellus compared to chlorpyrifos and buprofezin.Remarkably, only a 10-fold increase in resistance was observed in an L. striatellus variant with cycloxaprid resistance, contrasting with the increases of 106-fold and 332-fold for buprofezin and chlorpyrifos, respectively, under equivalent selection pressures.Dong et al. (2022) revealed that an A. gossypii variant with cycloxaprid resistance demonstrated minimal cross-resistance to six of other insecticides.The fold increases in resistance were 4.3 for imidacloprid, 2.9 for acetamiprid, 3.7 for thiamethoxam, 6.1 for nitenpyram, 2.2 for flupyradifurone, and 4.5 for sulfoxaflor.Similarly, Cui et al. (2016) and Zhang et al. (2019) suggested no cross-resistance between cycloxaprid and imidacloprid in A. gossypii and Nilaparvata lugens.Some recent research has investigated the molecular of cycloxaprid on various sucking insect pests (Yang et al., 2016;Jin et al., 2017).For example, Jin et al. (2021) identified 22 genes (e.g., genes in the families of CYP450 and GST)

FIGURE 6
Comprehensive KEGG enrichment analysis of DMs at three distinct time points post cycloxaprid exposure.The size of each dot represents the number of genes involved in the corresponding KEGG pathways.p-values, calculated using hypergeometric tests, demonstrate the statistical significance of the enrichment results.

FIGURE 7
Nine-quadrant diagrams showed correlations between genes and metabolites in Solenopsis invicta at 12 (A), 24 (B), and 48 h (C) post cycloxaprid exposure.These diagrams display dynamic relationships between gene expression and metabolic changes in response to cycloxaprid exposure over time.
potentially associated with Sogatella furcifera's resistance to cycloxaprid using transcriptomic analysis.Zhang et al. (2018) found that the resistance to cycloxaprid in L. striatellus was linked to the modified expression of nicotinic acetylcholine receptor subunits.Moreover, they reported that a knockdown of the Lsa1 gene significantly amplified the susceptibility of the cycloxaprid-resistant strain to cycloxaprid.While these insights shed light on the resistance mechanisms of sucking insects, it is crucial to note that successful control of eusocial pests, such as S. invicta, primarily depends on the slow-action and non-repellency of cycloxaprid.
In our study, transcriptomic analyses revealed alterations in various physiological processes of S. invicta post cycloxaprid exposure over different durations.At 12 h, most DEGs were upregulated and enriched pathways involved in fatty acid metabolism, fatty acid biosynthesis, AMPK signaling pathway, insulin signaling pathway, and acyltransferase activity.Fatty acids are essential for energy storage and release, cell membrane assembly, and the production of signaling molecules in insects (Canavoso et al., 2001;Arrese and Soulages, 2010).The insulin signaling pathway assumes a crucial role in sustaining energy homeostasis through regulating carbohydrate, lipid, protein, and energy metabolism in insects (Wu and Brown, 2006;Géminard et al., 2009;Nässel and Broeck, 2016), which promotes glucose uptake and glycolysis, stimulates lipogenesis, and suppresses lipolysis and gluconeogenesis (Broughton et al., 2005).Acyltransferases are critical in synthesizing and degrading various lipids, including fatty acids, glycerolipids, and phospholipids (Van der Horst, 2003;Arrese and Soulages, 2010).In addition, the AMPK is activated under energy-deprived conditions (e.g., low ATP levels) and can initiate various downstream processes that help maintain energy balance by regulating fatty acid synthesis, oxidation, and other related processes (Pan and Hardie, 2002;Hardie, 2007;Sinnett and Brenman, 2016).Changes in these pathways collectively point to an enhancement in lipid and energy metabolism in S. invicta.Neonicotinoid insecticides can bind to nAChRs located in the central nervous systems, leading to excitation, convulsions, and death of insects, accompanied by increased energy consumption (Matsuda et al., 2001;Tomizawa and Casida, 2003;Jeschke et al., 2011).It is probable that S. invicta actively enhanced its lipid metabolism to generate more energy to compensate for the increased energy expenditure after cycloxaprid exposure, thereby maintaining survival and recovery.
At 24 h, most DEGs were upregulated and enriched in pathways involved with protein digestion and absorption, pancreatic secretion, proteolysis, and amino sugar and nucleotide sugar metabolism.The upregulation of genes involved in protein digestion and absorption may suggest increased efficiency in breaking down dietary proteins and assimilating amino acids (Holtof et al., 2019).This process could be essential for

FIGURE 8
The correlation network showed connections between genes and metabolites in Solenopsis invicta post cycloxaprid exposure.The network indicated the complex relationships and potential interactions between the transcriptomic and metabolomic changes post cycloxaprid exposure.
Frontiers in Physiology frontiersin.orgmaintaining energy balance because amino acids can be an energy source when required (MacFarlane, 2018).Similarly, the upregulation of genes involved in pancreatic secretion and proteolysis indicates the enhanced capacity to produce and secrete digestive enzymes, which are crucial for breaking down proteins into smaller peptides and amino acids (van Hoef et al., 2011).These results strongly correlate with findings at 12 h, indicating that S. invicta may enhance its protein and carbohydrate metabolism to obtain more energy to support the heightened energy consumption post cycloxaprid exposure.At 48 h, most DEGs were downregulated and enriched in pathways of fatty acid metabolism, fatty acid biosynthesis, AMPK signaling, and insulin signaling.This may indicate a significant reduction in lipid and energy metabolism of S. invicta at the later stage after cycloxaprid exposure.Furthermore, the downregulation of genes involved in iron ion and heme binding suggests a decreased ability to bind and utilize iron and heme.This could impact various biological processes such as oxygen transport, energy production, and DNA synthesis (Nichol et al., 2002;Walter-Nuno et al., 2013).Upregulated pathways mainly include oxidoreductase activity, transmembrane transport, extracellular region, and cAMP signaling, implying the increased oxidative stress response that may alter membrane transport, enhance extracellular communication, and modulate intracellular signaling cascades (Kunieda et al., 2006;Holtof et al., 2019;Wu et al., 2019;Liu et al., 2021).
Interestingly, 22 DEGs were enriched in pathways related to odorant perception at 12 h.Solenopsis invicta typically detects odor molecules through olfactory receptors located on their antennae.The alterations in these genes suggest a disturbance of olfactory perception in S. invicta, which may ultimately cause reduced sensitivity to cycloxaprid.Additional researches are required in the future to examine the potential involvement of these genes in the nonrepellent effect of cycloxaprid on S. invicta.We also investigated the expression of detoxification-related genes, and identified 22 CYP450 family genes and 2 GST family genes.At 12 h, LOC105207882 and LOC120359407 were upregulated, suggesting that ants may actively adjust their detoxification mechanisms to counteract the effects of cycloxaprid during early exposure.At 24 h, the detoxification response seems to weaken because only one gene (LOC105206318) was upregulated, whereas four genes (LOC105198272, LOC105203688, LOC120357488, and LOC120359317) were downregulated.This shift indicates that S. invicta may experience increased stress due to continued exposure to cycloxaprid, thereby hindering their ability to maintain effective detoxification processes.At 48 h, 18 detoxification DEGs were downregulated, and only two genes (LOC105205661 and LOC105206584) were upregulated, indicating that S. invicta may be unable to sustain its detoxification capacity after prolonged exposure, which may lead to impaired physiological function and eventually death.It is worth noting that the expression changes in the CYP450 family and GST family genes could have interconnected impacts on other physiological processes, such as lipid metabolism, energy metabolism, and immune function.For instance, at 12 h, the activation of these genes may necessitate energy, thus promoting lipid and energy metabolism.However, at 48 h, the decline in lipid and energy metabolism may contribute to the inability of these genes to maintain the function of detoxification.
Metabolomic analyses showed the diversity and functions of DMs in S. invicta following exposure to cycloxaprid over different durations.These DMs are various fatty acids, nucleotides, hormones, and amines, which participate in numerous physiological processes, including lipid and energy metabolism, DNA synthesis, hormone signaling, neurotransmission, and overall metabolic activity (Weyel and Wegener, 1996;Gäde et al., 2003;Kaczmarek and Boguś, 2021;Barbero et al., 2023).Correlation analysis showed that 226 DEGs exhibited significant correlations with 54 DMs (|R| > 0.8, p < 0.05), indicating the tight interplay between metabolites and genes in S. invicta treated with cycloxaprid.The correlation network (Figure 8) showed that the nicotinamide riboside (NR, MEDP1791) was strongly interconnected with numerous metabolites and genes.NR is a pyridine-nucleoside and a form of vitamin B 3 and serves as an essential precursor to nicotinamide adenine dinucleotide (NAD + ) (Bogan and Brenner, 2008), which is a cofactor for enzymes involved in energy metabolism and various metabolic pathways such as glycolysis and fatty acid oxidation (Xie et al., 2020).Therefore, NR might contribute to the adaptive responses of S. invicta to cycloxaprid, and future studies are needed to verify this.
Notably, we found that changes in fatty acids were inversely correlated with the expression levels of genes implicated in fatty acid synthesis and metabolism pathways.Transcriptomic data revealed upregulation of DEGs associated with fatty acid pathways at 12 h and downregulation at 48 h.In contrast, metabolomic analyses detected only two upregulated fatty acids at 12 h, and eight carnitines were still upregulated at 48 h.Carnitine is an endogenous compound carrier molecule in lipid metabolism and energy homeostasis (Bremer, 1983;Bieber, 1988;Pande and Murthy, 1989).The diminished lipid and energy metabolism in S. invicta may disrupt fatty acid cycling, resulting in carnitine accumulation at 48 h.Numerous studies involving humans and animals have demonstrated the beneficial effects of carnitine accumulation in treating diverse health conditions, such as arrhythmia, gastrointestinal discomfort, and epilepsy (Longo et al., 2006;Fukuda et al., 2015;Wang et al., 2018b).Concurrently, recent research suggests that the regulation of specific fatty acids might play a significant role in determining the potential of Leptinotarsa decemlineata to develop resistance to neonicotinoid insecticide (Clements et al., 2020).Therefore, the alterations in fatty acid (including carnitine) might also enhance the capacity of S. invicta to adapt to cycloxaprid.It is important to note that the present study was conducted using S. invicta workers, and future studies would be valuable to assess the impacts of cycloxaprid on different developmental stages of S. invicta.

Conclusion
In conclusion, our study provides a comprehensive understanding of the slow-acting and non-repellent properties of cycloxaprid against S. invicta.The slow-acting nature of cycloxaprid can be ascribed to the dynamic physiological responses observed in S. invicta post exposure, encompassing the enhancement of lipid and energy metabolism, alterations in immune function, and adjustments in detoxification mechanisms.These enhancements enable S. invicta to maintain their survival and recovery during the initial stages of exposure to cycloxaprid.However, they are insufficient to completely escape the lethal threat posed by the insecticide.As exposure duration increases, cycloxaprid progressively impairs physiological function, ultimately leading to the ants' death.Furthermore, the non-repellent characteristic of cycloxaprid is likely associated with changes in olfactory perception-related genes, which may induce a disorder in S. invicta's olfactory perception, resulting in reduced sensitivity to cycloxaprid.Consequently, the ants do not avoid their counterparts contaminated with cycloxaprid, thereby enhancing the effectiveness of this insecticide in controlling S. invicta at colony levels.

FIGURE 1
FIGURE 1 Transcriptomic analysis of Solenopsis invicta at three time points post cycloxaprid exposure.(A) The number of upregulated and downregulated differentially expressed genes (DEGs) at three time points; (B) Venn diagram illustrating the overlap of DEGs at the three time points, highlighting the shared and unique responses to cycloxaprid exposure.

FIGURE 2
FIGURE 2Comprehensive Gene Ontology (GO) enrichment analysis of DEGs at three distinct time points post cycloxaprid exposure.The size of each dot represents the number of genes involved in the corresponding GO term, while the shape symbolizes the specific functional classification of the respective GO term.p-values, calculated using hypergeometric tests, demonstrate the statistical significance of the enrichment results.

FIGURE 3
FIGURE 3Comprehensive Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of DEGs at three distinct time points post cycloxaprid exposure.The size of each dot represents the number of genes involved in the corresponding KEGG pathways.p-values, calculated using hypergeometric tests, demonstrate the statistical significance of the enrichment results.

FIGURE 4
FIGURE 4Expression patterns of 8 selected DEGs in Solenopsis invicta across three time points post cycloxaprid exposure.(A) Expression levels of the selected genes based on Fragments Per Kilobase of transcript per Million (FPKM) values obtained from RNA-seq data."DE" denotes differential expression, which has been validated through the use of "DESeq2" or "EdgeR" tests, indicating a statistically significant difference in gene expression between the control and experimental groups at the corresponding time point.(B) Relative expression levels of the selected genes were obtained via quantitative realtime PCR (qRT-PCR) analysis.Error bars represent the standard error of the mean (SEM) from three replicates."*" marks instances where the paired t-test yielded a p-value <0.05, suggesting a significant difference in gene expression between the control and experimental groups at the corresponding time point.CK and, EG indicate the control group and the experimental group (fire ants treated with cycloxaprid), respectively, followed by the time point of sampling.The consistency between the RNA-seq and qRT-PCR data highlights the reliability of our transcriptomic analysis.

FIGURE 5
FIGURE 5 Metabolomic analysis of Solenopsis invicta at three time points post cycloxaprid exposure.(A) Classification of detected metabolites; (B) Partial Least Squares Discriminant Analysis (PLS-DA) results for six sample groups, illustrating the progressive separation between the experimental and control groups across the three time points; (C) The number of upregulated and downregulated differentially expressed metabolites (DMs) at three time points; (D) Venn diagram illustrating the overlap of DMs at the three time points, highlighting the shared and unique responses to cycloxaprid exposure.