Next Article in Journal
Multiple Myeloma as a Bone Disease? The Tissue Disruption-Induced Cell Stochasticity (TiDiS) Theory
Next Article in Special Issue
Neoadjuvant Systemic Treatment of Primary Angiosarcoma
Previous Article in Journal
Brain Metastases from Ovarian Cancer: Current Evidence in Diagnosis, Treatment, and Prognosis
Previous Article in Special Issue
Genomic Database Analysis of Uterine Leiomyosarcoma Mutational Profile
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

MDM2-Dependent Rewiring of Metabolomic and Lipidomic Profiles in Dedifferentiated Liposarcoma Models

1
Department of Biomedical Informatics, The Ohio State University, Columbus, OH 43210, USA
2
Division of Preclinical Innovation, National Center for Advancing Translational Sciences, NIH, 9800 Medical Center Dr., Rockville, MD 20892, USA
3
Biomedical Sciences Graduate Program, The Ohio State University, Columbus, OH 43210, USA
4
Department of Internal Medicine, The Ohio State University, Columbus, OH 43210, USA
5
The Huck Institutes of Life Sciences, Penn State University, State College, PA 16802, USA
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
These authors contributed equally to this work.
Cancers 2020, 12(8), 2157; https://doi.org/10.3390/cancers12082157
Submission received: 26 June 2020 / Revised: 22 July 2020 / Accepted: 30 July 2020 / Published: 4 August 2020
(This article belongs to the Special Issue Management of Soft Tissue Sarcomas and GIST)

Abstract

:
Dedifferentiated liposarcoma (DDLPS) is an aggressive mesenchymal cancer marked by amplification of MDM2, an inhibitor of the tumor suppressor TP53. DDLPS patients with higher MDM2 amplification have lower chemotherapy sensitivity and worse outcome than patients with lower MDM2 amplification. We hypothesized that MDM2 amplification levels may be associated with changes in DDLPS metabolism. Six patient-derived DDLPS cell line models were subject to comprehensive metabolomic (Metabolon) and lipidomic (SCIEX 5600 TripleTOF-MS) profiling to assess associations with MDM2 amplification and their responses to metabolic perturbations. Comparing metabolomic profiles between MDM2 higher and lower amplification cells yielded a total of 17 differentially abundant metabolites across both panels (FDR < 0.05, log2 fold change < 0.75), including ceramides, glycosylated ceramides, and sphingomyelins. Disruption of lipid metabolism through statin administration resulted in a chemo-sensitive phenotype in MDM2 lower cell lines only, suggesting that lipid metabolism may be a large contributor to the more aggressive nature of MDM2 higher DDLPS tumors. This study is the first to provide comprehensive metabolomic and lipidomic characterization of DDLPS cell lines and provides evidence for MDM2-dependent differential molecular mechanisms that are critical factors in chemoresistance and could thus affect patient outcome.

1. Introduction

Dedifferentiated liposarcoma (DDLPS) is a highly morbid, adipocytic tumor accounting for approximately 20% of all soft-tissue sarcomas [1]. Liposarcomas are generally thought to arise spontaneously rather than from preexisting benign lesions, and most patients lack recognized causative factors. Although complete surgical resection can be curative, DDLPS often develops in deep anatomic locations, such as the retroperitoneum or mediastinum, where its propensity to encase vital structures typically renders a complete surgical resection difficult or impossible. In fact, the five year survival of patients with these abdominal liposarcomas is only 20% [2]. Unfortunately, chemotherapy has limited efficacy in the treatment of DDLPS, with single-agent response rates of up to 30% [3]. Systemic therapeutic regimens improve survival only modestly when complete surgical resection is not feasible [3]. Improved treatments are critically needed for this highly morbid disease.
At the molecular level, DDLPS is characterized by an amplification of the 12q portion of the chromosome resulting in excess copies of the mouse double minute 2 homolog (MDM2) [4,5,6,7]. MDM2 inhibits the tumor suppressor TP53. Thus, an amplification in MDM2 results in a shift towards pro-growth pathways. Our previous work demonstrated that higher levels of MDM2 amplification are associated with worsened overall survival and resistance to DNA-damaging chemotherapy in liposarcomas [8,9]. Interestingly, although DDLPS are of fat origin, they produce scant fat themselves and resemble undifferentiated pleomorphic or spindle cell sarcoma, typically showing moderate or high cellularity, with moderate to marked pleomorphism [10]. Given these observations, we hypothesized that central to the oncogenicity of DDLPS is its ability to alter fat metabolism and other key metabolic pathways in an MDM2-dependent manner. We thus performed a comprehensive metabolomic and lipidomic profiling of patient-derived DDLPS cell lines with varying degrees of MDM2 amplification. Taken together, this is the first attempt to characterize the metabolome and lipidome of DDLPS in light of each cell line’s genomic amplification milieu.

2. Results

To test our hypothesis, that treatment impact on DDLPS tumors differs for tumors with higher vs. lower MDM2 amplification, we measured the metabolomic and lipidomic effects of MDM2 amplification, MDM2 inhibition, cholesterol inhibition, and ceramide treatment in six patient-derived DDLPS cell lines (Table 1). DDLPS cell lines were categorized as MDM2 lower and higher amplification based on prior published PCR analyses [11]. Lipo-224, LPS141, and Lipo-246 cells have higher MDM2 amplification, and Lipo-815, Lipo-224B, and Lipo-863 have lower MDM2 amplification (Table 1, Figure S1a). MDM2 amplification levels were verified by DNA copy number and RNA transcription levels, showing high concordance (Pearson’s r = 0.92, P = 0.03) [11]. Treatment of the cell lines with varying doses of doxorubicin confirmed that MDM2 higher cell lines had lower chemosensitivity than MDM2 lower cells (Figure 1a).

2.1. Metabolomic Changes Associated with MDM2 Amplification Levels

Metabolomic profiles were generated using the Metabolon platform in at least triplicate measurements. Table S1 provides a summary of the conditions used. A total of 541 metabolites were measured (including 62 Metabolomics Standard Initiative (MSI) Level 1 identified metabolites [13]). Measured metabolites include lipids, amino acids, nucleotides, carbohydrates, peptides, cofactors and vitamins, xenobiotics, and energy metabolites, as defined by the Metabolon “superpathway” designation (Figure 1b). Unsupervised clustering of samples by their metabolomic profiles (Figure S1b) showed that biological replicate samples cluster closely together, providing confidence in the reproducibility of the measurements.
When comparing metabolite levels between MDM2 higher and lower cells, we identified 17 altered metabolites (FDR-adjusted p value < 0.05, log2 fold change (FC) > 0.75), of which 15 were elevated in MDM2 higher cells and two were elevated in MDM2 lower cells (Figure S1b Supplementary Materials). Of these, we noted that lipids constitute 41% (seven species) of the metabolites altered, and the rest comprise three amino acids, three peptides, two nucleotides, one vitamin- and cofactor-associated metabolite, and one energy-associated metabolite (Figure 1c and Figure S1d). Interestingly, pathway analysis revealed that sphingolipid metabolism/de novo biosynthesis were the top enriched pathways. Additional pathways prioritized included synthesis of prostaglandins and thromboxanes, beta-oxidation, and bile acid synthesis/secretion (Figure S1c and Figure 1d).

