- Open Access
Aberrant activation of CYR61 enhancers in colorectal cancer development
Journal of Experimental & Clinical Cancer Research volume 38, Article number: 213 (2019)
High expression of secreted matricellular protein cysteine-rich 61 (CYR61) correlates with poor prognosis in colorectal cancer (CRC). Aberrant enhancer activation has been shown to correlate with expression of key genes involved in cancer progression. However, such mechanisms in CYR61 transcription regulation remain unexplored.
Expression of CYR61 was determined by immunohistochemistry (IHC), quantitative real-time PCR (qRT-PCR) and western blotting (WB) in CRC patients paraffin specimens and colon cell lines. ChIP-seq data of enhancer-characteristic histone modifications, in CRC tissues from the Gene Expression Omnibus (GEO) database, were reanalyzed to search for putative enhancers of CYR61. Dual-luciferase reporter assay was used to detected enhancer activity. Physical interactions between putative enhancers and CYR61 promoter were detected by chromosome conformation capture (3C) assay. Histone modification and transcription factors (TFs) enrichment were detected by ChIP-qPCR. Additionally, biological function of enhancers was investigated by transwell migration assays.
CRC tissues and cell lines expressed higher level of CYR61 than normal colon mucosa. Three putative enhancers located downstream of CYR61 were found in CRC tissues by ChIP-seq data reanalysis. Consistent with the ChIP-seq analysis results in the GEO database, the normal colon mucosal epithelial cell line NCM460 possessed no active CYR61 enhancers, whereas colon cancer cells exhibited different patterns of active CYR61 enhancers. HCT116 cells had an active Enhancer3, whereas RKO cells had both Enhancer1 and Enhancer3 active. Pioneer factor FOXA1 promoted CYR61 expression by recruiting CBP histone acetyltransferase binding and increasing promoter-enhancer looping frequencies and enhancer activity. CBP knockdown attenuated H3K27ac enrichment, promoter-enhancer looping frequencies, and enhancer activity. Small molecule compound 12-O-tetradecanoyl phorbol-13-acetate (TPA) treatment, which stimulated CYR61 expression, and verteporfin (VP) treatment, which inhibited CYR61 expression, confirmed that the enhancers regulated CYR61 expression. Knockdown and ectopic expression of CYR61 rescued cell migration changes induced by over-expressing and knockdown of FOXA1, respectively.
CYR61 enhancer activation, mediated by FOXA1 and CBP, occurs during CRC progression to up-regulate CYR61 expression and promote cell migration in CRC, suggesting inhibition of recruitment of FOXA1 and/or CBP to CYR61 enhancers may have therapeutic implications.
CYR61 (cysteine-rich 61/CCN1) belongs to the CYR61/CTGF/NOV (CCN) protein family . As a secreted matricellular protein, CYR61 can bind directly to various integrin receptors and heparan sulfate proteoglycans to regulate many cellular functions in a cell type- and context-dependent manner [1,2,3]. High expression of CYR61 is observed in colon cancer tissues and is closely related to shorter survival in colon cancer patients, and has been reported to promote cancer metastasis and cell migration [4,5,6,7,8,9]. Transcription factors (TFs), such as SOX4  and FOXK1 , can up-regulate CYR61 expression by binding to the CYR61 promoter in colon cancer cells. However, the research of CYR61 transcriptional regulation in CRC is limited.
Enhancers are important cis-regulatory elements and function as integrated TF docking platforms . Dysregulation of enhancers induces aberrant gene expression that drives the uncontrolled proliferation of human cancers, including colon cancer [11, 12]. Generally, active enhancers are loaded with lineage-specific TFs (sequence-dependent) and coactivator proteins (which lack sequence-specific DNA-binding) . Mediator complex, composed of coactivator proteins, can simultaneously bind to different TFs to mediate enhancer-promoter looping [14, 15] and facilitate delivery of important accessory factors to the promoter to potentiate transcription . Enhancer RNAs (eRNAs), which are transcribed from the enhancers, are considered to be reliable markers for active enhancer activity [17, 18]. Epigenetic modifications, such as histone modifications, can influence enhancer activity, with some histone modifications being considered hallmarks of enhancers and enhancer activity . In particular, DNA elements decorated with monomethylated H3 lysine 4 (H3K4me1) alone are considered to be primed enhancers and, when combinatorially deposited with acetylated H3 lysine 27 (H3K27ac), are considered to be active enhancers [16, 19].
CBP is closely related to its paralogue p300. As a co-activator, CBP can bind to TFs and bridge them to large protein complexes, such as Mediator complex . Moreover, CBP can act as a lysine acetyltransferase to acetylate TFs and histones to increase the accessibility of chromatin [20, 21]. Forkhead box A1 (FOXA1) has been reported to be a pioneer factor in that it can bind to and open compacted chromatin to facilitate the binding of other TFs, and its activity is dysregulated in many physiological and pathological conditions [22, 23]. Other forkhead transcription factors, such as FOXK1  and FOXO3a , have been reported to regulate CYR61 expression by binding to the CYR61 promoter in a sequence-dependent manner. Thus, FOXA1 might also regulate CYR61 expression in a sequence-dependent manner.
Although enhancers play important roles in gene transcriptional regulation, the role of enhancers in regulating CYR61 transcription in human colon cancer remains unexplored. By analyzing histone modification hallmarks of enhancers in colon cancer, we identified three putative enhancers located downstream of CYR61. In this study, we demonstrate that two of these three enhancers were aberrantly activated in CRC and, combined with FOXA1 and CBP, play key roles in the activation of CYR61 expression in colorectal cancer.
Patients and tissues
A total of 42 cases of colonic adenocarcinoma (along with corresponding matched normal colonic mucosa), obtained from colon cancer patients who underwent surgical treatment from 2015 to 2016 at Shantou Central Hospital (Guangdong, China), were available for examination in this study. Tissues were embedded in paraffin. All collected tissues had pathologic diagnoses by two independent pathologists at the Shantou Central Hospital. Patients who received radiotherapy or chemotherapy prior to surgery were excluded.
Specimen slides were incubated with anti-CYR61 primary rabbit polyclonal antibody (1:200; ThermoFisher, PA116579) overnight at 4 °C. Anti-mouse/rabbit HRP-labeled secondary antibody (Maixin Biological Technology Co. Ltd., KIT-5010) was applied for 20 min at RT. Slides were then incubated with diaminobenzidine (DAB) chromogen substrate and counterstained with hematoxylin. Image-Pro Plus v.6.0 (IPP 6.0) software (Media Cybernetics, Inc., USA) was used to assess the area and the integrated optical density (IOD) value of the stained region. Mean density = IOD/area. The average mean density for five random fields at 100× magnification was used for CYR61 statistical analysis.
The mRNA expression levels of CYR61 in normal colonic samples and colon cancer samples were analyzed by ONCOMINE microarray datasets (https://www.oncomine.org), with a cut-off p-value of 0.001, and fold change of 2.0 [25, 26].
The normal human colon mucosal epithelial cell line NCM460 (INCELL Corporation, USA) was cultured in M3D medium (INCELL Corporation, USA). The human colon carcinoma cell line LoVo (Cell Bank of the Chinese Academy of Sciences, China) was cultured in basic DMEM/F-12, HEPES (Gibco, USA); C2BBe1 [Caco-2 cell clone] (ATCC, USA) was cultured in DMEM (Gibco, USA) with 0.01 mg/ml human transferrin; the HCT116 cell line (Cell Bank of the Chinese Academy of Sciences, China) was cultured in McCoy’s 5A (modified) medium (McCoy’s 5A) (Gibco, USA); SW480 (ATCC, USA) was cultured in Leibovitz’s L-15 Medium (Gibco, USA); RKO (ATCC, USA) was cultured in minimum essential medium (MEM) (Gibco, USA). All cell culture media was supplemented with 10% fetal bovine serum. All cell lines were cultured in a humidified atmosphere at 37 °C. SW480 was cultured with atmospheric air, and other cell lines were incubated with 5% CO2.
Compounds used in this research were as follows: 12-O-tetradecanoyl phorbol-13-acetate (TPA) (Sigma, P8139) and verteporfin (VP) (Sigma, B1583–5).
RNA extraction, RT-PCR and quantitative PCR
Total RNA was isolated using Trizol (Takara, No.9109) according to the manufacturer’s instructions. Reverse transcription was performed using a PrimeScript™ RT reagent Kit with gDNA Eraser (Takara, RR047A). Quantitative real-time PCR was performed using AceQ qPCR SYBR Green Master Mix (Low ROX Premixed) (Vazyme, Q131–02) on a QuantStudio 12 K Flex Real-Time PCR System (ThermoFisher, USA) according to the manufacturer’s protocols. Primers we designed are shown in Additional file 1: Table S1.
Western blot assay (WB)
Antibodies used for western blot were as follows: anti-CYR61 (CST, #14479), anti-CBP (ThermoFisher, PA1–847), anti-FOXA1 (ThermoFisher, PA5–27157), and anti-β-actin (Santa Cruz, sc-130,656).
ChIP-seq and GRO-seq data analyses
ChIP-seq data analyses were used to search for putative enhancers of CYR61, and GRO-seq data analysis was used to detect expression of eRNAs at putative enhancer regions. ChIP-seq data of colon adenocarcinoma samples and normal colon tissues were downloaded from the GEO database, Data for ChIP-seq and GRO-seq were downloaded from the GEO database and listed in Additional file 1: Table S2. Raw data were aligned using Bowtie2  (version 2.2.9) to the human genome (build hg38, GRCh38) with parameter -p 15 -x after quality filtering. MACS  (version 1.4.2) was used with the parameter “-f BAM -g hs -n sample -p 1e-5 -B” for peak calling. The ChIP-seq and GRO-seq figures were visualized in the UCSC Genome Brower.
Chromatin immunoprecipitation (ChIP)
ChIP-qPCR was used to detect the histone modification and TF binding at the promoter and enhancer regions of CYR61. As we described before , the ChIP assay was performed using a Magna ChIP™ G Chromatin Immunoprecipitation Kit (Millipore, 17–611) according to the manufacturer’s protocol. ChIP-grade antibodies were as follows: 2 μg anti-RNA polymerase II (Millipore, 05–623), 2 μg anti-H3K27ac (ThermoFisher, 720,096), anti-H3K4me1 (Abcam, ab8895), 2 μg anti-H3K4me3 (Millipore, 17–614), 3 μg anti-CBP (ThermoFisher, PA1–847), and 3 μg anti-FOXA1 (ThermoFisher, PA5–27157). Immunoprecipitated DNA was detected by qPCR and normalized with input DNA. The sequences of the primers are listed in Additional file 1: Table S3.
RNA interference (RNAi) assay
All siRNAs were designed by GenePharma (Shanghai, China). Cells were transfected with siRNA using Lipofectamine RNAiMAX Reagent (Invitrogen, 13,778,100) according to the manufacturer’s protocol. The sequences of the siRNAs are listed in Additional file 1: Table S4.
Dual-luciferase reporter assay
Dual-luciferase reporter assays were used to detect putative enhancer activity. The pGL3-basic Vector (Promega, E1751) and pRL-SV40 Vector (Promega, E2231) were purchased. The pGL3-basic-promoter plasmid was constructed by inserting the CYR61 promoter region into the pGL3-basic plasmid, and the putative enhancer fragment or negative control region was inserted into the pGL3-basic-promoter plasmid to construct the pGL3-basic-promoter-enhancer/NC plasmids. Primers and insertion sites for dual-luciferase reporter assays are shown in Additional file 1: Table S5. All constructed plasmids were verified by sequencing. Plasmid transfection was performed with FuGENE HD transfection reagent (Promega, E2312) according to the manufacturer’s protocol. Forty-eight to seventy-two hours after transfection, cells were harvested. Luciferase activity was determined by the dual-luciferase reporter assay system (Promega, E1910).
Chromosome conformation capture (3C) assay
Chromosome conformation capture (3C) assay was used to detect physical interactions between enhancer and promoter [30, 31]. The 3C assay was performed as described previously [32, 33]. HindIII was used for genomic DNA digestion. Plasmids used in the 3C assays were BAC clones for ERCC3 (Invitrogen, CTD-3251 N23) and CYR61 (Invitrogen, RP11-963G4). Primers used in this assay are shown in Additional file 1: Table S6.
Cell migration assays were performed in a transwell chamber (24-well, 8 μm pore size; Corning). Plasmids GV219-mCYR61 and GV219-mFOXA1 (GeneChem, China) were used to over-express CYR61 and FOXA1. Cells were transfected with siRNA at 25 pmol and/or GV219 plasmid, 2.8 μg each, in a 6-well plate for 48~72 h, and then 1.0 × 105 cells were transplanted into a transwell chamber. After 48 h, migrated cells were stained with crystal violet. Five random fields at 200× magnification were used for cell counting for each membrane.
FOXA1 overexpression and CYR61 shRNA lentivirus vectors transduction
Lentiviral vector GV344 (LV-firefly_Luciferase-puromycin/shCYR61) was used to knockdown CYR61, and lentiviral vector CV572 (LV-Cherry-neomycin/FOXA1) was used to overexpress FOXA1. All lentiviral vectors were constructed by the Shanghai GenePharma Corporation (Shanghai, China). The targeting sequence of shCYR61 was 5′- GCATCCTATACAACCCTTT − 3′. A six-well plate was inoculated with 2 × 105 HCT116 cells. Twenty-two hours after inoculation, LV-firefly_Luciferase-puromycin/shNC and LV-firefly_Luciferase-puromycin/shCYR61 lentivirus were used to infect cultured cells. Fourteen hours later, the medium in each well was changed with fresh normal medium. Seventy-two hours after infection, 2 μg/ml puromycin (Sangon Biotech, China, # A610593) was used to screen stable cell clones. The transduction of FOXA1 overexpression lentivirus was performed after puromycin screening for 7 days and the transduction procedure was as shCYR61‘s. Seventy-two hours after infection with LV-Cherry-neomycin/con and LV-Cherry-neomycin/FOXA1, 800 μg/ml G418 (Sangon Biotech, China, # B540723) was used to select stable cell clones.
Six-week old male NOD SCID mice (Beijing Vital River Laboratory Animal Technology Co., Ltd., China, #406) weighing between 19 g to 22 g were used in our animal study. HCT116 cells (3 × 106 cells) transduced with lentivirus were tail vein injected into the mice (n = 5 per group). Tumor metastasis were monitored weekly by an in vivo imaging system (Perkin Elmer IVIS Lumina III) from the third week after injection. Mice were intraperitoneally injected with 150 μg D-luciferin (Solarbio, China, #D8390) solution per gram body weight for 10 min before subjected to bioluminescent imaging. Tumor size and metastasis were quantified using Living Image software (Perkin Elmer IVIS Lumina III).
Statistical analyses were conducted in SPSS Statistics 19.0. All data shown were determined for three independent experiments unless otherwise stated, and presented as the mean ± S.D., *P < 0.05, **P < 0.01, ***P < 0.001.
CYR61 is elevated in primary colonic adenocarcinoma tissue and colon cancer cell lines
We initially performed immunohistochemical analysis (IHC) to examine CYR61 protein expression in specimens from 42 cases of colonic adenocarcinomas and matched para-carcinoma tissues. The results showed that the expression levels of CYR61 protein in carcinoma tissue were significantly higher than those in para-carcinoma tissue (Fig. 1a, b). A significant association was noted between CYR61 protein expression and histological tumor grade (p < 0.001) (Fig. 1c). No significant correlation was seen between CYR61 expression and other clinicopathological parameters examined (Additional file 1: Figure S1, Table S7). ONCOMINE microarray dataset analysis from the Hong Colorectal dataset and Skrzypczak Colorectal 2 dataset showed that CYR61 mRNA was significantly up-regulated in primary colon adenocarcinoma tissues compared with normal controls (Fig. 1d, e), consistent with our tissue specimen analysis, indicating that CYR61 expression is up-regulated in primary colonic adenocarcinoma tissue. Furthermore, higher mRNA and protein expression of CYR61 was observed in several colon cancer cell lines (LoVo, C2BBe1, HCT116, SW480, and RKO), compared with the normal human colon mucosal epithelial cell line NCM460 (Fig. 1f, g). Thus, these results indicate that high expression of CYR61 is associated with tumor development and aggressiveness in colonic adenocarcinoma tissue and cell lines.
CYR61 expression is regulated by enhancers in colon cancer
Several recent studies have shown that enhancers promote gene transcription via a long-range interaction with their cognate promoters and distinct histone modifications in surrounding nucleosomes [16, 34], so we determined whether there were putative enhancers that interacted with the CYR61 promoter and promoted CYR61 expression in colon cancer cell lines. Firstly, we analyzed ChIP-seq data, of colon adenocarcinoma samples and normal colon tissues, downloaded from the GEO database and found three putative enhancers located at 26.5 kb, 31.5 kb, 48.6 kb downstream of the CYR61 transcriptional start site (Fig. 2a), referred to as Enhancer1 (chr1:85605899–85,608,570), Enhancer2 (chr1: 85610948–85,613,635) and Enhancer3 (chr1: 85628077–85,630,633), respectively. GRO-seq data of HCT116 showed all three regions transcribed eRNAs, denoted E1-eRNA, E2-eRNA, and E3-eRNA. These features are consistent with characteristics of an activated enhancer .
We detected eRNA expression at the three putative enhancer regions in NCM460 and colon cancer cell lines. In agreement with the results of the GRO-seq analysis, real-time quantitative PCR analysis showed that colon cancer cells expressed dramatically higher eRNA levels than NCM460 (Fig. 2b).
To investigate CYR61 enhancer activity in the above cell lines, we detected histone modification and RNA polymerase II (Pol II) enrichment in NCM460, HCT116 and RKO cell lines by ChIP-qPCR. As shown in Fig. 2c, H3K27ac, H3K4me1 and RNA polymerase II were not enriched at the three putative enhancer regions in NCM460 cells, but were enriched at Enhancer3 in HCT116 cells. H3K27ac ChIP-seq data, downloaded from GEO, showed high H3K27ac enrichment at the Enhancer3 region in HCT116 cells, but low H3K27ac enrichment at Enhancer1 and Enhancer2 regions (Additional file 1: Figure S2A), consistent with our ChIP-qPCR assay. In RKO cells, the H3K27ac and Pol II enrichment was observed at both Enhancer1 and Enhancer3 regions. Next, a dual-luciferase reporter assay was used to further determine the enhancer activity. The fragments inserted into the plasmid are shown in Fig. 2a. The inserted CYR61 promoter showed significantly higher promoter activity in HCT116 and RKO cells than that in NCM460 cells (Additional file 1: Figure S2B). Furthermore, all enhancers, inserted in both forward and reverse directions, showed enhancer activity in the above colon cancer cell lines (Fig. 2d).
As mentioned above, typically, enhancers contact with the promoters of their cognate gene through long-range interactions. We next determined whether the putative enhancers interacted with the CYR61 promoter by using a chromatin conformation capture (3C) assay. As shown in Fig. 2e, a strong interaction between Enhancer3 and the CYR61 promoter was identified in colon cancer cells HCT116 and RKO. The ninth test fragment (9F in Fig. 2a), which overlaps with the Enhancer3 region, displayed significantly higher interaction frequency with the CYR61 promoter compared to the neighboring DNA fragments (1F in Fig. 2a). Therefore, the 2.5 kb region located 48.6 kb away from the CYR6 promoter has all characteristics of an active enhancer (high level of H3K27ac and H3K4me1, transcription enhancement, and interaction with the CYR61 promoter). Interestingly, the interaction between Enhancer1 and the CYR61 promoter occurred at high frequency in RKO cells, but not in HCT116 cells. In NCM460 cells, no interaction between these putative enhancers and the CYR61 promoter was observed. These results imply that colon cancer cells exhibit different CYR61 regulation patterns by different enhancers. Moreover, the differential activity of enhancers between these cells further underscores the concept that these enhancers are only activated in cancer cells and promoted CYR61 transcription.
FOXA1 regulates enhancer activity to promote CYR61 expression
Most enhancer activators function by binding sequence-specifically to enhancers and promoters and forming protein-protein interactions with RNA polymerase and other general TFs. So, we determined the sequence-dependent TFs that mediate CYR61 promoter and enhancer interaction. By analyzing the CYR61 promoter sequence (chr1:85579761–85,580,961), Enhancer1 sequence (chr1:85605899–85,608,570) and Enhancer3 sequence (chr1:85628077–85,630,633) using transcription factor affinity prediction (TRAP) web tools (http://trap.molgen.mpg.de/cgi-bin/trap_multi_seq_form.cgi) , we found 6 TFs with statistically significant p-values (Fig. 3a). FOXA1, a pioneer factor in many physiological and pathological conditions was among the TFs listed and showed enrichment at TRAP-predicted loci (Fig. 3b). Furthermore, pioneer factors can bind to compacted chromatin to facilitate the binding of other TFs [22, 23]. Based on this, we reasoned that FOXA1 may regulate CYR61 enhancer activity and promote CYR61 expression.
We determined the enrichment of FOXA1 at the TRAP-predicted loci. FOXA1 was not enriched at TRAP-predicted loci in NCM460 cells, but enriched at the CYR61 promoter, Enhancer1 and Enhancer3 loci in RKO cells, and enriched at the CYR61 promoter and Enhancer3 loci in HCT116 cells (Additional file 1: Figure S3A-C). We next determined whether knockdown of FOXA1 could decrease enhancer activity and lower CYR61 expression. Supporting a role for FOXA1 in enhancing CYR61 expression, knockdown of FOXA1 expression with FOXA1-specific siRNA markedly reduced CYR61 expression both at the mRNA and protein levels in RKO and SW480 cells, and the eRNA levels were also reduced in RKO cells as well (Fig. 3c, d). Moreover, ChIP-qPCR results showed that FOXA1 knockdown reduced FOXA1 enrichment at both Enhancer1 and Enhancer3 loci, and reduced H3K27ac enrichment at the target loci (Fig. 3e). In parallel with the decrease of FOXA1 occupancy at the CYR61 promoter, Enhancer1 and Enhancer3 loci, and the CYR61 promoter/Enhancer1 and CYR61 promoter/Enhancer3 looping frequencies were decreased (Fig. 3f). To distinguish the contribution of FOXA1 to the individual CYR61 promoter, Enhancer1, and Enhancer3 activities, we knocked down FOXA1 by siRNA and performed a dual-luciferase assay. Both promoter-luc-Enhancer1+ and promoter-luc-Enhancer3+ showed decreased activity after FOXA1 knockdown, but the activity of promoter-luc-enhNC plasmid (which contained an NC fragment instead of an enhancer fragment) tended to be slightly decreased, but was not statistically different from the negative control group (Fig. 3g).
To further confirm the above results, we over-expressed FOXA1 by transfecting a GV219-mFOXA1 plasmid in HCT116 and SW480 cells, then detected CYR61 expression. Increasing FOXA1 expression markedly increased CYR61 expression at both the mRNA and protein levels, and this paralleled an increase in eRNA levels (Fig. 3h, i). In addition, FOXA1 over-expression increased both FOXA1 enrichment and H3K27ac enrichment at the promoter and Enhancer3 loci in HCT116 cells (Fig. 3j), and the CYR61 promoter/Enhancer3 looping frequencies were increased by FOXA1 over-expression (Fig. 3k). The dual-luciferase assay showed both promoter-luc-enhNC and promoter-luc-Enh3+ plasmid activities were increased after FOXA1 over-expression, with promoter-luc-Enh3+ showing a greater increase than promoter-luc-enhNC (Fig. 3l). Together, these results demonstrate that FOXA1 regulates CYR61 expression mainly by regulating its enhancer activity.
CBP collaborates with FOXA1 to regulate CYR61 enhancer activity and CYR61 expression
In the above experiments, we also found that with up- or down-regulation of FOXA1, the enrichment of CREB-binding protein (CBP) at Enhancer1, Enhancer3 and promoter loci changed accordingly (Figs. 3e, j, 4a, b, Additional file 1: Figure S3D-F). To validate an involvement of the histone acetyltransferase CBP in the regulation of enhancer activity, we knocked down CBP by siRNA and detected both CYR61 expression and enhancer activity. The results showed that, following CBP knockdown, both mRNA and protein levels of CYR61 were decreased in RKO and SW480 cells (Fig. 4c, d), and the expressions of E1-eRNA and E3-eRNA were likewise decreased in RKO cells (Fig. 4c). ChIP-qPCR analysis showed that, after CBP knockdown, CBP enrichment was statistically reduced at the Enhancer1 locus, and caused a decrease of H3K27ac enrichment at the CYR61 promoter, Enhancer1 and Enhancer3 loci (Fig. 4e). Parallel to the decreased CBP enrichment, CYR61 promoter/Enhancer1 looping frequencies and promoter-luc-Enhancer1+ plasmid activity were decreased after CBP knockdown (Fig. 4f, g). These data indicate that CBP is required for the maintenance of H3K27Ac histone modification and promoter-enhancer 3D structure. Additionally, CBP and FOXA1 knockdown did not affect each other’s expression level (Fig. 4h-k). Interestingly, FOXA1 knockdown decreased the CBP enrichment at the target loci (Fig. 4a), while CBP knockdown did not affect FOXA1 enrichment (Fig. 4l), suggesting that enrichment of FOXA1 is necessary for the recruitment of CBP at target loci.
Small molecule compounds regulate CYR61 expression through affecting enhancer activity
12-O-tetradecanoyl phorbol-13-acetate (TPA) is reported to up-regulate expression of CYR61 in breast cancer cells . Here, TPA treatment promoted the expression of CYR61 in both HCT116 and RKO cells (Fig. 5a, b, Additional file 1: Figure S4A-C), but not in NCM460 cells. To investigate whether the enhancers were affected, we detected eRNA expression and found that E3-eRNA in HCT116 cells, and both E1-eRNA and E3-eRNA in RKO cells were elevated by 100 ng/ml TPA treatment for 2 h (Fig. 5b). Further, our ChIP-qPCR results showed that TPA treatment caused H3K27ac, H3K4me3 and CBP enrichment at the CYR61 promoter and Enhancer3 region in HCT116 cells (Fig. 5c), indicating that TPA treatment induced enhancer activity. In addition, the above findings were confirmed by 3C analysis showing that CYR61 promoter/Enhancer3 looping frequencies were statistically increased in HCT116 after TPA treatment (Fig. 5d).
In addition to TPA, we characterized the effects of another small molecule drug, verteporfin (VP), a benzoporphyrin derivative, which has been reported to inhibit cell growth in cancer cells through down-regulating CYR61 expression [37, 38]. Therefore, we tested whether VP treatment could inhibit expression of CYR61 and suppress enhancer activity in colon cancer cells. VP treatment inhibited the expression of CYR61 mRNA and protein in both RKO and HCT116 cells, E3-eRNA in HCT116 cells, and both E1-eRNA and E3-eRNA in RKO cells (Fig. 5e, f, Additional file 1: Figure S4D-F). Since CYR61 mRNA levels were decreased in a time-dependent manner after treatment with 10 μg/ml VP for various times and bottomed out at 2 h (Additional file 1: Figure S4 E), we inferred that H3K4me1, H3K27ac, and CBP enrichment would bottomed out at 2 h also, but in fact, only H3K4me1 and H3K27ac enrichment at CYR61 promoter region showed significant decreases after VP treatment for 2 h. Considering that TF enrichment and histone modifications may change earlier than the target gene does, we tested the TF enrichment and histone modifications after VP treatment for 1 h, and found that treatment with VP for 1 h significantly inhibited H3K4me1, H3K27ac and CBP enrichment at the CYR61 promoter, Enhancer1 and Enhancer3 regions in RKO cells (Fig. 5g). These results suggested that the CBP enrichment and active histone modifications changed earlier than the target gene. CYR61 promoter/Enhancer1 and CYR61 promoter/Enhancer3 looping frequencies were statistically reduced after treatment of RKO cells with VP, as shown by 3C assay (Fig. 5h). These data suggest that TPA and VP regulate CYR61 expression at least partially through regulating CYR61 enhancer activity.
FOXA1 promotes colon cancer cell metastasis by increasing CYR61 enhancer activity
Since CYR61 is a key regulator of migration [7, 8, 39, 40], a transwell assay was performed to explore the functional roles of the enhancers in cancer progression. Ectopic expression of CYR61 increased the migration of HCT116 cells (Figs. 6a, 7b). Conversely, knockdown of CYR61 decreased the migration of RKO cells (Fig. 6c, d). Ectopic FOXA1 expression up-regulated CYR61 expression and increased HCT116 cell migration. Conversely, CYR61 knockdown, in cells ectopically expressing FOXA1, decreased cell migration (compared with ectopic FOXA1 expressing cells) (Fig. 6e, f). Accordingly, FOXA1 knockdown down-regulated CYR61 expression and decreased RKO cell migration, whereas ectopic expression of CYR61 in FOXA1-knockdown cells increased cell migration (Fig. 6g, h).
To verify whether FOXA1 could promote colon cancer cells metastasis by up-regulating CYR61 expression in vivo, we constructed mouse models by tail vein injection of cancer cells and observing lung metastases. We found that over-expression of FOXA1 up-regulated CYR61 expression levels and increased the metastatic potential of HCT116 cells; Conversely, CYR61 knockdown, in cells ectopically expressing FOXA1, decreased cell metastasis (Fig. 6i, j). The above data suggest that the identified enhancers, which are mediated by FOXA1, can regulate CYR61 expression to enhance the metastasis of colon cancer cells.
Enhancer regulation is an emerging mechanism of cancer progression [41, 42]. However, the functional mapping of cancer-specific enhancers is still in its early stages. Herein, we identify enhancers in the regulation of CYR61 expression in colon cancer and the underlying mechanism. Our findings are as follows: (i) selectively active enhancers of CYR61 are formed in colon cancer, but not in normal colon mucosa, that up-regulate CYR61 expression to promote the migration of colon cancer cells; (ii) FOXA1 up-regulates the expression of CYR61 by recruiting the lysine acetyltransferase CBP to the CYR61 promoter and enhancer regions, increasing enhancer activity and forming a promoter-enhancer loop structure; and (iii) treatment with TPA and VP, compounds previously shown to up- and down-regulate CYR61 expression, respectively, also respectively up- and down-regulate CYR61 expression by affecting enhancer activity.
Active CYR61 enhancers, which are enriched in H3K27Ac, are not present in normal colon tissue and cells, but are found in CRC tissue and colon cancer cell lines (Fig. 2a, b, Additional file 1: Figure S2A), indicating that CYR61 enhancers are activated in cancer cells. A recent report of super-enhancers at the c-MYC oncogene locus, which are not active in normal colon cells but are active in HCT116 cells, represents another example of cancer-specific enhancers . Moreover, different cancer cell lines show different enhancer regulation patterns, as reflected by one active enhancer (Enhancer3) in HCT116 cells, but two independent active enhancers (Enhancer1 and Enhancer3) in RKO cells. There was also a predicted Enhancer 2 (Fig. 2). However, cell experiments did not support Enhancer2 to be an active enhancer of CYR61. As Fig. 2c shows, Enhancer2 is decorated with low levels of H3K27ac and does not loop to the CYR61 promoter (Fig. 2e). Moreover, its activity remains unchanged after TPA and VP stimulation (Fig. 5). Likewise, Enhancer1 in HCT116 cells is also decorated with low levels of H3K27ac (Additional file 1: Figure S2A) and does not form a loop with the CYR61 promoter. Concomitantly, TF enrichment at Enhancer1 is also low (Additional file 1: Figure S3C, F).
Pioneer factor FOXA1 and acetyltransferase CBP play a strategic role in the regulation of CYR61 expression. On one hand, both FOXA1 and CBP are needed to maintain H3K27ac enrichment at the CYR61 promoter and enhancer regions. On the other hand, these two TFs are also required for the maintenance of the promoter-enhancer 3D structure. Decreasing FOXA1 binding reduced CBP enrichment at the enhancer regions. In contrast, decreasing CBP levels did not affect FOXA1 binding. These observations suggest that CBP is recruited after the binding of FOXA1, consistent with FOXA1 being a pioneer factor . We, therefore, infer that FOXA1 functions as a pioneer factor to anchor to chromatin and recruit additional TFs, such as CBP, to establish the open state of chromatin and mediate promoter-enhancer looping.
Although H3K27ac , H3K4me1  and H3K4me3  have been reported to be associated with gene transcription and enhancer activity, in our case, H3K27ac is the critical epigenetic modification involved in CYR61 regulation. In both the up-regulation of CYR61 by TPA stimulation or down-regulation of CYR61 by VP stimulation, only H3K27ac enrichment at CYR61 promoter and enhancer regions shows statistical changes that parallel CYR61 expression levels.
Our clinical data demonstrate that CYR61 expression levels are up-regulated in primary colonic adenocarcinoma tissue and positively correlate with tumor grade. Migration experiments and mouse metastasis experiments show that CYR61 promotes colon cancer cells migration and metastasis. Similar observations have been documented in earlier studies . By increasing enhancer activity, FOXA1 can positively regulate CYR61 expression to promote colon cancer cell migration and metastasis.
Based on our current study, we present a schematic model to illustrate the role of enhancers in the regulation of CYR61 expression (Fig. 7). In this model, abnormal enhancer activity (Enhancer3 in HCT116, and Enhancer1 and Enhancer3 in RKO cells) of CYR61 occurs in CRC and up-regulates the expression of CYR61. Based on the crucial role of Mediator in the organization of genomic DNA topological structure  and a recent study showing that FOXA1 is essential for the recruitment of Mediator , we speculate that, under our conditions, the Mediator complex is likely to contribute to the role of FOXA1 and CBP in long-range looping. By mediating binding to the Mediator complex, FOXA1 and CBP bring the CYR61 promoter and enhancers to the same transcriptional foci to enhance CYR61 transcription in an active chromatin region. However, the molecular mechanisms by which FOXA1 is activated and induced to bind to DNA elements in specific colon cancer subtypes remain to be fully understood.
In summary, our study reveals previously undocumented enhancers that are not active in normal colon mucosa, but become activated in CRC. This CRC-related selective enhancer regulation suggests that inhibition of recruitment of FOXA1 and/or CBP to enhancers and the promoter of CYR61 may have therapeutic implications in CRC.
chromosome conformation capture assay
forkhead box A1
Gene Expression Omnibus
acetylated H3 lysine 27
monomethylated H3 lysine 4
quantitative real-time PCR
Jun J, Lau LF. Taking aim at the extracellular matrix: CCN proteins as emerging therapeutic targets. Nat Rev Drug Discov. 2011;10(12):945–63.
Kim K, Chen C, Alpini G, Lau LF. CCN1 induces hepatic ductular reaction through integrin αvβ5–mediated activation of NF-κB. J Clin Investig. 2015;125(5):1886–900.
Kuonen F, Secondini C, Ruegg C. Molecular pathways: emerging pathways mediating growth, invasion, and metastasis of tumors progressing in an irradiated microenvironment. Clin Cancer Res. 2012;18(19):5196–202.
Wu G, Zhu Y, Zhang J. Sox4 up-regulates Cyr61 expression in Colon Cancer cells. Cell Physiol Biochem. 2014;34(2):405–12.
Jeong D, Heo S, Sung Ahn T, Lee S, Park S, Kim H, et al. Cyr61 expression is associated with prognosis in patients with colorectal cancer. BMC Cancer. 2014;14(1):164.
Ladwa R, Pringle H, Kumar R, West K. Expression of CTGF and Cyr61 in colorectal cancer. J Clin Pathol. 2010;64(1):58–64.
Huang X, Xiang L, Li Y, Zhao Y, Zhu H, Xiao Y, et al. Snail/FOXK1/Cyr61 signaling Axis regulates the epithelial–mesenchymal transition and metastasis in colorectal Cancer. Cell Physiol Biochem. 2018;47(2):590–603.
Yan J, Yang B, Lin S, Xing R, Lu Y. Downregulation of miR-142-5p promotes tumor metastasis through directly regulating CYR61 expression in gastric cancer. Gastric Cancer. 2018:1–12.
Tanaka F, Rizqiawan A, Higashikawa K, Tobiume K, Okui G, Shigeishi H, et al. Snail promotes Cyr61 secretion to prime collective cell migration and form invasive tumor nests in squamous cell carcinoma. Cancer Lett. 2013;329(2):243–52.
Shen H, Xu W, Guo R, Rong B, Gu L, Wang Z, et al. Suppression of enhancer Overactivation by a RACK7-histone demethylase complex. Cell. 2016;165(2):331–42.
Schuijers J, Manteiga JC, Weintraub AS, Day DS, Zamudio AV, Hnisz D, et al. Transcriptional dysregulation of MYC reveals common enhancer-docking mechanism. Cell Rep. 2018;23(2):349–60.
Li R, Grimm SA, Chrysovergis K, Kosak J, Wang X, Du Y, et al. Obesity, rather than diet, drives epigenomic alterations in colonic epithelium resembling cancer progression. Cell Metab. 2014;19(4):702–11.
Pelish HE, Liau BB, Nitulescu II, Tangpeerachaikul A, Poss ZC, Da Silva DH, et al. Mediator kinase inhibition further activates super-enhancer-associated genes in AML. Nature. 2015;526(7572):273–6.
Thomas-Claudepierre AS, Robert I, Rocha PP, Raviram R, Schiavo E, Heyer V, et al. Mediator facilitates transcriptional activation and dynamic long-range contacts at the IgH locus during class switch recombination. J Exp Med. 2016;213(3):303–12.
Lai F, Gardini A, Zhang A, Shiekhattar R. Integrator mediates the biogenesis of enhancer RNAs. Nature. 2015;525(7569):399–403.
Calo E, Wysocka J. Modification of enhancer chromatin: what, how, and why? Mol Cell. 2013;49(5):825–37.
Li W, Notani D, Rosenfeld MG. Enhancers as non-coding RNA transcription units: recent insights and future perspectives. Nat Rev Genet. 2016;17(4):207–23.
Bose DA, Donahue G, Reinberg D, Shiekhattar R, Bonasio R, Berger SL. RNA binding to CBP stimulates histone acetylation and transcription. Cell. 2017;168(1–2):135–149.e22.
Seumois G, Chavez L, Gerasimova A, Lienhard M, Omran N, Kalinke L, et al. Epigenomic analysis of primary human T cells reveals enhancers associated with TH2 memory cell differentiation and asthma susceptibility. Nat Immunol. 2014;15(8):777–88.
Raisner R, Kharbanda S, Jin L, Jeng E, Chan E, Merchant M, et al. Enhancer activity requires CBP/P300 Bromodomain-dependent histone H3K27 acetylation. Cell Rep. 2018;24(7):1722–9.
Dancy BM, Cole PA. Protein lysine acetylation by p300/CBP. Chem Rev. 2014;115(6):2419–52.
Wang X, Srivastava Y, Jankowski A, Malik V, Wei Y, Del RR, et al. DNA-mediated dimerization on a compact sequence signature controls enhancer engagement and regulation by FOXA1. Nucleic Acids Res. 2018;46(11):5470–86.
Roe J, Hwang C, Somerville TDD, Milazzo JP, Lee EJ, Da Silva B, et al. Enhancer Reprogramming Promotes Pancreatic Cancer Metastasis. Cell. 2017;170(5):875–888.e20.
Lee H, Chung J, Youn S, Kim J, Park K, Koo B, et al. Forkhead transcription factor FOXO3a is a negative regulator of Angiogenic immediate early gene CYR61, leading to inhibition of vascular smooth muscle cell proliferation and Neointimal hyperplasia. Circ Res. 2007;100(3):372–80.
Hong Y, Downey T, Eu KW, Koh PK, Cheah PY. A ‘metastasis-prone’ signature for early-stage mismatch-repair proficient sporadic colorectal cancer patients and its implications for possible therapeutics. Clin Exp Metastasis. 2010;27(2):83–90.
Skrzypczak M, Goryca K, Rubel T, Paziewska A, Mikula M, Jarosz D, et al. Modeling oncogenic signaling in Colon tumors by multidirectional analyses of microarray data directed for maximization of analytical reliability. PLoS One. 2010;5(10):e13091.
Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9(4):357–9.
Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008;9(9):R137.
Li Q, Su Z, Xu X, Liu G, Song X, Wang R, et al. AS1DHRS4, a head-to-head natural antisense transcript, silences the DHRS4 gene cluster in cis and trans. Proc Natl Acad Sci. 2012;109(35):14110–5.
Whitaker JW, Nguyen TT, Zhu Y, Wildberg A, Wang W. Computational schemes for the prediction and annotation of enhancers from epigenomic assays. Methods. 2015;72:86–94.
Zhu Y, Chen Z, Zhang K, Wang M, Medovoy D, Whitaker JW, et al. Constructing 3D interaction maps from 1D epigenomes. Nat Commun. 2016;7:10812.
Hagege H, Klous P, Braem C, Splinter E, Dekker J, Cathala G, et al. Quantitative analysis of chromosome conformation capture assays (3C-qPCR). Nat Protoc. 2007;2(7):1722–33.
Yang Y, Su Z, Song X, Liang B, Zeng F, Chang X, et al. Enhancer RNA-driven looping enhances the transcription of the long noncoding RNA DHRS4-AS1, a controller of the DHRS4 gene cluster. Sci Rep. 2016;6:20961.
Sun F, Chronis C, Kronenberg M, Chen X, Su T, Lay FD, et al. Promoter-enhancer communication occurs primarily within insulated neighborhoods. Mol Cell. 2018.
Thomas-Chollier M, Hufton A, Heinig M, O'Keeffe S, Masri NE, Roider HG, et al. Transcription factor binding predictions using TRAP for the analysis of ChIP-seq data and regulatory SNPs. Nat Protoc. 2011;6(12):1860–9.
Tsai M, Bogart DF, Li P, Mehmi I, Lupu R. Expression and regulation of Cyr61 in human breast cancer cell lines. Oncogene. 2002;21(6):964–73.
Al-Moujahed A, Brodowska K, Stryjewski TP, Efstathiou NE, Vasilikos I, Cichy J, et al. Verteporfin inhibits growth of human glioma in vitro without light activation. Sci Rep. 2017;7(1):7602.
Liang N, Zhang C, Dill P, Panasyuk G, Pion D, Koka V, et al. Regulation of YAP by mTOR and autophagy reveals a therapeutic target of tuberous sclerosis complex. J Exp Med. 2014;211(11):2249–63.
Wang G, Gu J, Gao Y. MicroRNA target for MACC1 and CYR61 to inhibit tumor growth in mice with colorectal cancer. Tumor Biol. 2016;37(10):13983–93.
Liu R, Huang S, Lei Y, Zhang T, Wang K, Liu B, et al. FGF8 promotes colorectal cancer growth and metastasis by activating YAP1. Oncotarget. 2015;6(2):935–52.
Sur I, Taipale J. The role of enhancers in cancer. Nat Rev Cancer. 2016;16(8):483–93.
Kovar H. Selective enhancer changes in osteosarcoma lung metastasis. Nat Med. 2018;24(2):126–7.
Hnisz D, Schuijers J, Lin CY, Weintraub AS, Abraham BJ, Lee TI, et al. Convergence of developmental and oncogenic signaling pathways at transcriptional super-enhancers. Mol Cell. 2015;58(2):362–70.
Ye Z, Chen Z, Sunkel B, Frietze S, Huang THM, Wang Q, et al. Genome-wide analysis reveals positional-nucleosome-oriented binding pattern of pioneer factor FOXA1. Nucleic Acids Res. 2016;44(16):7540–54.
Spitz F, Furlong EE. Transcription factors: from enhancer binding to developmental control. Nat Rev Genet. 2012;13(9):613–26.
Lee J, Wang C, Xu S, Cho Y, Wang L, Feng X, et al. H3K4 mono- and di-methyltransferase MLL4 is required for enhancer activation during cell differentiation. eLife. 2013;2:e01503.
Chen K, Chen Z, Wu D, Zhang L, Lin X, Su J, et al. Broad H3K4me3 is associated with increased transcription elongation and enhancer activity at tumor-suppressor genes. Nat Genet. 2015;47(10):1149–57.
Allen BL, Taatjes DJ. The mediator complex: a central integrator of transcription. Nat Rev Mol Cell Biol. 2015;16(3):155–66.
Fournier M, Bourriquen G, Lamaze FC, Côté MC, Fournier É, Joly-Beauparlant C, et al. FOXA and master transcription factors recruit mediator and Cohesin to the core transcriptional regulatory circuitry of cancer cells. Sci Rep. 2016;6(1):34962.
We thank Stanley L. Lin for critical reading of the manuscript and Guanwu Li, Xiaofen Zhan, Huiying Zhang, Zhongjing Su, Yuchun Lin, Xina Xie, Yingying Yang and Xiaoting Hong for technological and methodology expertise.
This work was supported by the National Natural Science Foundation of China (NSFC) (31770876, 81470441).
Availability of data and materials
The microarray datasets analyzed during the current study are available in the oncomine repository, https://www.oncomine.org/resource/main.html#a%3A6873%3Bd%3A156636578%3Bdso%3AgeneOverex%3Bdt%3ApredefinedClass%3Bec%3A%5B2%5D%3Bepv%3A150001.151078%2C2937%2C3508%3Bet%3Aover%3Bf%3A1896490%3Bg%3A3491%3Bp%3A200008978%3Bpg%3A1%3Bpvf%3A3020%2C9809%2C150004%3Bscr%3Adatasets%3Bss%3Aanalysis%3Bth%3Ag10.0%2Cp0.001%2Cfc2.0%3Bv%3A18.
The ChIP-seq datasets analyzed during the current study are available in the GEO repository, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM915331, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM883681,
https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM883660, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM956024, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM883688, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM2719751, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM2712746, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM2719748, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM2719749, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM1124062.
Ethics approval and consent to participate
This study was approved by local ethics committees (Shantou University Medical College, China, Permission number: SUMC-2015-42) and patients provided written informed consent for tissue storage and analysis.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Correlation of CYR61 expression of CRC tissues in IHC and clinicopathological parameters. Figure S2. Identification of CYR61 enhancers. Figure S3. Enrichment of FOXA1 and CBP at target loci. Figure S4. TPA and VP regulate CYR61 expression. Table S1. Primers sequences for real-time RT- qPCR. Table S2. Data for ChIP-seq and GRO-seq data downloaded from GEO. Table S3. Primers sequences for ChIP. Table S4. Sequences of siRNA oligonucleotide. Table S5. Plasmids and primers used in dual-luciferase reporter assays. Table S6. Primer sequences for 3C assay. Table S7. Histopathological features and clinical data of the patients. (DOCX 522 kb)