Skip to main content

Microarray based gene expression analysis of Sus Scrofa duodenum exposed to zearalenone: significance to human health



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 [35]. 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 [68]. 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 [1315]. 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 [1620].

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 [21], 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 [23], 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).

Fig. 1

qRT-PCR data validation of the microarray data

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.

Fig. 2

Protein expression quantification by ELISA at duodenum level for IL-1B, IL-8 and IL-4 for the control group and the group exposed to ZEA

Network analysis

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

Table 1 Top 4 canonical pathway targeted by ZEA at duodenum level

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.

Fig. 3

Gene network interaction generated using IPA, as response to ZEA exposure at duodenul level, as key central molecule is FOS, being relate to oxidative stress response and TGFB1 pathways. The altered expressed genes are displayed in red and green, based on fold change expression level, meanwhile the gene in grey are retrieved in literature to interact

Fig. 4

Gene network interaction generated using IPA, as response to ZEA exposure at duodenum level, as key central molecule is RB, being relate to alteration of the cell cycle or NFkB pathways. The altered expressed genes are displayed in red and green, based on fold change expression level, meanwhile the gene in grey are retrieved in literature to interact

Table 2 Top disease and function altered as response to ZEA exposure


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 [26]. 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 [27]. 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 [28], a mechanism that is altered by cytotoxic and carcinogenic agents [29], 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 [30]. The capacity to activate oxidative stress was proved also in an in vitro study [31]. It was reported that ZEA has pro-apoptotic effects, by increasing the Bax/Bcl-2 ration [32]. ZEN was observed to be involved in hepatic injury via Nrf2/ARE signaling pathway [33], mechanism altered also at intestinal level. Therefore, continuous exposure might be related to the activation of the tumor-promoting inflammatory mechanisms [30], 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 [34]. In the case of HCT116 cells it was observed that even at low concentration ZEA increase cell proliferation, colony formation and cell migration [35].

ZEA compounds were observed to exert estrogenic effects [7, 36], even at nM concentration [36], 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 [36]. ZEN could promote the autocrine action or expression of the ghrelin gene in piglet ovary [37].

Other studies confirmed the activation of MAPK and oxidative stress related genes as response to ZEA in MCF-7 breast cancer cells, [38] and macrophages [31], 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 [40]. 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 [36].

Alteration of Stathmin1 pathways are connected with invasion and migration capacity in digestive cancers [41]. The oncogenic role of the Stathmin1 pathways was proved to be related via miR-223 [42]. 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) [43]. 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 [44]. 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 [45].

The results remains to be further investigated particular on ZEA derivatives, a recent study present antitumoral effects of the semisynthetic zearalenone analogues [46]. Therefore the mechanisms of transformation and detoxification in parallel with the combinatorial effect need to be taken in account [47]. The procarcinogenic effect of ZEA was also studied by Ouanes et al., 2003 [48], and should take into account the accumulation in organism and the time required for metabolisation of this toxins [49]. Novel monitoring method based on urinary excretion of ZEA metabolites should be developed [50].

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 [52]. 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 [51], 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 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. [53].

The information related to homology annotation for all the transcripts sequence for S. scrofa and humans was used NCBI database ( 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://

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.

Table 3 qPCR primer sequences and characteristics

Cytokine concentration evaluation (ELISA)

Organ lysate was prepared as described by Taranu et al., (2014) [10] 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


  1. 1.

    Caldwell RW, Tuite J, Stob M, Baldwin R. Zearalenone production by Fusarium species. Appl Microbiol. 1970;20(1):31–4.

    CAS  PubMed  PubMed Central  Google Scholar 

  2. 2.

    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.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Yazar S, Omurtag GZ. Fumonisins, trichothecenes and zearalenone in cereals. Int J Mol Sci. 2008;9(11):2062–90.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  4. 4.

    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.

    CAS  PubMed  Google Scholar 

  5. 5.

    Iqbal SZ, Rabbani T, Asi MR, Jinap S. Assessment of aflatoxins, ochratoxin A and zearalenone in breakfast cereals. Food Chem. 2014;157:257–62.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    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.

    Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    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.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Diekman MA, Green ML. Mycotoxins and reproduction in domestic livestock. J Anim Sci. 1992;70(5):1615–27.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Hueza IM, Raspantini PC, Raspantini LE, Latorre AO, Gorniak SL. Zearalenone, an estrogenic mycotoxin, is an immunotoxic compound. Toxins. 2014;6(3):1080–95.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    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.

    Article  PubMed  Google Scholar 

  11. 11.

    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.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    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.

    Article  PubMed  PubMed Central  Google Scholar 

  13. 13.

    Tenaillon O, Skurnik D, Picard B, Denamur E. The population genetics of commensal Escherichia coli. Nat Rev Microbiol. 2010;8(3):207–17.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Sekirov I, Russell SL, Antunes LC, Finlay BB. Gut microbiota in health and disease. Physiol Rev. 2010;90(3):859–904.

    CAS  Article  Google Scholar 

  15. 15.

    Huffnagle G, Noverr MC. GI microbiota and regulation of the immune system. Preface. Adv Exp Med Biol. 2008;635:v–vi.

    PubMed  Google Scholar 

  16. 16.

    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.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    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.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    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.

    Article  Google Scholar 

  20. 20.

    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.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Bennett JW, Klich M. Mycotoxins. Clin Microbiol Rev. 2003;16(3):497–516.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    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.

    Google Scholar 

  24. 24.

    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.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    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.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Grenier B, Applegate TJ. Modulation of intestinal functions following mycotoxin ingestion: meta-analysis of published experiments in animals. Toxins. 2013;5(2):396–430.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    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.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    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.

    Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Chipman JK, Mally A, Edwards GO. Disruption of gap junctions in toxicity and carcinogenicity. Toxic Sci. 2003;71(2):146–53.

    CAS  Article  Google Scholar 

  30. 30.

    West AC, Jenkins BJ. Inflammatory and Non-Inflammatory Roles for Toll-Like Receptors in Gastrointestinal Cancer. Curr Pharm Des. 2015;21(21):2968–77.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    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.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    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.

  33. 33.

    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.

  34. 34.

    Grivennikov SI, Greten FR, Karin M. Immunity, Inflammation, and Cancer. Cell. 2010;140(6):883–99.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  35. 35.

    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.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    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.

    CAS  Article  PubMed  Google Scholar 

  37. 37.

    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.

    CAS  Article  PubMed  Google Scholar 

  38. 38.

    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.

    CAS  Article  PubMed  Google Scholar 

  39. 39.

    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.

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    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.

    CAS  Article  PubMed  Google Scholar 

  41. 41.

    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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  42. 42.

    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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  43. 43.

    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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  44. 44.

    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.

  45. 45.

    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.

    Article  PubMed  PubMed Central  Google Scholar 

  46. 46.

    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.

  47. 47.

    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.

  48. 48.

    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.

    CAS  Article  PubMed  Google Scholar 

  49. 49.

    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.

    CAS  Article  PubMed  Google Scholar 

  50. 50.

    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.

    CAS  Article  PubMed  Google Scholar 

  51. 51.

    Betts KS. Characterizing exposomes: tools for measuring personal environmental exposures. Environ Health Perspect. 2012;120(4):A158–163.

    Article  PubMed  PubMed Central  Google Scholar 

  52. 52.

    Buck Louis GM, Sundaram R: Exposome: Time for Transformative Research. Statistics in medicine 2012, 31(22):10.1002/sim.5496.

  53. 53.

    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.

    Article  PubMed  PubMed Central  Google Scholar 

  54. 54.

    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.

    CAS  Article  PubMed  Google Scholar 

  55. 55.

    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.

    Article  PubMed  Google Scholar 

  56. 56.

    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.

    CAS  Article  PubMed  Google Scholar 

Download references


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.

Authors’ contributions

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.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Author information



Corresponding authors

Correspondence to Cornelia Braicu or Ioana Berindan-Neagoe.

Additional files

Additional file 1: Figure S1.

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)

Additional file 2: Figure S2.

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)

Additional file 3: Figure S3.

Activation of colorectal carcinogenic mechanism as response to ZEA exposure; pathway generated using IPA. (JPG 5802 kb)

Additional file 4: Figure S4.

Alterations in GAP junction signaling with signification on the increasing the susceptibility to bacterial infection; pathway generated using IPA. (JPG 2711 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

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

Download citation


  • Zearalenone
  • Duodenum
  • Sus scrofa
  • Extrapolation to human
  • Health
  • Microarray