2.2. Raising the MDM2 Levels in DDLPS Lower Cell Lines Results in Induction of Key Metabolites Resembling Those of MDM2 Higher Cells

We previously demonstrated that MDM2-TP53 binding inhibitors temporarily result in further amplification of MDM2 by allowing TP53 to induce MDM2 expression [11]. Thus, MDM2-TP53 binding inhibitors temporarily increase MDM2 levels in both MDM2 higher and lower amplification DDLPS cell lines (Figure S2a). We evaluated the metabolomic profiles of DDLPS cell lines (LPS141, Lipo-246, Lipo-863, and Lipo-815) with and without MDM2-TP53 binding inhibitor RG7112. When comparing the effect of MDM2 inhibitor treatment in all cell lines (treated vs. untreated) regardless of MDM2 status, no significant metabolites were identified (Figure S2b).
When we considered the effect of MDM2 inhibitor treatment in MDM2 higher and lower cells separately, we observed six altered metabolites when comparing MDM2 lower cells before and after treatment. Of these six, four were upregulated by treatment and two were downregulated (Figure 2a,b and Figure S2c). No metabolites were altered in MDM2 higher cells (Figure 2c). The six metabolites altered by MDM2 binding inhibition in MDM2 lower cells comprised three lipids (1-(1-enyl-palmitoyl)-2-oleoyl-Glycerophosphoethanolamine (P-16:0/18:1), myristate (14:0), palmitoylcholine) and three nucleotides (adenylosuccinate, thymidine 5 -monophosphate, uridine). Pathways associated with these metabolites in the RaMPpathway database [12] included pyrimidine salvage, pyrimidine catabolism, and nucleotide salvage.

2.3. Independent Lipidomic Profiling Confirms Lipids Found to Be Altered by MDM2 Amplification from Metabolomic Analysis and Identifies Additional Relevant MDM2-Dependent Lipids

Given the larger numbers of lipids altered by MDM2 with our metabolomic analyses, we performed an independent lipidomic analysis (SCIEX 5600 TripleTOF-MS) in the same DDLPS cell line models grown analogously to the prior metabolomic profiling experiments. The quality of the experiment was assessed through PCA analysis and showed that pooled QC samples, blanks, and biological replicates clustered tightly together (Figure S3a). The lipid coverage in our metabolomic and lipidomic profiling was carefully assessed to identify common lipids covered by both platforms. Of all lipids measured, twenty-three percent were only identified in the Metabolon platform, 67% only in the lipidomic analysis, and 10% in both (Figure S3b). The correlation between the alterations of the ratios of abundance in MDM2 higher vs. lower cells in identical lipid species from both platforms was high (Pearson ρ = 0.67, Figure S3c,d), demonstrating that the measurements are reproducible and robust.
A total of 433 lipids were measured and identified, and lipid classes included a wide variety of glycerophospholipids such as glycerophosphocholines, glycerophosphoethanolamines, and glycerophosphoinositols, glycerolipids such as diradylglycerols and triacylglycerols, as well as sphingolipids such as sphingomyelins and ceramides (Figure 3a). When comparing lipid levels between MDM2 higher and lower amplification cells, seven lipids were altered, six of which were elevated in MDM2 low cells (monogalactosyldiacylglycerol (16:0_22:6), plasmenyl-Phosphatidylcholine (PC) (P-18:0/22:5), plasmenyl-PC(P-18:1/20:4), Sphingomyelin (SM) (d18:0/24:0), SM(d18:1/24:1), SM(d22:1/22:1)) and one of which was elevated in MDM2 high cells.
The results of the MDM2 higher vs. lower comparison in the lipidomic panel provided further evidence for perturbation of the sphingolipid metabolism pathway. Three sphingomyelins were elevated in MDM2 lower cells compared to MDM2 higher cells, in addition to two plasmalogens and one glycosyldiradylglycerol. Additionally, one hexosylceramide non-hydroxyfatty acid-sphingosine (HexCer_NS), (HexCer_NS (d18:1/16:0), p = 0.02) was significantly upregulated in MDM2 higher cells. This mirrored the results of the metabolomic analysis, where the same glycosylated ceramide (glycosyl-N-palmitoyl-sphingosine, p = 0.03) was also upregulated in MDM2 high compared to low. Three other glycosylated ceramides were borderline significant, showing log2 fold changes > 0.75, but unadjusted P-values < 0.05. When specifically assessing glycosylated ceramides and ceramides, we show that elevation of glycosylated ceramides in MDM2 higher cells is confirmed in both platforms (Figure S3c).

2.4. Induction of the Sphingolipid Pathway in DDLPS Models Using Atorvastatin Resulted in Chemoresistance

In previous literature, 3-hydroxy-3-methyl-glutaryl-CoA reductase reductase inhibitor (statin) treatment was noted to result in ceramide elevation [14,15,16]. Given the importance of the ceramide metabolites in our findings, we wanted to explore whether induction of ceramides altered chemosensitivity. To this end, we used the commonly prescribed HMG-CoA reductase inhibitor, atorvastatin. We observed that the growth of MDM2 higher cells was not affected by atorvastatin treatment, while MDM2 lower cells exhibited slower growth (Figure 3b). All DDLPS cell lines showed a decrease in cell viability, although the effect was more pronounced in MDM2 lower amplification cells (Figure S4a).
We then tested the effects of HMG-CoA reductase inhibition on the lipidome of both MDM2 higher and MDM2 lower cell lines (Figure 3c,d and Figure S4b). As expected, very few lipids were altered in MDM2 higher cells following treatment. Conversely, we observed dramatic shifts in the lipidome of MDM2 lower cells following atorvastatin treatment, with 52 lipids upregulated and 11 downregulated in response to treatment. Altered lipids included 11 ceramides, three glycosylated ceramides, and two sphingomyelins; all increased following atorvastatin treatment. Previous studies have demonstrated that statins can exert off-target effects on the sphingolipid metabolism pathway [17], which could explain the shifts observed in our data. Given that we noted elevated glycosylated ceramides, we hypothesized that this atorvastatin-driven shift could alter the chemosensitivity of the DDLPS models.
As a proof of concept, we examined the effects of the addition of atorvastatin on the cytotoxic effects of the anthracycline doxorubicin. Specifically, we tested a series of different concentrations of both atorvastatin and doxorubicin and then measured their effect on cell viability in the Lipo-246 DDLPS cell line model (MDM2 higher amplification). The cell viability in the combinations was more potent together than what we would have expected. Using the Chou–Talalay method [18], we note that the combination indices were nearly all above one indicating antagonistic agents. While both atorvastatin and doxorubicin both had inhibitory effects on the cell line viability as individual agents, when combined, we observed a cooperativity index > 1 for nearly all combinations, consistent with antagonism (Figure S4c).

