Next Article in Journal
Integrated Framework of the Immune-Defense Transcriptional Signatures in the Arabidopsis Shoot Apical Meristem
Next Article in Special Issue
Role of Adipose Tissue-Derived Autotaxin, Lysophosphatidate Signaling, and Inflammation in the Progression and Treatment of Breast Cancer
Previous Article in Journal
Premature Activation of Immune Transcription Programs in Autoimmune-Predisposed Mouse Embryonic Stem Cells and Blastocysts
Previous Article in Special Issue
Aspirin Modifies Inflammatory Mediators and Metabolomic Profiles and Contributes to the Suppression of Obesity-Associated Breast Cancer Cell Growth
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 21
Issue 16
10.3390/ijms21165744
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Intratumoral Adipocyte-High Breast Cancer Enrich for Metastatic and Inflammation-Related Pathways but Associated with Less Cancer Cell Proliferation

by
Yoshihisa Tokumaru
1,2,
Masanori Oshi
1,3,
Eriko Katsuta
1,
Li Yan
4,
Jing Li Huang
1,
Masayuki Nagahashi
5,
Nobuhisa Matsuhashi
2,
Manabu Futamura
2,
Kazuhiro Yoshida
2 and
Kazuaki Takabe
1,3,5,6,7,8,*
1
Breast Surgery, Department of Surgical Oncology, Roswell Park Comprehensive Cancer Center, Buffalo, NY 14263, USA
2
Department of Surgical Oncology, Graduate School of Medicine, Gifu University, 1-1 Yanagido, Gifu 501-1194, Japan
3
Department of Gastroenterological Surgery, Yokohama City University Graduate School of Medicine, Yokohama 236-0004, Japan
4
Department of Biostatistics & Bioinformatics, Roswell Park Comprehensive Cancer Center, Buffalo, NY 14263, USA
5
Department of Surgery, Niigata University Graduate School of Medical and Dental Sciences, Niigata 951-8510, Japan
6
Department of Surgery, University at Buffalo Jacobs School of Medicine and Biomedical Sciences, The State University of New York, Buffalo, NY 14263, USA
7
Department of Breast Oncology and Surgery, Tokyo Medical University, 6-7-1 Nishishinjuku, Shinjuku, Tokyo 160-8402, Japan
8
Department of Breast Surgery, Fukushima Medical University School of Medicine, Fukushima 960-1295, Japan
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(16), 5744; https://doi.org/10.3390/ijms21165744
Submission received: 30 July 2020 / Revised: 8 August 2020 / Accepted: 10 August 2020 / Published: 11 August 2020
(This article belongs to the Special Issue Unraveling the Involvement of the Adipose Tissue in Breast Cancer Progression)
Browse Figures
Review Reports Versions Notes

Abstract

:
Cancer-associated adipocytes are known to cause inflammation, leading to cancer progression and metastasis. The clinicopathological and transcriptomic data from 2256 patients with breast cancer were obtained based on three cohorts: The Cancer Genome Atlas (TCGA), GSE25066, and a study by Yau et al. For the current study, we defined the adipocyte, which is calculated by utilizing a computational algorithm, xCell, as “intratumoral adipocyte”. These intratumoral adipocytes appropriately reflected mature adipocytes in a bulk tumor. The amount of intratumoral adipocytes demonstrated no relationship with survival. Intratumoral adipocyte-high tumors significantly enriched for metastasis and inflammation-related gene sets and are associated with a favorable tumor immune microenvironment, especially in the ER+/HER2- subtype. On the other hand, intratumoral adipocyte-low tumors significantly enriched for cell cycle and cell proliferation-related gene sets. Correspondingly, intratumoral adipocyte-low tumors are associated with advanced pathological grades and inversely correlated with MKI67 expression. In conclusion, a high amount of intratumoral adipocytes in breast cancer was associated with inflammation, metastatic pathways, cancer stemness, and favorable tumor immune microenvironment. However, a low amount of adipocytes was associated with a highly proliferative tumor in ER-positive breast cancer. This cancer biology may explain the reason why patient survival did not differ by the amount of adipocytes.

1. Introduction

Breast cancer is one of the leading causes of cancer-related deaths in developed countries, including the United States [1,2]. Obesity is one of the known risk factors linked to the development and progression of breast cancer, particularly in the estrogen receptor-positive (ER+) subtype [3,4]. Excess adipose tissue secretes inflammatory cytokines and mediators [3,5]. Our group has indeed previously reported that a lipid mediator, sphingosine-1-phosphate, is more abundantly produced in the obese body and links obesity, inflammation, breast cancer metastasis [6], and treatment resistance [7]. In addition to adipose tissue, the role of adipocytes in the tumor microenvironment has lately drawn considerable attention [8].
Numerous investigators, including our group, have reported that not only cancer cells but also the tumor microenvironment plays a critical role in cancer progression and metastasis, as well as therapeutic responses [9,10]. Adipocytes, together with vascular endothelial cells [6,11,12], fibroblasts [13,14], and immune cells [15,16,17,18,19] constitute the tumor microenvironment [20,21,22]. Cancer-associated adipocytes interact with cancer cells and secrete the inflammatory cytokines, such as IL-6 and TNF-α, which contribute to pro-cancer inflammation that is known to aggravate cancer progression [23,24,25,26,27]. IL-6 promotes angiogenesis, invasion, and metastasis by activating the JAK/STAT3 pathway [28,29]. TNF-α promotes the invasion and metastasis of breast cancer by activating epithelial–mesenchymal transition (EMT) signaling [30,31,32]. In addition, cancer-associated adipocytes secrete chemokines such as Chemokine (C–C motif) ligand 2 (CCL2), which activates NOTCH1 signaling pathway and induces the stemness of the cancer cells [22,33]. Leptin, which is also produced by cancer-associated adipocytes, collaborates with IL-1, and promotes angiogenesis by upregulating expression of VEGF and VEGF receptors [22,34,35,36].
Although pre-clinical studies have shown how adipocytes may affect the progression of cancer, little is known about the clinical relevance of the amount of adipocytes in breast cancer. Recently, our group has been assessing the clinical relevance of the tumor microenvironment utilizing in silico approaches [10,14,17,18,37,38,39,40,41].
We demonstrated that factors affecting favorable tumor immunity may lead to better clinical outcomes. For instance, we demonstrated that pancreatic cancer with a higher amount of fibroblasts is associated with anti-cancer immunity and a higher rate of complete resection without residual tumor (R0) [14]. In addition, we reported that breast cancers with high expression of microRNA-143 were associated with a favorable tumor immune microenvironment and led to better survival [18]. In a similar fashion, we estimated the amount of adipocytes utilizing xCell, a bioinformatics algorithm that defined adipocytes by the expression pattern of 21 genes, as reported by Aran D et al. (Supplementary Table S1) [42].
Here, we hypothesized that breast cancer with a high amount of intratumoral adipocytes is associated with inflammation, metastasis, and worse survival. The aim of this study was to test this hypothesis utilizing in silico analysis of publicly available large databases.

2. Results

2.1. Transcriptionally Defined Adipocytes Significantly Correlated with Markers of Mature adipocytes

First, we assessed whether transcriptionally defined adipocytes appropriately reflected the characteristics of differentiated adipocytes by analyzing the correlation with their key markers: Adiponectin, leptin, lipoprotein lipase, and perilipin 1 (gene names are ADIPOQ, LEP, LPL, and PLIN1, respectively). Based on the TCGA cohort, transcriptionally defined adipocytes were strongly correlated with adiponectin, leptin, lipoprotein lipase, and perilipin 1 (Figure 1A; r = 0.867, r = 0.794, r = 0.672, and r = 0.875, respectively). We used two other cohorts to validate these results. Adipocytes were either moderately or strongly correlated with adiponectin, leptin, lipoprotein lipase, and perilipin 1 in the GSE25066 cohort (Figure 1B; r = 0.718, r = 0.472, r = 0.479, and r = 0.718, respectively) and the cohort of Yau et al. (Figure 1C; r = 0.758, r = 0.477, r = 0.615). No data for perilipin 1 was available in the latter cohort due to a lack of probe for PLIN1. These results implicated that transcriptionally defined adipocytes reasonably reflected mature adipocytes in a bulk tumor.

2.2. Breast Cancer Has High Level of Intratumoral Adipocytes Compared to the Other Cancers

Our methodology to analyze gene signatures to quantify adipocytes allowed us to measure the amount of adipocytes in any gene expression dataset, including the entire TCGA cohort. Since samples cataloged in TCGA were collected from a cancer cell-rich portion of the bulk tumors, we specifically termed our findings as “intratumoral adipocytes” to differentiate from infiltrating adipocytes in the peritumoral region. Utilizing our method, we compared the amount of intratumoral adipocyte among the different types of cancer in TCGA. Breast cancer (BRCA) possessed the second highest amount of adipocytes among the various cancers in TCGA, trailing only to hepatocellular carcinoma (LIHC) (Figure 1D). The other cancers with a relatively high amount of adipocytes included brain cancer (LGG, GBM), kidney cancer (KIRC, KICH, KIRP), mesothelioma (MESO), and cholangiocarcinoma (CHOL). On the contrary, almost no adipocytes were detected in melanoma (SLCM, UVM) or leukemia (LAML).

2.3. Intratumoral Adipocyte-High Tumors Were Not Associated with Worse Survival in Any of the Subtypes, nor with Pathological Complete Response (pCR) to Neoadjuvant Chemotherapy (NAC)

Given the detrimental role of cancer-associated adipocytes [26], we investigated whether the amount of intratumoral adipocyte was associated with the clinical outcomes in breast cancer. We divided the patients into high and low groups by the amount of adipocytes using a median cutoff. There was no statistically significant association between the amount of intratumoral adipocyte and patient survival, although intratumoral adipocyte-high tumors demonstrated a tendency for better overall survival (OS) in the whole TCGA cohort and the estrogen receptor positive/HER2 negative (ER+/HER2-) subtype (Figure 2; p = 0.085 and p = 0.088, respectively). There was no association between intratumoral adipocyte-high tumors and survival in either the HER2 positive (HER2+) or the triple negative (TNBC) subtype (Figure 2).
A previous in vitro study reported that the co-culture of adipocyte with human and murine cell lines promoted the resistance to doxorubicin, paclitaxel, and 5-fluorouracil (5-FU) independently from the subtypes [43]. Thus, we expected that intratumoral adipocyte-high tumors were associated with decreased pathological response (pCR) rate. We analyzed three independent neoadjuvant chemotherapy (NAC) cohorts that underwent different chemotherapy regimens; GSE25066 (n = 508; taxane and anthracycline) [44], GSE32646 (n = 115; 5-FU, epirubicin, cyclophosphamide, and paclitaxel) [45], and GSE20194 (n = 278; paclitaxel, 5-FU, cyclophosphamide and doxorubicin) [46]. Unexpectedly, intratumoral adipocyte-high tumors were not associated with decreased pCR rate in both ER+/HER2- and TNBC subtypes in any of NAC cohorts we analyzed (Supplementary Figure S1).

2.4. Intratumoral Adipocyte-High Tumors were not Associated with Advanced Cancer Staging

Cancer-associated adipocytes are known to promote tumor growth, invasion, and migration of breast cancer cells in vitro [25]. To this end, we hypothesized that tumors that were high in intratumoral adipocytes would be associated with the advanced American Joint Commission on Cancer (AJCC) clinical cancer staging in the TCGA cohort. Contrary to our expectation, less aggressive stage 1 tumors tended to associate with increased intratumoral adipocytes compared to tumors of other stages in the whole cohort and in any of the subtypes; however, none of these associations were statistically significant (Supplementary Figure S2; all p > 0.05).

2.5. Intratumoral Adipocyte-High Tumors Enriched for Inflammation, Metastasis, and Immune Response-Related Gene Sets

In order to examine whether intratumoral adipocytes demonstrated similar biological functions as previously reported in vitro [25], a gene set enrichment assay (GSEA) was performed on the TCGA cohort. In the whole cohort and in both ER+/HER2- and TNBC subtypes, intratumoral adipocyte-high tumors enriched for inflammation-related gene sets: TNF-α signaling via NFκB, IL-6/JAK/STAT3 signaling, and inflammatory response (Figure 3A, Supplementary Figures S3A and S4A). This was concordant with previous reports [25,26]. The aforementioned gene sets were not enriched in the GSE25066 cohort, which had fewer patients compared with TCGA. We also found intratumoral adipocyte-high tumors enriched for epithelial–mesenchymal transition (EMT), TGFβ signaling, angiogenesis, WNTβ catenin signaling, and Hedgehog signaling in the whole cohort and the ER+/HER2 subtype in TCGA (Figure 3B, Supplementary Figures S3B and S4B). These gene sets reflected on mechanisms of metastasis and cancer stemness. These were also enriched in intratumoral adipocyte-high tumors in the GSE25066 cohort. Only in the ER+/HER2 subtype was intratumoral adipocytes-high tumors enriched for Notch signaling. In the TNBC subtype, only TGFβ signaling was enriched among the metastasis and cancer stemness-related gene sets (Supplementary Figure S4B). On the other hand, intratumoral adipocyte-high tumors enriched for gene sets related to immune response such as IL2 signaling and IFN-γ only in ER+/HER2- subtype in TCGA cohort (Figure 3C, Supplementary Figures S3C and S4C). In the GSE25066 cohort, only IL2 STAT5 signaling was enriched among the immune response gene sets. These results suggest that intratumoral adipocyte-high tumors were associated with inflammation, metastasis, and immune-related gene sets, especially in the ER+/HER2- subtype.

2.6. Intratumoral Adipocyte-High Tumors Were Associated with Favorable Tumor Immune Microenvironment in ER+/HER2- Subtype

Given the data from the gene set enrichment analysis, we further explored the infiltration of immune cells within the tumor microenvironment using the whole TCGA breast cancer cohort. Regarding favorable immune cells, intratumoral adipocyte high tumors were associated with a high proportion of dendritic cells (DC), but with less with type 1 helper T cells (Th1) and M1 macrophages (M1) in the whole cohort (Figure 4A; all p < 0.001). In terms of unfavorable immune cells, less type 2 helper T cells (Th2), regulatory T cells (Treg), and M2 macrophages (M2) were found in intratumoral adipocyte-high tumors, compared to adipocyte-low tumors (Figure 4B; all p < 0.001). A similar pattern of results was found in the ER+/HER2- subtype, whereas only Th1 and Treg were associated with intratumoral adipocyte-high tumors in the TNBC subtype (Figure 4A,B). The cytolytic activity score (CYT) reflects the cytolytic activity of all the immune cells in the bulk tumor. High CYT was associated with intratumoral adipocyte-high tumors in the whole cohort and in the ER+/HER2- subtype, but not in TNBC subtype (Figure 4C). These results indicated that adipocyte-high tumors were associated with overall favorable tumor immune microenvironment in ER+/HER2- subtype, but not in TNBC.

2.7. Intratumoral Adipocyte-Low Tumor Enriched Cell Cycle- and Cell Proliferation-Related Gene Sets

We further investigated gene sets that may be associated with intratumoral adipocyte-low tumors, given that tumors with high adipocytes were enriched for metastasis and inflammation-related gene sets but did not translate to a survival difference. We performed GSEA on intratumoral adipocyte-low tumors in the TCGA and GSE25066 cohorts. We found that intratumoral adipocyte-low tumors significantly enriched for cell cycle and cell proliferation gene sets such as E2F targets, G2M checkpoint, MYC targets V1, and MYC targets V2 in the whole breast cancer cohort of TCGA. Similar results were found in patients with the ER+/HER2 subtype in the GSE25066 cohort (Figure 5, Supplementary Figures S5 and S6). These results demonstrate that intratumoral adipocyte-low tumors are associated with increased cell cycle and cell proliferation.

2.8. Adipocyte Low Tumor Correlated with High MKI67 Expression, and Was Associated with Advanced Nottingham Pathological Grade

Given that intratumoral adipocyte-low tumor strongly enriched most of the cell proliferation gene sets, we hypothesized that intratumoral adipocyte low tumors were a clinically highly proliferative cancer. We examined the association between the amount of intratumoral adipocyte and MKI67 expression, the most commonly used marker of cancer cell proliferation in the clinical setting [37,47,48]. As we expected, intratumoral adipocyte low tumors demonstrated significantly higher MKI67 expression consistently across all three cohorts TCGA, GSE25066, and the cohort from Yau et al. (Figure 6A, p < 0.001, p < 0.001, and p < 0.001, respectively). Furthermore, the amount of intratumoral adipocytes was weakly inversely correlated with MKI67 expression in TCGA (Figure 6B; r = −0.235, p < 0.01). This result was validated with two other cohorts, GSE25066 and the cohort reported by Yau et al. (Figure 6B; r = −2.09 and r = −0.322, respectively).
In clinical practice, cell proliferation was assessed morphologically by Nottingham pathological grade [18]. We found that the intratumoral adipocyte-low tumors were associated with advanced grade, grade 2 or 3, compared with grade 1 in the whole cohort and in the ER+/HER2- subtype (Figure 6C; p < 0.001 and p < 0.001 respectively). This result was also true in GSE25066 cohort (Figure 6C).

3. Discussion

In the current study, transcriptionally defined adipocytes reasonably correlated with previously reported markers of mature adipocytes. Breast cancer contained the second highest amount of intratumoral adipocytes among the cancers listed in TCGA. There was no survival difference nor pCR rate after NAC by the amount of intratumoral adipocytes in any of the subtypes in the TCGA cohort, nor in any of the NAC cohort. Intratumoral adipocyte-high tumors significantly enriched for inflammation and metastasis-related gene sets as well as favorable tumor immune microenvironment, particularly in the ER+/HER2- subtype from TCGA and in the GSE25066 cohort. On the other hand, intratumoral adipocyte-low tumors significantly enriched for cell cycle and cell proliferation-related gene sets in all subtypes in both the TCGA and GSE25066 cohorts. Intratumoral adipocyte-low tumors were significantly associated with a higher MKI67 expression and demonstrated an inverse correlation. Intratumoral adipocyte-low tumors were associated with advanced pathological grades, which is a clinical parameter of aggressive, highly proliferative breast cancer.
We have found that intratumoral adipocyte-high breast cancer enriched for the inflammation, metastasis, and cancer stemness related gene sets, which is in agreement with the previous reports [25,49,50,51,52]. This was more significant with the ER+/HER2- subtype compared with TNBC. Manabe et al. reported that the co-culture of mature adipocyte and ER+ breast cancer cells, such as MCF7, promoted the growth of these cancer cells [52]. In addition, a study by Lee et al. demonstrated that co-culture of adipocytes with cancer cells induced the upregulation of EMT-related genes, which led to the promotion of cancer cell invasion and migration [51]. Our study of clinical samples from large databases reflected an association with mechanisms of EMT that has been suggested by these in vitro experiments. Furthermore, intratumoral adipocyte-high tumors enriched for the gene sets associated with metastasis, such as TGFβ Signaling, Notch Signaling, Angiogenesis, WNTβ Catenin Signaling, and Hedgehog Signaling. This was more evident with the ER+/HER2- subtype compared with TNBC, which indicates that the role of adipocytes is more significant with the ER+/HER2- subtype. Our result is in perfect alignment with the previous studies that leptin, one of the well-known adipokines which enhances the cell proliferation and metastasis [53,54], and leptin receptor signaling crosstalk with ER signal pathway in ER+ subtypes but not in ER- subtypes [55,56].
Intratumoral adipocyte-high tumors are also enriched for inflammation-related gene sets, such as TNF-α signaling via NFκB, IL-6/JAK/STAT3 signaling, and inflammatory response in the TCGA cohort. Studies have shown the detrimental role of inflammatory cytokines, TNF-α and IL-6, in the promotion of cancer progression [6,7,12,17,40]. It has been reported that cancer associated adipocytes secrete IL-6 and TNF-α [25,26]. IL-6 activates the JAK/STAT3 pathway and promotes angiogenesis, invasion, metastasis [29,57], and TNF-α activate EMT and promote cancer progression [32]. Our result is consistent with these previous studies that a high amount of intratumoral adipokines associates with inflammation that aggravates cancer in the clinical setting.
Intratumoral adipocyte-high tumors also enriched for gene sets related to immune response such as IL2 signaling and IFN-γ and are associated with a smaller number of unfavorable immune cells such as Th2, Treg, and M2 macrophages, as well as with higher CYT. Overall, these associations suggest a favorable tumor immune microenvironment, especially in the ER+/HER2- subtype. Intratumoral adipocytes may contribute to breast cancer through inflammation, metastasis, and cancer stemness, but also associated with a favorable immune response that slows the cancer aggressiveness down, particularly in ER+/HER2- subtype.
It was surprising that the survival outcome did not differ between intratumoral adipocyte-high and low tumors in all the subtypes analyzed in the current study, especially when we found that intratumoral adipocyte associated with inflammation, metastasis, and cancer stemness, which have been previously reported to aggravate cancer [25,26,29]. Indeed, while the favorable immune response was associated with adipocyte-high tumors, there was no strong association with anti-cancer immune cell infiltration. Therefore, it is unlikely that favorable immune cell accumulation cancels out the detrimental effect of adipocytes. Conversely, we found was that intratumoral adipocyte-low tumors significantly enriched for four out of five designated hallmark cell proliferation gene sets (Figure 5), suggesting that they are extremely highly proliferative cancer. It is possible that the dense cellularity of highly proliferative cancer may not allow much space for adipocytes to infiltrate by physical pressure. Indeed, adipocytes are seldomly observed in the histological section of breast cancers, especially in the central portion of cancer where desmoplasia takes place [8]. The TCGA sample protocol mandated that samples be collected from the location within the bulk tumor that contains an average of 60% tumor cell nuclei with less than 20% necrosis [58]. Thus, it needs to be noted that the adipocytes analyzed in the current study are intratumoral adipocytes located closer to the center of the tumor. These are distinctively different from adipocytes from the peri-tumoral location, where adipocytes surrounding but not infiltrating the tumor may be included. Intratumoral adipocyte-low tumors enriched for G2M checkpoint gene set, which we recently reported to associate with aggressive characteristics in breast cancer [59]. The tumors with enhanced G2M checkpoint pathway demonstrated more advanced AJCC stage as well as more advanced Nottingham pathological grade, and thus increased risk of distant metastasis. We also reported that enhanced E2F signaling associated with clinically aggressive breast cancer [60]. Taken together, we speculate that the aggressive characteristics of adipocyte high tumors that cause inflammation and metastasis may be counterbalanced by the highly proliferative biology of adipocyte low tumors, leading to the no survival difference of intratumoral adipocyte-high versus low tumors.
This report is a retrospective study that mainly used publicly open databases, such as TCGA. There are several caveats to this study. TCGA and other publicly available breast cancer cohorts used in this study lack information on social history, which prevented us from analyzing the association between intra-tumoral adipocyte and social history, including smoking. One of the main strengths, as well as a limitation of this study, is that we transcriptomically defined the adipocyte. This approach allows us to assess the amount of intratumoral adipocytes in any bulk tumor with gene expression data, even when there was no intention to measure adipocytes at the time of the initial sample collection. However, because this approach utilizes the publicly available cohorts, our analyses are limited to the clinical parameters, the quality, and exact spatial location of where the sample was taken by the original authors. Further, although we have demonstrated that our methodology appropriately correlated with markers of mature adipocyte, there is a possibility that we may be picking up signals other than adipocytes. In addition, the mainstays of analyzing immune cells are immunohistochemistry and flow cytometry. However, these analyses were not performed due to the lack of access to the patient samples. We also did not conduct any in vitro and in vivo experiments, and thus are reliant on current literature to understand the underlying mechanisms.
In conclusion, a high amount of intratumoral adipocyte was associated with inflammation, metastatic pathways, cancer stemness, and favorable tumor immune microenvironment; whereas, a low amount of intratumoral adipocytes was associated with a highly proliferative tumor in ER-positive breast cancer. This cancer biology may explain the reason why patient survival did not differ by the amount of adipocytes.

4. Materials and Methods

4.1. Data Acquisition

A total of 1090 patient clinicopathological and gene expression data of The Cancer Genome Atlas (TCGA) Pan-Cancer study (TCGA PanCancer Atlas) was obtained through cBioPortal as previously described [17,19,40,41,48,60,61].
Clinical and normalized microarray-based gene expression data were obtained from the following cohorts that underwent neoadjuvant chemotherapies; Symmans et al. (GSE25066; n = 508, taxane and anthracycline) [44], Noguchi et al. (GSE32646; n = 115; 5-FU, epirubicin, cyclophosphamide, and paclitaxel) [45], and Shi et al. (GSE20194; n = 278; paclitaxel, 5-FU, cyclophosphamide and doxorubicin) [46].
The data of a study by Yau et al. (n = 683) [62] was obtained through UCSC Xena (http://xena.ucsc.edu, University of California Santa Cruz, Santa Cruz, CA, USA) [63]. Given that all the cohorts used in this study were in a publicly available de-identified database, approval from the Institutional Review Board was waived for this study.

4.2. Gene Set Enrichment Analysis (GSEA)

Gene set enrichment analysis (GSEA) was a publicly available software (GSEA version 4.0) provided by the Broad Institute (Cambridge, MA, USA) (http://software.broadinstitute.org/gsea/index.jsp). For GSEA, we utilized Hallmark gene sets used for this study, as previously described [14,18,19,40,41,47,64,65]. False discovery rate (FDR) of 0.25 was utilized to define the statistical significance of GSEA throughout the study, as the Broad Institute recommends.

4.3. Immune Cell Composition and Scores Related with Immune Activity

Intratumoral adipocytes were defined with a computational algorithm, xCell, which utilized the transcriptome data as described in the journal, Genome Biology, in 2017 by Aran D et al. [42]. In addition, xCell was used to analyze the cell composition of all the other immune cells examined in this study. Cytolytic activity (CYT) was calculated using the geometric mean of granzyme A and Perforin 1 expression values, as described previously [39,59,66].

4.4. Statistical Analysis

All the statistical analysis was performed using R software (version 4.0.2; http:///www.r-project.org/). Kaplan–Meier survival analysis was performed with greyzoneSurv packages in R for the survival analysis. One-way ANOVA or Fisher’s exact test was used to determine the significance of differences for groups. Other box plots were analyzed with one-way ANOVA. A two-sided p-value  <  0.05 was considered statistically significant. All boxplots are of the Tukey type, and the boxes depicted medians and inter-quartile ranges.

Supplementary Materials

The following are available online at https://www.mdpi.com/1422-0067/21/16/5744/s1, Figure S1: Relationships between pathological complete response (pCR) rate to neoadjuvant chemotherapy (NAC) and amount of intratumoral adipocyte, Figure S2: Stage was not associated with the amount of adipocytes, Figure S3: Intratumoral-high adipocyte tumors enriched the gene sets related with inflammation (A), metastasis (B) and immune response (C) in ER+/HER2- subtype, Figure S4: Intratumoral-high adipocyte tumors are enriched for gene sets related with inflammation (A), metastasis (B) and immune response (C) in the TNBC subtype, Figure S5: GSEA demonstrated that the intratumoral adipocyte-low tumors enriched the gene sets associated with cell cycle and cell proliferation in ER+/HER2- subtype, Figure S6: GSEA demonstrated that the intratumoral adipocyte-low tumors enriched the gene sets associated with cell cycle and cell proliferation in TNBC subtype, Table S1: The genes used to define adipocytes in xCell, the computational algorithm reported by Aran et al.

Author Contributions

Conceptualization, Y.T.; methodology, M.O., E.K.; formal analysis, Y.T.; resources, L.Y., N.M., M.F., K.Y.; writing—original draft preparation, Y.T.; writing—review and editing, J.L.H., M.N., K.T.; supervision, K.T.; project administration, K.T. All authors have read and agree to the published version of the manuscript.

Funding

This work was supported by NIH grant R01CA160688 to K.T., and National Cancer Institute (NCI) grant P30CA016056 involving the use of Roswell Park Cancer Institute’s Resources.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

EREstrogen receptor
HER2Human epidermal growth factor receptor 2
TNBCTriple negative breast cancer
TCGAThe Cancer Genome Atlas
METABRICMolecular Taxonomy of Breast Cancer International Consortium
GSEAGene Set Enrichment Analysis
CYTCytolytic Activity

References

  1. Siegel, R.L.; Miller, K.D.; Jemal, A. Cancer statistics, 2020. CA Cancer J. Clin. 2020, 70, 7–30. [Google Scholar] [CrossRef] [PubMed]
  2. Bray, F.; Ferlay, J.; Soerjomataram, I.; Siegel, R.L.; Torre, L.A.; Jemal, A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 2018, 68, 394–424. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Picon-Ruiz, M.; Morata-Tarifa, C.; Valle-Goffin, J.J.; Friedman, E.R.; Slingerland, J.M. Obesity and adverse breast cancer risk and outcome: Mechanistic insights and strategies for intervention. CA Cancer J. Clin. 2017, 67, 378–397. [Google Scholar] [CrossRef] [PubMed]
  4. Li, J.; Han, X. Adipocytokines and breast cancer. Curr. Probl. Cancer 2018, 42, 208–214. [Google Scholar] [CrossRef]
  5. Quail, D.F.; Dannenberg, A.J. The obese adipose tissue microenvironment in cancer development and progression. Nat. Rev. Endocrinol. 2019, 15, 139–154. [Google Scholar] [CrossRef]
  6. Nagahashi, M.; Yamada, A.; Katsuta, E.; Aoyagi, T.; Huang, W.C.; Terracina, K.P.; Hait, N.C.; Allegood, J.C.; Tsuchida, J.; Yuza, K.; et al. Targeting the SphK1/S1P/S1PR1 Axis That Links Obesity, Chronic Inflammation, and Breast Cancer Metastasis. Cancer Res. 2018, 78, 1713–1725. [Google Scholar] [CrossRef] [Green Version]
  7. Katsuta, E.; Yan, L.; Nagahashi, M.; Raza, A.; Sturgill, J.L.; Lyon, D.E.; Rashid, O.M.; Hait, N.C.; Takabe, K. Doxorubicin effect is enhanced by sphingosine-1-phosphate signaling antagonist in breast cancer. J. Surg. Res. 2017, 219, 202–213. [Google Scholar] [CrossRef]
  8. Choi, J.; Cha, Y.J.; Koo, J.S. Adipocyte biology in breast cancer: From silent bystander to active facilitator. Prog. Lipid Res. 2018, 69, 11–20. [Google Scholar] [CrossRef]
  9. Hinshaw, D.C.; Shevde, L.A. The Tumor Microenvironment Innately Modulates Cancer Progression. Cancer Res. 2019, 79, 4557–4566. [Google Scholar] [CrossRef] [Green Version]
  10. Katsuta, E.; Rashid, O.M.; Takabe, K. Clinical relevance of tumor microenvironment: Immune cells, vessels, and mouse models. Hum. Cell 2020. [Google Scholar] [CrossRef]
  11. Nojiri, T.; Hosoda, H.; Tokudome, T.; Miura, K.; Ishikane, S.; Otani, K.; Kishimoto, I.; Shintani, Y.; Inoue, M.; Kimura, T.; et al. Atrial natriuretic peptide prevents cancer metastasis through vascular endothelial cells. Proc. Natl. Acad. Sci. USA 2015, 112, 4086–4091. [Google Scholar] [CrossRef] [Green Version]
  12. Nagahashi, M.; Abe, M.; Sakimura, K.; Takabe, K.; Wakai, T. The role of sphingosine-1-phosphate in inflammation and cancer progression. Cancer Sci. 2018, 109, 3671–3678. [Google Scholar] [CrossRef] [Green Version]
  13. Sidaway, P. Breast cancer: Fibroblast subtypes alter the microenvironment. Nat. Rev. Clin. Oncol. 2018, 15, 264–265. [Google Scholar] [PubMed]
  14. Katsuta, E.; Rashid, O.M.; Takabe, K. Fibroblasts as a Biological Marker for Curative Resection in Pancreatic Ductal Adenocarcinoma. Int. J. Mol. Sci. 2020, 21, 3890. [Google Scholar] [CrossRef] [PubMed]
  15. Cha, Y.J.; Kim, E.S.; Koo, J.S. Tumor-associated macrophages and crown-like structures in adipose tissue in breast cancer. Breast Cancer Res. Treat. 2018, 170, 15–25. [Google Scholar] [CrossRef] [PubMed]
  16. Gajewski, T.F.; Schreiber, H.; Fu, Y.X. Innate and adaptive immune cells in the tumor microenvironment. Nat. Immunol. 2013, 14, 1014–1022. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Tokumaru, Y.; Oshi, M.; Katsuta, E.; Yan, L.; Satyananda, V.; Matsuhashi, N.; Futamura, M.; Akao, Y.; Yoshida, K.; Takabe, K. KRAS signaling enriched triple negative breast cancer is associated with favorable tumor immune microenvironment and better survival. Am. J. Cancer Res. 2020, 10, 897–907. [Google Scholar]
  18. Tokumaru, Y.; Asaoka, M.; Oshi, M.; Katsuta, E.; Yan, L.; Narayanan, S.; Sugito, N.; Matsuhashi, N.; Futamura, M.; Akao, Y.; et al. High Expression of microRNA-143 is Associated with Favorable Tumor Immune Microenvironment and Better Survival in Estrogen Receptor Positive Breast Cancer. Int. J. Mol. Sci. 2020, 21, 3213. [Google Scholar] [CrossRef]
  19. Takeshita, T.; Yan, L.; Asaoka, M.; Rashid, O.; Takabe, K. Late recurrence of breast cancer is associated with pro-cancerous immune microenvironment in the primary tumor. Sci. Rep. 2019, 9, 16942. [Google Scholar] [CrossRef] [Green Version]
  20. Cassim, S.; Pouyssegur, J. Tumor Microenvironment: A Metabolic Player that Shapes the Immune Response. Int. J. Mol. Sci. 2019, 21, 157. [Google Scholar] [CrossRef] [Green Version]
  21. DeNardo, D.G.; Coussens, L.M. Inflammation and breast cancer. Balancing immune response: Crosstalk between adaptive and innate immune cells during breast cancer progression. Breast Cancer Res. BCR 2007, 9, 212. [Google Scholar] [CrossRef] [PubMed]
  22. Wu, Q.; Li, B.; Li, Z.; Li, J.; Sun, S.; Sun, S. Cancer-associated adipocytes: Key players in breast cancer progression. J. Hematol. Oncol. 2019, 12, 95. [Google Scholar] [CrossRef] [PubMed]
  23. Gyamfi, J.; Eom, M.; Koo, J.S.; Choi, J. Multifaceted Roles of Interleukin-6 in Adipocyte-Breast Cancer Cell Interaction. Transl. Oncol. 2018, 11, 275–285. [Google Scholar] [CrossRef] [PubMed]
  24. Balkwill, F. Tumour necrosis factor and cancer. Nat. Rev. Cancer 2009, 9, 361–371. [Google Scholar] [CrossRef] [PubMed]
  25. Rybinska, I.; Agresti, R.; Trapani, A.; Tagliabue, E.; Triulzi, T. Adipocytes in Breast Cancer, the Thick and the Thin. Cells 2020, 9, 560. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Chu, D.T.; Phuong, T.N.T.; Tien, N.L.B.; Tran, D.K.; Nguyen, T.T.; Thanh, V.V.; Quang, T.L.; Minh, L.B.; Pham, V.H.; Ngoc, V.T.N.; et al. The Effects of Adipocytes on the Regulation of Breast Cancer in the Tumor Microenvironment: An Update. Cells 2019, 8, 857. [Google Scholar] [CrossRef] [Green Version]
  27. Zhu, Z.; Zhu, X.; Yang, S.; Guo, Z.; Li, K.; Ren, C.; Zhou, Y.; Dou, J. Yin-yang effect of tumour cells in breast cancer: From mechanism of crosstalk between tumour-associated macrophages and cancer-associated adipocytes. Am. J. Cancer Res. 2020, 10, 383–392. [Google Scholar]
  28. Yu, H.; Lee, H.; Herrmann, A.; Buettner, R.; Jove, R. Revisiting STAT3 signalling in cancer: New and unexpected biological functions. Nat. Rev. Cancer 2014, 14, 736–746. [Google Scholar] [CrossRef]
  29. Deng, T.; Lyon, C.J.; Bergin, S.; Caligiuri, M.A.; Hsueh, W.A. Obesity, Inflammation, and Cancer. Annu. Rev. Pathol. 2016, 11, 421–449. [Google Scholar] [CrossRef] [Green Version]
  30. Leibovich-Rivkin, T.; Liubomirski, Y.; Bernstein, B.; Meshel, T.; Ben-Baruch, A. Inflammatory factors of the tumor microenvironment induce plasticity in nontransformed breast epithelial cells: EMT, invasion, and collapse of normally organized breast textures. Neoplasia (New York N.Y.) 2013, 15, 1330–1346. [Google Scholar] [CrossRef] [Green Version]
  31. Li, C.W.; Xia, W.; Huo, L.; Lim, S.O.; Wu, Y.; Hsu, J.L.; Chao, C.H.; Yamaguchi, H.; Yang, N.K.; Ding, Q.; et al. Epithelial-mesenchymal transition induced by TNF-α requires NF-κB-mediated transcriptional upregulation of Twist1. Cancer Res. 2012, 72, 1290–1300. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Chua, H.L.; Bhat-Nakshatri, P.; Clare, S.E.; Morimiya, A.; Badve, S.; Nakshatri, H. NF-kappaB represses E-cadherin expression and enhances epithelial to mesenchymal transition of mammary epithelial cells: Potential involvement of ZEB-1 and ZEB-2. Oncogene 2007, 26, 711–724. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Tsuyada, A.; Chow, A.; Wu, J.; Somlo, G.; Chu, P.; Loera, S.; Luu, T.; Li, A.X.; Wu, X.; Ye, W.; et al. CCL2 mediates cross-talk between cancer cells and stromal fibroblasts that regulates breast cancer stem cells. Cancer Res. 2012, 72, 2768–2779. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Crean-Tate, K.K.; Reizes, O. Leptin Regulation of Cancer Stem Cells in Breast and Gynecologic Cancer. Endocrinology 2018, 159, 3069–3080. [Google Scholar] [CrossRef]
  35. Newman, G.; Gonzalez-Perez, R.R. Leptin-cytokine crosstalk in breast cancer. Mol. Cell. Endocrinol. 2014, 382, 570–582. [Google Scholar] [CrossRef] [Green Version]
  36. Gonzalez-Perez, R.R.; Lanier, V.; Newman, G. Leptin’s Pro-Angiogenic Signature in Breast Cancer. Cancers 2013, 5, 1140–1162. [Google Scholar] [CrossRef] [Green Version]
  37. Takahashi, H.; Asaoka, M.; Yan, L.; Rashid, O.M.; Oshi, M.; Ishikawa, T.; Nagahashi, M.; Takabe, K. Biologically Aggressive Phenotype and Anti-cancer Immunity Counterbalance in Breast Cancer with High Mutation Rate. Sci. Rep. 2020, 10, 1852. [Google Scholar] [CrossRef] [Green Version]
  38. Asaoka, M.; Patnaik, S.K.; Zhang, F.; Ishikawa, T.; Takabe, K. Lymphovascular invasion in breast cancer is associated with gene expression signatures of cell proliferation but not lymphangiogenesis or immune response. Breast Cancer Res. Treat. 2020, 181, 309–322. [Google Scholar] [CrossRef]
  39. Takahashi, H.; Kawaguchi, T.; Yan, L.; Peng, X.; Qi, Q.; Morris, L.G.T.; Chan, T.A.; Tsung, A.; Otsuji, E.; Takabe, K. Immune Cytolytic Activity for Comprehensive Understanding of Immune Landscape in Hepatocellular Carcinoma. Cancers 2020, 12, 1221. [Google Scholar] [CrossRef]
  40. Okano, M.; Oshi, M.; Butash, A.L.; Katsuta, E.; Tachibana, K.; Saito, K.; Okayama, H.; Peng, X.; Yan, L.; Kono, K.; et al. Triple-Negative Breast Cancer with High Levels of Annexin A1 Expression Is Associated with Mast Cell Infiltration, Inflammation, and Angiogenesis. Int. J. Mol. Sci. 2019, 20, 4197. [Google Scholar] [CrossRef] [Green Version]
  41. Okano, M.; Oshi, M.; Butash, A.L.; Asaoka, M.; Katsuta, E.; Peng, X.; Qi, Q.; Yan, L.; Takabe, K. Estrogen Receptor Positive Breast Cancer with High Expression of Androgen Receptor has Less Cytolytic Activity and Worse Response to Neoadjuvant Chemotherapy but Better Survival. Int. J. Mol. Sci. 2019, 20, 2655. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Aran, D.; Hu, Z.; Butte, A.J. xCell: Digitally portraying the tissue cellular heterogeneity landscape. Genome Biol. 2017, 18, 220. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Lehuédé, C.; Li, X.; Dauvillier, S.; Vaysse, C.; Franchet, C.; Clement, E.; Esteve, D.; Longué, M.; Chaltiel, L.; Le Gonidec, S.; et al. Adipocytes promote breast cancer resistance to chemotherapy, a process amplified by obesity: Role of the major vault protein (MVP). Breast Cancer Res. BCR 2019, 21, 7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Hatzis, C.; Pusztai, L.; Valero, V.; Booser, D.J.; Esserman, L.; Lluch, A.; Vidaurre, T.; Holmes, F.; Souchon, E.; Wang, H.; et al. A genomic predictor of response and survival following taxane-anthracycline chemotherapy for invasive breast cancer. JAMA 2011, 305, 1873–1881. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Miyake, T.; Nakayama, T.; Naoi, Y.; Yamamoto, N.; Otani, Y.; Kim, S.J.; Shimazu, K.; Shimomura, A.; Maruyama, N.; Tamaki, Y.; et al. GSTP1 expression predicts poor pathological complete response to neoadjuvant chemotherapy in ER-negative breast cancer. Cancer Sci. 2012, 103, 913–920. [Google Scholar] [CrossRef] [PubMed]
  46. Shi, L.; Campbell, G.; Jones, W.D.; Campagne, F.; Wen, Z.; Walker, S.J.; Su, Z.; Chu, T.M.; Goodsaid, F.M.; Pusztai, L.; et al. The MicroArray Quality Control (MAQC)-II study of common practices for the development and validation of microarray-based predictive models. Nat. Biotechnol. 2010, 28, 827–838. [Google Scholar]
  47. Tokumaru, Y.; Katsuta, E.; Oshi, M.; Sporn, J.C.; Yan, L.; Le, L.; Matsuhashi, N.; Futamura, M.; Akao, Y.; Yoshida, K.; et al. High Expression of miR-34a Associated with Less Aggressive Cancer Biology but Not with Survival in Breast Cancer. Int. J. Mol. Sci. 2020, 21, 3045. [Google Scholar] [CrossRef]
  48. Takahashi, H.; Oshi, M.; Asaoka, M.; Yan, L.; Endo, I.; Takabe, K. Molecular Biological Features of Nottingham Histological Grade 3 Breast Cancers. Ann. Surg. Oncol. 2020. [Google Scholar] [CrossRef]
  49. Osman, A.; Afify, S.M.; Hassan, G.; Fu, X.; Seno, A.; Seno, M. Revisiting Cancer Stem Cells as the Origin of Cancer-Associated Cells in the Tumor Microenvironment: A Hypothetical View from the Potential of iPSCs. Cancers 2020, 12, 879. [Google Scholar] [CrossRef] [Green Version]
  50. Zaoui, M.; Morel, M.; Ferrand, N.; Fellahi, S.; Bastard, J.P.; Lamazière, A.; Larsen, A.K.; Béréziat, V.; Atlan, M.; Sabbah, M. Breast-Associated Adipocytes Secretome Induce Fatty Acid Uptake and Invasiveness in Breast Cancer Cells via CD36 Independently of Body Mass Index, Menopausal Status and Mammary Density. Cancers 2019, 11, 2012. [Google Scholar] [CrossRef] [Green Version]
  51. Lee, Y.; Jung, W.H.; Koo, J.S. Adipocytes can induce epithelial-mesenchymal transition in breast cancer cells. Breast Cancer Res. Treat. 2015, 153, 323–335. [Google Scholar] [CrossRef] [PubMed]
  52. Manabe, Y.; Toda, S.; Miyazaki, K.; Sugihara, H. Mature adipocytes, but not preadipocytes, promote the growth of breast carcinoma cells in collagen gel matrix culture through cancer-stromal cell interactions. J. Pathol. 2003, 201, 221–228. [Google Scholar] [CrossRef] [PubMed]
  53. Shin, E.; Koo, J.S. The Role of Adipokines and Bone Marrow Adipocytes in Breast Cancer Bone Metastasis. Int. J. Mol. Sci. 2020, 21, 4967. [Google Scholar] [CrossRef] [PubMed]
  54. Maroni, P. Leptin, Adiponectin, and Sam68 in Bone Metastasis from Breast Cancer. Int. J. Mol. Sci. 2020, 21, 1051. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Khanal, T.; Kim, H.G.; Do, M.T.; Choi, J.H.; Won, S.S.; Kang, W.; Chung, Y.C.; Jeong, T.C.; Jeong, H.G. Leptin induces CYP1B1 expression in MCF-7 cells through ligand-independent activation of the ERα pathway. Toxicol. Appl. Pharmacol. 2014, 277, 39–48. [Google Scholar] [CrossRef]
  56. Strong, A.L.; Ohlstein, J.F.; Biagas, B.A.; Rhodes, L.V.; Pei, D.T.; Tucker, H.A.; Llamas, C.; Bowles, A.C.; Dutreil, M.F.; Zhang, S.; et al. Leptin produced by obese adipose stromal/stem cells enhances proliferation and metastasis of estrogen receptor positive breast cancers. Breast Cancer Res. BCR 2015, 17, 112. [Google Scholar] [CrossRef] [Green Version]
  57. Johnson, D.E.; O’Keefe, R.A.; Grandis, J.R. Targeting the IL-6/JAK/STAT3 signalling axis in cancer. Nat. Rev. Clin. Oncol. 2018, 15, 234–248. [Google Scholar] [CrossRef]
  58. The Cancer Genome Atlas Network. Comprehensive molecular portraits of human breast tumours. Nature 2012, 490, 61–70. [Google Scholar] [CrossRef] [Green Version]
  59. Oshi, M.; Takahashi, H.; Tokumaru, Y.; Yan, L.; Rashid, O.M.; Matsuyama, R.; Endo, I.; Takabe, K. G2M Cell Cycle Pathway Score as a Prognostic Biomarker of Metastasis in Estrogen Receptor (ER)-Positive Breast Cancer. Int. J. Mol. Sci. 2020, 21, 2921. [Google Scholar] [CrossRef] [Green Version]
  60. Oshi, M.; Takahashi, H.; Tokumaru, Y.; Yan, L.; Rashid, O.M.; Nagahashi, M.; Matsuyama, R.; Endo, I.; Takabe, K. The E2F Pathway Score as a Predictive Biomarker of Response to Neoadjuvant Therapy in ER+/HER2- Breast Cancer. Cells 2020, 9, 1643. [Google Scholar] [CrossRef]
  61. Katsuta, E.; Yan, L.; Takeshita, T.; McDonald, K.A.; Dasgupta, S.; Opyrchal, M.; Takabe, K. High MYC mRNA Expression Is More Clinically Relevant than MYC DNA Amplification in Triple-Negative Breast Cancer. Int. J. Mol. Sci. 2019, 21, 217. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Yau, C.; Esserman, L.; Moore, D.H.; Waldman, F.; Sninsky, J.; Benz, C.C. A multigene predictor of metastatic outcome in early stage hormone receptor-negative and triple-negative breast cancer. Breast Cancer Res. BCR 2010, 12, R85. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Goldman, M.J.; Craft, B.; Hastie, M.; Repečka, K.; McDade, F.; Kamath, A.; Banerjee, A.; Luo, Y.; Rogers, D.; Brooks, A.N.; et al. Visualizing and interpreting cancer genomics data via the Xena platform. Nat. Biotechnol. 2020, 38, 675–678. [Google Scholar] [CrossRef] [PubMed]
  64. Oshi, M.; Katsuta, E.; Yan, L.; Ebos, J.M.L.; Rashid, O.M.; Matsuyama, R.; Endo, I.; Takabe, K. A Novel 4-Gene Score to Predict Survival, Distant Metastasis and Response to Neoadjuvant Therapy in Breast Cancer. Cancers 2020, 12, 1148. [Google Scholar] [CrossRef]
  65. Gandhi, S.; Elkhanany, A.; Oshi, M.; Dai, T.; Opyrchal, M.; Mohammadpour, H.; Repasky, E.A.; Takabe, K. Contribution of Immune Cells to Glucocorticoid Receptor Expression in Breast Cancer. Int. J. Mol. Sci. 2020, 21, 4635. [Google Scholar] [CrossRef]
  66. Rooney, M.S.; Shukla, S.A.; Wu, C.J.; Getz, G.; Hacohen, N. Molecular and genetic properties of tumors associated with local immune cytolytic activity. Cell 2015, 160, 48–61. [Google Scholar] [CrossRef] [Green Version]
Ijms 21 05744 g001 550
Figure 1. Transcriptionally defined adipocytes have a moderate to strong correlation with representative markers of mature adipocytes. The amount of intratumoral adipocytes in various cancers is also represented. (A) The Cancer Genome Atlas (TCGA) cohort (B) GSE25066 (C) Study by Yau et al. (D) The distribution of the amount of intratumoral adipocytes in various cancer types. ADIPOQ, adiponectin; LEP, leptin; LPL, lipoprotein lipase; PLIN1, perilipin 1; ACC, Adrenocortical carcinoma; BLCA, Bladder Urothelial Carcinoma; BRCA, Breast invasive carcinoma; CESC, Cervical squamous cell carcinoma and endocervical adenocarcinoma; CHOL, Cholangiocarcinoma; COAD, Colon adenocarcinoma; DLBC, Lymphoid Neoplasm Diffuse Large B-cell Lymphoma; ESCA, Esophageal carcinoma; GBM, Glioblastoma multiforme; HNSC, Head and Neck squamous cell carcinoma; KICH, Kidney Chromophobe; KIRC, Kidney renal clear cell carcinoma; KIRP, Kidney renal papillary cell carcinoma; LAML, Acute Myeloid Leukemia; LGG, Brain Lower Grade Glioma; LIHC, Liver hepatocellular carcinoma; LUAD, Lung adenocarcinoma; LUSC, Lung squamous cell carcinoma; MESO, Mesothelioma; OV, Ovarian serous cystadenocarcinoma; PAAD, Pancreatic adenocarcinoma; PCPG, Pheochromocytoma and Paraganglioma; PRAD, Prostate adenocarcinoma; READ, Rectum adenocarcinoma; SARC, Sarcoma; SKCM, Skin Cutaneous Melanoma; STAD, Stomach adenocarcinoma; TGCT, Testicular Germ Cell Tumors; THCA, Thyroid carcinoma; THYM, Thymoma; UCEC, Uterine Corpus Endometrial Carcinoma; UCS, Uterine Carcinosarcoma; UVM, Uveal Melanoma. p-value  <  0.05 was considered statistically significant. r represents a Spearman’s correlation coefficient.
Figure 1. Transcriptionally defined adipocytes have a moderate to strong correlation with representative markers of mature adipocytes. The amount of intratumoral adipocytes in various cancers is also represented. (A) The Cancer Genome Atlas (TCGA) cohort (B) GSE25066 (C) Study by Yau et al. (D) The distribution of the amount of intratumoral adipocytes in various cancer types. ADIPOQ, adiponectin; LEP, leptin; LPL, lipoprotein lipase; PLIN1, perilipin 1; ACC, Adrenocortical carcinoma; BLCA, Bladder Urothelial Carcinoma; BRCA, Breast invasive carcinoma; CESC, Cervical squamous cell carcinoma and endocervical adenocarcinoma; CHOL, Cholangiocarcinoma; COAD, Colon adenocarcinoma; DLBC, Lymphoid Neoplasm Diffuse Large B-cell Lymphoma; ESCA, Esophageal carcinoma; GBM, Glioblastoma multiforme; HNSC, Head and Neck squamous cell carcinoma; KICH, Kidney Chromophobe; KIRC, Kidney renal clear cell carcinoma; KIRP, Kidney renal papillary cell carcinoma; LAML, Acute Myeloid Leukemia; LGG, Brain Lower Grade Glioma; LIHC, Liver hepatocellular carcinoma; LUAD, Lung adenocarcinoma; LUSC, Lung squamous cell carcinoma; MESO, Mesothelioma; OV, Ovarian serous cystadenocarcinoma; PAAD, Pancreatic adenocarcinoma; PCPG, Pheochromocytoma and Paraganglioma; PRAD, Prostate adenocarcinoma; READ, Rectum adenocarcinoma; SARC, Sarcoma; SKCM, Skin Cutaneous Melanoma; STAD, Stomach adenocarcinoma; TGCT, Testicular Germ Cell Tumors; THCA, Thyroid carcinoma; THYM, Thymoma; UCEC, Uterine Corpus Endometrial Carcinoma; UCS, Uterine Carcinosarcoma; UVM, Uveal Melanoma. p-value  <  0.05 was considered statistically significant. r represents a Spearman’s correlation coefficient.
Ijms 21 05744 g001
Ijms 21 05744 g002 550
Figure 2. Kaplan Meier survival analysis (Disease-free survival (DFS) and overall survival (OS)) of adipocyte high (red line) vs. low (blue line) breast cancer in whole TCGA breast cancer cohort, ER+/HER2- subtype, Her2+ subtype, and TNBC subtype. The number of patients at risk is shown below the X-axis of each panel. DFS, disease-free survival; OS, overall survival; ER+/HER2-, ER-positive, and HER2 negative; HER2+, HER2 positive; TNBC, triple-negative breast cancer. p-value  <  0.05 was considered statistically significant.
Figure 2. Kaplan Meier survival analysis (Disease-free survival (DFS) and overall survival (OS)) of adipocyte high (red line) vs. low (blue line) breast cancer in whole TCGA breast cancer cohort, ER+/HER2- subtype, Her2+ subtype, and TNBC subtype. The number of patients at risk is shown below the X-axis of each panel. DFS, disease-free survival; OS, overall survival; ER+/HER2-, ER-positive, and HER2 negative; HER2+, HER2 positive; TNBC, triple-negative breast cancer. p-value  <  0.05 was considered statistically significant.
Ijms 21 05744 g002
Ijms 21 05744 g003a 550Ijms 21 05744 g003b 550
Figure 3. The gene set enrichment analysis (GSEA) of adipocyte high tumors in the whole cohort of TCGA breast cancer and GSE25066. (A) GSEA of Hallmark inflammation-related gene sets. (B) GSEA of Hallmark metastasis- and cancer stemness-related gene sets. (C) GSEA of Hallmark immune response-related gene sets. NES, normalized enrichment score; FDR, false discovery rate. FDR < 0.25 was considered statistically significant.
Figure 3. The gene set enrichment analysis (GSEA) of adipocyte high tumors in the whole cohort of TCGA breast cancer and GSE25066. (A) GSEA of Hallmark inflammation-related gene sets. (B) GSEA of Hallmark metastasis- and cancer stemness-related gene sets. (C) GSEA of Hallmark immune response-related gene sets. NES, normalized enrichment score; FDR, false discovery rate. FDR < 0.25 was considered statistically significant.
Ijms 21 05744 g003aIjms 21 05744 g003b
Ijms 21 05744 g004 550
Figure 4. Immune cell composition by high and low adipocytes in the whole TCGA cohort, ER+/Her2-, and TNBC subtypes. (A) Favorable immune cells: CD8+ (CD8+ T cell), CD4+ (CD4+ T cell), Th1 (type 1 helper T cell), M1, (M1 macrophage), DC (Dendritic Cell). (B) Unfavorable immune cells: Th2 (type 2 helper T cell), Treg (regulatory T cell), M2 (M2 macrophage. (C) Cytolytic activity score (CYT). ER+/HER2-, ER-positive and HER2 negative; HER2+, HER2 positive; TNBC, triple-negative breast cancer. p-value  <  0.05 was considered statistically significant.
Figure 4. Immune cell composition by high and low adipocytes in the whole TCGA cohort, ER+/Her2-, and TNBC subtypes. (A) Favorable immune cells: CD8+ (CD8+ T cell), CD4+ (CD4+ T cell), Th1 (type 1 helper T cell), M1, (M1 macrophage), DC (Dendritic Cell). (B) Unfavorable immune cells: Th2 (type 2 helper T cell), Treg (regulatory T cell), M2 (M2 macrophage. (C) Cytolytic activity score (CYT). ER+/HER2-, ER-positive and HER2 negative; HER2+, HER2 positive; TNBC, triple-negative breast cancer. p-value  <  0.05 was considered statistically significant.
Ijms 21 05744 g004
Ijms 21 05744 g005 550
Figure 5. GSEA of adipocyte low tumors enriched Hallmark cell cycle and cell proliferation-related gene sets in whole TCGA cohort. NES, normalized enrichment score; FDR, false discovery rate. FDR < 0.25 was considered statistically significant.
Figure 5. GSEA of adipocyte low tumors enriched Hallmark cell cycle and cell proliferation-related gene sets in whole TCGA cohort. NES, normalized enrichment score; FDR, false discovery rate. FDR < 0.25 was considered statistically significant.
Ijms 21 05744 g005
Ijms 21 05744 g006 550
Figure 6. The amount of intratumoral adipocytes correlate with MKI67 expression and Nottingham pathological grade. (A) MKI67 gene expression by intratumoral adipocyte High (shaded box) vs. Low (open box) in whole TCGA, GSE25066, and Yau et al.’s cohort. (B) Pearson correlation curve of MKI67 expression and adipocyte amount in the whole TCGA, GSE25066 and Yau et al.’s cohort. (C) Nottingham pathological grade (grade 1: Open box, grade 2: Shaded box, grade 3: Closed box) and amount of adipocytes in whole cohort and ER+/HER2-, Her2+, TNBC subtype in TCGA and GSE25066 cohort. ER+/HER2-, ER-positive and HER2 negative; HER2+, HER2 positive; TNBC, triple-negative breast cancer. p-value  <  0.05 was considered statistically significant. r represents a Spearman’s correlation coefficient.
Figure 6. The amount of intratumoral adipocytes correlate with MKI67 expression and Nottingham pathological grade. (A) MKI67 gene expression by intratumoral adipocyte High (shaded box) vs. Low (open box) in whole TCGA, GSE25066, and Yau et al.’s cohort. (B) Pearson correlation curve of MKI67 expression and adipocyte amount in the whole TCGA, GSE25066 and Yau et al.’s cohort. (C) Nottingham pathological grade (grade 1: Open box, grade 2: Shaded box, grade 3: Closed box) and amount of adipocytes in whole cohort and ER+/HER2-, Her2+, TNBC subtype in TCGA and GSE25066 cohort. ER+/HER2-, ER-positive and HER2 negative; HER2+, HER2 positive; TNBC, triple-negative breast cancer. p-value  <  0.05 was considered statistically significant. r represents a Spearman’s correlation coefficient.
Ijms 21 05744 g006

Share and Cite

MDPI and ACS Style

Tokumaru, Y.; Oshi, M.; Katsuta, E.; Yan, L.; Huang, J.L.; Nagahashi, M.; Matsuhashi, N.; Futamura, M.; Yoshida, K.; Takabe, K. Intratumoral Adipocyte-High Breast Cancer Enrich for Metastatic and Inflammation-Related Pathways but Associated with Less Cancer Cell Proliferation. Int. J. Mol. Sci. 2020, 21, 5744. https://doi.org/10.3390/ijms21165744

AMA Style

Tokumaru Y, Oshi M, Katsuta E, Yan L, Huang JL, Nagahashi M, Matsuhashi N, Futamura M, Yoshida K, Takabe K. Intratumoral Adipocyte-High Breast Cancer Enrich for Metastatic and Inflammation-Related Pathways but Associated with Less Cancer Cell Proliferation. International Journal of Molecular Sciences. 2020; 21(16):5744. https://doi.org/10.3390/ijms21165744

Chicago/Turabian Style

Tokumaru, Yoshihisa, Masanori Oshi, Eriko Katsuta, Li Yan, Jing Li Huang, Masayuki Nagahashi, Nobuhisa Matsuhashi, Manabu Futamura, Kazuhiro Yoshida, and Kazuaki Takabe. 2020. "Intratumoral Adipocyte-High Breast Cancer Enrich for Metastatic and Inflammation-Related Pathways but Associated with Less Cancer Cell Proliferation" International Journal of Molecular Sciences 21, no. 16: 5744. https://doi.org/10.3390/ijms21165744

APA Style

Tokumaru, Y., Oshi, M., Katsuta, E., Yan, L., Huang, J. L., Nagahashi, M., Matsuhashi, N., Futamura, M., Yoshida, K., & Takabe, K. (2020). Intratumoral Adipocyte-High Breast Cancer Enrich for Metastatic and Inflammation-Related Pathways but Associated with Less Cancer Cell Proliferation. International Journal of Molecular Sciences, 21(16), 5744. https://doi.org/10.3390/ijms21165744

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