Histone deacetylase inhibitors: A new mode for inhibition of cholesterol metabolism

Background Eukaryotic gene expression is a complex process involving multiple cis and trans activating molecules to either facilitate or inhibit transcription. In recent years, many studies have focused on the role of acetylation of histone proteins in modulating transcription, whereas deacetylation of these same proteins is associated with inactivation or repression of gene expression. This study explores gene expression in HepG2 and F9 cell lines treated with Trichostatin A (TSA), a potent histone deacetylase inhibitor. Results These experiments show that TSA treatment results in clear repression of genes involved in the cholesterol biosynthetic pathway as well as other associated pathways including fatty acid biosynthesis and glycolysis. TSA down regulates 9 of 15 genes in this pathway in the F9 embryonal carcinoma model and 11 of 15 pathway genes in the HepG2 cell line. A time course study on the effect of TSA on gene expression of various enzymes and transcription factors involved in these pathways suggests that down regulation of Srebf2 may be the triggering factor for down regulation of the cholesterol biosynthesis pathway. Conclusion Our results provide new insights in the effects of histone deacetylases on genes involved in primary metabolism. This observation suggests that TSA, and other related histone deacetylase inhibitors, may be useful as potential therapeutic entities for the control of cholesterol levels in humans.


Background
Histone deacetylases (HDACs) are important chromatin remodeling enzymes that are generally involved in transcriptional repression [1]. Mammalian HDACs are classified into three main categories depending on their primary homology to Saccharomyces cerevisiae HDACs (RPD3, HDA1 and SIR2). Histone deacetylase inhibitors (HDACIs) tend to show equal effects on gene activation and repression [2][3][4]. HDACIs have been shown to induce differentiation, apoptosis or growth arrest in a variety of transformed cell lines [5]. This is generally attributed to the ability of these inhibitors to induce an open chromatin conformation facilitating transcription of regulatory genes like p21 which inhibit tumor cell growth [6]. These qualities make HDACIs promising targets for chemotherapeutic intervention.
Recently many different types of HDAC inhibitors have been discovered (Figure 1). These include short chain fatty acids (sodium butyrate, phenylbutyrate, valproic acid) [7], hydroxamic acids (trichostatin A (TSA), suberoylanilide hydromaxic acid (SAHA), pyroxamide, cyclic hydroxamic acid-containing peptides (CHAPs), cinnamic acid bishydroxamic acid (CBHA) and scriptaid) [8,9], cyclic tetrapeptides (trapoxin, apicidin, depsipeptide) [10][11][12][13]13], and benzamides (MS-275) [14,15]. Most HDAC inhibitors (HDACIs) developed to date inhibit both Class I and II HDACs equally with the exceptions being valproic acid (5 fold more selective for HDAC1 vs HDACs 5 and 6) and FK-228 (Class I selective). Class I and II HDACs are inhibited by trichostatin A (TSA) and related compounds whereas Class III HDACs are not. As noted, HDACIs have been shown to promote cell cycle arrest, differentiation, and apoptosis in many transformed cultured cell types. In animal models, HDACIs have been shown to inhibit growth of breast, prostate, lung and stomach cancers, as well as neuroblastomas and leukemias, with little toxicity [16,17]. In a previous study looking at the combination regimen of all trans retinoic acid (RA) with the HDACI, Trichostatin A (TSA), we identified several new targets for HDACIs [18]. We also identified critical differences in gene regulation subsequent to treatment with these two agents and a novel promoter module associated with the regulation of a subset of these differentially regulated genes. These analyses focused on the anticancer therapeutic potential of these compounds alone or in combination. Recent analysis of these data identified certain crucial metabolic pathways that have not previously been shown to respond to HDACI treatment and which may be critical in identifying new therapies for cardiovascular health. In this report we discuss the possible role of HDAC inhibition on cholesterol metabolism.

Microarray results from F9 cell treatments
Of the 12,451 mouse genes on the Affymetrix MU74Av2 microarray, 1248 genes (upregulated expression of 489 genes and decreased expression of 759 genes) were found to be significantly differentially expressed following TSA treatment. Of these, only 463 genes were found to be differentially expressed at an arbitrary two-fold or greater level of expression (226 genes up; 237 genes down) (Tables 1 &2, Additional file 1). The raw CEL files for the microarray data are available for download at the Gene Expression Omnibus under series GSE1437. Genes for which up regulated expression was noted were involved in retinoid binding and/or metabolism (e.g., Crabp2, Rbp1,  Mybl2, Tal1). Interestingly, we also found down regulated genes involved in pyrimidine biosynthesis (Dhodh) and in the cholesterol metabolism pathway (Mvk, Lss, Hmgcr, Fasn and Sqle) ( Figure 2). This latter finding was intriguing and we decided to extend our investigations beyond the pluripotent mouse EC cells. These experiments were repeated using the more relevant human hepatocarcinoma derived HepG2 cells, since these are hepatic in origin and the liver is the primary source for cholesterol and fatty acid metabolism in humans. While these are not primary hepatocytes, this cell line offers the ability to both explore the hereto unknown effects of TSA on cholesterol metabolism and also look at the previously known targets of this drug.

Microarray results from HepG2 experiments
Of the 54,613 human genes on the Affymetrix HU133 plus 2.0 array, only 6,513 showed significant differential expression following TSA treatment (p value < 0.05). The raw CEL files for the microarray data are available for download at the Gene Expression Omnibus under series GSE4465. TSA treatment of this cell line resulted in 1561 genes being up regulated and 4952 genes being down regulated at this level of significance ( Figure 3). This observation was surprising since the current paradigm of HDAC inhibition suggests equal effects on gene activation and repression after treatment with an HDAC inhibitor [3,4]. Furthermore, using a two-fold cutoff, the genelist was reduced to 3229 genes with 254 up and 2975 down regulated (see Tables 3 and 4; Additional file 1) further emphasizing the extreme extent of gene expression down regulation following HDACI treatment. Genes that showed up regulated expression include phospholipid transfer protein (Pltp), tissue inhibitors of metalloproteinase 1 and 2 (Timp1, Timp2) and transforming growth factor beta 1 (Tgfβ1). Perhaps more importantly, the interesting downregulated genes included thymidylate synthetase (Tyms), formyltetrahydrofolate dehydrogenase (Fthfd), dihydroorotate dehydrogenase (Dhodh) and CTP synthase II (CTPSII) (all of which are related to pyrimidine biosynthesis) as well as genes related to lipid transport and fatty acid metabolism including low density lipoprotein receptor (ldlr), enoyl-Coenzyme A hydratase, acyl-Coenzyme A dehydrogenase (Acadm), apolipoproteins A5, C3, L1, high density lipoprotein binding protein (vigilin); and 3-hydroxy-3-methylglutaryl-Coenzyme A synthase 1 (Hmgcs1), farnesyl-diphosphate farnesyltransferase 1 (Fdft1), squalene epoxidase (Sqle), sterol regulatory element binding transcription factor 2 (Srebf2) and 7dehydrocholesterol reductase (Dhcr7). Notably 11 of these genes are involved in cholesterol metabolism (Figure 2).

Quantitative PCR results
Sybr green qPCR was used to validate microarray expression data for a subset of the differentially expressed genes. The expression patterns of 10 genes from the F9 microarray data set and 21 from the HepG2 microarray data set ( Figure 4) were all confirmed by qPCR. Furthermore we decided to examine the levels of gene expression at early and late time points for 11 of these genes that have a role in cholesterol and lipid metabolism. The relative gene expression was obtained for these genes at 3 h, 6 h, 9 h, 12

Discussion
In a previous study we had used microarray analyses to examine the effects of RA and TSA on embryonal carcinoma cell growth and differentiation using the prototypical EC cell line F9 [18]. Results from these studies identified several important genes and pathways differentially regulated by these compounds. In this report we identify new target pathways for TSA treatment based on further analysis of this data. Most importantly, the regulatory pathways that are affected include pyrimidine metabolism and cholesterol biosynthesis. The pyrimidine pathway is of interest because one of the rate limiting enzymes in this pathway, dihydroorotate dehydrogenase (dhodh), has been targeted for inhibition in murine models of rheumatoid arthritis as well as in the human T-lymphoblastoma cell line (A3.01) [21,22]. Dhodh catalyzes the fourth committed step in the de novo biosynthesis of pyrimidines. Activated lymphocytes expand their pyrimidine pool by eightfold during proliferation [23]. In rheumatoid arthritis, inflammation and degradation of synovial tissue are initiated by the influx of lymphocytes (B cells, CD4+, CD8+ and T cells) [24]. Thus, inhibiting activated T-cells by decreasing their supply of pyrimidines via TSA treatment could provide an attractive alternative method for treating rheumatoid arthritis. Interestingly, HDACIs (TSA and phenylbutyrate) were used as treatments in a rat model of rheumatoid arthritis, and resulted in reduced inflammation, and inhibition both of synovial hyperplasia and bone or cartilage destruction [25]. The authors also found that HDACIs inhibited the expression of tumor necrosis factor-α, which functions to stimulate matrix degradation in rheumatoid arthritis [26], therefore suggesting a mechanism by which HDACIs may alleviate some effects of rheumatoid arthritis. Further extending and supporting these results, in this study, we found that TSA itself could significantly inhibit the expression of dhodh even in non-lymphatic cells (60% in F9 and 25% in HepG2), providing an alternative (or synergistic) mechanism by which HDACI might suppress rheumatoid arthritis in both mice and men. Moreover, the mRNA levels of thymidylate synthetase (Tyms), another key enzyme in this pathway, were decreased 8 fold in HepG2 cells (2.4 fold in F9 cells) further potentiating this effect of TSA treatment. A previous study using chondrocytes showed that HDACIs such as TSA and sodium butyrate, blocked the induction of matrix metalloproteinases (MMP-1, MMP-13) as well as aggrecan-degrading enzymes (Adamts4, Adamts5 and Adamts9) [27]. Both of these enzyme families mediate cartilage destruction. In our The pathway most significantly affected by TSA treatment in F9 EC cells is that of cholesterol biosynthesis, most specifically those steps involved in the synthesis of low density lipoprotein. There are two main types of lipoproteins that transport cholesterol in the blood: low density lipoproteins (LDL) and high density lipoproteins (HDL). HDL particles are generally considered to be "good cholesterol", while LDL is considered "bad cholesterol" [29]. Several genes encoding essential enzymes in the LDL synthesis pathway are down regulated by these treatments. Pathway analysis of microarray data using genes showing statistically significant (p < 0.05) differential gene expression indicated that expression levels of 9 enzymes out of the 15 in the cholesterol biosynthesis pathway are decreased following TSA treatment. They include HMG CoA reductase (Hmgcr), mevalonate kinase (Mvk), di-pmevalonate decarboxylase (Mvd), isopentenyl-PP isomerase (Idi1), squalene synthatase (Fdft1), squalene epoxidase (Sqle), lanosterol synthase (Lss) and lanosterol oxidase (Sc4mol) and NAD(P)-dependent steroid dehydrogenase (Nsdhl) (Figure 2). Of these only 6 of these 9, including Hmgcr, Mvd, Idi1, Sqle, Lss, Sc4mol and Nsdhl, showed a greater than two fold differential expression on TSA treatment. The z-score assigned to each category by MAPPFinder reflects the degree to which the expression of genes in that category was greater than that expected by chance. A positive z-score indicates that a large number of genes in that category are differentially expressed between the compared conditions while a negative Z score indicates that the there are fewer genes meeting the criterion than would be expected by random chance. If the MAPP-Finder data truly obeyed the assumptions of the hypergeometric distribution, then a Z score or 1.96 or -1.96 would correlate with a p value of 0.05. The z-score for this pathway was highly significant with values of 4.85 associated with TSA treatment. Real-time qPCR analyses verified decreased expression of 3 of 5 genes (Hmgcr, Mvd, Lss) in this pathway but the expression of neither Mvd nor Lss was significantly down regulated by TSA. Decreasing high cholesterol levels using TSA treatments may work well since repression of few of the detected genes may be sufficient to induce the response, i.e., a reduction in cholesterol intermediates and synthesis.
Following the analysis of these pluripotent EC cells, we decided to investigate these effects in the HepG2 cell line which arose from a carcinoma of the human liver, the primary organ for cholesterol and fatty acid metabolic processes. While we realize that primary hepatocytes would be a better model to evaluate this pathway, we choose the HepG2 cell line as an means to evaluate this phenomenon but allowing for use of the known anti-cancer effects of TSA as a control. Expression data from HepG2 cells also indicated that multiple enzymes in cholesterol biosynthesis and fatty acid synthesis pathways were significantly down regulated (Figure 3). The mRNA transcript levels that were repressed at a greater than 2 fold level of significance included HMG CoA synthase (Hmgcs1), HMG CoA reductase (Hmgcr), sterol receptor binding factor-2 (Srebf2) and lanosterol 14 α-demethylase (Cyp51a1) (involved in cholesterol metabolism), and others including fatty acid synthase (Fasn), fatty acid binding protein (Fabp), farnesyl diphosphate synthase (Fdps), acetyl-coA carboxylase (Acaca), acetyl-coA dehydrogenase (Acadm), acetyl-coA acetyl transferase (Acat2), peroxisome proliferative activated receptor, gamma (Pparγ) and a variety of apolipoproteins that are involved in fatty acid and triglyceride metabolism. Quantitative PCR studies verified that TSA treatment reduced expression of Hmgcr, Hmgcs1, Srebf2, Fabp Fasn, Fdps, Acaca, Acadm, Acat2, ApoA5, C1, E and L1 as well as Cyp27a1, Ldlr, Pparγ and Tyms (Figure 4). The down regulation seems to be a complex phenomenon involving genes that regulate these pathways at different levels. Most evident is the down regulation of Srebf2 which in turn acts as a transcription factor regulating the expression of enzymes like Hmgcr (the target for the statin class of drugs) and Mvd. It is known that Srebf2 overexpression induces all 12 enzymes in the cholesterol biosynthesis pathway and inhibition of Srebf2 by TSA might inhibit the expression of these enzymes [30]. In fact, our microarray data demonstrates that the levels of almost all these enzymes are down (10 genes pass all cutoff filters) following TSA treatment. The repression of Srebf2 occurs at a early time point (around 3 hrs) and continually repressed over 24 hrs (2.4 fold repression at 24 h). This effect is probably responsible for the down regulation of the cholesterol pathway since expression of Hmgcs (repressed at 6 h), Hmgcr (repressed at 12 h) and Ldlr (repressed at 24 h) are all known to be induced by Srebf2 [31]. Srebf-1a and -1c are more involved with regulation of fatty acid synthesis and lipogenesis [30,32]. While we were unable to detect the levels of Srebf1 expression in our microarray experiments, TSA treatment modestly (1.3 fold) decreased expression levels of cytosolic NADPdependent isocitrate dehydrogenase1 (Idh1) which pro-vides the cytosolic NADPH required for proper functioning of both the cholesterol and fatty acid biosynthetic pathways. The Idh1 promoter is activated by Srebf1a and Srebf2 in human hepatoma cells [33]. Thus there appears to be a concerted down regulation of both pathways through a synergistic effect. Cyp51a1 (lanosterol 14αdemethylase) is another important intermediate in cholesterol metabolism and has in recent years gained importance as a target for the development of hypocholesterolemic agents [34]. Our microarray data showed that levels of Cyp51a1 were also down regulated further adding support to the possible use of this HDACI as a way to decrease plasma cholesterol levels.
Finally, atherosclerosis is the underlying disorder associated with most cardiovascular disease [35]. This disorder is characterized by deposits of fatty substances, cholesterol, cellular waste products, calcium and other substances in the inner lining of an artery (collectively known as plaques) [36]. Cholesterol has been implicated as the major contributor to this condition as atherosclerosis is strongly correlated with an increase in serum cholesterol levels [37,38]. Generally, serum levels should be between 140 and 200 mg per deciliter (mg/dl) whereas high levels surpassing 240 mg/dl indicate one is at high risk for cardiovascular disease [39]. Thus, atherosclerosis is characterized by elevated levels of LDL [40]. The activity of the hepatic LDL receptor (Ldlr) is the primary determinant of plasma LDL cholesterol levels and Ldlr transcription is in turn regulated by Srebf2. When the levels of hepatocellular sterols drop, Srebf2 is activated and this process restores the normal levels by concurrent activation of de novo cholesterol synthesis and increased uptake of plasma cholesterol through Ldlr. LDL receptor is also post transcriptionally regulated by proprotein convertase subtilisn/kexin type 9a (Pcsk9) in an inverse manner [41]. While our microarray and qPCR data shows decreased Ldlr expression following TSA treatment, microarray gene expression levels of Pcsk9 are also down regulated. This suggests existence of a mechanism for potential compensatory increase in Ldlr levels or activity post-transcriptionally. Our time course experiments did not show a significant repression of Ldlr levels until 24 h further highlighting the complex nature of Srebf2 regulation.
Our data with TSA treatment also showed a decrease in the levels of gene expression for a variety of apolipoproteins including apoA1, apoA5, apoB, apoC1, apoE, apoL1. This observation highlights the complex relationship of apolipoprotein levels and lipoprotein metabolism. While elevated levels of apoB and reduced levels of apoA1 are associated with increased cardiac disease, serum levels of apoB100 associated VLDL are regulated in turn by Acat2 which stimulates cholesteryl ester secretion into apoBcontaining lipoproteins. Acat inhibitors are being devel-oped as a therapeutic means to lower LDL cholesterol without affecting cholesterol uptake [42,43]. Also apoE deficient mice show high levels of cholesterol and develop spontaneous atherosclerosis while mice with partial or complete deficiency of high-mobility group A2 protein (Hmga2) are able to resist diet-induced obesity [44]. Acat2 inhibition using antisense nucleotides was previously shown to alleviate atherosclerosis in apoB-Ldlr -/-mice [45]. This study also found Acat2 inhibition to be effective in reducing plasma cholesterol, increasing plasma triglycerides, and shifting LDL cholesteryl ester fatty acids to become mainly polyunsaturated. In our study, TSA treatment showed a modest decrease in apoB, apoE and apoA1 in addition to decreased levels of Acat2, Fasn and Hmga2. This indicates that triglyceride metabolism is perturbed by TSA and further studies may be necessary to evaluate the possibility of using TSA and other HDACIs for modulating triglyceride metabolism. Cyp27 has been reported to be regulated by the nuclear receptor subfamily of which PPARγ is a member and levels of both these genes have been found to be high in atherosclerotic lesions. Levels of Cyp27a1 (maximal repression at 24 h) and Pparγ (maximal repression between 9-12 h) were found to be repressed by TSA treatment in both the microarray and qPCR data. This observation adds credence to the potential for development of TSA like HDACIs for atherosclerosis.

Conclusion
Our results show that cholesterol metabolism is significantly down regulated by TSA both directly and indirectly and thus HDACI therapy may be a relatively novel tool to develop for use in controlling cholesterol levels. This study only addresses the effect of TSA treatment on transcript levels of the rate limiting enzymes and transcription factors and further studies evaluating protein expression levels are necessary to derive firm conclusions on regulation of this pathway. Additional studies exploring the different classes of HDACIs with respect to their effects on regulation of the genes in the cholesterol pathway would also help dissect the details of this innovative application for these drugs.

RNA extraction and purification
Both F9 and HepG2 cells were harvested with 4 mL of Trireagent (Molecular Research Center, Inc) and RNA isolation was carried out according to the manufacturer's protocol. Total RNA was purified using the RNeasy cleanup kit and protocol (Qiagen), quantified and then analyzed for degradation on a BioAnalyzer (Agilent).

Hybridization of sample to GeneChip Microarrays
RNA was converted to biotinylated cRNA (complimentary RNA) from oligo-dT-primed cDNA using standard Affymetrix protocols. Biotinylated cRNA was used to probe the MU74Av2 (F9 samples) or HU133 Plus 2.0 (HepG2 samples) Affymetrix GeneChip microarrays. A total of six samples (three controls and three TSA treated) for each cell line were analyzed.

Statistical Analysis
The raw data (CEL files) were imported into GeneSpring software (v7.2) for further analysis. A two-step normalization algorithm was implemented to select differential gene expression in response to TSA samples (ethanol treated samples as baseline). In the first normalization step, a global scaling per chip method was used in which the signal of each gene was divided by the mean intensity (50 th percentile) of the chip. This normalization step was followed by a per gene normalization which divides each gene by the average intensity of that gene in several control samples. Hierarchical clustering was used to organize the data in discrete expression profiles. Selection of statistically significant genes from each expression profile was done using a p-value cut off of ≤ 0.05 with the cross gene error model (CGEM) combined with Welch t-test. The multiple testing correction (Benjamini and Hochberg false discovery rate) was integrated within each test. Additionally we also analyzed the HepG2 data using both  MAS5 (Microarray Suite) as well as RMA (Robust Microarray Analysis) algorithms and genes that passed all criteria from both sets of analyses were used for follow up studies.

Pathway and Functional Cluster Analysis
The differentially expressed genes selected as described above were subjected to functional cluster analysis using MAPPFinder in conjunction with GenMAPP (Gene Microarray Pathway Profiler) 2.0 [19,20]

Quantitative Real-Time PCR
To verify the data obtained from microarrays, 5 μg of total RNA was taken from the same pool of RNA as used for the microarray experiments. The RNA was DNase treated (Ambion) and reverse transcribed to cDNA which served as the template for quantitative PCR (qPCR). Real-time relative qPCR (SYBR Green; Applied Biosystems) was performed in triplicate using a GeneAmp 5700 (F9 samples) or a HT7900 sequence detection system (HepG2 samples) according to the manufacturer's instructions. Primers were specifically designed using Primer express software (Applied Biosystems). 1 μg of cDNA was amplified in 1× SYBR green buffer. PCR conditions were: 10 min at 95°C for AmpliTaq Gold DNA polymerase activation, 45 thermal cycles of 15 sec at 95°C to denature and 1 min at 60°C to anneal and extend. Relative expression levels were analyzed using the method using the GAPDH expression level as a control. Samples treated with TSA were compared to the baseline expression value determined from ethanol treated samples at the respective time points and the fold change is shown. (For primer sequences see Table 5).
For quantitative analysis of the data, C T (threshold-cycle number) values were normalized to those of GAPDH, with use of the ΔΔC T method.