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

Xu, Xiaoxiang; Yan, Guorong; Chang, Juan; Wang, Ping; Yin, Qingqiang; Liu, Chaoqi; Zhu, Qun; Lu, Fushan

doi:https://doi.org/10.1155/2020/5974157

Oxidative Medicine and Cellular Longevity

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 5974157 | https://doi.org/10.1155/2020/5974157

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

Xiaoxiang Xu,¹Guorong Yan,²Juan Chang,¹Ping Wang,¹Qingqiang Yin,¹Chaoqi Liu,¹Qun Zhu,³and Fushan Lu⁴

Academic Editor: Víctor M. Mendoza-Núñez

Received24 Jul 2020

Revised17 Sept 2020

Accepted05 Oct 2020

Published24 Oct 2020

Abstract

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 () and superoxide dismutase and catalase activities () and decreased lactate dehydrogenase release, the contents of malonaldehyde, IL-8, and NF-κB (), the relative mRNA abundances of IL-6, IL-8, TNF-α, COX-2, NF-κB, Bax, and caspase 3 (), 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.

1. 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 DON-contaminated 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–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 stress-activated 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–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.

2. Materials and Methods

2.1. Materials and Reagents

Glycyrrhinic acid was obtained from Luoyang Lansealy Technology Co., Ltd. (Luoyang, Henan, China). DON was purchased from Sigma-Aldrich (St. Louis, MO, USA), dissolved in dimethylsulfoxide (DMSO) (Beijing Solarbio Biotechnology Co., Ltd., Beijing, China), formulated into 1 mg/mL stock solution, and stored at -20°C. Phosphate-buffered saline (PBS), penicillin- (10000 U/mL) streptomycin (10 mg/mL), 0.25% pancreatin with or without ethylenediaminetetraacetic acid (EDTA), and 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) were purchased from Beijing Solarbio Biotechnology Co., Ltd., Beijing, China. High-glucose Dulbecco’s Modified Eagle’s Medium (H-G DMEM) and fetal bovine serum (FBS) were purchased from Biological Industries (Kibbutz Beit-Haemek, Israel). The Annexin V-FITC/PI kit was purchased from Shanghai Qihai Futai Biotechnology Co., Ltd., Shanghai, China. The IL-8, NF-κB, and caspase 3 concentrations assay kits were purchased from Jiangsu Meimian Industrial Co., Ltd., Jiangsu, China. Trizol reagent was purchased from Invitrogen (Grand Island, NY, USA). The reverse transcription kit and TB GREEN kit were purchased from TaKaRa Bio Inc. (Dalian, China). All cell culture materials were purchased from Corning Costar (Cambridge, NY, USA). Rabbit polyclonal antibody of Bax (abs119724) and TNF-α (abs123966) and goat anti-rabbit antibody of IgG were purchased from Absin Bioscience Inc., Shanghai, China. The β-actin was purchased from Bioworld Technology Inc., Nanjing, China.

2.2. 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₂. 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.

2.3. Cell Viability

The effect of GA on alleviating DON cytotoxicity was determined using MTT assay. The cells ( 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).

2.4. Annexin V-FITC/PI Apoptosis Assay

The apoptosis cells were examined using the Annexin V-FITC apoptosis detection kit. The cells ( cells/well) were seeded in a 6-well plate and incubated for 24 h and then subsequently 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. After four treatments, the cells were digested by pancreatin without EDTA and centrifuged at 1000 rpm for 5 min to remove the medium. Then, the precipitation was washed twice with PBS and resuspended in 100 μL 1x binding buffer and stained with 5 μL Annexin V-FITC and 5 μL PI at room temperature in the dark for 15 min. Finally, 400 μL 1x binding buffer was added to each tube, and cell apoptosis levels were analyzed by flow cytometry (Em 530 nm, Ex 488 nm, Becton Dickinson Company, NJ, USA).

2.5. Determinations of IL-8, NF-κB, and Caspase 3

