- Research article
- Open Access
Microarray based gene expression analysis of Sus Scrofa duodenum exposed to zearalenone: significance to human health
BMC Genomics volume 17, Article number: 646 (2016)
Zearalenone (ZEA) is a secondary metabolite produced by Fusarium species. ZEA was proved to exert a wide range of unwanted side effects, but its mechanism of action, particularly at duodenum levels, remains unclear. In our study based on the microarray technology we assessed the alteration of gene expression pattern Sus scrofa duodenum which has been previously exposed to ZEA. Gene expression data was validated by qRT-PCR and ELISA. The gene expression data were further extrapolated the results to their human orthologues and analyzed the data in the context of human health using IPA (Ingenuity Pathways Analysis).
Using Agilent microarray technology, we found that gene expression pattern was significantly affected by ZEA exposure, considering a 2-fold expression difference as a cut-off level and a p-value < 0.05. In total, we found 1576 upregulated and 2446 downregulated transcripts. About 1084 genes (764 downregulated and 751 overexpressed) were extrapolated to their human orthologues. IPA analysis showed various altered key cellular and molecular pathways. As expected, we observed a significant alteration of immune response related genes, MAPK (mitogen activate protein kinases) pathways or Toll-Like Receptors (TLRs). What captured our attention was the modulation of pathways related to the activation of early carcinogenesis.
Our data demonstrate that ZEA has a complex effect at duodenum level. ZEA is able to activate not only the immune response related genes, but also those relate to colorectal carcinogenesis. The effects can be more dramatic when connected with the exposure to other environmental toxic agents or co-occurrence with different microorganisms.
Zearalenone (ZEA) is a secondary metabolite produced by certain species that belong to the genus Fusarium, a common group of fungus species [1, 2]. These fungi are common contaminants present is almost all types of crops, including those used as feed for farm animal diets as well as cereals consumed by humans [3–5]. The ZEA mycotoxin has a particular way of exerting its effects by interfering in the physiological estrogen signaling pathways, since its chemical structure is similar to that of the estrogenic hormone 17-β-estradiol [6–8]. Among the different farm animals, pigs proved to be highly sensitive to the effects of ZEA mycotoxin via contaminated feed. They are known to display various deleterious effects caused by this mycotoxin, but, regardless of the estrogenic activity of ZEA, many of the observed damages are, presumably, not always mediated by the estrogen receptor [9, 10]. Also, the histology, anatomy, physiology – and consequently the pathology – of this species are very similar to that of humans, making it a suitable model for the study of various diseases and physio-pathological processes [11, 12].
The intestine is the interface between the organism and the whole exposome, so it interacts with the existent microflora and with the pathogenic agents that may occur [13–15]. These particular tissues were chosen because, as part of the gastro-intestinal (GI) tract, they are involved in the immune response and inflammation on the one hand, and they are among the most exposed to the action of the contaminants on the other hand. This makes these organs suitable subjects for this microarray profiling study [16–20].
The animal model used for this study is the pig, due to the previously mentioned genetic similarity between Sus scrofa and Homo sapiens. At the same time, pigs consume high quantities of maize, a cereal which is very prone to Fusarium mycotoxin contamination , and which is also a raw material in the diet of humans. Such studies are of particular interest because humans also consume high quantities of different cereals which are, many times, contaminated with various mycotoxins such as ZEA [9, 22]. Even though there are regulations regarding the maximum tolerated values of food contaminating agents, including ZEA , some such toxins have shown increased resistance to decontaminating procedures, and a wide range of side effects [24, 25]. Therefore it is important to acquire a better understanding of the influence that ZEA has on the health of humans and their surrounding environment.
Evaluation of duodenal gene expression pattern
A significant supposition in many toxicological investigations is that the molecular states of organisms indicate their biological responses to a particular toxin, like in our case ZEA mycotoxin. Using the Agilent microarray technology, we found that gene expression was considerably affected by ZEA at duodenum level, considering 2-fold expression difference as a cut-off level and p-value < 0.05. In total, we found 4023 transcripts with an altered expression level, from which 1576 were upregulated and 2446 downregulated. Sus scrofa is a key mammalian model system for studying complex human diseases and therefore it is useful to assess the impact of the toxin on this model, then to extrapolate the gene expression profile and to analyze the data in the context of human health. We were able to obtain the human orthologues for about 1084 genes (764 downregulated and 751 overexpressed).
Microarray data validation by qRT-PCR
The obtained microarray data were validated by qRT-PCR. Therefore four genes were selected (NFKB1, IL-6, TNF-α and SOD2) and three housekeeping genes (ACTB, GAPDH, B2M) were used for the normalization of the data. qRT-PCR data confirm the microarray downregulation expression of these genes, moreover it shows a good correlation among them (Fig. 1).
The impact of ZEA exposure on cytokine protein expression at duodenum levels
The microarray data displayed previously showed an alteration of molecules involved in immune response. Therefore, we assessed the protein level for IL-1β, IL-8, IL-4 and TNF-α for cellular lysates. As it can be seen from Fig. 2, we observed a downregulation of IL-1β, IL-8, IL-4 at protein level as a response to ZEA exposure. The level of TNF-α was under the limit of detection provided by our method.
A primary goal of the study was to identify the possible implications of the altered genes at cellular and molecular level, as well as the related functions and pathways that might be mediated by ZEA (Table 1). To accomplish this objective in an unbiased manner, we performed IPA analysis for all the genes with and altered expression levels that were extrapolated to their human orthologues. This facilitated the assessment of the potential impact by analyzing the networks, biological functions, and canonical pathways. Additional file 1: Figure S1 presents the top canonical pathways, based on the overlap value, and displays the genes related to the alteration of Toll-Like Receptors (TLRs) and the activation of the inflammatory cytokine in parallel with the alteration of the expression level for the adhesion molecules. In Additional file 2: Figure S2 are emphasized the altered genes related to MAPK (mitogen activate protein kinases), an early event of carcinogenesis, fact demonstrated in Additional file 3: Figure S3. Also, we observed alterations in GAP junction signaling (Additional file 4: Figure S4).
The analysis of probe sets that are differentially expressed with a fold change greater than 2 revealed 24 altered pathways, specific for the toxin exposure. Most of the altered networks were related to the alteration of normal cell status, the most significant network containing 34 focus molecules (Fig. 3), with the core of this network being the FOS gene. Another interesting network is that represented by the one with the ID7 (Table 2), related to carcinogenesis, cell death and organism damage (Figs. 3 and 4 ), the core or this network being represented by ER. A curious aspect of the gene expression alteration at duodenum level is represented by the network with the ID20 (Table 1), where the core of the network is the RB1 gene; this gene is recognized to be involved in TGFB1 or p53 signaling cell cycle regulation, but also to be related to mitochondrial damage.
ZEA mycotoxin is absorbed up to 90 % in the upper part of the GIT (gastrointestinal tract), and it goes into enterohepatic circulation as other mycotoxins . One of the important roles of the GIT is its function as an immune barrier. This function is accomplished through a number of particularities. Firstly, it possesses its own immune system and it is estimated that up to 70 % of the immune defenses of the organism are located in the intestine. Secondly, its morphology plays a role as a physical barrier through tight junctions (TJs) formed mainly from occludin and claudin isoforms proteins, and gap junctions (GJs) that permit the transfer of ions, nucleotides and other small molecules between adjacent cells, and are formed mainly by connexins. Last but not least, the intestinal microbiota plays a very important role in protecting against pathogen invasion . Moreover, the physical barrier can be disrupted due to defective TJs and GJs. ZEA has been shown to reduce mRNA levels of occludin and claudin-4 and also the protein levels of connexin , a mechanism that is altered by cytotoxic and carcinogenic agents , fact demonstrated also by present study, via multiple mechanisms relate to cell-to-cell communication.
By quantifying the modified expression of tens of thousands of genes, we can recognize mechanistically-relevant genes, leading to a better understanding of the ZEA toxicological effects at duodenum level.
ZEA was proved to affect the alteration of several genes involves in inflammatory response, oxidative stress, proteinases or other relevant genes . The capacity to activate oxidative stress was proved also in an in vitro study . It was reported that ZEA has pro-apoptotic effects, by increasing the Bax/Bcl-2 ration . ZEN was observed to be involved in hepatic injury via Nrf2/ARE signaling pathway , mechanism altered also at intestinal level. Therefore, continuous exposure might be related to the activation of the tumor-promoting inflammatory mechanisms , expressed by the inflammatory cytokines or by the antioxidant enzyme, these being validated by qRT-PCR. It is clear that ZEA leads to an inflammatory microenvironment, and this can be correlated with increased mutation rates, all these representing early carcinogenic events. This fact sustained by other environmental toxic agents can be an essential in the activation of tumor initiation . In the case of HCT116 cells it was observed that even at low concentration ZEA increase cell proliferation, colony formation and cell migration .
ZEA compounds were observed to exert estrogenic effects [7, 36], even at nM concentration , fact confirmed by the present study. The estrogenic effects seem to be altered via Erk1/2 phosphorylation, mechanism, in parallel with the sustaining of cell proliferation . ZEN could promote the autocrine action or expression of the ghrelin gene in piglet ovary .
Other studies confirmed the activation of MAPK and oxidative stress related genes as response to ZEA in MCF-7 breast cancer cells,  and macrophages , this findings being confirmed by this present in vivo study at duodenum level. The oxidative species seems to be generated via p53, JNK or p38 pathways [31, 39]. In our study NFKB1 was found to be down-regulated, also it was connected with a proapoptotic effect in colon carcinoma cells. c-Myc is an early response gene with important regulatory functions in the cell cycle, and known to be regulated in colon cancer. c-Myc was downregulated by ZEA exposure. In the same time, the IPA network analysis showed the complex action of ZEA, the analysis being more focused on the alteration of the mitochondrial pathways connected with endocrine dysfunction . There are several classes of altered genes, such as Ras-related genes, kinases/phosphatases and their binding proteins, genes retrieved also in study. These Ras pathways contribute to cell growth and proliferation .
Alteration of Stathmin1 pathways are connected with invasion and migration capacity in digestive cancers . The oncogenic role of the Stathmin1 pathways was proved to be related via miR-223 . In a recent review paper Reuter et al., (2010) present a clear connection between oxidative stress and the activation of the transcription factors (NFkB, AP-1, HIF-1α, PPAR-γ, Wnt and Nrf2) . The present study confirms that the activation of these transcription factors is related to the alteration of other multiple genes like growth factors, cytokines/chemokines, cell cycle related genes, cell proliferation or differentiation anti-inflammatory molecules, or those in invasion and migration. From the top 10 genes that were highly expressed in the case of duodenum exposure to ZEA, a representative one is CTNNA1. CTNNA1, along with other MAPK related genes [31, 39], have been proved to be involved in gastric cancer . NFAT-signaling pathway activation was proved to be implicated in aggressive forms of inflammatory pathologies. Moreover, this pathway was proved to have a significant role in inflammatory breast cancer .
The results remains to be further investigated particular on ZEA derivatives, a recent study present antitumoral effects of the semisynthetic zearalenone analogues . Therefore the mechanisms of transformation and detoxification in parallel with the combinatorial effect need to be taken in account . The procarcinogenic effect of ZEA was also studied by Ouanes et al., 2003 , and should take into account the accumulation in organism and the time required for metabolisation of this toxins . Novel monitoring method based on urinary excretion of ZEA metabolites should be developed .
Policies and regulations of the mycotoxin need to be assessed in the context of the exposome exposure, not only guided by compound exposure and evaluation of the dose-response analyses in risk assessment [51, 52]. The impact of mycotoxins on human health needs to be considered in a transformative mode to take into account the whole exposome exposure . Even in the case of a single toxic exposure (ZEA), the effect on gene expression is significant at duodenum level. There are several clues that environmental toxic exposure implies , including the case of ZEA. This toxin can interact and affect more dramatically health status that is presented by microarray data. ZEA is only one component of the exposome, therefore ZEA toxic effects can be more dramatic in the context of multiple exposures.
As ZEA caused alteration of gene expression, chronic exposure can create a metabolic environment that allows mutated cells to acquire a selective advantage. Our data demonstrate that ZEA activates not only the immune related genes response, but also those related to carcinogenesis. The effect can be more dramatic when is connected with other environmental toxic agents and this should be taken into account for further implementation of novel norms for the control of this toxin.
Animals and treatments
Two groups of weaned TOPIG piglets [(Landrace × Large White) × (Duroc × Pietrain)], n = 6 per group/pen) with an initial average body weight of 9.88 ± 0.2 Kg and derived from the experimental farm of the National Research Institute for Biology and Animal Nutrition, Balotesti Romania were used in this study. The animals individually identified by ear tag were housed in pens and fed with experimental diets: a control diet without mycotoxin or a diet contaminated with ZEA (100 ppb) for 30 days. Assigned diet and water were provided ad libitum every day of the experiment. Pigs were observed twice daily and no clinical signs or death was recorded throughout the entire experimental period. At the end of the experiment, animals were stunned and slaughtered in an EU-licensed abattoir according with the EU Council Directive 2010/63/CE. Organ samples were collected on ice from all animals, weighed and were stored at –80 °C until the analyses.
Animals were raised in agreement with the Romanian Law 43/2014 for handling and protection of animals used for experimental purposes and the EU Council Directive 98/58/EC concerning the protection of farmed animals. The study protocol was accepted by the Ethical Committee of the National Research-Development Institute for Animal Nutrition and Biology, Balotesti, Romania (Ethical Committee no. 52/2014). At the end of the research period, animals were stunned and slaughtered in an EU-licensed abattoir in accordance with the EU Council Directive 2010/63/CE.
RNA extraction and quality control
The duodenum tissue collected from weaned pigs exposed to ZEA for 33d was disrupted using a mechanical homogenizer, in the presence of TRI Reagent (Sigma-Aldrich), and the total RNA was extracted following the phase separation protocol with chloroform and isopropanol, according to the manufacturer’s protocol. The RNA thus obtained underwent a further purification step with the RNeasy Micro Kit from Qiagen, which uses silicagel spin columns to generate high quality RNA from small quantities of tissue. The precise concentrations and integrity of extracted nucleic acids, necessary for the subsequent microarray step were then evaluated using the NanoDrop2000 (Thermo Scientific) and the Bioanalyzer 2100 (Agilent Technologies).
The microarray probes were synthesized from equal quantities of 100 ng of total RNA, by using Agilent Low Input Quick Amp Labeling Kit (5190-2305) according to the manufacturer’s protocol. Subsequent to this step, the hybridization products were purified using RNeasy Mini kit (Qiagen). Probe quality control was executed using the NanoDrop2000 spectrophotometer (Thermo Scientific), all the samples having a specific activity higher than 6 pmol/μl Cy3/μg cRNA (specific activity > 8 pmol Cy3/μg cRNA). The fragmentation and hybridization were performed based on the Agilent one color protocol.
The samples were hybridized (17 h at 65 °C) on a custom microarray slide using SureDesign, containing up of 60mer oligos for over 59,000 probes representative for a high number of Sus scrofa transcripts (AMADID 056850). The microarray slides were scanned with the SureScan Microarray Scanner (1x60k array slides with 61 × 21 mm, resolution 3 μM) from Agilent and the images were processed with Feature Extraction 18.104.22.168 software. GeneSpring GX v.13.0 from Agilent was used for data analysis, which started with the removal of control probes and quantile normalization. A threshold of 2.0 was applied for the fold change, and then the moderated t test was used, together with the FDR (False Discovery Rate) correction method and extrapolated to human as described by Petric et al. .
The information related to homology annotation for all the transcripts sequence for S. scrofa and humans was used NCBI database (www.ncbi.nlm.nih.gov) and BLAST. The analysis in the context of human health for the altered genes as response to ZEA exposure was done by ingenuity pathway analysis (IPA; htp://www.ingenuity.com).
qRT-PCR data validation
For the validation of the gene profiles provided by microarray analysis, the expression levels of five randomly selected genes: TNF-α (Tumor Necrosis Factor α), IL-8 (Interleukin 8), IL-6 (Interleukin 6), IL-1β (Interleukin 1 beta), IL-10 (Interleukin 10) and NFkB1 (nuclear factor of kappa light polypeptide gene enhancer in B-cells 1) were measured by real time RT-qPCR in all samples used for microarray analysis. 1 μg of total RNA extracted from each sample was used to generate cDNA using M-MLV Revers Trascriptase kit (Invitrogen, Life Technologies), according to the manufacturer’s recommendations. Quantitative fluorescent real-time PCR reactions were set up in a total volume of 20 μl using 10 μl of cDNA (diluted 1:10 with nuclease-free water), SYBR Green PCR Master Mix (Applied Biosystems, Life Technologies), 0.3 μM each of gene-specific primer and performed in the Rotor- Gene-Q (QIAGEN GmbH, Germany) machine. The primer pairs used in the present study, listed in Table 3, were obtained from Eurogentec (San Diego, USA). The cycling conditions were set up according to manufacturer’s recommendation, as follows: UDG pre-treatment at 50 °C for 2 min, initial denaturation step at 95 °C for 15 s, followed by 40 cycles of 95 °C for 15 s, 60 °C for 60s with a single fluorescence measurement. The specificity of the PCR products was confirmed by analysis of the dissociation curve. The melting curve programme consisted of temperatures between 60 and 95 °C with a heating rate of 0.1 °C/s and continuous fluorescence measurement. Negative controls were used for each primer pair. The relative quantification of gene expression changes were quantified using the comparative method. The expression levels of two reference genes, GAPDH and β-actin were used for data normalisation. These reference genes were experimentally validated for duodenum tissues and the lack of treatment effect and expression variation was the criteria for reference gene choice. The results were expressed as relative fold change (Fc) in comparison with control samples.
Cytokine concentration evaluation (ELISA)
Organ lysate was prepared as described by Taranu et al., (2014)  and cytokine concentrations in the supernatant were determined by ELISA. Briefly, the purified fraction of anti-swine cytokines (IL-1β, IL-8, IL-4) were used as capture antibodies in conjunction with biotinylated anti-swine antibodies (R&D system). Streptavidin-HRP (Sigma) and TMB (Sigma) were used for detection. Absorbance was read at 450 nm using an ELISA plate reader (Tecan, Sunrise, Austria). Recombinant swine cytokine proteins were used as standard and results were expressed as picograms of cytokine/mL.
ACTB, Beta-actin; AP-1, Activator protein 1; ARE, antioxidant response element; B2M, Beta-2-Microglobulin; Bax, BCL2 Associated X Protein; Bcl-2, B-cell lymphoma 2; CTNNA1, Catenin Alpha 1; ER, Estrogen Recptor; Erk1/2, extracellular signal–regulated kinases 1/2; FDR, False Discovery Rate; GAPDH , Glyceraldehyde-3-Phosphate Dehydrogenase; GAPDH, Glyceraldehyde-3-Phosphate Dehydrogenase; GIT, gastrointestinal tract; GJs, gap junctions; HIF-1α , Hypoxia-inducible factor 1-alpha; IL-10, Interleukin 10; IL-1β, Interleukin 1 beta; IL-4, Interleukin 4; IL-6, Interleukin 6; IL-8, Interleukin 8; IPA, Ingenuity Pathways Analysis; JNK, Jun amino-terminal kinases; MAPK, mitogen activate protein kinases; NFAT, Nuclear factor of activated T-cells; NFKB1, nuclear factor of kappa light polypeptide gene enhancer in B-cells 1; Nrf2, nuclear erythroid 2-related factor 2; PPAR-γ, Peroxisome Proliferator-activated Receptor γ; RB1, retinoblastoma; SOD2, Superoxide dismutase 2; TGFB1, Transforming Growth Factor Beta 1; TJs, tight junctions; TLRs - Toll-Like Receptors; TNF-α, Tumor Necrosis Factor α; ZEA, Zearalenone
Caldwell RW, Tuite J, Stob M, Baldwin R. Zearalenone production by Fusarium species. Appl Microbiol. 1970;20(1):31–4.
Hidy PH, Baldwin RS, Greasham RL, Keith CL, McMullen JR. Zearalenone and some derivatives: production and biological activities. Adv Appl Microbiol. 1977;22:59–82.
Yazar S, Omurtag GZ. Fumonisins, trichothecenes and zearalenone in cereals. Int J Mol Sci. 2008;9(11):2062–90.
Fazekas B, Tar A. Determination of zearalenone content in cereals and feedstuffs by immunoaffinity column coupled with liquid chromatography. J AOAC Int. 2001;84(5):1453–9.
Iqbal SZ, Rabbani T, Asi MR, Jinap S. Assessment of aflatoxins, ochratoxin A and zearalenone in breakfast cereals. Food Chem. 2014;157:257–62.
So MY, Tian Z, Phoon YS, Sha S, Antoniou MN, Zhang J, Wu RS, Tan-Un KC. Gene expression profile and toxic effects in human bronchial epithelial cells exposed to zearalenone. PLoS One. 2014;9(5):e96404.
Minervini F, Giannoccaro A, Cavallini A, Visconti A. Investigations on cellular proliferation induced by zearalenone and its derivatives in relation to the estrogenic parameters. Toxicol Lett. 2005;159(3):272–83.
Diekman MA, Green ML. Mycotoxins and reproduction in domestic livestock. J Anim Sci. 1992;70(5):1615–27.
Hueza IM, Raspantini PC, Raspantini LE, Latorre AO, Gorniak SL. Zearalenone, an estrogenic mycotoxin, is an immunotoxic compound. Toxins. 2014;6(3):1080–95.
Taranu I, Braicu C, Marin DE, Pistol GC, Motiu M, Balacescu L, Beridan Neagoe I, Burlacu R. Exposure to zearalenone mycotoxin alters in vitro porcine intestinal epithelial cells by differential gene expression. Toxicol Lett. 2014;232(1):310–25.
Kollarczik B, Gareis M, Hanelt M. In vitro transformation of the Fusarium mycotoxins deoxynivalenol and zearalenone by the normal gut microflora of pigs. Nat Toxins. 1994;2(3):105–10.
Piotrowska M, Slizewska K, Nowak A, Zielonka L, Zakowska Z, Gajecka M, Gajecki M. The effect of experimental fusarium mycotoxicosis on microbiota diversity in porcine ascending colon contents. Toxins. 2014;6(7):2064–81.
Tenaillon O, Skurnik D, Picard B, Denamur E. The population genetics of commensal Escherichia coli. Nat Rev Microbiol. 2010;8(3):207–17.
Sekirov I, Russell SL, Antunes LC, Finlay BB. Gut microbiota in health and disease. Physiol Rev. 2010;90(3):859–904.
Huffnagle G, Noverr MC. GI microbiota and regulation of the immune system. Preface. Adv Exp Med Biol. 2008;635:v–vi.
Salminen S, Bouley C, Boutron-Ruault MC, Cummings JH, Franck A, Gibson GR, Isolauri E, Moreau MC, Roberfroid M, Rowland I. Functional food science and gastrointestinal physiology and function. Br J Nutr. 1998;80 Suppl 1:S147–171.
Tait Wojno ED, Artis D. Innate lymphoid cells: balancing immunity, inflammation, and tissue repair in the intestine. Cell Host Microbe. 2012;12(4):445–57.
Chen X, Feng BS, Zheng PY, Liao XQ, Chong J, Tang SG, Yang PC. Fc gamma receptor signaling in mast cells links microbial stimulation to mucosal immune inflammation in the intestine. Am J Pathol. 2008;173(6):1647–56.
Maier E, Anderson RC, Roy NC. Understanding how commensal obligate anaerobic bacteria regulate immune functions in the large intestine. Nutrients. 2015;7(1):45–73.
Marin DE, Taranu I, Burlacu R, Tudor DS. Effects of zearalenone and its derivatives on the innate immune response of swine. Toxicon. 2010;56(6):956–63.
Oswald IP, Desautels C, Laffitte J, Fournout S, Peres SY, Odin M, Le Bars P, Le Bars J, Fairbrother JM. Mycotoxin fumonisin B1 increases intestinal colonization by pathogenic Escherichia coli in pigs. Appl Environ Microbiol. 2003;69(10):5870–4.
Bennett JW, Klich M. Mycotoxins. Clin Microbiol Rev. 2003;16(3):497–516.
European Food Safety Authority. Scientific opinion on the risk for public health related to the presence of zearalenone in food. Eur Food Saf Auth J. 2011;9:1–124.
Zinedine A, Soriano JM, Molto JC, Manes J. Review on the toxicity, occurrence, metabolism, detoxification, regulations and intake of zearalenone: an oestrogenic mycotoxin. Food Chem Toxicol. 2007;45(1):1–18.
Ryu D, Hanna MA, Eskridge KM, Bullerman LB. Heat stability of zearalenone in an aqueous buffered model system. J Agric Food Chem. 2003;51(6):1746–8.
Grenier B, Applegate TJ. Modulation of intestinal functions following mycotoxin ingestion: meta-analysis of published experiments in animals. Toxins. 2013;5(2):396–430.
Pinton P, Braicu C, Nougayrede JP, Laffitte J, Taranu I, Oswald IP. Deoxynivalenol impairs porcine intestinal barrier function and decreases the protein expression of claudin-4 through a mitogen-activated protein kinase-dependent mechanism. J Nutr. 2010;140(11):1956–62.
Liu M, Gao R, Meng Q, Zhang Y, Bi C, Shan A. Toxic effects of maternal zearalenone exposure on intestinal oxidative stress, barrier function, immunological and morphological changes in rats. PLoS One. 2014;9(9):e106412.
Chipman JK, Mally A, Edwards GO. Disruption of gap junctions in toxicity and carcinogenicity. Toxic Sci. 2003;71(2):146–53.
West AC, Jenkins BJ. Inflammatory and Non-Inflammatory Roles for Toll-Like Receptors in Gastrointestinal Cancer. Curr Pharm Des. 2015;21(21):2968–77.
Yu JY, Zheng ZH, Son YO, Shi X, Jang YO, Lee JC. Mycotoxin zearalenone induces AIF- and ROS-mediated cell death through p53- and MAPK-dependent signaling pathways in RAW264.7 macrophages. Toxicol in vitro. 2011;25(8):1654–63.
Hu J, Xu M, Dai Y, Ding X, Xiao C, Ji H, Xu Y. Exploration of Bcl-2 family and caspases-dependent apoptotic signaling pathway in Zearalenone-treated mouse endometrial stromal cells. Biochem Biophys Res Commun. 2016;476(4):553-9.
Long M, Yang SH, Han JX, Li P, Zhang Y, Dong S, Chen X, Guo J, Wang J, He JB. The Protective Effect of Grape-Seed Proanthocyanidin Extract on Oxidative Damage Induced by Zearalenone in Kunming Mice Liver. Int J Mol Sci. 2016;17(6). doi:10.3390/ijms17060808.
Grivennikov SI, Greten FR, Karin M. Immunity, Inflammation, and Cancer. Cell. 2010;140(6):883–99.
Abassi H, Ayed-Boussema I, Shirley S, Abid S, Bacha H, Micheau O. The mycotoxin zearalenone enhances cell proliferation, colony formation and promotes cell migration in the human colon carcinoma cell line HCT116. Toxicol Lett. 2016;254:1–7.
Parveen M, Zhu Y, Kiyama R. Expression profiling of the genes responding to zearalenone and its analogues using estrogen-responsive genes. FEBS Lett. 2009;583(14):2377–84.
Dai M, Jiang S, Yuan X, Yang W, Yang Z, Huang L. Effects of zearalenone-diet on expression of ghrelin and PCNA genes in ovaries of post-weaning piglets. Anim Reprod Sci. 2016;168:126–37.
Mike N, Papp G, Certik M, Czibulya Z, Kunsagi-Mate S, Ember I, Vagvolgyi C, Pesti M, Gazdag Z. Regulation of cytotoxic, non-estrogenic, oxidative stress-induced processes of zearalenone in the fission yeast Schizosaccharomyces pombe. Toxicon. 2013;73:130–43.
Lu J, Yu JY, Lim SS, Son YO, Kim DH, Lee SA, Shi X, Lee JC. Cellular mechanisms of the cytotoxic effects of the zearalenone metabolites alpha-zearalenol and beta-zearalenol on RAW264.7 macrophages. Toxicol in vitro. 2013;27(3):1007–17.
Busk OL, Ndossi D, Verhaegen S, Connolly L, Eriksen G, Ropstad E, Sorlie M. Relative quantification of the proteomic changes associated with the mycotoxin zearalenone in the H295R steroidogenesis model. Toxicon. 2011;58(6-7):533–42.
Jeon TY, Han ME, Lee YW, Lee YS, Kim GH, Song GA, Hur GY, Kim JY, Kim HJ, Yoon S, et al. Overexpression of stathmin1 in the diffuse type of gastric cancer and its roles in proliferation and migration of gastric cancer cells. Br J Cancer. 2010;102(4):710–8.
Kang W, Tong JH, Chan AW, Lung RW, Chau SL, Wong QW, Wong N, Yu J, Cheng AS, To KF. Stathmin1 plays oncogenic role and is a target of microRNA-223 in gastric cancer. PLoS One. 2012;7(3):e33919.
Reuter S, Gupta SC, Chaturvedi MM, Aggarwal BB. Oxidative stress, inflammation, and cancer: How are they linked? Free Radic Biol Med. 2010;49(11):1603–16.
Sugimoto S, Komatsu H, Morohoshi Y, Kanai T. Recognition of and recent issues in hereditary diffuse gastric cancer. J Gastroenterol. 2015;50(8):831-43.
Remo A, Simeone I, Pancione M, Parcesepe P, Finetti P, Cerulo L, Bensmail H, Birnbaum D, Van Laere SJ, Colantuoni V, et al. Systems biology analysis reveals NFAT5 as a novel biomarker and master regulator of inflammatory breast cancer. J Transl Med. 2015;13:138.
Tadpetch K, Kaewmee B, Chantakaew K, Kantee K, Rukachaisirikul V, Phongpaichit S: Synthesis and cytotoxic activities of semisynthetic zearalenone analogues. Bioorg Med Chem Lett. 2016;26(15):3612-6.
Warth B, Braun D, Ezekiel CN, Turner PC, Degen GH, Marko D: Bio-monitoring of mycotoxins in human breast milk: Current state and future perspectives. Chem Res Toxicol. 2016;29(7):1087-97.
Ouanes Z, Abid S, Ayed I, Anane R, Mobio T, Creppy EE, Bacha H. Induction of micronuclei by Zearalenone in Vero monkey kidney cells and in bone marrow cells of mice: protective effect of Vitamin E. Mutat Res. 2003;538(1-2):63–70.
Gajecka M, Slawuta P, Nicpon J, Kolacz R, Kielbowicz Z, Zielonka L, Dabrowski M, Szweda W, Gajecki M, Nicpon J. Zearalenone and its metabolites in the tissues of female wild boars exposed per os to mycotoxins. Toxicon. 2016;114:1–12.
Thanner S, Czegledi L, Schwartz-Zimmermann HE, Berthiller F, Gutzwiller A. Urinary deoxynivalenol (DON) and zearalenone (ZEA) as biomarkers of DON and ZEA exposure of pigs. Mycotoxin Res. 2016;32(2):69–75.
Betts KS. Characterizing exposomes: tools for measuring personal environmental exposures. Environ Health Perspect. 2012;120(4):A158–163.
Buck Louis GM, Sundaram R: Exposome: Time for Transformative Research. Statistics in medicine 2012, 31(22):10.1002/sim.5496.
Petric RC, Braicu C, Bassi C, Pop L, Taranu I, Dragos N, Dumitrascu D, Negrini M, Berindan-Neagoe I. Interspecies Gene Name Extrapolation--A New Approach. PLoS One. 2015;10(9):e0138751.
Grenier B, Bracarense AP, Schwartz HE, Trumel C, Cossalter AM, Schatzmayr G, Kolf-Clauw M, Moll WD, Oswald IP. The low intestinal and hepatic toxicity of hydrolyzed fumonisin B(1) correlates with its inability to alter the metabolism of sphingolipids. Biochem Pharmacol. 2012;83(10):1465–73.
von der Hardt K, Kandler MA, Brenn G, Scheuerer K, Schoof E, Dotsch J, Rascher W. Comparison of aerosol therapy with different perfluorocarbons in surfactant-depleted animals. Crit Care Med. 2004;32(5):1200–6.
Meadus WJ, MacInnis R, Dugan ME. Prolonged dietary treatment with conjugated linoleic acid stimulates porcine muscle peroxisome proliferator activated receptor gamma and glutamine-fructose aminotransferase gene expression in vivo. J Mol Endocrinol. 2002;28(2):79–86.
This work was supported by funds from the Romanian Ministry of Research and Technology: PNII-PCCA-102/2011-2016 with title ‘Impact of feed co-contamination and mitigating solutions to increase feed safety, animal health and food quality”.
Availability of data and materials
All the data will be available on ArrayExpress, accession number: E-MTAB-3885.
CB, IT and IB-N: designed and coordinated the study; CB, AJ, DG carried out the microarray study; CB, IB-N, RC-P performed microarray statistical analysis; A-MG, D-EM, IT, GP carried out in vivo study, microarray validation and ELISA cytokine assessment. All the authors participated for the writing of the paper and agreed to the final variant of the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Alteration of Toll-Like Receptors (TLRs) and the activation of the inflammatory cytokine in parallel with the alteration of the expression level for the adhesion molecules as response to ZEA exposure; pathway generated using IPA (Ingenuity Pathway Analysis). (JPG 5492 kb)
Alteration of genes related to MAPK (mitogen activate protein kinases), an early event of carcinogenesis, fact demonstrated as response to ZEA exposure; pathway generated using IPA. (JPG 2327 kb)
Activation of colorectal carcinogenic mechanism as response to ZEA exposure; pathway generated using IPA. (JPG 5802 kb)
Alterations in GAP junction signaling with signification on the increasing the susceptibility to bacterial infection; pathway generated using IPA. (JPG 2711 kb)
About this article
Cite this article
Braicu, C., Cojocneanu-Petric, R., Jurj, A. et al. Microarray based gene expression analysis of Sus Scrofa duodenum exposed to zearalenone: significance to human health. BMC Genomics 17, 646 (2016). https://doi.org/10.1186/s12864-016-2984-8
- Sus scrofa
- Extrapolation to human