2.5. Glycosylated Ceramides Are Consistently Elevated in MDM2 Higher Cells

Noting the elevation of glycosylated ceramides across panels in MDM2 higher cells compared to lower, we evaluated the sensitivity of the cell lines to treatment with non-glycosylated ceramides and found that the cell lines’ viability decreased, regardless of MDM2 status (Figure 4a). Glycosyl-N-palmitoyl-sphingosine, a glycosylated ceramide, showed elevated levels in MDM2 high cells. This was the same glycosylated ceramide that was also highly upregulated in MDM2 high cells in the lipidomic panel (HexCer_NS (d18:1/16:0), log2 FC = 2.36, FDR = 0.02, Figure 4b). Further, Lipid Ontology enrichment (LION/Web) [19] analysis of lipids perturbed between untreated MDM2 higher and MDM2 lower cells returned “ceramide phosphocholines (sphingomyelins)” as a significantly enriched term (Figure 4c). Lastly, to further validate differences in glycosylated ceramide levels between MDM2 higher and lower cells, we performed unsupervised clustering of cell lines in the lipidomic panel by abundance of all glycosylated ceramides (Figure 4d) and found that cells clustered completely by MDM2 amplification status.

3. Discussion

Prior studies have identified upregulation of nucleotide and serine synthesis [20] in liposarcoma compared to healthy cells, as well as the activity of the nucleoside salvage pathway [21]. We observed two nucleotides, three peptides, and two amino acids that were altered between MDM2 higher and lower cells, which could indicate that MDM2 higher and lower tumors could be differentially sensitive to the previously reported therapeutic targets of amino acid and nucleotide salvage and synthesis. We further observed significant differences in lipid metabolism between MDM2 higher and lower cells and MDM2-dependent shifts in cellular growth of DDLPS cells in response to modulation of lipid metabolism via statin treatment. Sphingolipid metabolism was consistently observed and is a known important therapeutic target in a variety of cancer types. For example, sphingomyelin synthase is frequently inhibited to prevent the conversion of ceramide to sphingomyelin to maintain its pro-apoptotic activity [22]. In the metabolomic panel, two lipids altered between MDM2 higher and lower cells mapped to the sphingolipid metabolism/glycosphingolipid metabolism pathways. We identified three sphingomyelins that were depleted in MDM2 higher cells through our lipidomic assay. Interestingly, while glycosylated ceramides were strong predictors of MDM2 amplification status, ceramides were not (results not shown), meaning that the ratio of glycosylated/nonglycosylated ceramides was shifted between MDM2 higher and lower cells. Glycosylated ceramides are known to be able to drive drug resistance without a change in the levels of their associated nonglycosylated form [23], counteracting the pro-apoptotic signaling exerted by ceramides [24]. This suggests a possible mechanism of resistance to atorvastatin in MDM2 higher cells driven by elevated ratios of glycosylated:nonglycosylated ceramides. These results are concordant with previous findings showing the elevation of ceramides in response to stress (e.g., chemotherapy) [25]. Sphingolipid metabolism is an important therapeutic target in a variety of cancer types. For example, sphingomyelin synthase is frequently inhibited to prevent the conversion of ceramide to sphingomyelin to maintain its pro-apoptotic activity [22]. Both our lipid and metabolite data provided evidence for the perturbation of sphingolipid metabolism in MDM2 higher cells compared to lower. These results are concordant with previous findings showing the elevation of ceramides in response to stress (e.g., chemotherapy), thereby leading to cell death [25].
We also note that MDM2 higher cells showed an increase in metabolites involved in de novo fatty acid synthesis. Many cancer types, such as breast, colorectal, and ovarian cancer, show an upregulation of fatty acid synthesis [26,27,28]. Upregulation of lipogenic enzymes that produce fatty acids has been associated with poor prognosis and resistance to chemotherapy [29]. Further, three of the lipids upregulated in MDM2 higher compared to lower cells were saturated fatty acids: myristate (14:0), stearate (18:0), and nonadecanoate (19:0). Saturated fatty acids are known drivers of chronic inflammation through TLR4/nf κ B-dependent signaling [30], which is a hallmark of cancers with increased risk [31]. This hypothesis is further supported by the identification of beta-oxidation and alpha-linoleic metabolism as perturbed pathways between MDM2 higher and lower cells, which is evidence of increased lipid metabolism/synthesis. Upregulation of beta oxidation is increasingly being recognized as a consistent feature of the cancer metabolic landscape [32,33]. Previous studies have indicated that DDLPS tumors have a distinct fatty acid composition compared to other liposarcoma cell types and benign lipomas [34].
In light of the different atorvastatin treatment responses displayed by MDM2 high and low cells (Figure 3a), the shift in bile acid pathways observed between MDM2 high and low cells was a significant finding. Bile acids are steroids generated from cholesterol in the liver, undergoing further metabolism in the gut by enzymes derived from intestinal bacteria [35]. In addition to the key role they play in intestinal uptake of lipids and vitamins, bile acids act as signaling molecules that regulate cell growth, as well as inflammation [36]. As such, perturbation of bile acid synthesis is known to be a factor in a variety of cancer types, such as colorectal cancer [37], hepatocellular carcinoma [38], and cholangiocarcinoma [39].

4. Materials and Methods

4.1. In Vitro Models

The culture of human LPS cell lines (Lipo-246, Lipo-863, Lipo-815, Lipo-224, Lipo-224B) has been previously reported [40]. Dr. Jonathan Fletcher (Boston, MA) generously provided us with the LPS141 cell line. All cells were cultured in Dulbecco’s Modified Eagle’s Medium (DMEM) and supplemented with 10% fetal bovine serum (FBS), 100 U/mL penicillin, and 100 U/mL streptomycin. These cells were cultured in a humidified chamber delivering 5% CO2 at 37 ° C.

4.2. Chemical Reagents

Doxorubicin was purchased from Cayman Biochemicals. The MDM2 inhibitor SAR405838 was purchased from Selleckchem. Atorvastatin was purchased from Sigma-Aldrich. All drugs were prepared per the manufacturers’ instructions. Serial dilutions were made to obtain final concentrations for cellular assays of DMSO not exceeding 0.01%.

4.3. Cell Proliferation via the MTT Assay and Cooperativity Evaluation

Exponentially growing DDLPS cell lines were seeded into 96 well plates and treated with the indicated compounds. After 24 h, treatment was added to the cells in the plate. The treatment conditions were as follows: doxorubicin only, atorvastatin only, constant atorvastatin doses paired with increasing doxorubicin doses, and constant doxorubicin doses paired with increasing atorvastatin doses. Doxorubicin was added in the following doses: 0.1 μ M, 0.3 μ M, 0.7 μ M, 1.2 μ M, 2 μ M. Atorvastatin was added in the following doses: 5 μ M, 10 μ M, 15 μ M, 20 μ M, 30 μ M. After 72 h, the XTT cell proliferation kit from Roche Applied Science (11465015001) was used to assess cell viability following the manufacturer’s instructions. The absorbance was measured at 470 nm. The calculation of effect size was done via Calcusyn software. Live cell imaging was performed using the Incucyte Zoom system.

