Genome-wide changes accompanying knockdown of fatty acid synthase in breast cancer

Background The lipogenic enzyme fatty acid synthase (FAS) is up-regulated in a wide variety of cancers, and is considered a potential metabolic oncogene by virtue of its ability to enhance tumor cell survival. Inhibition of tumor FAS causes both cell cycle arrest and apoptosis, indicating FAS is a promising target for cancer treatment. Results Here, we used gene expression profiling to conduct a global study of the cellular processes affected by siRNA mediated knockdown of FAS in MDA-MB-435 mammary carcinoma cells. The study identified 169 up-regulated genes (≥ 1.5 fold) and 110 down-regulated genes (≤ 0.67 fold) in response to knockdown of FAS. These genes regulate several aspects of tumor function, including metabolism, cell survival/proliferation, DNA replication/transcription, and protein degradation. Quantitative pathway analysis using Gene Set Enrichment Analysis software further revealed that the most pronounced effect of FAS knockdown was down-regulation in pathways that regulate lipid metabolism, glycolysis, the TCA cycle and oxidative phosphorylation. These changes were coupled with up-regulation in genes involved in cell cycle arrest and death receptor mediated apoptotic pathways. Conclusion Together these findings reveal a wide network of pathways that are influenced in response to FAS knockdown and provide new insight into the role of this enzyme in tumor cell survival and proliferation.


Background
Up-regulation of fatty acid synthase (FAS), the enzyme responsible for the endogenous synthesis of palmitate, is increasingly recognized as a hallmark of cancer [1,2]. While normal cells obtain most fatty acids from circulating lipids, tumor cells have developed an increased reliance on endogenous fatty acid synthesis to satisfy their metabolic needs [2]. This elevation of FAS occurs early in human cancer, is associated with aggressive forms of the disease, and is linked to poor prognosis [3][4][5][6][7]. Conse-quently, strategies to target FAS are becoming increasingly exploited as attractive approaches for cancer therapy.
FAS is critically important for tumor cell survival and function; providing the necessary fatty acids for membrane formation and signal transduction [2]. The targeted knockdown of tumor FAS by small molecule inhibitors or small interfering RNA (siRNA) has been shown to induce both cell cycle arrest and apoptosis in cultured cells and suppresses tumor growth in xenograft bearing mice [8][9][10][11][12][13]. This anti-tumor activity is linked to increased expres-sion of p27 kip1 [9,14] and decreased Akt phosphorylation [15]. FAS inhibition can also lead to the transcriptional suppression of the Her2/neu gene suggesting an active role for FAS in gene regulation [16]. However, detailed understanding of the regulatory mechanisms linking inhibition of FAS to these anti-tumor effects remains elusive.
In the present study, we utilized siRNA, BeadArray technology, and pathway analysis to define the genome-wide changes that take place following knockdown of FAS (12-48 h). We have identified a core set of 279 genes representing the FAS knockdown signature in the MDA-MB-435 mammary carcinoma cell model. Functional classification of these target genes, combined with quantitative pathway analysis, revealed extensive changes in metabolism, cell survival/proliferation, DNA replication/transcription, and ubiquitin dependent protein degradation as a consequence of FAS inhibition. Taken together, our results provide a detailed overview of the anti-tumorigenic signaling network induced in tumor cells by the targeted knockdown of FAS.

Changes in gene expression resulting from knockdown of FAS
MDA-MB-435 mammary carcinoma cells were selected as the model for defining the FAS knockdown signature. For our experiment, four independent siRNA duplexes targeting FAS (FAS #1-#4) were chosen based on there ability to knockdown the enzyme and to induce tumor cell apoptosis after 72 h [17]. Inhibition of FAS by each duplex was verified using 25 nM of siRNA as demonstrated by a decrease in FAS mRNA, protein, and fatty acid biosynthesis after 48 h relative to non-silencing control siRNA (Figure 1a-c). The non-silencing control siRNA was selected based on minimal cross reactivity with known targets and had no impact on FAS expression or activity when compared against Lipofectamine 2000 transfection alone (data not shown). Gene expression profiles were examined on two separate occasions following transfection with FAS siRNA. Treatment times (12, 24, 36 and 48 h) were chosen to capture early and late gene changes associated with the block in cell cycle progression and the advent of apoptosis occurring in response to knockdown of FAS. Abrogation of FAS (>70%) was verified within 12 h of transfection and persisted throughout the 48 h experiment. Significant changes in cell viability were not observed 48 h post-transfection indicating cells were viable during the time of gene analysis (data not shown).
The first approach we used in analyzing the effect of knockdown of FAS on gene expression was to identify target genes modified by at least 3 of the siRNA duplexes. The expression signature of each FAS siRNA was determined by identifying genes significantly changed over the time course using MB statistics and two way ANOVA analysis. For each siRNA duplex only genes with a p-value ≤ 0.05 and at least 1.2 fold changes in both biological replicates at a given time point were maintained for expression analysis, and are published as supporting information online (see Additional Files 1, 2, 3, and 4). The FAS knockdown signature was defined as the overlap of significant gene changes occurring in response to at least 3 of the siRNA duplexes with an average gene expression change 1.5 fold that observed in controls. This allowed us to identify gene changes specifically associated with knockdown of FAS, while at the same time eliminate potential off-target effects unique to the individual siRNA sequences. Using this approach, we identified 279 genes whose expression changed in response to knockdown of FAS (169 genes upregulated and 110 down-regulated; see Additional File 5). Alterations in gene expression occurred as early as 12 h after knockdown of FAS with the majority of genes affected by 24 h (Figure 1d). These results show that inactivation of FAS has a profound affect on gene transcription.
Knockdown of FAS affected genes that regulate a variety of biological processes including cell proliferation, DNA replication, transcription and apoptosis (Table 1; see Additional File 5). We found that the abrogation of FAS enhanced the expression of several anti-proliferative genes including CCNG1, CDKN1A (p21cip1), SPRY2 and SPRY4. This coincided with up-regulation of the cell surface apoptosis gene ANXA1, the mitochondrial apoptosis pathway gene APAF1, and the death receptors TNFRSF10B (TNF-related apoptosis-inducing ligand (TRAIL) receptor) and TNFRSF21 (TNF-α receptor). In addition, results showed knockdown of FAS up-regulates SGPL1, a key enzyme in sphingosine-1-phosphate metabolism whose increased expression also leads to accumulation of the pro-apoptotic signaling molecule ceramide [18]. This finding is consistent with reports that FAS inhibition upregulates ceramide synthesis, which has been found to be a necessary step leading to the induction of tumor cell apoptosis following the loss of FAS [17].
Several of the genes significantly altered by knockdown of FAS are known to play a role in the regulation of lipid metabolism. Inactivation of FAS increased the expression of INSIG1, and the leptin receptors, LEPR and LEPROTL1, which block transcription, proteolytic cleavage and transcriptional activation of the sterol regulatory elementbinding protein (SREBP) family of transcription factors that promote lipid biosynthesis (Table 1) [19,20]. Additionally, loss of FAS led to down-regulation of ID2 and ID3, which are dominant negative helix-loop-helix (HLH) proteins that bind SREBP 1c and functionally repress FAS promoter activity [21]. ID2 and ID3 also repress expression of the cyclin-dependent kinase inhibitor p21 cip1 FAS knockdown generates specific and time-dependent gene expression patterns     [22,23], thus, as expected the suppression of these genes in response to knockdown of FAS was coupled with upregulation of the p21 cip1 gene (Table 1). These data provide evidence of a direct link between regulation of FAS and control of cell cycle progression.
Another important biological category affected by knockdown of FAS is regulation of protein ubiquitination (Table 1). We found that knockdown of FAS altered the expression of several E2 ubiquitin conjugation enzymes (2 up-regulated) and E3 ubiquitin ligases (4 down-regulated and 6 up-regulated) which function to target proteins for degradation by the proteosome [24]. In addition, we identified up-regulation of 2 splice variants of UBE2V1, which are catalytically inactive E2 enzymes [25]. UBE2V1 variants result from the co-transcription of UBE2V1 with the neighboring upstream gene, Kua, which shares sequence consensus to fatty acid hydrolases (vari- . This targeting of ubiquitination enzymes indicates that the knockdown of FAS affects the tumor cell proteome not only through changes in transcription, but also on the post-translational level.

Knockdown of FAS leads to widespread changes in pathways involved in metabolism
The second approach we used to analyze the gene array results was to identify gene pathways coordinately up-regulated or down-regulated by knockdown of FAS using Gene Set Enrichment Analysis (GSEA). GSEA is a unique application that coordinates small changes in gene expression from a large number of functionally related genes in order to identify pathways that are significantly changed [26]. As expected, fatty acid metabolism (NES = -1.808, NOM p-val = 0.002, FDR q-val = 0.085) was significantly suppressed by knockdown of FAS ( Figure 2). Fatty acid elongation (NES = -1.537, NOM p-val = 0.037, FDR q-val = 0.394) was also reduced, which is consistent with the loss of available palmitate.

Knockdown of FAS up-regulates pathways involved in cell cycle arrest and apoptosis
The anti-tumorigenic effects of knockdown of FAS were traced to an up-regulation in genes that modify tumor proliferation ( Figure 3) pathways consistent with our previous finding that FAS regulates the G 1 checkpoint through these pathways [9]. Additionally, knockdown of FAS elevated the expression of several cyclin-dependent kinase inhibitors (p16INK4 (CDKN2A), p15INK4 (CDKN2B) and p21 Cip1 ) which function as negative regulators of the Rb pathway [27], demonstrating FAS induces a high degree of regulation over the G 1 /S transition (Figure 3). Results also showed that knockdown of FAS elevated genes involved in the G 2 pathway including the checkpoint regulators CHEK2, WEE1 and PLK1 (Table 1; Figure 3). These data indicate that FAS controls tumor proliferation by regulating several aspects of the cell division cycle.
The apoptotic effects of knockdown of FAS were recently shown to be blocked by co-treatment with the pan-caspase inhibitor z-VAD-fmk [28], suggesting a role for caspases in mediating FAS induced tumor cell death. Here, we found up-regulation of caspase 7 and caspase 8 (receptor-mediated apoptosis) in response to knockdown of FAS ( Figure 3). This was accompanied by enhanced expression of pathways and genes that activate caspase 8

Discussion
The present study provides a comprehensive overview of the genomic changes accompanying knockdown of tumor FAS. To our knowledge, this is the first report demonstrating that targeted knockdown of a metabolic enzyme can influence overall energy producing pathways within the tumor cell. The major findings of this work are as follows: first, siRNA mediated knockdown of FAS affects the transcription of genes involved in tumor cell energy metabo-

Down-regulation of metabolic pathways by targeted knockdown of FAS
Up-regulation of cell cycle arrest and cell death pathways in response to knockdown of FAS lism as demonstrated by down-regulation in lipid metabolism, glycolysis, krebs-TCA cycle and oxidative phosphorylation pathways. Second, we show that loss of FAS is anti-tumorigenic by up-regulating cell cycle arrest and death receptor mediated apoptosis pathways. Third, the inhibition of FAS leads to changes in genes that regulate transcription and ubiquitin-dependent protein degradation. Altogether, these findings provide novel insight into the cellular processes affected by knockdown of FAS in tumor cells.
A vital assumption in siRNA-mediated gene silencing is that knockdown of the target gene is exclusive, and without off-target consequences. Yet, increasing evidence has suggested siRNA can induce off-target effects that are unique to the individual siRNA duplex, but unrelated to the identity of the target gene [29]. These siRNA duplex specific off-target effects are thought to arise from unintended sequence homology with the off-target transcript. In this study, we addressed the concern of potential siRNA off-target effects by examining gene expression signatures in response to four siRNA duplexes targeting different regions of the FAS gene and focused our analysis on the common gene changes that were associated with knockdown of the enzyme. This eliminated off-target effects and enabled us to study gene changes specifically associated with knockdown of FAS. Altogether, evidence supporting the role of FAS inhibition in inducing the described genomic effects is as follows: first, we showed that FAS siRNA strongly abrogates the expression of the FAS gene, reduces FAS protein levels and blocks the biosynthetic activity of the enzyme, as indicated by reduced [ 14 C]malonyl-CoA incorporation into palmitate. Second, pathway analysis demonstrated that fatty acid biosynthesis and metabolism are among the most significant pathways down-regulated in response to FAS siRNA. Third, up-regulation of the p27 and Rb gene pathways is consistent with our prior work showing knockdown of FAS elicits a G 1 /S arrest through these pathways [9]. We also found increased expression of the p21 Cip1 gene which is up-regulated at the protein level in response to knockdown of FAS [14,30]. Like p27, p21 Cip1 acts as a negative regulator of the G 1 /S transition by inactivating the cyclin-dependent kinases that function to phosphorylate Rb and promote S phase entry [31].
Knockdown of FAS by siRNA or the small molecule inhibitor Orlistat induces apoptosis by activating caspase 8 (Knowles LM and Smith JW, unpublished observations). This is consistent with the current observation that siRNA targeting FAS up-regulates TNF-α and TRAIL death receptor mediated apoptosis pathways. In addition to TRAIL, cell death through FAS inhibition has also been linked to up-regulation of BNIP3 and DAP Kinase 2 [17]. These two proapoptotic genes were detected by microarray analysis 72 h after FAS inhibition, which explains why we did not observe similar changes as our study was limited to 48 h. A possible mechanism accounting for the ability of FAS inhibition to induce these apoptotic effects is the accumulation of ceramide, which has been shown to induce these death receptor pathways [32]. Ceramide accumulation is the result of malonyl-CoA build-up following inhibition of FAS [17]. Importantly, we have found that knockdown of FAS up-regulates SGPL1. Ceramide generation and promotion of cellular apoptosis have been observed in response to increased expression of SGPL1 [18], suggesting this enzyme may represent a leverage point by which knockdown of FAS influences ceramide production.
In addition to changes in proliferation and survival, knockdown of FAS also suppresses genes in the glycolysis, krebs-TCA cycle and oxidative phosphorylation metabolic pathways suggesting an overall reduction in cellular energy metabolism. Consistent with the down-regulation of oxidative phosphorylation gene expression, Chajès et al. [33] report impaired mitochondrial function and generation of reactive oxygen species in response to FAS inhibition. Targeted inactivation of ATP citrate lyase, a lipogenic enzyme which functions immediately upstream FAS, has the opposite effect by stimulating mitochondrial function [34]. This is interesting as we have found that knockdown of FAS, but not of ATP citrate lyase, induces apoptosis (Knowles LM and Smith JW, unpublished observations); raising the question of whether the changes in cell survival and metabolism are linked together. Indeed, growing evidence suggests TNF-α may influence metabolism, as it can induce the down-regulation of genes involved in oxidative phosphorylation [35], inhibition of electron transport chain activity [36] and generation of reactive oxygen species [37]. Identifying whether induction of apoptosis and reduced energy metabolism are connected (via TNF-α or another molecular target) will be an important direction for future investigations, and may prove critical for understanding the anti-tumorigenic benefits of a knockdown in tumor FAS.
The overall goal of our study was to define the regulatory changes that take place following inhibition of tumor FAS. To this end, we have identified a number of target genes specifically modified by knockdown of FAS. Of particular interest among these are the dominant negative HLH transcription repressors ID2 and ID3, which function to bind and inactivate the basic HLH family of transcription factors [38], and therefore, are recognized for their ability to block cell differentiation and promote proliferation [39]. Like FAS, ID over-expression correlates with tumor development in an array of cancers [40,41], and knockdown of ID2/ID3 induces growth arrest and apoptosis [39,42]. Our finding that knockdown of FAS leads to the suppression of transcripts encoding ID2 and ID3 provides evidence of a link between FAS and ID expression. This link is further supported by observations that ID2/ID3 proteins can also suppress FAS gene expression by inactivating the basic HLH transcription factor SREBP1 c that controls FAS expression [21]. Further characterization of the molecular interplay between FAS and ID is needed to more fully understand the role of these in tumorigenesis.

Conclusion
In conclusion, this study provides a comprehensive database of genomic changes that occur in response to knockdown of FAS and confirms that FAS is central not only for tumor cell metabolism but also for tumor cell signaling. Further characterization of these signals on the protein level will offer a more complete understanding of the cellular targets and pathways affected by knockdown of FAS and may lead to new therapeutic strategies for tumor prevention and treatment.

Real-time quantitative PCR
Total RNA was isolated from adherent cells using the Qiagen RNeasy kit (Qiagen, Inc, Valencia, CA, USA). Genomic DNA contamination was removed using RNase-Free DNase (Qiagen). RNA purity was assessed by A 260 / A 280 absorption and RNA integrity was verified by agarose gel electrophoresis. cDNA was synthesized from 4 μg RNA using the Superscript III First-Strand Synthesis System (Invitrogen). Real-time quantitative PCR reactions were performed using Power SYBR ® Green PCR Master Mix (Applied Biosystems, Foster City, CA, USA) on a Stratagene Mx3000p QPCR System (La Jolla, CA, USA). cDNA was amplified using forward and reverse primers for FAS (5'-AACTCCATGTTTGGTGTTTG-3' and 5'-CACATGCG-GTTTAATTGTG-3') and normalized to the housekeeping gene P0 (5'CAAGACTGGAGACAAAGTGG-3' and 5'AATCTGCAGACAGACACTGG-3'). The formation of a single PCR product was verified using automated melting curve analysis.

Western blot analysis
Adherent cell populations were lysed in 2× SDS sample buffer, separated by SDS-PAGE and transferred to nitrocellulose as previously described [9]. Membranes were probed for FAS and β-tubulin protein expressions using anti-FAS (BD Transduction Laboratories, San Jose, CA, USA) and anti-β-tubulin (Upstate Biotechnology, Lake Placid, NY, USA) antibodies. Immunoreactivity was detected using anti-mouse IgG conjugated peroxidase and visualized by enhanced chemiluminescence.

FAS activity
FAS activity was determined by measuring the incorporation of [ 14 C]malonyl-CoA into cellular fatty acids as previously described [9]. Briefly, MDA-MB-435 cells were lysed in 20 mM Tris (pH 7.5), 1 mM dithiothreitol, and 1 mM EDTA by sonication. Insoluble material was removed by centrifugation (14,000 rpm) for 15 min at 4°C. The resulting lysates (70 ug) were incubated in reaction buffer containing 500 μM NADPH, 166.6 μM acetyl CoA, 100 mM KCl and 0.4 μCi [ 14 C]malonyl-CoA (GE Healthcare, Piscataway, NJ, USA) for 25 min at 37°C. Reactions were chased with 25 nM cold malonyl-CoA for 15 min and terminated by the addition of chloroform:methanol (1:1). The chloroform extract was dried under N 2 and extracted with water-saturated butanol. The butanol extract was evaporated under N 2 , and labeled fatty acids were detected using a Beckman model LS 5000CE Scinillation Counter (Beckman Coulter, Inc., Fullerton, CA, USA).

BeadArray gene expression
MDA-MB-435 cells were transfected with siRNA on two separate occasions to generate independent measurements for each treatment, and monitored for whole genome expression at 12 h intervals over a 48 h time period. Total RNA was isolated and examined for purity/ integrity as described above. Labeled cRNA was prepared from 500 ng RNA using the Illumina ® RNA Amplification Kit from Ambion (Austin, TX, USA). The labeled cRNA (1500 ng) was hybridized overnight at 55°C to the Sentrix ® HumanRef-6 Expression BeadChip (>46,000 gene transcripts; Illumina, San Diego, CA, USA) according to the manufacturer's instructions. BeadChips were subsequently washed and developed with fluorolink streptavidin-Cy3 (GE Healthcare). BeadChips were scanned with an Illumina BeadArray Reader and hybridization efficiency was monitored using BeadStudio software (Illumina) to measure internal controls built into the Illumina system.

Expression data analysis
Expression data was filtered to identify genes expressed on the BeadChip at >0.99 confidence (16,585 genes), and normalized per chip using non-linear Normalize.quantiles or VSN normalization to remove array to array variability (Bioconductor Project [43]). The efficiency of both normalization methods were validated by measuring gene to gene correlations across BeadChips according to the methods of Ploner et al. [44] using R 2.0.1 software (R Development Core Team, Vienna, Austria).
Differentially expressed genes that changed over time in response to each of the four FAS siRNA duplex were identified using two independent methods. In the first method, per chip normalized Normalize.quantiles data was imported into the GeneSpring GX 7.3.1 software package (Agilent Technologies, Palo Alto, CA, USA), and per gene normalized to the appropriate biological control for each time point. Genes that significantly changed over time in response to each FAS siRNA duplex were compared against control using a two-way ANOVA with a Benjamini & Hochberg FDR of 0.05. In the second method, VSN normalized data was imported into the Bioconductor Timecourse package [43]. Genes regulated over time by each FAS siRNA duplex in both biological replicates were determined using mb. long() statistics, with an mb>0 used as criterion for the identification of significantly regulated genes. Genes identified as significant using these two complementary approaches were combined, and genes with a p-value ≤ 0.05 and with 1.2 fold changes in both the biological replicates of each duplex at a given time point, were maintained for expression analysis. Genes identified as significantly modified by each of the 4 FAS siRNA duplexes were then overlapped using venn diagrams to identify a set of genes modified by 3 or more of the siRNA treatments. Expression profiles of the core gene list were averaged, and genes modified by 1.5 fold were visually examined and ordered by hierarchical clustering using the Pearson Correlation similarity measurement (GeneSpring).

Gene set enrichment analysis
For pathway analysis, data from the 16,585 genes expressed in MDA-MB-435 cells was exported from Gene-Spring, filtered for duplicate symbols and analyzed using GSEA software (Broad Institute) according to published methods [26]. Briefly, data was overlapped on 522 gene sets (s2.symbols.gmt) downloaded with the GSEA package and measured for the enrichment of genes at the top or bottom of the gene list to determine their correlation with the gene set's phenotype. The GSEA parameters used included: metric for ranking genes, signal2noise; enrichment statistic, weighted; permutation type, phenotype; permutation number, 1000; and gene set size restrictions, 4 minimum, 500 maximum. Gene sets significantly modified by FAS siRNA treatment were identified using a nominal p-value < 0.05 and a multiple hypothesis testing FDR < 0.25. NES represents the enrichment of genes in the designated GSEA gene set, ranked according to the overrepresentation of genes at the top or bottom of the list, normalized to gene set size.