- Open Access
Differential effects of class I isoform histone deacetylase depletion and enzymatic inhibition by belinostat or valproic acid in HeLa cells
Molecular Cancer volume 7, Article number: 70 (2008)
Histone acetylation is an epigenetic modification involved in the regulation of gene expression, balanced by histone acetyl transferases and histone deacetylase (HDAC) enzymes. HDAC inhibitors (HDACi) induce growth arrest and cell death in transformed cells, and are currently in many clinical cancer trials. The transcriptional response to HDACi is complex, as is the response to HDAC isoform knockdown (KD). Here, we describe for the first time in a human cancer cell line, a transcriptional comparison of treatment by two structurally unrelated HDACi; belinostat and valproic acid with the KD of HDAC1, 2 and 3 isoforms.
HDAC KD showed anti-proliferative effects, although to a lesser extent than HDACi treatment. Moreover, we found a 2-fold increased resistance of HDAC1 knockdown cells to belinostat, suggesting this isoenzyme as a selective target. While both HDACi treatment and individual class I HDAC KD produce significant transcriptional effects, three-times higher for HDACi, the gene-expression profiles of class I HDAC KD compared with that obtained by HDACi treatment exhibited less than 4% of altered genes in common between the two modes of inhibition. Further, cell-specific effects of HDAC KD are evident by comparison with a recent study in a different cell line.
The increased resistance to belinostat in response to HDAC1 depletion indicates the possibility of using this isoform as a predictive biomarker of response to HDACi treatment. Further, the transcriptional response to chemical inhibition of HDACs is very different from that of KD of individual class I HDAC isoforms. These data suggest that the anti-tumor effect of HDACi is indeed linked to class I inhibition, but may be more complex than simply targeting individual HDAC enzymes.
The transcription of genes is highly regulated by epigenetic chromatin modifications, including the acetylation of lysine residues protruding from nucleosomal histones. Thus, histone acetylation status is maintained by the opposing actions of histone acetyl transferase and histone deacetylase (HDAC) enzymes [1, 2]. HDACs modify gene expression via multiple mechanisms. The deacetylation of histones causes general chromosome condensation, and also plays a role in transcriptional regulation by forming a combinatorial 'histone code' that regulates downstream responses [2, 3]. Additionally, a variety of non-histone targets such as transcription factors, structural and chaperone proteins are targeted by HDAC enzymes . The Zn2+-dependent mammalian HDAC isoenzymes are divided into three classes based on their homology to yeast deacetylase proteins. Class I HDAC isoforms include HDAC1, -2 and -3 that are ubiquitously expressed as well as the low-abundance HDAC8. Class II (HDAC4, 5, 6, 7, 9, 10) and IV (HDAC11) isoforms display a more restricted tissue pattern of expression . A number of cofactors are required for HDAC activity; indeed, they reside in multi-protein complexes including co-regulators and other chromatin-modifying enzymes .
Recent advances into the biology of HDAC enzymes reveal a substantial division of labor between HDAC subtypes [2, 5]. Modulating HDAC expression demonstrates that class I HDACs are essential for proliferation and survival. Hence, HDAC1 and HDAC3 are believed to be important for proliferation [6–9], whereas HDAC2 is likely involved in the regulation of apoptosis [10, 11]. HDAC8 has been implicated in smooth muscle cell contractility , though its knockdown (KD) also affects proliferation in tumor cells . Class II HDACs are mainly involved in cell differentiation and development , while selective HDAC6 inhibition by tubacin also induced cytotoxicity without accompanying gene-expression changes . Aberrant expression of HDAC1, 2, 3 and 6 has been observed in various tumor types [16–21], and HDAC2-mutant mice display reduced tumor development . Further, the transformed epigenome of neoplastic cells includes specific hypo-acetylation of histone H4 . Together, these findings provide the rationale for the targeted inhibition of HDAC enzymes. HDACi treatment increases global acetylation levels, which ultimately results in cell cycle arrest, apoptosis or terminal differentiation of transformed cells. A considerable variation in the gene-expression response to HDACi depending on cell line and structural class of drug has been demonstrated, and because HDACi treatment potentially affects the entire transcriptome, it is interesting that pan-HDAC inhibition changes the expression of a relatively small percentage of genes [24, 25]. There are several structurally distinct HDACi currently in clinical trials for the treatment of solid and hematological cancers, of which the hydroxamate Zolinza (vorinostat, SAHA), recently gained approval for the treatment of cutaneous T-cell lymphoma .
Despite several reports on the effects of HDAC KD in human and other species, a direct comparison of global gene-expression changes between individual class I HDAC KD and HDACi treatment has not previously been performed on human cancer cell lines. In this report, we examined viability parameters and transcriptional profiles of human HDAC1, 2 and 3 KD, and directly compared expression profiles with treatment of near-IC50 doses of two structurally distinct HDACi; the pan-inhibitory hydroxamate belinostat (PXD101) and the class I selective short-chain fatty acid valproic acid (VPA) . Further, we compared HeLa class I HDAC KD microarray data with that obtained in a recent similar study on U2OS cells.
Depletion of HDAC1, 2 and 3 affect viability
Efficient and specific down-regulation of HDAC1, -2 and -3 was obtained in HeLa cells at both protein (Fig. 1A) and mRNA levels (Table 1, bold), by using the siRNA technology. Viability, as measured by metabolically active cells present in culture, was consistently reduced by 20, 23 and 16% following HDAC1, -2 and -3 KD, respectively (Fig. 1B). A similar effect was seen in HCT116 and MCF-7 cells (data not shown). In HDAC1+2 double KD cells, proliferation was reduced by 35% and 25% when compared with single HDAC1 KD and HDAC2 KD cells, respectively (see additional file 4). Apoptotic effector caspase-3/7 activity was significantly increased for HDAC1, -2 and combination KD (1.65-, 1.45- and 1.63-fold respectively, p < 0.01), but not for HDAC3 KD alone. Further, a dose-response of 1.4, 1.8 and 2.3-fold increased apoptosis at 0.1, 1.0 and 10.0 μM at 24 hours is evident for belinostat treatment (Fig. 1C). No indication of cell cycle deregulation was observed for class I HDAC KD in HeLa at 48 hours post-transfection. However, an increase in the subdiploid population corresponding to fragmented cells was observed for especially HDAC2 and to some extent in HDAC3 KD cells, though not for HDAC1 KD cells. In comparison, belinostat treatment showed marked cell cycle alterations and cell debris (Fig. 1D).
HDAC1 knockdown reduces sensitivity to the HDACi belinostat
Next, we examined how HeLa cells respond to HDACi treatment following individual class I HDAC enzyme down-regulation (Fig. 2). Interestingly, HDAC1 KD significantly increased IC50 values almost 2-fold towards the hydroxamate belinostat (p < 0.01), which was not seen in response to either HDAC2 or 3 depletion. When examining VPA, no significance was observed for either HDAC KD condition.
Gene-expression profiles of belinostat and VPA treatment
Global gene-expression analysis has previously been performed following HDACi treatment regimens primarily in human cell lines [24, 25], but only once recently for individual human class I HDAC KD . However, a direct comparison of gene-expression profiles for each has not been reported. To determine first the transcriptional responses to the HDACi used, DNA chip analyses were performed in independent triplicates for each condition; mock treated control, belinostat and VPA treatment in HeLa cells (0.5 μM and 3.0 mM respectively for 24 hours). The doses chosen were close to the IC50 values in HeLa: 0.76 μM and 3.3 mM for belinostat and VPA respectively, and induced histone H3 and H4 hyper-acetylation (data not shown). Differential gene-expression patterns were detected between each experimental condition versus control. Fig. 3 summarizes the number of non-redundant genes significantly deregulated in response to drug treatment, at an arbitrary 2.0 fold-change cut-off value. The number of genes deregulated by belinostat or VPA is 5.3 and 6.0% respectively. Further, a greater proportion of genes are induced by both drugs (63 and 66% respectively). The relationship of differentially expressed genes between belinostat and VPA was illustrated by Venn diagrams (Fig. 4), and demonstrated that approximately 30% of altered genes responded identically between drug treatments.
In the literature, certain genes have been identified whose expression is affected by HDACi treatment [24, 25]. For instance, Glaser et al. recognized a common set of 13 core genes whose expression were universally altered in response to various HDACi in multiple cell types . Some of the commonly affected genes are listed in Table 1 for belinostat and VPA treatment, and these include up-regulation of CDKN1A (p21) and FUCA and down-regulation of TYMS, CTPS and KPNB1 (importin β) by both drugs. Hence, this study confirms a subset of 5 of the 13 core genes besides other known HDACi target genes. Comprehensive gene-lists of all conditions at above 2-fold changes are accessible (see additional file 2). Data were validated by qRT-PCR analysis on 9 selected genes on RNA samples used in microarray analysis plus independent ones, and an overall good correlation to the microarray data was observed (see additional file 1).
Gene-expression profiles of class I HDAC depletion
Further, we conducted a genome-wide analysis of the transcriptional response to siRNA-mediated depletion of three class I HDACs in HeLa cells. As for HDACi samples, DNA chip analyses were carried out in independent triplicates for each condition; scrambled siRNA control, HDAC1, -2 and -3 KD at 48 hours post-transfection. Differential gene-expression patterns between each knockdown condition and scrambled control were identified by statistical analysis. Efficient KD was confirmed by the microarray data, as each HDAC isoform was specifically down-regulated 7–10 fold (Table 1, bold). The proportion of non-redundant significant transcripts affected by the down-regulation of each HDAC enzyme at 2.0-fold change or more is in the range 1.4–2.0% (shown in Fig. 3), and highest for HDAC1 KD. As for HDACi drugs, a slight overweight of transcripts was induced for HDAC1 and -2 KD samples. In contrast, HDAC3 KD was the only condition not showing this pattern. The proportion of genes with identical expression between KD conditions was in the order of 19–27% (Fig. 4), with HDAC1 KD displaying the least overlap with the other two KD conditions. Further, the individual HDAC isoenzyme targeted by its specific siRNA was the most down-regulated gene in each respective sample, and siRNA targeting one HDAC did not affect expression levels of other class I HDACs. The complete gene lists for all conditions ≥ 2.0-fold changes are accessible (see additional file 2). Also, data were validated by qRT-PCR analysis on 6 selected genes on RNA samples used in microarray analysis plus independent ones, which overall had good correlations to the microarray data (see additional file 1). We further addressed the effect of a combined HDAC1+2 KD by analyzing mRNA expression of three genes (HRASLS3, CCND1 and THBS) found to be affected by each individual HDAC KD. For the CCND1 and THBS genes, KD of either HDAC1 or -2 reduced expression by approximately 50–75% compared to control; an effect not observed in double HDAC1+2 KD (Supplemental figure 2B). For the expression of the HRASLS3 gene, an increase of approximately 50% is seen with single HDAC1 or -2 KD, which increases to approximately 200% in HDAC1+2 double KD cells. Together, these data indicate a degree of redundancy of the HDAC1 and -2 proteins.
Cell-specific effects of individual class I HDAC depletion
A previous report by Senese et al.  studied the transcriptional effect of HDAC1, 2 and 3 KD in the human U2OS osteosarcoma cell line, by microarray analysis. In a direct comparison study, we find very little overlap between the results obtained in the present study and the data recently published by Senese et al. . As discussed below, this apparent discrepancy can be attributed to both methodological and biological differences between the two studies. First, while the experimental design in the Senese study relies on 2 technical replicates of a biological pool on each array, which is then scanned twice (before and after amplification), in our study, we have chosen the more traditional approach with 3 independent biological replicates for each experimental condition and array. Because the biological variation between individual samples most often is much larger than the assay variation , it appears more efficient to perform single arrays on independent biological samples, rather than replicate arrays on a limited number of samples.
Second, we had the opportunity to reanalyze the data (CEL files kindly provided by Dr. Susanna Chiocca, European Institute of Oncology, Milan, Italy) from the Senese study. In our hands, the number of significantly regulated genes following HDAC siRNA inhibition is much lower than reported in the original study (Additional file 5, slide 19). This disparity is probably due to the stringent filtering criteria applied in our analysis, where we require the absolute difference between differentially expressed genes to be 50 or more, as otherwise, the risk of obtaining false positive results due to genes that are close to the background level, is high. Finally, we have made a direct comparison between our knockdown experiments and those performed by Senese et al. From this analysis, it is evident that the overriding differences observed between the two studies are due to cell type specific effects. Both the PCA plot based on 1965 probes that vary between experimental conditions (Additional file5 , slide 29) and the unsupervised hierarchical clustering (Additional file 5, slide 30) imply that the two data sets are clearly separated, and only have a few genes in common when looking at each specific HDAC. That it is indeed the biology that is the main variable between the two studies, and not a technical artifact, is supported by looking at the expression levels of the three HDACs, which are all among the most down-regulated genes, when their corresponding siRNA is transfected into the cells. Thus, we do see HDAC specific effects common for the two studies, but the predominant effects are due to the diverse biology of the two systems. Interestingly, HDAC3 is only number 122 among the most down-regulated genes in the HDAC3 siRNA experiment in the Senese et al. paper. This can be explained by the already low HDAC3 expression in the control cells, so going from a low basal expression to a lower siRNA inhibited expression level will not give as large a fold-change as when looking at a gene with a high basal expression. In summary, the low degree of concordance between the two studies can be attributed mainly to cell specific differences, which further emphasizes the importance and impact of the model system chosen for a particular study (see additonal file 6).
Dissimilar transcriptional profiles between individual HDAC depletion and HDACi treatment
We thus wanted to directly compare individual class I HDAC enzyme knockdown with the treatment of two structurally distinct classes of HDACi in the same cell line, as this has not been done previously in human cells. Approximately three times more genes are deregulated by enzymatic HDACi treatment than by individual class I HDAC depletion (Fig. 3). Between treatment types, a slightly greater proportion of genes are in common for the two structurally distinct HDACi, than are in common for KD of individual HDACs from the same class. However, the biggest difference occurs when comparing single class I HDAC KD with HDACi treatment, as merely 1.6–3.4% of altered genes overlap between these conditions, depending on enzyme and drug assayed (Fig. 4). Thus, only a small proportion of genes respond identically when either knocking down single class I HDAC enzymes or inhibiting multiple HDACs enzymatically. A complete list of genes whose expression is altered in the same direction between one or more HDAC KDs and one or both HDACi treatments are listed (see additional file 3), including genes with ≥ 2.0-fold changes only. This comprises merely 109 genes altered in response to both types of HDAC inhibition, and includes genes involved in processes such as transcriptional regulation, gametogenesis and development, metabolism and intracellular trafficking. Of the 109 genes, 7 are up-regulated simultaneously in all 5 conditions; 3 with unknown function besides HRASLS3 involved in negative cell cycle regulation, CGA involved in cell-cell signaling, PDK4 involved in carbohydrate metabolism and the ion transporter ATP10D.
Targeting cancer through epigenetic control mechanisms is an area of growing interest. While HDACi show promise in clinical trials, the contribution of each HDAC isoenzyme in the anti-proliferative response of HDACi is unknown. In the present study, we directly compared gene-expression profiles between the two modes of HDAC inhibition; single class I HDAC protein depletion by siRNA and enzymatic HDACi treatment in a human cancer cell line. It is recognized, that HDACs function in multi-protein complexes and their depletion therefore might have a dissimilar outcome to HDACi treatment , however this has not been directly addressed previously.
The reduced viability that we observe upon individual HDAC1, -2 and -3 knockdown has been published on class I HDAC KD in cancer cells, especially via proliferation for HDAC1 and -3, and via apoptosis for HDAC2 [6, 8, 11, 19]. We also detected an increased subdiploid population of HDAC2 and less for -3 KD cells, whereas caspase activity was increased for HDAC1 and -2 KD cells. Thus, mediators of apoptosis following HDAC KD might be dissimilar between the isoforms examined. Caspase 3 as a mediator of apoptosis in HDAC1 KD cells was recently reported , as was an increased subdiploid population for HDAC2 KD [10, 19] and HDAC3 KD [6, 19, 28], but not in HDAC1 KD , thus supporting our results. Further, we found no major alterations in cell cycle distribution in response to class I HDAC KD, which is in agreement with other reports [10, 11, 28]. To conclude, class I HDAC KD causes a reduction in viability and an increase in apoptosis, however at much lower levels than detected for HDACi treatment, as this is not transferred to alterations in cell cycle distributions.
Published data suggest a wide range in the proportion of genes deregulated in response to HDACi treatment; between 1–22%. This depends on factors such as class of compound, dosage, incubation time and choice of cell line [24, 25, 29–31]. Hence, our data on belinostat and VPA in HeLa cells are within this broad range. Between belinostat and VPA, the shared proportion of genes of 30% probably correspond to the overlapping functions as HDAC inhibitors as both drugs affect some typical HDACi-induced genes, whereas differences are attributed to structural dissimilarities, HDAC class specificity, and non-HDACi functions of VPA. Other reports comparing the transcriptional response of different HDACi compounds find approximately 45% similarities between trichostatin A (TSA) and either tributyrate or vorinostat and 77% identical genes between tributyrate and vorinostat treatment, when examining three cancer cell lines , while vorinostat and depsipeptide had very similar responses in one cell line, especially in the first hours of treatment . Further, of the limited 'core' set of 13 genes universally affected by HDACi treatment , 5 were reproduced by both drugs in this study. In response to single class I HDAC down-regulation, none of these 13 genes were altered, however the expression of a considerable amount of genes were altered that included genes involved in proliferation, apoptosis or adhesion. For HDAC1, this corresponds to data on C. elegans in which 2.2% were altered by ≥ 1.8-fold , albeit lower than the 7% observed in HDAC1 knockout of untransformed murine embryonic stem cells at 2-fold or more , probably due to the complete abrogation of HDAC1 in this system. HDACi treatment and individual HDAC KD have been shown to cause both up- and down-regulation of multiple gene targets [5, 32, 33]. The knockdown of class I HDAC enzymes in this report showed that near equal proportions of genes were induced as were repressed by HDAC KD, with a slight overweight of induced genes for HDAC1 and -2 KD and a slight overweight of down-regulated genes for HDAC3, possibly separating this isoform as mainly a transcriptional activator. As HDAC1 and -2 reside in the same co-repressor complexes, the disruption of these might have more similar outcomes. Moreover, we found that HDAC1 KD altered the greatest number of genes, and hence might affect gene transcription to a larger extent than HDAC2 and -3. Between the three KD conditions, we found most genes (73–80%) to be uniquely deregulated upon individual HDAC KD, with HDAC1 having the least degree of overlap. This suggests distinctive transcriptional targets for HDAC enzymes from the same class, and could thus provide the basis for discrete functions between class I HDACs . In comparison with genes affected by HDAC1, -2 or -3 KD by siRNA in human U2OS cells in a recent study , the majority were not reproduced herein, and generally point to cell-line specific responses to HDAC depletion. This emphasizes the importance of comparing HDAC KD with HDACi treatment in the same cell line.
Finally, we compared individual KD of class I HDAC members with two dissimilar HDACi compounds at near-IC50 doses. At the treatment regimens chosen, three times more genes were deregulated by HDACi treatment than by individual class I HDAC KD. As these drugs target multiple HDACs, this is not unexpected. The overlap of genes between HDACi treatments and between individual HDAC KD was in a similar range; 20–30%. When looking into the genes whose expression overlapped between HDACi treatment and individual KD of the target HDACs of these compounds, a surprisingly low degree of similarity was observed, namely less than 4% of regulated genes. The reason for the low degree of overlap could have several explanations. First, some degree of redundancy might occur after individual HDAC KD. A prior study in Drosophila showed an overlapping proportion of 20% (469 of 2347 genes regulated in total) between DHDAC1 KD and TSA treatment, each for 5 days post-treatment. However, reducing TSA treatment to 6 hours also reduced the overlap to 4.5% (52/1151), thus differences in experimental set-up probably account for a large variation in these numbers. For DHDAC3 KD, the overlap with TSA treatment was 2%, and the authors conclude that especially DHDAC1 affected gene expression in a similar manner to TSA . The closer resemblance between DHDAC1 and TSA profiles might be because Drosophila has fewer HDAC enzymes and DHDAC1 is orthologous to both human HDAC1 and -2. Second, depleting HDAC levels most likely interferes with the multi-protein complexes in which they reside in a different manner than by enzymatic drug inhibition of HDAC, causing differential cellular responses. It has previously been shown in Drosophila, that DHDAC1 deficiency and point mutations had dissimilar phenotypic outcomes, the latter presumably by altering HDAC complexes rather than disrupting them . Third, we showed that the transcriptional profile obtained by individual HDAC KD is not simply elaborated by inhibiting multiple HDAC enzymes but altered altogether, and thus other mechanisms might contribute to the HDACi effects other than targeting individual class I HDAC enzymes. These differences might explain why single class I HDAC KD is not as toxic as pan-inhibitory HDACi treatment and fails to produce identical phenotypic effects, despite the probable effects of HDACi mainly via class I HDAC enzymes.
Cell culture and drugs
Human cervix cancer cells HeLa, CCL-2 (American Type Culture Collection, Manassas, VA) and mammary cancer cells MCF-7 (HTB-22) were propagated in DMEM + glutamax media supplemented with penicillin and streptomycin and 10% FBS; the colon cancer cell line HCT116 (CCL-247) was maintained in RPMI-1640 media supplemented with glutamine, penicillin, streptomycin and 10% FBS (Invitrogen, Carlsbad, CA). All were grown in a humidified atmosphere of 5% CO2 at 37°C and passaged twice a week. Belinostat was synthesized as described in recent patent applications (International publication number US 6,888,027), and valproic acid was purchased from Sigma-Aldrich (St. Louis, MO). Drugs were dissolved in sterile water, aliquoted and stored at -20°C until use.
Transfection of siRNA
Pre-designed targeting siRNA SmartPOOL was purchased from Dharmacon (Lafayette, CO) (non-targeting siRNA D-001206-13, HDAC1 M-003493-02, HDAC2 M-003495-01, HDAC3 M-003496-00). Cells were plated in 6-well plates, 250,000/well in complete media and incubated overnight prior to aspiration of media and replacement with OPTI-MEM (Invitrogen) with a final concentration of 50 nM siRNA complexed with oligofectamine (Invitrogen). Cells were incubated 4–6 hours before addition of 1 ml growth medium with 20% FCS.
Cells were plated in 6-well plates at 250,000/well and left overnight before the transfection procedure, or replaced with fresh media with drug at 0.5 μM for belinostat or 3.0 mM for VPA. 48 hours post-transfection and 24 hours after drug treatment, total RNA was extracted with Trizol according to the manufacturers protocol (Invitrogen). For microarray samples, 5 of the 6 wells pr condition were pooled to minimize well-to-well variation. The 6th well was lysed directly in SDS sample buffer (Invitrogen), and used for protein analysis.
DNA microarray analysis
RNA integrity was quality checked on the Agilent 2100 Bioanalyser (Agilent Technologies, Santa Clara, CA), then processed and hybridized onto Affymetrix arrays according to the manufacturer's protocol. Briefly, 5 μg RNA pr sample was used to to generate biotin-labeled antisense cRNA. After fragmentation, the labeled cRNA samples were hybridized to Affymetrix HG-U133 Plus 2.0 arrays (Affymetrix, Santa Clara, CA), washed and stained with phycoerytrin conjugated streptavidin, and finally scanned in the Affymetrix GeneArray® scanner to generate fluorescent images, as described in the Affymetrix GeneChip® protocol. All statistical analyses, including pre-processing of data was carried out in R, version 2.3 (R development core team. R: A language and Environment for Statistical Computing; http://www.R-project.org). DNA chips were checked for quality assurance parameters such as visual image inspection, replicate scatter plots and RNA degradation plots, before normalization for mean overall expression using the gcrma package (GC robust multiarray algorithm) . Agglomerative hierachical clustering showed that the biological replicates clustered together as expected (data not shown). The statistical linear model based method of Limma (linear models for microarray data package)  was found to be most sensitive at identifying genes with differential expression between control and each condition. Raw p-values were adjusted for multiple testing using the Benjamini & Hochberg method  to reduce the number of false positives, and a 5% significance threshold applied. Comparisons of gene lists between conditions was performed using VennMapper .
Quantitative reverse transcription polymerase chain reaction
RNA was reverse transcribed by RT-PCR using 1 volume diluted RNA (100 ng/μl) and 1 volume 2× RT-master mix (High Capacity cDNA archive kit, Applied Biosystems, Foster City, CA) exposed to 25°C 10 minutes and 37°C 2 hours. Samples were analyzed for gene-expression levels by qRT-PCR, performed on ABI PRISM™ 7500 Sequence Detection System (Applied Biosystems). Experiments were done in triplicate by mixing 1 μL probe, 10 μL 2× Taqman master mix and 9 μL cDNA diluted 1:50, and subjecting samples to 40 cycles of amplification (15 seconds denaturation at 95°C, 1 minute annealing and elongation at 60°C), using GAPDH (HDAC KD) or β-Actin (drug treatments) as endogenous controls. All pre-validated FAM-labeled probes were purchased from Applied Biosystems. Subsequent data analysis was performed using DART-PCR version 1.0 .
Cells were lysed directly in SDS sample buffer (Invitrogen) and electrophoretically separated and transferred to nitrocellulose paper, following the manufacturer's instructions (Invitrogen pre-cast Novex gels). Blots were blocked in 10% non-fat skimmed milk, incubated overnight with primary antibody, washed in Tris-buffered saline with Tween-20 and visualized by HRP-conjugated secondary antibodies (1:2,000) and ECL plus reagent (GE Healthcare, Chalfont St. Giles, Buckinghamshire, UK). Antibodies uses were: rabbit anti-HDAC1 and -3 (1:1,000, Cell Signaling Technology, Danvers, MA), mouse anti-HDAC2 (1:10,000, Abcam, Cambridge, UK), rabbit anti-Actin (1:5,000 Sigma-Aldrich).
Scrambled control and HDAC KD cells were plated in triplicate at 10,000/well in 96-well format 24 hours post-transfection. Cells were incubated 48 hours, without drug for viability measurements, and within expected belinostat toxicity limits for determination of IC50 values. Cells were lysed directly with CellTiter-Glo luminescent viability assay (Promega, Madison, WI), and luminescence proportional to ATP present hence metabolically active cells was measured (HTS 7000 plus Bioassay reader, PerkinElmer Life and Analytical Sciences, Waltham, MA). Data were normalized to the scrambled control, and IC50 values determined in Prism 4 by generation of a sigmoidal dose-response curve with variable slope (GraphPad Software, San Diego, CA). Significant changes in mean viability or IC50 values in the four groups (control, HDAC1, -2 or -3 KD) were calculated by ANOVA one-way analysis of variance repeated measures test and Dunnett's multiple comparisons test in Prism (GraphPad Software).
Caspase-Glo 3/7 assay
Cells were plated at 104/well, in quadroplicates for each HDAC KD condition and control 48 hours post-transfection, or in triplicates for drug treatments. Plates were incubated for 24 hours prior to direct lysis by Caspase-Glo 3/7 reagent (Promega), and luminescence reading according to caspase 3/7 activity.
Cell cycle analysis
Transfected cells were incubated for 48 hours before analysis. For drug treatments, cells were plated in 6-well format 2.5 × 105/well, incubated overnight, treated with drug for 24 hours and processed as follows. Cells were permeabilized by incubation in ice-cold 70% ethanol, rehydrated in PBS supplemented with Tween-20 and FBS, RNAse treated and DNA stained with propidium iodide. Cells were analyzed using a FACSCalibur instrument and the CellQuest software (BD Biosciences, Mountain View, CA).
histone deacetylase inhibitor
half maximal inhibitory concentration
quantitative reverse transcriptase polymerase chain reaction
small interfering ribonucleic acid
Bolden JE, Peart MJ, Johnstone RW: Anticancer activities of histone deacetylase inhibitors. Nat Rev Drug Discov. 2006, 5: 769-784. 10.1038/nrd2133
Gallinari P, DiMarco S, Jones P, Pallaoro M, Steinkuhler C: HDACs, histone deacetylation and gene transcription: from molecular biology to cancer therapeutics. Cell Res. 2007, 17: 195-211.
Strahl BD, Allis CD: The language of covalent histone modifications. Nature. 2000, 403: 41-45. 10.1038/47412
Minucci S, Pelicci PG: Histone deacetylase inhibitors and the promise of epigenetic (and more) treatments for cancer. Nat Rev Cancer. 2006, 6: 38-51. 10.1038/nrc1779
Bernstein BE, Tong JK, Schreiber SL: Genomewide studies of histone deacetylase function in yeast. Proc Natl Acad Sci USA. 2000, 97: 13708-13713. 10.1073/pnas.250477697
Glaser KB, Li J, Staver MJ, Wei RQ, Albert DH, Davidsen SK: Role of class I and class II histone deacetylases in carcinoma cells using siRNA. Biochem Biophys Res Commun. 2003, 310: 529-536. 10.1016/j.bbrc.2003.09.043
Khabele D, Son DS, Parl AK, Goldberg GL, Augenlicht LH, Mariadason JM, Rice VM: Drug-Induced Inactivation or Gene Silencing of Class I Histone Deacetylases Suppresses Ovarian Cancer Cell Growth: Implications for Therapy. Cancer Biol Ther. 2007, 6: 795-801.
Lagger G, O'Carroll D, Rembold M, Khier H, Tischler J, Weitzer G, Schuettengruber B, Hauser C, Brunmeir R, Jenuwein T, Seiser C: Essential function of histone deacetylase 1 in proliferation control and CDK inhibitor repression. EMBO J. 2002, 21: 2672-2681. 10.1093/emboj/21.11.2672
Senese S, Zaragoza K, Minardi S, Muradore I, Ronzoni S, Passafaro A, Bernard L, Draetta GF, Alcalay M, Seiser C, Chiocca S: A role for histone deacetylase 1 in human tumor cell proliferation. Mol Cell Biol. 2007, 27: 4784-4795. 10.1128/MCB.00494-07
Huang BH, Laban M, Leung CH, Lee L, Lee CK, Salto-Tellez M, Raju GC, Hooi SC: Inhibition of histone deacetylase 2 increases apoptosis and p21Cip1/WAF1 expression, independent of histone deacetylase 1. Cell Death Differ. 2005, 12: 395-404. 10.1038/sj.cdd.4401567
Zhu P, Huber E, Kiefer F, Gottlicher M: Specific and redundant functions of histone deacetylases in regulation of cell cycle and apoptosis. Cell Cycle. 2004, 3: 1240-1242.
Waltregny D, Glenisson W, Tran SL, North BJ, Verdin E, Colige A, Castronovo V: Histone deacetylase HDAC8 associates with smooth muscle alpha-actin and is essential for smooth muscle cell contractility. FASEB J. 2005, 19: 966-968.
Vannini A, Volpari C, Filocamo G, Casavola EC, Brunetti M, Renzoni D, Chakravarty P, Paolini C, De FR, Gallinari P, Steinkuhler C, DiMarco S: Crystal structure of a eukaryotic zinc-dependent histone deacetylase, human HDAC8, complexed with a hydroxamic acid inhibitor. Proc Natl Acad Sci USA. 2004, 101: 15064-15069. 10.1073/pnas.0404603101
Verdin E, Dequiedt F, Kasler HG: Class II histone deacetylases: versatile regulators. Trends Genet. 2003, 19: 286-293. 10.1016/S0168-9525(03)00073-8
Haggarty SJ, Koeller KM, Wong JC, Grozinger CM, Schreiber SL: Domain-selective small-molecule inhibitor of histone deacetylase 6 (HDAC6)-mediated tubulin deacetylation. Proc Natl Acad Sci USA. 2003, 100: 4389-4394. 10.1073/pnas.0430973100
Choi JH, Kwon HJ, Yoon BI, Kim JH, Han SU, Joo HJ, Kim DY: Expression profile of histone deacetylase 1 in gastric cancer tissues. Jpn J Cancer Res. 2001, 92: 1300-1304.
Sakuma T, Uzawa K, Onda T, Shiiba M, Yokoe H, Shibahara T, Tanzawa H: Aberrant expression of histone deacetylase 6 in oral squamous cell carcinoma. Int J Oncol. 2006, 29: 117-124.
Sasaki H, Moriyama S, Nakashima Y, Kobayashi Y, Kiriyama M, Fukai I, Yamakawa Y, Fujii Y: Histone deacetylase 1 mRNA expression in lung cancer. Lung Cancer. 2004, 46: 171-178. 10.1016/j.lungcan.2004.03.021
Wilson AJ, Byun DS, Popova N, Murray LB, L'Italien K, Sowa Y, Arango D, Velcich A, Augenlicht LH, Mariadason JM: 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-13558. 10.1074/jbc.M510023200
Zhang Z, Yamashita H, Toyama T, Sugiura H, Omoto Y, Ando Y, Mita K, Hamaguchi M, Hayashi S, Iwase H: HDAC6 expression is correlated with better survival in breast cancer. Clin Cancer Res. 2004, 10: 6962-6968. 10.1158/1078-0432.CCR-04-0455
Zhu P, Martin E, Mengwasser J, Schlag P, Janssen KP, Gottlicher M: Induction of HDAC2 expression upon loss of APC in colorectal tumorigenesis. Cancer Cell. 2004, 5: 455-463. 10.1016/S1535-6108(04)00114-X
Zimmermann S, Kiefer F, Prudenziati M, Spiller C, Hansen J, Floss T, Wurst W, Minucci S, Gottlicher M: Reduced body size and decreased intestinal tumor rates in HDAC2-mutant mice. Cancer Res. 2007, 67: 9047-9054. 10.1158/0008-5472.CAN-07-0312
Fraga MF, Ballestar E, Villar-Garea A, Boix-Chornet M, Espada J, Schotta G, Bonaldi T, Haydon C, Ropero S, Petrie K, Iyer NG, Perez-Rosado A, Calvo E, Lopez JA, Cano A, Calasanz MJ, Colomer D, Piris MA, Ahn N, Imhof A, Caldas C, Jenuwein T, Esteller M: Loss of acetylation at Lys16 and trimethylation at Lys20 of histone H4 is a common hallmark of human cancer. Nat Genet. 2005, 37: 391-400. 10.1038/ng1531
Glaser KB, Staver MJ, Waring JF, Stender J, Ulrich RG, Davidsen SK: Gene expression profiling of multiple histone deacetylase (HDAC) inhibitors: defining a common gene set produced by HDAC inhibition in T24 and MDA carcinoma cell lines. Mol Cancer Ther. 2003, 2: 151-163.
Gray SG, Qian CN, Furge K, Guo X, Teh BT: Microarray profiling of the effects of histone deacetylase inhibitors on gene expression in cancer cell lines. Int J Oncol. 2004, 24: 773-795.
Khan N, Jeffers M, Kumar S, Hackett C, Boldog F, Khramtsov N, Qian X, Mills E, Berghs SC, Carey N, Finn PW, Collins LS, Tumber A, Ritchie JW, Jensen PB, Lichenstein HS, Sehested M: Determination of the class and isoform selectivity of small molecule HDAC inhibitors. Biochem J. 2007, 409: 581-589. 10.1042/BJ20070779
Glaser KB: HDAC inhibitors: Clinical update and mechanism-based potential. Biochem Pharmacol. 2007, 74: 659-671. 10.1016/j.bcp.2007.04.007
Takami Y, Nakayama T: N-terminal region, C-terminal region, nuclear export signal, and deacetylation activity of histone deacetylase-3 are essential for the viability of the DT40 chicken B cell line. J Biol Chem. 2000, 275: 16191-16201. 10.1074/jbc.M908066199
de Ruijter AJ, Meinsma RJ, Bosma P, Kemp S, Caron HN, van Kuilenburg AB: Gene expression profiling in response to the histone deacetylase inhibitor BL1521 in neuroblastoma. Exp Cell Res. 2005, 309: 451-467. 10.1016/j.yexcr.2005.06.024
Moore PS, Barbi S, Donadelli M, Costanzo C, Bassi C, Palmieri M, Scarpa A: Gene expression profiling after treatment with the histone deacetylase inhibitor trichostatin A reveals altered expression of both pro- and anti-apoptotic genes in pancreatic adenocarcinoma cells. Biochim Biophys Acta. 2004, 1693: 167-176. 10.1016/j.bbamcr.2004.07.001
Peart MJ, Smyth GK, van Laar RK, Bowtell DD, Richon VM, Marks PA, Holloway AJ, Johnstone RW: Identification and functional significance of genes regulated by structurally different histone deacetylase inhibitors. Proc Natl Acad Sci USA. 2005, 102: 3697-3702. 10.1073/pnas.0500369102
Whetstine JR, Ceron J, Ladd B, Dufourcq P, Reinke V, Shi Y: Regulation of tissue-specific and extracellular matrix-related genes by a class I histone deacetylase. Mol Cell. 2005, 18: 483-490. 10.1016/j.molcel.2005.04.006
Zupkovitz G, Tischler J, Posch M, Sadzak I, Ramsauer K, Egger G, Grausenburger R, Schweifer N, Chiocca S, Decker T, Seiser C: Negative and positive regulation of gene expression by mouse histone deacetylase 1. Mol Cell Biol. 2006, 26: 7913-7928. 10.1128/MCB.01220-06
Foglietti C, Filocamo G, Cundari E, De RE, Lahm A, Cortese R, Steinkuhler C: Dissecting the biological functions of Drosophila histone deacetylases by RNA interference and transcriptional profiling. J Biol Chem. 2006, 281: 17968-17976. 10.1074/jbc.M511945200
Mottus R, Sobel RE, Grigliatti TA: Mutational analysis of a histone deacetylase in Drosophila melanogaster: missense mutations suppress gene silencing associated with position effect variegation. Genetics. 2000, 154: 657-668.
Wu Z, Irizarry RA, Gentleman R, Murillo FM, Spencer F: A Model Based Background Adjustment for Oligonucleotide Expression Arrays. Technical Report, John Hopkins University, Department of Biostatistics Working Papers. 2003.
Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article 3.
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc B. 1995, 57: 289-300.
Smid M, Dorssers LC, Jenster G: Venn Mapping: clustering of heterologous microarray data based on the number of co-occurring differentially expressed genes. Bioinformatics. 2003, 19: 2065-2071. 10.1093/bioinformatics/btg282
Peirson SN, Butler JN, Foster RG: Experimental validation of novel and conventional approaches to quantitative real-time PCR data analysis. Nucleic Acids Res. 2003, 31: e73- 10.1093/nar/gng073
Simon RM, Korn EL, McShane LM, Radmacher MD, Wright GW, Zhao Y: Design and analysis of DNA microarray investigations. 2003, New York: Springer.
We thank Sanne Jul Christiansen, Mette Frandsen and Annette Nielsen for valuable technical skills, and Daniel Llorente Tabor for helping with IT (all are TopoTarget A/S employees).
Declaration of competing interests
Marielle Dejligbjerg was partly sponsored by a grant from the Danish Ministry of Science, Technology and Innovation, and the biotechnology company TopoTarget A/S. Morten Grauslund is a shareholder and full-time employee at TopoTarget A/S. Thomas Litman is an employee of Exiqon A/S. Laura Collins is a previous employee of TopoTarget A/S. Xiaozhong Qian and Michael Jeffers are previous employees of CuraGen Corp. Henri Lichenstein is a shareholder and full-time employee at CuraGen Corp. Peter B Jensen is CEO of and shareholder in TopoTarget A/S. Maxwell Sehested is CSO of and shareholder in TopoTarget A/S.
MD designed and performed siRNA transfections, qRT-PCR, western blotting, viability assays and flow cytometry analyses, helped with secondary microarray analysis and drafted the manuscript. MG participated in the design of the study and drafting of the manuscript. TL and LC performed microarray data analysis. XQ assisted in siRNA transfections, western blotting and flow cytometry analyses. MJ, HL, PBJ and MS conceived of the study, participated in its design and coordination, and provided intellectual discussions and ideas regarding the content of manuscript. All participating authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Additional file A: Supplemental Figure 1. Validation of gene-expression changes from microarray analysis, by qRT-PCR. For each gene, data from microarray analysis and corresponding values from qRT-PCR analysis are included. Validation was shown for both one of the same RNA purifications as were used for microarray analysis ("A"), and independent ones in HeLa ("B"). A) Genes affected by HDACi treatment, B) Genes affected by HDAC KD (DOC 426 KB)
Additional file 2: Additional file B: Supplemental Table A: Complete gene lists at ≥ 2-fold changes. Complete gene lists of significantly altered genes (5%) in all five conditions: HDAC 1, 2 or 3 KD, VPA or belinostat treatments, at 2-fold changes or more. (XLS 576 KB)
Additional file 3: Additional file C: Supplemental Table B: Transcriptional overlap between conditions. The overlapping proportion of genes regulated by one or more HDACi treatments and by one or more class I HDAC depletion, by at least 2-fold. Thus, these genes respond identically to both types of HDAC inhibition. (XLS 42 KB)
Additional file 4: Additional file D: Supplemental Figure 2. Effect of HDAC1, 2 and 1+2 depletion using the siRNA technology in HeLa cells. A). Effect on cell proliferation by measuring intracellular ATP levels. B) Effect on gene expression by measuring mRNA levels of selected genes by qRT-PCR. (DOC 30 KB)
Additional file 5: Additional file E: Detailed comparison of microarray data from this study to the Senese et al. study . The Power Point file contains 40 slides showing how the analysis was carried out and the results obtained herein (PPT 5 MB)
Additional file 6: Additional file F: Supplemental Table C: Combined analysis of the Senese and Dejligbjerg data. The Excel table allows direct comparison and sorting of the genes that were more than 1.5-fold differentially expressed following HDAC knockdown in the two studies. (XLS 4 MB)
About this article
Cite this article
Dejligbjerg, M., Grauslund, M., Litman, T. et al. Differential effects of class I isoform histone deacetylase depletion and enzymatic inhibition by belinostat or valproic acid in HeLa cells. Mol Cancer 7, 70 (2008). https://doi.org/10.1186/1476-4598-7-70
- Individual Class
- HDACi Treatment
- HDAC Enzyme
- Individual HDAC