4.4. Western Blotting

Western blots were performed using Odyssey CLx (Li-Cor) and ECL (PerkinElmer). The antibodies were used as indicated per experiment: p53, p21, β -actin (Santa Cruz); MDM2 (Abcam); cleaved caspase-3 (Cell Signaling).

4.5. Metabolomic and Lipidomic Data Acquisition

DDLPS were plated in 60 mm dishes, drug treated with atorvastatin or SAR405838, and cultured for 72 h. Because these treatments are not cytotoxic agents, peak changes are seen later as metabolic and genomic shifts take more time to manifest, requiring a longer time point for treatment. After treatment, cells were collected using 1 mL of cold (−20 ° C) methanol and a cell scraper and then immediately stored at −80 ° C. Cell suspensions were sent to Metabolon Inc. [41] for comprehensive metabolomic profiling and to the Penn State University Metabolomics Facility for lipidomic profiling.
For lipidomic profiling, cells were pelleted and processed using a chloroform:methanol homogenization followed by an isopropanol:acetonitrile extraction as previously described [42]. Samples were separated by reverse phase HPLC using a Prominence 20 UFLCXR system (Shimadzu, Columbia MD) with a Waters (Milford, MA) CSH C18 column (100 mm × 2.1 mm 1.7 μ m particle size) maintained at 55 ° C and a 20 m aqueous/acetonitrile/isopropanol gradient, at a flow rate of 225 μ L/min. For electrospray ionization positive mode, Solvent A was 40% water, 60% acetonitrile with 10 mM ammonium formate and 0.1% formic acid, and Solvent B was 90% isopropanol, 10% acetonitrile with 10 mM ammonium formate and 0.1% formic acid. For electrospray ionization negative mode, Solvent A was 40% water, 60% isopropanol with 10 mM ammonium acetate, and Solvent B was 90% isopropanol, 10% acetonitrile with 10 mM ammonium acetate. The initial conditions were 60% A and 40% B, increasing to 43% B at 2 min, 50% B at 2.1 min, 54% B at 12 min, 70% B at 12.1 min, and 99% B at 18 min, held at 99% B until 20.0 min, before returning to the initial conditions. The eluate was delivered into a 5600 (QTOF) TripleTOF using a Duospray™ ion source (all AB Sciex, Framingham, MA, USA). The capillary voltage was set at 5.5 kV in positive ion mode and 4.5 kV in negative ion mode, with a declustering potential of 80 V. The mass spectrometer was operated in IDA (information-dependent acquisition) mode with a 100 ms survey scan from 100 to 1200 m/z and up to 20 MS/MS product ion scans (100 ms) per duty cycle using a collision energy of 50 V with a 20 V spread.

4.6. Metabolomic Data Preprocessing

Metabolite levels were pre-processed by Metabolon, including sample normalization by Bradford protein concentration, median scaling, and missing value imputation to minimum values. There were a total of 541 identified metabolites. Metabolites with 50% missing values or greater in all samples were removed, eliminating 33 metabolites. Metabolites with a high coefficient of variation (greater than 150) were also removed, as extremely variable metabolites may be a result of technical error rather than actual biological variation. After applying these filters, four-hundred eighty metabolites remained for further analysis.

4.7. Lipidomic Data Preprocessing

The Proteowizard software suite [43] was used to convert the wiff/wiff.scan raw data files to mzML/ms2 file formats, using the following (1) peakpicking filters: pickertype=cwt, signal-to-noise ratio = 0.1, and peakSpace = 0.1 and 2) conversion filters: mslevel 1-1 was used for .mzl conversion, and mslevel 2-2 was used for .ms2 conversion. Converted files were then input into MZmine2 Version 2.42 [44] for peak calling and chromatogram alignment. Mass detection with a centroid mass detector on MS Level 1 was performed. Next, the ADAP chromatogram builder was used, with a minimum scan span of 3 above a group intensity threshold of 1000, a minimum highest intensity of 1000, and an m/z tolerance of 0.005 m/z or 10 ppm. Chromatogram deconvolution was performed with the ADAP wavelet algorithm, with a minimum feature height of 1000, a coefficient/area threshold of 10, a peak duration range of 0–2, and a retention time wavelet range of 0–0.1. Next, isotopic peak grouping was performed, with an m/z tolerance of 0.005 m/z or 5 ppm, a retention time tolerance of 1% relative, and a maximum charge of 3. Next, a duplicate feature filter was applied, with a filter mode of “NEW AVERAGE”, an m/z tolerance of 0.005 m/z or 5 ppm, and a retention time tolerance of 1% relative. Join alignment was performed with an m/z tolerance of 0.009 m/z or 10 ppm, a weight for m/z of 85, a retention time tolerance of 1% relative, and a weight for RT of 85. Finally, the same RT and m/z range gap filler was applied, with an m/z tolerance of 0.005 m/z or 5 ppm. This generated csv feature tables for the positive and negative ion modes, which were used as input for peak identification using the LipidMatch software [45].
For LipidMatch lipid identification, the default Sciex parameters were used. A retention time window of 0.14 m was used. An m/z search tolerance of 0.005 Da was used for MS1, and an M/Z search window of 10 ppm was used for MS2. Isolation window was set to 1 Da. Minimum scans for fragments were set to 1. The intensity threshold for MS2 was set to 1000. The All-ion-fragmentation (AIF)minimum number of scans was set to 5. The AIF correlation cutoff (adjusted R2 correlation between precursors and fragments) was set to 0.6. The features tables were configured so that they matched the default order of columns and rows expected by LipidMatch (Column 1 is row ID; Column 2 is M/Z; Column 3 is retention time; numeric data starts on Row 2). Abundance values were normalized for each ionization mode independently using the total ion current method [46]. Normalized abundances were log2 transformed, and data from both ionization modes were merged. Further analysis was restricted to lipids identified from MS2 data. The highest abundance feature was selected for further analysis when multiple features had the same identification. If a lipid was identified in both modes, negative adducts were preferred over positive adducts. Sodium adducts were only used if matching identifications with other adducts were not found. A total of 430 unique lipids were identified and quantified.

4.8. Statistical Analysis

Differentially abundant metabolites and lipids between groups (e.g., MDM2 status, treatment status, etc.) were assessed using linear mixed effect modeling with cell line of origin incorporated as a random effect. Resulting P-values were corrected for multiple comparisons using the method by Benjamini-Hochberg to correct the false discovery rate (FDR) [47]. Metabolites were considered statistically significant when their associated FDR-adjusted P-values were greater than 0.05 and absolute log2 fold change > 0.75. Results of all statistical testing can be found in Supplementary Materials Files 2 and 3.

4.9. Pathway Analysis

We developed a novel network-based pathway enrichment strategy that integrates pathway annotations from the RaMP database [12], as well as chemical structure similarity. Metabolites associated with phenotype (e.g., MDM2 amplification or treatment) that mapped to KEGG pathways [48,49,50] were used as the input. Ten out of 18 metabolites altered by MDM2 amplification mapped to KEGG pathways. We obtained a list of all KEGG pathways that contained at least one of these altered metabolites. We then obtained a list of all metabolites involved in this list of pathways. We used these metabolites to build a mutual pathway participation network, where nodes were metabolites and edges encoded pathway similarity between metabolites. Pathway similarity was quantified using a Jaccard index of pathway annotations shared by metabolites connected by an edge. A chemical similarity network was constructed between this same set of metabolites, where edges represented the chemical fingerprint overlap between two metabolites using a Tanimoto score. The chemical similarity network was then binarized, with edges in the 90th percentile or higher of chemical similarity receiving an edge weight of 1 and all other edge weights assigned a 0. The two separate similarity networks were then merged into a consensus network model by summing edge weights. We then searched for metabolites proximal to our set of ten or nine altered metabolites using a random walk with restarts strategy [51]. To control for the connectivity of the network structure, we ran the algorithm with 10,000 random seed sets of ten or nine and compared the proximity score of nodes with the real seed set to the distribution of random sets. Metabolites that were scored very high compared to their randomized distribution (97th percentile or higher) were extracted to form a subnetwork with the seed nodes. Lastly, we applied the Louvain clustering algorithm to identify highly interconnected clusters (e.g., red, blue, and green clusters in Figure S1c). Metabolites in clusters were tested separately for pathway enrichment using RaMP (FDR-adjusted p-value < 0.05).

4.10. Code and Data Availability

Preprocessed data along with the R code used for the analysis can be found at https://github.com/andyptt21/MDM2_reprograms_DDLPS_metabolism. Raw lipidomics data were submitted to Metabolomics Workbench under study track ST001405.

5. Conclusions

To our knowledge, this is the first detailed examination of the metabolome/lipidome of dedifferentiated liposarcoma cell lines. As such, our experimental approach allowed us to interrogate the metabolite and lipid differences that could underlie the heightened chemoresistance of MDM2 higher tumor cells as compared to lower. Our metabolomic/lipidomic analyses uncovered a wide variety of differences between MDM2 higher and lower cells, and in response atorvastatin treatment. Many of the pathway level shifts we observed suggest that MDM2 plays a role in a wide variety of cancer-related metabolic processes, including ceramide metabolism, amino acid synthesis, lipogenesis, and inflammation. The mirrored and complementary results we observed between the metabolomic and lipidomic datasets served as the validation of our findings. These data serve as a key basis for the further development of metabolic targets in this morbid disease, as well as potential clinical work. Metabolically-driven anti-cancer therapies based on these observations may potentially spare untargeted toxicities to the patient in this difficult to treat disease. Future research focus will include examination of mouse xenograft models of DDLPS and of patient tissue for validation of our in vitro findings of perturbation in ceramide metabolism.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6694/12/8/2157/s1, Figure S1: MDM2 status and metabolomic (Metabolon) profiling in DDLPS cell lines, Figure S2: Comparing the effect of MDM2 inhibitor in DDLPS cells, Figure S3: Quality assessment of lipidomic data, Figure S4: Effect of MDM2 amplification/Atorvastatin treatment in DDLPS cells, Table S1: Sarcoma cell line characteristics for metabolomic and lipidomic profiling.

Author Contributions

Conceptualization, J.L.C. and E.A.M.; formal analysis, A.P. (Andrew Patt), B.D., C.S., K.-L.B., P.S., A.V., A.P. (Andrew Patterson), J.H., M.H., and E.A.M.; funding acquisition, J.L.C. and E.A.M.; investigation, B.D., P.S., J.L.C., and E.A.M.; project administration, J.L.C.; resources, P.S., A.P. (Andrew Patterson), J.L.C., and E.A.M.; visualization, A.P. (Andrew Patt), and B.D.; writing, original draft, A.P. (Andrew Patt) and J.L.C.; writing, review and editing, A.P. (Andrew Patt), B.D., C.S., K.-L.B., A.V., A.P. (Andrew Patterson), J.H., M.H., and E.A.M. All authors read and agreed to the published version of the manuscript.

Funding

This research was supported in part by the Intramural/Extramural research program of the NCATS, NIH. Funding for this project was also provided in part by Pelotonia, a non-profit cancer charity affiliated with Ohio State University. Lastly, this work was supported by the National Institute Of General Medical Sciences of the National Institutes of Health under Award Number 2T32GM068412-11A1, as well as the National Library of Medicine grant 4T15LM011270-05. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Acknowledgments

We would like to thank Jeremy Koelmel, University of Florida, for technical support in using LipidMatch.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jo, V.Y.; Fletcher, C.D. WHO classification of soft tissue tumours: An update based on the 2013 (4th) edition. Pathology 2014, 46, 95–104. [Google Scholar] [CrossRef] [PubMed]
  2. Fabre-Guillevin, E.; Coindre, J.M.; Somerhausen, N.d.S.A.; Bonichon, F.; Stoeckle, E.; Bui, N.B. Retroperitoneal liposarcomas: Follow-up analysis of dedifferentiation after clinicopathologic reexamination of 86 liposarcomas and malignant fibrous histiocytomas. Cancer 2006, 106, 2725–2733. [Google Scholar] [CrossRef] [PubMed]
  3. Schöffski, P.; Ray-Coquard, I.L.; Cioffi, A.; Bui, N.B.; Bauer, S.; Hartmann, J.T.; Krarup-Hansen, A.; Grünwald, V.; Sciot, R.; Dumez, H.; et al. Activity of eribulin mesylate in patients with soft-tissue sarcoma: A phase 2 study in four independent histological subtypes. Lancet Oncol. 2011, 12, 1045–1052. [Google Scholar] [CrossRef]
  4. Beird, H.C.; Wu, C.C.; Ingram, D.R.; Wang, W.L.; Alimohamed, A.; Gumbs, C.; Little, L.; Song, X.; Feig, B.W.; Roland, C.L.; et al. Genomic profiling of dedifferentiated liposarcoma compared to matched well-differentiated liposarcoma reveals higher genomic complexity and a common origin. Cold Spring Harb. Mol. Case Stud. 2018, 4, a002386. [Google Scholar] [CrossRef] [Green Version]
  5. Nilbert, M.; Rydholm, A.; Mitelman, F.; Meltzer, P.S.; Mandahl, N. Characterization of the 12q13-15 amplicon in soft tissue tumors. Cancer Genet. Cytogenet. 1995, 83, 32–36. [Google Scholar] [CrossRef]
  6. Binh, M.B.; Sastre-Garau, X.; Guillou, L.; de Pinieux, G.; Terrier, P.; Lagacé, R.; Aurias, A.; Hostein, I.; Coindre, J.M. MDM2 and CDK4 immunostainings are useful adjuncts in diagnosing well-differentiated and dedifferentiated liposarcoma subtypes: A comparative analysis of 559 soft tissue neoplasms with genetic data. Am. J. Surg. Pathol. 2005, 29, 1340–1347. [Google Scholar] [CrossRef]
  7. Crago, A.M.; Singer, S. Clinical and molecular approaches to well differentiated and dedifferentiated liposarcoma. Curr. Opin. Oncol. 2011, 23, 373–378. [Google Scholar] [CrossRef] [Green Version]
  8. Ricciotti, R.W.; Baraff, A.J.; Jour, G.; Kyriss, M.; Wu, Y.; Liu, Y.; Li, S.C.; Hoch, B.; Liu, Y.J. High amplification levels of MDM2 and CDK4 correlate with poor outcome in patients with dedifferentiated liposarcoma: A cytogenomic microarray analysis of 47 cases. Cancer Genet. 2017, 218–219, 69–80. [Google Scholar] [CrossRef]
  9. Bill, K.L.J.; Seligson, N.D.; Hays, J.L.; Awasthi, A.; Demoret, B.; Stets, C.W.; Duggan, M.C.; Bupathi, M.; Brock, G.N.; Millis, S.Z.; et al. Degree of MDM2 Amplification Affects Clinical Outcomes in Dedifferentiated Liposarcoma. Oncologist 2019, 24, 989–996. [Google Scholar] [CrossRef]
  10. Thway, K. Well-differentiated liposarcoma and dedifferentiated liposarcoma: An updated review. Semin Diagn. Pathol. 2019, 36, 112. [Google Scholar] [CrossRef]
  11. Bill, K.L.; Garnett, J.; Meaux, I.; Ma, X.; Creighton, C.J.; Bolshakov, S.; Barriere, C.; Debussche, L.; Lazar, A.J.; Prudner, B.C.; et al. SAR405838: A Novel and Potent Inhibitor of the MDM2:p53 Axis for the Treatment of Dedifferentiated Liposarcoma. Clin. Cancer Res. 2016, 22, 1150–1160. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Zhang, B.; Hu, S.; Baskin, E.; Patt, A.; Siddiqui, J.K.; Mathé, E.A. RaMP: A Comprehensive Relational Database of Metabolomics Pathways for Pathway Enrichment Analysis of Genes and Metabolites. Metabolites 2018, 8, 16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Sumner, L.W.; Amberg, A.; Barrett, D.; Beale, M.H.; Beger, R.; Daykin, C.A.; Fan, T.W.; Fiehn, O.; Goodacre, R.; Griffin, J.L.; et al. Proposed minimum reporting standards for chemical analysis Chemical Analysis Working Group (CAWG) Metabolomics Standards Initiative (MSI). Metabolomics 2007, 3, 211–221. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Al Mamun Bhuyan, A.; Nüßle, S.; Cao, H.; Zhang, S.; Lang, F. Simvastatin, a Novel Stimulator of Eryptosis, the Suicidal Erythrocyte Death. Cell. Physiol. Biochem. 2017, 43, 492–506. [Google Scholar] [CrossRef]
  15. Binnington, B.; Nguyen, L.; Kamani, M.; Hossain, D.; Marks, D.L.; Budani, M.; Lingwood, C.A. Inhibition of Rab prenylation by statins induces cellular glycosphingolipid remodeling. Glycobiology 2016, 26, 166–180. [Google Scholar] [CrossRef] [Green Version]
  16. Goulitquer, S.; Croyal, M.; Lalande, J.; Royer, A.L.; Guitton, Y.; Arzur, D.; Durand, S.; Le Jossic-Corcos, C.; Bouchereau, A.; Potin, P.; et al. Consequences of blunting the mevalonate pathway in cancer identified by a pluri-omics approach. Cell Death Dis. 2018, 9, 745. [Google Scholar] [CrossRef]
  17. Bergheanu, S.C.; Reijmers, T.; Zwinderman, A.H.; Bobeldijk, I.; Ramaker, R.; Liem, A.H.; van der Greef, J.; Hankemeier, T.; Jukema, J.W. Lipidomic approach to evaluate rosuvastatin and atorvastatin at various dosages: Investigating differential effects among statins. Curr. Med. Res. Opin. 2008, 24, 2477–2487. [Google Scholar] [CrossRef]
  18. Chou, T.C. Drug combination studies and their synergy quantification using the Chou-Talalay method. Cancer Res. 2010, 70, 440–446. [Google Scholar] [CrossRef] [Green Version]
  19. Molenaar, M.R.; Jeucken, A.; Wassenaar, T.A.; van de Lest, C.H.A.; Brouwers, J.F.; Helms, J.B. LION/web: A web-based ontology enrichment tool for lipidomic data analysis. Gigascience 2019, 8, giz061. [Google Scholar] [CrossRef] [Green Version]
  20. Cisse, M.Y.; Pyrdziak, S.; Firmin, N.; Gayte, L.; Heuillet, M.; Bellvert, F.; Fuentes, M.; Delpech, H.; Riscal, R.; Arena, G.; et al. Targeting MDM2-dependent serine metabolism as a therapeutic strategy for liposarcoma. Sci. Transl. Med. 2020, 12, eaay2163. [Google Scholar] [CrossRef]
  21. Braas, D.; Ahler, E.; Tam, B.; Nathanson, D.; Riedinger, M.; Benz, M.R.; Smith, K.B.; Eilber, F.C.; Witte, O.N.; Tap, W.D.; et al. Metabolomics strategy reveals subpopulation of liposarcomas sensitive to gemcitabine treatment. Cancer Discov. 2012, 2, 1109–1117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Pashikanti, S.; Afrin, F.; Meldrum, T.C.; Stegelmeier, J.L.; Pavek, A.; Habashi, Y.A.; Fatema, K.; Barrott, J.J. Quantifying Fluorescently Labeled Ceramide Levels in Human Sarcoma Cell Lines in Response to a Sphingomyelin Synthase Inhibitor. Methods Protoc. 2019, 2, 76. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Inokuchi, J.; Mason, I.; Radin, N.S. Antitumor activity via inhibition of glycosphingolipid biosynthesis. Cancer Lett. 1987, 38, 23–30. [Google Scholar] [CrossRef] [Green Version]
  24. Liu, Y.Y.; Hill, R.A.; Li, Y.T. Ceramide glycosylation catalyzed by glucosylceramide synthase and cancer drug resistance. Adv. Cancer Res. 2013, 117, 59–89. [Google Scholar] [PubMed] [Green Version]
  25. Ogretmen, B. Sphingolipid metabolism in cancer signalling and therapy. Nat. Rev. Cancer 2018, 18, 33–50. [Google Scholar] [CrossRef] [PubMed]
  26. Kuo, C.Y.; Ann, D.K. When fats commit crimes: Fatty acid metabolism, cancer stemness and therapeutic resistance. Cancer Commun. 2018, 38, 47. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Jafari, N.; Drury, J.; Morris, A.J.; Onono, F.O.; Stevens, P.D.; Gao, T.; Liu, J.; Wang, C.; Lee, E.Y.; Weiss, H.L.; et al. De Novo Fatty Acid Synthesis-Driven Sphingolipid Metabolism Promotes Metastatic Potential of Colorectal Cancer. Mol. Cancer Res. 2019, 17, 140–152. [Google Scholar] [CrossRef] [Green Version]
  28. Zhao, J.; Zhang, X.; Gao, T.; Wang, S.; Hou, Y.; Yuan, P.; Yang, Y.; Yang, T.; Xing, J.; Li, J.; et al. SIK2 enhances synthesis of fatty acid and cholesterol in ovarian cancer cells and tumor growth through PI3K/Akt signaling pathway. Cell Death Dis. 2020, 11, 25. [Google Scholar] [CrossRef] [Green Version]
  29. Menendez, J.A.; Lupu, R. Fatty acid synthase (FASN) as a therapeutic target in breast cancer. Expert Opin. Ther. Targets 2017, 21, 1001–1016. [Google Scholar] [CrossRef]
  30. Li, B.; Leung, J.C.K.; Chan, L.Y.Y.; Yiu, W.H.; Tang, S.C.W. A global perspective on the crosstalk between saturated fatty acids and Toll-like receptor 4 in the etiology of inflammation and insulin resistance. Prog. Lipid Res. 2020, 77, 101020. [Google Scholar] [CrossRef]
  31. Taniguchi, K.; Karin, M. NF-κB, inflammation, immunity and cancer: Coming of age. Nat. Rev. Immunol. 2018, 18, 309–324. [Google Scholar] [CrossRef] [PubMed]
  32. Wang, T.; Fahrmann, J.F.; Lee, H.; Li, Y.J.; Tripathi, S.C.; Yue, C.; Zhang, C.; Lifshitz, V.; Song, J.; Yuan, Y.; et al. JAK/STAT3-Regulated Fatty Acid-Oxidation Is Critical for Breast Cancer Stem Cell Self-Renewal and Chemoresistance. Cell Metab. 2018, 27, 136–150. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Ma, Y.; Temkin, S.M.; Hawkridge, A.M.; Guo, C.; Wang, W.; Wang, X.Y.; Fang, X. Fatty acid oxidation: An emerging facet of metabolic transformation in cancer. Cancer Lett. 2018, 435, 92–100. [Google Scholar] [CrossRef] [PubMed]
  34. Singer, S.; Millis, K.; Souza, K.; Fletcher, C. Correlation of lipid content and composition with liposarcoma histology and grade. Ann. Surg. Oncol. 1997, 4, 557–563. [Google Scholar] [CrossRef]
  35. Ridlon, J.M.; Kang, D.J.; Hylemon, P.B. Bile salt biotransformations by human intestinal bacteria. J. Lipid Res. 2006, 47, 241–259. [Google Scholar] [CrossRef] [Green Version]
  36. Li, T.; Apte, U. Bile Acid Metabolism and Signaling in Cholestasis, Inflammation, and Cancer. Adv. Pharmacol. 2015, 74, 263–302. [Google Scholar]
  37. Jia, W.; Xie, G.; Jia, W. Bile acid-microbiota crosstalk in gastrointestinal inflammation and carcinogenesis. Nat. Rev. Gastroenterol. Hepatol. 2018, 15, 111–128. [Google Scholar] [CrossRef] [Green Version]
  38. El-Mir, M.Y.; Badia, M.D.; Luengo, N.; Monte, M.J.; Marin, J.J. Increased levels of typically fetal bile acid species in patients with hepatocellular carcinoma. Clin. Sci. 2001, 100, 499–508. [Google Scholar] [CrossRef]
  39. Alsaleh, M.; Sithithaworn, P.; Khuntikeo, N.; Loilome, W.; Yongvanit, P.; Chamadol, N.; Hughes, T.; O’Connor, T.; Andrews, R.H.; Holmes, E.; et al. Characterisation of the Urinary Metabolic Profile of Liver Fluke-Associated Cholangiocarcinoma. J. Clin. Exp. Hepatol. 2019, 9, 657–675. [Google Scholar] [CrossRef] [Green Version]
  40. Peng, T.; Zhang, P.; Liu, J.; Nguyen, T.; Bolshakov, S.; Belousov, R.; Young, E.D.; Wang, X.; Brewer, K.; López-Terrada, D.H.; et al. An experimental model for the study of well-differentiated and dedifferentiated liposarcoma; deregulation of targetable tyrosine kinase receptors. Lab. Investig. 2011, 91, 392–403. [Google Scholar] [CrossRef]
  41. Evans, A.M.; DeHaven, C.D.; Barrett, T.; Mitchell, M.; Milgram, E. Integrated, nontargeted ultrahigh performance liquid chromatography/electrospray ionization tandem mass spectrometry platform for the identification and relative quantification of the small-molecule complement of biological systems. Anal. Chem. 2009, 81, 6656–6667. [Google Scholar] [CrossRef] [PubMed]
  42. Folch, J.; Lees, M.; Sloane Stanley, G.H. A simple method for the isolation and purification of total lipides from animal tissues. J. Biol. Chem. 1957, 226, 497–509. [Google Scholar] [PubMed]
  43. Chambers, M.C.; Maclean, B.; Burke, R.; Amodei, D.; Ruderman, D.L.; Neumann, S.; Gatto, L.; Fischer, B.; Pratt, B.; Egertson, J.; et al. A cross-platform toolkit for mass spectrometry and proteomics. Nat. Biotechnol. 2012, 30, 918–920. [Google Scholar] [CrossRef] [PubMed]
  44. Katajamaa, M.; Miettinen, J.; Orešič, M. MZmine: Toolbox for processing and visualization of mass spectrometry based molecular profile data. Bioinformatics 2006, 22, 634–636. Available online: http://xxx.lanl.gov/abs/https://academic.oup.com/bioinformatics/article-pdf/22/5/634/539891/btk039.pdf (accessed on 20 May 2020). [CrossRef] [Green Version]
  45. Koelmel, J.P.; Kroeger, N.M.; Ulmer, C.Z.; Bowden, J.A.; Patterson, R.E.; Cochran, J.A.; Beecher, C.W.W.; Garrett, T.J.; Yost, R.A. LipidMatch: An automated workflow for rule-based lipid identification using untargeted high-resolution tandem mass spectrometry data. BMC Bioinform. 2017, 18, 331. [Google Scholar] [CrossRef]
  46. Fung, E.T.; Enderwick, C. ProteinChip clinical proteomics: Computational challenges and solutions. BioTechniques 2002, 34–38. [Google Scholar] [CrossRef] [Green Version]
  47. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B Methodol. 1995, 57, 289–300. [Google Scholar] [CrossRef]
  48. Kanehisa, M.; Furumichi, M.; Tanabe, M.; Sato, Y.; Morishima, K. KEGG: New perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017, 45, D353–D361. [Google Scholar] [CrossRef] [Green Version]
  49. Kanehisa, M.; Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef]
  50. Kanehisa, M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019, 28, 1947–1951. [Google Scholar] [CrossRef]
  51. Dirmeier, S. Diffusr: Network Diffusion Algorithms; R Package Version 0.1.4; 2018; Available online: https://rdrr.io/cran/diffusr/ (accessed on 4 August 2020).
Figure 1. DDLPS cell lines with higher or lower MDM2 amplification show distinct chemosensitivity and metabolomic profiling. (a) Cell viability of DDLPS cell lines is elevated in MDM2 higher cells after doxorubicin treatment, (b) Distribution of metabolite classes (Metabolon superpathway) represented in the Metabolon panel. (c) Metabolites that are altered between MDM2 higher and lower amplification cell lines (FDR-adjusted p-value < 0.05 and |log2 fold change| > 0.75). (d) Over-represented pathways (FDR-adjusted p-value < 0.05) in metabolites that are altered between cells with higher vs. lower MDM2 amplification in the blue module of the final network enrichment model (Figure S1c). Pathways are grouped by cluster as determined by the relational database of metabolomics pathways (RaMP) pathway clustering algorithm [12].
Figure 1. DDLPS cell lines with higher or lower MDM2 amplification show distinct chemosensitivity and metabolomic profiling. (a) Cell viability of DDLPS cell lines is elevated in MDM2 higher cells after doxorubicin treatment, (b) Distribution of metabolite classes (Metabolon superpathway) represented in the Metabolon panel. (c) Metabolites that are altered between MDM2 higher and lower amplification cell lines (FDR-adjusted p-value < 0.05 and |log2 fold change| > 0.75). (d) Over-represented pathways (FDR-adjusted p-value < 0.05) in metabolites that are altered between cells with higher vs. lower MDM2 amplification in the blue module of the final network enrichment model (Figure S1c). Pathways are grouped by cluster as determined by the relational database of metabolomics pathways (RaMP) pathway clustering algorithm [12].
Cancers 12 02157 g001
Figure 2. p53 reactivation using an MDM2 binding inhibitor causes lipid upregulation in MDM2 lower cell lines. (a) Metabolites altered in DDLPS MDM2 lower amplification cells that are untreated or treated with MDM2 binding inhibition (RG7112) (FDR-adjusted p-val < 0.05 and |log2 fold change| > 0.75). (b) Volcano plot demonstrating shifts in metabolite levels caused by MDM2 inhibitor treatment in MDM2 lower cells. (c) No statistically significant metabolites were identified in MDM2 higher cells.
Figure 2. p53 reactivation using an MDM2 binding inhibitor causes lipid upregulation in MDM2 lower cell lines. (a) Metabolites altered in DDLPS MDM2 lower amplification cells that are untreated or treated with MDM2 binding inhibition (RG7112) (FDR-adjusted p-val < 0.05 and |log2 fold change| > 0.75). (b) Volcano plot demonstrating shifts in metabolite levels caused by MDM2 inhibitor treatment in MDM2 lower cells. (c) No statistically significant metabolites were identified in MDM2 higher cells.
Cancers 12 02157 g002
Figure 3. Effects of statin treatment on DDLPS cell lines depends on MDM2 status. (a) Coverage of lipids, categorized by their LipidMaps superclass. (b) Live cell imaging of DDLPS cell lines demonstrates differential growth patterns stratified by MDM2 status after atorvastatin treatment (average triplicate experiments). (c,d) We observed a more prominent dysregulation of lipids in response to atorvastatin treatment in MDM2 lower cells (c), compared to MDM2 higher cells (d).
Figure 3. Effects of statin treatment on DDLPS cell lines depends on MDM2 status. (a) Coverage of lipids, categorized by their LipidMaps superclass. (b) Live cell imaging of DDLPS cell lines demonstrates differential growth patterns stratified by MDM2 status after atorvastatin treatment (average triplicate experiments). (c,d) We observed a more prominent dysregulation of lipids in response to atorvastatin treatment in MDM2 lower cells (c), compared to MDM2 higher cells (d).
Cancers 12 02157 g003
Figure 4. Ceramides are an important differentiator of DDLPS subtypes. (a) Treatment with ceramides is lethal to DDLPS cells, regardless of MDM2 status. (b) Glycosylated ceramides and lipids in the sphingolipid pathways were amongst the most altered lipids between MDM2 high and low cells in the metabolomic panel. (c) LIONenrichment analysis of terms associated with lipids found different between MDM2 higher and lower cells, untreated. (d) Hierarchical clustering of cell lines by abundance of all glycosylated ceramides in the lipidomic panel. For column colors, red is MDM2 higher and blue is MDM2 lower.
Figure 4. Ceramides are an important differentiator of DDLPS subtypes. (a) Treatment with ceramides is lethal to DDLPS cells, regardless of MDM2 status. (b) Glycosylated ceramides and lipids in the sphingolipid pathways were amongst the most altered lipids between MDM2 high and low cells in the metabolomic panel. (c) LIONenrichment analysis of terms associated with lipids found different between MDM2 higher and lower cells, untreated. (d) Hierarchical clustering of cell lines by abundance of all glycosylated ceramides in the lipidomic panel. For column colors, red is MDM2 higher and blue is MDM2 lower.
Cancers 12 02157 g004
Table 1. Clinical and molecular characteristics of DDLPS cell lines.
Table 1. Clinical and molecular characteristics of DDLPS cell lines.
Cell LineMDM2 mRNA LevelMDM2 Amplification LevelGenderAge
LPS141473.4HighM80
Lipo-246583.1HighM60
Lipo-224A345.3HighF81
Lipo-224B169.9LowF81
Lipo-815106.1LowM66
Lipo-86379.4LowM74

Share and Cite

MDPI and ACS Style

Patt, A.; Demoret, B.; Stets, C.; Bill, K.-L.; Smith, P.; Vijay, A.; Patterson, A.; Hays, J.; Hoang, M.; Chen, J.L.; et al. MDM2-Dependent Rewiring of Metabolomic and Lipidomic Profiles in Dedifferentiated Liposarcoma Models. Cancers 2020, 12, 2157. https://doi.org/10.3390/cancers12082157

AMA Style

Patt A, Demoret B, Stets C, Bill K-L, Smith P, Vijay A, Patterson A, Hays J, Hoang M, Chen JL, et al. MDM2-Dependent Rewiring of Metabolomic and Lipidomic Profiles in Dedifferentiated Liposarcoma Models. Cancers. 2020; 12(8):2157. https://doi.org/10.3390/cancers12082157

Chicago/Turabian Style

Patt, Andrew, Bryce Demoret, Colin Stets, Kate-Lynn Bill, Philip Smith, Anitha Vijay, Andrew Patterson, John Hays, Mindy Hoang, James L. Chen, and et al. 2020. "MDM2-Dependent Rewiring of Metabolomic and Lipidomic Profiles in Dedifferentiated Liposarcoma Models" Cancers 12, no. 8: 2157. https://doi.org/10.3390/cancers12082157

APA Style

Patt, A., Demoret, B., Stets, C., Bill, K. -L., Smith, P., Vijay, A., Patterson, A., Hays, J., Hoang, M., Chen, J. L., & Mathé, E. A. (2020). MDM2-Dependent Rewiring of Metabolomic and Lipidomic Profiles in Dedifferentiated Liposarcoma Models. Cancers, 12(8), 2157. https://doi.org/10.3390/cancers12082157

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop