Abstract
Interferon regulatory factor 6 (IRF6) acts as a tumor suppressor and controls cell differentiation in ectodermal and craniofacial tissues by regulating expression of target genes. However, its function in gastric cancer (GC) remains unknown to date. In this study, we found that the IRF6 expression was significantly downregulated in GC. And the decreased expression of IRF6 was clinically correlated with poor prognosis of GC. Moreover, loss-of-function and gain-of-function studies showed that IRF6 was negatively regulated by ZEB1 but positively regulated by ELF3. Additionally, transcription factor ZEB1 and ELF3 could directly bind on IRF6 promoter, which suggested that transcription factor IRF6 is transcriptionally regulated by ZEB1 and ELF3. Nevertheless, we found that IRF6 expression was negatively related to its promoter methylation in TCGA stomach cancer cohorts. The downregulation of IRF6 in GC might be due to the overexpression of ZEB1 and the DNA methylation of IRF6 promoter.
Introduction
Gastric cancer (GC) is the third leading cause of cancer-related death worldwide, especially in East Asia (1). Despite a decrease in its incidence in some regions of the world, the 5-years overall survival rate remains unsatisfactory because most of GC patients are diagnosed in an advanced stage, with a poor prognosis and limited treatment options (2, 3). Although there is a great advancement on the gastric carcinogenesis, the molecular mechanisms underlying GC progression remains unclear (4). Hence, better understanding of the GC progression is essential to identify new effective diagnostic markers and novel effective therapies for GC patients.
Interferon regulatory factor 6 (IRF6) belongs to a family of nine transcription factors that share a highly conserved helix–turn–helix DNA-binding domain and a less conserved protein-binding domain (5). Unlike other IRF family members, IRF6 is not involved in IFN gene regulation, but instead have an essential role in skin development and keratinocyte differentiation (6, 7). Recent studies have proved that IRF6, regulated by TP63, plays a tumor suppressor role in squamous cell carcinomas through a Notch-dependent mechanism (8–10). However, the IRF6 gene expression regulation in gastrointestinal cancer types are not yet reported.
Previous studies have reported that IRF6 was downregulated during EMT process of breast cancer and prostate cancer (11–13). ZEB1 and ELF3 (E74 like ETS transcription factor 3) are two of transcription factors that involved in epithelial-mesenchymal transition (EMT) process. Transcription factor ZEB1 is known to be a master regulator of EMT process in varied cancer types (14, 15). Unlike the promotion of EMT by ZEB1, ELF3 is proved to be a negative regulator of EMT in ovarian cancer cells (16). Besides, it's recently reported that ELF3 functions as an antagonist of oncogenic-signaling-induced expression of EMT-TF ZEB1 in colorectal cancer (17). These results suggested that ELF3 and ZEB1 might play opposite roles in the EMT process of tumor cells.
In this study, we identified that the expression level of transcription factor IRF6 protein and mRNA was significantly downregulated in GC. Besides, the lower expression of IRF6 predicts poorer prognosis of GC. In addition, our data highlights that ZEB1 and ELF3 are two transcriptional regulators of IRF6 in GC; IRF6 is negatively regulated by ZEB1 but positively regulated by ELF3 in GC.
Materials and Methods
Expression Data Analysis
Six public gastric cancer microarray gene profiling datasets (GSE26942, GSE35809, GSE54129, GSE62254, GSE63089, and GSE79973) were downloaded from the Gene Expression Omnibus (GEO) in the NCBI web server. The clinical information of GSE62254 gastric cancer cohort was download from Cristescu et al. (18). The expression data of normal stomach tissue and stomach cancer tissues was obtained from GTEx and TGCA, respectively by using GEPIA (http://gepia.cancer-pku.cn/index.html) and UCSC (http://genome.ucsc.edu/). The chip-seq data of ELF3 and ZEB1 download from NCBI was analyzed by Cistrome (http://cistrome.org/db/#/). The chip-seq data of Dnase I, H3K27Ac, EP300, SRF, POLR2A, MED1, CREB1 and BRD4 were obtained from http://cistrome.org/browser/?genome=hg38wugb&datahub=http://dc2.cistrome.org/data5/browser/1535255763.json&gftk=refGene,full&coordinate = chr7:27066839-27266927. The possible binding sites of ELF3 and ZEB1 in the 2000-length IRF6 promoter was predicted by JASPAR (http://jaspar.genereg.net/). The correlation of IRF6 promoter methylation and IRF6 expression was analyzed by MEXPRESS (https://mexpress.be/) and LinkedOmics (http://www.linkedomics.org/login.php).
Cell Culture
The human gastric cancer cell line (N87, BGC823, AGS, SGC7901, MGC803, and HGC27) and a normal gastric epithelium cell line (GES-1) were purchased from the Shanghai Cell Bank of Chinese Academy of Sciences (Shanghai, China). All the GC cell lines and the GES-1 Cell line were cultured in DMEM medium containing 10% fetal bovine serum (FBS), 100 U/mL penicillin,100 U/mL streptomycin and 0.03% glutamine at 37°C in 5% CO2.
Cell Transfection
The siRNAs (si-ELF3#1: 5′-GCUGCAACCUGUGAGAUUA-3′, si-ELF3#2: 5′-CC-UCUGCAAUUGUGCCCUU-3′, si-ELF3#3: 5′-CCAUGAGGUACUACUACAA-3′; si-ZEB1: 5′-UGAUCAGCCUCAAUCUGCA-3′, si-NC: 5′-UUCUCCGAACGUGUCACG U-3′) used to knock down ELF3 and ZEB1 expression in GC cell lines were designed and synthesized by Genepharma (Shanghai, China). The ORF region of ZEB1 cDNA was cloned into pcDNA3.1 (+). Gastric cancer cells were grown in 6-well plates and transfected by Lipofectamine 2000 (Invitrogen) according to the manufacturer's instructions. At 48 h post-transfection, cells were harvested for qPCR analysis. For ELF3 overexpression experiment, ELF3 and controlled scrambled plasmids were purchased from GeneChem (Shanghai, China). Transfection was performed with Lipofectamine 2000 (Invitrogen) according to the manufacturer's instructions.
RNA Extraction and Quantitative RT-PCR
For RNA extraction, gastric cancer cells were grown in 6-well plates and transfected with siRNAs. After 48 h, remove the medium and directly add 800 μL Trizol into the 6-well plate to harvest samples. Total RNA was extracted using Trizol reagent (Invitrogen, USA) according to the manufacturer's instructions. The isolated RNA was treated with RNase-free DNase I (Roche) for 15–30 min as we described before (19). PCR was performed to ensure removal of genomic DNA by using RNA samples as templates. Reverse transcription was performed to obtain cDNA using the PrimeScript™ RT reagent Kit (Perfect Real Time, Takara) according to the manufacturer's instructions. For Quantitative RT-PCR, all the cDNA samples were diluted 5 times. The qPCR protocol was using One Step TB Green PrimeScript™ RT-PCR Kit II (Takara) according to the manufacturer's instructions. The qPCR analysis was conducted on Bio-Rad CFX Manager 3.1 real-time PCR system. The Cycling conditions were as follows: 95°C for 30 s, 95°C for 5 s and 60°C for 30 s. The reaction was performed for 40 cycles. The mRNA expression of ELF3, IRF6, and ZEB1 were determined by using the specific primers (ELF3-F: 5′-CACTGATGGCAAGCTCTTC-3′, ELF3-R: 5′-GGAGCG-CAGGAACTTGAAG-3′; IRF6-F: 5′-CCAGTAGTGGCTCGGATGAT-3′, IRF6-R: 5′-CAGCTCTCCTGGGTTTG AAG-3′; ZEB1-F: 5′-ACCTCTTCACAGGTTGCTCCT-3′, ZEB1-R: 5′-AGTGC AGGAGCTGAGAGTCA-3′, ACTIN-F: 5′-ATCGTCCACCGCA-AATGCTTCTA-3′, ACTIN-R: 5′-AGCCATGCCAATCTCATCTTGTT-3′). The experiment was performed with three replicates and average values were presented. Comparative quantification was determined using the 2−ΔΔCt Method.
Dual Luciferase Reporter Assay
For IRF6 promoter cloning, total genomic DNAs were extracted from AGS cell lines using TIANamp Genomic DNA Kit (TIANGEN, China). The 1,000 bp-length of IRF6 promoters were amplified by PCR (forward primer: GTGACCCATGCCTATATTT, reverse primer: GGGCGCCTGGCTCTACCCAA) and then clone into pMD18-T vectors (Takara) which were next used as templates for wildtype and mutant IRF6 promoters cloning. For reporter vectors construction (WT), the wildtype promoters of IRF6 (NM_006147) were amplified by PCR (forward primer: actagtacgcgtatttGAAAGAGAAAAAAGCAAACA, reverse primer: agagtttaaacgtcgacatttGGGCGCCTGGCTCTACCCAA). After digesting the pEZX-FR01 vector (GeneCopoeia, USA) into a linear fragment by SwaI, the IRF6 promoters and the linear vector fragments are recombined by seamless cloning (Yeasen Biotech, China) according to the manufacturer's instructions. For reporter vectors construction (Mut), two DNA fragments (mutant IRF6 promoter) was amplified by using the templates that mentioned above. The primers are as follows: F1, actagtacgcgtatttTTGAACTGGGTGCCA; R1, GATCGATCTTCTTTTTTGTTGTTG; F2, AAAAGAAGATCGATCTTATTATTCTCATTG, R2, agagtttaaacgtcgacatttTCTCCCCGTCCCGCAC. After digesting the pEZX-FR01 vector into a linear fragment by SwaI, the linear vector fragments and the mutant IRF6 promoters are recombined by seamless cloning according to the manufacturer's instructions. AGS and SGC7901 cells were seeded into 12-well-tissue plates 24 h before transfection, and then co-transfected with 5 ng siRNA and 1 mg plasmid using the Lipofectamine 2000 Reagent (Invitrogen), according to the manufacturer's instructions. After another 48 h, cells were assayed using the Dual-Luciferase reporter assay system kit (GeneCopoeia, USA). All experiments were performed in triplicate and data were pooled from three independent experiments.
Chromatin Immunoprecipitation Assay
Chromatin immunoprecipitation (CHIP) assays were performed using CHIP Assay Kit (56383S, Cell Signal Technology, USA) according to the manufacturer's protocol. ELF3 antibody were purchase from Santa Cruz Biotechnology (sc-376055, USA); ZEB1 antibody were purchase from Santa Cruz Biotechnology (sc-25388, USA). Briefly, SGC7901 cells were collected and fixed for 10 min at 37°C with 1% formaldehyde, followed in sequence with SDS lysis and DNA shearing, protein and DNA immunoprecipitation, cross-linked DNA reversal and DNA purification, and finally the immunoprecipitated DNA fragments were detected by PCR assays. The normal rabbit IgG was used as the negative control. The sequence of primers for CHIP (anti-ELF3) were used as follows (forward: TCTGAACTCCCAGTCGCTTC; reverse: TATGACACTCCG CGTTTCTG); the sequence of primers for CHIP (anti-ZEB1) were used as follows (forward: CGGGCGGATGCGAAGGCT; reverse: GGGCGCCTGGCTCTACCCAA).
Statistical Analysis
Data from at least three independent experiments performed in triplicate were presented as the mean ± S.D. The differences between two groups were determined by student's t-test. Differences among multiple groups were determined by one-way ANOVA. Comparisons were performed using the Spearman's correlation test; p < 0.05 was considered statistically significant.
Results
Reduced IRF6 Expression Was Clinically Correlated With Poor Prognosis in GC
The Human Protein Atlas (HPA) is an effort to map the subcellular location of all human proteins (http://www.proteinatlas.org/), which contains a large number of immunohistochemistry (IHC) images of sections from 46 different normal human tissue types and 20 different cancer types (20). To evaluate subcellular antibody staining patterns of IRF6 in normal and cancer tissue of stomach, the corresponding IHC images of IRF6 were downloaded from HPA web server. The results showed that IRF6 proteins were mainly located in cytoplasm and nucleus of glandular cells in normal stomach tissues and rarely expressed in cancerous tissues (Figure 1A). Furthermore, we analyzed the expression of IRF6 in two microarray gene profiling data (GSE54129 and GSE79973) that contained expression data of both normal tissues and cancer tissues. The results showed that the transcripts level of IRF6 in the gastric cancer tissues was lower than in the corresponding normal tissues (Figure 1B).
Figure 1
In order to understand the significance of decreased IRF6 expression in gastric cancer, we further determined the potential associations between IRF6 expression level and clinicopathological features in GSE62254 GC cohorts (n = 300). The results showed that diffuse GC tissues possessed lower IRF6 expression than intestinal GC tissues (Figure 1C). Moreover, IRF6 expression level was negatively correlated with advanced clinical TNM stages (Figures 1D–F). In addition, GC patients with higher expression of IRF6 tended to possess longer overall survival time and disease free time than patients with lower expression of IRF6 (Figures 1G,H). These results together suggested that lower IRF6 expression predicted poorer prognosis in GC.
IRF6 Was Directly Regulated by ZEB1 in GC
Previous study has reported that gastric cancer should be divided into four subtypes, including MSI, MSS/TP53+, MSS/TP53–, and MSS/EMT (18). During analysis of correlation between IRF6 expression level and clinicopathological features in GSE62254 GC cohorts, we noted that IRF6 expression in MSS/EMT subtype was remarkably lower than other three subtypes (Figure 2A). As ZEB1 was a master regulator of EMT process in many cancer types, we firstly considered to analyze the expression correlation of ZEB1 and IRF6 in GC. The results showed that IRF6 expression level was negatively associated with ZEB1 expression in both GC tissues and normal stomach tissues (Figures 2B–D).
Figure 2
To identify if ZEB1 could directly regulate IRF6 expression in GC, we firstly analyzed the available chip-seq data of ZEB1 using the online website of Cistrome. The result showed that an obvious ZEB1 peak was observed in the promoter region of IRF6 in different cell lines (Figure S1). Here, we took the chip-seq of ZEB1 in pancreatic cancer cell line as an example for further analysis. As shown in the bottom of Figure 2E, the most abundant section in the IRF6 promoter region pulled down by the ZEB1 antibody was located in 110 bp upstream of TSS of IRF6 (NM_006147). Besides, this section contained two ZEB1 binding sites predicted by JASPAR web server (Figures 2E–G), which suggested that IRF6 might be transcriptionally regulated by ZEB1. In order to confirm this possibility, we examined the IRF6 expression level after knockdown and overexpression of ZEB1 in AGS cell line (Figures 2H,I). The results showed that knockdown expression of ZEB1 significantly increased IRF6 expression; while overexpression of ZEB1 significantly decreased IRF6 expression (Figures 2J,K). Additionally, the result of CHIP assay showed that ZEB1 could directly bind on the IRF6 promoters in gastric cancer cells (Figure 2L). These results together indicated that transcription factor ZEB1 negatively regulated IRF6 expression via directly binding on IRF6 promoter.
IRF6 Was Positively Regulated by ELF3 in GC
IRF6 was reported to be positively regulated by TP63 in squamous cell carcinomas and normal skin tissues. To verify if the regulation of TP63 on IRF6 expression exists in other tissues, we downloaded the expression data of IRF6 and TP63 in varies kinds of normal tissues from GTEx Portal. Based on comparative analysis, we found that TP63 and IRF6 showed a very similar expression profile in most kinds of tissues, especially in human skin tissue. However, a certain degree of IRF6 expression can still be detected in the normal human tissues that rarely express TP63, such as stomach, colon and small intestine (Figure 3A). Therefore, there must be other genes that can regulate IRF6 expression in the human gastrointestinal tissues.
Figure 3
Transcription factors tends to possess similar expression patterns to their target genes. To identify which gene regulated IRF6 expression in human gastrointestinal tissues, we analyzed expression of all the transcription factors that co-expressed with IRF6. And the result indicated that transcription factor ELF3 showed a very similar expression profile with IRF6 in most kinds of tissues (Figure 3A). Furthermore, by comparing the correlation of between ELF3 and IRF6 expression in stomach, colon and small intestine tissues, we found that the IRF6 expression was highly relevant to ELF3 in the normal human gastrointestinal tissues (Figures 3B–D). Moreover, no obvious expression correlation was observed between TP63 and IRF6 in the normal human gastrointestinal tissues (Figures 3E–G).
Considering highly similar expression profile between ELF3 and IRF6 in normal human stomach tissues, we wanted to know if ELF3 and IRF6 were still co-expressed in the gastric cancer tissues. To explore the ELF3 and IRF6 expression profile in human GC tissues, five microarray gene profiling data were downloaded from GEO datasets. The GSE26942 dataset consists of 217 GC samples; GSE35809 consists of 70 GC samples; GSE54129 consists of 132 GC samples; GSE62254 consists of 300 GC samples; GSE63089 consists of 45 paired samples. It was identified that ELF3 was co-expressed with IRF6 in all five gastric cancer GSE datasets (Figure 4A). Furthermore, we investigated and analyzed the expression of ELF3 and IRF6 in GC tissues from TCGA by using the bioinformatics tool GEPIA, a web server for cancer and normal gene expression profiling and interactive analyses (21). The results showed that ELF3 and IRF6 were highly co-expressed in GC tissues form TCGA (Figure 4B). In addition, after examination of ELF3 and IRF6 expression in six GC cell lines and the normal gastric epithelial cells GES-1 by performing qRT-PCR analysis, we found that ELF3 and IRF6 were both overexpressed in N87, BGC823, SGC7901 and AGS, but both hardly expressed in HGC27 compared to their expression level in GES-1 (Figures 4C,D). These results together indicated that the expression of ELF3 and IRF6 in GC tissues and cell lines were highly relevant.
Figure 4
However, it is still unclear whether ELF3 could regulate IRF6 or not. Thus, loss-of-function and gain-of-function studies were conducted in GC cell lines. Here, we chose SGC7901 and AGS gastric cancer cell lines to perform the RNA interference experiment and ELF3 overexpression experiment. The results showed that knockdown of ELF3 expression in GC cell lines significantly decreased IRF6 expression (Figures 4E,F), while overexpression of ELF3 in GC cell lines remarkably increased IRF6 expression (Figures 4G,H). These results suggested that IRF6 was positively regulated by transcription factor ELF3 in GC.
ELF3 Directly Binds on the Promoter of IRF6
To evaluate if ELF3 directly regulate IRF6 expression or not, we firstly analyzed the only available chip-seq data of ELF3 antibody in the pancreatic ductal carcinoma, which obtained from the online website of Cistrome (22). As expected, an obvious ELF3 peak was observed in the 700 bp-length of IRF6 promoter (Figure S2A). Coincidentally, bioinformatics analysis shows that the 2,000 bp-length promoter contained a possible ELF3 binding site, which was located at 660–672 bp upstream of IRF6 transcription start site (Figure S2B).
Furthermore, our unexpected findings showed that ELF3 might regulate IRF6 expression via binding on the enhancers near the IRF6 gene. Recent genome-wide studies have established that enhancers can be defined as DNA sequences that bind the transcriptional co-activator p300/CREB1/SRF/MED1/BRD4, that bind H3K27AC, and then recruit RNA polymerase II to start transcription (23). Therefore, we analyze the available chip-seq data of these transcription factors according to the online website of Cistrome. The results showed that there are five enhancers located near the IRF6 gene (from 50 kb downstream of IRF6 gene to 50 kb upstream of IRF6 gene). An obvious ELF3 binding peak was found in all the five enhancer regions, which suggested that ELF3 might be a transcriptional regulator of IRF6 (Figure 5A). Moreover, it is worth mentioning that the enhancer 3 was located at 700 bp-length promoter of IRF6.
Figure 5
To detect the ELF3 binding activity on the IRF6 promoter, we designed a mutant IRF6 promoter that the predicted ELF3 binding sites by JASPAR (Figure 5B). After the two IRF6 promoter sequences were cloned into the luciferase reporter vector, the luciferase activity was examined in the AGS and SGC7901 cell lines that transfected with ELF3 siRNAs pool and negative control siRNA. As expected, a significant decrease was observed in the both AGS and SGC7901 cell that transfected with ELF3 siRNAs pool and wildtype IRF6 promoter luciferase reporter, compared to the cell transfected with negative control siRNA and wildtype IRF6 promoter luciferase reporter. However, no obvious changes were detected in the AGS and SGC7901 cell that transfected with mutant IRF6 promoter luciferase reporter (Figure 5C). In addition, chromatin immunoprecipitation (CHIP) result in SGC7901 further confirmed that ELF3 could indeed directly bind to the promoter of IRF6 (Figure 5D). These results together indicated that ELF3 was also a transcriptional regulator of IRF6.
Discussion
Interferon regulatory factors (IRFs) family, consisting of nine members in mammals, commonly possessed a novel helix-turn-helix DNA-binding domain and was proved to be involved in regulation of immunity and oncogenesis (24). Transcription factor IRF6, as a member of IRF family, was reported to play an essential role in keratinocyte differentiation, craniofacial development, salivary Glands and pancreas development (7, 25, 26). While the function of IRF6 in cancers has only been reported in squamous cell carcinomas and breast cancer. In squamous cell carcinomas, IRF6 expression is downregulated and knockdown of IRF6 promotes invasive behavior of tumor cells. In breast cancer, IRF6 overexpression results in a significant reduction of breast cancer cell numbers through arresting cell cycles (27). In other words, IRF6 plays a tumor suppressor role in squamous cell carcinomas and breast cancer. However, in addition to high expression in the breast and skin tissues, IRF6 is also highly expressed in gastrointestinal tissues according to analysis of the RNA-seq data in GTEx dataset (Figure 3A). Up to now, the expression profile and function of IRF6 in gastrointestinal cancer has not been reported yet.
In this study, we determined to investigate the function of transcription factor IRF6 in gastric cancer. Firstly, we identified that the expression levels of IRF6 protein and transcript are both significantly reduced in GC compared with normal stomach tissues (Figures 1A,B). Besides, the decreased expression level of IRF6 was closely related to malignant progression and poor prognosis of gastric cancer (Figures 1C–H). These results strongly suggested that IRF6 might function as a tumor suppressor in GC, which was similar to the role of IRF6 in squamous cell carcinomas and breast cancer. Therefore, we speculated that the mechanism of IRF6 downregulation and anti-cancer effect might be conservative at least in these cancers. Due to the reason why IRF6 was downregulated in these 3 cancer types remains largely unclear, our following studies focused on identifying the possible transcriptional regulators of IRF6 in gastric cancer.
Previous studies have reported that IRF6 was downregulated during EMT process of breast cancer and prostate cancer (11–13). Our results also showed that IRF6 was lowly expressed in MSS/EMT subtype of GC. Besides, a significant negative correlation between expression of ZEB1 and IRF6 was observed in GC. And the CHIP assay in gastric cancer cells showed that ZEB1 could bind on IRF6 promoter. More importantly, loss-of-function and gain-of-function studies showed that IRF6 was negatively regulated by ZEB1. Taken together, we firstly confirmed that transcription factor ZEB1 was a negative transcriptional regulator of IRF6 in GC. It is worth noting that lots of evidences show that transcription factor ZEB1 was overexpressed in gastrointestinal cancers (28–30). Therefore, the downregulation of IRF6 in GC or in EMT process of GC, SCC and breast cancer may be due to upregulation of ZEB1.
In addition, previous studies showed that IRF6 is a downstream target gene of the NOTCH signaling pathway and induced by the NOTCH signaling pathway in breast cancer and keratinocytes (9, 31, 32). NOTCH signaling pathway has been reported to play critical roles in the development and progression of human cancers through regulating ZEB1 expression and EMT pathway (33–38). For instance, knockdown of NOTCH3 upregulated ZEB1 expression in esophageal squamous cell carcinoma (39). Interestingly, analysis of microarray data of NOTCH3 knockdown (GSE27424) showed that silencing NOTCH3 results in a significant increase of ZEB1 and a significant decrease of IRF6 in esophageal squamous cell carcinoma. Based on our current findings, it is highly probable that the NOTCH signaling pathway regulates IRF6 expression through regulation of ZEB1.
On the other hand, we tried to identify possible transcriptional regulators that positively regulate IRF6 expression. Due to IRF6 was known to be a target gene of TP63 in skin development and squamous cell carcinomas, we firstly attempted to analyze the co-expression of IRF6 and TP63 in human gastrointestinal tissues. However, we found that a certain degree of IRF6 expression can still be detected in the human gastrointestinal tissues that hardly express TP63, such as stomach, colon and small intestine (Figures 3A,E–G). Therefore, there must be other genes that can regulate the expression of IRF6 in the human gastrointestinal tissues.
Inspired of this, we considered to analyze all the transcription factors that highly co-expressed with IRF6 and finally identified ELF3 as a possible positive transcriptional regulator of IRF6 in GC due to ELF3 was highly co-expressed with IRF6 in normal stomach tissues, GC tissues and GC cell lines (Figures 3C, 4A–D). In addition, knockdown of ELF3 significantly decreased IRF6 expression and overexpression of ELF3 significantly increased IRF6 expression in GC cell lines (Figures 4E–H). Besides, it's identified that there are at least five enhancers located near the IRF6 gene and transcription factor ELF3 could bind on all the enhancers (Figure 5A). Additionally, the dual luciferase reporter and CHIP assay showed that ELF3 could directly bind on the promoter of IRF6 (Figures 5B–D). These results together indicate that ELF3 is a transcriptional regulator of IRF6 in GC.
Nevertheless, previous studies has reported that promoter DNA methylation at CpG islands is a common mechanism used by cancer cells to repress expression of tumor suppressor genes (40). Besides, it's reported that IRF6 is aberrantly silenced by DNA methylation of the 5′ IRF6 CGI in melanoma and squamous cell carcinomas (10, 41). Therefore, we also analyzed the correlation between IRF6 expression and promoter DNA methylation in GC tissues from TCGA by bioinformatics method. As expected, a significant negative correlation between IRF6 expression and promoter DNA methylation was observed, suggested promoter DNA methylation was also a potent reason of IRF6 downregulation in GC (Figure S3).
In conclusion, transcription factor IRF6 was downregulated in GC. The reduced expression of IRF6 predicted poor prognosis of GC. In addition, ZEB1 is identified to be a negative transcriptional factor of IRF6 gene, while ELF3 is identified to be a positive transcriptional factor of IRF6 gene in GC. The downregulation of IRF6 in GC might be owing to overexpression of ZEB1 and the DNA methylation of IRF6 promoter.
Statements
Author contributions
The research was designed by YL and SQ. The experiments were performed by SQ, DL, PC, JW, XQ, XZ, and LX, and the data were analyzed by SQ. The manuscript was written by SQ.
Acknowledgments
This study was supported by grants from the National Natural Science Foundation of China (grant no. 81802375); the Natural Science Foundation of Hubei Provincial Department of Education (grant nos. Q20172105 and Q20182103); the Foundation of Health and Family planning Commission of Hubei Province (grant no. WJ2017F066); the Faculty Development Grants from Hubei University of Medicine (grant nos. 2016QDJZR07 and 2017QDJZR05). We also want to thank for the help of University of British Columbia China Scholarship Council 6-Month Visiting Scholar Overseas Study Program.
Conflict of interest
The 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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2019.00220/full#supplementary-material
Figure S1The chip-seq data of ZEB1 in different cancer cell lines was analyzed by Cistrome web server. An obvious peak located in the promoter region of IRF6 was observed in blood, kidney, pancreas, lung, and breast tissues.
Figure S2Bioinformatics analysis showed that ELF3 might be a transcriptional regulator of IRF6. (A) The chip-seq data of ELF3 in pancreatic cancer cell lines was analyzed by Cistrome web server. An obvious peak located in the 700 bp-length of promoter region of IRF6 was observed (B) The possible ELF3 binding sites in 2,000 bp-length of IRF6 promoter was analyzed by JASPAR.
Figure S3The IRF6 expression was negatively related to IRF6 promoter methylation. (A) An obvious CpG islands was located in the promoter region of IRF6. (B) The association between IRF6 expression and IRF6 promoter methylation in the TCGA stomach cancer cohort was analyzed by MEXPRESS web tool. (C) Genes that negatively correlated with IRF6 methylation was analyzed by using LinkedOmics web server. (D) The IRF6 expression and IRF6 promoter methylation in 372 stomach cancer samples of TCGA were analyzed by using LinkedOmics web server.
References
1.
TorreLABrayFSiegelRLFerlayJLortettieulentJJemalA. Global cancer statistics, 2012. CA Cancer J Clin. (2015) 65:87–108. 10.3322/caac.21262
2.
LiFXDuHLiSXLiuJT. The association between metabolic syndrome and gastric cancer in Chinese. Front Oncol. (2018) 8:326. 10.3389/fonc.2018.00326
3.
ShimizuDKandaMKoderaY. Review of recent molecular landscape knowledge of gastric cancer. Histol Histopathol. (2018) 33:11–26. 10.14670/HH-11-898
4.
ShahBKKhanalAHewettY. Second primary malignancies in adults with gastric cancer–a US population-based study. Front Oncol. (2016) 6:82. 10.3389/fonc.2016.00082
5.
XuLHuangTJHuHWangMYShiSMYangQet al. The developmental transcription factor IRF6 attenuates ABCG2 gene expression and distinctively reverses stemness phenotype in nasopharyngeal carcinomaCancer Lett. (2017) 431:230–43. 10.1016/j.canlet.2017
6.
IngrahamCKinoshitaAKondoSYangBSajanSTroutKet al. Abnormal skin, limb and craniofacial morphogenesis in mice deficient for interferon regulatory factor 6 (Irf6). Nat Genet. (2006) 38:1335–40. 10.1038/ng1903
7.
RichardsonRJJillDSaimonMHardmanMJLynnetteKBoot-HandfordRPet al. Irf6 is a key determinant of the keratinocyte proliferation-differentiation switch. Nat Genet. (2006) 38:1329–34. 10.1038/ng1894
8.
GritlilindeA. p63 and IRF6: brothers in arms against cleft palate. J. Clin. Invest. (2010) 120:1386. 10.1172/JCI42821
9.
†RestivoG†NguyenBC†DziunyczPRistorcelliERyanRJHZuysalYet al. IRF6 is a mediator of Notch pro-differentiation and tumour suppressive function in keratinocytes. Embo J. (2011) 30:4571. 10.1038/emboj.2011.325
10.
BottiESpalloneGMorettiFMarinariBPinettiVGalantiSet al. Developmental factor IRF6 exhibits tumor suppressor activity in squamous cell carcinomas. Proc Natl Acad Sci USA. (2011) 108:13710. 10.1073/pnas.1110931108
11.
BaileyCMZhilaKEShinjiKMargaryanNVSeftorREBWheatonWWet al. Mammary serine protease inhibitor (Maspin) binds directly to interferon regulatory factor 6: identification of a novel serpin partnership. J Biol Chem. (2005) 280:34210–7. 10.1074/jbc.M503523200
12.
EladKSylvieDMSimsAHPhilippeGHelenCMeehanRRet al. An in vitro model that recapitulates the epithelial to mesenchymal transition (EMT) in human breast cancer. PLoS ONE. (2011) 6:e17083. 10.1371/journal.pone.0017083
13.
RocaHHernandezJWeidnerSMceachinRCFullerDSudSet al. Transcription factors OVOL1 and OVOL2 induce the mesenchymal to epithelial transition in human cancer. PLoS ONE. (2013) 8:e76773. 10.1371/journal.pone.0076773
14.
Moreno-BuenoGPortilloFCanoA. Transcriptional regulation of cell polarity in EMT and cancer. Oncogene. (2008) 27:6958–69. 10.1038/onc.2008.346
15.
HuangYDuJMiYLiTGongYOuyangHet al. Long non-coding RNAs contribute to the inhibition of proliferation and EMT by Pterostilbene in human breast cancer. Front Oncol. (2018) 8:629. 10.3389/fonc.2018.00629
16.
YeungTLLeungCSWongKKGutierrez-HartmannAKwongJGershensonDMet al. ELF3 is a negative regulator of epithelial-mesenchymal transition in ovarian cancer cells. Oncotarget. (2017) 8:16951–63. 10.18632/oncotarget.15208
17.
LiuDSkomorovskaYSongJBowlerEHarrisRRavaszMet al. ELF3 is an antagonist of oncogenic-signalling-induced expression of EMT-TF ZEB1. Cancer Biol Ther. (2018) 20:90–100. 10.1080/15384047.2018.1507256
18.
CristescuRLeeJNebozhynMKimKMTingJCWongSSet al. Molecular analysis of gastric cancer identifies subtypes associated with distinct clinical outcomes. Nat Med. (2015) 21:449–56. 10.1038/nm.3850
19.
QinSSTangYHChenYPWuPZLiMRWuGJet al. Overexpression of the starch phosphorylase-like gene (PHO3) in Lotus japonicus has a profound effect on the growth of plants and reduction of transitory starch accumulation. Front Plant Sci. (2016) 7:1315. 10.3389/fpls.2016.01315
20.
Issac NiwasSAndreasKRSVirginieUPalanisamyPCarolineKMartinSet al. Automated classification of immunostaining patterns in breast tissue from the human protein atlas. J Pathol Inform. (2013) 4:S14. 10.4103/2153-3539.109881
21.
TangZLiCKangBGaoGLiCZhangZ. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. (2017) 45:W98–102. 10.1093/nar/gkx247
22.
LiuTOrtizJATaingLMeyerCALeeBZhangYet al. Cistrome: an integrative platform for transcriptional regulation studies. Genome Biol. (2011) 12:R83. 10.1186/gb-2011-12-8-r83
23.
KimTKHembergMGrayJMCostaAMBearDMWuJet al. Widespread transcription at neuronal activity-regulated enhancers. Nature. (2010) 465:182–7. 10.1038/nature09033
24.
SavitskyDTamuraTYanaiHTaniguchiT. Regulation of immunity and oncogenesis by the IRF transcription factor family. Cancer Immunol Immunother. (2010) 59:489–510. 10.1007/s00262-009-0804-6
25.
FakhouriWDMetwalliKNajiABakhietSQuispesalcedoANitschkeLet al. Intercellular genetic interaction between Irf6 and Twist1 during craniofacial development. Sci Rep. (2017) 7:7129. 10.1038/s41598-017-06310-z
26.
MetwalliKADoMANguyenKMallickSKinKFarokhniaNet al. Interferon regulatory factor 6 is necessary for salivary glands and pancreas development. J Dent Res. (2017) 97:226–36. 10.1177/0022034517729803
27.
BaileyCMAbbottDEMargaryanNVZhilaKEHendrixMJC. Interferon regulatory factor 6 promotes cell cycle arrest and is regulated by the proteasome in a cell cycle-dependent manner. Mol Cell Biol. (2008) 28:2235–43. 10.1128/MCB.01866-07
28.
JiaBLiuHKongQLiB. Overexpression of ZEB1 associated with metastasis and invasion in patients with gastric carcinoma. Mol Cell Biochem. (2012) 366:223–9. 10.1007/s11010-012-1299-6
29.
OkugawaYToiyamaYTanakaKMatsusitaKFujikawaHSaigusaSet al. Clinical significance of zinc finger E-box binding homeobox 1 (ZEB1) in human gastric cancer. J Surg Oncol. (2012) 106:280–5. 10.1002/jso.22142
30.
ZhangGJZhouTTianHPLiuZLXiaSS. High expression of ZEB1 correlates with liver metastasis and poor prognosis in colorectal cancer. Oncol Lett. (2013) 5:564–8. 10.3892/ol.2012.1026
31.
ZenginT. Role of IRF6 as a Mediator of Notch Signaling in Breast Cancer Cell Lines. Izmir: Izmir Institute of Technology (2014).
32.
ZenginTEkinciBKucukkoseCYalcin-OzuysalO. IRF6 is involved in the regulation of cell proliferation and transformation in MCF10A cells downstream of Notch signaling. PLoS ONE. (2015) 10:e0132757. 10.1371/journal.pone.0132757
33.
WangZLiYKongDSarkarFH. The role of Notch signaling pathway in epithelial-mesenchymal transition (EMT) during development and tumor aggressiveness. Curr Drug Targets. (2010) 11:745–51. 10.2174/138945010791170860
34.
NatsuizakaMOhashiSNaganumaSKalmanRAOhyamaARhoadesBet al. Notch regulates squamous differentiation, cell plasticity and tumor heterogeneity in esophageal carcinoma. Cancer Res. (2011) 71:5194. 10.1158/1538-7445.AM2011-5194
35.
KangHAnHJSongJYKimTHHeoJHAhnDHet al. Notch3 and Jagged2 contribute to gastric cancer development and to glandular differentiation associated with MUC2 and MUC5AC expression. Histopathology. (2012) 61:576–86. 10.1111/j.1365-2559.2012.04274.x
36.
EspinozaIMieleL. Notch inhibitors for cancer treatment. Pharmacol Ther. (2013) 139:95–110. 10.1016/j.pharmthera.2013.02.003
37.
LiYMaJQianXWuQXiaJMieleLet al. Regulation of EMT by Notch signaling pathway in tumor progression. Curr Cancer Drug Targets. (2013) 13:957–62. 10.2174/15680096113136660101
38.
NowellCSRadtkeF. Notch as a tumour suppressor. Nat Rev Cancer. (2017) 17:145–59. 10.1038/nrc.2016.145
39.
OhashiSNatsuizakaMNaganumaSKagawaSKimuraSItohHet al. A Notch3-mediated squamous cell differentiation program limits expansion of EMT-competent cells that express the ZEB transcription factors. Cancer Res. (2011) 71:6836–47. 10.1158/0008-5472.CAN-11-0846
40.
JonesPALairdPW. Cancer epigenetics comes of age. Nat Genet. (1999) 21:163–7. 10.1038/5947
41.
NobeyamaYNakagawaH. Silencing of interferon regulatory factor gene 6 in melanoma. PLoS ONE. (2017) 12:e0184444. 10.1371/journal.pone.0184444
Summary
Keywords
IRF6, ELF3, ZEB1, transcriptional regulation, gastric cancer
Citation
Li D, Cheng P, Wang J, Qiu X, Zhang X, Xu L, Liu Y and Qin S (2019) IRF6 Is Directly Regulated by ZEB1 and ELF3, and Predicts a Favorable Prognosis in Gastric Cancer. Front. Oncol. 9:220. doi: 10.3389/fonc.2019.00220
Received
09 July 2018
Accepted
12 March 2019
Published
04 April 2019
Volume
9 - 2019
Edited by
Rupert Langer, University of Bern, Switzerland
Reviewed by
Toru Furukawa, School of Medicine, Tohoku University, Japan; Jens Neumann, Ludwig Maximilian University of Munich, Germany
Updates
Copyright
© 2019 Li, Cheng, Wang, Qiu, Zhang, Xu, Liu and Qin.
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: Ying Liu ying_liu1002@163.comShanshan Qin qinss77@163.com
This article was submitted to Gastrointestinal Cancers, a section of the journal Frontiers in Oncology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.