Comparative Transcriptome Analysis Reveals the Protective Mechanism of Glycyrrhinic Acid for Deoxynivalenol-Induced Inflammation and Apoptosis in IPEC-J2 Cells

Deoxynivalenol (DON) is the most common mycotoxin that frequently contaminates human food and animal feed, resulting in intestinal diseases and systemic immunosuppression. Glycyrrhinic acid (GA) exhibits various pharmacological activities. To investigate the protective mechanism of GA for DON-induced inflammation and apoptosis in IPEC-J2 cells, RNA-seq analysis was used in the current study. The IPEC-J2 cells were treated with the control group (CON), 0.5 μg/mL DON, 400 μg/mL GA, and 400 μg/mL GA+0.5 μg/mL DON (GAD) for 6 h. Results showed that 0.5 μg/mL DON exposure for 6 h could induce oxidative stress, inflammation, and apoptosis in IPEC-J2 cells. GA addition could specifically promote the proliferation of DON-induced IPEC-J2 cells in a dose- and time-dependent manner. In addition, GA addition significantly increased Bcl-2 gene expression (P < 0.05) and superoxide dismutase and catalase activities (P < 0.01) and decreased lactate dehydrogenase release, the contents of malonaldehyde, IL-8, and NF-κB (P < 0.05), the relative mRNA abundances of IL-6, IL-8, TNF-α, COX-2, NF-κB, Bax, and caspase 3 (P < 0.01), and the protein expressions of Bax and TNF-α. Moreover, a total of 1576, 289, 1398, and 154 differentially expressed genes were identified in CON vs. DON, CON vs. GA, CON vs. GAD, and DON vs. GAD, respectively. Transcriptome analysis revealed that MAPK, TNF, and NF-κB signaling pathways and some chemokines played significant roles in the regulation of inflammation and apoptosis induced by DON. GA may alleviate DON cytotoxicity via the TNF signaling pathway by downregulating IL-15, CCL5, and other gene expressions. These results indicated that GA could alleviate DON-induced oxidative stress, inflammation, and apoptosis via the TNF signaling pathway in IPEC-J2 cells.


Introduction
The gastrointestinal tract of animals plays an important role in the digestion and absorption of various nutrients and provides a pivotal barrier against food contaminants and other harmful irritants such as toxins, stress, and pathogens [1,2]. The injury gastrointestinal tract will cause barrier dysfunction, leading to inflammatory bowel diseases, necrotizing enterocolitis, and other intestinal diseases [3]. Deoxynivalenol (DON), called colloquially vomitoxin, is one of the most common mycotoxins mainly produced by Fusarium graminearum and Fusarium culmorum to contaminate human food and animal feed [4,5]. In general, the primary target of DON is the intestinal tract. Consumption of DONcontaminated food or feed brings about the compromised intestinal epithelial barrier function [6], the augment of intestinal permeability [7], disordered intestinal structure [8], and lower nutrient absorption [9], which leads to intestinal diseases and systemic immunosuppression. In addition, swine is the most sensitive domestic animal to DON. The IPEC-J2 cell line, isolated from the jejunal epithelium of the neonatal piglet, was widely used as an in vitro model for studying intestinal functions [10].
According to previous studies, DON exposure could induce cytotoxicity, oxidative stress, inflammation, and apoptosis and influence cell growth and functions [11][12][13]. DON triggers epithelial inflammation and oxidative stress by stimulating proinflammatory cytokine production in different epithelial cells as well as activates prototypical signaling pathways linked to immunity and inflammation such as the nuclear factor-(NF-) κB and mitogen-activated protein kinase (MAPK) [14,15]. The MAPK pathway plays an important role in regulating cell proliferation, inflammation, and immune reactions in the intestine, which consists of extracellular signal-regulated kinase (ERK), p38, and stressactivated protein kinase (JNK) pathways [16]. NF-κB is essential for the immune system, which can modulate cytokine expression and effector enzymes, as well as bind to various receptors involved in immunization [17]. External factors such as toxins and pathogens could stimulate the activation of the MAPK and NF-κB pathways to cause cell apoptosis [18]. Many researches have uncovered the toxic effects of DON on intestinal cells; however, the molecular mechanism underlying DON-induced inflammation and apoptosis of intestinal epithelial cells is not completely clear. Therefore, it is very imperative to fully reveal the toxic effects of DON on intestinal cells and to find effective nutritional and protective strategies for alleviating DON-induced intestinal damage to maintain animal health.
Traditional Chinese medicines have been recognized as good feed additives in improving animal production performance, immunity, and disease resistance [19,20]. Licorice, a perennial herb, has well-known detoxifying effects [21]. Glycyrrhinic acid (GA) is a botanical extract in licorice with various pharmacological activities including anti-inflammatory, immune regulation, antioxidation, antivirus, anticancer, and hypolipidemic effects [22][23][24]. Licorice extracts have been reported for clinical use in the treatment of liver injury [25] as well as for preventing diseases and promoting animal growth and meat quality [26]. Although many studies have confirmed the anti-inflammatory and liver protection effects of GA, the intestinal protective effect of GA against DON-induced toxicity, inflammation, and apoptosis has not been found.
To further compare the differential gene expression profiles and elucidate the repairing mechanisms of GA against inflammation and apoptosis in DON-stimulated IPEC-J2 cells, the high-throughput RNA sequencing (RNA-seq) was used in this study. The results will contribute to the identification of functional genes to provide a more comprehensive understanding of GA for suppressing inflammation and apoptosis in DON-stimulated IPEC-J2 cells. It will lay a theoretical basis for the application of GA as a feed additive for alleviating DON cytotoxicity to improve animal health and production performance.

Cell Culture and Treatments.
The IPEC-J2 cell line was obtained from the College of Animal Science and Technology, Jiangxi Agricultural University (Jiangxi, China) and cultured in H-G DMEM supplemented with 10% FBS, 100 U/mL penicillin, and 100 μg/mL streptomycin at 37°C in an incubator with 5% CO 2 . When the confluence of cells reached 80-90%, cells were seeded into 96-well or 6-well plates and cultured for 24 h prior to different treatments. Based on the previous research in our laboratory (Xu et al., 2020), 0.5 μg/mL DON treated for 6 h was selected to make the cell damage model. GA and DON were diluted with H-G DMEM without serum and antibiotics.

Cell
Viability. The effect of GA on alleviating DON cytotoxicity was determined using MTT assay. The cells (1 × 10 4 cells/well) were seeded in 96-well plates and cultured for 24 h and subsequently incubated with an increasing concentration of GA (0, 50, 100, 200, 400, and 800 μg/mL) for 3, 6, 12, and 24 h, respectively. After the optimal GA concentration was selected, the cells were incubated with 0.5 μg/mL DON and the optimal GA concentration for 6 h. Then, each well was added with 10 μL 5 mg/mL MTT and incubated for 4 h at 37°C. Thereafter, the supernatant was removed and 150 μL DMSO was added to each well to solubilize the formazan. The plate was shaken at room temperature for 10 min. The absorbance was measured at 490 nm wavelength with a reference wavelength of 630 nm by an ELx800 microplate reader (BioTek Instruments Inc., Winooski, VT, USA).

Data
Analysis of RNA Sequencing. The quality of raw reads was evaluated and cleaned by trimming the 3′ adapter sequence using cutadapt [27]. The raw reads with an average sequencing accuracy below 99.9% were removed. The clean reads were aligned to the reference porcine genome assembly Sus scrofa 11.1 using HISAT2 with default parameters [28]. The mapping results of different genomic regions were eval-uated. Read counts for each gene were calculated using HTSeq with union strategy [29]. In order to avoid the effects of sequencing coverage and gene length, read counts were normalized into fragments per kilobase of exon model per million mapped fragments (FPKM), which were used as input for the following analysis. Pearson correlation coefficient was calculated between different samples using the expression values of all detected genes. Moreover, principal component (PC) analysis was also performed among all samples using all detected genes. Differential expression analysis was performed on DESeq2 (Bioconductor version 1.6.2), and genes with the threshold of P < 0:05 and the absolute value of log 2 ðfold changeÞ > 1 among two groups were regarded as the differentially expressed genes (DEGs).

