- Open Access
Nrf2 status affects tumor growth, HDAC3 gene promoter associations, and the response to sulforaphane in the colon
Clinical Epigenetics volume 7, Article number: 102 (2015)
The dietary agent sulforaphane (SFN) has been reported to induce nuclear factor erythroid 2 (NF-E2)-related factor 2 (Nrf2)-dependent pathways as well as inhibiting histone deacetylase (HDAC) activity. The current investigation sought to examine the relationships between Nrf2 status and HDAC expression in preclinical and translational studies.
Wild type (WT) and Nrf2-deficient (Nrf2−/+) mice were treated with the colon carcinogen 1,2-dimethylhydrazine (DMH) followed by 400 ppm SFN in the diet (n = 35 mice/group). WT mice were more susceptible than Nrf2−/+ mice to tumor induction in the colon. Tumors from WT mice had higher HDAC levels globally and locally on genes such as cyclin-dependant kinase inhibitor 2a (Cdkn2a/p16) that were dysregulated during tumor development. The average tumor burden was reduced by SFN from 62.7 to 26.0 mm3 in WT mice and from 14.6 to 11.7 mm3 in Nrf2−/+ mice. The decreased antitumor activity of SFN in Nrf2−/+ mice coincided with attenuated Cdkn2a promoter interactions involving HDAC3. HDAC3 knockdown in human colon cancer cells recapitulated the effects of SFN on p16 induction. Human subjects given a broccoli sprout extract supplement (200 μmol SFN equivalents), or reporting more than five cruciferous vegetable servings per week, had increased p16 expression that was inversely associated with HDAC3 in circulating peripheral blood mononuclear cells (PBMCs) and in biopsies obtained during screening colonoscopy.
Nrf2 expression varies widely in both normal human colon and human colon cancers and likely contributes to the overall rate of tumor growth in the large intestine. It remains to be determined whether this influences global HDAC protein expression levels, as well as local HDAC interactions on genes dysregulated during human colon tumor development. If corroborated in future studies, Nrf2 status might serve as a biomarker of HDAC inhibitor efficacy in clinical trials using single agent or combination modalities to slow, halt, or regress the progression to later stages of solid tumors and hematological malignancies.
Histone deacetylase (HDAC) enzymes have emerged as important regulators of cancer development [1, 2]. Downregulation of specific HDACs can increase global histone acetylation, turn-on epigenetically silenced genes, and trigger cell cycle arrest, apoptosis, or differentiation in cancer cells [3–5]. Pan-HDAC inhibitors are currently undergoing clinical evaluation as anticancer agents, but the quest continues for more specific HDAC inhibitors with improved efficacy towards hematological and solid tumors .
We reported that a natural compound, sulforaphane (SFN), targets HDAC3 for protein turnover in human colon cancer cells [6–12]. SFN is obtained from broccoli and other cruciferous vegetables that are rich sources of the precursor, glucoraphanin [13, 14]. SFN was first identified as an inducer of phase 2 detoxification enzymes, acting via the nuclear factor erythroid 2 (NF-E2)-related factor 2 (Nrf2) signaling pathway . However, Nrf2 in certain circumstances can preserve rather than attenuate cancer phenotypes [16, 17].
The current investigation sought to examine the relationships between Nrf2 status and HDAC expression in a widely used model of colon carcinogenesis [18, 19]. Preclinical experiments included post-initiation SFN treatment and highlighted a role for HDAC3 in regulating cyclin-dependant kinase inhibitor 2a (Cdkn2a/p16) expression. The findings were extended to human subjects according to their cruciferous vegetable consumption or supplement intake.
Dosing schedule and HDAC3 levels dictate anticancer outcomes in the colon
1,2-Dimethylhydrazine (DMH) was administered for 10 weeks, and 1 week later SFN was given in the diet by continuous or alternating daily dosing schedules (Fig. 1a). There was a significant reduction in tumor multiplicity and tumor burden only after continuous SFN treatment (Fig. 1b, c). In colon tumors, daily SFN lowered HDAC activity (Fig. 2a) and HDAC3 protein expression (Fig. 2b, c) while increasing global histone H4 acetylation (Fig. 2b, c). The findings in vivo recapitulate prior observations on HDAC3 protein turnover by SFN in cell-based assays . Several SFN metabolites have been examined in the context of the HDAC turnover mechanism [11, 20, 21], namely, SFN-glutathione (SFN-GSH), SFN-cysteine-glycine (SFN-CG), SFN-cysteine (SFN-Cys), and SFN-n-acetylcysteine (SFN-NAC). These metabolites were detected, as reported , in tissues of SFN-treated mice (data not shown).
Nrf2 genetic background influences HDAC protein levels in colon tumors
Continuous SFN treatment was used in a modified protocol that employed fewer DMH doses and a longer post-initiation phase in wild type (WT, Nrf2+/+) and Nrf2−/+ mice (Fig. 3a). The longer study duration of 35 weeks resulted in an average tumor burden of 62.7 mm3 in WT mice (Fig. 3b and Table 1, part B), compared with 18.4 mm3 at 25 weeks (Table 1, part A). Notably, in the 35-week study, DMH produced an average tumor burden of 14.6 mm3 in Nrf2−/+ mice (Fig. 3b and Table 1, part B), compared with 62.7 mm3 in the Nrf2+/+ controls. In SFN-fed mice, the average tumor burden was reduced from 62.7 to 26.0 mm3 in WT mice and from 14.6 to 11.7 mm3 in Nrf2−/+ mice (Fig. 3b and Table 1, part B). In colon tumors, Nrf2 levels were significantly lower in Nrf2−/+ vs. Nrf2+/+ mice, regardless of SFN treatment (Fig. 3c).
Immunoblotting of tissue lysates from tumor and adjacent normal colon revealed that HDAC3 expression was reduced by SFN treatment in WT mice (Fig. 4a), but not in mice on the Nrf2−/+ background (Fig. 4b). In general, HDAC expression was lower in colon tissues of Nrf2−/+ animals (compare y-axes of densitometry plots in Fig. 4a, b). No treatment-related differences were detected for HDAC3 mRNA levels or for other class I HDACs (data not shown).
Nrf2 status impacts HDAC3 levels on p16 in mouse colon tumors
Gene expression arrays differentiated between the genes most altered by SFN treatment when comparing tumor with adjacent normal colon (Fig. 5a). For the complete list of genes and the respective fold changes, see Additional file 1: Table S1. Scatter plots were generated with an arbitrary cutoff of fivefold in either direction. When tumors were compared with adjacent normal colon, cyclin-dependant kinase inhibitor 2A (Cdkn2a/P16), a well-known tumor suppressor, was surprisingly the most highly overexpressed gene in the tumors of WT mice (Fig. 5b). Estrogen receptor 1, alpha (Esr1), a nuclear hormone receptor involved in the regulation of eukaryotic gene expression that affects cellular proliferation and differentiation in target tissues, also was highly expressed, whereas known tumor suppressors cyclin-dependant kinase 1A (Cdkn1a/p21), serine peptidase inhibitor (serpin b5), and oncogenes such as Met proto-oncogene (Met), and protein kinase C, alpha (Prkca) were underexpressed in the tumor. On the Nrf2−/+ background (Fig. 5c), tumor suppressor genes such as Cdkn2a and transformation-related protein 53 (Tp53), as well as reported oncogenic factors E26 avian leukemia oncogene 1,5′ domain (Ets1) and insulin-like growth factor 2 receptor (Igf2r) were overexpressed, whereas underexpressed genes included tumor suppressor genes cyclin-dependent kinase inhibitor 2B (Cdkn2b/p15) and serpin b5. mutL homolog 1 (Mlh1), which is known to be silenced in human colon cancer was also attenuated in mouse colon tumors. SFN treatment compressed the spread of altered genes (Fig. 5d, e), although Cdkn2a/p16 levels remained consistently elevated, especially in WT mice given SFN.
Candidates from the PCR arrays were validated by quantitative reverse transcription PCR (qRT-PCR), confirming Cdkn2a/p16 as being among the most highly overexpressed genes (Additional file 2: Figure S1). SFN altered Cdkn2a/p16 mRNA expression in colon tumor and adjacent normal-looking colon of WT and Nrf2−/+ mice (compare red arrows in Additional file 2: Figure S1, A vs. B, C vs. D, E vs. F, and G vs. H). Interestingly, SFN increased Cdkn2a/p16 mRNA levels in colon tumors of WT mice but had the opposite effect in Nrf2−/+ mice (Fig. 6a). These findings were supported by immunohistochemical analysis of p16 protein expression (Fig. 6b–e). Thus, high levels of p16 protein were detected in tumors compared with adjacent normal colon, especially in WT mice fed with SFN (Fig. 6c, inset).
Based on prior reports linking Cdkn2a and HDAC3 [23–25], we performed chromatin immunoprecipitation (ChIP) assays in vivo and observed HDAC3 interactions on the p16 proximal promoter region to be lower in colon tumors of Nrf2−/+ mice compared with WT (Fig. 6f, black vs. grey bars), especially after SFN treatment (green bars). A runt-related transcription factor 1 (Runx1) binding site, ~10 Kb upstream of the p16 promoter , had minimal HDAC3 interactions and served as a negative control for the ChIP assays (Fig. 6f, region “4”). Notably, the findings for HDAC3-dependent regulation of p16 did not represent a generic response of all highly dysregulated genes; HDAC3 interactions on Esr1, for example, were unaffected by SFN (Additional file 2: Figure S2).
p16 is regulated directly by HDAC3, but not Nrf2, in human colon cancer cells
In human colon cancer cells, Cdkn2a mRNA expression was induced 24 h after SFN treatment (Fig. 7a), and immunoblotting corroborated the increased expression of p16 protein (Fig. 7b). RNAi-mediated knockdown of HDAC3 also increased p16 levels, whereas the knockdown of Kelch-like ECH-associated protein 1 (Keap1), which releases Nrf2 from the cytoplasm to the nucleus for gene activation, did not increase p16 expression. Keap1 knockdown induced the Nrf2 target gene heme oxygenase-1 (HO1, a positive control for Nrf2 activation) to a similar extent as SFN treatment, whereas HDAC3 knockdown had no effect on HO1 (Fig. 7c). Knockdown of Keap1 and HDAC3 was confirmed at the mRNA and protein level (Fig. 7d, e). We conclude that Cdkn2a/p16 is an HDAC3-regulated gene that does not involve direct Nrf2 interactions, whereas the reverse is true for HO1.
HDAC3 and p16 are reciprocally regulated in humans after SFN intake
To examine these relationships in humans, healthy subjects consumed a broccoli sprout extract (BSE) supplement or placebo once a day for 7 days, and blood was drawn at 0, 1, 3, and 6 h post-consumption on day 1 (first day of study) and day 7 (to study effects due to dose accumulation, if any), and on days 8, 9, and 14 as a follow-up to ensure systemic clearance of SFN and metabolites (Fig. 8a). After consuming BSE (200 μmol SFN equivalents), plasma levels of SFN metabolites peaked at 1–3 h, rapidly decreased at 6 h, and were undetectable at 24 h (Fig. 8b, black symbols). The specific SFN metabolites were in agreement with prior studies  and were higher on day 1 as compared to day 7 (Additional file 3: Table S2). In a recent report , half the dose (100 μmol SFN equivalents) given every 12 h produced lower SFN tissue concentrations than in the present investigation. We infer that a single daily dose of 200 μmol SFN equivalents might be necessary to elicit HDAC inhibitory responses in vivo.
In mice, a single oral dose of SFN reduced HDAC3 protein expression and increased p16 protein level in splenocytes, supporting the use of these end points as potential biomarkers in systemic tissues (Additional file 2: Figure S3). Thus, HDAC3 and p16 protein expression changes also were examined in circulating peripheral blood mononuclear cells (PBMCs) from human volunteers (Fig. 8c). HDAC3 was reduced as early as 1 h after BSE consumption and continued to remain low up to 6 h later. Notably, the loss of HDAC3 coincided with increased p16 expression during this time period. On day 7, HDAC3 levels returned to baseline, whereas p16 remained elevated (Fig. 8c). This implies that cross-talk between HDAC3 and p16 can be uncoupled at later times, possibly via interactions of p16 with alternative HDACs, histone acetyltransferases, and/or their co-regulators [20, 21]. For example, HDAC3/SMRT inhibition and turnover has been observed to precede changes in other HDACs, such as HDAC6 . Lower SFN metabolite levels at day 7 vs. day 1 hinted at a possible compensatory mechanism following repeated daily SFN intake. Further studies are needed to examine the possible induction of pathways that favor enhanced SFN metabolism and/or excretion.
Based on findings from the short-term intervention trial with BSE, HDAC3 and p16 were examined in the context of more typical human dietary intake patterns (Fig. 8d). Patients scheduled for a screening colonoscopy were stratified as high vs. low cruciferous vegetable consumers based on a validated questionnaire . In PBMCs obtained immediately prior to colonoscopy, p16 mRNA levels were significantly higher for subjects reporting >5 vs. 0–1 servings of cruciferous vegetables per week (Fig. 8e). We next examined the corresponding colon biopsies for selected proteins of interest, i.e., p16, HDACs, and histone acetylation (Fig. 8f). With higher cruciferous vegetable intake, histone acetylation and p16 expression were increased, whereas HDAC1 levels were unchanged. Based on densitometry measurements of immunoblots, a significant inverse association was observed for p16 and HDAC3 in colon biopsies (Fig. 8g), whereas p16 and acetyl histone H4 lysine 12 (AcH4K12) (Fig. 8h) and HDAC3 and Nrf2 (Fig. 8i) were positively correlated.
SFN was first identified as a potent inducer of phase 2 enzymes, acting via the Nrf2 pathway to induce detoxification pathways that favor carcinogen excretion and elimination from the body [14, 15]. However, SFN also can be effective in post-initiation protocols , independent of carcinogen exposure. This is exemplified by the tumor suppression observed for SFN in genetic models, such as the adenomatous polyposis coli/multiple intestinal neoplasia (ApcMin/+) mouse, in which HDAC inhibition was identified as a contributing mechanism . Therefore, we examined the interplay between two key mechanisms implicated in the SFN antitumor activity, namely, Nrf2 induction and HDAC inhibition. Although pharmacological HDAC inhibitors have been tested in preclinical models of colon cancer , the role of Nrf2 was not examined. We specifically sought to test the hypothesis that Nrf2 status might affect HDAC3 protein expression levels in colon tumors, and thus the inhibitory response to SFN acting preferentially on HDAC3 [20, 21].
As in the ApcMin/+ mouse [10, 31], SFN suppressed tumorigenesis in the DMH model, and this was accompanied by reduced HDAC activity and HDAC3 protein expression in the colon tumors. These findings are in accordance with data from human colon cancer cells showing HDAC inhibition by SFN and its metabolites [11, 20], and by other dietary isothiocyanates , to involve HDAC3 protein turnover. Metabolites implicated in the HDAC3 turnover mechanism [11, 20, 21], such as SFN-Cys and SFN-NAC, were detected in tissues of SFN-treated mice, as reported . The corroboration of HDAC3 as a target of SFN in vivo is important considering the critical role of this HDAC in regulating colon cancer growth and tumorigenesis . Loss of tumor suppression in mice fed with SFN on alternating days (Fig. 1) might be related to the inability to sustain high enough SFN metabolite levels for effective HDAC3 inhibition. Indeed, SFN metabolites are cleared within 24 h in mice , as in human subjects consuming BSE (Fig. 8b).
Nrf2-deficient mice are generally more sensitive to carcinogens and agents that trigger chronic inflammation [32–35]. However, Nrf2−/+ mice treated with DMH and observed for up to 35 weeks had a significantly reduced tumor burden and lower HDAC protein expression compared with WT animals. The differential response to the carcinogen likely was not attributable to genes that regulate DMH metabolism (Cyp2E1) or DNA repair (Mgmt), since their expression was similar in WT and Nrf2−/+ mouse colon (Additional file 2: Figure S4). In a recent study using urethane to initiate lung tumors , resistance to tumor growth was observed in Nrf2-deficient mice compared with WT. The enhanced tumor growth in WT mice adds to the discussion on pros vs. cons of Nrf2 signaling in different stages of cancer development [17, 36].
From the gene expression arrays, qRT-PCR and immunoblotting experiments, tumor suppressors p21 and p15 were induced only marginally, if at all, by SFN in the preclinical model reported here, in marked contrast to p16. For example, in gene expression arrays (Additional file 1: Table S1), p21 and p15 were attenuated slightly in the colon tumors from SFN-treated mice, whereas p16 was induced greater than 20-fold. Although p16 more typically is considered a tumor suppressor protein, high Cdkn2a/p16 levels have been detected in mouse colon tumors induced by azoxymethane, a metabolite of DMH , and in human benign tumors and high-grade malignancies . Overexpression of p16 has been linked to the so-called oncogene-induced senescence (OIS) in benign tumors or as a mechanism to arrest uncontrolled proliferation in advanced cancers . Early upregulation of p16 in some tumors might represent an attempt to correct for one or more dysregulated signaling pathways. We focused on p16 as a major target dysregulated in DMH-induced colon tumors and noted that SFN increased p16 expression in tumors of WT mice but had the opposite effect in Nrf2−/+ mice. Based on evidence that HDAC inhibition activates p16 [39–41], we confirmed HDAC3 interactions on p16 to be higher in colon tumors of WT vs. Nrf2−/+ mice (Fig. 6f), and inversely associated with p16 mRNA levels (Fig. 6a). In human colon cancer cells, HDAC3 knockdown increased p16 levels to a similar extent as SFN treatment, whereas Keap1 knockdown had no effect on p16 (Fig. 7). This suggested that p16 is regulated by HDAC3 but is not a direct target of Nrf2.
What, then, connects Nrf2 genetic background to altered HDAC3 levels on p16? A working model can be proposed (Additional file 2: Figure S5). We speculate that Nrf2 deficiency in mice, through mechanisms that remain to be clarified, attenuates the overall rate of colon tumor growth and global HDAC levels within the tumor. This in turn diminishes HDAC interactions on key genes dysregulated during tumor development, regardless of whether or not they are directly regulated by Nrf2 binding. Mechanisms affecting promoter methylation and transcription factor access [42–45] might influence which genes are most altered, and the ultimate response to SFN treatment.
The BSE supplement is being evaluated in several human trials (clinicaltrials.gov). In the current investigation, SFN metabolites were detected in plasma shortly after BSE intake, and in circulating PBMCs, there was increased p16 and transiently decreased HDAC3 expression. As predicted from the preclinical model, in human colon biopsies HDAC3 was inversely associated with p16 (Fig. 8g), and positively correlated with Nrf2 (Fig. 8i). We were interested in the variability of Nrf2 expression in normal colon biopsies (Fig. 8f). In The Cancer Genome Atlas (TCGA) database , NFE2L2 (Nrf2) levels also varied markedly in normal colon and in colon tumors and were inversely correlated with Cdkn2a/p16 (Fig. 9), in accordance with data from cell-based, preclinical, and translational studies reported here. We are now recruiting additional human volunteers to obtain both polyps and adjacent normal colon biopsies, seeking to corroborate the proposed cross-talk between Nrf2, HDAC3, and p16.
Nrf2 status in the mouse colon appears to serve as an arbiter of overall colon tumor growth, HDAC protein expression in colon tumors, and the response to HDAC inhibitor treatment mediated by downstream molecular targets such as p16. A key issue will be the extent to which Nrf2 status influences HDAC levels and p16 expression at different stages of human colon cancer development and the ultimate response to pan-HDAC or HDAC-selective inhibitors. This could have implications beyond the treatment of colorectal cancer, for example, in other solid tumors and hematological malignancies currently undergoing clinical evaluation with HDAC inhibitors . In the broad context of these various clinical trials and their overall aims, we believe that Nrf2 status is worthy of further investigation as a possible determinant of tumor growth and HDAC inhibitor responsiveness.
Animals and diets
Male WT or Nrf2−/+ mice at 8–10 weeks of age were randomized to 25 mice/group (pilot study) or 35 mice/group (main study). Tumors were induced by 1,2-dimethylhydrazine (DMH, Sigma-Aldrich), as reported . After DMH treatment, mice continued on regular AIN93 diet for an additional week before administering AIN93 diet or AIN93 diet supplemented with 400 ppm d,l-SFN (Toronto Research Chemicals, Inc.), either continuously or on alternating days. At end of the study, the colon was removed, and tumors were scored for number, size, and position by individuals blinded to the treatment. Tumor volume was calculated using the formula xy 2 * 0.5 (x = long diameter, y = short diameter). Total tumor burden per animal was calculated by adding the individual tumor volumes. Tumor and adjacent normal-looking tissue was removed and one portion was fixed in 10 % buffered formalin, while the other portion was flash-frozen in liquid nitrogen and stored at −80 °C. The work was approved by the Institutional Animal Care and Use Committee.
Frozen samples of colon tumors and adjacent tissue were thawed and subjected to immunoblotting using the methodology described previously [20, 21, 48]. Antibodies were for HDAC1 and HDAC3 (Santa Cruz), p16 (Proteintech), acetyl histone H4K12 and histone H4 (Cell Signaling), and β-actin (Sigma-Aldrich).
PCR arrays and qPCR
RT2 Profiler arrays were run as per manufacturer’s instructions (SA Biosciences, Qiagen, Valencia, CA, USA) on colon tumor samples (pooled, n = 6) and matched controls (pooled, n = 6). From the threshold cycle (C t) value, the relative gene expression of each target was normalized to human β-actin gene (ACTB) or murine β-actin gene (Actb) (β-actin gene in human and mouse, respectively). Two or more separate qRT-PCR experiments were performed to validate targets of interest, as reported .
The ChIP-IT Express Enzymatic kit (Active Motif, Carlsbad, CA) was used, as reported . Frozen tissue (100 mg) was cut into pieces, cross-linked with formaldehyde, and homogenized in order to isolate the nuclear fraction. DNA fragmentation was performed via enzymatic shearing, using a proprietary cocktail (Active Motif) that randomly cleaves between nucleosomes, generating DNA fragments of ~500 bp. Ten microliters of fragmented chromatin was kept as input while the remaining was immunoprecipitated (IP) with anti-HDAC3 antibody (Santa Cruz). After reversing the cross-linking and proteinase treatment, DNA was purified using QIAquick PCR Purification kits (Qiagen). PCR was run on a Roche LightCycler 480 II with preincubation for 5 min at 95 °C, then 45 cycles at 95 °C for 10 s, 60 °C for 10 s, and 72 °C for 10 s. Each experiment was repeated at least twice.
Formalin-fixed paraffin-embedded mouse colon tumor and adjacent normal tissue were processed for immunohistochemistry as reported . Slides containing 5 μm sections were rehydrated and placed in an Autostainer (Dako). After primary antibody to p16 (Proteintech) for 30 min and One-Step HRP Polymer anti-IgG (ImmunoBioscience) for 7 min, Nova Red (Vector Labs) was applied for 5 min followed by hematoxylin (Dako). Images were acquired on a Nikon E400 microscope equipped with a CCD camera.
HCT116 cells were obtained from American Type Culture Collection (Manassas, VA, USA) and validated as reported . Cells were transfected for 48 h with HDAC3 siRNA (Trilencer-27, OriGene), Keap1 siRNA (Sigma-Aldrich), control siRNA (OriGene), or Lipofectamine 2000 alone, using the manufacturer’s protocol (Invitrogen).
Phase “0” trial
A pilot intervention study was performed based on a prior protocol . Ten healthy subjects avoided cruciferous vegetables, starting 1 week before day 1 of the study and continuing through day 14. After an overnight fast, volunteers ate a standardized breakfast together with a broccoli sprout extract (BSE) supplement (IND #111736, 200 μmol SFN equivalents, n = 5) or a placebo (n = 5), once a day for 7 days. BSE and placebo were obtained from Johns Hopkins University (Baltimore, MD, USA), with SFN content validated by LC-MS/MS, as reported . The Institutional Review Board (IRB) approved the protocol, and all participants provided written consent.
Whole blood was collected into ethylene diamine tetraacetic acid (EDTA) Vacutainers (VWR, Radnor, PA, USA) at 0, 1, 3, and 6 h post-consumption on days 1 and 7 and on days 8, 9, and 14. After centrifuging at 2000 rpm for 30 min, plasma was removed, acidified with trifluoroacetic acid, and stored at −80 °C. Samples were analyzed for SFN metabolites as previously described . The remaining whole blood was processed to isolate peripheral blood mononuclear cells (PBMC), as described previously , and frozen at −80 °C.
Screening colonoscopy study
Men and women aged >50 years and scheduled for a screening colonoscopy were recruited based on cruciferous vegetable consumption (n = 28). Recruitment and data collection were implemented through the Oregon Clinical and Translational Research Institute and the Oregon Health & Science University Cancer Institute, with IRB approval and written consent from each participant. Subjects completed a validated  cruciferous vegetable food frequency questionnaire (CVFFQ) and had three 24-h dietary recalls over a 3-week period. The CVFFQ assessed intake over the previous 12 months, relating to number of servings, serving size, intake of raw and cooked vegetables, method of cooking, and use of condiments. Data from the CVFFQ were analyzed by the Arizona Diet, Behavioral, and Quality of Life Assessment Center, University of Arizona, Tucson, AZ, USA. Volunteers were stratified into low (0–1 serving/week, n = 5) and high (≥5 servings/week, n = 23) consumers. Dietary recalls were used to provide information on possible changes in diet between the CVFFQ and the clinic visit. Blood was obtained and processed as in the BSE trial. In addition, two biopsies of rectal colon and two of proximal colon were taken with standard forceps and placed in formalin or flash-frozen in liquid nitrogen.
Results were expressed as mean ± SD. Analysis of variance (ANOVA) was used for group comparisons, followed by Bonferroni’s multiple comparison test (GraphPad Prism v 5.04). Student’s t test was used for paired comparisons, with P < 0.05 considered as significant.
acetyl histone H4 lysine 12
- ACTB :
human β-actin gene
- Actb :
murine β-actin gene
- ApcMin/+ :
adenomatous polyposis coli/multiple intestinal neoplasia
broccoli sprout extract
cyclin-dependant kinase inhibitor 1A
- Cdkn2a/p16 :
cyclin-dependant kinase inhibitor 2A
cyclin-dependent kinase inhibitor 2B
cruciferous vegetable food frequency questionnaire
cytochrome P450 isoform 2E1
ethylene diamine tetraacetic acid
- Esr1 :
estrogen receptor 1
- Ets1 :
E26 avian leukemia oncogene 1,5′ domain
histone deacetylase 1
histone deacetylase 3
- HO1 :
insulin-like growth factor 2 receptor
Kelch-like ECH-associated protein 1
O6-alkylguanine DNA methyltransferase
mutL homolog 1
nuclear factor erythroid 2 (NF-E2)-related factor 2
transformation-related protein 53
peripheral blood mononuclear cells
parts per million
protein kinase C, alpha
runt-related transcription factor 1
- Serpin b5:
serine peptidase inhibitor
The Cancer Genome Atlas
Fredly H, Gjertsen B, Bruserud Ø. Histone deacetylase inhibition in the treatment of acute myeloid leukemia: the effects of valproic acid on leukemic cells, and the clinical and experimental evidence for combining valproic acid with other antileukemic agents. Clin Epigenetics. 2013;5:12.
Mack GS. Epigenetic cancer therapy makes headway. J Natl Cancer Inst. 2006;98:1443–4.
Wilson AJ, Byun DS, Popova N, Murray LB, L'Italien K, Sowa Y, et al. Histone deacetylase 3 (HDAC3) and other class I HDACs regulate colon cell maturation and p21 expression and are deregulated in human colon cancer. J Biol Chem. 2006;281:13548–58.
Marks PA, Xu W-S. Histone deacetylase inhibitors: potential in cancer therapy. J Cell Biochem. 2009;107:600–8.
Delcuve GP, Khan DH, Davie JR. Roles of histone deacetylases in epigenetic regulation: emerging paradigms from studies with inhibitors. Clin Epigenetics. 2012;4:5.
Dashwood RH, Myzak MC, Ho E. Dietary HDAC inhibitors: time to rethink weak ligands in cancer chemoprevention? Carcinogenesis. 2006;27:344–9.
Rajendran P, Ho E, Williams DE, Dashwood RH. Dietary phytochemicals, HDAC inhibition, and DNA damage/repair defects in cancer cells. Clin Epigenetics. 2011;3:4.
Rajendran P, Williams DE, Ho E, Dashwood RH. Metabolism as a key to histone deacetylase inhibition. Crit Rev Biochem Mol Biol. 2011;46:181–99.
Myzak MC, Dashwood RH. Chemoprotection by sulforaphane: keep one eye beyond Keap1. Cancer Lett. 2006;233:208–18.
Myzak MC, Dashwood WM, Orner GA, Ho E, Dashwood RH. Sulforaphane inhibits histone deacetylase in vivo and suppresses tumorigenesis in Apcmin mice. FASEB J. 2006;20:506–8.
Myzak MC, Karplus PA, Chung F-L, Dashwood RH. A novel mechanism of chemoprotection by sulforaphane: inhibition of histone deacetylase. Cancer Res. 2004;64:5767–74.
Myzak MC, Tong P, Dashwood W-M, Dashwood RH, Ho E. Sulforaphane retards the growth of human PC-3 xenografts and inhibits HDAC activity in human subjects. Exp Biol Med. 2007;232:227–34.
Fahey JW, Zhang Y, Talalay P. Broccoli sprouts: an exceptionally rich source of inducers of enzymes that protect against chemical carcinogens. Proc Natl Acad Sci USA. 1997;94:10367–72.
Zhang Y, Talalay P, Cho CG, Posner GH. A major inducer of anticarcinogenic protective enzymes from broccoli: isolation and elucidation of structure. Proc Natl Acad Sci USA. 1992;89:2399–403.
Kensler TW, Egner PA, Agyeman AS, Visvanathan K, Groopman JD, Chen JG, et al. Keap1-nrf2 signaling: a target for cancer prevention by sulforaphane. Top Curr Chem. 2013;329:163–77.
Satoh H, Moriguchi T, Takai J, Ebina M, Yamamoto M. Nrf2 prevents initiation but accelerates progression through the Kras signaling pathway during lung carcinogenesis. Cancer Res. 2013;73:4158–68.
Kensler TW, Wakabayashi N. Nrf2: friend or foe for chemoprevention? Carcinogenesis. 2010;31:90–9.
Corpet DE, Pierre F. How good are rodent models of carcinogenesis in predicting efficacy in humans? A systematic review and meta-analysis of colon chemoprevention in rats, mice and men. Eur J Cancer. 2005;41:1911–22.
Sumiyoshi H, Wargovich MJ. Chemoprevention of 1,2-dimethylhydrazine-induced colon cancer in mice by naturally occurring organosulfur compounds. Cancer Res. 1990;50:5084–7.
Rajendran P, Delage B, Dashwood WM, Yu TW, Wuth B, Williams DE, et al. Histone deacetylase turnover and recovery in sulforaphane-treated colon cancer cells: competing actions of 14-3-3 and Pin1 in HDAC3/SMRT corepressor complex dissociation/reassembly. Mol Cancer. 2011;10:68.
Rajendran P, Kidane AI, Yu TW, Dashwood WM, Bisson WH, Löhr CV, et al. HDAC turnover, CtIP acetylation and dysregulated DNA damage signaling in colon cancer cells treated with sulforaphane and related dietary isothiocyanates. Epigenetics. 2013;8:612–23.
Clarke JD, Hsu A, Williams DE, Dashwood RH, Stevens JF, Yamamoto M, et al. Metabolism and tissue distribution of sulforaphane in Nrf2 knockout and wild-type mice. Pharm Res. 2011;28:3171–9.
Feng Y, Wang X, Xu L, Pan H, Zhu S, Liang Q, et al. The transcription factor ZBP-89 suppresses p16 expression through a histone modification mechanism to affect cell senescence. FEBS J. 2009;276:4197–206.
Wang X, Feng Y, Xu L, Chen Y, Zhang Y, Su D, et al. YY1 restrained cell senescence through repressing the transcription of p16. Biochim Biophys Acta. 1783;2008:1876–83.
Zheng S, Li Q, Zhang Y, Balluff Z, Pan Y-X. Histone deacetylase 3 (HDAC3) participates in the transcriptional repression of the p16 (INK4a) gene in mammary gland of the female rat offspring exposed to an early-life high-fat diet. Epigenetics. 2012;7:183–90.
Linggi B, Müller-Tidow C, van de Locht L, Hu M, Nip J, Serve H, et al. The t(8;21) fusion protein, AML1 ETO, specifically represses the transcription of the p14(ARF) tumor suppressor in acute myeloid leukemia. Nat Med. 2002;8:743–50.
Atwell LL, Hsu A, Wong CP, Stevens JF, Bella D, Yu TW, et al. Absorption and chemopreventive targets of sulforaphane in humans following consumption of broccoli sprouts or a myrosinase-treated broccoli sprout extract. Mol Nutr Food Res. 2015;59:424–33.
Thomson CA, Rock CL, Caan BJ, Flatt SW, Al-Delaimy WA, Newman VA, et al. Increase in cruciferous vegetable intake in women previously treated for breast cancer participating in a dietary intervention trial. Nutr Cancer. 2007;57:11–9.
Chung FL, Conaway CC, Rao CV, Reddy BS. Chemoprevention of colonic aberrant crypt foci in Fischer rats by sulforaphane and phenethyl isothiocyanate. Carcinogenesis. 2000;21:2287–91.
Glauben R, Batra A, Stroh T, Erben U, Fedke I, Lehr HA, et al. Histone deacetylases: novel targets for prevention of colitis-associated cancer in mice. Gut. 2008;57:613–22.
Hu R, Khor TO, Shen G, Jeong WS, Hebbar V, Chen C, et al. Cancer chemoprevention of intestinal polyposis in ApcMin/+ mice by sulforaphane, a natural product derived from cruciferous vegetables. Carcinogenesis. 2006;27:2038–46.
Cheung KL, Lee JH, Khor TO, Wu TY, Li GX, Chan J, et al. Nrf2 knockout enhances intestinal tumorigenesis in Apc(min/+) mice due to attenuation of anti-oxidative stress pathway while potentiates inflammation. Mol Carcinog. 2014;53:77–84.
Khor TO, Huang M-T, Kwon KH, Chan JY, Reddy BS, Kong A-N. Nrf2-deficient mice have an increased susceptibility to dextran sulfate sodium-induced colitis. Cancer Res. 2006;66:11580–4.
Khor TO, Huang MT, Prawan A, Liu Y, Hao X, Yu S, et al. Increased susceptibility of Nrf2 knockout mice to colitis-associated colorectal cancer. Cancer Prev Res. 2008;1:187–91.
Osburn WO, Karim B, Dolan PM, Liu G, Yamamoto M, Huso DL, et al. Increased colonic inflammatory injury and formation of aberrant crypt foci in Nrf2-deficient mice upon dextran sulfate treatment. Int J Cancer. 2007;121:1883–91.
Sporn MB, Liby KT. NRF2 and cancer: the good, the bad, and the importance of context. Nat Rev Cancer. 2012;12:564–71.
Wang QS, Papanikolaou A, Nambiar PR, Rosenberg DW. Differential expression of p16(INK4a) in azoxymethane-induced mouse colon tumorigenesis. Mol Carcinog. 2000;28:139–47.
Romagosa C, Simonetti S, López-Vicente L, Mazo A, Lleonart ME, Castellvi J, et al. p16(Ink4a) overexpression in cancer: a tumor suppressor gene associated with senescence and high-grade tumors. Oncogene. 2011;30:2087–97.
Chen M, Voeller D, Marquez VE, Kaye FJ, Steeg PS, Giaccone G, et al. Enhanced growth inhibition by combined DNA methylation/HDAC inhibitors in lung tumor cells with silenced CDKN2A. Int J Oncol. 2010;37:963–71.
Yamaguchi J, Sasaki M, Sato Y, Itatsu K, Harada K, Zen Y, et al. Histone deacetylase inhibitor (SAHA) and repression of EZH2 synergistically inhibit proliferation of gallbladder carcinoma. Cancer Sci. 2010;101:355–62.
Valentini A, Gravina P, Federici G, Bernardini S. Valproic acid induces apoptosis, p16INK4A upregulation, and sensitization to chemotherapy in human melanoma cells. Cancer Biol Ther. 2007;6:185–91.
Chien WW, Ffrench M. Regulation of p16INK4a, senescence and oncogenesis. Med Sci. 2006;22:865–71.
Ohtani N, Zebedee Z, Huot TJ, Stinson JA, Sugimoto M, Ohashi Y, et al. Opposing effects of Ets and Id proteins on p16INK4a expression during cellular senescence. Nature. 2001;409:1067–70.
Wang W, Pan K, Chen Y, Huang C, Zhang X. The acetylation of transcription factor HBP1 by p300/CBP enhances p16INK4A expression. Nucleic Acids Res. 2012;40:981–95.
Borinstein SC, Conerly M, Dzieciatkowski S, Biswas S, Washington MK, Trobridge P, et al. Aberrant DNA methylation occurs in colon neoplasms arising in the azoxymethane colon cancer model. Mol Carcinog. 2010;49:94–103.
Meyer LR, Zweig AS, Hinrichs AS, Karolchik D, Kuhn RM, Wong M, et al. The UCSC Genome Browser database: extensions and updates 2013. Nucleic Acids Res. 2013;41:D64–9.
Tuma RS. Epigenetic therapies move into new territory, but how exactly do they work? J Natl Cancer Inst. 2009;101(19):1300–1.
Wang R, Dashwood WM, Nian H, Löhr CV, Fischer KA, Tsuchiya N, et al. NADPH oxidase overexpression in human colon cancers and rat colon tumors induced by 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine (PhIP). Int J Cancer. 2011;128:2581–90.
Nian H, Delage B, Pinto JT, Dashwood RH. Allyl mercaptan, a garlic-derived organosulfur compound, inhibits histone deacetylase and enhances Sp3 binding on the P21WAF1 promoter. Carcinogenesis. 2008;29:1816–24.
Parasramka M, Dashwood WM, Wang R, Saeed HH, Williams DE, Ho E, et al. A role for low-abundance miRNAs in colon cancer: the miR-206/Krüppel-like factor 4 (KLF4) axis. Clin Epigenetics. 2012;4:16.
We thank Rong Wang, Hui Nian, Chris Larsen, Kate Cleveland, Hassaan Saeed, Tian-Wei Yu, and Soyoun Ahn for their technical assistance. Daniel Sudakin, M.D., Mary Garrard, R.N., Karen Hardin, Joshua Hay, and Matthew French assisted with the clinical trials. BSE supplement and placebos were from Drs. P. Talalay and J. Fahey (Johns Hopkins University, Baltimore, MD, USA). This work is supported by a grant CA090890 from the National Cancer Institute, P30 grants ES00210 and ES02351 from the National Institute of Environmental Health Sciences, and by a Chancellor’s Research Initiative from Texas A&M University. The open access publishing fees for this article have been covered by the Texas A&M University Online Access to Knowledge (OAK) Fund, supported by the University Libraries and the Office of the Vice President for Research.
The authors declare that they have no competing interests.
PR, EH, DAL, and RHD were responsible for conception and design of the studies. PR, WMD, YK, LL, KAF, and CVL developed the methodology. YY provided the animals, and DAL was responsible for the clinical component, including patient recruitment. Analysis and interpretation of data (statistical analysis, online data mining) were performed by PR, EK, and GJ. Writing, review, and/or revision of the manuscript was performed by PR, WMD, EK, GJ, CVL, DEW, EH, DAL, and RHD. All authors read and approved the final manuscript.
Gene expression changes in mouse colon. The Table shows gene expression data of '84 genes' in the RT2 Profiler array as described in Methods. Each colum (as labeled) provides a list of the different genes, gene names, and the gene expression data (relative to ACTB) normalized to normal colon in WT vehicle control mice calculated using the RT2 Profier PCR Array Data Analysis software. Each row provides data for a single gene, positive values indicating higher expression whereas negative values indicating lower gene expression compared to normal colon in WT vehicle controls. (DOCX 30 kb)
Validation of gene array data in mouse colon tumors and normal colon. a–d Relative mRNA expression of selected genes normalized to the β-actin gene (Actb) in normal colonic mucosa, and e–h in DMH-induced colon tumors. Data = mean ± SD (n = 3). Red arrow, p16. Figure S2. HDAC3 interactions with Esr1 are unaffected by SFN. Primers were designed to interrogate the promoter region adjacent to the transcriptional start site of Esr1, a gene strongly overexpressed in DMH-induced colon tumors (see Fig. 5b). ChIP data = mean ± SD from three independent ChIP assays. Esr1 served as a negative control for the ChIP assays with p16 (see Fig. 6f). Figure S3. SFN altered HDAC3 and p16 expression in mouse splenocytes. Mice (n = 3) were administered a single oral gavage of 200 μmol SFN, and mononuclear cells were isolated from mouse spleen at 6 h, as reported earlier . a Whole cell lysates were immunoblotted for HDAC3, p16, and AcH4K12 with densitometry data normalized to β-actin. b A colon tumor from DMH-treated mice was included as a reference control. Figure S4. No change in the expression of genes involved in DMH activation (Cyp2E1) and DNA repair (Mgmt). Relative mRNA expression levels were normalized to the β-actin gene, Actb. Data = mean ± SD (n = 3). Figure S5. Working model for Nrf2 status impacting tumor growth, HDAC3 levels, and p16 induction. (PPTX 276 kb)
SFN and its metabolites in human plasma following ingestion of a BSE supplement. The levels of "SFN and its metabolites" in human plasma was measured using Liquid Chromatrography-Mass Spectrometry (LC-MS) methods, as described previously . Calculated pharmacokinetic parameters - Tmax, Cmax, AUC, and half-life (mean) values are shown in 'rows' for each SFN metabolite, measured on Days 1 and 7, shown in respective columns. (DOCX 19 kb)