Skip to main content

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

Abstract

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–7]. Consequently, 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–13]. This anti-tumor activity is linked to increased expression of p27kip1 [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.

Results

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).

Figure 1
figure 1

FAS knockdown generates specific and time-dependent gene expression patterns. (a-c) Target mRNA, protein and fatty acid biosynthesis knockdown by FAS siRNA duplexes. MDA-MB-435 tumor cells were exposed to 25 nM of four different siRNA duplexes targeted against FAS (FAS #1-#4) or non-silencing control siRNA for 48 h. Efficiency of FAS knockdown was determined by measuring FAS mRNA (a), FAS protein (b) and fatty acid biosynthesis (c). β-tubulin served as a loading control for FAS protein expression. Values are the mean ± SE of two replicates per treatment. (d) Expression profile represents 279 genes differentially modified by 1.5 fold in response to knockdown of FAS (shown in hours on the X-axis). Horizontal lines represent the average expression of individual genes modified by at least 3 of the FAS siRNA duplexes. Red and green indicate increased and decreased expression, respectively, relative to non-silencing control siRNA.

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 up-regulated 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 (p21cip 1), 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 up-regulates 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].

Table 1 Selected FAS target genes involved in proliferation/apoptosis, metabolism, transcription and protein ubiquitination.

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 element-binding 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 p21cip1 [22, 23], thus, as expected the suppression of these genes in response to knockdown of FAS was coupled with up-regulation of the p21cip 1gene (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 (variant 2, Kua-UEV fusion gene [GI_40806191-I]; variant 4, Kua gene [GI_40806192-I]). 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.

Figure 2
figure 2

Down-regulation of metabolic pathways by targeted knockdown of FAS. The figure displays select pathways found to be down-regulated in response to FAS siRNA treatment compared to non-silencing control siRNA. Significance was determined using a nominal p-value < 0.05 or FDR < 0.250. The expression levels of the genes significantly modified in the pathway are coded colorimetrically: red, high expression; blue, low expression. FAS siRNA treatments for each time point are ordered as follows (a and b indicate different biological replicates): FAS #1a, FAS #1b, FAS #2a, FAS #2b, FAS #3a, FAS #3b, FAS #4a, FAS #4b. For a complete list of all pathways down-regulated by knockdown of FAS, see Additional File 6.

The effects of knockdown of FAS on the expression of metabolic genes extended to pathways other than fatty acid metabolism (Figure 2). Inactivation of FAS down-regulated both the glycolysis/gluconeogenesis (NES = -1.560, NOM p-val = 0.029, FDR q-val = 0.377) and krebs-TCA cycle (NES = -1.737, NOM p-val = 0.014, FDR q-val = 0.094) pathways. Down-regulation of the krebs-TCA cycle was consistent with an overall suppression in mitochondrial genes involved in energy metabolism and oxidative phosphorylation (GO 005739; NES = -1.515, NOM p-val = 0.039, FDR q-val = 0.365). Furthermore, this down-regulation in genes regulating glucose utilization was coupled with a reduction in insulin signaling (NES = -1.635, NOM p-val = 0.005, FDR q-val = 0.280). This suggests that inhibition of FAS may restrict the ability of tumor cells to produce the energy necessary to thrive.

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). GSEA analysis showed that knockdown of FAS up-regulated genes involved in DNA damage signaling (NES = 1.809, NOM p-val = 0.000, FDR q-val = 0.068) and cell cycle arrest (NES = 1.504, NOM p-val = 0.036, FDR q-val = 0.163). Induction of cell cycle arrest was coupled with up-regulation of the p27 (NES = 1.582, NOM p-val = 0.009, FDR q-val = 0.128) and Rb (NES = 1.581, NOM p-val = 0.010, FDR q-val = 0.126) pathways consistent with our previous finding that FAS regulates the G1 checkpoint through these pathways [9]. Additionally, knockdown of FAS elevated the expression of several cyclin-dependent kinase inhibitors (p16INK4 (CDKN2A), p15INK4 (CDKN2B) and p21Cip 1) which function as negative regulators of the Rb pathway [27], demonstrating FAS induces a high degree of regulation over the G1/S transition (Figure 3). Results also showed that knockdown of FAS elevated genes involved in the G2 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.

Figure 3
figure 3