GO Functional Annotation and KEGG Pathway
Enrichment Analyses. In order to deeply gain insights into the functions of DEGs obtained from different treatments, Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were carried out. GO annotation including biological process (BP), molecular function (MF), and cellular component (CC) was performed using the topGO package.
The significance of genes annotated in a GO term was calculated by hypergeometric distribution. The enrichment of genes involved in the same biological functions facilitated the understanding of gene biological functions. KEGG was a pathway-related database for understanding highlevel functions and utilities of the biological systems, and a cluster profile package was used for this enrichment.

PPI Network Construction and Module Analyses.
To establish protein-protein interaction (PPI) network of interested DEGs, the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING, version 11.0, https://string-db.org/) database was used, which covered interactive relationships of proteins obtained from experiments, data mining, and homological prediction for various species [30]. The minimum required interaction score was set as high (0.90) for two different DEGs. Then, the Cytoscape tools (v3.6.1, http://cytoscape .org/) was used for the visualization of PPI networks. For seeking hub genes and significant network, the cytoHubba and Molecular Complex Deletion (MCODE) plug-ins for Cytoscape were downloaded to determine hub genes using the Maximal Clique Centrality (MCC) ranking method and significant network modules with default parameters. 2.11. Western Blot Analysis. After four treatments, total proteins of the cells were extracted with RIPA buffer (EpiZyme Biotechnology, Shanghai, China), and concentrations were determined by the BCA protein assay kit (EpiZyme Biotechnology, Shanghai, China). The equal amounts of proteins were separated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and transferred onto the polyvinylidene difluoride (PVDF) membranes. The membranes were blocked with TBST buffer containing 5% skim milk powder for 2 h, followed by incubation with the primary antibodies at 4°C overnight, and then further incubated with the secondary antibodies for 2 h at room temperature. The immunopositive bands were visualized by enhanced chemiluminescence.
2.12. Statistical Analysis. The data were expressed as means ± SD. There were 3 replications for each treatment. Statistical analysis was analyzed using the one-way ANOVA method (Duncan's multiple comparison test) with the SPSS 20.0 software (Sishu Software, Shanghai Co., Ltd., Shanghai, China). All graphs were generated using GraphPad Prism 8. Statistical differences were considered significant at P < 0:05.

GA Alleviated Cytotoxicity and Oxidative Stress in the DON-Induced IPEC-J2 Cells.
To study the effects of GA on IPEC-J2 cells, the cell viability evaluation was performed using MTT assay. Figure 1(a) shows that GA promoted cell proliferation in a dose-and time-dependent manner, and cell viabilities were significantly increased when cells were treated with GA for 6 h (P < 0:01). Therefore, 6 h was set as the optimal reaction time to confirm the effects of GA on alleviating DONinduced cells. As shown in Figure 1(b), the treatment with 0.5 μg/mL DON for 6 h significantly decreased cell viability (P < 0:01); however, 400 μg/mL GA addition could markedly enhance cell viability induced by DON exposure (P < 0:01). Therefore, experimental condition with 400 μg/mL GA and 6 h reaction time was employed in the subsequent experiments.
The results in Figures 1(c)-1(f) showed that DON at a concentration of 0.5 μg/mL for 6 h significantly increased LDH release and malonaldehyde (MDA) levels and decreased superoxide dismutase (SOD) and catalase (CAT) activities, compared with the control group (P < 0:01). However, GA addition significantly decreased LDH release and MDA level (P < 0:01) and increased the activities of SOD (P < 0:01) and CAT (P < 0:05). The above results suggested that DON exposure could cause serious cell damage, and GA could efficiently alleviate the cytotoxicity and oxidative stress in DON-induced cells.

GA Alleviated IPEC-J2 Cell Apoptosis and Inflammation
Induced by DON. The results of Annexin V/FITC/PI staining in Figures 2(a) and 2(b) showed that the total apoptotic rate of IPEC-J2 cells induced by DON was significantly increased compared with the control group (P < 0:01), whereas GA addition significantly decreased it (P < 0:01). Figure 2(c) shows that DON significantly decreased the viable cell rates and increased the late and early apoptotic cell rates (P < 0:01); however, GA addition significantly increased the viable cell rates (P < 0:01) and decreased the late apoptotic cell rates (P < 0:01) and early apoptotic cell rates (P < 0:05).
3.3. RNA-seq Data Analyses. In order to deeply uncover the molecular mechanism of GA for alleviating IPEC-J2 cell apoptosis and inflammation induced by DON, RNA-seq was performed in four groups including CON, DON, GA, and GAD. The sequencing quality of Q20 (base error < 1%) and Q30 (base error < 0:1%) was 97.96% and 94.51%, respectively. After removing reads containing poly-N and adapters and filtering low-quality reads, approximately 42 million clean reads (93.59%) were kept, in which the reads with an average of 96.64% were mapped to the reference genome (Supplementary Table 2). The correlation plot showed distinct groups for different samples using unsupervised hierarchical clustering. After PCA analysis, the first two PC were used for evaluating the relationship of 12 samples according to their gene expression profiles. PC1 and PC2 explained 86% and 10% of the total variations, respectively, and the samples in different groups were gathered tightly (Figure 4(a)). These results showed a satisfactory data quality of RNA-seq, and three biological replications in each treatment displayed a homogeneity regarding whole gene expressions.
An expression matrix including 21278 genes from RNA-seq was obtained. A total of 1576 DEGs (668 upregulated and 908 downregulated), 289 DEGs (211 upregulated and 78 downregulated), and 1398 DEGs (734 upregulated and 664 downregulated) were identified in the DON, GA, and GAD groups, compared with the control group, respectively. Figure 4(b) shows the upregulated and downregulated genes in CON vs. DON, CON vs. GA, CON vs. GAD, and DON vs. GAD, respectively. Subsequently, the DEG clusters were analyzed in a heat map, which showed a more similarity among samples from the same treatment based on those genes (Figure 4(c)). Three biological replications in each treatment were clustered together with almost identical distance using DEGs by hierarchical clustering, indicating good reproducibility of each treatment. In  (Table 1).

GO Annotation and KEGG Enrichment.
To understand the biological processes and pathway of DEGs, they were functionally annotated and enriched into GO and KEGG databases, respectively. GO analysis were performed regarding three main ontologies: biological process (BP), molecular function (MF), and cellular component (CC). The top 10 GO terms for DEGs identified from CON vs. DON, CON vs. GA, CON vs. GAD, and DON vs. GAD are summarized in Figure 5(a). In the BP level, these four clusters of DEGs were mainly annotated as system development and cell differentiation in the top 10 terms; however, some DEGs from CON vs. GA were annotated as a cellular response to a stimulus, and some DEGs from CON vs. GAD were annotated as the immune system process. In the MF level, almost all the DEGs were annotated as  Oxidative Medicine and Cellular Longevity sequence-specific DNA binding except for DON vs. GAD. In the CC level, DEGs were intensively annotated as intracellular organelle in CON vs. DON and CON vs. GAD; however, almost all the DEGs were annotated as an extracellular region or cell membrane in CON vs. GA, which showed a similar annotation as DON vs. GAD to some extent.
KEGG pathway enrichments were conducted for DEGs in four different comparisons mentioned above to identify pathways that changed significantly under different experimental conditions. All the DEGs were enriched according to cellular processes, environmental information processing, human diseases, metabolism, and organismal systems (Figures 5(b) and 5(c)). The DEGs in the CON vs. DON group were mainly involved in the following three pathways: (1) many signal transduction pathways including TNF, Notch, MAPK, Ras, and NF-κB; (2) some immune-related signal pathways including Th1 and Th2 cell differentiation, B cell receptor, and Toll-like receptor; (3) three metabolic pathways including steroid biosynthesis, biosynthesis of unsaturated fatty acids, and terpenoid backbone biosynthesis (Supplementary Table 3). In CON vs. GA, the signal pathways included cell adhesion molecules (CAMs), leukocyte transendothelial migration, complement and coagulation cascades, and PI3K-  Table 4). There were more than 900 DEGs to be shared between CON vs. DON and CON vs. GAD; their enriched pathways were similar, which was related to signal transduction and immune systems (Supplementary  Table 5). Interestingly, some same pathways such as TNF, cAMP, and leukocyte transendothelial migration were identified in DON vs. GAD and CON vs. DON, in which the Jak-STAT signaling pathway was unique. In addition, immune systems including the chemokine signaling pathway, Fc epsilon RI signaling pathway, complement and coagulation cascades, and Fc gamma R-mediated phagocytosis were also identified (Supplementary Table 6).
The functional annotation and enrichment among the 79 coexpressed DEGs in CON vs. DON and DON vs. GAD were investigated. GO annotation showed that most annotated BP ontology focused on immune cell migration such as myeloid   Oxidative Medicine and Cellular Longevity Oxidative Medicine and Cellular Longevity leukocyte, leukocyte, neutrophil, mononuclear, and granulocyte. Furthermore, some DEGs were annotated as the external stimulus, inflammatory response, Notch binding, receptor binding, chemokine, chemokine receptor, IL-8 (CXCL8) receptor binding, and so on in the MF aspect (Supplementary Figure 1). After KEGG enrichment, DEGs including CXCL8, CCL5, and IL-15 were enriched in the chemokine and TNF signaling pathway and cytokine-cytokine receptor interaction pathway (Supplementary Table 7).