IPEC-J2 cells ( cells/well) were seeded in 6-well plates and cultured for 24 h. Following four treatments (control, 0.5 μg/mL DON, 400 μg/mL GA, and 400 μg/mL GA+0.5 μg/mL DON) for 6 h, the cell supernatants were collected to determine IL-8, NF-κB, and caspase 3 contents, and the cells were washed twice with PBS to extract the proteins for further use. The IL-8, NF-κB, and caspase 3 concentration assays were quantified by enzyme-linked immunosorbent assay (ELISA) in accordance with the manufacturer’s instructions.

2.6. RNA Extraction, Library Construction, and Sequencing

IPEC-J2 cells ( cells/well) were seeded in 6-well plates and cultured for 24 h. Following four treatments (control, 0.5 μg/mL DON, 400 μg/mL GA, and 400 μg/mL GA+0.5 μg/mL DON) for 6 h, the cells were washed twice with PBS and then harvested. Total RNA was isolated using Trizol (Invitrogen Corporation, Grand Island, NY, USA) according to the manufacturer’s instructions and stored at -80°C. The RNA quantity and purity were assessed on an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA), and the RNA integrity was checked by agarose gel electrophoresis. The purified mRNA was enriched by Oligo(dT) beads and fragmented into 300 bp by Mg²⁺ ion fragmentation buffer. Then, the fragmented mRNA was reversely transcripted into cDNA with 6 bp random primers. After that, the second-strand cDNA was synthesized using the first-strand cDNA as a template. The synthesized product was purified with a PCR purification kit (Qiagen, Venlo, The Netherlands), followed by end repair, poly(A) addition, and Illumina sequencing adapter ligation. The sequencing libraries were made, and the Agilent 2100 Bioanalyzer was used to check their sizes. Finally, sequencing libraries were sequenced using the Illumina HiSeq 4000 platform (Illumina, San Diego, CA, USA) with a paired-end (PE 150 bp) sequencing strategy.

2.7. 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 evaluated. 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 and the absolute value of among two groups were regarded as the differentially expressed genes (DEGs).

2.8. 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 high-level functions and utilities of the biological systems, and a cluster profile package was used for this enrichment.

2.9. 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.10. Real-Time Quantitative PCR

Four treatments and the RNA extraction were carried out as above. About 1 μg total RNA was used to erase gDNA at 42°C for 2 min and then reversely transcripted into cDNA using the TB GREEN kit (TaKaRa Bio Inc., Dalian, China). Real-time PCR was performed with the CFX Connect™ Real-Time PCR Detection System (Bio-Rad, Hercules, CA, USA) in a 20 μL PCR reaction system containing 2 μL cDNA, 1 μL reverse primer,1 μL forward primer, 10 μL SYBR mixture, and 6 μL deionized water. The GAPDH was used as a housekeeping gene to normalize the gene levels, and the RT-PCR data were analyzed using the 2^-ΔΔCT method (Livak and Schmittgen, 2001). All the primers used in this study were synthesized by Shanghai Shenggong Biotech Co., Ltd., Shanghai, China (Supplementary Table 1).

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 . 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 .

3. Results

3.1. 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 (). Therefore, 6 h was set as the optimal reaction time to confirm the effects of GA on alleviating DON-induced cells. As shown in Figure 1(b), the treatment with 0.5 μg/mL DON for 6 h significantly decreased cell viability (); however, 400 μg/mL GA addition could markedly enhance cell viability induced by DON exposure (). Therefore, experimental condition with 400 μg/mL GA and 6 h reaction time was employed in the subsequent experiments.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 1

Cell viability and status affected by GA and DON. (a) Effects of different GA concentrations on cell viability at different incubation times (). (b) Effects of 0.5 μg/mL DON and different GA concentrations on cell viability for 6 h incubation (). (c–f) Effects of GA on LDH release and antioxidant indexes in DON-induced cells (). DON: the single DON (0.5 μg/mL) group; GA: the single GA (400 μg/mL) group; GA+DON: 400 μg/mL GA+0.5 μg/mL DON group. indicates the significant difference, compared with the control group; ^# and ^## indicate the significant difference, compared with the DON group; “ns” indicates the insignificant difference, compared with the DON group.

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 (). However, GA addition significantly decreased LDH release and MDA level () and increased the activities of SOD () and CAT (). 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.

3.2. 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 (), whereas GA addition significantly decreased it (). Figure 2(c) shows that DON significantly decreased the viable cell rates and increased the late and early apoptotic cell rates (); however, GA addition significantly increased the viable cell rates () and decreased the late apoptotic cell rates () and early apoptotic cell rates ().

(a)

(b)

(c)

Figure 2

GA alleviating cytotoxicity of DON-induced IPEC-J2 cells. (a) Flow cytometry analysis of Annexin V/FITC/PI staining cells; LACR, NCR, VCR, and EACR represent the late apoptotic cell rates, necrotic cell rates, viable cell rates, and early apoptotic cell rates, respectively. (b) Quantification of the total cell apoptotic rates (early apoptotic cell rates+late apoptotic cell rates). (c) Quantification and analysis of the late apoptotic cell rates, necrotic cell rates, viable cell rates, and early apoptotic cell rates in four groups. DON: the single DON (0.5 μg/mL) group; GA: the single GA (400 μg/mL) group; GA+DON: 400 μg/mL GA+0.5 μg/mL DON group. indicates the significant difference, compared with the control group; ^# and ^## indicate the significant difference, compared with the DON group.

As shown in Figures 3(a) and 3(b), DON exposure significantly increased the contents and the relative mRNA abundances of IL-8, caspase 3, and NF-κB (), upregulated the relative mRNA abundances of Bax, IL-6, TNF-α, and COX-2 (), and downregulated the Bcl-2 mRNA abundances (). However, compared with the DON group, GA addition significantly decreased the contents of IL-8 () and NF-κB (), downregulated the mRNA abundances of IL-8, NF-κB, Bax, IL-6, TNF-α, and COX-2 (), and upregulated the Bcl-2 mRNA abundance (). Figure 3(c) shows that the protein expressions of Bax and TNF-α were increased by DON exposure but decreased by GA addition, corresponding with the changes of mRNA abundances. It was summarized that GA addition could alleviate apoptosis and inflammation in DON-induced IPEC-J2 cells.

(a)

(b)

(c)

Figure 3

GA alleviating IPEC-J2 cell apoptosis and inflammation induced by DON. (a) The contents of IL-8, caspase 3, and NF-κB. (b) Effects of GA on regulating the relative mRNA abundances of apoptotic and inflammatory genes in DON-induced IPEC-J2 cells. (c) The relative protein expressions of Bax and TNF-α in four groups. DON: the single DON (0.5 μg/mL) group; GA: the single GA (400 μg/mL) group; GA+DON: 400 μg/mL GA+0.5 μg/mL DON group. and indicate the significant difference, compared with the control group; ^# and ^## indicate the significant difference, “ns” indicates , compared with the DON group.

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 () and Q30 () 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.

(a)

(b)

(c)

(d)

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 addition, the coexpressed DEGs in DON vs. CON, GA vs. CON, GAD vs. CON, and GAD vs. DON were analyzed. As shown in the Venn diagram (Figure 4(d)), there were 79 coexpressed DEGs in CON vs. DON and DON vs. GAD. Interestingly, there were 37 downregulated DEGs in CON vs. DON, but they were upregulated in DON vs. GAD; there were 26 upregulated DEGs in CON vs. DON, but they were downregulated in DON vs. GAD (Table 1).

3.4. 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 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.

(a)

(b)

(c)

Figure 5

Functional annotation and pathway enrichment results. (a) Gene Ontology (GO) annotation of DEGs in four different comparisons. The top 10 functional classified GO terms of DEGs annotated by the subontology of GO analysis including biological processes (BP), molecular function (MF), and cellular components (CC). (b, c) Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment of DEGs for four different comparisons. (b) Top 10 classifications regarding different KEGG pathways. (c) Scatterplot of top 20 pathways in KEGG enrichment.

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-Akt (Supplementary 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 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).

