Loss of HSulf-1 promotes altered lipid metabolism in ovarian cancer
Cancer & Metabolism volume 2, Article number: 13 (2014)
Loss of the endosulfatase HSulf-1 is common in ovarian cancer, upregulates heparin binding growth factor signaling and potentiates tumorigenesis and angiogenesis. However, metabolic differences between isogenic cells with and without HSulf-1 have not been characterized upon HSulf-1 suppression in vitro. Since growth factor signaling is closely tied to metabolic alterations, we determined the extent to which HSulf-1 loss affects cancer cell metabolism.
Ingenuity pathway analysis of gene expression in HSulf-1 shRNA-silenced cells (Sh1 and Sh2 cells) compared to non-targeted control shRNA cells (NTC cells) and subsequent Kyoto Encyclopedia of Genes and Genomics (KEGG) database analysis showed altered metabolic pathways with changes in the lipid metabolism as one of the major pathways altered inSh1 and 2 cells. Untargeted global metabolomic profiling in these isogenic cell lines identified approximately 338 metabolites using GC/MS and LC/MS/MS platforms. Knockdown of HSulf-1 in OV202 cells induced significant changes in 156 metabolites associated with several metabolic pathways including amino acid, lipids, and nucleotides. Loss of HSulf-1 promoted overall fatty acid synthesis leading to enhance the metabolite levels of long chain, branched, and essential fatty acids along with sphingolipids. Furthermore, HSulf-1 loss induced the expression of lipogenic genes including FASN, SREBF1, PPARγ, and PLA2G3 stimulated lipid droplet accumulation. Conversely, re-expression of HSulf-1 in Sh1 cells reduced the lipid droplet formation. Additionally, HSulf-1 also enhanced CPT1A and fatty acid oxidation and augmented the protein expression of key lipolytic enzymes such as MAGL, DAGLA, HSL, and ASCL1. Overall, these findings suggest that loss of HSulf-1 by concomitantly enhancing fatty acid synthesis and oxidation confers a lipogenic phenotype leading to the metabolic alterations associated with the progression of ovarian cancer.
Taken together, these findings demonstrate that loss of HSulf-1 potentially contributes to the metabolic alterations associated with the progression of ovarian pathogenesis, specifically impacting the lipogenic phenotype of ovarian cancer cells that can be therapeutically targeted.
Metabolic re-programming has recently emerged as a new hallmark of cancer. Alteration of cellular metabolism in cancer cells is proposed to increase the availability of essential building blocks that support uncontrolled cellular proliferation. Most cancer cells, although diversified by etiology and type, reprogram their metabolism to accumulate metabolic intermediates as sources of building blocks. The Warburg effect is one of the most important metabolic alteration in cancer, in which neoplastic cells exhibit higher glucose uptake and utilization by altering glucose metabolism even in the presence of oxygen to produce lactate from glucose and thereby decouples glycolysis from mitochondrial oxidation. Although alterations in fatty acids (FAs) and lipid metabolism have received less attention, recently, their importance in cancer metabolism is being increasingly recognized. The total lipid pools required for membrane synthesis of dividing cells are derived mainly from FAs and in part from acetyl CoA. Cancer cells meet their FAs demand mainly by increasing de novo FA synthesis rather than from exogenous sources and partly from acetyl CoA[5, 6]. The activated FAs can then be utilized in the synthesis of membrane phospholipids (PLs) including phosphatidylcholine (PC), phosphatidylethanolamine (PE) in addition to sterols, sphingolipids, and lysolipids to meet the energy demands and proliferation. The rest of the activated FAs then function as signaling molecules or esterified as triglycerols or sterol esters and stored in lipid droplets (LDs)[7, 8].
Accumulating evidence suggest that activation of oncogenes such as MYC, NF-kB, K-RAS[2, 9] and loss of tumor suppressor genes (P53, LKB1/AMPK)[10–12] forms a basis for altered metabolism of cancer cells. The oncogenic activation of c-MYC turns on the Ras-Raf-MAPK signaling pathway along with HIF1α and PI3K-Akt-mTOR axis which transcriptionally stimulate the expression of most glycolytic and glutaminolytic genes and subsequently activate lipid metabolism. Similarly, loss of tumor suppressors forms the basis of the Warburg effect leading to carcinogenesis. For example, p53 activates TIGAR to reduce the cellular accumulation of fructose 2,6 bisphosphate, an allosteric activator of phosphofructo kinase, a critical control point in glycolysis leads to decreased glycolysis. Moreover, growth factor-mediated phosphorylation of pyruvate kinase isoform 2 (PKM2) and mutations of genes in IDH1/2 in the metabolic pathways such as glycolysis and the TCA cycle respectively have been identified that support growth of the rapidly proliferating cells and survive metabolic stress[16, 17].
Human Sulfatase 1 (HSulf-1), an endosulfatase established as a putative tumor suppressor in ovarian cancer, has been shown to modulate the signaling of growth factors and cytokines in tumor microenvironment. Our previous work demonstrated that loss of HSulf-1 modulates heparin-binding growth factors such as bFGF, VEGF, HGF, PDGF, and heparin binding EGF (HB-EGF) signaling, which plays an important role in tumor progression, metastasis, and angiogenesis[19–22]. Moreover, serous tumors with moderate to high levels of HSulf-1 had better prognosis in terms of overall survival, implicating its critical role in the progression of ovarian cancer. Our more recent data demonstrated that HSulf-1 knockdown clones in the OV202 ovarian cancer cell line (OV202Sh1 and Sh2 cells) have significantly increased the ability to form anchorage-independent colonies in vitro and enhanced tumorigenicity in vivo. Also, in breast cancer, HSulf-1 is negatively regulated by HIF1α, but positively by von Hippel-Lindau tumor suppressor gene[23, 25, 26]. These findings led us to hypothesize that HSulf-1 might play a unique role to alter tumor microenvironment raising the possibility that its loss might alter the cellular metabolism and levels of the resulting metabolites as a downstream effect of altered growth factor signaling.
In the current study, global changes in metabolism were investigated in HSulf-1 silenced OV202 cells by microarray, metabolic data analysis, and Western blotting. Here, we for the first time report that loss of HSulf-1 promotes overall fatty acid synthesis and oxidation leading to a lipogenic phenotype to promote cancer growth in ovarian cancer.
OV202 cell line was low-passage primary line established at the Mayo Clinic. OV202NTC, Sh1, Sh2, and Cl 11 cells were cultured in 5% CO2-95% air humidified atmosphere at 37°C with minimal essential medium supplemented with 20% fetal bovine serum and 1 μg/ml puromycin, with nonessential amino acids. All cell lines were tested using a PCR-based assay and found to be free of Mycoplasma contamination.
HSulf-1 short-hairpin (sh) RNA1 (Sh1- AGCTACCCTGGGTTCCTTTGT) which targets the 3′-untranslated region (UTR) was cloned into lentiviral vector pLKO.1-puro as described previously. HSulf-1 shRNA2 (Sh2- CGTCGAATTTGAAGGTGAAAT) and nontargeted control shRNAs (NTC shRNA- ACTTACGAGTGACAGTAGATT) cloned into the lentivirus vector pLKO.1-puro were chosen from the human library (MISSION TRC-Hs 1.0) and purchased as glycerol stock from Sigma. Transfection with Fugene (Roche) was performed according to the manufacturer's instructions. Transducted cells were selected with 1 μg/ml puromycin.
Rescue of HSulf-1 in Sh1 cells
pcDNA-HSulf-1 plasmid was cloned as described previously. Since Sh1 shRNA targeted the 3′UTR of HSulf-1, we rescued the expression of HSulf-1 in this cell line with CMV-driven WT expression construct and selected stable clone C11 as previously described. Vector only transfected cells served as controls.
Microarray expression data analysis
OV202 NTC, Sh1, and Sh2 cells in triplicates were profiled using Illumina Human HT-12 3.0 Expression Beadchip array as previously described. Microarray expression data were analyzed on the log2 scale. Data quality was assessed via box and whisker plots along with residual and pair-wise MVA plots before and after normalization[29, 30]. All arrays were normalized together using fastlo, a non-linear normalization similar to cyclic loess which runs in a fraction of the time. Both supervised and unsupervised analyses were performed. Supervised analysis to determine differentially expressed genes was performed using Significance Analysis of Microarrays (SAM). For SAM analysis for unpaired samples, Biometric Research Branch (BRB)-ArrayTools (Version 3.7.0, developed by Dr. Richard Simon and Amy Peng Lam.) was used with Delta set to 0.822, resulting in false discovery rate <5%. Unsupervised clustering was performed using the one minus correlation metric with average linkage. Heat maps were generated for visualization. Pathway analysis was performed using Ingenuity Pathway Analysis (Ingenuity® Systems,http://www.ingenuity.com).
Liquid chromatography/mass spectrometry (LC/MS, LC/MS2)
The LC/MS portion of the platform was based on a Waters ACQUITY UPLC and a Thermo-Finnigan LTQ mass spectrometer, which consisted of an electrospray ionization (ESI) source and linear ion-trap (LIT) mass analyzer. The sample extract was split into 2 aliquots, dried, then reconstituted in acidic or basic LC-compatible solvents, each of which contained 11 or more injection standards at fixed concentrations. One aliquot was analyzed using acidic positive ion optimized conditions and the other using basic negative ion optimized conditions in two independent injections using separate dedicated columns. Extracts reconstituted in acidic conditions were gradient eluted using water and methanol both containing 0.1% Formic acid, while the basic extracts, which also used water/methanol, contained 6.5 mM ammonium bicarbonate. The MS analysis alternated between MS and data-dependent MS2 scans using dynamic exclusion.
Gas chromatography/mass spectrometry (GC/MS)
The samples destined for GC/MS analysis were re-dried under vacuum desiccation for a minimum of 24 h prior to being derivatized under dried nitrogen using bistrimethyl-silyl-triflouroacetamide (BSTFA). The GC column was 5% phenyldimethyl silicone and the temperature ramp is from 40°C to 300°C in a 16-min period. Samples were analyzed by a Thermo-Finnigan Trace DSQ fast-scanning single-quadrupole mass spectrometer using electron impact ionization. The instrument was tuned and calibrated for mass resolution and mass accuracy on a daily basis. The information output from the raw data files was automatically extracted as discussed below.
Data extraction and compound identification
Peaks were identified using Metabolon's proprietary peak integration software. Compounds were identified by comparison to library entries of purified standards or recurrent unknown entities. Identification of known chemical entities was based on comparison to metabolomic library entries of purified standards. As of this writing, more than 2,600 commercially available purified standard compounds had been identified and registered into LIMS for distribution to both the LC and GC platforms for determination of their analytical characteristics. The combination of chromatographic properties and mass spectra gave an indication of a match to the specific compound or an isobaric entity. Metabolon data analysts use proprietary visualization and interpretation software to confirm the consistency of peak identification among the various samples. Library matches for each compound were checked for each sample and corrected if necessary.
Raw data from each sample was normalized to protein concentration as measured by Bradford assay prior to statistical analysis.
Each sample received was accessioned into the Metabolon LIMS system and was assigned by the LIMS, a unique identifier, which was associated with the original source identifier only. This identifier was used to track all sample handling, tasks, and results. The samples (and all derived aliquots) were bar-coded and tracked by the LIMS system. All portions of any sample were automatically assigned their own unique identifiers by the LIMS when a new task was created; the relationship of these samples was also tracked. All samples were maintained at -80°C until processed.
The sample preparation process was carried out using the automated MicroLab STAR® system from Hamilton Company (Reno, NV, USA). Recovery standards were added prior to the first step in the extraction process for QC purposes. Sample preparation was conducted using a proprietary series of organic and aqueous extractions to remove the protein fraction while allowing maximum recovery of small molecules. The resulting extract was divided into four fractions; two for analysis by LC and one for analysis by GC and a forth as a spare. Samples were placed briefly on a TurboVap® (Zymark, Hopkinton, MA, USA) to remove the organic solvent. Each sample was then frozen and dried under vacuum. Samples were then prepared for the appropriate instrument, either LC/MS or GC/MS.
For QA/QC purposes, a number of additional samples are included with each day's analysis. Furthermore, a selection of QC compounds is added to every sample, including those under test. These compounds are carefully chosen so as not to interfere with the measurement of the endogenous compounds. These QC samples are primarily used to evaluate the process control for each study as well as aiding in the data curation.
Differently expressed metabolites between Sh1/Sh2 and baseline conditions were firstly mapped to KEGG metabolites IDs according to Human Metabolome Database (HMDB; URL:http://www.hmdb.ca/). Then, pathway analysis was performed to highlight relevant metabolic pathways defined in KEGG database (http://www.genome.jp/kegg/), using an on-line tool named MetaboAnalyst (http://www.metaboanalyst.ca/)[34, 35]. Specifically, two types of pathway analysis were done: one is over-representation analysis using hypergeometric test, asking if differentially expressed metabolites are particularly enriched in a same pathway; the other is pathway topology analysis summarizing relative-betweeness centrality, investigating potential pathway impact of observed metabolite changes based on known pathway topology relationships.
Western blot analysis
Western blot analysis was performed as described previously. Whole cell lysates were analyzed with the following antibodies: FASN, ASCL1 (Cell signaling), SREBP1c, PLA2G3, HSulf-1 (Abcam, AB96533), CPT1A, HSL, DAGLA, β-tubulin (GeneTex) and β-actin (Sigma-Aldrich).
Quantitative real-time PCR (qRT-PCR) was carried out using SYBR-Green PCR Master Mix (Applied Biosystems, Foster City, CA, USA), with specific primers for the genes shown in this study. GAPDH or 18S ribosomal subunit (Applied Biosystems) were used as internal control in a Light Cycler kit (BioRad Chromo 4). Normalization across samples was carried out using the average of the constitutive human gene 18S and/or GAPDH primers and calculated as previously described. Binding efficiencies of primer sets for both target and reference genes were similar.
Cells (50,000) were seeded on a coverslip in a 24-well plate and were grown for 24 hours in the presence of complete growth medium. Cells were washed and fixed in 4% paraformaldehyde for 10 min at room temperature before staining with 1 μg/ml BODIPY (493/503; Sigma, St. Louis, MI, USA) in PBS for 10 min at room temperature. Coverslips were washed with PBS and mounted in a slide with Prolong Gold Antifade Reagent (Invitrogen). BODIPY stained cells were examined under inverted confocal fluorescence microscope (Zeiss).
To determine the effect of PLA2G3 on lipid droplet biogenesis OV202 NTC cells were transiently transfected with plasmids containing empty vector or cDNA encoding PLA2G3. After 24 h of transfection, we performed BODIPY staining to visualize lipid droplets. PLA2G3 plasmid was obtained on a MTA from Addgene.
Fatty acid synthesis
Cells were washed twice in cold PBS and resuspend in lysis buffer (50 mM Tris-HCl, pH 7.4, 1 mM EDTA, 150 mM NaCl and PMSF. Cells were sonicated and homogenized by dounce homogenizer followed by centrifugation at 13,000 rpm for 15 min at 4°C. The supernatant was collected and proteins were measured by Bradford assay, and 100 μg of protein was used to conduct FASN activity assay. FASN activity was measured by protocol described by Vazquez-Martin et al.. Briefly, 100 μg of protein was incubated with 240 μM NADPH, 30 μM acetyl CoA and 50 μM malonyl CoA in assay buffer (200 mM potassium phosphate, pH 6.6, 1 mM DTT, 1 mM EDTA) and oxidation of NADPH was measured by monitoring the absorbance at 340 nM over the period of the time. Results presented here compares FASN activity monitored for 10 min. Values are presented as nanomolar NADPH oxidized per minute per milligram of protein.
Fatty acid oxidation
Oxygen consumption rate was measured using a Seahorse Bioscience XF24 flux analyzer. 5 × 104 cells were seeded per well in triplicates in MEM-α containing 20% FBS in an XF24 well culture microplates and incubated overnight in a 37°C/10% CO2 incubator. The assay medium for FAO is low-buffered KHB buffer (110 mM NaCl, 4.7 mM KCl, 2 mM MgSO4, 1.2 mM Na2HPO4, 2.5 mM glucose adjusted to pH 7.4) supplemented with 0.5 mM carnitine. For induction of FAO, BSA conjugated palmitate was injected to a final concentration of 50 μM. XF analyses were performed in the XF Extracellular Flux Analyzer (Seahorse Bioscience, Billerica, MA, USA). Three basal rates were measured prior to automated injection of palmitic acid (50 μM) coupled to BSA vehicle or BSA vehicle alone. After treatment for 55 min, the carnitine palmitoyl transferase-1 inhibitor, Etomoxir (ETO, 50 μM), was added. Oxygen consumption rates were measured by using time-resolved method (Seahorse Bioscience XF24) (21). Data were normalized to protein content (assayed after completion of measurements).
Equal number of cells (1 × 105) was plated in triplicate in 12-well plates. OV202NTC, Sh1, and Sh2 cells were counted after 24, 48, and 72 h using a cellometer (Nexelom, Lawrence, MA, USA). For Etomoxir treatment, Equal number of cells (1 × 105) were seeded in 12-well plates in triplicate and treated with increasing concentration of Etomoxir (0 to 100 μM) for 24 h and total cell numbers were counted using cellometer.
Loss of HSulf-1 comprehensively altered major metabolic pathways
We recently reported that HSulf-1 knockdown clones in the OV202 ovarian cancer cell line (OV202Sh1 and Sh2 cells) have significantly increased ability to form anchorage-independent colonies in vitro and enhanced tumorigenicity in vivo. Consistent with these observations, our growth assays showed enhanced growth rate in Sh1 and Sh2 cells compared to NTC cells (Additional file1: Figure S1). Here, to elucidate the function of putative tumor suppressor HSulf-1 in the metabolism of ovarian cancer, we performed gene expression profiling of stably knockdown HSulf-1 clonal lines OV202 Sh1 and Sh2 cells (referred to from hereon as Sh1 and Sh2) compared to HSulf-1 expressing non-targeted control cells (OV202NTC, referred to as NTC) in triplicates using Illumina HumanHT 12 v3 platform. Unbiased hierarchical clustering and heat maps showed that genes in several different pathways were differentially expressed in Sh1 and Sh2 compared to NTC cells (Figure 1A). We found that over 1,645 and 780 genes were differentially expressed in Sh1 and Sh2 cells, respectively, compared to NTC at 2.6 FC (p and FDR <0.0001). We identified 500 and 280 altered genes in Sh1 and Sh2 cells respectively by significance analysis of microarrays (SAM) from the comprehensive list of 2,752 genes which encoded all known human metabolic enzymes and transporters reported by Possemato et al.. Ingenuity pathway analysis (http://www.ingenuity.com) for these genes showed that most genes were differentially regulated in the fatty acid/lipid pathways in Sh1 and Sh2 cells compared to NTC cells (Additional file2: Table S1, Figure 1B). We next explored the 271 genes in the lipid related pathways from our microarray data by unsupervised clustering (Figure 1C) and found that 26% (73 of 271) genes were differentially expressed in Sh1 and Sh2 cells compared to NTC cells (Figure 1D), indicating that tumor suppressor HSulf-1 possibly regulates the lipid metabolism in ovarian cancer cells.
Loss of HSulf-1 altered global metabolic profile in ovarian cancer cells
To determine whether loss of HSulf-1 has effect in cellular metabolism, we performed unbiased global metabolic profiling using the Metabolon platform (Metabolon Inc, Durham, NC, USA) in Sh1 and Sh2 cells compared to NTC. The samples were extracted using Metabolon’s standard solvent extraction method from cells in logarithmic phase with 5 biological replicates for each sample and distributed into equal parts for analysis on the GC/MS and LC/MS/MS platforms. Hierarchical clustering revealed that metabolite levels of Sh1 and Sh2 cells cluster together separately from control NTC (Figure 2A). Our initial principal component analysis (PCA) revealed that Sh1 and Sh2 cells had specific group of metabolites which were different from NTC cells (Figure 2B). Additionally, t test in each principal component dimension showed that PCA differences between Sh1-Sh2 and NTC cells mainly reside in the first and second principal components, with statistically significant p values of 1.8e - 10 and 4.5e - 6, respectively. Consistent with this, two-dimensional PCA plots showed differences between Sh1-Sh2 and NTC are mainly in the first and second PC dimensions (Figure 2C,D,E). We also performed linear regression analysis to evaluate association of each principal component vs. major metabolite classes, to determine if lipid class is one of the dominant factors determining the first and second principal components. The association was evaluated in a multivariate regression model (Additional file3: Figure S2), where coefficients of first/second PCs were treated as dependent variable, and metabolite class labels were regarded as independent variables. This analysis showed that the lipid class has statistically significant association with both first and second PC dimensions, with p values = 5.0e - 3 and 1.4e - 3, respectively. Interestingly, the peptide class, despite its relatively small size, also has significant associations with both first and second PC dimensions, while the amino acid class has significant association with the first principal component.
Moreover, this analysis showed a total of 338 known metabolites altered by the loss of HSulf-1 in Sh1 and Sh2 cells. Among them, Sh1 and Sh2 cells had a total of 193 and 188 biochemical, respectively, which were significantly altered (p < 0.05, Welch's t test) compared to NTC cells (Figure 2F, Additional file4: Table S2, and Additional file5: Table S3). Additionally, with the loss of HSulf-1, 129 and 114 metabolites were downregulated, whereas 64 and 74 metabolites were upregulated in Sh1 and Sh2 cells, respectively. Also, 40% of these biochemicals were altered in the same direction in Sh1 and Sh2 cells compared to NTC cells (p = <0.05, with false discovery rate (FDR) <0.05, Additional file6: Table S4). Interestingly, we found that among all the metabolites altered, 28% were lipids, 22% were amino acids, whereas sugar and dipeptide comprised 7% (Figure 2G), implicating that loss of HSulf-1 mediates a global metabolic alteration in ovarian cancer with changes in the lipid class being a major contributor.
However, some metabolites were differentially altered in the Sh1 and Sh2 cells. The differentially changed metabolites between these ShRNAs could also be a function of the level of knockdown and where the ShRNAs could have integrated. The extent of knockdown in Sh1 targeting the 3'UTR was close to 100%. However, with Sh2 RNA targeting the open reading frame, there was still some level of HSulf-1 present (Figure 1A,).
Alteration of lipid metabolites upon HSulf-1 loss
The lipogenic phenotype characterized by the activation of lipid metabolism is recognized as a universal feature of most cancers[40, 41]. Apart from the fatty acid (FA) uptake, cancer cells requires de novo FA biosynthesis to synthesize new membranes, to store energy in lipid droplets and to form the lipidic platform for signaling in membrane level in lipid rafts for increased signaling of cell growth receptors[42, 43]. Moreover, circulating lipids also play a significant role in cancer cell growth, migration, and invasion[44, 45]. Of the 156 common metabolites altered in Sh1 and Sh2 cells, 44 (28%) metabolites belonged to the lipid class including long-chain FAs, lysolipids, sphingolipids, glycerolipids, eicosanoids, and carnitine (Additional file6: Table S4), in which 20 metabolites were upregulated (45%). Using the Kyoto Encyclopedia of Genes and Genomics (KEGG) database, we mapped these metabolites to major pathways impacted with alterations in the key junctions[34, 35]. These included linoleic acid, glycerophospholipid, arachidonic acid (Additional file7: Figure S3), and sphingolipid pathways[46, 47]. Consistent with these data, we found that among the 20 long-chain fatty acid metabolites detected, 13 metabolites were significantly increased by 2–6-fold including palmitate, stearate, and oleate, while the remaining 7 long-chain FAs did not show any significant differences in both Sh1 and Sh2 cells (Figure 3A). Of note, docosadienoate (22:2n6), 10-nonadeconoate and eicosenoate (20:1n9) levels were increased approximately sevenfold compared to NTC. Additionally, all the detected branched fatty acid metabolites were augmented (Figure 3B) while five of six essential fatty acids detected were increased 2–5.5-fold (Figure 3C).
The key structural lipids in cell membranes are the glycerol-phospholipids including phosphatidyl-choline, phosphatidyl-ethanolamine, phosphatidyl-serine, phosphatidyl-inositol, and phosphatidic acid in addition to other lipids, such as sterols, sphingolipids, and lyso-phospholipids. Sphingolipids such as ceramide, sphingosine, and sphingosine-1-phosphate are bioactive lipids which can dictate the signaling including growth factor responses, inflammation, apoptosis, and proliferation. Our metabolite analysis revealed that HSulf-1 deficiency was closely associated with increased levels of sphinganine (FC = 21.79, p = <0.001), sphingosine (FC = 5.08, p = <0.001), palmitoyl sphingomyelin (FC = 1.62, p = <0.0047, and FC = 1.52, p = <0.001), and stearoyl sphingomyelin (FC = 6.24, p = <0.001 and FC = 5.30, p = <0.001) (Figure 3D). Additionally, we found that knockdown of HSulf-1 expression also largely affected the metabolite levels of choline/inositol pathway (Figure 3E) and lysophospholipids (Figure 3F) including reduced glycerophosphorylcholine (GPC) (FC = -1.5, p = 0.0940 in Sh1 cells only) and glycerol 3-phosphate (G3P) (FC = -1.6, p = 0.0154 and FC = -1.9, p = <0.001). Altogether, these results implicated a key role of HSulf-1 in increased lipid metabolism and signaling.
Deficiency of HSulf-1 in ovarian cancer induced higher expression of ‘lipogenic genes’
All the lipid molecules in cells are derived in part from acetyl CoA, and many contain FAs. These FA building blocks come from either exogenous sources or from de novo FA synthesis. Thus, malignant cells synthesize their own FA de novo and thereby exhibit a preference over exogenous FA uptake, while most normal human cells prefer exogenous sources. Our microarray analysis showed that mRNA levels of lipogenic enzymes fatty acid synthase (FASN), sphingosine kinase 1 (SPHK1), PLA2G4A, PLA2G3, sterol regulatory element-binding transcription factor 1 (SREBF1), and peroxisome proliferator-activated receptor (PPARγ) are also upregulated upon loss of HSulf-1 (Figure 4A). Consistent with the increased levels of these genes and lipid metabolites, qRT PCR and immunoblotting showed enhanced mRNA (Figure 4B) and protein (Figure 4C) expressions of FASN, SPHK1, PLA2G4A, PLA2G3, SREBF1, and PPARγ in the Sh1 and Sh2 cells compared to the NTC cells. The increase in the mRNA and protein expression of FASN along with the increased production of several long chain FAs in Sh1 and Sh2 cells indicated that fatty acid synthesis was enhanced in HSulf-1 silenced cells. To confirm, we also measured the enzymatic activity of FASN and found almost twofold higher activity of FASN in Sh1 and Sh2 cells compared to NTC cells (Figure 4D). These data suggest that loss of HSulf-1 can increase the activity of FASN to enhance FA synthesis.
Interestingly, FASN, SPHK1, PLA2G4A, PLA2G3, SREBF1, and PPARγ are the key enzymes which are localized and involved in the biogenesis of lipid droplets (LD)[50–54]. The FAs are activated by covalent modification by CoA via fatty-acyl-CoA synthetases and esterified to glycerol generating triglycerides or sterol esters that are stored in lipid droplets. Accumulation of cytoplasmic LDs forms a basis of increased growth and chemoresistance in neoplastic cells[40, 55]. To explore whether increased expression of above enzymes resulted in LD formation, we imaged the LDs by bodipy staining which clearly showed increased LDs in Sh1 and Sh2 cells. As expected, there were very few LDs in the NTC cells (Figure 4E). To confirm whether HSulf-1 was directly involved in LD biogenesis, we also rescued the expression of HSulf-1 in Sh1 cells by stable transfection of CMV-driven HSulf-1 expression construct (clone 11) with vector-transfected Sh1 cells served as controls (Figure 4F,G). Bodipy staining revealed that rescue of HSulf-1 significantly reduced the number of lipid droplets in Sh1 (Sh1 clone 11) compared to Sh1 vector control with many LDs (Figure 4G) demonstrating that loss of HSulf-1 promotes LD biogenesis. Additionally, transient transfection of PLA2G3 into OV202NTC cells showed increased LDs in these cells compared to vector-transfected controls (Figure 4H), suggesting a role of PLA2G3 in loss of HSulf-1 mediated LD biogenesis.
Loss of HSulf-1 facilitates enhanced β-oxidation and lipolysis
The analysis of metabolite data additionally revealed the increased levels of acetylcarnitine (FC = 3.74, p = <0.001 and FC = 3.22, p = <0.001), butyrylcarnitine, hexanoylcarnitine and octanoylcarnitine in HSulf-1 silenced cells. Moreover, increased levels of oleolycarnitine were 16.84- and 6.64-fold in Sh1 and Sh2, respectively (Figure 5A). Carnitine is important for shuttling FAs across mitochondrial membranes for oxidation and so the escalation of carnitine level enhances β-oxidation support the increased growth of ovarian cancer. To determine if β-oxidation was enhanced in Sh1 and Sh2 cells, we measured the mRNA and protein level of carnitine palmitoyltransferase 1 (CPT1A) and the results showed a higher level of CPT1A in both in Sh1 and Sh2 cells (Figure 5B). The increase in the expression of CPT1A mRNA along with increased production of several of the long chain FAs in Sh1 and Sh2 cells indicated that these cells may utilize long chain FAs by β-oxidation to generate more ATP to accommodate increase proliferation. Next, to confirm this notion, we determined the fatty acid oxidation (FAO), and results showed that FAO was higher both in Sh1 and Sh2 compared to NTC cells upon addition of palmitate (Figure 5C). Subsequently, we demonstrated that Etomoxir, a specific inhibitor of CPT1A, decreased FAO more significantly in Sh1 and Sh2 cells compared to NTC cells. We also demonstrated that Sh1 and Sh2 cells were more sensitive toward Etomoxir treatment than NTC confirming a major role of FAO in their survival (Additional file8: Figure S4).Since FAs are stored in lipid droplets and released by the action of lipolytic enzymes, we determined the expression of monoacylglycerol lipase (MAGL), diacylglycerol lipase alpha (DAGLA), long-chain acyl-CoA synthetase (ACSL1), and hormone-sensitive lipase (HSL). Immunoblot analysis showed that the expression of these enzymes was upregulated in Sh1 and Sh2 cells compared to NTC cells (Figure 5D), demonstrating that loss of HSulf-1 activates beta oxidation and lipolysis.
Altered metabolism in amino acid and peptide super pathways in cells with loss of HSulf-1
In addition to the altered lipid metabolism, the metabolomic analysis revealed that amino acids and their derivatives were the second most class of metabolites altered in HSulf-1 knockdown cells (Additional file4: Table S2). Most of the amino acids including serine, threonine, aspartate, asparagine, alanine, phenylalanine, tyrosine, tryptophan, arginine, ornithine, proline, and methionine were decreased in the range of -1.36- to -3.33-fold (p < 0.05) in both OV202 Sh1 and Sh2 cells. In contrast, glycine was increased marginally in Sh1 and Sh2 cells (range 1.08- to 1.1-fold; p < 0.001 for each). Interestingly, we found significant upregulation of cysteine (eight to tenfold, p < 0.001) and its derivatives, hypotaurine (FC = 2.38, p = <0.001 and FC = 1.99, p = <0.0244) and taurine (FC = 9.91, p = <0.0032 and FC = 5.38, p = <0.0090) in both Sh1 and Sh2 cells (Additional file4: Table S2). We also observed alterations of N-acetylated amino acids including N-acetyl-alanine, serine, and threonine resulted from the action of N-acetyltransferases on acetyl-CoA and L-amino acids where all the six N-acetylated amino acid tested were downregulated. The amino acid-derived antioxidant, both reduced and oxidized glutathione levels were lower in HSulf-1-deficient cells. Additionally, among different dipeptide molecules, gamma-glutamylglutamate (FC = 3.73, 2.77 in Sh1 and Sh2, p = <.0001) and gamma-glutamylmethionine (FC = 7.31, 6.23 in Sh1 and Sh2 respectively, p = <.0001) levels were increased upon loss of HSulf-1. This enhancement may be particularly relevant as increased expression of gamma-glutamylglutamate is reported to be associated with tumor progression and drug resistance observed in human malignancies. Additionally, a significant increase of pro-hydroxy-proline, (FC = 2.05, p = 0.0057 and FC = 1.92, p = 0.0088) a dipeptide, also demonstrated in Sh1 and Sh2 cells. Pro-hydroxy-proline is a marker of collagen and extracellular matrix degradation (Additional file4: Table S2).
Aberrant cellular metabolism in cancer is now well known and is directly related to tumorigenesis in most of the cancers[57, 58]. Multiple signaling pathways and several molecules are involved in the synthesis and degradation of the lipids and also the activities of lipid metabolizing enzymes are regulated by a complex interplay between metabolic, tumor suppressor, and oncogenic signaling. Though loss of HSulf-1, a putative tumor suppressor gene, was well known to promote tumorigenesis, angiogenesis[22, 23], and invasion in breast[18, 21] and ovarian[23, 24] cancers, the role of HSulf-1 was never elucidated in altered metabolism of ovarian cancer cells so far. Thus, in the present study, we have shown that loss of the putative tumor suppressor, HSulf-1 promotes altered metabolic pathways including lipid, amino acid, and nucleotide. Of the several pathways altered by loss of HSulf-1, the lipid and amino acid pathway-related metabolites accounted for 50% of the total altered metabolites identified. In ovarian cancer, altered lipid metabolism was detected in patients during early and late stages of disease compared to healthy controls[60, 61]. In contrast to lipid pathway-related metabolites, there are very few reports on the alteration in the amino acid levels in OVCa. We found a significant downregulation of most of the amino acids with loss of HSulf-1 consistent with Zhang et al.'s report where they showed presence of lower levels of amino acids in the serum of patients with esophageal adenocarcinomas. This downregulation could be due to differences in uptake from the media, increased utilization, and/or catabolism since increased demand for utilization of amino acids has been reported in other cancers[63, 64]. Other studies have shown an increase in amino acid levels in the serum of patients with colon and breast cancers[65, 66]. While the majority of amino acids were downregulated, we found significant upregulation of cysteine and its derivatives hypotaurine and taurine in Sh1 and Sh2 cells. Significant upregulation in cysteine, taurine, and hypotaurine could be the result of increased levels of cysteine dioxygenase type I (CDO1)[67–69] in Sh1 and Sh2 cells (data not shown). While taurine is reported to be downregulated in cancer, it is also reported to be increased in tumors of the prostate, squamous cell carcinoma and liver metastasis[71, 72]. Additionally, levels of myo-inositol and taurine concentrations both in vivo and in vitro are correlated with cell density of the tumors.
Additionally, we also saw an increase in the dipeptide pro-hydroxy-pro in Sh1 and Sh2 cells. Increase in pro-hydroxy-proline is consistent with the report that major cartilage ECM proteins type II collagen and aggrecan were significantly lower in HSulf-/- chondrocytes suggesting that loss of HSulf-1 may regulate the overall balance of cartilage matrix synthesis and degradation. Consistent with this, significant alterations in several UDP-glycosylation moieties, including UDP-acetylglucosamine/UDP-acetylgalactosamine UDP-glucose, UDP-glucuronate, and UDP-galactose are also supportive of changes in extracellular matrix remodeling with reduced HSulf-1 expression.
The major pathway identified by both ingenuity pathway analysis and metabolic profiling was the alterations in the lipid pathway. Lipids are as important building blocks as carbohydrates to form the basic skeleton of rapidly dividing cells, and therefore, large amounts of FAs are required to accommodate high rates of proliferation in cancer cells. Additionally, the source of FAs may determine the phospholipid composition of membranes. Moreover, it was reported that in high-grade ovarian cancer, long-chain fatty acids were elevated. Our results indicate that with the loss of HSulf-1, there was a significant increase in the long-chain FAs along with branched and essential FAs in ovarian cancer. These high levels of FAs in ovarian cancer with the absence of HSulf-1 indicate a major role of HSulf-1 in FA synthesis. This enhanced production of FAs might meet the structural needs of highly proliferative ovarian cancer cells' requirement, when the tumor suppressor HSulf-1 is lost.
Sphingolipids and lysolipid or lysophospholipids are bioactive lipid molecules which play pivotal role in cancer pathogenesis. The most important sphingolipids are ceramide and sphingosine-1-phosphate and the balance between these two define the cell's fate. In cancer, sphingolipid metabolism is altered and includes changes in the levels of sphingolipids and the enzymes involved in their metabolism. Lysolipids with their various FA side chains are also bioactive lipids which mainly function as growth-stimulating factor and induce cell proliferation, differentiation, and cell migration. Herein, we report an increase in sphingolipids including sphinganine, sphingosine, palmitoylsphingomyelin along with altered choline/inositol and lysolipid content. Interestingly, according to the KEGG database, metabolite mapping and glycerophospholipid, arachidonic, and sphingolipid metabolisms were identified as key junctions. These results clearly indicate that with the loss of HSulf-1, the cell remodeled its own lipid synthesis in a manner to supply both structural and signaling lipids to malignant cells.
In cancer cells, the increased lipid synthesis is due to the higher expression and activity of lipogenic enzymes[78, 79]. Changes in the expression and activity of enzymes involved in lipid metabolism are regulated by metabolic and oncogenic signalling pathways[80, 81]. We identified enhanced expression of FASN, SREBF1, SPHK1, and PPARγ related to lipid biogenesis. Upregulation of these key lipogenic enzymes with the loss of HSulf-1 were further corroborated by their mRNA and protein expressions. Moreover higher enzymatic activity of FASN strongly supports the enhanced biogenesis of FAs with HSulf-1 loss. Additionally, phospholipase-related genes PLA2G3, PLA2G4A, PNPLA2, PLA2G4C, and PNPLA6 involved in lipolysis were also found to be higher in Sh1 and Sh2. Interestingly, PLA2G3 is an important enzyme in lipid metabolism and implicated in LD biogenesis[50, 52]. Recently, it was reported that cancer cells contain increased numbers of lipid droplets compared with normal tissue which are storage sites for triglycerides and cholesterol to be used as energy source. Loss of HSulf-1 induced an increase in PLA2G3 expression that may lead to the accumulation of LD in ovarian cancer. Furthermore, reduction in the number of LDs with the rescue of HSulf-1 expression in Sh1 cells establishes a direct connection to the loss of a putative tumor suppressor HSulf-1 and enhanced LD biogenesis in ovarian cancer.
In accordance with these results, we also demonstrated lower level of carnitine along with higher level of carnitine derivatives suggestive of enhanced transport of FAs through the mitochondrial membrane. More importantly, the augmented CPT1A level in Sh1 and Sh2 cells strongly support the possibility of enhanced β-oxidation. Increased and altered FA synthesis in cancer is a well-accepted phenomenon; however, enhanced expression and activity of lipolytic enzymes in tumor cells is a recent observation. In prostate cancer, malignant cells solely rely on FAO as their energy resource. Here, we report that loss of HSulf-1 also promotes enhanced FAO in Sh1 and Sh2 cells compared to NTC. Consistent with increased lipolysis, our results indicated a significant increase in the expression of important lipolytic enzymes including MAGL, DAGLA, HSL, and ASCL1 in Sh1 and Sh2 cells compared to NTC cells. Higher MAGL expression has been shown to stimulate pro-tumorigenic signals and promotes survival, tumor growth, and migration. In this regard, the relationship between synthesis, storage, and utilization of free fatty acids (FFA) through beta oxidation for energy production is very poorly understood. While the prevailing belief that cells that have increased FA synthesis usually do not undergo beta oxidation, it has also been suggested that newly synthesized FFA are immediately converted into neutral lipids and stored in lipid droplets. The FFAs are released by the action of lipases such as mono-acyl glycerol (MAGL) which is then used for new membrane synthesis, lipid signaling and beta oxidation for energy production for the anabolic reactions. Our data seems to support this hypothesis since the Sh1/ Sh2 cells have both increased the rate of FA synthesis and beta oxidation compared to NTC cells. The FFA may be immediately converted to neutral lipids and stored in lipid droplets to be released under stress conditions.
Increasing evidence indicates that metabolic alterations induced by loss of tumor suppressor genes are common in cancer. These alterations are critical for growth and survival of cancer cells. This is the first comprehensive report on metabolic alterations induced by loss of HSulf-1, a putative tumor suppressor gene in ovarian and breast cancer. Our results clearly indicate that loss of HSulf-1 remodels lipid metabolism in ovarian cancer. It is well accepted that HSulf-1 is a major regulator of growth factor-mediated signaling and altered the tumor microenvironment in ovarian cancer; however, how HSulf-1 directly reconstructs the lipid metabolism in terms of FA synthesis, lipolysis, and enhanced LD formation is yet to be unraveled. A better understanding of the function of HSulf-1, its role in energy metabolism (glycolysis, TCA cycles, PPP) and its regulation on target genes involved in altered cellular metabolism could lead to a better understanding of tumorigenesis as well as development of new targeted therapy.
HSulf-1 is reported to be an important tumor-suppressor gene and its expression is lost in a majority of ovarian tumors. This study demonstrates a significant alteration of cellular metabolism upon loss of HSulf-1 in OV202 cells. Microarray analysis and metabolite profiling is performed to ascertain the impact of HSulf-1 on the overall metabolic changes including alterations in lipid and amino acid pathways. The major finding from our study shows the metabolic reprogramming of cells toward an enhanced lipid metabolism upon absence of HSulf-1. Cells adopt a lipogenic phenotype which is manifested with an excess fatty acid synthesis and an upregulated beta-oxidation. Furthermore, HSulf-1-deficient cells accumulate a huge amount of cytoplasmic lipid droplets to accommodate the excessive fatty acid syntheses. The present findings are supported by the increasing evidences of enhanced lipid metabolism in cancer cells facilitating cell survival, proliferation, and signaling. Our results indicate that, loss of HSulf-1 is enabling the cells to synthesize more lipids to expedite the high proliferation rate and survival.
long-chain acyl-CoA synthetase
carnitine palmitoyltransferase 1
diacylglycerol lipase alpha
false discovery rate
fatty acid oxidation
fatty acid synthase
LC/MS/MS: PCA: principal component analysis
hormone sensitive lipase
Kyoto Encyclopedia of Genes and Genomics
OV202: non-targeted control cells
peroxisome proliferator-activated receptor gamma
significance analysis of microarray
stably knockdown HSulf-1 clonal line 1 in OV202
stably knockdown HSulf-1 clonal line 2 in OV202
sphingosine kinase 1
sterol regulatory element-binding transcription factor 1.
Hanahan D, Weinberg RA: Hallmarks of cancer: the next generation. Cell. 2011, 144 (5): 646-674.
Fan TW, Lane AN, Higashi RM, Farag MA, Gao H, Bousamra M, Miller DM: Altered regulation of metabolic pathways in human lung cancer discerned by (13)C stable isotope-resolved metabolomics (SIRM). Mol Cancer. 2009, 8: 41-
Warburg O: On the origin of cancer cells. Science. 1956, 123 (3191): 309-314.
Menendez JA, Lupu R: Fatty acid synthase and the lipogenic phenotype in cancer pathogenesis. Nat Rev Cancer. 2007, 7 (10): 763-777.
Medes G, Thomas A, Weinhouse S: Metabolism of neoplastic tissue. IV A study of lipid synthesis in neoplastic tissue slices in vitro. Cancer Res. 1953, 13 (1): 27-29.
Ookhtens M, Kannan R, Lyon I, Baker N: Liver and adipose tissue contributions to newly formed fatty acids in an ascites tumor. Am J Physiol. 1984, 247 (1 Pt 2): R146-153.
Tauchi-Sato K, Ozeki S, Houjou T, Taguchi R, Fujimoto T: The surface of lipid droplets is a phospholipid monolayer with a unique Fatty Acid composition. J Biol Chem. 2002, 277 (46): 44507-44512.
Wilfling F, Wang H, Haas JT, Krahmer N, Gould TJ, Uchida A, Cheng JX, Graham M, Christiano R, Frohlich F, Liu X, Buhman KK, Coleman RA, Bewersdorf J, Farese RV, Walther TC: Triacylglycerol synthesis enzymes mediate lipid droplet growth by relocalizing from the ER to lipid droplets. Dev Cell. 2013, 24 (4): 384-399.
Yun J, Rago C, Cheong I, Pagliarini R, Angenendt P, Rajagopalan H, Schmidt K, Willson JK, Markowitz S, Zhou S, Diaz LA, Velculescu VE, Lengauer C, Kinzler KW, Vogelstein B, Papdopoulos N: Glucose deprivation contributes to the development of KRAS pathway mutations in tumor cells. Science. 2009, 325 (5947): 1555-1559.
Bensaad K, Tsuruta A, Selak MA, Vidal MN, Nakano K, Bartrons R, Gottlieb E, Vousden KH: TIGAR, a p53-inducible regulator of glycolysis and apoptosis. Cell. 2006, 126 (1): 107-120.
Matoba S, Kang JG, Patino WD, Wragg A, Boehm M, Gavrilova O, Hurley PJ, Bunz F, Hwang PM: p53 regulates mitochondrial respiration. Science. 2006, 312 (5780): 1650-1653.
Shackelford DB, Shaw RJ: The LKB1-AMPK pathway: metabolism and growth control in tumour suppression. Nature Rev. 2009, 9 (8): 563-575.
Levine AJ, Puzio-Kuter AM: The control of the metabolic switch in cancers by oncogenes and tumor suppressor genes. Science. 2010, 330 (6009): 1340-1344.
Luo W, Semenza GL: Emerging roles of PKM2 in cell metabolism and cancer progression. Trends Endocrinol Metab. 2012, 23 (11): 560-566.
Reitman ZJ, Jin G, Karoly ED, Spasojevic I, Yang J, Kinzler KW, He Y, Bigner DD, Vogelstein B, Yan H: Profiling the effects of isocitrate dehydrogenase 1 and 2 mutations on the cellular metabolome. Proc Natl Acad Sci USA. 2011, 108 (8): 3270-3275.
Jones RG, Thompson CB: Tumor suppressors and cell metabolism: a recipe for cancer growth. Genes Dev. 2009, 23 (5): 537-548.
Vogelstein B, Kinzler KW: Cancer genes and the pathways they control. Nat Med. 2004, 10 (8): 789-799.
Khurana A, Liu P, Mellone P, Lorenzon L, Vincenzi B, Datta K, Yang B, Linhardt RJ, Lingle W, Chien J, Baldi A, Shridhar V: HSulf-1 modulates FGF2- and hypoxia-mediated migration and invasion of breast cancer cells. Cancer Res. 2011, 71 (6): 2152-2161.
Lai J, Chien J, Staub J, Avula R, Greene EL, Matthews TA, Smith DI, Kaufmann SH, Roberts LR, Shridhar V: Loss of HSulf-1 up-regulates heparin-binding growth factor signaling in cancer. J Biol Chem. 2003, 278 (25): 23107-23117.
Lai JP, Chien J, Strome SE, Staub J, Montoya DP, Greene EL, Smith DI, Roberts LR, Shridhar V: HSulf-1 modulates HGF-mediated tumor cell invasion and signaling in head and neck squamous carcinoma. Oncogene. 2004, 23 (7): 1439-1447.
Narita K, Chien J, Mullany SA, Staub J, Qian X, Lingle WL, Shridhar V: Loss of HSulf-1 expression enhances autocrine signaling mediated by amphiregulin in breast cancer. J Biol Chem. 2007, 282 (19): 14413-14420.
Narita K, Staub J, Chien J, Meyer K, Bauer M, Friedl A, Ramakrishnan S, Shridhar V: HSulf-1 inhibits angiogenesis and tumorigenesis in vivo. Cancer Res. 2006, 66 (12): 6025-6032.
Liu P, Khurana A, Rattan R, He X, Kalloger S, Dowdy S, Gilks B, Shridhar V: Regulation of HSulf-1 expression by variant hepatic nuclear factor 1 in ovarian cancer. Cancer Res. 2009, 69 (11): 4843-4850.
He X, Khurana A, Roy D, Kaufmann SH, Shridhar V: Loss of HSulf-1 expression enhances tumorigenicity through inhibiting Bim expression in ovarian cancer. Int J Cancer. 2014, 135 (8): 1783-1789.
He X, Khurana A, Maguire JL, Chien J, Shridhar V: HtrA1 sensitizes ovarian cancer cells to cisplatin-induced cytotoxicity by targeting XIAP for degradation. Int J Cancer. 2012, 130 (5): 1029-1035.
Khurana A, Beleford D, He X, Chien J, Shridhar V: Role of heparan sulfatases in ovarian and breast cancer. Am J Cancer Res. 2013, 3 (1): 34-45.
Staub J, Chien J, Pan Y, Qian X, Narita K, Aletti G, Scheerer M, Roberts LR, Molina J, Shridhar V: Epigenetic silencing of HSulf-1 in ovarian cancer: implications in chemoresistance. Oncogene. 2007, 26 (34): 4969-4978.
Chien J, Fan JB, Bell DA, April C, Klotzle B, Ota T, Lingle WL, Gonzalez Bosquet J, Shridhar V, Hartmann LC: Analysis of gene expression in stage I serous tumors identifies critical pathways altered in ovarian cancer. Gynecol Oncol. 2009, 114 (1): 3-11.
Cunningham JM, Oberg AL, Borralho PM, Kren BT, French AJ, Wang L, Bot BM, Morlan BW, Silverstein KA, Staggs R, Zeng Y, Lamblin AF, Hilker CA, Fan JB, Steer CJ, Thibodeau SN: Evaluation of a new high-dimensional miRNA profiling platform. BMC Med Genomics. 2009, 2: 57-
Eckel JE, Gennings C, Therneau TM, Burgoon LD, Boverhof DR, Zacharewski TR: Normalization of two-channel microarray experiments: a semiparametric approach. Bioinformatics. 2005, 21 (7): 1078-1083.
Ballman KV, Grill DE, Oberg AL, Therneau TM: Faster cyclic loess: normalizing RNA arrays via linear models. Bioinformatics. 2004, 20 (16): 2778-2786.
Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001, 98 (9): 5116-5121.
Wishart DS, Jewison T, Guo AC, Wilson M, Knox C, Liu Y, Djoumbou Y, Mandal R, Aziat F, Dong E, Bouatra S, Sinelnikov I, Arndt D, Xia J, Liu P, Yallou F, Bjorndahl T, Perez-Pinero R, Eisner R, Allen F, Neveu V, Greiner R, Scalbert A: HMDB 3.0--the human metabolome database in 2013. Nucleic Acids Res. 2013, 41 (Database issue): D801-807.
Aittokallio T, Schwikowski B: Graph-based methods for analysing networks in cell biology. Brief Bioinform. 2006, 7 (3): 243-255.
Goeman JJ, Buhlmann P: Analyzing gene expression data in terms of gene sets: methodological issues. Bioinformatics. 2007, 23 (8): 980-987.
He X, Ota T, Liu P, Su C, Chien J, Shridhar V: Downregulation of HtrA1 promotes resistance to anoikis and peritoneal dissemination of ovarian cancer cells. Cancer Res. 2010, 70 (8): 3109-3118.
Vazquez-Martin A, Colomer R, Brunet J, Lupu R, Menendez JA: Overexpression of fatty acid synthase gene activates HER1/HER2 tyrosine kinase receptors in human breast epithelial cells. Cell Prolif. 2008, 41 (1): 59-85.
Possemato R, Marks KM, Shaul YD, Pacold ME, Kim D, Birsoy K, Sethumadhavan S, Woo HK, Jang HG, Jha AK, Chen WW, Barrett FG, Stransky N, Tsun ZY, Cowley GS, Barretina J, Kalaany NY, Hsu PP, Ottina K, Chan AM, Yuan B, Garraway LA, Root DE, Mino-Kenudson M, Brachtel EF, Driggers EM, Sabatini DM: Functional genomics reveal that the serine synthesis pathway is essential in breast cancer. Nature. 2011, 476 (7360): 346-350.
Hervé Abdi ALJW: Principal component analysis. Wiley Int Rev: Comput Stat. 2010, 2 (4): 433-459.
Hopperton KE, Duncan RE, Bazinet RP, Archer MC: Fatty acid synthase plays a role in cancer metabolism beyond providing fatty acids for phospholipid synthesis or sustaining elevations in glycolytic activity. Exp Cell Res. 2014, 320 (2): 302-310.
Swinnen JV, Brusselmans K, Verhoeven G: Increased lipogenesis in cancer cells: new players, novel targets. Curr Opin Clin Nutr Metabol Care. 2006, 9 (4): 358-365.
Currie E, Schulze A, Zechner R, Walther TC, Farese RV: Cellular fatty acid metabolism and cancer. Cell Metab. 2013, 18 (2): 153-161.
Swierczynski J, Sledzinski T: [Metabolic and regulatory function of fatty acid synthase]. Postepy Biochem. 2012, 58 (2): 175-185.
Zhang F, Du G: Dysregulated lipid metabolism in cancer. World J Biol Chem. 2012, 3 (8): 167-174.
Zaidi N, Lupien L, Kuemmerle NB, Kinlaw WB, Swinnen JV, Smans K: Lipogenesis and lipolysis: the pathways exploited by the cancer cells to acquire fatty acids. Prog Lipid Res. 2013, 52 (4): 585-589.
Xia J, Mandal R, Sinelnikov IV, Broadhurst D, Wishart DS: MetaboAnalyst 2.0--a comprehensive server for metabolomic data analysis. Nucleic Acids Res. 2012, 40 (Web Server issue): W127-133.
Xia J, Psychogios N, Young N, Wishart DS: MetaboAnalyst: a web server for metabolomic data analysis and interpretation. Nucleic Acids Res. 2009, 37 (Web Server issue): W652-660.
Kraveka YAH JM: Bioactive sphingolipids: an overview on ceramide, ceramide 1-phosphate dihydroceramide, sphingosine, sphingosine 1-phosphate. Handbook of Neurochemistry and Molecular Neurobiology. Edited by: Lajtha A, Tettamanti G, Goracci G. 2010, New York: Springer, 373-383.
Mashima T, Seimiya H, Tsuruo T: De novo fatty-acid synthesis and related pathways as molecular targets for cancer therapy. Br J Cancer. 2009, 100 (9): 1369-1372.
Gubern A, Casas J, Barcelo-Torns M, Barneda D, de la Rosa X, Masgrau R, Picatoste F, Balsinde J, Balboa MA, Claro E: Group IVA phospholipase A2 is necessary for the biogenesis of lipid droplets. J Biol Chem. 2008, 283 (41): 27369-27382.
McDonough PM, Agustin RM, Ingermanson RS, Loy PA, Buehrer BM, Nicoll JB, Prigozhina NL, Mikic I, Price JH: Quantification of lipid droplets and associated proteins in cellular models of obesity via high-content/high-throughput microscopy and automated image analysis. Assay Drug Dev Technol. 2009, 7 (5): 440-460.
Sato H, Kato R, Isogai Y, Saka G, Ohtsuki M, Taketomi Y, Yamamoto K, Tsutsumi K, Yamada J, Masuda S, Ishikawa Y, Ishii T, Kobayashi T, Ikeda K, Taguchi R, Hatakeyama S, Hara S, Kudo I, Itabe H, Murakami M: Analyses of group III secreted phospholipase A2 transgenic mice reveal potential participation of this enzyme in plasma lipoprotein modification, macrophage foam cell formation, and atherosclerosis. J Biol Chem. 2008, 283 (48): 33483-33497.
Schadinger SE, Bucher NL, Schreiber BM, Farmer SR: PPARgamma2 regulates lipogenesis and lipid accumulation in steatotic hepatocytes. Am J Physiol Endocrinol Metab. 2005, 288 (6): E1195-1205.
Wang H, Maechler P, Antinozzi PA, Herrero L, Hagenfeldt-Johansson KA, Bjorklund A, Wollheim CB: The transcription factor SREBP-1c is instrumental in the development of beta-cell dysfunction. J Biol Chem. 2003, 278 (19): 16622-16629.
Cheng YX, Hu M, Chen L, Huang JL, Xia LB, Li BS, Zhou LM, Hong L: The mechanism of lipid raft mediating chemoresistance of cervical cancer. Saudi Med J. 2012, 33 (5): 508-514.
Corti A, Franzini M, Paolicchi A, Pompella A: Gamma-glutamyltransferase of cancer cells at the crossroads of tumor progression, drug resistance and drug targeting. Anticancer Res. 2010, 30 (4): 1169-1181.
Merida I, Avila-Flores A: Tumor metabolism: new opportunities for cancer therapy. Clin Transl Oncol. 2006, 8 (10): 711-716.
Sheng H, Niu B, Sun H: Metabolic targeting of cancers: from molecular mechanisms to therapeutic strategies. Curr Med Chem. 2009, 16 (13): 1561-1587.
Tania M, Khan MA, Song Y: Association of lipid metabolism with ovarian cancer. Curr Oncol. 2010, 17 (5): 6-11.
Xiao Y, Chen Y, Kennedy AW, Belinson J, Xu Y: Evaluation of plasma lysophospholipids for diagnostic significance using electrospray ionization mass spectrometry (ESI-MS) analyses. Ann N Y Acad Sci. 2000, 905: 242-259.
Xu Y, Shen Z, Wiper DW, Wu M, Morton RE, Elson P, Kennedy AW, Belinson J, Markman M, Casey G: Lysophosphatidic acid as a potential biomarker for ovarian and other gynecologic cancers. JAMA. 1998, 280 (8): 719-723.
Zhang J, Bowers J, Liu L, Wei S, Gowda GA, Hammoud Z, Raftery D: Esophageal cancer metabolite biomarkers detected by LC-MS and NMR methods. PLoS One. 2012, 7 (1): e30181-
Wu H, Xue R, Lu C, Deng C, Liu T, Zeng H, Wang Q, Shen X: Metabolomic study for diagnostic model of oesophageal cancer using gas chromatography/mass spectrometry. J Chromatogr B Analyt Technol Biomed Life Sci. 2009, 877 (27): 3111-3117.
Lai HS, Lee JC, Lee PH, Wang ST, Chen WJ: Plasma free amino acid profile in cancer patients. Semin Cancer Biol. 2005, 15 (4): 267-276.
Poschke I, Mao Y, Kiessling R, de Boniface J: Tumor-dependent increase of serum amino acid levels in breast cancer patients has diagnostic potential and correlates with molecular tumor subtypes. J Transl Med. 2013, 11: 290-
Denkert C, Budczies J, Weichert W, Wohlgemuth G, Scholz M, Kind T, Niesporek S, Noske A, Buckendahl A, Dietel M, Fiehn O: Metabolite profiling of human colon carcinoma–deregulation of TCA cycle and amino acid turnover. Mol Cancer. 2008, 7: 72-
Stipanuk MH, Ueki I: Dealing with methionine/homocysteine sulfur: cysteine metabolism to taurine and inorganic sulfur. J Inherit Metab Dis. 2011, 34 (1): 17-32.
Stipanuk MH, Ueki I, Dominy JE, Simmons CR, Hirschberger LL: Cysteine dioxygenase: a robust system for regulation of cellular cysteine levels. Amino Acids. 2009, 37 (1): 55-63.
Ueki I, Roman HB, Valli A, Fieselmann K, Lam J, Peters R, Hirschberger LL, Stipanuk MH: Knockout of the murine cysteine dioxygenase gene results in severe impairment in ability to synthesize taurine and an increased catabolism of cysteine to hydrogen sulfide. Am J Physiol Endocrinol Metab. 2011, 301 (4): E668-684.
Fong MY, McDunn J, Kakar SS: Identification of metabolites in the normal ovary and their transformation in primary and metastatic ovarian cancer. PLoS One. 2011, 6 (5): e19963-
El-Sayed S, Bezabeh T, Odlum O, Patel R, Ahing S, MacDonald K, Somorjai RL, Smith IC: An ex vivo study exploring the diagnostic potential of 1H magnetic resonance spectroscopy in squamous cell carcinoma of the head and neck region. Head Neck. 2002, 24 (8): 766-772.
Moreno A, Lopez LA, Fabra A, Arus C: 1H MRS markers of tumour growth in intrasplenic tumours and liver metastasis induced by injection of HT-29 cells in nude mice spleen. NMR Biomed. 1998, 11 (3): 93-106.
Lehtimaki KK, Valonen PK, Griffin JL, Vaisanen TH, Grohn OH, Kettunen MI, Vepsalainen J, Yla-Herttuala S, Nicholson J, Kauppinen RA: Metabolite changes in BT4C rat gliomas undergoing ganciclovir-thymidine kinase gene therapy-induced programmed cell death as studied by 1H NMR spectroscopy in vivo, ex vivo, and in vitro. J Biol Chem. 2003, 278 (46): 45915-45923.
Otsuki S, Hanson SR, Miyaki S, Grogan SP, Kinoshita M, Asahara H, Wong CH, Lotz MK: Extracellular sulfatases support cartilage homeostasis by regulating BMP and FGF signaling pathways. Proc Natl Acad Sci U S A. 2010, 107 (22): 10202-10207.
Behrouz Z, Justin Bottsford M, Wei H, Rebecca Stone L, Alpa Nick M, Xinna Z, De Koen G, Sood SK: Metabolic distinctions in high grade epithelial ovarian cancer. In Proceedings of the AACR 103rd Annual Meeting, 2012; Chicago, IL. Cancer Res. 2012, 72 (8):
Ryland LK, Fox TE, Liu X, Loughran TP, Kester M: Dysregulation of sphingolipid metabolism in cancer. Cancer Biol Ther. 2011, 11 (2): 138-149.
Sutphen R, Xu Y, Wilbanks GD, Fiorica J, Grendys EC, LaPolla JP, Arango H, Hoffman MS, Martino M, Wakeley K, Griffin D, Blanco RW, Cantor AB, Xiao Y, Krischer JP: Lysophospholipids are potential biomarkers of ovarian cancer. Cancer Epidemiol Biomarkers Prev. 2004, 13 (7): 1185-1191.
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.
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.
Scaglia N, Chisholm JW, Igal RA: Inhibition of stearoylCoA desaturase-1 inactivates acetyl-CoA carboxylase and impairs proliferation in cancer cells: role of AMPK. PLoS One. 2009, 4 (8): e6812-
Yang YA, Han WF, Morin PJ, Chrest FJ, Pizer ES: Activation of fatty acid synthesis during neoplastic transformation: role of mitogen-activated protein kinase and phosphatidylinositol 3-kinase. Exp Cell Res. 2002, 279 (1): 80-90.
Accioly MT, Pacheco P, Maya-Monteiro CM, Carrossini N, Robbs BK, Oliveira SS, Kaufmann C, Morgado-Diaz JA, Bozza PT, Viola JP: Lipid bodies are reservoirs of cyclooxygenase-2 and sites of prostaglandin-E2 synthesis in colon cancer cells. Cancer Res. 2008, 68 (6): 1732-1740.
Zha S, Ferdinandusse S, Hicks JL, Denis S, Dunn TA, Wanders RJ, Luo J, De Marzo AM, Isaacs WB: Peroxisomal branched chain fatty acid beta-oxidation pathway is upregulated in prostate cancer. Prostate. 2005, 63 (4): 316-323.
Nomura DK, Long JZ, Niessen S, Hoover HS, Ng SW, Cravatt BF: Monoacylglycerol lipase regulates a fatty acid network that promotes cancer pathogenesis. Cell. 2010, 140 (1): 49-61.
The work is supported in part by the grants from the National Institutes of Health P50CA136393, CA106954, CA123249, Department of Defense Ovarian Cancer Research Program (W81XWH-13-1-0119, OC120250) and Mayo Clinic CCaTS grant number UL1TR000135 (VS) and CNPQ/Brazil and the Korea Science and Engineering Foundation (KOSEF) grant funded by the Korea government (MEST 2012–0005755) to JCP and SHK, respectively. The authors wish to acknowledge Dr. Andre Terzic (Mayo Clinic) for the use of seahorse equipment.
EDK and LNB are employees of Metabolon, Inc. and, as such, have affiliations with or financial involvement with Metabolon, Inc. These authors have no other relevant affiliations or financial involvement with any organization or entity with a financial interest in or financial conflict with the subject matter or materials discussed in the manuscript apart from those disclosed. The rest of the authors declare that they have no competing interests.
DR technically performed the experiments in Figures 4C,E,G,H and5B,C,D. He also wrote portions of the manuscript and the ‘Methods’ section for FASN and beta oxidation section and was involved in the discussions. SM wrote a major part of the Introduction and Results and discussion, provided figures related to lipid metabolites, western blots in Figure 2E, and performed beta oxidation-related studies in Figure 5C. CW did the bioinformatics analysis of the KEGG pathway of the metabolites, the ingenuity pathway analysis of the microarray data and provided the following figures: Figure 1B,D, Additional file3: Figure S2, Figure 2, Additional file7: Figure S3, and wrote the corresponding section in the MS. XH generated the OV202NTC, Sh1, and Sh2 cells and the rescue clone Sh1 Cl11. AK performed the FASN activity assay and edited the manuscript. SG edited the manuscript and provided input regarding the metaboanalyst analysis. RH and DBJ extracted RNA and provided cDNAs for real-time PCR analysis. SK edited the manuscript and provided financial support for visiting student DBJ. EC and JCP were involved in the discussions and editing of the manuscript. CF trained DR and SM in the use of Seahorse-related experiments for beta oxidation. AM, SD, and JBG edited the manuscript and provided funds to support a technician involved in the study. SR and AO designed the microarray study and normalized the microarray data, provided the write-up related to these methods in the ‘Methods’ section, and reviewed the final manuscript. EDK edited the manuscript. LNB provided the initial write-up on the analysis of the overall metabolomics data and edited the manuscript. JC did the SAM analysis and was involved in the discussion and editing of the manuscript. VS was involved in the design, execution and analysis and did the QPCR of the lipid-related genes and wrote the manuscript with SM and DR. All authors read and approved the final manuscript.
Debarshi Roy, Susmita Mondal, Chen Wang contributed equally to this work.
An erratum to this article is available at http://dx.doi.org/10.1186/2049-3002-2-24.
Electronic supplementary material
Additional file 1: Figure S1: Loss of HSulf-1 induced enhanced proliferation in OV202 cells. Equal number of cells of NTC, Sh1, and Sh2 were plated in triplicates (1 × 105) and counted after 24, 48, and 72 h. The increase in cell count in both Sh1 and Sh2 were statistically significant in 48 h (*p < 0.05) and 72 hr (**p < 0.001). These experiments were repeated twice. (TIFF 91 KB)
Additional file 3: Figure S2: Multivariate regression results between major metabolite classes. The eight major metabolite classes are amino acid, carbohydrate, cofactors and vitamins, energy, lipid, nucleotide, peptide, and xenobiotics. ‘Estimate’ and ‘Std. Error’ are estimate of regression coefficient and estimation of standard-deviation error; ‘t value’ is the t statistics for each coefficient estimate, and ‘Pr(>|t|)’ is the corresponding p value for each coefficient estimate. (TIFF 249 KB)
Additional file 4: Table S2: Changes in major metabolic pathways including amino acids, lipids and nucleotides. List of metabolites identified through mass spectrometry and the super-pathway and sub-pathway for them are shown. The green- and red-shaded boxes are metabolites that are downregulated and/or upregulated in Sh1 and Sh2 cells compared to NTC cells respectively. The p values and the platform used to identify these metabolites are also shown. (XLS 198 KB)
Additional file 5: Table S3: Global metabolic changes upon HSulf-1 loss. List of significantly altered metabolites in Sh1 and Sh2 cells compared to NTC cells in different pathways (XLS). Heat map of statistically significant biochemical profiles in this study. By paired comparisons, shaded cells indicate p ≤ 0.05 (red indicates that the mean values are significantly higher compared untreated control; green values significantly lower). Blue-bolded text indicates 0.05 < p < 0.10. All data were normalized using Bradford protein concentration (red—upregulated, green—downregulated, and blue—approaching significance). (XLS 59 KB)
Additional file 6: Table S4: Alterations in lipid metabolites. Lipid metabolic pathways in Sh1 and Sh2 cells compared to NTC cells. Heat map of statistically significant (p ≤ 0.05, FDR < 0.05) biochemicals are shown. Red indicates that the mean values are significantly higher by comparison with untreated control; green values significantly lower. All data were normalized using Bradford protein concentration (red—upregulated, green—downregulated. (XLSX 10 KB)
Additional file 7: Figure S3: Metaboanalyst pathway analysis. (A) Statistics for pathways with major change based on high impact (linoleic acid metabolism) or p value (pathways glycerophospholipid, arachidonic, and sphingolipid metabolic pathways). (B) Of the 12 highly significant KEGG pathways plotted according to global test p value (intensity of color in the vertical axis) and impact factor (size of the circles in the horizontal axis), all 4 belong to the lipid pathways. (TIFF 272 KB)
Additional file 8: Figure S4: Effect of etomoxir on cellular growth with increasing concentration of Etomoxir treatment (0 to 100 μM) in NTC, Sh1, and Sh2 cells (n = 2). At 60 μM and onwards, the cell growth inhibition was statistically significant (p < 0.05) in both Sh1 and Sh2 cells compared to NTC. (TIFF 100 KB)
About this article
Cite this article
Roy, D., Mondal, S., Wang, C. et al. Loss of HSulf-1 promotes altered lipid metabolism in ovarian cancer. Cancer Metab 2, 13 (2014). https://doi.org/10.1186/2049-3002-2-13