PPI Network Construction and Module Analysis.
To further study the interaction among the DEGs, PPI networks were constructed in CON vs. DON. As shown in Figure 6(a), a network including 209 nodes and 572 relationship pairs was constructed, and genes including CXCL8, CCL4, IL-15, IL-6, IL1A, NFKB2, NFKB1A, and MAPK-related were highly interacted. The top 30 highly connected genes (hub genes) were searched out by cyto-Hubba, and the first 28 hub genes formed two independent networks, which were the top 2 significant modules including 15 nodes and 105 edges as well as 13 nodes and 78 edges extracted from Figure 6(a) by MCODE with scores of 15 and 13, respectively (Figure 6(b)). In order to investigate the interaction among genes involved in the alleviation of GA, the 79 DEGs coexpressed in both CON vs. DON and DON vs. GAD were used to construct the PPI network, and one network was obtained (Figure 6(c)). CXCL8 was the most hub gene interacted with a range of genes including CCL5, PTGS2, IL-15, HCAR1, RAC2, S1PR1, PECAM1, and LTB4R.

Discussion
Deoxynivalenol is the most prevalent mycotoxin, which is widespread and frequently exists in most human foods and animal feeds, resulting in intestinal diseases, systemic immunosuppression, and other diseases [31]. Pigs are more sensitive to DON than other kinds of animals. DON can lead to emesis, diarrhea, low nutrient absorption and growth, immune system disorders, and economical loss [32,33]. It has been indicated that DON exposure can induce cytotoxicity, oxidative stress, apoptosis, and intestinal barrier dysfunction [34,35]. It is also a ribotoxic mycotoxin to provoke oxidative stress and inflammatory responses [36]. Although the cytotoxic effects of DON have been well documented, few studies are conducted to alleviate its toxicity. Therefore, there is an urgent need to develop effective nutritional and protective strategies to alleviate the damage induced by DON and to improve the intestinal health of animals. GA has characters of anti-inflammation and liver protection, which is widely used to promote the growth and development in animals, but the protective effects of GA on alleviating DONinduced intestinal inflammation, apoptosis, and its active mechanism are rarely reported. Hence, this study provided the potential mechanism of GA for alleviating inflammation and apoptosis in DON-stimulated IPEC-J2 cells using molecular biochemistry techniques and RNA-seq analysis, which supported the application of GA as one kind of feed additive for alleviating DON cytotoxicity in animal feeding.

11
Oxidative Medicine and Cellular Longevity The recent study has revealed that DON can induce cellular oxidative stress and decrease the activities of antioxidant enzymes [37]. SOD and CAT are primary enzymes participating in repairing damage caused by oxidative stress. SOD can effectively eliminate O 2to protect the body from the influence of superoxide anion, and CAT can decompose ROS produced in the metabolic process [38]. It was reported that DON could alter membrane integrity, cellular redox signaling, and antioxidant status of the cells [39]. MDA is the important product of membrane lipid peroxide to aggravate the damage of cell membrane, which can indirectly reflect the membrane integrity of cells [40]. This study showed that DON exposure significantly increased LDH release and MDA level and decreased the SOD and CAT activities, indicating DON cytotoxicity, in agreement with the previous report [11]; however, GA addition significantly increased the activities of these antioxidant enzymes and decreased the LDH release and MDA level, inferring that GA could alleviate the cytotoxicity and oxidative stress induced by DON.
In addition, the related proapoptotic and proinflammatory factors such as Bax, caspase 3, IL-6, IL-8, TNF-α, COX-2, and NF-κB were significantly upregulated after DON stimulation in IPEC-J2 cells, which is in accordance with the previous study [11]. Proinflammatory cytokines such as IL-6, IL-8, and TNF-α are associated with the severity of inflammation and intestinal disease, which can directly affect intestinal epithelial cells. Excessive expression of them can result in damage to the epithelial barrier and secretion of chemokines such as COX-2 and NF-κB, even causing epithelial cell apoptosis [41]. COX-2 is an immediate-early response gene, which is regarded as an inflammatory marker associated with a variety of inflammatory conditions and cancer [42]. NF-κB plays a critical role in regulating inflammatory and immune responses as well as controlling cell proliferation and cell death [43]. Moreover, apoptosis is a genetically programmed procedure due to its central role in normal cell development and homeostasis at various physiological and pathological conditions [44,45]. Bax, Bcl-2, and caspase 3 are the important mediators of cell apoptosis. Therefore, the upregulation of these genes in this study clearly indicated that DON exposure induced inflammation and apoptosis in IPEC-J2 cells. However, these proinflammatory and proapoptosis genes were downregulated by GA addition. In addition, the contents of IL-8, caspase 3, and NF-κB in cell supernatant were increased after DON exposure, while GA could reduce their contents except for caspase 3. Furthermore, Annexin V-FITC/PI staining further con-firmed that GA significantly decreased cell apoptosis induced by DON. The consistencies among the cytokine concentrations in cell supernatant, cytokine mRNA abundances, and protein expressions in this study prove the cytoprotective role of GA. Based on the above results, it was preliminary speculated that GA could act as a cytoprotective agent against DON-induced cytotoxicity, oxidative stress, inflammation, and apoptosis in IPEC-J2 cells.
To further characterize the potential mechanisms of GA for alleviating cell damage induced by DON, the RNA-seq was adopted to compare the differential gene expression profiles of different treatment groups in this study. The results showed that there were different numbers of DEGs during comparisons among 4 groups. GO and KEGG enrichment analyses of the DEGs were performed to analyze the functions of genes. It was found that the DEGs are primarily related to immune and inflammatory response. Cytokines and chemokines play critical parts in inflammation and immune response. Cytokines can directly promote or limit intestinal cell proliferation, permeability, and apoptosis, thereby destroying the intestinal epithelial barrier [46]. Chemokines are responsible for chemotactic cell migration and play important roles in the humoral and cellular immune responses to induce the production of corresponding inflammatory factors [47,48]. Except for the proinflammatory factors, DON can also upregulate several chemokines. In the current study, the differentially expressed chemokines such as CCL4, CCL5, and CXCL8 were identified in the CON vs. DON. In general, CCL4 acts as a chemoattractant for a variety of immune cells [49], and CCL5 is a key proinflammatory chemokine [50], which can be produced by many kinds of cells including CD4 + T lymphocytes,