3.5. 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 cytoHubba, 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.

(a)

(b)

(c)

Figure 6

Protein-protein interaction (PPI) network construction of DEGs. (a) The PPI network of DEGs identified from CON vs. DON. (b) The top 30 highly connected DEGs with 19 edges of the PPI network and the top 2 significant modules extracted from the PPI network (MCODE scores are 15 and 13). (c) PPI network constructed from 79 coexpressed genes between CON vs. DON and DON vs. GAD. Red represents gene expression upregulated, and blue represents gene expression downregulated in (a) and (b). Node size indicates the number of genes interacted in (a) and (c). Line thickness indicates the strength of data support in (c). The interaction score was set to 0.9 in (a) and (b) and 0.4 for (c).

4. 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 DON-induced 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.

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^2- to 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 confirmed 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, 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.

5. 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.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Xiaoxiang Xu helped in the investigation, methodology, data curation, and writing of the original draft. Guorong Yan contributed in the data curation, formal analysis, and software. Juan Chang performed project administration, investigation, and methodology. Ping Wang, Chaoqi Liu, and Qun Zhu helped in the methodology and investigation. Qingqiang Yin performed funding acquisition, project administration, methodology, and investigation and wrote, reviewed, and edited the manuscript. Fushan Lu helped in the methodology and data curation. Xiaoxiang Xu and Guorong Yan contributed equally to this work.

Acknowledgments

This work was funded by the Natural Science Foundation of Henan Province (182300410029) and the Xinxiang Key Scientific and Technological Project (ZD19005).

Supplementary Materials

Supplementary Table S1: the detailed summary of sequencing data for each sample. Supplementary Table S2: the significant enriched KEGG pathways of the DEGs in CON vs. DON. Supplementary Table S3: the significant enriched KEGG pathways of the DEGs in CON vs. GA. Supplementary Table S4: the significant enriched KEGG pathways of the DEGs in CON vs. GAD. Supplementary Table S5: the significant enriched KEGG pathways of the DEGs in DON vs. GAD. Supplementary Table S6: the significant enriched KEGG pathways of the DEGs in CON vs. DON and DON vs. GAD. Supplementary Table S7: primer sequences of some genes for quantitative real-time PCR. Supplementary Figure S1: Gene Ontology (GO) annotation of DEGs in CON vs. DON and DON vs. GAD. The top 10 functional classified GO terms of DEGs annotated by subontology of GO analysis; BP: biological processes; MF: molecular function; CC: cellular components (CC). (Supplementary Materials)

References

J. D. Soderholm and M. H. Perdue, “II. Stress and intestinal barrier function,” American Journal of Physiology. Gastrointestinal and Liver Physiology, vol. 280, no. 1, pp. G7–G13, 2001.
View at: Publisher Site | Google Scholar
H. Clevers, “The intestinal crypt, a prototype stem cell compartment,” Cell, vol. 154, no. 2, pp. 274–284, 2013.
View at: Publisher Site | Google Scholar
M. Camilleri, “Leaky gut: mechanisms, measurement and clinical implications in humans,” Gut, vol. 68, no. 8, pp. 1516–1526, 2019.
View at: Publisher Site | Google Scholar
Z. Ren, C. Guo, S. Yu et al., “Progress in mycotoxins affecting intestinal mucosal barrier function,” International Journal of Molecular Sciences, vol. 20, no. 11, p. 2777, 2019.
View at: Publisher Site | Google Scholar
Q. Wu, K. Kuča, H.-U. Humpf, B. Klímová, and B. Cramer, “Fate of deoxynivalenol and deoxynivalenol-3-glucoside during cereal-based thermal food processing: a review study,” Mycotoxin Research, vol. 33, no. 1, pp. 79–91, 2017.
View at: Publisher Site | Google Scholar
C. Nossol, P. Landgraf, S. Kahlert et al., “Deoxynivalenol affects cell metabolism and increases protein biosynthesis in intestinal porcine epithelial cells (IPEC-J2): DON increases protein biosynthesis,” Toxins (Basel), vol. 10, no. 11, p. 464, 2018.
View at: Publisher Site | Google Scholar
P. Pinton, C. Braicu, J. P. Nougayrede, J. Laffitte, I. Taranu, and I. P. Oswald, “Deoxynivalenol impairs porcine intestinal barrier function and decreases the protein expression of claudin-4 through a mitogen-activated protein kinase-dependent mechanism,” The Journal of Nutrition, vol. 140, no. 11, pp. 1956–1962, 2010.
View at: Publisher Site | Google Scholar
M. Kolf-Clauw, J. Castellote, B. Joly et al., “Development of a pig jejunal explant culture for studying the gastrointestinal toxicity of the mycotoxin deoxynivalenol: histopathological analysis,” Toxicology In Vitro, vol. 23, no. 8, pp. 1580–1584, 2009.
View at: Publisher Site | Google Scholar
K. Ghareeb, W. A. Awad, J. Bohm, and Q. Zebeli, “Impacts of the feed contaminant deoxynivalenol on the intestine of monogastric animals: poultry and swine,” Journal of Applied Toxicology, vol. 35, no. 4, pp. 327–337, 2015.
View at: Publisher Site | Google Scholar
A. J. Brosnahan and D. R. Brown, “Porcine IPEC-J2 intestinal epithelial cells in microbiological investigations,” Veterinary Microbiology, vol. 156, no. 3-4, pp. 229–237, 2012.
View at: Publisher Site | Google Scholar
R. Kang, R. Li, P. Dai, Z. Li, Y. Li, and C. Li, “Deoxynivalenol induced apoptosis and inflammation of IPEC-J2 cells by promoting ROS production,” Environmental Pollution, vol. 251, pp. 689–698, 2019.
View at: Publisher Site | Google Scholar
S. Adesso, G. Autore, A. Quaroni, A. Popolo, L. Severino, and S. Marzocco, “The food contaminants nivalenol and deoxynivalenol induce inflammation in intestinal epithelial cells by regulating reactive oxygen species release,” Nutrients, vol. 9, no. 12, p. 1343, 2017.
View at: Publisher Site | Google Scholar
M. L. Y. Wan, P. C. Turner, V. A. Co, M. F. Wang, K. M. A. Amiri, and H. El-Nezami, “Schisandrin A protects intestinal epithelial cells from deoxynivalenol-induced cytotoxicity, oxidative damage and inflammation,” Scientific Reports, vol. 9, no. 1, p. 19173, 2019.
View at: Publisher Site | Google Scholar
C. Peng, X. Ding, L. Zhu et al., “β-Conglycinin-induced intestinal porcine epithelial cell damage via the nuclear factor κB/mitogen-activated protein kinase signaling pathway,” Journal of Agricultural and Food Chemistry, vol. 67, no. 32, pp. 9009–9021, 2019.
View at: Publisher Site | Google Scholar
Y. Tang, J. Li, F. Li et al., “Autophagy protects intestinal epithelial cells against deoxynivalenol toxicity by alleviating oxidative stress via IKK signaling pathway,” Free Radical Biology & Medicine, vol. 89, pp. 944–951, 2015.
View at: Publisher Site | Google Scholar
L. Osaki and P. Gama, “MAPKs and signal transduction in the control of gastrointestinal epithelial cell proliferation and differentiation,” International Journal of Molecular Sciences, vol. 14, no. 5, pp. 10143–10161, 2013.
View at: Publisher Site | Google Scholar
M. S. Hayden and S. Ghosh, “Signaling to NF-B,” Genes & Development, vol. 18, no. 18, pp. 2195–2224, 2004.
View at: Publisher Site | Google Scholar
Y. Sun, W. Z. Liu, T. Liu, X. Feng, N. Yang, and H. F. Zhou, “Signaling pathway of MAPK/ERK in cell proliferation, differentiation, migration, senescence and apoptosis,” Journal of Receptor and Signal Transduction Research, vol. 35, no. 6, pp. 600–604, 2015.
View at: Publisher Site | Google Scholar
C. W. Huang, T. T. Lee, Y. C. Shih, and B. Yu, “Effects of dietary supplementation of Chinese medicinal herbs on polymorphonuclear neutrophil immune activity and small intestinal morphology in weanling pigs,” Journal of Animal Physiology and Animal Nutrition, vol. 96, no. 2, pp. 285–294, 2012.
View at: Publisher Site | Google Scholar
A. Abdallah, P. Zhang, Q. Zhong, and Z. Sun, “Application of traditional Chinese herbal medicine by-products as dietary feed supplements and antibiotic replacements in animal production,” Current Drug Metabolism, vol. 20, no. 1, pp. 54–64, 2019.
View at: Publisher Site | Google Scholar
L. Gu, X. Wang, Z. Liu et al., “A study of Semen Strychni-induced renal injury and herb–herb interaction of Radix Glycyrrhizae extract and/or Rhizoma Ligustici extract on the comparative toxicokinetics of strychnine and brucine in rats,” Food and Chemical Toxicology, vol. 68, pp. 226–233, 2014.
View at: Publisher Site | Google Scholar
G. L. Bentz, A. J. Lowrey, D. C. Horne et al., “Using glycyrrhizic acid to target sumoylation processes during Epstein-Barr virus latency,” PLoS One, vol. 14, no. 5, article e0217578, 2019.
View at: Publisher Site | Google Scholar
A. Afkhami-Poostchi, M. Mashreghi, M. Iranshahi, and M. M. Matin, “Use of a genetically engineered E. coli overexpressing β-glucuronidase accompanied by glycyrrhizic acid, a natural and anti-inflammatory agent, for directed treatment of colon carcinoma in a mouse model,” International Journal of Pharmaceutics, vol. 579, article 119159, 2020.
View at: Publisher Site | Google Scholar
K. Akutagawa, T. Fujita, K. Ouhara et al., “Glycyrrhizic acid suppresses inflammation and reduces the increased glucose levels induced by the combination of Porphyromonas gulae and ligature placement in diabetic model mice,” International Immunopharmacology, vol. 68, pp. 30–38, 2019.
View at: Publisher Site | Google Scholar
X. Huo, X. Meng, J. Zhang, and Y. Zhao, “Hepatoprotective effect of different combinations of 18α-and 18β-glycyrrhizic acid against CCl4-induced liver injury in rats,” Biomedicine & Pharmacotherapy, vol. 122, article 109354, 2020.
View at: Publisher Site | Google Scholar
M. M. Alfajaro, H.-J. Kim, J.-G. Park et al., “Anti-rotaviral effects of glycyrrhiza uralensis extract in piglets with rotavirus diarrhea,” Virology Journal, vol. 9, no. 1, p. 310, 2012.
View at: Publisher Site | Google Scholar
M. Martin, “Cutadapt removes adapter sequences from high-throughput sequencing reads,” EMBnet Journal, vol. 17, no. 1, 2011.
View at: Publisher Site | Google Scholar
D. Kim, B. Langmead, and S. L. Salzberg, “HISAT: a fast spliced aligner with low memory requirements,” Nature Methods, vol. 12, no. 4, pp. 357–360, 2015.
View at: Publisher Site | Google Scholar
S. Anders, P. T. Pyl, and W. Huber, “HTSeq--a Python framework to work with high-throughput sequencing data,” Bioinformatics, vol. 31, no. 2, pp. 166–169, 2015.
View at: Publisher Site | Google Scholar
D. Szklarczyk, A. L. Gable, D. Lyon et al., “STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets,” Nucleic Acids Research, vol. 47, no. D1, pp. D607–D613, 2019.
View at: Publisher Site | Google Scholar
Y. Liao, Z. Peng, L. Chen, A. K. Nussler, L. Liu, and W. Yang, “Deoxynivalenol, gut microbiota and immunotoxicity: a potential approach?” Food and Chemical Toxicology, vol. 112, pp. 342–354, 2018.
View at: Publisher Site | Google Scholar
I. Alassane-Kpembi, O. Puel, P. Pinton, A. M. Cossalter, T. C. Chou, and I. P. Oswald, “Co-exposure to low doses of the food contaminants deoxynivalenol and nivalenol has a synergistic inflammatory effect on intestinal explants,” Archives of Toxicology, vol. 91, no. 7, pp. 2677–2687, 2017.
View at: Publisher Site | Google Scholar
Z. Ren, H. Deng, Y. Deng et al., “Combined effects of deoxynivalenol and zearalenone on oxidative injury and apoptosis in porcine splenic lymphocytes in vitro,” Experimental and Toxicologic Pathology, vol. 69, no. 8, pp. 612–617, 2017.
View at: Publisher Site | Google Scholar
Z. Q. Zhang, S. B. Wang, R. G. Wang, W. Zhang, P. L. Wang, and X. O. Su, “Phosphoproteome analysis reveals the molecular mechanisms underlying deoxynivalenol-induced intestinal toxicity in IPEC-J2 cells,” Toxins (Basel), vol. 8, no. 10, p. 270, 2016.
View at: Publisher Site | Google Scholar
A. Springler, S. Hessenberger, G. Schatzmayr, and E. Mayer, “Early activation of MAPK p44/42 is partially involved in DON-induced disruption of the intestinal barrier function and tight junction network,” Toxins (Basel), vol. 8, no. 9, p. 264, 2016.
View at: Publisher Site | Google Scholar
R. Krishnaswamy, S. N. Devaraj, and V. V. Padma, “Lutein protects HT-29 cells against deoxynivalenol-induced oxidative stress and apoptosis: prevention of NF-κB nuclear localization and down regulation of NF-κB and cyclo-oxygenase – 2 expression,” Free Radical Biology & Medicine, vol. 49, no. 1, pp. 50–60, 2010.
View at: Publisher Site | Google Scholar
Q. Wu, X. Wang, E. Nepovimova et al., “Antioxidant agents against trichothecenes: new hints for oxidative stress treatment,” Oncotarget, vol. 8, no. 66, pp. 110708–110726, 2017.
View at: Publisher Site | Google Scholar
J. He, C. Xia, Y. He et al., “Proteomic responses to oxidative damage in meat from ducks exposed to heat stress,” Food Chemistry, vol. 295, pp. 129–137, 2019.
View at: Publisher Site | Google Scholar
X. Zhang, L. Jiang, C. Geng, J. Cao, and L. Zhong, “The role of oxidative stress in deoxynivalenol-induced DNA damage in HepG2 cells,” Toxicon, vol. 54, no. 4, pp. 513–518, 2009.
View at: Publisher Site | Google Scholar
J. Xie, H. Ye, M. Du, Q. Yu, Y. Chen, and M. Shen, “Mung bean protein hydrolysates protect mouse liver cell line nctc-1469 cell from hydrogen peroxide-induced cell injury,” Foods, vol. 9, no. 1, p. 14, 2020.
View at: Publisher Site | Google Scholar
M. Yang, H.-B. Lin, S. Gong et al., “Effect of astragalus polysaccharides on expression of TNF-α, IL-1β and NFATc4 in a rat model of experimental colitis,” Cytokine, vol. 70, no. 2, pp. 81–86, 2014.
View at: Publisher Site | Google Scholar
A. Rostom, C. Dubé, G. Lewin et al., “Nonsteroidal anti-inflammatory drugs and cyclooxygenase-2 inhibitors for primary prevention of colorectal cancer: a systematic review prepared for the U.S. Preventive Services Task Force,” Annals of Internal Medicine, vol. 146, no. 5, pp. 376–389, 2007.
View at: Publisher Site | Google Scholar
B. B. Aggarwal, Y. Takada, S. Shishodia et al., “Nuclear transcription factor NF-kappa B: role in biology and medicine,” Indian Journal of Experimental Biology, vol. 42, no. 4, pp. 341–353, 2004.
View at: Google Scholar
S. Elmore, “Apoptosis: a review of programmed cell death,” Toxicologic Pathology, vol. 35, no. 4, pp. 495–516, 2016.
View at: Publisher Site | Google Scholar
L. Portt, G. Norman, C. Clapp, M. Greenwood, and M. T. Greenwood, “Anti-apoptosis and cell survival: a review,” Biochimica et Biophysica Acta, vol. 1813, no. 1, pp. 238–259, 2011.
View at: Publisher Site | Google Scholar
C. Andrews, M. H. McLean, and S. K. Durum, “Cytokine tuning of intestinal epithelial function,” Frontiers in Immunology, vol. 9, 2018.
View at: Publisher Site | Google Scholar
C. E. Hughes and R. J. B. Nibbs, “A guide to chemokines and their receptors,” The FEBS Journal, vol. 285, no. 16, pp. 2944–2971, 2018.
View at: Publisher Site | Google Scholar
J. W. Griffith, C. L. Sokol, and A. D. Luster, “Chemokines and chemokine receptors: positioning cells for host defense and immunity,” Annual Review of Immunology, vol. 32, no. 1, pp. 659–702, 2014.
View at: Publisher Site | Google Scholar
R. S. Bystry, V. Aluvihare, K. A. Welch, M. Kallikourdis, and A. G. Betz, “B cells and professional APCs recruit regulatory T cells via CCL4,” Nature Immunology, vol. 2, no. 12, pp. 1126–1132, 2001.
View at: Publisher Site | Google Scholar
R. E. Marques, R. Guabiraba, R. C. Russo, and M. M. Teixeira, “Targeting CCL5 in inflammation,” Expert Opinion on Therapeutic Targets, vol. 17, no. 12, pp. 1439–1460, 2013.
View at: Publisher Site | Google Scholar
V. Appay and S. L. Rowland-Jones, “RANTES: a versatile and controversial chemokine,” Trends in Immunology, vol. 22, no. 2, pp. 83–87, 2001.
View at: Publisher Site | Google Scholar
Q. Liu, A. Li, S. Yu et al., “DACH1 antagonizes CXCL8 to repress tumorigenesis of lung adenocarcinoma and improve prognosis,” Journal of Hematology & Oncology, vol. 11, no. 1, p. 53, 2018.
View at: Publisher Site | Google Scholar
H. Zhang, X. Deng, C. Zhou, W. Wu, and H. Zhang, “Deoxynivalenol induces inflammation in IPEC-J2 cells by activating P38 Mapk and Erk1/2,” Toxins (Basel), vol. 12, no. 3, p. 180, 2020.
View at: Publisher Site | Google Scholar
S.-C. Lin, P.-Y. Chu, W.-T. Liao et al., “Glycyrrhizic acid induces human MDA-MB-231 breast cancer cell death and autophagy via the ROS-mitochondrial pathway,” Oncology Reports, vol. 39, no. 2, pp. 703–710, 2017.
View at: Publisher Site | Google Scholar
D. Girard, M. E. Paquet, R. Paquin, and A. D. Beaulieu, “Differential effects of interleukin-15 (IL-15) and IL-2 on human neutrophils: modulation of phagocytosis, cytoskeleton rearrangement, gene expression, and apoptosis by IL-15,” Blood, vol. 88, no. 8, pp. 3176–3184, 1996.
View at: Publisher Site | Google Scholar
M. Pelletier, C. Ratthe, and D. Girard, “Mechanisms involved in interleukin-15-induced suppression of human neutrophil apoptosis: role of the anti-apoptotic Mcl-1 protein and several kinases including Janus kinase-2, p38 mitogen-activated protein kinase and extracellular signal-regulated kinases-1/2,” FEBS Letters, vol. 532, no. 1-2, pp. 164–170, 2002.
View at: Publisher Site | Google Scholar
J. V. Dzimianski, F. E. M. Scholte, E. Bergeron, and S. D. Pegan, “ISG15: it’s complicated,” Journal of Molecular Biology, vol. 431, no. 21, pp. 4203–4216, 2019.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Xiaoxiang Xu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

693

Downloads

864

Citations