Epistatic interaction between the lipase-encoding genes and causes liposarcoma in mice

Liposarcoma is an often fatal adult-onset tumor of fat tissue. Lipolysis, the central pathway of fat tissue metabolism, has been implicated in cancer. We generated mice that were deficient in two key enzymes of lipolysis, adipose triglyceride lipase (ATGL) and hormone-sensitive lipase (HSL). Strikingly, all mice with combined ATGL and HSL deficiency developed liposarcoma by 11–14 months of age. No liposarcoma occurred in single knockout or normal controls. Transcriptome analysis revealed that a subset of genes is dysregulated by 3 months of age. Our study reveals a novel epistatic interaction in fat cells between these two lipase genes and that causes a unique form of liposarcoma in mice. The double knockout mice provide a novel tool to study the early stages of liposarcoma development, prognostic markers and preventive treatments.

Published in the journal: . PLoS Genet 13(5): e32767. doi:10.1371/journal.pgen.1006716
Category: Research Article
doi: 10.1371/journal.pgen.1006716


Liposarcoma is an often fatal adult-onset tumor of fat tissue. Lipolysis, the central pathway of fat tissue metabolism, has been implicated in cancer. We generated mice that were deficient in two key enzymes of lipolysis, adipose triglyceride lipase (ATGL) and hormone-sensitive lipase (HSL). Strikingly, all mice with combined ATGL and HSL deficiency developed liposarcoma by 11–14 months of age. No liposarcoma occurred in single knockout or normal controls. Transcriptome analysis revealed that a subset of genes is dysregulated by 3 months of age. Our study reveals a novel epistatic interaction in fat cells between these two lipase genes and that causes a unique form of liposarcoma in mice. The double knockout mice provide a novel tool to study the early stages of liposarcoma development, prognostic markers and preventive treatments.


Liposarcomas are malignant tumors of fat. It is the commonest soft tissue sarcoma and the annual incidence is 2.5 cases per million [1]. Specific biomarkers for early diagnosis and specific curative treatments are not available for the disease. The WHO classification of liposarcomas identifies four categories: well-differentiated, myxoid, dedifferentiated (DDLS) and pleomorphic [2]. Well-differentiated liposarcoma is now considered to be a precursor to dedifferentiated liposarcoma [3, 4]. Distinct molecular changes are found in each category. For instance, myxoid liposarcoma has been associated with a specific t12;16 translocation and expression of a FUS-CHOP fusion protein [5]. Alternative lengthening of telomeres and other telomere maintenance mechanisms are active in pleomorphic, dedifferentiated and other liposarcomas [68]. Amplification of cyclin-dependent kinase 4 (CDK4) and murine double minute 2 (MDM2) occurs in well-differentiated and dedifferentiated liposarcomas [9], as does down-regulation of PTEN [10]. Some dedifferentiated liposarcomas show amplification of STAT6 [2]. Activation of the PI3K/AKT pathway has been implicated in liposarcomas [1113] and loss of HIF-2a promotes liposarcoma growth [14]. Despite these major advances, important questions remain about the causes and the biology of liposarcomas.

Lipolysis, the pathway by which triglycerides are degraded, has been implicated in cancer, both as a source of fatty acids for tumor growth and as a mechanism of cancer-associated wasting [1517]. Lipolysis has been studied biochemically in greatest detail in adipose tissue. Two major adipose lipases are adipose triglyceride lipase (ATGL) and hormone-sensitive lipase (HSL). ATGL, encoded by the PNPLA2 gene on chromosome 11p15.5, is the main TG hydrolase of adipose tissue. It catalyzes the first step of lipolysis, i.e., the cleavage of triglyceride to diglyceride. Of note, deletion of PNPLA2 is reported in well differentiated liposarcoma and sarcoma [18, 19]. HSL, encoded by LIPE, can catalyze the second step of lipolysis, the cleavage of diglycerides to monoglycerides. Deletions of the chromosome 19p13 region containing LIPE are frequent in DDLS and correlate with poor prognosis [18]. Although PNPLA2 deletion has been reported in well differentiated liposarcoma and LIPE deletion, in DDLS with poor outcome [18], neither ATGL deficiency nor HSL deficiency reportedly causes liposarcoma in mice.

Increasing evidence suggests that cancer development may be driven by tissue-specific epistatic interactions due to mutations in multiple functionally-related genes [20, 21]. We generated three mouse lines with lipase deficiencies in adipose tissues, two single adipose knockout mice with deficiency of either ATGL or HSL plus double adipose knockout (DAKO) mice deficient in both ATGL and HSL. Strikingly, all DAKO mice developed a unique form of liposarcoma.


Combined deletion of Pnpla2 and Lipe in adipose tissue causes liposarcomas in mice