Up-regulation of cell cycle arrest and cell death pathways in response to knockdown of FAS. The figure displays select pathways found to be up-regulated in response to FAS siRNA treatment compared to non-silencing control siRNA. Significance was determined using a nominal p-value < 0.05 or FDR < 0.250. The expression levels of the genes significantly modified in the pathway are coded colorimetrically: red, high expression; blue, low expression. FAS siRNA treatments for each time point are ordered as follows (a and b indicate different biological replicates): FAS #1a, FAS #1b, FAS #2a, FAS #2b, FAS #3a, FAS #3b, FAS #4a, FAS #4b. For a complete list of all pathways up-regulated by knockdown of FAS, see Additional File 7.

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 such as tumor necrosis factor (tnfr1, NES = 1.408, NOM p-val = 0.038, FDR q-val = 0.230; tnfr2, NES = 1.788, NOM p-val = 0.002, FDR q-val = 0.070; Figure 3). GSEA results also demonstrated compensatory up-regulation of NF-κB (NF-κB Induced, NES = 1.5, NOM p-val = 0.022, FDR q-val = 0.166) and the MAP kinase (NES = 1.539, NOM p-val = 0.016, FDR q-val = 0.133) survival pathways in line with previous reports (see Additional File 7; [14]).

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 metabolism 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 [14C]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 G1/S arrest through these pathways [9]. We also found increased expression of the p21Cip 1gene which is up-regulated at the protein level in response to knockdown of FAS [14, 30]. Like p27, p21Cip1 acts as a negative regulator of the G1/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 up-stream 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.

Methods

Cell line and culture conditions

The MDA-MB-435 mammary tumor cell line was obtained from Janet Price at the University of Texas Southwestern. MDA-MB-435 cells were maintained in Minimum Essential Medium Eagle (MEM) with Earl's Salts (Mediatech, Inc., Herndon, VA, USA) supplemented with 10% fetal bovine serum (Irvine Scientific, Santa Ana, CA, USA), 2 mM L-glutamine (Invitrogen Life Technologies, Inc., Carlsbad, CA, USA), MEM vitamins (Invitrogen Life Technologies, Inc.), nonessential amino acids (Mediatech, Inc.) and antibiotics (Omega Scientific, Inc., Tarzana, CA, USA). Cells were grown at 37°C under a humidified, 5% CO2 atmosphere.

FAS gene silencing using siRNA