13
Oxidative Medicine and Cellular Longevity epithelial cells, fibroblasts, and platelets [51]. CXCL8 can be induced by proinflammatory cytokines TNF-α [52]. It is evident that DON can activate the NF-κB signaling pathway to trigger MAPK signaling pathways, which are essential for controlling inflammation by mediating the production of inflammatory factors [13]. A recent study demonstrated that inflammation occurred by activating the P38 MAPK and Erk1/2 pathway by comparing transcriptomes with and without DON treatment (2 μg/mL for 2 h) in IPEC-J2 cells [53]. Although the same MAPK, TNF, NF-κB signaling pathways, and cytokine-cytokine receptor interaction were enriched in this study, the number of shared DEGs in two studies was small when using the same DEG identification threshold. A total of 138 shared DEGs including CXCL8, IL1A, IL-6, FOS, and MAP3K5 were identified, in which only eight genes including FOS and IL1A showed an inconsistent expression state, mirroring the validity of results in both studies. The low proportion of overlapped DEGs suggests that there are some different pathways involved in inflammation induced by DON. Therefore, these findings highlight that TNF, MAPK, and NF-κB signaling pathways and chemokines are involved in the inflammation and immune response of IPEC-J2 cells induced by DON.
GA has many functions of stimulating the antioxidant status and reinforcing the immune system function [54], which is also validated by RNA-seq analysis in this study. Most DEGs in CON vs. GA were annotated as the GO terms of cell development and the regulation of stimulus. The KEGG pathway enrichment analysis demonstrated that DEGs were significantly enriched into cell adhesion molecules and the immune system. Some pathways such as TNF and MAPK were insignificantly enriched in CON vs. GA. The similar pathways enriched in CON vs. DON were also enriched in CON vs. GAD, indicating that GA was involved in some of the same pathways to alleviate the damage from DON-induced IPEC-J2 cells. To uncover the potential pathway involved in the protection effect of GA, the DEGs in DON vs. GAD were investigated. Some of the 154 DEGs including CXCL8, IL-15, CCL5, and RAC2 in DON vs. GAD were enriched in chemokine, cytokine-cytokine receptor interaction, TNF, cAMP, and Jak-STAT signaling pathways, which were also included in the pathways enriched from CON vs. DON. It is demonstrated that chemokine, cytokine, and TNF signaling pathways may act as major contributors for the protection effect of GA.
This research showed that 79 DEGs coexpressed in both CON vs. DON and DON vs. GAD might be involved in the protection effect of GA. Among these DEGs, CXCL8, CCL5, and IL-15 interacted with each other in a high strength of data support were enriched to cytokine-cytokine receptor interaction and the TNF signaling pathway. This result further indicated that the TNF signaling pathway may participate the protection effect of GA. The 26 DEGs upregulated in the DON group and downregulated in the GAD group belonging to IL-15 and CCL5 enriched in the TNF signaling pathway and cytokine-cytokine receptor interaction. Furthermore, CCL5 and RAC2 are enriched in the chemokine signaling pathway, and RAC2 and HCAR1 are enriched in the cAMP signaling pathway. Particularly, CCL5 interacted with IL-15, ISG15, PTGS2, HCAR1, and CXCL8. CCL5 and IL-15 may act as key genes in the whole alleviation process of GA. IL-15 is a human neutrophil agonist known to induce RNA and de novo protein synthesis, phagocytosis, and apoptosis, which is also often regarded as a proinflammatory cytokine [55]. A previous study demonstrated that IL-15 can delay the apoptosis of neutrophil by recovering the expression of the antiapoptotic MCL1 protein via several kinases including Jak-2, p38, and ERK1/2 [56]. ISG15, a member of the family of ubiquitin-like proteins, is an important component of host responses to microbial infection [57]. For other DEGs, RAC2 and TXNIP are involved in oxidative stress, and LTB4R takes part in inflammatory and immune response. These results show that GA may alleviate the cell damage from DON via the TNF signaling pathway by downregulating IL-15 and CCL5 gene expressions, and other immunization-related and inflammation-related genes including ISG15, RAC2, TXNIP, and LTB4R also play important roles in the process.

Conclusion
DON exposure can induce oxidative stress, inflammation, and apoptosis in IPEC-J2 cells, which can be alleviated by GA addition. MAPK, TNF, and NF-κB signaling pathways and some chemokines serve as major parts to participate in the inflammation and apoptosis induced by DON; however, GA shows alleviating effect via the TNF signaling pathway by downregulating IL-15, CCL5, and other immunization-related and inflammation-related gene expressions (Figure 7). This study provides a new insight into the protective mechanism of GA against the damage induced by mycotoxin exposure in intestinal epithelial cells of pigs, which also lays a theoretical foundation for the development and application of GA as a potential feed additive to alleviate DON-induced cytotoxicity for improving animal health and productions.

Data Availability
The data used to support the findings of the present study are available from the corresponding authors upon request. edited the manuscript. Fushan Lu helped in the methodology and data curation. Xiaoxiang Xu and Guorong Yan contributed equally to this work.