Complete deficiencies of ATGL and/or HSL in adipose tissue were confirmed in both white and brown adipose tissues by genomic Southern blotting (S1A Fig), Western blotting (S1B Fig) and real-time PCR of mRNA (S1C Fig). In isolated adipocytes, the rate of maximum beta adrenergic-stimulated lipolysis was markedly decreased in each single knockout line and was nearly unmeasurable in DAKO adipocytes (S2A Fig). Physiologically, DAKO mice were unable to maintain their blood glucose values with a normal postprandial fasting (S2B Fig), consistent with rapid depletion of carbohydrate reserves in the absence of lipid-derived energy from white adipose tissue (WAT). Cold tolerance, a measure of brown adipose tissue (BAT) function, was markedly reduced only in DAKO mice (S2C Fig).

ATGL and HSL catalyze the first two steps of lipolysis (Fig 1A). To test whether combined deficiency of ATGL plus HSL might promote liposarcoma (Fig 1A), 24 mice of each genotype were followed longitudinally. Both WAT and BAT were examined regularly by palpation and also directly at postmortem. No malignant tumors were found in WAT. In contrast, BAT of DAKO mice was hypertrophic and developed firm irregular tumors which invaded surrounding skin and tissue (Fig 1B). Histological evaluation of these tumors revealed liposarcoma (Fig 1B). Liposarcomas were detected between 11 and 14 months (24/24 mice, 100%) (Fig 1C). In contrast, neither single knockout line (ATGLAKO, HSLAKO) developed cancer (Fig 1C).

Characterization of lipase-deficient mice.
Fig. 1. Characterization of lipase-deficient mice.
(A) Experimental plan for generation of three different knockout mouse lines with lipolytic deficiency in adipose tissues. The positions of ATGL and HSL in the standard representation of the lipolytic pathway are shown below. (B) Liposarcoma in the left interscapular region of a 12-month-old DAKO mouse, gross and microscopic aspects. (C) Kaplan Meier plot displaying liposarcoma onset in the three knockout mouse lines and controls. (D) Principal component analysis (PCA) of whole transcriptome data, showing genotype-phenotype correlation for the early premalignant (3 months) and late (14 months) tumor stage. TG, triglycerides; DG, diglycerides; MG, monoglycerides; CTL, control.

The clear genotype-phenotype correlation was confirmed with expression microarrays. Principle component analysis revealed that by 3 months, DAKO mice showed a clearly distinct transcriptional profile when compared to either of the single knockouts or to normal controls (Fig 1D). DAKO tumors were strikingly different from all other samples including those from pretumor DAKO mice.

Premalignant DAKO BAT is morphologically abnormal and lacks UCP-1

Morphologically, by 3 months, BAT of the single and double knockout strains clearly differed from that of normal controls (Fig 2A). In all three knockout lines, nearly all brown adipocytes contained a single large lipid droplet in contrast to the multilocular morphology of normal BAT. The interscapular BAT depots of each knockout line were hypertrophic, particularly in DAKO mice (Fig 2B). Also, uncoupling protein-1 (UCP-1), which normally is expressed strongly and specifically in BAT, was undetectable in BAT of DAKO mice, in contrast with normal controls and with each single knockout line (Fig 2C).

Tissue specific response related to genotype.
Fig. 2. Tissue specific response related to genotype.
(A) Histology of WAT and BAT, showing controls, single (ATGLAKO, HSLAKO) and double knockout (DAKO) mice. (B) Mass of interscapular brown adipose tissue in different genotypes (3 months). (C) Western blot of the brown adipose-specific marker, UCP-1 (3 months). *, p < 0.05; **, p < 0.01; ***, p < 0.001.

Expression profiling of premalignant DAKO BAT

In order to study the mechanism of tumor initiation and progression over time, whole-transcriptome profiling was performed for all genotypes illustrated in Fig 1A. Differential expression analysis between normal controls and DAKO mice at three time points (3, 7 and 12 months) revealed that 21 genes were consistently differentially expressed at each point (Fig 3A). The greatest difference occurred at 3 months, when 499 genes were found to be differentially expressed between normal control and DAKO BAT. Gene set enrichment analysis (GSEA) using reactome pathways [22, 23] identified six pathways which were significantly enriched at an adjusted p value cutoff of p < 0.05 (Fig 3B). Down-regulation occurred for the gene sets of fatty acid, triacylglycerol and ketone body metabolism, the tricarboxylic acid (TCA) cycle and respiratory chain and genes of lipid metabolism. In contrast, genes involved in the immune response, particularly adaptive immunity, were up-regulated. Together, these results in premalignant DAKO BAT suggest a global down-regulation of oxidative energy metabolism and up-regulation of the immune response.

Early transcriptome analysis.
Fig. 3. Early transcriptome analysis.
(A) Venn diagram of intersections between differentially-expressed genes, comparing DAKO (n = 4) and WT mice (n = 4) at three different time points (3, 7 and 12 months). (B) Gene Set Enrichment Analysis for Reactome pathway, highlighting up- and down-regulated pathways for the differentially-expressed genes in BAT of 3-month-old mice. (C) Heatmap for all four genotypes, displaying the 21 genes that are consistently differentially-expressed between BAT of DAKO and control mice. Highlighted in red are Gpnmb and G0S2, the two most up- and down-regulated of these genes at all time points. (D) Boxplots indicating the expression levels of Gpnmb and G0s2 in the four mouse lines studied, at 3, 7 and 12 months.

Heat-map clustering of the 21 genes that are consistently differentially-expressed in DAKO and control BAT reveals a clear genotype-phenotype correlation, with marked differences between DAKO and normal control BAT and intermediate profiles in each single knockout line (Fig 3C). Among these genes, the two with the greatest up- and down-regulation in liposarcoma were glycoprotein nonmetastatic melanoma protein B (Gpnmb) and G0–G1 switch gene 2 (G0S2), respectively.

Gpnmb, a circulating glycoprotein, has been identified as a therapeutic target for several cancers [24]. Therapeutic suppression of Gpnmb is currently in phase 2 clinical trials for several types of cancer (ClinicalTrials.gov Identifier: NCT01997333, NCT01156753, NCT02713828, NCT02487979, NCT02302339, NCT02280785, NCT00921570, NCT02363283). Gpnmb is the top up-regulated gene in malignant and also premalignant DAKO BAT, compared to each of the other three genotypes (Fig 3C and 3D). Of note, plasma Gpnmb levels were also significantly (~3 fold) elevated in malignant and also premalignant DAKO mice compared to single knockout and normal controls, as early as 3 months of age, long before detectable liposarcoma development (S3 Fig).

Because Gpnmb is expressed in macrophages [25, 26], we explored whether macrophage expression of Gpnmb might account for its high level in DAKO BAT. Macrophage-related transcripts were found to be higher in HSLAKO BAT than in BAT of DAKO or in the other control lines (S4A Fig), but DAKO plasma shows higher levels of Gpnmb and DAKO BAT shows higher levels of Gpnmb expression than do corresponding samples from HSLAKO mice or from the other genotypes (S3 Fig). Therefore Gpnmb expression is high in DAKO BAT, and more so than anticipated from the observed level of macrophage infiltration. To explore whether adipocytes express Gpnmb, we performed in vitro differentiation of cultured brown adipocytes isolated from intrascapular DAKO BAT and also of cultured NIH 3T3-L1 white preadipocytes. Gpnmb was detectable after differentiation in both DAKO brown adipocytes and in 3T3-L1 cells (S4B Fig). Taking these observations together, it is likely that the high levels of Gpnmb in DAKO mice arise directly from adipocytes, although a fraction may derive from macrophages.

G0s2, the endogenous inhibitor of ATGL [27], is epigenetically silenced in several types of cancer [28, 29]. Interestingly, G0s2 was the transcript showing the greatest down-regulation in DAKO mice compared with normal controls and single knockout mice (Fig 3C and 3D).

Also, well-known liposarcoma markers such as MDM2 and p16/cdkn2 were detected. MDM2 is highly expressed in DAKO BAT at all ages studied (S5 Fig). The expression of two inhibitors of Cdk4, Cdkn2a and Cdkn2b, was significantly increased, as a late event, present only in DAKO liposarcomas (S5 Fig). Other cancer-associated genes such as Top2a, Ptn, Ptk7, Tnc and Mmp14 were highly expressed in DAKO mice (S5 Fig). Therefore, the expression of some known liposarcoma-related genes coincides with DAKO liposarcoma development.

High expression of angiogenesis-related genes in DAKO tumors

To further explore the late onset of tumor formation, we compared premalignant and malignant DAKO BAT. Morphologically, after 7 months, extensive necrosis is observed in BAT of DAKO mice (S6A Fig). Compared to late stage non-tumor BAT samples of DAKO mice, Gene Set Enrichment Analysis for Reactome pathway in DAKO liposarcomas revealed that the pathway with the greatest relative up-regulation of gene expression in DAKO liposarcomas was angiogenesis (S6B Fig), In this pathway, the highest up-regulation was seen for Esm1, encoding Endocan, expression of which has been reported to be increased up to 30-fold in human liposarcoma (SW872) cell lines [30, 31]. Other genes of angiogenesis that show upregulation in DAKO liposarcoma are Tnc, Ptn and Tnfrsf12a (S6B Fig). Together, these results suggest that angiogenesis is a key feature of DAKO liposarcoma development.

Differences in gene expression between WAT and BAT were reduced in DAKO mice

Morphologically, BAT samples from single knockout and DAKO mice show monolocular lipid droplets, as in WAT (Fig 2A). We compared expression data from WAT and BAT of three-month-old mice (S7 Fig). The number of genes that are differentially-expressed between brown and white adipose tissues is smaller in DAKO mice (114) than in normal controls (435), ATGLAKO (228) or HSLAKO mice (350) (S7A Fig). A subset of 48 genes is differentially expressed between BAT and WAT in all genotypes studied (S7B Fig). Based on these markers we performed time-course analyses. A subset of these BAT markers (Ucp1, Cidea, Cox8b) decreased progressively in DAKO BAT (S8 Fig). Of note, Pparγ, a critical regulator of identity and differentiation in white and brown adipocytes [32], is slightly but significantly down-regulated in 7- and 12-month-old DAKO BAT (S8 Fig).

Gene expression and copy number changes in DAKO liposarcoma and in human cancers

In order to explore whether some of the changes observed in DAKO liposarcoma might also be present in some human liposarcomas and other tumors, gene expression data from 58 DDLS patients from the TCGA data [33] were analyzed. Analysis of the orthologues of the 21 genes that were consistently differentially-expressed in DAKO BAT, revealed that LIPE, PNPLA2 and G0S2 are also among the five most down-regulated genes in human liposarcoma (Figs 3C and 4). Therefore, tissue specific inactivation of Pnpla2 and Lipe causes liposarcoma in mice, and down-regulation of PNPLA2 and LIPE expression occurs in human DDLS.

Comparison of gene expression in DAKO adipose tissue and in human dedifferentiated liposarcomas.
Fig. 4. Comparison of gene expression in DAKO adipose tissue and in human dedifferentiated liposarcomas.
Heatmap displaying gene expression patterns for all human-mouse orthologous genes in the publicly-available data set of 58 dedifferentiated liposarcomas. The vital status line indicates dead (red) or alive (green).

We also examined the expressions of GPNMB and PCOLCE in the 58 DDLS samples. No significant difference was found in the small number of samples, although there was high variability among the samples. Interestingly, among the 497 different types of cancer-related expression profiles in TCGA, GPNMB and PCOLCE were highly expressed in the two available pleomorphic liposarcoma samples (S9A Fig) and in sarcoma samples in general (S9B Fig).

In addition, copy number analysis for 265 patients with soft tissue sarcoma from the TCGA dataset was performed to test the hypothesis that PNPLA2 and LIPE might be tumor suppressors in some tissues. The loss of chromosome 11p15.5 harboring PNPLA2 was significantly associated (p < 0.01, chi-square) with soft tissue sarcomas (S10 Fig). In contrast, copy number changes of the region on chromosome 19q13 containing LIPE were not observed (S10 Fig).


Our study showed that the combined deficiency of the two major lipases of adipose tissue, ATGL and HSL, causes a unique form of liposarcoma with complete penetrance in mice. This clearly reveals a positive synergistic epistatic interaction between the Pnpla2 and Lipe genes, because neither single lipase deficiency causes liposarcoma.

There is growing interest in epistasis in cancer. Most cancers cannot be solely explained by single, highly-penetrant genetic changes [20, 21]. Large scale genetic analysis of human cancers has revealed co-enrichment (~5 to 16 fold) of genes whose products are known to physically or functionally interact in a tissue specific manner [20]. However, to date, few examples of functionally-proven biochemical interactions between members of the same metabolic pathway have been directly linked to epistatic mechanisms of cancer [21].

Each lipase of the lipolytic pathway shown in Fig 1 has clearly been related to cancer. Activation of monoacylglycerol lipase (MGL), which catalyzes the last step of lipolysis, promotes the progression of several cancers [34]. MGL has multiple physiological roles in addition to the cleavage of monoacylglycerols [35]; the mechanism by which MGL activation promotes cancer is unclear. ATGL is now linked to cancer in many ways. Deficiency of CGI-58, the coactivator of ATGL encoded by ABHD5, results in increased tumorigenesis and malignant transformation in mice and ABHD5 deletions occur in several human cancers [36]. G0S2, the endogenous inhibitor of ATGL, is epigenetically methylated in several cancers [28], has been considered to function as a tumor suppressor [37] and is markedly down-regulated in DAKO adipose tissue and DAKO liposarcoma. Al-Zougbi et al [19] recently reported lung cancer in 25% of ATGL-deficient mice that had prolonged survival. This group also studied ATGL copy number in a variety of human malignancies, and tested ATGL expression histochemically in human lung adenocarcinoma and squamous cell carcinoma, pancreatic intraductal neoplasia and adenocarcinoma and uterine leiomyosarcoma. They found lower ATGL expression in these cancers than in the corresponding normal tissues. The authors hypothesized that reduction of ATGL might reduce fatty acid availability and oxidation by the tumor cells, increasing their dependence on other energy sources like glucose and glutamine. A role for HSL in cancer development is less well documented than for MGL and ATGL, but a study of patients with well-differentiated and dedifferentiated liposarcoma showed that loss of the chromosome 19q13 region containing LIPE occurred frequently in DDLS, and that this was associated with decreased survival [18]. HSL deficiency is clearly essential for the development of DAKO liposarcomas. Taken together, these observations suggest that disruption of lipolytic pathway enzymes may occur widely in cancer and sometimes can be the primary and sufficient cause of cancer, as in DAKO liposarcoma.

ATGL and HSL are closely related metabolically. Functional overlap between these two lipases, as directly shown by lipolysis in isolated DAKO adipocytes and by intolerance to short fasting in DAKO mice, provides a potential mechanism for the epistatic phenomenon that occurs when both are deficient. During revision of this article we noted the description by Bi et al [38] of a mouse model of liposarcoma caused by adipocyte-specific overexpression of Notch1 intracellular protein (N1ICP). The authors conclude that deficiency of fatty acids within adipocytes and reduction of PPAR-γ signaling are critical events in the development of that form of liposarcoma. Intriguingly, although N1ICP-related liposarcomas occur in white adipocytes whereas DAKO liposarcoma occurs in brown adipose tissue, both the DAKO and N1ICP-related liposarcomas occur with complete penetrance, and in both, intracellular fatty acid pools and PPAR-γ signaling are predicted to be reduced. The discovery of the precise mechanism(s) by which combined deficiencies of ATGL and HSL cause liposarcoma will require further study. In DAKO BAT, major changes in gene expression are present by 3 months. This premalignant state may provide clues to the conditions under which DAKO liposarcomas develop.

In conclusion, our findings suggest that a synergistic epistatic interaction between two functionally related genes of the lipolysis pathway plays an important role in cancer development. The DAT mice described here provide a novel system for deciphering the early stages of development of liposarcoma and for testing prognostic markers and preventive treatment strategies.

Materials and methods

Ethics statement

Experiments were approved by Animal Facility Committee of CHU Sainte-Justine Hospital (protocol 620) according to the guidelines of the Canadian Council on Animal Care (http://www.ccac.ca/en_/).

Creation of ATGLAKO, HSLAKO and DAKO mice

Four mouse strains were used: normal controls, mice deficient in adipose tissue for ATGL (ATGLAKO), HSL (HSLAKO) or both (DAKO). The creation of ATGLAKO mice was described previously [39]. HSLAKO mice were homozygous for a floxed allele of Lipe and harbored a transgene, from which Cre is expressed from a Fabp4 promoter, which is active in adipose tissues. Further details about the creation of HSLAKO mice are provided in S1 Text. DAKO mice were created by breeding of the two single knockout lines. Mice expressing only the Fabp4-Cre transgene were used as controls. Each targeted allele had been transferred for at least 8 generations to a C57BL/6J background.

Mouse breeding and handling

Mice were raised with a light-dark cycle with lights on at 7am and off at 7pm, and free access to water and food (no. 2019, Teklad Global Rodent, 9% fat; Harlan Laboratories, Madison, WI). Mice were housed at 21° C unless otherwise specified. For sacrifice, mice were anesthetized with sodium pentobarbital (Somnotol; MTC Pharmaceuticals, Hamilton, Ontario, Canada). After cardiac puncture, organs were rapidly removed, weighed, and frozen. Experiments were approved by the Canadian Council on Animal Care-accredited animal facility committee of CHU Sainte-Justine Hospital.

Western blotting & plasma GPNMB measurement

Western blotting was performed as described [40]. The following antibodies were used: ATGL (#10006409–1, Cayman Chemical), HSL [40], UCP-1 (# RB-10599-P0, Lab vision corporation) and GPNMB (#AF2330-SP, R&D Systems). Plasma levels of GPNMB soluble protein were measured using the commercially available ELISA kit from ORIGENE (EA100800).

Lipolysis in isolated adipocytes

Isolation of adipocytes and measurement of lipolysis were performed as described [39].

Fasting tests in lipase-deficient mice

Six mice of each genotype, aged 3 months, were fasted starting at 7am. Plasma glucose was measured at 12 noon.


Six mice of each genotype, aged 2.5 months, were placed at 4°C. Five hours after food removal, body temperature was measured using a rectal thermometer (BAT-10 Thermometer; Physitemp Instruments, Inc., Clifton, NJ).

Whole-transcriptome analysis

Interscapular BAT from 3, 7 and 12-month old mice for tumor and non-tumor BAT and perigonadal WAT was rapidly obtained from anesthetized 5h fasted mice, then was snap frozen in liquid nitrogen. RNA was extracted using the RNeasy Lipid Tissue Mini Kit (QIAGEN, Cat. No.74804). Illumina mouse Ref-8 v2.0 Expression BeadChips were used for whole-transcriptome profiling. http://support.illumina.com/array/array_kits/mouseref-8_v2_expression_beadchip_kit/questions.html. Raw expression data was obtained as output from the BeadArrray system via the Beadstudio software interface. Quality control and preprocessing of the raw data was performed using the R & Bioconductor software suite. Linear models for microarray data (limma) were used to identify significant expression differences of genes between biological replicates. Probe identifiers were mapped to HGNC symbols and Refseq identifiers for subsequent pathway and Gene Set Enrichment Analysis. For differential gene expression analysis, all possible comparisons between the three genotypes and the normal controls were designed. Genes designated as having significant differential expression had a p-value threshold of 0.01 after correction for multiple testing by the false discovery rate method (FDR) [41]. Genes were sorted according to the empirical Bayes log odds of differential expression (B value) and the adjusted p values. Hierarchical clustering of genes based on Pearson’s correlation coefficients revealed groups of genes with similar expression patterns across genotypes. Groups of differentially-expressed genes were displayed as heatmaps or Venn diagrams. Gene Set Enrichment analysis was performed using the standard parameters for the Reactome pathway hallmark data set (www.broadinstitute.org/gsea/msigdb/collections) [22, 23].

Expression and CNV data from sarcoma patients

Gene expression and CNV data from 265 soft sarcoma patients (SARC cohort) and 972 breast cancer patients was obtained from the TCGA dataset [33] via the firebrowse data repository (www.firebrowse.org). Expression data was analyzed for a subset of 58 patients with dedifferentiated liposarcoma based on human-mouse ortholog comparisons in the DAKO mouse model. CNV data was processed using GISTIC2 [41] and estimated values were assigned to the following thresholds of -2, -1, 0, 1, 2, representing homozygous deletion, single copy deletion, diploid normal copy, low-level copy number amplification, or high-level copy number amplification, respectively. The gene visible data source (www.genevisible.com), containing expression profiles of 497 cancers, and the TCGA dataset were used to visualize cancer-specific expression profiles for PCOLCE and GPNMB. CNV data from 265 soft tissue sarcoma patients from the TCGA resource was analyzed to identify gene deletions and losses associated with chromosomal regions harboring PNPLA2 and LIPE.

Culture and differentiation of interscapular brown adipocytes and NIH 3T3-L1 cells

Interscapular brown adipocytes from 3-day-old mice were isolated and cultured, then differentiation was induced as described [42]. NIH3T3-L1 cells were cultured and differentiation was induced as described [43].

Web resources

Reactome pathway hallmark data set; www.broadinstitute.org/gsea/msigdb

Broad portal for TCGA cancer data set; www.firebrowse.org

Genevisible expression data set; www.genevisible.com

Supporting Information

Attachment 1


1. Morrison BA. Soft tissue sarcomas of the extremities. Proc (Bayl Univ Med Cent). 2003;16(3):285–90. PubMed Central PMCID: PMCPMC1200782.

2. Doyle LA, Tao D, Marino-Enriquez A. STAT6 is amplified in a subset of dedifferentiated liposarcoma. Mod Pathol. 2014;27(9):1231–7. doi: 10.1038/modpathol.2013.247 24457460

3. Thway K, Jones RL, Noujaim J, Zaidi S, Miah AB, Fisher C. Dedifferentiated Liposarcoma: Updates on Morphology, Genetics, and Therapeutic Strategies. Adv Anat Pathol. 2016;23(1):30–40. doi: 10.1097/PAP.0000000000000101 26645460

4. Iwasa Y, Nakashima Y. Dedifferentiated liposarcoma with lipoma-like well-differentiated liposarcoma: clinicopathological study of 30 cases, with particular attention to the comingling pattern of well- and dedifferentiated components: a proposal for regrouping of the present subclassification of well-differentiated liposarcoma and dedifferentiated liposarcoma. Int J Surg Pathol. 2013;21(1):15–21. doi: 10.1177/1066896912449040 22674916

5. Knight JC, Renwick PJ, Dal Cin P, Van den Berghe H, Fletcher CD. Translocation t(12;16)(q13;p11) in myxoid liposarcoma and round cell liposarcoma: molecular and cytogenetic analysis. Cancer Res. 1995;55(1):24–7. 7805034

6. Jeyapalan JN, Mendez-Bermudez A, Zaffaroni N, Dubrova YE, Royle NJ. Evidence for alternative lengthening of telomeres in liposarcomas in the absence of ALT-associated PML bodies. Int J Cancer. 2008;122(11):2414–21. doi: 10.1002/ijc.23412 18311780

7. Lee JC, Jeng YM, Liau JY, Tsai JH, Hsu HH, Yang CY. Alternative lengthening of telomeres and loss of ATRX are frequent events in pleomorphic and dedifferentiated liposarcomas. Mod Pathol. 2015;28(8):1064–73. doi: 10.1038/modpathol.2015.67 26022452

8. Johnson JE, Gettings EJ, Schwalm J, Pei J, Testa JR, Litwin S, et al. Whole-genome profiling in liposarcomas reveals genetic alterations common to specific telomere maintenance mechanisms. Cancer Res. 2007;67(19):9221–8. doi: 10.1158/0008-5472.CAN-07-1133 17909028

9. Abbas Manji G, Singer S, Koff A, Schwartz GK. Application of molecular biology to individualize therapy for patients with liposarcoma. Am Soc Clin Oncol Educ Book. 2015:213–8. doi: 10.14694/EdBook_AM.2015.35.213 25993159

10. Smith KB, Tran LM, Tam BM, Shurell EM, Li Y, Braas D, et al. Novel dedifferentiated liposarcoma xenograft models reveal PTEN down-regulation as a malignant signature and response to PI3K pathway inhibition. Am J Pathol. 2013;182(4):1400–11. PubMed Central PMCID: PMCPMC3620414. doi: 10.1016/j.ajpath.2013.01.002 23416162

11. Guo S, Lopez-Marquez H, Fan KC, Choy E, Cote G, Harmon D, et al. Synergistic effects of targeted PI3K signaling inhibition and chemotherapy in liposarcoma. PLoS One. 2014;9(4):e93996. PubMed Central PMCID: PMCPMC3973642. doi: 10.1371/journal.pone.0093996 24695632

12. Gutierrez A, Snyder EL, Marino-Enriquez A, Zhang YX, Sioletic S, Kozakewich E, et al. Aberrant AKT activation drives well-differentiated liposarcoma. Proc Natl Acad Sci U S A. 2011;108(39):16386–91. PubMed Central PMCID: PMCPMC3182699. doi: 10.1073/pnas.1106127108 21930930

13. Demicco EG, Torres KE, Ghadimi MP, Colombo C, Bolshakov S, Hoffman A, et al. Involvement of the PI3K/Akt pathway in myxoid/round cell liposarcoma. Mod Pathol. 2012;25(2):212–21. PubMed Central PMCID: PMCPMC3266971. doi: 10.1038/modpathol.2011.148 22020193

14. Nakazawa MS, Eisinger-Mathason TS, Sadri N, Ochocki JD, Gade TP, Amin RK, et al. Epigenetic re-expression of HIF-2alpha suppresses soft tissue sarcoma growth. Nat Commun. 2016;7:10539. PubMed Central PMCID: PMCPMC4742834. doi: 10.1038/ncomms10539 26837714

15. 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–9. PubMed Central PMCID: PMCPMC4002264. doi: 10.1016/j.plipres.2013.08.005 24001676

16. Arner P, Langin D. Lipolysis in lipid turnover, cancer cachexia, and obesity-induced insulin resistance. Trends Endocrinol Metab. 2014;25(5):255–62. doi: 10.1016/j.tem.2014.03.002 24731595

17. Nieman KM, Kenny HA, Penicka CV, Ladanyi A, Buell-Gutbrod R, Zillhardt MR, et al. Adipocytes promote ovarian cancer metastasis and provide energy for rapid tumor growth. Nat Med. 2011;17(11):1498–503. PubMed Central PMCID: PMCPMC4157349. doi: 10.1038/nm.2492 22037646

18. Crago AM, Socci ND, DeCarolis P, O'Connor R, Taylor BS, Qin LX, et al. Copy number losses define subgroups of dedifferentiated liposarcoma with poor prognosis and genomic instability. Clin Cancer Res. 2012;18(5):1334–40. PubMed Central PMCID: PMCPMC3294014. doi: 10.1158/1078-0432.CCR-11-2820 22241790

19. Al-Zoughbi W, Pichler M, Gorkiewicz G, Guertl-Lackner B, Haybaeck J, Jahn SW, et al. Loss of adipose triglyceride lipase is associated with human cancer and induces mouse pulmonary neoplasia. Oncotarget. 2016;7(23):33832–40. PubMed Central PMCID: PMCPMC5085122. doi: 10.18632/oncotarget.9418 27213586

20. Park S, Lehner B. Cancer type-dependent genetic interactions between cancer driver alterations indicate plasticity of epistasis across cell types. Mol Syst Biol. 2015;11(7):824. PubMed Central PMCID: PMCPMC4547852. doi: 10.15252/msb.20156102 26227665

21. Wang X, Fu AQ, McNerney ME, White KP. Widespread genetic epistasis among cancer genes. Nat Commun. 2014;5:4828. doi: 10.1038/ncomms5828 25407795

22. Liberzon A, Birger C, Thorvaldsdottir H, Ghandi M, Mesirov JP, Tamayo P. The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Syst. 2015;1(6):417–25. PubMed Central PMCID: PMCPMC4707969. doi: 10.1016/j.cels.2015.12.004 26771021

23. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50. PubMed Central PMCID: PMCPMC1239896. doi: 10.1073/pnas.0506580102 16199517

24. Maric G, Rose AA, Annis MG, Siegel PM. Glycoprotein non-metastatic b (GPNMB): A metastatic mediator and emerging therapeutic target in cancer. Onco Targets Ther. 2013;6:839–52. PubMed Central PMCID: PMCPMC3711880. doi: 10.2147/OTT.S44906 23874106

25. Ripoll VM, Irvine KM, Ravasi T, Sweet MJ, Hume DA. Gpnmb is induced in macrophages by IFN-gamma and lipopolysaccharide and acts as a feedback regulator of proinflammatory responses. J Immunol. 2007;178(10):6557–66. 17475886

26. Kumagai K, Tabu K, Sasaki F, Takami Y, Morinaga Y, Mawatari S, et al. Glycoprotein Nonmetastatic Melanoma B (Gpnmb)-Positive Macrophages Contribute to the Balance between Fibrosis and Fibrolysis during the Repair of Acute Liver Injury in Mice. PLoS One. 2015;10(11):e0143413. PubMed Central PMCID: PMCPMC4657955. doi: 10.1371/journal.pone.0143413 26599547

27. Yang X, Lu X, Lombes M, Rha GB, Chi YI, Guerin TM, et al. The G(0)/G(1) switch gene 2 regulates adipose lipolysis through association with adipose triglyceride lipase. Cell Metab. 2010;11(3):194–205. PubMed Central PMCID: PMCPMC3658843. doi: 10.1016/j.cmet.2010.02.003 20197052

28. Kusakabe M, Kutomi T, Watanabe K, Emoto N, Aki N, Kage H, et al. Identification of G0S2 as a gene frequently methylated in squamous lung cancer by combination of in silico and experimental approaches. Int J Cancer. 2010;126(8):1895–902. doi: 10.1002/ijc.24947 19816938

29. Barreau O, Assie G, Wilmot-Roussel H, Ragazzon B, Baudry C, Perlemoine K, et al. Identification of a CpG island methylator phenotype in adrenocortical carcinomas. J Clin Endocrinol Metab. 2013;98(1):E174–84. doi: 10.1210/jc.2012-2993 23093492

30. Almog N, Ma L, Raychowdhury R, Schwager C, Erber R, Short S, et al. Transcriptional switch of dormant tumors to fast-growing angiogenic phenotype. Cancer Res. 2009;69(3):836–44. doi: 10.1158/0008-5472.CAN-08-2590 19176381

31. Delehedde M, Devenyns L, Maurage CA, Vives RR. Endocan in cancers: a lesson from a circulating dermatan sulfate proteoglycan. Int J Cell Biol. 2013;2013:705027. PubMed Central PMCID: PMCPMC3625564. doi: 10.1155/2013/705027 23606845

32. Rosen ED, Spiegelman BM. Molecular regulation of adipogenesis. Annu Rev Cell Dev Biol. 2000;16:145–71. doi: 10.1146/annurev.cellbio.16.1.145 11031233

33. Cancer Genome Atlas Research N, Weinstein JN, Collisson EA, Mills GB, Shaw KR, Ozenberger BA, et al. The Cancer Genome Atlas Pan-Cancer analysis project. Nat Genet. 2013;45(10):1113–20. PubMed Central PMCID: PMCPMC3919969. doi: 10.1038/ng.2764 24071849

34. 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. PubMed Central PMCID: PMCPMC2885975. doi: 10.1016/j.cell.2009.11.027 20079333

35. Shimizu T, Tanaka K, Shimizu S, Higashi Y, Yawata T, Nakamura K, et al. Possible inhibitory role of endogenous 2-arachidonoylglycerol as an endocannabinoid in (+/-)-epibatidine-induced activation of central adrenomedullary outflow in the rat. Neuropharmacology. 2015;95:278–89. doi: 10.1016/j.neuropharm.2015.03.034 25882827

36. Ou J, Miao H, Ma Y, Guo F, Deng J, Wei X, et al. Loss of abhd5 promotes colorectal tumor development and progression by inducing aerobic glycolysis and epithelial-mesenchymal transition. Cell Rep. 2014;9(5):1798–811. PubMed Central PMCID: PMCPMC4268306. doi: 10.1016/j.celrep.2014.11.016 25482557

37. Yim CY, Sekula DJ, Hever-Jardine MP, Liu X, Warzecha JM, Tam J, et al. G0S2 Suppresses Oncogenic Transformation by Repressing a MYC-Regulated Transcriptional Program. Cancer Res. 2016;76(5):1204–13. PubMed Central PMCID: PMCPMC4775337. doi: 10.1158/0008-5472.CAN-15-2265 26837760

38. Bi P, Yue F, Karki A, Castro B, Wirbisky SE, Wang C, et al. Notch activation drives adipocyte dedifferentiation and tumorigenic transformation in mice. J Exp Med. 2016;213(10):2019–37. PubMed Central PMCID: PMCPMC5030803. doi: 10.1084/jem.20160157 27573812

39. Wu JW, Wang SP, Casavant S, Moreau A, Yang GS, Mitchell GA. Fasting energy homeostasis in mice with adipose deficiency of desnutrin/adipose triglyceride lipase. Endocrinology. 2012;153(5):2198–207. doi: 10.1210/en.2011-1518 22374972

40. Fortier M, Soni K, Laurin N, Wang SP, Mauriege P, Jirik FR, et al. Human hormone-sensitive lipase (HSL): expression in white fat corrects the white adipose phenotype of HSL-deficient mice. J Lipid Res. 2005;46(9):1860–7. doi: 10.1194/jlr.M500081-JLR200 15961788

41. Du P, Kibbe WA, Lin SM. lumi: a pipeline for processing Illumina microarray. Bioinformatics. 2008;24(13):1547–8. doi: 10.1093/bioinformatics/btn224 18467348

42. Klein J, Fasshauer M, Ito M, Lowell BB, Benito M, Kahn CR. beta(3)-adrenergic stimulation differentially inhibits insulin signaling and decreases insulin-induced glucose uptake in brown adipocytes. J Biol Chem. 1999;274(49):34795–802. 10574950

43. Song Z, Cheng J, Yang H, Li Y, Gao Q, Shi X, et al. Differentiation of 3T3-L1 preadipocytes is inhibited under a modified ceiling culture. Cell Biol Int. 2015;39(5):638–45. doi: 10.1002/cbin.10428 25572439

Genetika Reprodukční medicína
Zapomenuté heslo

Nemáte účet?  Registrujte se

Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.


Nemáte účet?  Registrujte se