Next Article in Journal
Control of Aflatoxigenic Molds by Antagonistic Microorganisms: Inhibitory Behaviors, Bioactive Compounds, Related Mechanisms, and Influencing Factors
Next Article in Special Issue
T-2 Toxin Induces Oxidative Stress, Apoptosis and Cytoprotective Autophagy in Chicken Hepatocytes
Previous Article in Journal
Inflammasome Activation Induced by a Snake Venom Lys49-Phospholipase A2 Homologue
Previous Article in Special Issue
Reduced Toxicity of Trichothecenes, Isotrichodermol, and Deoxynivalenol, by Transgenic Expression of the Tri101 3-O-Acetyltransferase Gene in Cultured Mammalian FM3A Cells
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxins
Volume 12
Issue 1
10.3390/toxins12010023
toxins-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptome Analysis of Ochratoxin A-Induced Apoptosis in Differentiated Caco-2 Cells

by
Xue Yang
1,2,3,4,†,
Yanan Gao
1,2,3,4,†,
Qiaoyan Yan
1,2,3,4,†,
Xiaoyu Bao
1,2,3,4,
Shengguo Zhao
1,2,3,4,
Jiaqi Wang
1,2,3,4 and
Nan Zheng
1,2,3,4,*
1
Key Laboratory of Quality & Safety Control for Milk and Dairy Products of Ministry of Agriculture and Rural Affairs, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China
2
Laboratory of Quality and Safety Risk Assessment for Dairy Products of Ministry of Agriculture and Rural Affairs, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China
3
Milk and Dairy Product Inspection Center of Ministry of Agriculture and Rural Affairs, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China
4
State Key Laboratory of Animal Nutrition, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Toxins 2020, 12(1), 23; https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12010023
Submission received: 4 November 2019 / Revised: 21 December 2019 / Accepted: 24 December 2019 / Published: 31 December 2019
(This article belongs to the Special Issue Toxicological Effects of Mycotoxins on Target Cells)
Browse Figures
Review Reports Versions Notes

Abstract

:
Ochratoxin A (OTA), an important mycotoxin that occurs in food and animal feed, has aroused widespread concern in recent years. Previous studies have indicated that OTA causes nephrotoxicity, hepatotoxicity, genotoxicity, immunotoxicity, cytotoxicity, and neurotoxicity. The intestinal toxicity of OTA has gradually become a focus of research, but the mechanisms underlying this toxicity have not been described. Here, differentiated Caco-2 cells were incubated for 48 h with different concentrations of OTA and transcriptome analysis was used to estimate damage to the intestinal barrier. Gene expression profiling was used to compare the characteristics of differentially expressed genes (DEGs). There were altogether 10,090 DEGs, mainly clustered into two downregulation patterns. The Search Tool for Retrieval of Interacting Genes (STRING), which was used to analyze the protein–protein interaction network, indicated that 24 key enzymes were mostly responsible for regulating cell apoptosis. Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) analysis was used to validate eight genes, three of which were key genes (CASP3, CDC25B, and EGR1). The results indicated that OTA dose-dependently induces apoptosis in differentiated Caco-2 cells. Transcriptome analysis showed that the impairment of intestinal function caused by OTA might be partly attributed to apoptosis, which is probably associated with downregulation of murine double minute 2 (MDM2) expression and upregulation of Noxa and caspase 3 (CASP3) expression. This study has highlighted the intestinal toxicity of OTA and provided a genome-wide view of biological responses, which provides a theoretical basis for enterotoxicity and should be useful in establishing a maximum residue limit for OTA.
Keywords:
ochratoxin A; differentiated Caco-2 cells; cell apoptosis; transcriptome analysis
Key Contribution: Transcription analysis indicated that OTA-induced intestinal toxicity may be induced by apoptosis through the regulation of MDM2 and CASP3 which provide a theoretical basis for toxicological evaluation of OTA.

1. Introduction

Ochratoxin A (OTA) is a fungal secondary metabolite produced by certain Penicillium and Aspergillus species, including Penicillium verrucosum, Aspergillus ochraceus, and Aspergillus niger [1]. OTA was first isolated from A. ochraceus in 1965 and was found to contaminate the food chain worldwide [2]. It is widely found in various grains and vegetables [3,4,5], as well as in food products of animal origin, such as meat, eggs, and milk [6,7,8,9]. Milk, which has high bioavailability and is an abundant source of nutrients, is widely recognized to be an important component of the human diet. As the consumption of milk has increased over recent years, the mycotoxins found in milk have received increasing attention. A provisional tolerable weekly intake of 100 ng/kg.bw/week has been established for OTA by the Joint FAO/WHO Expert Committee on Food Additives (JECFA) [10], JECFA, 2001), although a maximum residue limit (MRL) for OTA in milk has not been agreed upon internationally. A study in Italy detected OTA concentrations of 70–110 ng/L in organic milk [11]. In China, Huang et al. [12] measured levels of OTA in raw, powdered, and liquid cow’s milk, and found mean concentrations of 56.7, 27.0, and 26.8 ng/kg, respectively. In Sudan, the level of OTA in a contaminated milk sample was as high as 2730 ng/L [13]. OTA has been classified in Group 2B (possible carcinogens in humans) by The International Agency for Research on Cancer, because of evidence of carcinogenicity in animals, but not in humans [14].
The gastrointestinal tract (GIT) is essential for human health and provides a barrier between the external environment and the tightly regulated internal environment [15]. The GIT can be exposed to numerous contaminated foods and high doses of some mycotoxins [16]. Early studies on OTA focused mainly on the diversity of toxic effects in different animal species [17,18,19]. Recent studies, however, have reported the toxic effect of OTA on the intestine [20,21]. OTA-induced intestinal damage has been reported in both animals and in vitro intestinal models [22,23]. It has been shown to damage the intestinal epithelium in chickens and rats [24,25,26] and also shows toxicity in intestinal epithelial cells, including a porcine intestinal cell line (IPEC-J2) and human intestinal epithelial lines (HT-29-D4 cells and Caco-2 cells) [20,27,28]. Previously published studies have shown that cell apoptosis is one of the ways by which OTA exerts intestinal toxicity [20,26,29]. Wang et al. [20] suggested that apoptosis induced by OTA may play a major role in the intestinal toxicity of this mycotoxin, and Bouaziz et al. [30] also suggested that OTA causes toxicity through apoptosis. OTA has been shown to induce apoptosis in different cell lines [31,32,33,34,35], which may be one of the main cellular mechanisms underlying the toxic effects. However, the molecular mechanisms responsible for cell apoptosis, which leads to intestinal toxicity, are still inadequately understood. It is, therefore, important to investigate the mechanism of apoptosis of intestinal epithelial cells following exposure to OTA. In the present study, we used differentiated Caco-2 cells as they can form polarized apical/mucosal and basolateral/serosal membranes that are similar to those formed by epithelial cells in the small intestine [36]. Moreover, it has been acknowledged by the Food and Drug Administration that differentiated Caco-2 cells are a suitable model for evaluating the impact of toxins on intestinal barrier function [37,38].
With large-scale transcription approaches, a comprehensive overview, provided by high throughput data, can easily reveal biological pathways and processes that have not been found before [39]. Using whole-genome transcriptome profiling, RNA sequencing (RNA-seq), an unbiased sequencing tool, has been used to detect changes of gene expression in tissue samples or cells [40,41]. The aim of this study was to investigate the mechanism of OTA-induced apoptosis in differentiated Caco-2 cells and to use RNA-seq technology to evaluate changes in gene expression and profile to clarify the mechanism of OTA-induced apoptosis in differentiated Caco-2 cells.

2. Results

2.1. OTA Induces Apoptosis in Differentiated Caco-2 Cells in a Dose-Dependent Manner

After treatment with OTA (0.0005, 0.005, and 4 μg/mL) for 48 h, Annexin V-PI dual staining and flow cytometry were used to detect apoptosis of differentiated Caco-2 cells (Figure 1a). As depicted in the histogram (Figure 1b), the percentage of apoptotic cells was increased by treatment for 48 h with 0.0005 μg/mL (2.9% ± 0.67%), 0.005 μg/mL (5.23% ± 0.21%), and 4 μg/mL (20.9% ± 2.49%) OTA, compared with the control (2.3% ± 0.08%). The percentage of cell apoptosis was significantly increased in differentiated Caco-2 cells (p < 0.05) when the concentration of OTA was 4 μg/mL (Figure 1b). The number of living cells in the Q3 regions of the flow cytometry plots also decreased as the concentration of OTA increased (Figure 1a).

2.2. Effect of OTA on Gene Expression Patterns

There were 18,654 unigene annotations in the control group (56.56% of all the 32,938 reference unigene sequences) and 18,475 (56.56%), 18,898 (57.30%), and 19,853 (55.61%) unigene annotations in the 0.0005, 0.005, and 4 μg/mL OTA groups., respectively. Using a p-value < 0.05 and a 2-fold change (FC) as the conditions for discrimination, we identified 503, 2139, and 9402 differentially expressed genes (DEGs) when the cells were incubated for 48 h with 0.0005, 0.005, and 4 μg/mL OTA, respectively. At these three concentrations of OTA, 35, 258, and 1437 genes were upregulated and 468, 1881, and 7965 genes were downregulated, respectively (Figure 2). At all concentrations of OTA, downregulation of genes was the main trend (93%, 88%, and 85% of the total DEGs, respectively). The number of DEGs in differentiated Caco-2 cells increased with increasing OTA concentration. Using a Venn diagram, we found that 3 upregulated DEGs and 169 downregulated DEGs were commonly modulated by all three concentrations of OTA (Figure 3a,b). Compared with the control, there were 257 overlapping DEGs between the 0.0005 and 0.005 μg/mL OTA groups and 1532 overlapping DEGs between the 0.005 and 4 μg/mL groups (Supplementary Figure S1), showing that the common DEGs were increased in a dose-dependent manner. To assess the expression patterns of mRNAs at different concentrations of OTA, we used heatmaps to analyze the overall transcriptome differences (Figure 3c). The same differently expressed transcripts were present in three separate runs of the control group (CTL1, CTL2, CTL3), the 0.0005 μg/mL treatment group (0.0005-1,2,3), the 0.005 μg/mL treatment group (0.005-1,2,3), and the 4 μg/mL treatment group (4-1,2,3) (Figure 3c). The heatmap (Figure 3c) shows accurate repeatability and high reliability.

2.3. Gene Ontology (GO) Annotation and KEGG Enrichment Analysis of DEGs

Using the FC and p-value thresholds described in the Methods section above, we obtained 10,090 DEGs (compared with the control, the sum of DEGs of the OTA treatment at different concentrations), which could be clustered into eight profiles using the short time-series expression miner to obtain dynamic expression patterns of the DEGs. Among these eight profiles, two classic downregulated profiles (profiles 0 and 3) were significantly enriched (Supplementary Figure S2). A total of 7199 DEGs were mainly clustered into these two downregulated profiles, which contained 1883 and 5226 genes, respectively (Supplementary Figure S3). To discover the functions of the DEGs and the associated biological processes altered by OTA treatment of differentiated Caco-2 cells, we carried out a gene ontology (GO) (http://www.geneontology.org/) enrichment analysis, in which the DEGs were divided into three independent GO categories. Single DEGs could be annotated to more than two GO terms and the most enriched GO terms are displayed in Figure S3. The three main categories of the GO classification were assessed for enriched genes.
All the DEGs were mapped in the Kyoto Encyclopedia of Genes and Genomes (KEGG, https://www.kegg.jp) database to detect the response pathways altered by treatment with OTA. Of the DEGs that could be annotated to the KEGG pathway, most were associated with metabolism and signal transduction pathways. In profile 0, the enriched pathways were pyrimidine metabolism, the TNF signaling pathway, ribosome biogenesis in eukaryotes, and the Wnt signaling pathway. In profile 3, the enriched pathways were endocytosis, cell cycle, ubiquitin-mediated proteolysis, pyrimidine metabolism, and AMPK signaling (Figure 4).

2.4. Key Pathways and DEGs Related to Cell Apoptosis

Using the GO annotation and KEGG analysis, we can select the key pathways that are associated with the toxic effects of OTA in differentiated Caco-2 cells. DEGs from the key significantly enriched pathways, particularly those shared between the different concentrations of OTA, may be thought of as the key gene expression regulators that respond to treatment with OTA. These 10 key pathways, together with some key genes that participate in or regulate the cell cycle and cell apoptosis, are shown in Table 1.
As shown in Figure 5, the 50 DEGs that occurred frequently in these key pathways, or participated in individual key pathways, were combined and employed to construct a protein–protein interaction (PPI) network using STRING. Cytoscape 3.1 was then used for further filtering (Edge score > 0.8). The 10 key pathways and 24 key DEGs associated with cell apoptosis were identified mainly from those shared between the different concentrations of OTA in the GO and KEGG analysis, and these key pathways were largely associated with the regulation of cell apoptosis and the cell cycle. According to the degree of connectivity of each node, we selected the following 24 key enzymes: murine double minute 2 (MDM2), v-akt murine thymoma viral oncogene homolog 1 (AKT1), tumor protein p53 (TP53), caspase 3 (CASP3), caspase 9 (CASP9), hras proto-oncogene gtpase (HRAS), mechanistic target of rapamycin kinase (MTOR), epidermal growth factor receptor (EGFR), bcl-2-like protein 1 (BCL2L1), tumor necrosis factor (TNF), atm serine/threonine kinase (ATM), cytochrome C somatic (CYCs), mitogen-activated protein kinase 1 and 14 (MAPK1, MAPK14), ribosomal protein s6 kinase B1 (RPS6KB1), x-linked inhibitor of apoptosis (XIAP), nuclear factor kappa B subunit 1 (NFKB1), cyclin B1 (CCNB1), E1A binding protein p300 (EP300), inhibitor of nuclear factor kappa B kinase subunit beta (IKBKB), TSC complex subunit 2 (TSC2), phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha (PIK3CA), cyclin dependent kinase 1 (CDK1), and forkhead box O1 (FOXO1).

2.5. Validation of RNA-Seq Results by qRT-PCR

Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) analysis was performed to validate the high throughput data. We determined the expression of the following eight randomly selected genes: caspase 3 (CASP3), C–X–C motif chemokine ligand 2 (CXCL2), cell division cycle 25B (CDC25B), early growth response 1 (EGR1), FRY like transcription coactivator (FRYL), H2B histone family members (H2BFS), sedoheptulokinase (SHPK), and transcription factor EC (TFEC). Expression levels of CDC25B, SHPK, ERG1, FRYL, and TFEC significantly decreased after OTA treatment. In contrast, expression of CASP3, CXCL2, and H2BFS significantly increased (Figure 6a). These expression patterns were in good agreement with the data obtained by RNA-Seq, confirming the reliability of the two methods. There was a strong correlation (r = 0.803) between RNA-seq and qRT-PCR (Figure 6b). The protein levels of two randomly selected genes, MDM2 and CASP3, were measured by Western blotting. MDM2 levels significantly decreased and CASP3 levels increased after treatment with OTA (Figure 7).

3. Discussion

OTA has deleterious effects on humans and animals, resulting in worldwide illness and economic losses. Exposure to OTA occurs from ingestion of contaminated food and feed [42]. The small intestine is the main site of OTA absorption, with the largest absorption in the proximal jejunum [43]. Because of their location and function, intestinal epithelial cells are a potential target for the toxic effects of OTA. In the present study, flow cytometry showed that the number of apoptotic cells increased with increasing OTA concentration (Figure 1), which is consistent with previous studies demonstrating that OTA can induce apoptosis of intestinal cells [20,44,45]. Apoptosis, which is a form of programmed cell death, may disrupt the integrity of the intestinal barrier [46,47]. OTA induces apoptosis by modulating BcL-2 family members [30,48], and activation of MEK/ERK1-2 signaling has been shown to be crucial for OTA-induced apoptosis in HK-2 cells [49]. Although OTA-induced apoptosis has been studied previously, there has been a tendency to focus mainly on a particular pathway [20,49,50]. In this study, we chose instead to search for related pathways from a holistic perspective because transcriptome data can quickly and economically supply accurate genome information, particularly in identifying the effects of biological pathways and processes [51,52]. Elena et al. (2016) used transcriptomic analysis to reveal the different toxicity mechanisms of OTA and citrinin [53]. Hence, we used a transcriptomic study to provide a genome-wide biological view of changes that occurred after treatment of differentiated Caco-2 cells with different concentrations of OTA.
For the GO annotation, the main GO terms of the DEGs were cell processes and metabolic processes. Another microarray study demonstrated that OTA decreased mRNA levels of genes involved in liver cell metabolism [54]. KEGG enrichment analysis showed that there are key pathways, such as the cell cycle, the MAPK signaling pathway, the TNF signaling pathway, and the p53 signaling pathway, that are highly associated with cell apoptosis (Figure 4). These results indicate that OTA-induced cell apoptosis probably involves changes to the enriched pathways, such as the cell cycle, and other important signaling pathways. From the results obtained in this study, we speculate that cell apoptosis happens in an orderly manner and can be regulated by many mechanisms. Liang et al. [55] pointed out that OTA suppresses the cell cycle, particularly DNA replication, leading to cell cycle arrest and apoptosis. OTA has been shown to regulate apoptosis, the cell cycle, and cell fate following DNA damage by stimulating MAPK signaling pathways, such as p38 MAPK, ERK1/2, and JNK [56,57,58]. Apoptosis of liver and/or kidney cells has also been shown to be influenced by activation of certain signal transduction pathways, such as MAPK, ERK, p38, and JNK [56,59,60]. Another study obtained data that were of great importance in characterizing the modulatory effect of the TNF signaling pathway on apoptotic signaling [61]. Activation of the p53 pathway has also been shown to play a key role in OTA-induced apoptosis of human and monkey kidney cells [33]. These results are in agreement with reported research demonstrating that OTA triggers a p53-dependent apoptotic pathway in human hepatoma cells [62]. Using knockout mice, Hibi et al. showed that the p53 pathway plays a key role in OTA-induced genotoxicity [63]. Because previous research has highlighted the role of p53 signaling in OTA-induced apoptosis, in this study, we concentrated on the p53 pathway and analyzed its specific regulatory mechanism.
The p53 signaling pathway, which is amplified by different concentrations of OTA, plays a key role in regulating DNA repair mechanisms, oxidative stress, cell cycle arrest, and cell apoptosis [64,65]. In the renal outer medulla, OTA causes genotoxicity by deregulating molecular functions such as DNA double-strand break repair, cell cycle arrest in response to DNA damage, and p53-associated factors [66]. Alterations in oxidative stress and calcium homeostasis have also been associated with OTA-induced toxicity [67]. The proto-oncogene, MDM2, a nuclear protein that is the main cellular antagonist of p53, can regulate the potentially lethal activities of p53 by tightly combining with a p53 tumor suppressor protein and negatively regulating its stability and transcriptional activity [68]. Furthermore, MDM2 can not only inhibit the transactivation of p53 activity, but also stabilize p53 and reduce p53-mediated apoptosis [69]. The p53-MDM2 feedback loop may be important in regulating cell apoptosis. Previous studies have demonstrated that apoptotic pathways are regulated by proteins such as the tumor suppressor p53 and inhibitor of apoptosis (IAP) proteins, which are highly regulated by MDM2 through an autoregulatory feedback loop [70,71]. OTA may induce apoptosis in differentiated Caco-2 cells through the feedback loop interaction of MDM2 and the p53 signaling pathway. A study by Gu et al. [72] showed that gambogic acid induced apoptosis of wild-type p53-expressing cancer cells through downregulation of MDM2. Another study found that violacein decreased the expression of MDM2 and caused apoptosis of human breast cancer cells by activating PARP, CDKN1A, TNF-α, and p53 cleavage [73]. Our transcriptome analysis showed that OTA induced apoptosis in differentiated Caco-2 cells by activating key downstream genes in the p53 signaling pathway. In these key pathways, some DEGs, such as CASP3, AKT1, MAPK1, HRAS, PIK3CA, and MDM2, were highly involved (Figure 5), which appeared in the group 4 µg/mL specially instead of common DEGs (Figure 3). Several genes (CASP3, CDC25B, CXCL, and EGR1) were randomly selected from these key DEGs and common DEGs (FRYL, TFEC, SHPK, and H2BFS) were used for gene verification by qRT-PCR analysis. These key genes, and other DEGs, showed similar results in the RNA-seq analysis (Figure 6), and further study is needed to elucidate the relationships between the genes and regulation of cell apoptosis. In contrast to our own study, Kuroda et al. found that expression of the CDK1 gene was increased by treatment with OTA [74]. We believe that the difference between the studies is because OTA increases DNA damage in gpt delta rats and CDK1 could accelerate the resection of DNA broken ends during the homologous recombination process [74,75,76]. Through analysis of the enrichment pathway, we found that p53 is a target signaling pathway and speculate that OTA probably induces cell apoptosis by downregulating MDM2 expression. Western blotting confirmed that OTA-induced cell apoptosis in differentiated Caco-2 cells is linked to perturbation of MDM2 (Figure 7).
In addition to the downregulation of MDM2, there are several important genes in the p53 signaling pathway that are upregulated in the apoptosis of differentiated Caco-2 cells. These are Noxa (phorbol-12-myristate-13-acetate-induced Protein1, PMAIP1), tumor protein p53 regulated apoptosis-inducing protein 1 (P53AIP1), and caspase 3 (cysteinyl aspartate proteases, CASP3). P53 transactivates noxa, bax, puma, and other apoptotic response genes, and their products trigger mitochondrial apoptosis pathways [50]. Noxa is a member of the pro-apoptotic B-cell lymphoma 2 (BCL-2) family and can not only promote activation of caspases and apoptosis, but also promote changes in the mitochondrial membrane and outflow of apoptogenic proteins from the mitochondria [77,78]. Functional studies have shown that P53AIP is a pro-apoptosis molecule, and it is believed to play a vital role in mediating p53-dependent apoptosis. P53AIP1 can lead to CASP3 activation and downstream activities by decreasing mitochondrial membrane potential and inducing release of cytochrome c. CASP3 plays a vital role in many events involved in cell apoptosis, including activation of the caspase cascade and the execution of apoptosis [79,80,81]. OTA-induced downregulation of MDM2 could regulate the downstream upregulation of Noxa and P53AIP1, and eventually activate CASP3 in the apoptosis of Caco-2 cells.
To conclude, we have shown that OTA-induced apoptosis in differentiated Caco-2 cells may involve downregulation of MDM2 and upregulation of CASP3, through activation of the p53-mediated cell apoptosis signaling pathway. A more exhaustive study is, however, needed to fully elucidate all of the underlying mechanisms. We have also demonstrated the intestinal toxicity of OTA and provided a genome-wide view of biological responses, which provides a theoretical basis for enterotoxicity and should be useful in establishing an MRL for OTA.

4. Materials and Methods

4.1. Chemicals and Reagents

OTA powder (C20H18ClNO6; molecular weight, 403) was purchased from Pribolab (Qingdao, China). Human colon adenocarcinoma Caco-2 cells (passage number 18) were acquired from the American Type Culture Collection (Manassas, VA, USA). Fetal bovine serum (FBS) and Dulbecco’s modified Eagle’s medium (DMEM) were obtained from Gibco (Carlsbad, CA, USA). Nonessential amino acids (NEAA), trypsin (2.5%), antibiotics (100 units/mL penicillin, 100 µg/mL streptomycin), phosphate-buffered saline (PBS), Western and IP Cell Lysis Buffer, and an Annexin V-FITC apoptosis detection kit were supplied by Beyotime Biotechnology (Shanghai, China). A stock solution of OTA (1000 μg/mL) was obtained by dissolving OTA in methanol, and stored at −20 °C for later use. Rabbit anti-β-actin (58169S) and rabbit anti-MDM2 antibodies (86934S) were obtained from Cell Signaling Technology (Boston, MA, USA), and goat anti-rabbit IgG conjugated to horseradish peroxidase was obtained from Bioss (Beijing, China).

4.2. Cell Culture and Treatments

Differentiated Caco-2 cells were seeded into six-well Transwell chambers (Corning, NY, USA) at a density of 1 × 105 cells per well in DMEM containing 4.5 g/L glucose, 1% NEAA, 10% FBS, and antibiotics. The cells were incubated at 37 °C in a humidified atmosphere containing 5% CO2. A polarized epithelial monolayer of differentiated Caco-2 cells was formed by replacing the medium every other day for 21 days. Because an MRL for OTA in milk has not yet been established and previous studies have indicated that OTA and aflatoxin M1 (AFM1) have the same cytotoxicity in human intestinal Caco-2 cells [82], we chose the concentration of 0.0005 μg/mL for OTA, based on the Chinese (0.5 μg/kg) limits for AFM1 in milk. We established a concentration gradient when measuring cell viability and transepithelial electrical resistance [83] and chose the 4 μg/mL concentration to investigate the effect of an upper concentration limit on intestinal injury, without affecting cell survival. After 48 h of treatment with OTA (0.0005, 0.005, and 4 μg/mL), serum-free medium containing the same concentration of methanol was added to the control group. Differentiated Caco-2 cells were collected after 48 h for subsequent cell apoptosis studies, transcriptomics analysis, and qRT-PCR.

4.3. Cell Apoptosis Assay by Annexin V-FITC/PI FACS

Following the instructions of the Annexin V-FITC apoptosis detection kit, the differentiated Caco-2 cells were washed twice with cold PBS after treatment with different concentrations of OTA for 48 h. Cell samples were trypsinized, transferred to centrifuge tubes, and centrifuged at 1000× g for 5 min. The supernatant was discarded and the cells were resuspended in PBS (1 mL). The cell suspension was transferred to a 1.5 mL centrifuge tube and centrifuged again at 1000× g for 5 min. The supernatant was discarded and the cells were gently resuspended by adding Annexin V-FITC (5 μL), propidium iodide (10 μL), and Annexin V-FITC binding solution (195 μL). The cells were gently vortexed and incubated for 10–20 min at room temperature. An FC 500 MCL flow cytometer (Becton Dickinson, Mountain View, CA, USA) was used to analyze the cell samples within 1 h.

4.4. RNA Extraction, Library Construction, and Transcriptome Sequencing

4.4.1. RNA Extraction

After treatment with OTA (0.0005, 0.005, and 4 μg/mL) for 48 h, Trizol (Invitrogen, Camarillo, CA, USA) was used to extract total RNA from the differentiated Caco-2 cells in accordance with the manufacturer’s instructions. Remaining DNA was then removed by treatment with RNase-free DNase I (Takara Bio, Kusatsu, Shiga, Japan) for 30 min at 37 °C. Triple replicates of each treatment of differentiated Caco-2 cells were combined into a sample to reduce sample variability. After purification using an RNeasy Mini Kit (Qiagen, Dusseldorf, Germany), a NanoDrop 2000 spectrophotometer (Thermo Scientific, Wilmington, DE, USA) was used to assess RNA quality and RNase-free agarose gel electrophoresis was used to assess quantity.

4.4.2. Construction of cDNA Library and Illumina Sequencing

RNA-seq was performed on RNA samples from differentiated Caco-2 cells treated with each of the three concentrations of OTA described above. High-quality RNA from each sample was enriched using magnetic Oligo (dT) beads and used to construct and sequence a cDNA library. Suitable cDNA fragments were selected as templates for PCR amplification using index primers and NEB universal PCR primers. When constructed, the cDNA library was sequenced using a HiSeqTM 2500 RNA sequencing system (Illumina, San Diego, CA, USA).

4.4.3. De Novo Assembly and Quantification of Gene Abundance

Before data analysis, the quality of the original data was controlled and noise was reduced by data filtering. In order to obtain high-quality clean reads for subsequent information analysis, the clean reads were filtered more rigorously to remove reads containing adapters or more than 10% of unknown nucleotides (N), as well as all reads with more than 50% low-quality sequence (q-value ≤ 20). The transcriptome assembly program Trinity was used for de novo transcriptome reconstruction. The sequencing error rate was used to evaluate the quality of sequenced reads, the saturation of the library, and the randomness of sequencing. Cuffmerge was used to combine the transcripts from different replicates of one group into a comprehensive set of transcripts for subsequent differential analysis and to filter the unique annotation files for manually introduced assembly errors. Annotation of unigenes was obtained using Blast and compared with data obtained using GO, KEGG, KOG, NR, COG, and Swissprot databases. Gene expression levels were then calculated using both the “raw counts” mode and “FPKM” [84].

4.4.4. Differentially Expressed Genes and their Dynamic Expression Profile

The edge R package (http://www.rproject.org/) was used to determine the DEGs between different treatment groups [62]. Firstly, we used the general filtering standard (FDR (False Discovery Rate) < 0.05 and |log2FC|>1) to identify significant DEGs. All DEGs were analyzed using the short time-series expression miner [85]. GO annotation was analyzed using Blast2GO software [86]. Blast software was used against the KEGG database to determine the KEGG pathway annotation [75]. GO/KEGG functional enrichment analysis was carried out for the genes in each trend, and the p-value was calculated by the hypothesis test. After the p-value was corrected by FDR, the GO term and path satisfying this condition were assigned a q-value ≤ 0.05 threshold.

4.5. Validation of RNA-Seq Results Using qRT-PCR

A Prime Script™ II 1st Strand cDNA Synthesis Kit and TB Green™ Premix Ex Taq™ II (Takara, Kusatsu, Shiga, Japan) were used to reverse transcribe cDNAs. The qRT-PCR conditions were as follows: 95 °C for 30 s; 40 cycles at 95 °C for 5 s; 60 °C for 30 s; and, finally, 72 °C for 20 s. The sequences of the primers, which were synthesized by a commercial company (Sangon Biotech Co., Ltd., Shanghai, China), are shown in Table 2. The data were analyzed and 7500 Software v.2.0.1 (Applied Biosystems, Foster City, CA, USA) was then used to determine the cycle threshold (Ct). The GAPDH gene (house-keeping gene) was selected as an internal control to normalize the expression data. The 2−ΔΔCt method was used to calculate the relative expression of genes, and the mean and standard deviation of three biologic replicates are shown as the results [87].

4.6. Western Blotting Assays

Western blotting was used to verify differentially expressed proteins. Caco-2 cells were cultured in transwell chambers, with or without OTA (0.0005, 0.005, and 4 μg/mL), for 48 h. First, the cell samples were lysed with Western and IP Cell Lysis Buffer, and then equal amounts of protein were subjected to sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE). The samples were next transferred to polyvinylidene fluoride membranes by dry rotation and the membranes were blocked with 5% skim milk powder dissolved in Tris-buffered saline for 2 h at room temperature. The samples were then incubated for 2 h with specific primary antibodies (1:1000 in TBS), washed three times with TBST (TBS containing 0.1% Tween 20), and incubated for 2 h with secondary antibodies. The bands were imaged using a Tanon-5200 Chemiluminescent Imaging System (Tanon Science & Technology Co., Ltd.) and band densities were analyzed using Image J 2 × software (Version 2.1.0, National Institutes of Health, Bethesda, MD, USA, 2006). Intensity values were normalized to human β-actin.

4.7. Statistical Analysis

Analysis of data was performed using SPSS® Statistics version 19. Analysis of variance (ANOVA) followed by Tukey’s multiple comparison was used to test statistical differences between the OTA treatment groups and the control group. An asterisk (*) indicates a p-value < 0.05, which was regarded as statistically significant.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6651/12/1/23/s1, Figure S1: Venn diagram depicting the DEGs common to different dose of OTA exposure, Figure S2: DEGs expression profiles in ochratoxin A (OTA) treated differentiated Caco-2 cells, Figure S3: DEGs expression profile in the OTA treated differentiated Caco-2 cells.

Author Contributions

Conceptualization, X.Y., N.Z., Y.G., Q.Y. and J.W.; Formal analysis, X.Y., Y.G., S.Z. and X.B.; Project administration, N.Z. and J.W.; Supervision, N.Z. and Q.Y.; Writing—original draft, X.Y.; Writing—review & editing, X.Y., N.Z. and Y.G. All authors have read and agreed to the published version of the manuscript.

Funding

The study was supported by The National Key Research and Development Program (2017YFD0500502), The Scientific Research Project for Major Achievements of The Agricultural Science and Technology Innovation Program (ASTIP) (No.CAAS-ZDXT2019004), and Modern Agro-Industry Technology Research System of the PR China (CARS-36).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Creppy, E.E. Update of survey, regulation and toxic effects of mycotoxins in Europe. Toxicol. Lett. 2002, 127, 19–28. [Google Scholar] [CrossRef]
  2. Van der Merwe, K.J.; Steyn, P.S.; Fourie, L.; Scott, D.B.; Theron, J.J. Ochratoxin A, a toxic metabolite produced by Aspergillus ochraceus Wilh. Nature 1965, 205, 1112–1113. [Google Scholar] [CrossRef] [PubMed]
  3. Duarte, S.C.; Pena, A.; Lino, C.M. A review on ochratoxin A occurrence and effects of processing of cereal and cereal derived food products. Food Microbiol. 2010, 27, 187–198. [Google Scholar] [CrossRef] [PubMed]
  4. Clark, H.A.; Snedeker, S.M. Ochratoxin A: Its cancer risk and potential for exposure. J. Toxicol. Environ. Health B Crit. Rev. 2006, 9, 265–296. [Google Scholar] [CrossRef] [PubMed]
  5. Kononenko, G.P.; Burkin, A.A.; Zotova, E.V.; Soboleva, N.A. Ochratoxin A: Study of grain contamination. Prikl. Biokhim. Mikrobiol. 2000, 36, 209–213. [Google Scholar] [PubMed]
  6. Duarte, S.C.; Lino, C.M.; Pena, A. Food safety implications of ochratoxin A in animal-derived food products. Vet. J. 2012, 192, 286–292. [Google Scholar] [CrossRef]
  7. Ringot, D.; Chango, A.; Schneider, Y.J.; Larondelle, Y. Toxicokinetics and toxicodynamics of ochratoxin A, an update. Chem. Biol. Interact. 2006, 159, 18–46. [Google Scholar] [CrossRef]
  8. Capriotti, A.L.; Caruso, G.; Cavaliere, C.; Foglia, P.; Samperi, R.; Laganà, A. Multiclass mycotoxin analysis in food, environmental and biological matrices with chromatography/mass spectrometry. Mass Spectrom. Rev. 2012, 31, 466–503. [Google Scholar] [CrossRef]
  9. Flores-Flores, M.E.; Lizarraga, E.; de Cerain, A.L.; Gonzalez-Penas, E. Presence of mycotoxins in animal milk: A review. Food Control. 2015, 53, 163–176. [Google Scholar] [CrossRef]
  10. Joint FAO/WHO Expert Committee on Food Additives. Safety Evaluation of Certain Mycotoxins in Food; World Health Organization: Geneva, Switzerland, 2001. [Google Scholar]
  11. Pattono, D.; Gallo, P.F.; Civera, T. Detection and quantification of Ochratoxin A in milk produced in organic farms. Food Chem. 2011, 127, 374–377. [Google Scholar] [CrossRef]
  12. Huang, L.C.; Zheng, N.; Zheng, B.Q.; Wen, F.; Cheng, J.B.; Han, R.W.; Xu, X.M.; Li, S.L.; Wang, J.Q. Simultaneous determination of aflatoxin M1, ochratoxin A, zearalenone and alpha-zearalenol in milk by UHPLC-MS/MS. Food Chem. 2014, 146, 242–249. [Google Scholar] [CrossRef] [PubMed]
  13. Elzupir, A.O.; Makawi, S.Z.A.; Elhussein, A.M. Determination of Aflatoxins and Ochratoxin a in Dairy Cattle Feed and Milk in Wad Medani, Sudan. J. Anim. Vet. Adv. 2009, 8, 2508–2511. [Google Scholar]
  14. Cancer IAFO. Some Naturally Occurring Substances: Food Items and Constituents, Heterocyclic Aromatic Amines and Mycotoxins; International Agency for Research on Cancer: Geneva, Switzerland, 1993. [Google Scholar]
  15. Gambacorta, L.; Pinton, P.; Avantaggiato, G.; Oswald, I.P.; Solfrizzo, M. Grape Pomace, an Agricultural Byproduct Reducing Mycotoxin Absorption: In Vivo Assessment in Pig Using Urinary Biomarkers. J. Agric. Food Chem. 2016, 64, 6762–6771. [Google Scholar] [CrossRef] [PubMed]
  16. Bouhet, S.; Oswald, I.P. The effects of mycotoxins, fungal food contaminants, on the intestinal epithelial cell-derived innate immune response. Vet. Immunol. Immunopathol. 2005, 108, 199–209. [Google Scholar] [CrossRef] [PubMed]
  17. Pfohl-Leszkowicz, A.; Manderville, R.A. Ochratoxin A: An overview on toxicity and carcinogenicity in animals and humans. Mol. Nutr. Food Res. 2007, 51, 61–99. [Google Scholar] [CrossRef]
  18. Elling, F.; Hald, B.; Jacobsen, C.; Krogh, P. Spontaneous toxic nephropathy in poultry associated with ochratoxin A. Acta Pathol. Microbiol. Scand. A 1975, 83, 739–741. [Google Scholar] [CrossRef]
  19. Gekle, M.; Silbernagl, S. Renal toxicodynamics of ochratoxin A: A pathophysiological approach. Kidney Blood Press. Res. 1996, 19, 225–235. [Google Scholar] [CrossRef]
  20. Wang, H.; Chen, Y.; Zhai, N.; Chen, X.; Gan, F.; Li, H.; Huang, K. Ochratoxin A-Induced Apoptosis of IPEC-J2 Cells through ROS-Mediated Mitochondrial Permeability Transition Pore Opening Pathway. J. Agric. Food Chem. 2017, 65, 10630–10637. [Google Scholar] [CrossRef]
  21. Gao, Y.; Li, S.; Wang, J.; Luo, C.; Zhao, S.; Zheng, N. Modulation of Intestinal Epithelial Permeability in Differentiated Caco-2 Cells Exposed to Aflatoxin M1 and Ochratoxin A Individually or Collectively. Toxins 2018, 10, 13. [Google Scholar] [CrossRef] [Green Version]
  22. Maresca, M.; Fantini, J. Some food-associated mycotoxins as potential risk factors in humans predisposed to chronic intestinal inflammatory diseases. Toxicon 2010, 56, 282–294. [Google Scholar] [CrossRef]
  23. Schrickx, J.; Lektarau, Y.; Fink-Gremmels, J. Ochratoxin A secretion by ATP-dependent membrane transporters in Caco-2 cells. Arch. Toxicol. 2006, 80, 243–249. [Google Scholar] [CrossRef] [PubMed]
  24. Stoev, S.D.; Stefanov, M.; Denev, S.; Radic, B.; Domijan, A.M.; Peraica, M. Experimental mycotoxicosis in chickens induced by ochratoxin A and penicillic acid and intervention with natural plant extracts. Vet. Res. Commun. 2004, 28, 727–746. [Google Scholar] [CrossRef] [PubMed]
  25. Dortant, P.M.; Peters-Volleberg, G.W.; Van Loveren, H.; Marquardt, R.R.; Speijers, G.J. Age-related differences in the toxicity of ochratoxin A in female rats. Food Chem. Toxicol. 2001, 39, 55–65. [Google Scholar] [CrossRef]
  26. Solcan, C.; Pavel, G.; Floristean, V.C.; Chiriac, I.S.; Şlencu, B.G.; Solcan, G. Effect of ochratoxin A on the intestinal mucosa and mucosa-associated lymphoid tissues in broiler chickens. Acta Vet. Hung. 2015, 63, 30–48. [Google Scholar] [CrossRef] [Green Version]
  27. Cano-Sancho, G.; González-Arias, C.A.; Ramos, A.J.; Sanchis, V.; Fernández-Cruz, M.L. Cytotoxicity of the mycotoxins deoxynivalenol and ochratoxin A on Caco-2 cell line in presence of resveratrol. Toxicol. In Vitro 2015, 29, 1639–1646. [Google Scholar] [CrossRef] [Green Version]
  28. Revajová, V.; Levkut, M.; Levkutová, M.; Bořutová, R.; Grešaková, L.; Košiková, B.; Leng, L. Effect of lignin supplementation of a diet contaminated with Fusarium mycotoxins on blood and intestinal lymphocyte subpopulations in chickens. Acta Vet. Hung. 2013, 61, 354–365. [Google Scholar] [CrossRef]
  29. Ranaldi, G.; Caprini, V.; Sambuy, Y.; Perozzi, G.; Murgia, C. Intracellular zinc stores protect the intestinal epithelium from Ochratoxin A toxicity. Toxicol. In Vitro 2009, 23, 1516–1521. [Google Scholar] [CrossRef]
  30. Bouaziz, C.; Sharaf el dein, O.; Martel, C.; El Golli, E.; Abid-Essefi, S.; Brenner, C.; Lemaire, C.; Bacha, H. Molecular events involved in ochratoxin A induced mitochondrial pathway of apoptosis, modulation by Bcl-2 family members. Environ. Toxicol. 2011, 26, 579–590. [Google Scholar] [CrossRef]
  31. Rached, E.; Pfeiffer, E.; Dekant, W.; Mally, A. Ochratoxin A: Apoptosis and aberrant exit from mitosis due to perturbation of microtubule dynamics? Toxicol. Sci. 2006, 92, 78–86. [Google Scholar] [CrossRef] [Green Version]
  32. Gekle, M.; Sauvant, C.; Schwerdt, G. Ochratoxin A at nanomolar concentrations: A signal modulator in renal cells. Mol. Nutr. Food Res. 2005, 49, 118–130. [Google Scholar] [CrossRef]
  33. Li, J.; Yin, S.; Dong, Y.; Fan, L.; Hu, H. P53 activation inhibits ochratoxin A-induced apoptosis in monkey and human kidney epithelial cells via suppression of JNK activation. Biochem. Biophys. Res. Commun. 2011, 411, 458–463. [Google Scholar] [CrossRef] [PubMed]
  34. Cui, J.; Xing, L.; Li, Z.; Wu, S.; Wang, J.; Liu, J.; Wang, J.; Yan, X.; Zhang, X. Ochratoxin A induces G2 phase arrest in human gastric epithelium GES-1 cells in vitro. Toxicol. Lett. 2010, 193, 152–158. [Google Scholar] [CrossRef] [PubMed]
  35. Chopra, M.; Link, P.; Michels, C.; Schrenk, D. Characterization of ochratoxin A-induced apoptosis in primary rat hepatocytes. Cell Biol. Toxicol. 2010, 26, 239–254. [Google Scholar] [CrossRef] [PubMed]
  36. Artursson, P.; Palm, K.; Luthman, K. Caco-2 monolayers in experimental and theoretical predictions of drug transport. Adv. Drug Deliv. Rev. 2012, 64, 280–289. [Google Scholar] [CrossRef]
  37. Akbari, P.; Braber, S.; Varasteh, S.; Alizadeh, A.; Garssen, J.; Fink-Gremmels, J. The intestinal barrier as an emerging target in the toxicological assessment of mycotoxins. Arch. Toxicol. 2017, 91, 1007–1029. [Google Scholar] [CrossRef] [Green Version]
  38. Sambuy, Y.; De Angelis, I.; Ranaldi, G.; Scarino, M.L.; Stammati, A.; Zucco, F. The Caco-2 cell line as a model of the intestinal barrier: Influence of cell and culture-related factors on Caco-2 cell functional characteristics. Cell Biol. Toxicol. 2005, 21, 1–26. [Google Scholar] [CrossRef]
  39. Derks, K.W.; Misovic, B.; van den Hout, M.C.; Kockx, C.E.; Gomez, C.P.; Brouwer, R.W.; Vrieling, H.; Hoeijmakers, J.H.; van IJcken, W.F.; Pothof, J. Deciphering the RNA landscape by RNAome sequencing. RNA Biol. 2015, 12, 30–42. [Google Scholar] [CrossRef] [Green Version]
  40. Nagalakshmi, U.; Waern, K.; Snyder, M. RNA-Seq: A method for comprehensive transcriptome analysis. Curr. Protoc. Mol. Biol. 2010, 11, 1–13. [Google Scholar] [CrossRef]
  41. Mutz, K.O.; Heilkenbrinker, A.; Lönne, M.; Walter, J.G.; Stahl, F. Transcriptome analysis using next-generation sequencing. Curr. Opin. Biotechnol. 2013, 24, 22–30. [Google Scholar] [CrossRef]
  42. Peraica, M.; Flajs, D.; Domijan, A.M.; Ivić, D.; Cvjetković, B. Ochratoxin A Contamination of Food from Croatia. Toxins 2010, 2, 2098–2105. [Google Scholar] [CrossRef] [Green Version]
  43. Sergent, T. Differential modulation of ochratoxin A absorption across Caco-2 cells by dietary polyphenols, used at realistic intestinal concentrations. Toxicol. Lett. 2005, 159, 60–70. [Google Scholar] [CrossRef] [PubMed]
  44. Chen, R.; Deng, L.; Yu, X.; Wang, X.; Zhu, L.; Yu, T.; Zhang, Y.; Zhou, B.; Xu, W.; Chen, L.; et al. MiR-122 partly mediates the ochratoxin A-induced GC-2 cell apoptosis. Toxicol. In Vitro 2015, 30, 264–273. [Google Scholar] [CrossRef] [PubMed]
  45. Giromini, C.; Rebucci, R.; Fusi, E.; Rossi, L.; Saccone, F.; Baldi, A. Cytotoxicity, apoptosis, DNA damage and methylation in mammary and kidney epithelial cell lines exposed to ochratoxin A. Cell Biol. Toxicol. 2016, 32, 249–258. [Google Scholar] [CrossRef] [PubMed]
  46. Zhao, D.Y.; Zhang, W.X.; Qi, Q.Q.; Long, X.; Li, X.; Yu, Y.B.; Zuo, X.L. Brain-derived neurotrophic factor modulates intestinal barrier by inhibiting intestinal epithelial cells apoptosis in mice. Physiol. Res. 2018, 67, 475–485. [Google Scholar] [CrossRef] [PubMed]
  47. Rai, M.F.; Tycksen, E.D.; Sandell, L.J.; Brophy, R.H. Advantages of RNA-seq compared to RNA microarrays for transcriptome profiling of anterior cruciate ligament tears. J. Orthop. Res. 2018, 36, 484–497. [Google Scholar] [CrossRef] [PubMed]
  48. Assaf, H.; Azouri, H.; Pallardy, M. Ochratoxin A induces apoptosis in human lymphocytes through down regulation of Bcl-xL. Toxicol. Sci. 2004, 79, 335–344. [Google Scholar] [CrossRef] [Green Version]
  49. Ozcan, Z.; Gul, G.; Yaman, I. Ochratoxin A activates opposing c-MET/PI3K/Akt and MAPK/ERK 1-2 pathways in human proximal tubule HK-2 cells. Arch. Toxicol. 2015, 89, 1313–1327. [Google Scholar] [CrossRef]
  50. Bouaziz, C.; Sharaf El Dein, O.; El Golli, E.; Abid-Essefi, S.; Brenner, C.; Lemaire, C.; Bacha, H. Different apoptotic pathways induced by zearalenone, T-2 toxin and ochratoxin A in human hepatoma cells. Toxicology 2008, 254, 19–28. [Google Scholar] [CrossRef]
  51. Vettorazzi, A.; van Delft, J.; de Cerain López, A. A review on ochratoxin A transcriptomic studies. Food Chem. Toxicol. 2013, 59, 766–783. [Google Scholar] [CrossRef]
  52. Zhang, Y.; Qi, X.; Zheng, J.; Luo, Y.; Huang, K.; Xu, W. High-Throughput Tag-Sequencing Analysis of Early Events Induced by Ochratoxin A in HepG-2 Cells. J. Biochem. Mol. Toxicol. 2016, 30, 29–36. [Google Scholar] [CrossRef]
  53. Vanacloig-Pedros, E.; Proft, M.; Pascual-Ahuir, A. Different Toxicity Mechanisms for Citrinin and Ochratoxin A Revealed by Transcriptomic Analysis in Yeast. Toxins 2016, 8, 273. [Google Scholar] [CrossRef] [PubMed]
  54. Boesch-Saadatmandi, C.; Matzner, N.; Matzner, N.; Lang, F.; Blank, R.; Wolffram, S.; Blaschek, W.; Rimbach, G. Ochratoxin A lowers mRNA levels of genes encoding for key proteins of liver cell metabolism. Cancer Genom. Proteom. 2008, 5, 319–332. [Google Scholar]
  55. Liang, R.; Shen, X.L.; Zhang, B.; Li, Y.; Xu, W.; Zhao, C.; Luo, Y.; Huang, K. Apoptosis Signal-regulating Kinase 1 promotes Ochratoxin A-induced renal cytotoxicity. Sci. Rep. 2015, 5, 8078–8089. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Gekle, M.; Schwerdt, G.; Freudinger, R.; Mildenberger, S.; Wilflingseder, D.; Pollack, V.; Dander, M.; Schramek, H. Ochratoxin A induces JNK activation and apoptosis in MDCK-C7 cells at nanomolar concentrations. J. Pharmacol. Exp. Ther. 2000, 293, 837–844. [Google Scholar] [PubMed]
  57. Darif, Y.; Mountassif, D.; Belkebir, A.; Zaid, Y.; Basu, K.; Mourad, W.; Oudghiri, M. Ochratoxin A mediates MAPK activation, modulates IL-2 and TNF-alpha mRNA expression and induces apoptosis by mitochondria-dependent and mitochondria-independent pathways in human H9 T cells. J. Toxicol. Sci. 2016, 41, 403–416. [Google Scholar] [CrossRef]
  58. Kaminska, B. MAPK signalling pathways as molecular targets for anti-inflammatory therapy—From molecular mechanisms to therapeutic benefits. Biochim. Biophys. Acta 2005, 1754, 253–262. [Google Scholar] [CrossRef]
  59. Horvath, A.; Upham, B.L.; Ganev, V.; Trosko, J.E. Determination of the epigenetic effects of ochratoxin in a human kidney and a rat liver epithelial cell line. Toxicon 2002, 40, 273–282. [Google Scholar] [CrossRef]
  60. Sauvant, C.; Holzinger, H.; Gekle, M. The nephrotoxin ochratoxin A induces key parameters of chronic interstitial nephropathy in renal proximal tubular cells. Cell Physiol. Biochem. 2005, 15, 125–134. [Google Scholar] [CrossRef]
  61. Sharma, R.P.; He, Q.; Johnson, V.J.; Voss, K.A. Increased expression of CD95-ligand and other apoptotic signaling factors by fumonisin B1, a hepatotoxic mycotoxin, in livers of mice lacking tumor necrosis factor α. Cytokine 2003, 24, 226–236. [Google Scholar] [CrossRef]
  62. Wang, L.; Feng, Z.; Wang, X.; Zhang, X. DEGseq: An R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics 2010, 26, 136–138. [Google Scholar] [CrossRef]
  63. Hibi, D.; Kijima, A.; Suzuki, Y.; Ishii, Y.; Jin, M.; Sugita-Konishi, Y.; Yanai, T.; Nishikawa, A.; Umemura, T. Effects of p53 knockout on ochratoxin A-induced genotoxicity in p53-deficient gpt delta mice. Toxicology 2013, 304, 92–99. [Google Scholar] [CrossRef] [PubMed]
  64. Shangary, S.; Wang, S. Targeting the MDM2-p53 interaction for cancer therapy. Clin. Cancer Res. 2008, 14, 5318–5324. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Fridman, J.S.; Lowe, S.W. Control of apoptosis by p53. Oncogene 2003, 22, 9030–9340. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Hibi, D.; Kijima, A.; Kuroda, K.; Suzuki, Y.; Ishii, Y.; Jin, M.; Nakajima, M.; Sugita-Konishi, Y.; Yanai, T.; Nohmi, T.; et al. Molecular mechanisms underlying ochratoxin A-induced genotoxicity: Global gene expression analysis suggests induction of DNA double-strand breaks and cell cycle progression. J. Toxicol. Sci. 2013, 38, 57–69. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Arbillaga, L.; Vettorazzi, A.; Gil, A.G.; van Delft, J.H.; García-Jalón, J.A.; de Cerain, A.L. Gene expression changes induced by ochratoxin A in renal and hepatic tissues of male F344 rat after oral repeated administration. Toxicol. Appl. Pharmacol. 2008, 230, 197–207. [Google Scholar] [CrossRef] [PubMed]
  68. Vassilev, L.T.; Uesugi, M. In vivo activation of the p53 pathway by small-molecule antagonists of MDM2. Science 2004, 303, 844–848. [Google Scholar] [CrossRef] [Green Version]
  69. De Rozieres, S.; Maya, R.; Oren, M.; Lozano, G. The loss of mdm2 induces p53-mediated apoptosis. Oncogene 2000, 19, 1691–1697. [Google Scholar] [CrossRef] [Green Version]
  70. Silke, J.; Meier, P. Inhibitor of Apoptosis (IAP) Proteins-Modulators of Cell Death and Inflammation. Cold Spring Harb. Perspect. Biol. 2013, 5, 008730. [Google Scholar] [CrossRef]
  71. Picksley, S.M.; Lane, D.P. The p53-mdm2 autoregulatory feedback loop: A paradigm for the regulation of growth control by p53? Bioessays 1993, 15, 689–690. [Google Scholar] [CrossRef]
  72. Gu, H.; Wang, X.; Rao, S.; Wang, J.; Zhao, J.; Ren, F.L.; Mu, R.; Yang, Y.; Qi, Q.; Liu, W.; et al. Gambogic acid mediates apoptosis as a p53 inducer through down-regulation of mdm2 in wild-type p53-expressing cancer cells. Mol. Cancer Ther. 2008, 7, 3298–3305. [Google Scholar] [CrossRef] [Green Version]
  73. Alshatwi, A.A.; Subash-Babu, P.; Antonisamy, P. Violacein induces apoptosis in human breast cancer cells through up regulation of BAX, p53 and down regulation of MDM2. Exp. Toxicol. Pathol. 2016, 68, 89–97. [Google Scholar] [CrossRef] [PubMed]
  74. Kuroda, K.; Hibi, D.; Ishii, Y.; Takasu, S.; Kijima, A.; Matsushita, K.; Masumura, K.I.; Watanabe, M.; Sugita-Konishi, Y.; Sakai, H.; et al. Ochratoxin A induces DNA double-strand breaks and large deletion mutations in the carcinogenic target site of gpt delta rats. Mutagenesis 2014, 29, 27–36. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Chapman, J.R.; Taylor, M.R.G.; Simon, J. Boulton, Playing the End Game: DNA Double-Strand Break Repair Pathway Choice. Mol. Cell 2012, 47, 497–510. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Pfohl-Leszkowicz, A. MESNA protects rats against nephrotoxicity but not carcinogenity induced by Ochratoxin A, implicating two separate pathways. Med. Biol. 2002, 9, 6357–6361. [Google Scholar]
  77. Guikema, J.E.; Amiot, M.; Eldering, E. Exploiting the pro-apoptotic function of NOXA as a therapeutic modality in cancer. Expert Opin. Ther. Targets 2017, 21, 767–779. [Google Scholar] [CrossRef]
  78. Chen, H.C.; Kanai, M.; Inoue-Yamauchi, A.; Tu, H.C.; Huang, Y.; Ren, D.; Kim, H.; Takeda, S.; Reyna, D.E.; Chan, P.M.; et al. An interconnected hierarchical model of cell death regulation by the BCL-2 family. Nat. Cell Biol. 2015, 17, 1270–1281. [Google Scholar] [CrossRef] [Green Version]
  79. Comim, C.M.; Tatiana, B.; Denis, G.; Felipe, D.P.; Joao, Q.; Stephen, L. Caspase-3 mediates in part hippocampal apoptosis in sepsis. Mol. Neurobiol. 2013, 47, 394–398. [Google Scholar] [CrossRef]
  80. Larsen, B.D.; Sorensen, C.S. The caspase-activated DNase: Apoptosis and beyond. FEBS J. 2017, 284, 1160–1170. [Google Scholar] [CrossRef]
  81. D’Amelio, M.; Sheng, M.; Cecconi, F. Caspase-3 in the central nervous system: Beyond apoptosis. Trends Neurosci. 2012, 35, 700–709. [Google Scholar] [CrossRef]
  82. Gao, Y.N.; Wang, J.Q.; Li, S.L.; Zhang, Y.D.; Zheng, N. Aflatoxin M1 cytotoxicity against human intestinal Caco-2 cells is enhanced in the presence of other mycotoxins. Food Chem. Toxicol. 2016, 96, 79–89. [Google Scholar] [CrossRef]
  83. Gao, Y.; Li, S.L.; Bao, X.Y.; Luo, C.C.; Yang, H.G.; Wang, J.Q.; Zhao, S.G.; Zheng, N. Transcriptional and Proteomic Analysis Revealed a Synergistic Effect of Aflatoxin M1 and Ochratoxin A Mycotoxins on the Intestinal Epithelial Integrity of Differentiated Human Caco-2 Cells. J. Proteom. Res. 2018, 17, 3128–3142. [Google Scholar] [CrossRef] [PubMed]
  84. Trapnell, C.; Roberts, A.; Goff, L.; Pertea, G.; Kim, D.; Kelley, D.R.; Pimentel, H.; Salzberg, S.L.; Rinn, J.L.; Pachter, L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 2012, 7, 562–578. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Ernst, J.; Bar-Joseph, Z. STEM: A tool for the analysis of short time series gene expression data. BMC Bioinf. 2006, 7, 191. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Kanehisa, M.; Araki, M.; Goto, S.; Hattori, M.; Hirakawa, M.; Itoh, M.; Katayama, T.; Kawashima, S.; Okuda, S.; Tokimatsu, T.; et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2007, 36, 480–484. [Google Scholar] [CrossRef]
  87. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta, C.(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
Toxins 12 00023 g001 550
Figure 1. Effects of ochratoxin A (OTA) on apoptosis in differentiated Caco-2 cells. Incubation of cells for 48 h with OTA (0, 0.0005, 0.005, and 4 μg/mL) was followed by analysis using flow cytometry. (a) Representative flow cytometry plots are presented for the control group (CTL) and OTA groups (0.0005 μg/mL, 0.005 μg/mL, and 4 μg/mL). (b) Histogram showing number of differentially expressed genes after treatment with three concentrations of OTA compared with the control group. * represents a significant difference (p < 0.05).
Figure 1. Effects of ochratoxin A (OTA) on apoptosis in differentiated Caco-2 cells. Incubation of cells for 48 h with OTA (0, 0.0005, 0.005, and 4 μg/mL) was followed by analysis using flow cytometry. (a) Representative flow cytometry plots are presented for the control group (CTL) and OTA groups (0.0005 μg/mL, 0.005 μg/mL, and 4 μg/mL). (b) Histogram showing number of differentially expressed genes after treatment with three concentrations of OTA compared with the control group. * represents a significant difference (p < 0.05).
Toxins 12 00023 g001
Toxins 12 00023 g002 550
Figure 2. Differentially expressed genes (DEGs) identified in differentiated Caco-2 cells after exposure to 0.0005, 0.005, and 4 μg/mL OTA for 48 h. Histogram shows the number of DEGs compared with the control (CTL).
Figure 2. Differentially expressed genes (DEGs) identified in differentiated Caco-2 cells after exposure to 0.0005, 0.005, and 4 μg/mL OTA for 48 h. Histogram shows the number of DEGs compared with the control (CTL).
Toxins 12 00023 g002
Toxins 12 00023 g003 550
Figure 3. Gene expression profiles of differentiated Caco-2 cells treated with ochratoxin A (OTA). (a) Venn diagram depicting upregulated DEGs common to different doses of OTA. (b) Venn diagram depicting downregulated DEGs common to different doses of OTA. (c) Hierarchical clustering of common DEGs in differentiated Caco-2 cells based on log10-transformed expression values (fragments per kilobase of transcript per million fragments mapped, FPKM).
Figure 3. Gene expression profiles of differentiated Caco-2 cells treated with ochratoxin A (OTA). (a) Venn diagram depicting upregulated DEGs common to different doses of OTA. (b) Venn diagram depicting downregulated DEGs common to different doses of OTA. (c) Hierarchical clustering of common DEGs in differentiated Caco-2 cells based on log10-transformed expression values (fragments per kilobase of transcript per million fragments mapped, FPKM).
Toxins 12 00023 g003
Toxins 12 00023 g004 550
Figure 4. Scatter plots of the top 20 enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway terms (a), profile 0; (b), profile 3. Enriched items were measured by the rich factor, q value (q < 0.05), and number of genes. TNF, tumor necrosis factor.
Figure 4. Scatter plots of the top 20 enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway terms (a), profile 0; (b), profile 3. Enriched items were measured by the rich factor, q value (q < 0.05), and number of genes. TNF, tumor necrosis factor.
Toxins 12 00023 g004
Toxins 12 00023 g005 550
Figure 5. Protein–protein interaction (PPI) network for 50 important enzymes encoded by key differentially expressed genes (DEGs) of key pathways. Different nodes represent different enzymes. Interactions between these enzymes are represented by different size nodes; the larger the node, the stronger the connectivity.
Figure 5. Protein–protein interaction (PPI) network for 50 important enzymes encoded by key differentially expressed genes (DEGs) of key pathways. Different nodes represent different enzymes. Interactions between these enzymes are represented by different size nodes; the larger the node, the stronger the connectivity.
Toxins 12 00023 g005
Toxins 12 00023 g006 550
Figure 6. (a) Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) results for eight DEGs compared with RNA-seq results. The blue bars represent RNA-seq data and the orange bars represent q-PCR data. (b) Correlation analysis between RNA-seq and qPCR.
Figure 6. (a) Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) results for eight DEGs compared with RNA-seq results. The blue bars represent RNA-seq data and the orange bars represent q-PCR data. (b) Correlation analysis between RNA-seq and qPCR.
Toxins 12 00023 g006
Toxins 12 00023 g007 550
Figure 7. Effect of OTA (0.0005, 0.005, and 4 μg/mL) on cell apoptosis-related proteins in differentiated Caco-2 cells. Protein extracts were immunoblotted for candidate proteins (a) and band densities were quantified (b). Results are the mean of three separate experiments performed in triplicate ± SD. * p < 0.05, significantly different from control group.
Figure 7. Effect of OTA (0.0005, 0.005, and 4 μg/mL) on cell apoptosis-related proteins in differentiated Caco-2 cells. Protein extracts were immunoblotted for candidate proteins (a) and band densities were quantified (b). Results are the mean of three separate experiments performed in triplicate ± SD. * p < 0.05, significantly different from control group.
Toxins 12 00023 g007
Table
Table 1. Key pathways and important related differentially expressed genes (DEGs) that are associated with the toxic effect of ochratoxin A (OTA) in differentiated Caco-2 cells. mTOR, mechanistic target of rapamycin kinase; MAPK, mitogen-activated protein kinase; TNF, tumor necrosis factor.
Table 1. Key pathways and important related differentially expressed genes (DEGs) that are associated with the toxic effect of ochratoxin A (OTA) in differentiated Caco-2 cells. mTOR, mechanistic target of rapamycin kinase; MAPK, mitogen-activated protein kinase; TNF, tumor necrosis factor.
PathwayPathway IDKey DEGs from the Pathway
Cell cycleko04110MDM2, CDK1, TP53, EP300, ATM, CDKN1B, TGFB1, CHEK1, CDC25B, CCNB1, PDPK1
P53 signaling pathwayko04115MDM2, PMAIP1, CASP3, TP53AIP, CDK1, TP53, ATM, CYCS, TSC2, CHEK1, CASP9, CCNB1, BAX
HIF-1 signaling pathwayko04066HIF1A, NFKB1, CDKN1B, MTOR, EP300, AKT1, PIK3R, EGFR, NOX1, INSR, PIK3CA
PI3k-Akt signaling pathwayko04151CDKN1B, NFKB1, CCND1, BCL2L1, MYC, PIK3CA, EGFR, MDM2, TP53, IKBKB, AKT1, HRAS, TSC2, IL3RA, MAPK1, SOS1, CASP3, MTOR, INSR, BAD, CASP9, PDPK1
mTOR signaling pathwayko04150MTOR, TNF, IKBKB, PRKAA, PIK3CA, HRAS, TSC2, MAPK1, AKT1, SOS1, INSR, PDPK1
Apoptosisko04210XIAP, BCL2L1, NFKB1, PIK3CA, CASP3, CYC, BAX, BAD, IL3RA, CASP9, CFLAR, TP53, ATM, HRAS, TRADD, MAPK1, RIPK1, CDKN1B, PMAIP1, TP53AIP, AKT1, TNF, IKBKB
Foxo signaling pathwayko04068CDKN1B, PIK3CA, MDM2, PDPK1, EGFR1, FOXO1, AKT1, IRS1, INSR, MDM2, TGFB1, IKBKB, ATM, HRAS, MAPK1, SOS1, EGFR, EP300, CCNB1, TGFB1
Insulin signaling pathwayko04910MAPK1, SOS1, PIK3CA, AKT1, MTOR, IRS1, INSR, HRAS, TSC2, TRADD, MAPK1, CALM1, IKBKB, FOXO1
MAPK signaling pathwayko04010TNF, CASP3, NFKB1, TP53, MAPK14, TGFB1, PIK3CA, HRAS, SOS1, EGFR, MAPK1, TRAF2, CDC25B, AKT1, TNF, IKBKB
TNF signaling pathwayko04668XIAP, TNF, NFKB1, CASP3, IKBKB, MAPK1, MAPK14, CASP8, TRADD, FADD, MAPK1, RIPK1, TRAF2, AKT1, PIK3CA
Table
Table 2. Primer sequences for the quantification of genes by quantitative reverse transcription-polymerase chain reaction (qRT-PCR). CASP 3, caspase 3; CXCL2, C–X–C motif chemokine ligand 2; CDC25B, cell division cycle 25B; EGR1, early growth response 1 (EGR1), H2BFS, H2B histone family members; SHPK, sedoheptulokinase; TFEC, transcription factor EC.
Table 2. Primer sequences for the quantification of genes by quantitative reverse transcription-polymerase chain reaction (qRT-PCR). CASP 3, caspase 3; CXCL2, C–X–C motif chemokine ligand 2; CDC25B, cell division cycle 25B; EGR1, early growth response 1 (EGR1), H2BFS, H2B histone family members; SHPK, sedoheptulokinase; TFEC, transcription factor EC.
GenesProduct Length (bp)Forward Primer Sequence (5′–3′)Reverse Primer Sequence (5′–3′)
GAPDH235GGAGTCCACTGGCGTCTTGAGTCCTTCCACGATACCAAA
CASP3109TCCTGAGATGGGTTTATGTTGTTTCCCTGAGGTTTGC
CXCL2150CCAAACCGAAGTCATAGCGAACAGCCACCAATAAGC
CDC25B296GTAGACGGAAAGCACCAAGATCCCTGATGAAACGGCAC
EGR1229CACGAACGCCCTTACGCTCATCGCTCCTGGCAAACT
H2BFS119TGCTCGTCTCAGGCTCGTAGCTTCCTGCCGTCCTTCTTCT
SHPK58AGTAGATGCGGCAATGGTTTGGTAGGGATGGCTGTG
TEFC94GCACTGGAGGGATAAATGTAAAGACACCCGAAGGAT

Share and Cite

MDPI and ACS Style

Yang, X.; Gao, Y.; Yan, Q.; Bao, X.; Zhao, S.; Wang, J.; Zheng, N. Transcriptome Analysis of Ochratoxin A-Induced Apoptosis in Differentiated Caco-2 Cells. Toxins 2020, 12, 23. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12010023

AMA Style

Yang X, Gao Y, Yan Q, Bao X, Zhao S, Wang J, Zheng N. Transcriptome Analysis of Ochratoxin A-Induced Apoptosis in Differentiated Caco-2 Cells. Toxins. 2020; 12(1):23. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12010023

Chicago/Turabian Style

Yang, Xue, Yanan Gao, Qiaoyan Yan, Xiaoyu Bao, Shengguo Zhao, Jiaqi Wang, and Nan Zheng. 2020. "Transcriptome Analysis of Ochratoxin A-Induced Apoptosis in Differentiated Caco-2 Cells" Toxins 12, no. 1: 23. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12010023

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop