Guo-Liang Zhang1,2Jun-Lin Song3Chuan-Liang Ji2Yu-Long Feng2Jie Yu2Charles M. Nyachoti4
Gong-She Yang1*- 1College of Animal Science and Technology, Northwest A&F University, Yangling, China
- 2National Engineering Research Center for Gelatin-based Traditional Chinese Medicine, Dong-E-E-Jiao Co., Ltd., Liaocheng, China
- 3Central Laboratory, Qingdao Agricultural University, Qingdao, China
- 4Department of Animal Science, University of Manitoba, Winnipeg, MB, Canada
Zearalenone (ZEA) is a natural contaminant existing in food and feed products that exhibits a negative effect on domestic animals’ reproduction. Donkeys possess high economic value in China and are at risk of exposure to ZEA. However, few information is available on ZEA-induced toxicity and no report on toxicity in donkeys can be found in scientific literature. We investigated the biological effects of ZEA exposure on donkey granulosa cells (dGCs) by using RNA-seq analysis. ZEA at 10 and 30 μM were administered to GCs within 72 h of in vitro culture. ZEA at 10 μM significantly altered the tumorigenesis associated genes in dGCs. Exposure to 10 and 30 μM ZEA treatment significantly reduced mRNA expression of PTEN, TGFβ, ATM, and CDK2 genes, particularly, the ZEA treatment significantly increased the expression of PI3K and AKT genes. Furthermore, immunofluorescence, RT-qPCR, and Western blot analysis verified the gene expression of ZEA-exposed GCs. Collectively, these results demonstrated the deleterious effect of ZEA exposure on the induction of ovarian cancer related genes via the PTEN/PI3K/AKT signaling pathway in dGCs in vitro.
Introduction
Zearalenone (ZEA) is a mycotoxin produced by various Fusarium fungi (Bennett and Klich, 2003) that infects grains and maize worldwide. Similar to aflatoxins (AFs), ZEA is one of the most important and widespread trichothecenes that cause extensive and recurring economic damage in cereal grains and animal feedstuffs (Escriva et al., 2015). In domestic animals, ZEA causes porcine ovarian atrophy (Vanyi et al., 1994) and equine follicular hematomas (Cortinovis et al., 2013) and exhibits significant genotoxic potential and induces DNA damage (Zhang et al., 2017a) in experimental animals (mice and rats). Owing to its estrogenic activity, ZEA could cause reproductive disorders in a wide variety of species-specific organs in animals (Poor et al., 2015). Both low and high concentrations of ZEA can cause abortion and reproductive failure in livestock (Osweiler et al., 1990; Dacasto et al., 1995; Zwierzchowski et al., 2005). Moreover, ZEA is a potential carcinogen with a possible correlation to xenoestrogens and breast cancer risk (Yu et al., 2004).
The effects of ZEA and some of its metabolites on estrous mares have only been reported in a few cases. Juhasz et al. (2001) reported the effects of 10 day low dose ZEA exposure on the reproductive system of ovulating mares fed with 7 mg purified ZEA. No effect was observed on the length of the inter-ovulatory intervals, follicular phases, and ovarian follicular activity (the maximum size and number of the ovulatory follicles, growth rate, and the initial increase in the number of large follicles) (Juhasz et al., 2001). Another study showed that feeding mares with oats that were naturally contaminated with ZEA (12 mg/kg) had no relevant effects on the cycle length and the release of reproductive hormones (Cortinovis et al., 2013). However, the authors stated that mares fed with the naturally contaminated oats had a high incidence of follicular hematomas, which did not occur in the control (Cortinovis et al., 2013). Similar to the findings in estrous pig (Zhang et al., 2017b), ZEA and its derivatives induced in vitro apoptosis of the granulosa cells (GCs), which were collected from estrous mares’ ovaries (Minervini et al., 2006). Hence, the authors suggested that ZEA could induce follicular atresia in domestic animals. These effects could be due to the direct interaction with ERs and 3-α/β-HSD enzymes present in the GCs and ovary, which are responsible for the synthesis or metabolism of the endogenous steroid hormones (Minervini et al., 2006).
The mechanism of ZEA toxicity is not fully understood, but ZEA is known to possess acute and chronic toxic effects in animals. Ovaries are female reproductive organs comprising follicles of varying sizes. The early stages of follicular growth depend on the development of the GCs and oocytes, which are in constant communication with each other. The development of one cell type influences the other’s compartment. During follicle development, GCs replicate, secrete hormones, and support the growth of the oocyte (Hamel et al., 2008). Previous investigations demonstrated that ZEA may alter GC’s function in swine (Ranzenigo et al., 2008; Zhang et al., 2017b). This study aims to evaluate the in vitro toxicity of 10 and 30 μM ZEA in donkey ovarian GCs through transcriptome analysis.
Materials and Methods
Reagents
Zearalenone was purchased from Sigma Company (St. Louis, MO, United States). Stock solutions of ZEA were fixed by dissolving ZEA in dimethyl sulfoxide (DMSO). DMSO (D12345), fetal bovine serum (FBS; 10100147), M-199 medium (11150-059), penicillin and streptomycin were procured from Gibco Company (Carlsbad, CA, United States).
Animals
The mature donkeys’ ovaries used in the experiment were obtained from the Dong E Donkey Production Company (Qingdao, Shandong, China). The ovaries were collected from the slaughterhouse of the company and maintained at 32–35°C for the isolation of GCs. All procedures of animal handling in this study were reviewed and approved by the Ethical Committee (Agreement No. 2017-18) of Qingdao Agricultural University.
Isolation and Culture of dGCs
Donkey GCs (dGCs) were aspirated from the antral follicles using a 10 ml syringe (Zhu et al., 2012). Standing for 15–18 min, the dGCs were centrifuged at 300 g for 5 min in accordance with previous report (Qin et al., 2015). Then the dGCs were cultured in DMEM medium (HyClone, SH30022.01, Beijing, China) supplement with 10% FBS (10099-141, Gibco, Australia) and 1% penicillin-streptomycin (Hyclone, SV30010) in incubator with 5% CO2 at 37°C (Duda et al., 2011).
The primary GCs were passaged after culture 48–36 h. To avoid the stress of passage response, drug exposure was performed until 12 h later. The GCs were inoculation in a 6 cm petri dish (Corning, 430166, United States) at a density of 1 × 106 cells. ZEA was added to the cultured medium at final concentrations of the 10 or 30 μM, then the cells were incubated with ZEA for 72 h. The control and 10 μM ZEA group added the same dose of DMSO to 30 μM ZEA group for accuracy.
Immunofluorescence and Cell Counting
The GCs were collected and fixed in 4% paraformaldehyde for 2 h, then heated at 42°C for 2 h, finally attached onto a poly-lysine coated slide. For immunofluorescence, the sections were blocked with the BDT (10% goat serum in TBS, 3% BSA) for 35 min, and then incubated overnight with primary antibodies at 4°C (Supplementary Table S1). The sections were then incubated with Cy3/FITC-conjugated goat anti-rabbit secondary antibody (Beyotime, A0208, Nantong, China) at a dilution of 1:50 at 37°C for 1.5 h. Finally, the sections were incubated with Hoechst33342 (Beyotime, C1022) to visualize nuclei of GCs for 3 min at room temperature. The immunosignal was detected using a fluorescence microscope (Olympus, XB51, Japan), the images were captured and analyzed in accordance with cellSens Standard.
TUNEL Staining
TUNEL BrightRed Apoptosis Detection Kit was utilized to evaluate the GCs apoptosis (Vazyme, A11302, Nanjing, China). Briefly, the GCs were fixed with 4% paraformaldehyde for 2 h after ZEA treatment. Observation under fluorescence microscopy was carried out after TUNEL reaction following manufacturer’s recommendations [60 min at 37°C without light, and DNA staining with Hoechst 33342 (Beyotime, C1022) incubation at room temperature for 3 min]. The TUNEL positive cells were detected under the fluorescence microscope (Olympus, XB51). To analyze TUNEL positive cell ratio, more than 2,000 cells were counted in each group, and each group included three biological replicates at least.
RNA Extraction, Reverse Transcription, and RNA-seq
Total RNA was extracted from the cultured GCs by the use of an RNAprep pure MicroKit in line with the manufacturer’s instruction (Aidlab, RN07, Beijing, China). And the first-strand cDNA acquisition was utilized a cDNA Synthesis Kit (TransGen, AT311-03, Beijing, China) with reference to previous study (Pan et al., 2011). The reaction program was setup at: 42°C for 15 min, 85°C for 30 s, and 4°C for cooling. Then, RNA sequencing was performed with Hiseq 4000 platform by Novogene (Beijing, China). Three biological replicates were manipulated in each group. Primary RNA-seq data were uploaded to the SRA repository. The SRA index number was SRP139976. The RNA-seq data matrix used for DESeq2 analysis were uploaded to the GEO repository and the GEO accession number was GSE1169501.
Data Preprocessing and Identification of Differentially Expressed Genes
R Bioconductor/DESeq2 package was applied to analyze the dGCs groups (control, 10 and 30 μM ZEA) to identify the difference of gene expression. DESeq2, the negative binomial was applied as the reference distribution and taken own normalization approach for raw counts in differential expression analysis. In other methods to avoid possible biases, the data of differential expression analysis has to be previously normalized (Benjamini et al., 2001; Luo and Brouwer, 2013). Because the design of sequencing contains biological replicates in each group, log2|fold change| is not setup as the filter condition. So padj < 0.01 was really considered as statistical significance.
KEGG Enrichment Analysis
R Bioconductor/clusterProfiler package was applied to analyze functional profiles (KEGG) of differentially expressed genes (DEGs) (Yu et al., 2012). Furthermore, R Bioconductor/Pathview package was applied to visualize KEGG enrichment results. According to the DEGs log2|fold change| value shows enrichment signaling pathway active status. The p-value adjusted used the Benjamini & Hochberg (BH) method (Luo and Brouwer, 2013). Still, the padj < 0.05 was considered as statistical significance.
Gene Set Enrichment Analysis
Gene set enrichment analysis (GSEA) does not need to specify a clearly differential expression of gene threshold, the algorithm performs GSEA based on overall trend of the gene raw read count. The file has to be made suitable for GSEA software analysis with R script (Subramanian et al., 2005). The data for GSEA analysis contained nine samples within three groups. The GSEA software option “Collapse dataset to gene symbols” parameters was “false” and option “Permutation type” parameters were “gene set.” Other parameters were set as default as software manual references. The GSEA gene set with FDR q-value <0.05 were defined as significant difference.
Protein–Protein Interaction Network Construction and Modules Mining
Search Tool for the Retrieval of Interacting Genes/Proteins2 (STRING) is a database of protein–protein interaction (PPI). This database contains the direct and physically related interactions between known and predicted protein and genes. The R Bioconductor/STRINGdb was applied for PPI of interested DEGs (Franceschini et al., 2013). The Cytoscape software was applied to visualize PPI results (Shannon et al., 2003).
Quantitative Real-Time PCR
Total RNA and cDNA reverse transcription was extracted as described above. Quantitative real-time PCR (RT-qPCR) was performed using Light-Cycler® 480 SYBR Green I Master Kit (Roche, Germany) on LightCycler 480 real-time PCR instrument. The RT-qPCR reaction was set as: 10 min at 95°C, followed by 45 cycles of 95°C for 10 s, 60°C for 30 s, 72°C for 30 s, and cooling step at 4°C. Quantitative RT-PCR primers are listed in Supplementary Table S2 and the primer efficiency in the Supplementary Figure S1. The standard curve gene GAPDH in dGCs was used as the reference to normalize the related genes’ mRNA expression. The difference multiples = 2-ΔΔCT method was employed for relative quantification of PCR. Each gene was expressed as the mean ± standard deviation (SD), which was calculated from independent biological replicates at least three times.
Western Blotting
Protein lysates isolated from dGCs were used for western blotting according to the standard protocol (Zhang et al., 2010; Chao et al., 2012). Proteins from each ZEA treatment were separated by 10% SDS-PAGE, then transferred onto the PVDF membranes. After, the membranes were incubated with primary antibodies (Supplementary Table S1) at 4°C overnight. Then the membranes were incubated with secondary antibodies (Beyotime, A0208) at 37°C for 2 h in TBST after rinsing three times with TBST. The final detection of related genes was carried out by AlphaImager® HP (ProteinSimple, 92-13824-00, United States). The band intensity was quantified using GAPDH as internal control and measured with IPWIN software.
Statistical Methods
Data are presented by mean ± SD. Different effects between the control and ZEA treatment groups of donkey was statistically determined by One-way ANOVA for multiple comparisons. All analyses were conducted using Graph-Pad Prism analysis software (San Diego, CA, United States). All experiments were repeated at least three times unless otherwise noted. Results were considered statistically significant at p-value <0.05.
Results
The Apoptosis and Apoptosis-Related Gene Expression of dGCs Exposed to ZEA
Donkey granulosa cells were cultured in vitro and exposed to 10 or 30 μM ZEA for 72 h (Figure 1A). The percentages of TUNEL positive dGCs significantly increased as a result of exposure of ZEA (10 μM: 36.04 ± 1.52%; 30 μM: 44.30 ± 1.33%) compared to that of the control (0 μM: 9.83 ± 0.21%; P < 0.01; Figure 1B). As shown in Figure 1C, the ratios of BAX/BCL2 mRNA expression significantly increased in 10 and 30 μM ZEA exposed dGCs.
FIGURE 1. Zearalenone exposure increasing apoptosis and inducing the apoptosis-related gene mRNA abundance in cultured GCs. GCs were stained blue with Hoechst33342 Solution. TUNEL assay was performed using immunostaining. (A) Immunofluorescent staining of TUNEL and Hoechst33342 of GCs. Bar indicates 50 μm. (B) The percentages of TUNEL positive GCs. (C) The mRNA levels of BAX/BCL2 in cultured GCs exposed to ZEA, and the control determined by Q-PCR. The results are presented as mean ± SD. All experiments were repeated at least three times. ∗P < 0.05; ∗∗P < 0.01.
We performed RNA-seq to verify the effect of ZEA exposure on dGCs. Based on the criterion FDR <0.05, 14,506 DEGs was observed between the control and the ZEA-treated dGCs (Figure 2A). A total of 7,253 and 6,984 DEGs were noted between the control and 10 and 30 μM ZEA treatment groups, respectively. A total of 269 DEGs were observed between the 10 and 30 μM ZEA treatment groups. In this study, the DEGs between the control and ZEA treatment groups and among the treatment groups were obtained from: 0 μM vs. 10 or 30 μM and 10 vs. 30 μM ZEA group, respectively (Figure 3A).
FIGURE 2. The GCs gene expression profiling after ZEA treatment. (A) Scatterplot of gene expression after ZEA treatment. Control group vs. 10 μM ZEA-treatment group, control group vs. 30 μM ZEA-treatment group, 10 μM ZEA-treatment group vs. 30 μM ZEA-treatment group. Red and blue plot represents genes expressed differently. (B) The nine samples shown in the 2D plane spanned by the first three principal components. (C) Heatmap indicate the group difference of DEGs in the 10 and 30 μM ZEA-treated groups compared with the control group, and the repeatability within each group. The results are presented as mean ± SD. All experiments were repeated at three times.
FIGURE 3. Gene expression of RNA-seq analysis in granulosa cells between control and ZEA-treated groups. (A) Venn diagram showed that the different expression of 14,506 genes in control and ZEA-treated groups. (B) Dispersion plot. The dispersion estimate plot shows the gene-wise estimates (black), the fitted values (red), and the final maximum a posteriori estimates used in testing (blue). (C) KEGG biological processes involving DEGs in ZEA-treatment GCs.
We obtained three RNA-seq replicates for the 0, 10, and 30 μM ZEA-treated dGCs. The variation of the replicates from the control and ZEA treatment groups are shown in Figure 2B. A heat map was drawn from the DEG results (Figure 2C). In this study, we chose the DEGs from the 0 μM vs. 10 and 30 μM groups.
To explore the potential mechanism of ZEA exposure in dGCs, the STRING database was applied in annotating functional interactions DEGs between the control and ZEA-treatment group. Also PPI networks was visualized by Cytoscape (Supplementary Figures S2A,B). The center nodes of the networks, some interesting genes, were observed, such as PTEN, AKT, ATM, TGFβ, PI3K, CCND2, HEY2, CDK2, and FDX1. In addition, the cBioPortal was used to provide visualization, analysis of DEGs related to the ovarian cancer from large-scale cancer genomics data sets (Supplementary Figures S2C, S3).
DEGs Involved in KEGG Pathways
In order to attain functional insights of DEGs, the R package of clusterProfiler was carried out to establish extremely affected KEGG pathways. The Venn diagram were constructed by the results of DEGs (Figure 3A). A total of 5,863 DEGs containing in control and 30 μM ZEA groups were quantified (Figure 3B). The DEGs were significantly enriched in the PI3K-AKT signaling pathway (Count = 125, padj = 0.0049), endocytosis (Count = 91, padj = 0.049), regulation of actin cytoskeleton (Count = 90, padj = 0.0001), and proteoglycans in cancer (Count = 83, padj = 0.0009) in ZEA treated dGCs (Figure 3C).
The GSEA analysis showed that OVARIAN_CANCER_LMP and TNF_SIGNALING_VIA_NFKB gene set was enriched by 10 μM ZEA vs. the control of dGCs (Figures 4A,B), while CANCER_HEAD_AND_NECK_VS_CERVICAL and BREAST_CANCER_16Q24_AMPLICON gene set was enriched by the 30 μM ZEA vs. the control in dGCs (Figures 4C,D), respectively. The GSEA enrichment plot revealed a regulated tendency for most of enriched genes in ZEA-treatment compared with the control. The GSEA analysis also verified the KEGG analysis.
FIGURE 4. Differentially expressed genes in GSEA analysis. (A,B) GSEA for GCs treated with 10 μM ZEA. (C,D) GSEA for GCs treated with 30 μM ZEA.
In addition, six genes were identified by performing the KEGG and GSEA pathway analysis in ZEA treatment, the PI3K and AKT genes were involved in PI3K-AKT signaling pathway of dGCs upregulating. DGCs in ZEA treatment, another four genes involved in anti-oncogene process downregulating were also identified, such as PTEN, TGFβ, CDK2, and ATM. Finally, CCND2, CDK6, TNF, and TP53 genes in treatment dGCs, were identified involving in cancer process either upregulation or downregulation.
Specific Impact of ZEA Exposure to dGCs
It is hypothesis that 10 and 30 μM ZEA exposure may affect the proliferation and carcinogenesis of dGCs (Figure 4). RT-qPCR was conducted to verify the hypothesis by evaluating the expression of different transcript together with enzymes in the pathway of dGCs between the control and ZEA treatments. RT-qPCR analyses indicated that 10 and 30 μM ZEA exposure significantly downregulated mRNA abundance of PTEN while upregulated PI3K and AKT genes in dGCs (P < 0.05 or P < 0.01; Figure 5A). As illustrated in Figure 5B, ZEA-treated dGCs exhibited lower protein levels of PTEN but higher protein levels of PI3K and AKT compared with that of the control dGCs (P < 0.05 or P < 0.01). Moreover, 10 and 30 μM ZEA exposure significantly down-regulated the mRNA abundance and protein levels of CDK2, TGFβ, and ATM genes than the control dGCs (P < 0.05 or P < 0.01; Figures 6A,B). Interestingly, the ZEA treatments significantly decreased the number of PTEN and TGFβ but significantly increased the number of PI3K immunofluorescence positive dGCs, in comparison to the control (Figures 7A,B).
FIGURE 5. Zearalenone exposure affecting mRNA and protein abundance of PTEN-related genes in cultured donkey GCs. (A) Quantitative RT-PCR for PTEN, AKT, and PI3K transcription factors. The mRNA levels of all genes were normalized to GCs GAPDH gene. (B) Protein levels of PTEN/GAPDH, AKT/GAPDH, and PI3K/GAPDH by Western blotting. The protein levels were normalized to GAPDH. The results are presented as mean ± SD. All experiments were repeated at least three times. ∗P < 0.05; ∗∗P < 0.01.
FIGURE 6. Zearalenone exposure affecting mRNA and protein abundance of tumorigenesis related genes in cultured GCs. (A) Quantitative RT-PCR for CDK2, TGFβ, and ATM transcription factors. The mRNA levels of all genes were normalized to GCs GAPDH gene. (B) Protein levels of CDK2/GAPDH TGFβ/GAPDH, and ATM/GAPDH by Western blotting. The protein levels were normalized to GAPDH. The results are presented as mean ± SD. All experiments were repeated at least three times. ∗P < 0.05; ∗∗P < 0.01.
FIGURE 7. Immunofluorescence assay examining the expression of phosphor-PTEN, TGFβ (A), and PI3K (B) proteins. The fluorescence intensity and percentages of positive cells were analyzed, respectively. Bar indicates 50 μm. Data are presented as means ± SD. ∗P < 0.05; ∗∗P < 0.01.
Discussion
Numerous studies have proven ZEA’s cytotoxic effect on the reproductive (Kiang et al., 1978; Mehmood et al., 2000; Nikov et al., 2000), immune (Abbes et al., 2006a,b; Luongo et al., 2006), and endocrine (Mueller et al., 2004) systems, as well as on heredity (Kouadio et al., 2005). ZEA indirectly affects mammalian fertility by impairing the formation of primordial follicle (Zhang et al., 2017a), thereby causing changes in gene expression; it also induces DNA damage in ovarian GCs (Zhang et al., 2017b). Not much information is available on ZEA-induced alteration in dGC proliferation and development and impairment of mammalian fertility. This research is the first to describe the GC transcriptomes of donkeys.
The RNA-seq results showed that exposure to 10 μM ZEA significantly transformed the mRNA expression of thousands of genes in dGCs. In particular, PTEN genes were downregulated, suggesting that ZEA exposure is involved with the cancer-related PTEN signaling pathway. The tumor suppressor gene, PTEN, plays an essential role in cell growth, survival, and tumor formation. Decrease in PTEN expression, either at the protein or mRNA level, has been associated with many primary malignancies, including ovarian cancer (Kechagioglou et al., 2014). Hence, the downregulation of PTEN is a hallmark of tumors. The PI3K-AKT pathway is one of the two major signaling pathways that have been identified as important in cancer. Through phosphorylation, PI3K-AKT inhibits the activity of proapoptotic members while activating anti-apoptotic members. The reduction in PTEN expression indirectly stimulates PI3K-AKT activity, thereby contributing to oncogenesis in mammals (Zhang et al., 2017c). Our studies showed that exposure to 10 μM ZEA significantly upregulated the PI3K-AKT gene expression in the dGCs, implying that exposure to low concentrations of ZEA (10 μM) might increase the donkey’s risk for ovarian cancer via the PTEN/PI3K/AKT pathway by suppressing the expression of antitumor genes or by activating the expression of cancer-causing genes (Zhang et al., 2017c). Moreover, it was reported that decreased TGFβ-mediated signaling might predispose an individual to develop cancer (Pasche, 2001; Galliher et al., 2006; Jin et al., 2017). Previous studies have assessed the association between TGFβ and risk for various forms of cancer, and several meta-analyses have demonstrated that TGFβ is associated with risk for ovarian cancer (Akhurst and Derynck, 2001; Nilsson and Skinner, 2002; de la Cruz-Merino et al., 2009). Interestingly, TGFβ expression in the ZEA treatment samples was significantly decreased, indicating an increased risk for ovarian cancer. The direct effect of TGFβ is not solely responsible for influencing tumor behavior (Hagedorn et al., 2001; Tuxhorn et al., 2002). Lack of TGFβ in fibroblasts can result in mammary gland tumor progression (Bhowmick et al., 2004; Matise et al., 2012).
Bioinformatics analysis confirmed that ataxia-telangiectasia mutated (ATM) and cyclin-dependent kinase 2 (CDK2) genes were affected by exposure to ZEA. A previous study showed that the ATM gene plays an essential role in DNA double-strand breaks (DSBs) repair (Kurz and Lees-Miller, 2004). These DSBs if left unchecked can result in the development of cancer (Weber and Ryan, 2015). Several studies showed that suppression of ATM is associated with a variety of tumors (Weber and Ryan, 2015). Here, we provide evidence that exposure to 10 μM ZEA significantly decreased ATM gene expression, suggesting a common mechanism of action in donkeys. Other studies revealed that cell cycle dysregulation, resulting in uncontrolled proliferation, is also a hallmark of cancer (Vijayaraghavan et al., 2017). The CDK family is composed of proteins associated with cell cycle regulation and is frequently mutated or overexpressed in ovarian cancer. Deregulation of the CDK2/4/6 signaling pathway is among the most common aberrations found in ovarian cancer (D’Andrilli et al., 2004). Interestingly, exposure to 10 μM ZEA significantly downregulated the expression of CDK2 genes, which regulate the cell cycle and are involved in ovarian cancer. Finally, exposure to ZEA increased the dGCs’ apoptosis rate and elevated the expression of the BAX/BCL2 genes.
Conclusion
This research is the first study to investigate ZEA-induced impairment of dGCs. ZEA (10 or 30 μM), a potentially carcinogenic substance, can directly cause tumorigenesis and in vitro apoptosis of dGCs. This study developed an innovative, integrated, and low-cost approach to study GC exposure to ZEA. ZEA disrupts the endocrine and reproductive performance of domestic animals, and this study may help elucidate the mechanism of ZEA toxicity.
Ethics Statement
All donkeys were treated humanely during slaughter (No. 11002009000012, production license number: SCXK (Qingdao): 2012–0766, Dongeejiao, Qingdao, China). The procedures of animal handling in this study were reviewed and approved by the Ethical Committee of Qingdao Agricultural University (Protocol No. 2017-018).
Author Contributions
G-LZ and G-SY designed the study and wrote the manuscript. J-LS, C-LJ, Y-LF, JY, and CN conducted the experiments and analyzed the data.
Funding
This work was supported by National Science and Technology Major Project of China (2016ZX08006003), National Key Research and Development Program of China (2016YFD0501207), Construction of E-Jiao Standardization (ZYBZH-Y-SD-31), and the Agriculture Breeds Engineering of Shandong Province of China (2017LZGC020).
Conflict of Interest Statement
C-LJ, Y-LF, and JY were employed by company Dong-E-E-Jiao Co., Ltd.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors are grateful to Dr. Paul W. Dyce for his suggestions on the manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2018.00293/full#supplementary-material
FIGURE S1 | The standard curve and primer efficiency of all Q-PCR primers. (A–F) represent the standard curve and primer efficiency of PTEN, AKT, PI3K, CDK2, TGFβ, and ATM genes.
FIGURE S2 | Protein–protein interaction (PPI) network based on the STRING database to annotate functional interactions between DEGs in control and 10 μM ZEA-treatment groups (A), in control and 30 μM ZEA-treatment groups (B). Cancer type detailed summary of DEGs related to the ovarian cancer (C). The node degree ≥20 was selected as the threshold.
FIGURE S3 | Visualization of DEGs related to the ovarian cancer from large-scale cancer genomics data sets. The node degree ≥20 was selected as the threshold.
TABLE S1 | Primary antibodies.
TABLE S2 | Primers used for quantitative-PCR.
Footnotes
References
Abbes, S., Ouanes, Z., ben Salah-Abbes, J., Houas, Z., Oueslati, R., Bacha, H., et al. (2006a). The protective effect of hydrated sodium calcium aluminosilicate against haematological, biochemical and pathological changes induced by zearalenone in mice. Toxicon 47, 567–574. doi: 10.1016/j.toxicon.2006.01.016
Abbes, S., Salah-Abbes, J. B., Ouanes, Z., Houas, Z., Othman, O., Bacha, H., et al. (2006b). Preventive role of phyllosilicate clay on the immunological and biochemical toxicity of zearalenone in Balb/c mice. Int. Immunopharmacol. 6, 1251–1258. doi: 10.1016/j.intimp.2006.03.012
Akhurst, R. J., and Derynck, R. (2001). TGF-beta signaling in cancer–a double-edged sword. Trends Cell Biol. 11, S44–S51. doi: 10.1016/S0962-8924(01)02130-4
Benjamini, Y., Drai, D., Elmer, G., Kafkafi, N., and Golani, I. (2001). Controlling the false discovery rate in behavior genetics research. Behav. Brain Res. 125, 279–284. doi: 10.1016/S0166-4328(01)00297-2
Bennett, J. W., and Klich, M. (2003). Mycotoxins. Clin. Microbiol. Rev. 16, 497–516. doi: 10.1128/CMR.16.3.497-516.2003
Bhowmick, N. A., Chytil, A., Plieth, D., Gorska, A. E., Dumont, N., Shappell, S., et al. (2004). TGF-beta signaling in fibroblasts modulates the oncogenic potential of adjacent epithelia. Science 303, 848–851. doi: 10.1126/science.1090922
Chao, H. H., Zhang, X. F., Chen, B., Pan, B., Zhang, L. J., Li, L., et al. (2012). Bisphenol a exposure modifies methylation of imprinted genes in mouse oocytes via the estrogen receptor signaling pathway. Histochem. Cell Biol. 137, 249–259. doi: 10.1007/s00418-011-0894-z
Cortinovis, C., Pizzo, F., Spicer, L. J., and Caloni, F. (2013). Fusarium mycotoxins: effects on reproductive function in domestic animals–a review. Theriogenology 80, 557–564. doi: 10.1016/j.theriogenology.2013.06.018
Dacasto, M., Rolando, P., Nachtmann, C., Ceppa, L., and Nebbia, C. (1995). Zearalenone mycotoxicosis in piglets suckling sows fed contaminated grain. Vet. Hum. Toxicol. 37, 359–361.
D’Andrilli, G., Kumar, C., Scambia, G., and Giordano, A. (2004). Cell cycle genes in ovarian cancer: steps toward earlier diagnosis and novel therapies. Clin. Cancer Res. 10, 8132–8141. doi: 10.1158/1078-0432.CCR-04-0886
de la Cruz-Merino, L., Henao-Carrasco, F., Garcia-Manrique, T., Fernandez-Salguero, P. M., and Codes-Manuel de Villena, M. (2009). Role of transforming growth factor beta in cancer microenvironment. Clin. Transl. Oncol. 11, 715–720. doi: 10.1007/s12094-009-0433-8
Duda, M., Knet, M., Tabarowski, Z., and Slomczynska, M. (2011). Luteal macrophage conditioned medium affects steroidogenesis in porcine granulosa cells. Reprod. Biol. 11, 117–134. doi: 10.1016/S1642-431X(12)60049-7
Escriva, L., Font, G., and Manyes, L. (2015). In vivo toxicity studies of Fusarium mycotoxins in the last decade: a review. Food Chem. Toxicol. 78, 185–206. doi: 10.1016/j.fct.2015.02.005
Franceschini, A., Szklarczyk, D., Frankild, S., Kuhn, M., Simonovic, M., Roth, A., et al. (2013). String v9.1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 41, D808–D815. doi: 10.1093/nar/gks1094
Galliher, A. J., Neil, J. R., and Schiemann, W. P. (2006). Role of transforming growth factor-beta in cancer progression. Future Oncol. 2, 743–763. doi: 10.2217/14796694.2.6.743
Hagedorn, H. G., Bachmeier, B. E., and Nerlich, A. G. (2001). Synthesis and degradation of basement membranes and extracellular matrix and their regulation by TGF-beta in invasive carcinomas (Review). Int. J. Oncol. 18, 669–681. doi: 10.3892/ijo.18.4.669
Hamel, M., Dufort, I., Robert, C., Gravel, C., Leveille, M. C., Leader, A., et al. (2008). Identification of differentially expressed markers in human follicular cells associated with competent oocytes. Hum. Reprod. 23, 1118–1127. doi: 10.1093/humrep/den048
Jin, Y., Chen, W., Yang, H., Yan, Z., Lai, Z., Feng, J., et al. (2017). Scutellaria barbata D. Don inhibits migration and invasion of colorectal cancer cells via suppression 435 of PI3K/AKT and TGF-beta/Smad signaling pathways. Exp. Ther. Med. 14, 5527–5534. doi: 10.3892/etm.2017.5242
Juhasz, J., Nagy, P., Kulcsar, M., Szigeti, G., Reiczigel, J., and Huszenicza, G. (2001). Effect of low-dose zearalenone exposure on luteal function, follicular activity and uterine oedema in cycling mares. Acta Vet. Hung. 49, 211–222. doi: 10.1556/004.49.2001.2.11
Kechagioglou, P., Papi, R. M., Provatopoulou, X., Kalogera, E., Papadimitriou, E., Grigoropoulos, P., et al. (2014). Tumor suppressor PTEN in breast cancer: heterozygosity, mutations and protein expression. Anticancer Res. 34, 1387–1400.
Kiang, D. T., Kennedy, B. J., Pathre, S. V., and Mirocha, C. J. (1978). Binding characteristics of zearalenone analogs to estrogen receptors. Cancer Res. 38, 3611–3615.
Kouadio, J. H., Mobio, T. A., Baudrimont, I., Moukha, S., Dano, S. D., and Creppy, E. E. (2005). Comparative study of cytotoxicity and oxidative stress induced by deoxynivalenol, zearalenone or fumonisin B1 in human intestinal cell line Caco-2. Toxicology 213, 56–65. doi: 10.1016/j.tox.2005.05.010
Kurz, E. U., and Lees-Miller, S. P. (2004). DNA damage-induced activation of ATM and ATM-dependent signaling pathways. DNA Repair 3, 889–900. doi: 10.1016/j.dnarep.2004.03.029
Luo, W., and Brouwer, C. (2013). Pathview: an R/Bioconductor package for pathway-based data integration and visualization. Bioinformatics 29, 1830–1831. doi: 10.1093/bioinformatics/btt285
Luongo, D., Severino, L., Bergamo, P., De Luna, R., Lucisano, A., and Rossi, M. (2006). Interactive effects of fumonisin B1 and alpha-zearalenol on proliferation and cytokine expression in jurkat T cells. Toxicol. In Vitro 20, 1403–1410. doi: 10.1016/j.tiv.2006.06.006
Matise, L. A., Palmer, T. D., Ashby, W. J., Nashabi, A., Chytil, A., Aakre, M., et al. (2012). Lack of transforming growth factor-beta signaling promotes collective cancer cell invasion through tumor-stromal crosstalk. Breast Cancer Res. 14:R98. doi: 10.1186/bcr3217
Mehmood, Z., Smith, A. G., Tucker, M. J., Chuzel, F., and Carmichael, N. G. (2000). The development of methods for assessing the in vivo oestrogen-like effects of xenobiotics in CD-1 mice. Food Chem. Toxicol. 38, 493–501. doi: 10.1016/S0278-6915(00)00022-3
Minervini, F., Giannoccaro, A., Fornelli, F., Dell’Aquila, M. E., Minoia, P., and Visconti, A. (2006). Influence of mycotoxin zearalenone and its derivatives (alpha and beta zearalenol) on apoptosis and proliferation of cultured granulosa cells from equine ovaries. Reprod. Biol. Endocrinol. 4:62. doi: 10.1186/1477-7827-4-62
Mueller, S. O., Simon, S., Chae, K., Metzler, M., and Korach, K. S. (2004). Phytoestrogens and their human metabolites show distinct agonistic and antagonistic properties on estrogen receptor alpha (ERalpha) and ERbeta in human cells. Toxicol. Sci. 80, 14–25. doi: 10.1093/toxsci/kfh147
Nikov, G. N., Hopkins, N. E., Boue, S., and Alworth, W. L. (2000). Interactions of dietary estrogens with human estrogen receptors and the effect on estrogen receptor-estrogen response element complex formation. Environ. Health Perspect. 108, 867–872. doi: 10.1289/ehp.00108867
Nilsson, E. E., and Skinner, M. K. (2002). Role of transforming growth factor beta in ovarian surface epithelium biology and ovarian cancer. Reprod. Biomed. Online 5, 254–258. doi: 10.1016/S1472-6483(10)61828-7
Osweiler, G. D., Stahr, H. M., and Beran, G. W. (1990). Relationship of mycotoxins to swine reproductive failure. J. Vet. Diagn. Invest. 2, 73–75. doi: 10.1177/104063879000200115
Pan, B., Chao, H., Chen, B., Zhang, L., Li, L., Sun, X., et al. (2011). DNA methylation of germ-cell-specific basic helix-loop-helix (HLH) transcription factors. Mol. Hum. Reprod. 17, 550–561. doi: 10.1093/molehr/gar017
Pasche, B. (2001). Role of transforming growth factor beta in cancer. J. Cell. Physiol. 186, 153–168. doi: 10.1002/1097-4652(200002)186:2<153::AID-JCP1016>3.0.CO;2-J
Poor, M., Kunsagi-Mate, S., Sali, N., Koszegi, T., Szente, L., and Peles-Lemli, B. (2015). Interactions of zearalenone with native and chemically modified cyclodextrins and their potential utilization. J. Photochem. Photobiol. B 151, 63–68. doi: 10.1016/j.jphotobiol.2015.07.009
Qin, X., Cao, M., Lai, F., Yang, F., Ge, W., Zhang, X., et al. (2015). Oxidative stress induced by zearalenone in porcine granulosa cells and its rescue by curcumin in vitro. PLoS One 10:e0127551. doi: 10.1371/journal.pone.0127551
Ranzenigo, G., Caloni, F., Cremonesi, F., Aad, P. Y., and Spicer, L. J. (2008). Effects of Fusarium mycotoxins on steroid production by porcine granulosa cells. Anim. Reprod. Sci. 107, 115–130. doi: 10.1016/j.anireprosci.2007.06.023
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504. doi: 10.1101/gr.1239303
Subramanian, A., Tamayo, P., Mootha, V. K., Mukherjee, S., Ebert, B. L., Gillette, M. A., et al. (2005). Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. U.S.A. 102, 15545–15550. doi: 10.1073/pnas.0506580102
Tuxhorn, J. A., McAlhany, S. J., Yang, F., Dang, T. D., and Rowley, D. R. (2002). Inhibition of transforming growth factor-beta activity decreases angiogenesis in a human prostate cancer-reactive stroma xenograft model. Cancer Res. 62, 6021–6025.
Vanyi, A., Bata, A., Glavits, R., and Kovacs, F. (1994). Perinatal oestrogen syndrome in swine. Acta Vet. Hung. 42, 433–446.
Vijayaraghavan, S., Moulder, S., Keyomarsi, K., and Layman, R. M. (2017). Inhibiting CDK in cancer therapy: current evidence and future directions. Target. Oncol. 13, 21–38. doi: 10.1007/s11523-017-0541-2
Weber, A. M., and Ryan, A. J. (2015). ATM and ATR as therapeutic targets in cancer. Pharmacol. Ther. 149, 124–138. doi: 10.1016/j.pharmthera.2014.12.001
Yu, G., Wang, L. G., Han, Y., and He, Q. Y. (2012). Clusterprofiler: an R package for comparing biological themes among gene clusters. OMICS 16, 284–287. doi: 10.1089/omi.2011.0118
Yu, Z., Hu, D., and Li, Y. (2004). Effects of zearalenone on mRNA expression and activity of cytochrome P450 1A1 and 1B1 in MCF-7 cells. Ecotoxicol. Environ. Saf. 58, 187–193. doi: 10.1016/j.ecoenv.2004.03.001
Zhang, G. L., Sun, X. F., Feng, Y. Z., Li, B., Li, Y. P., Yang, F., et al. (2017a). Zearalenone exposure impairs ovarian primordial follicle formation via down-regulation of Lhx8 expression in vitro. Toxicol. Appl. Pharmacol. 317, 33–40. doi: 10.1016/j.taap.2017.01.004
Zhang, G. L., Zhang, R. Q., Sun, X. F., Cheng, S. F., Wang, Y. F., Ji, C. L., et al. (2017b). RNA-seq based gene expression analysis of ovarian granulosa cells exposed to zearalenone in vitro: significance to steroidogenesis. Oncotarget 8, 64001–64014. doi: 10.18632/oncotarget.19699
Zhang, Z., Chen, Q., Zhang, J., Wang, Y., Hu, X., Yin, S., et al. (2017c). Associations of genetic polymorphisms in PTEN/AKT/mTOR signaling pathway genes with cancer risk: a meta-analysis in Asian population. Sci. Rep. 7:17844. doi: 10.1038/s41598-017-17250-z
Zhang, P., Chao, H., Sun, X., Li, L., Shi, Q., and Shen, W. (2010). Murine folliculogenesis in vitro is stage-specifically regulated by insulin via the Akt signaling pathway. Histochem. Cell Biol. 134, 75–82. doi: 10.1007/s00418-010-0708-8
Zhu, L., Yuan, H., Guo, C., Lu, Y., Deng, S., Yang, Y., et al. (2012). Zearalenone induces apoptosis and necrosis in porcine granulosa cells via a caspase-3- and caspase-9-dependent mitochondrial signaling pathway. J. Cell. Physiol. 227, 1814–1820. doi: 10.1002/jcp.22906
Keywords: donkey, granulosa cells, tumorigenesis, gene expression, RNA-seq
Citation: Zhang G-L, Song J-L, Ji C-L, Feng Y-L, Yu J, Nyachoti CM and Yang G-S (2018) Zearalenone Exposure Enhanced the Expression of Tumorigenesis Genes in Donkey Granulosa Cells via the PTEN/PI3K/AKT Signaling Pathway. Front. Genet. 9:293. doi: 10.3389/fgene.2018.00293
Received: 08 April 2018; Accepted: 13 July 2018;
Published: 31 July 2018.
Edited by:
Pierre R. Bushel, National Institute of Environmental Health Sciences (NIEHS), United StatesReviewed by:
Xu Wang, Auburn University, United StatesClarLynda Williams-DeVane, North Carolina Central University, United States
Copyright © 2018 Zhang, Song, Ji, Feng, Yu, Nyachoti and Yang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Gong-She Yang, MTUzMzUzMjkxOTFAMTYzLmNvbQ==