Four individual FAS siRNA sequences corresponding to 5'-GAGCGUAUCUGUGAGAAACUU-3' (nucleotides 6241–6259; FAS#1), 5'-GACGAGAGCACCUUUGAUGUU-3' (nucleotides 1741–1749; FAS#2), 5'-UGACAUCGUCCAUUCGUUUUU-3' (nucleotides 1758–1776; FAS#3), 5'-UGACAUCGUCCAUUCGUUUUU-3' (nucleotides 6236–6254; FAS#4) were designed and synthesized by Dharmacon (Lafayette, CO, USA). MDA-MB-435 cells were plated at 7.81 × 103 cells/cm2 in 10-cm2 plates and grown for 24 h prior to transfection with 25 nM FAS#1, FAS#2, FAS#3, FAS#4 or non-silencing control duplex #2 (D-001210-02) siRNA in Opti-MEM medium (Invitrogen) using Lipofectamine 2000 reagent (Invitrogen). After 5 h, transfection medium was replaced with normal culture medium and cells were grown for 12–48 h post-transfection.

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 A260/A280 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'-CACATGCGGTTTAATTGTG-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 [14C]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 [14C]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 N2 and extracted with water-saturated butanol. The butanol extract was evaporated under N2, 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 GeneSpring, 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.

Abbreviations

FAS:

fatty acid synthase

GSEA:

Gene Set Enrichment Analysis

HLH:

helix-loop-helix

siRNA:

small interfering RNA

SREBP:

sterol regulatory element-binding protein

TRAIL:

TNF-related apoptosis-inducing ligand.

References

  1. Wakil SJ: Fatty acid synthase, a proficient multifunctional enzyme. Biochemistry. 1989, 28 (11): 4523-4530. 10.1021/bi00437a001.

    Article  CAS  PubMed  Google Scholar 

  2. Kuhajda FP: Fatty-acid synthase and human cancer: new perspectives on its role in tumor biology. Nutrition. 2000, 16 (3): 202-208. 10.1016/S0899-9007(99)00266-X.

    Article  CAS  PubMed  Google Scholar 

  3. Piyathilake CJ, Frost AR, Manne U, Bell WC, Weiss H, Heimburger DC, Grizzle WE: The expression of fatty acid synthase (FASE) is an early event in the development and progression of squamous cell carcinoma of the lung. Hum Pathol. 2000, 31 (9): 1068-1073. 10.1053/hupa.2000.9842.

    Article  CAS  PubMed  Google Scholar 

  4. Swinnen JV, Roskams T, Joniau S, Van Poppel H, Oyen R, Baert L, Heyns W, Verhoeven G: Overexpression of fatty acid synthase is an early and common event in the development of prostate cancer. Int J Cancer. 2002, 98 (1): 19-22. 10.1002/ijc.10127.

    Article  CAS  PubMed  Google Scholar 

  5. Alo PL, Visca P, Marci A, Mangoni A, Botti C, Di Tondo U: Expression of fatty acid synthase (FAS) as a predictor of recurrence in stage I breast carcinoma patients. Cancer. 1996, 77 (3): 474-482. 10.1002/(SICI)1097-0142(19960201)77:3<474::AID-CNCR8>3.0.CO;2-K.

    Article  CAS  PubMed  Google Scholar 

  6. Innocenzi D, Alo PL, Balzani A, Sebastiani V, Silipo V, La Torre G, Ricciardi G, Bosman C, Calvieri S: Fatty acid synthase expression in melanoma. J Cutan Pathol. 2003, 30 (1): 23-28. 10.1034/j.1600-0560.2003.300104.x.

    Article  CAS  PubMed  Google Scholar 

  7. Gansler TS, Hardman W, Hunt DA, Schaffel S, Hennigar RA: Increased expression of fatty acid synthase (OA-519) in ovarian neoplasms predicts shorter survival. Hum Pathol. 1997, 28 (6): 686-692. 10.1016/S0046-8177(97)90177-5.

    Article  CAS  PubMed  Google Scholar 

  8. Pizer ES, Chrest FJ, DiGiuseppe JA, Han WF: Pharmacological inhibitors of mammalian fatty acid synthase suppress DNA replication and induce apoptosis in tumor cell lines. Cancer Res. 1998, 58 (20): 4611-4615.

    CAS  PubMed  Google Scholar 

  9. Knowles LM, Axelrod F, Browne CD, Smith JW: A fatty acidsynthase blockade induces tumor cell-cycle arrest by down-regulating Skp2. J Biol Chem. 2004, 279 (29): 30540-30545. 10.1074/jbc.M405061200.

    Article  CAS  PubMed  Google Scholar 

  10. Kridel SJ, Axelrod F, Rozenkrantz N, Smith JW: Orlistat is a novel inhibitor of fatty acid synthase with antitumor activity. Cancer Res. 2004, 64 (6): 2070-2075. 10.1158/0008-5472.CAN-03-3645.

    Article  CAS  PubMed  Google Scholar 

  11. Pizer ES, Wood FD, Heine HS, Romantsev FE, Pasternack GR, Kuhajda FP: Inhibition of fatty acid synthesis delays disease progression in a xenograft model of ovarian cancer. Cancer Res. 1996, 56 (6): 1189-1193.

    CAS  PubMed  Google Scholar 

  12. Pizer ES, Thupari J, Han WF, Pinn ML, Chrest FJ, Frehywot GL, Townsend CA, Kuhajda FP: Malonyl-coenzyme-A is a potential mediator of cytotoxicity induced by fatty-acid synthase inhibition in human breast cancer cells and xenografts. Cancer Res. 2000, 60 (2): 213-218.

    CAS  PubMed  Google Scholar 

  13. De Schrijver E, Brusselmans K, Heyns W, Verhoeven G, Swinnen JV: RNA interference-mediated silencing of the fatty acid synthase gene attenuates growth and induces morphological changes and apoptosis of LNCaP prostate cancer cells. Cancer Res. 2003, 63 (13): 3799-3804.

    CAS  PubMed  Google Scholar 

  14. Menendez JA, Mehmi I, Atlas E, Colomer R, Lupu R: Novel signaling molecules implicated in tumor-associated fatty acid synthase-dependent breast cancer cell proliferation and survival: Role of exogenous dietary fatty acids, p53-p21WAF1/CIP1, ERK1/2 MAPK, p27KIP1, BRCA1, and NF-kappaB. Int J Oncol. 2004, 24 (3): 591-608.

    CAS  PubMed  Google Scholar 

  15. Wang HQ, Altomare DA, Skele KL, Poulikakos PI, Kuhajda FP, Di Cristofano A, Testa JR: Positive feedback regulation between AKT activation and fatty acid synthase expression in ovarian carcinoma cells. Oncogene. 2005, 24 (22): 3574-3582. 10.1038/sj.onc.1208463.

    Article  CAS  PubMed  Google Scholar 

  16. Menendez JA, Vellon L, Mehmi I, Oza BP, Ropero S, Colomer R, Lupu R: Inhibition of fatty acid synthase (FAS) suppresses HER2/neu (erbB-2) oncogene overexpression in cancer cells. Proc Natl Acad Sci USA. 2004, 101 (29): 10715-10720. 10.1073/pnas.0403390101.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  17. Bandyopadhyay S, Zhan R, Wang Y, Pai SK, Hirota S, Hosobe S, Takano Y, Saito K, Furuta E, Iiizumi M: Mechanism of apoptosis induced by the inhibition of fatty acid synthase in breast cancer cells. Cancer Res. 2006, 66 (11): 5934-5940. 10.1158/0008-5472.CAN-05-3197.

    Article  CAS  PubMed  Google Scholar 

  18. Reiss U, Oskouian B, Zhou J, Gupta V, Sooriyakumaran P, Kelly S, Wang E, Merrill AH, Saba JD: Sphingosine-phosphate lyase enhances stress-induced ceramide generation and apoptosis. J Biol Chem. 2004, 279 (2): 1281-1290. 10.1074/jbc.M309646200.

    Article  CAS  PubMed  Google Scholar 

  19. Yang T, Espenshade PJ, Wright ME, Yabe D, Gong Y, Aebersold R, Goldstein JL, Brown MS: Crucial step in cholesterol homeostasis: sterols promote binding of SCAP to INSIG-1, a membrane protein that facilitates retention of SREBPs in ER. Cell. 2002, 110 (4): 489-500. 10.1016/S0092-8674(02)00872-3.

    Article  CAS  PubMed  Google Scholar 

  20. Nogalska A, Sucajtys-Szulc E, Swierczynski J: Leptin decreases lipogenic enzyme gene expression through modification of SREBP-1c gene expression in white adipose tissue of aging rats. Metabolism. 2005, 54 (8): 1041-1047. 10.1016/j.metabol.2005.03.007.

    Article  CAS  PubMed  Google Scholar 

  21. Moldes M, Boizard M, Liepvre XL, Feve B, Dugail I, Pairault J: Functional antagonism between inhibitor of DNA binding (Id) and adipocyte determination and differentiation factor 1/sterol regulatory element-binding protein-1c (ADD1/SREBP-1c) trans-factors for the regulation of fatty acid synthase promoter in adipocytes. Biochem J. 1999, 344 (Pt 3): 873-880. 10.1042/0264-6021:3440873.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  22. Matsumura ME, Lobe DR, McNamara CA: Contribution of the helix-loop-helix factor Id2 to regulation of vascular smooth muscle cell proliferation. J Biol Chem. 2002, 277 (9): 7293-7297. 10.1074/jbc.M108986200.

    Article  CAS  PubMed  Google Scholar 

  23. Forrest ST, Taylor AM, Sarembock IJ, Perlegas D, McNamara CA: Phosphorylation regulates Id3 function in vascular smooth muscle cells. Circ Res. 2004, 95 (6): 557-559. 10.1161/01.RES.0000142735.67542.5a.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  24. Nakayama KI, Nakayama K: Ubiquitin ligases: cell-cycle control and cancer. Nat Rev Cancer. 2006, 6 (5): 369-381. 10.1038/nrc1881.

    Article  CAS  PubMed  Google Scholar 

  25. Thomson TM, Lozano JJ, Loukili N, Carrio R, Serras F, Cormand B, Valeri M, Diaz VM, Abril J, Burset M: Fusion of the human gene for the polyubiquitination coeffector UEV1 with Kua, a newly identified gene. Genome Res. 2000, 10 (11): 1743-1756. 10.1101/gr.GR-1405R.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  26. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102 (43): 15545-15550. 10.1073/pnas.0506580102.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  27. Sherr CJ, McCormick F: The RB and p53 pathways in cancer. Cancer Cell. 2002, 2 (2): 103-112. 10.1016/S1535-6108(02)00102-2.

    Article  CAS  PubMed  Google Scholar 

  28. Brusselmans K, De Schrijver E, Verhoeven G, Swinnen JV: RNA interference-mediated silencing of the acetyl-CoA-carboxylase-alpha gene induces growth inhibition and apoptosis of prostate cancer cells. Cancer Res. 2005, 65 (15): 6719-6725. 10.1158/0008-5472.CAN-05-0571.

    Article  CAS  PubMed  Google Scholar 

  29. Jackson AL, Bartz SR, Schelter J, Kobayashi SV, Burchard J, Mao M, Li B, Cavet G, Linsley PS: Expression profiling reveals off-target gene regulation by RNAi. Nat Biotechnol. 2003, 21 (6): 635-637. 10.1038/nbt831.

    Article  CAS  PubMed  Google Scholar 

  30. Li JN, Gorospe M, Chrest FJ, Kumaravel TS, Evans MK, Han WF, Pizer ES: Pharmacological inhibition of fatty acid synthase activity produces both cytostatic and cytotoxic effects modulated by p53. Cancer Res. 2001, 61 (4): 1493-1499.

    CAS  PubMed  Google Scholar 

  31. Rank KB, Evans DB, Sharma SK: The N-terminal domains of cyclin-dependent kinase inhibitory proteins block the phosphorylation of cdk2/Cyclin E by the CDK-activating kinase. Biochem Biophys Res Commun. 2000, 271 (2): 469-473. 10.1006/bbrc.2000.2648.

    Article  CAS  PubMed  Google Scholar 

  32. Bollinger CR, Teichgraber V, Gulbins E: Ceramide-enriched membrane domains. Biochim Biophys Acta. 2005, 1746 (3): 284-294. 10.1016/j.bbamcr.2005.09.001.

    Article  CAS  PubMed  Google Scholar 

  33. Chajes V, Cambot M, Moreau K, Lenoir GM, Joulin V: Acetyl-CoA carboxylase alpha is essential to breast cancer cell survival. Cancer Res. 2006, 66 (10): 5287-5294. 10.1158/0008-5472.CAN-05-1489.

    Article  CAS  PubMed  Google Scholar 

  34. Hatzivassiliou G, Zhao F, Bauer DE, Andreadis C, Shaw AN, Dhanak D, Hingorani SR, Tuveson DA, Thompson CB: ATP citrate lyase inhibition can suppress tumor cell growth. Cancer Cell. 2005, 8 (4): 311-321. 10.1016/j.ccr.2005.09.008.

    Article  CAS  PubMed  Google Scholar 

  35. Dahlman I, Forsgren M, Sjogren A, Nordstrom EA, Kaaman M, Naslund E, Attersand A, Arner P: Downregulation of Electron Transport Chain Genes in Visceral Adipose Tissue in Type 2 Diabetes Independent of Obesity and Possibly Involving Tumor Necrosis Factor-{alpha}. Diabetes. 2006, 55 (6): 1792-1799. 10.2337/db05-1421.

    Article  CAS  PubMed  Google Scholar 

  36. Schulze-Osthoff K, Bakker AC, Vanhaesebroeck B, Beyaert R, Jacob WA, Fiers W: Cytotoxic activity of tumor necrosis factor is mediated by early damage of mitochondrial functions. Evidence for the involvement of mitochondrial radical generation. J Biol Chem. 1992, 267 (8): 5317-5323.

    CAS  PubMed  Google Scholar 

  37. Lo YY, Cruz TF: Involvement of reactive oxygen species in cytokine and growth factor induction of c-fos expression in chondrocytes. J Biol Chem. 1995, 270 (20): 11727-11730. 10.1074/jbc.270.20.11727.

    Article  CAS  PubMed  Google Scholar 

  38. Benezra R, Davis RL, Lockshon D, Turner DL, Weintraub H: The protein Id: a negative regulator of helix-loop-helix DNA binding proteins. Cell. 1990, 61 (1): 49-59. 10.1016/0092-8674(90)90214-Y.

    Article  CAS  PubMed  Google Scholar 

  39. Barone MV, Pepperkok R, Peverali FA, Philipson L: Id proteins control growth induction in mammalian cells. Proc Natl Acad Sci USA. 1994, 91 (11): 4985-4988. 10.1073/pnas.91.11.4985.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  40. Sikder HA, Devlin MK, Dunlap S, Ryu B, Alani RM: Id proteins in cell growth and tumorigenesis. Cancer Cell. 2003, 3 (6): 525-530. 10.1016/S1535-6108(03)00141-7.

    Article  CAS  PubMed  Google Scholar 

  41. Coppe JP, Itahana Y, Moore DH, Bennington JL, Desprez PY: Id-1 and Id-2 proteins as molecular markers for human prostate cancer progression. Clin Cancer Res. 2004, 10 (6): 2044-2051. 10.1158/1078-0432.CCR-03-0933.

    Article  CAS  PubMed  Google Scholar 

  42. Asirvatham AJ, Schmidt MA, Chaudhary J: Non-redundant inhibitor of differentiation (Id) gene expression and function in human prostate epithelial cells. Prostate. 2006, 66 (9): 921-935. 10.1002/pros.20366.

    Article  CAS  PubMed  Google Scholar 

  43. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5 (10): R80-10.1186/gb-2004-5-10-r80.

    Article  PubMed Central  PubMed  Google Scholar 

  44. Ploner A, Miller LD, Hall P, Bergh J, Pawitan Y: Correlation test to assess low-level processing of high-density oligonucleotide microarray data. BMC Bioinformatics. 2005, 6: 80-10.1186/1471-2105-6-80.

    Article  PubMed Central  PubMed  Google Scholar 

Download references

Acknowledgements

We wish to thank Kang Liu from the Burnham Institute for Medical Research Microarray Facility for BeadChip hybridizations and Roy Williams for help with gene analysis. This work was supported by Department of Defense Grant W81XWH-04-1-0515, and National Institutes of Heath Grant R01 CA108959 (to J.W.S.).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Jeffrey W Smith.

Additional information

Authors' contributions

JS and LK conceived and designed the study. LK carried out the siRNA transfections, validated target knockdown, and prepared RNA for BeadArray expression profiling, performed data analysis and drafted the manuscript. JS supervised the coordination of the study and participated in manuscript preparation. All authors read and approved the final manuscript.

Electronic supplementary material

12864_2006_881_MOESM1_ESM.xls

Additional file 1: Genes significantly modified by FAS siRNA duplex #1 (≥1.2 or ≤0.83 fold). The fold change values represent the average from 2 replicates of FAS siRNA duplex #1 compared to non-silencing control siRNA. Significance was determined using a p-value of < 0.05, with a 1.2 fold change cutoff for both biological replicates at a given time point. (XLS 328 KB)

12864_2006_881_MOESM2_ESM.xls

Additional file 2: Genes significantly modified by FAS siRNA duplex #2 (≥1.2 or ≤0.83 fold). The fold change values represent the average from 2 replicates of FAS siRNA duplex #2 compared to non-silencing control siRNA. Significance was determined using a p-value of < 0.05, with a 1.2 fold change cutoff for both biological replicates at a given time point. (XLS 372 KB)

12864_2006_881_MOESM3_ESM.xls

Additional file 3: Genes significantly modified by FAS siRNA duplex #3 (≥1.2 or ≤0.83 fold). The fold change values represent the average from 2 replicates of FAS siRNA duplex #3 compared to non-silencing control siRNA. Significance was determined using a p-value of < 0.05, with a 1.2 fold change cutoff for both biological replicates at a given time point. (XLS 404 KB)

12864_2006_881_MOESM4_ESM.xls

Additional file 4: Genes significantly modified by FAS siRNA duplex #4 (≥1.2 or ≤0.83 fold). The fold change values represent the average from 2 replicates of FAS siRNA duplex #4 compared to non-silencing control siRNA. Significance was determined using a p-value of < 0.05, with a 1.2 fold change cutoff for both biological replicates at a given time point. (XLS 382 KB)

12864_2006_881_MOESM5_ESM.xls

Additional file 5: FAS knockdown signature of 279 genes modified by 1.5 fold. The table displays genes modified by 1.5 fold compared to non-silencing control siRNA. The fold change values represent the average from at least 3 of the FAS siRNA duplexes. (XLS 104 KB)

12864_2006_881_MOESM6_ESM.xls

Additional file 6: GSEA pathways down-regulated in response to knockdown of FAS. The table displays all pathways found to be down-regulated in response to FAS siRNA treatment compared to non-silencing control siRNA. Significance was determined using a nominal p-value < 0.05 or FDR < 0.250. (XLS 17 KB)

12864_2006_881_MOESM7_ESM.xls

Additional file 7: GSEA pathways up-regulated in response to knockdown of FAS. The table displays all pathways found to be up-regulated in response to FAS siRNA treatment compared to non-silencing control siRNA. Significance was determined using a nominal p-value < 0.05 or FDR < 0.250. (XLS 24 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Knowles, L.M., Smith, J.W. Genome-wide changes accompanying knockdown of fatty acid synthase in breast cancer. BMC Genomics 8, 168 (2007). https://doi.org/10.1186/1471-2164-8-168

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-8-168

Keywords