- 1Department of Nuclear Medicine, Seoul National University College of Medicine, Seoul, South Korea
- 2Department of Molecular Medicine and Biopharmaceutical Sciences, Graduate School of Convergence Science and Technology, Seoul National University, Seoul, South Korea
- 3Department of Thoracic and Cardiovascular Surgery, Seoul National University Hospital, Seoul, South Korea
Introduction: Tumor immune microenvironment (TIME) promotes immune escape, allowing for tumor progression and metastasis. In spite of the current evidence of the complicated role of immune cells in promoting or suppressing cancer progression, the heterogeneity of TIME according to the tumor site has been scarcely investigated. Here, we analyzed transcriptomic profiles of metastatic breast cancer to understand how TIME varies according to tumor sites.
Methods: Two gene expression datasets from metastatic breast cancer of various sites and a single-cell RNA sequencing dataset of primary breast cancer and metastatic lymph nodes were analyzed. The immune cell-type enrichment of each tumor was estimated. Immune cell types were identified by clustering analysis, and the proportions of cell types in TIME were assessed according to the tumor site.
Results: Metastatic bone lesions showed more neutrophils than breast lesions. Tumors clustered according to immune cell type were significantly associated with tumor site. In single-cell analyses, the TIMEs of metastatic lymph nodes showed fewer macrophages than those of primary tumors. Differentially expressed gene signatures in the primary tumor and metastatic lymph nodes were associated with macrophage activation.
Conclusion: We conclude that metastatic sites show variable enrichment patterns of immune cells, and that the TIME of metastatic lesions should be considered in precise immuno-oncology treatments.
Introduction
The tumor immune microenvironment (TIME) plays a key role in tumor progression and metastasis. Interactions between immune cells, stromal cells, and chemokines establish a pre-metastatic niche for circulating tumor cells (1–3). The cell types in TIME play different pro- and antitumor roles during tumor progression (4). For instance, neutrophils are recruited to the metastatic site to initiate metastasis by activating dormant tumor cells and enhancing migration of tumor cells (5, 6). Tumor-associated macrophages (TAMs), a specific subtype of macrophages, facilitate tumor progression to support invasion of tumor cells and suppress the tumoricidal immune response (7–9). In contrast, CD8+ T-cells and natural killer cells inhibit tumor growth by promoting antitumor immunity (10, 11).
The roles of immune regulation in different metastatic sites have recently been elucidated (12). A metastatic niche is established in the lymph nodes by suppressing autoreactive immunity mediated by cytokines such as IL-17 and regulatory T-cells that support immune evasion of tumor cells (13, 14). In the brain, one of the most common organs of hematogenous metastasis, metastasis is associated with different cellular populations such as astrocytes, microglia, and macrophages, which play a role in cancer cell infiltration and survival (15, 16). Protecting cancer stem cells from immune cells is an important mechanism promoting metastasis in the bone (17). Despite recent findings of specific cellular populations in TIME and the close interaction between cancer cells and immune cells, to the best of our knowledge, there has been no explorative interrogation of cell enrichment within TIME according to the metastatic site. Notably, as the cellular composition of TIME can vary across tumor types and is regarded as a key predictive biomarker for immunotherapy (18, 19), a comprehensive understanding of TIME according to the metastatic site has a clinical implication with regard to tumor-specific immunotherapy (20).
Therefore, this study aimed to assess the difference in the cellular composition of TIME in different metastatic sites. We used multiple datasets, including microarray and single-cell RNA-sequencing (scRNA-seq), from breast cancer patients to assess enrichment and proportions of immune cells based on the metastatic site. In addition, the proportions of cell types in TIME were assessed in the scRNA-seq dataset. Differential gene expression analysis was conducted to compare primary breast tumors and metastatic lymph nodes.
Methods
Data Sources
We used four publicly available datasets in this study. For bulk tissue analyses, we obtained two microarray datasets of metastatic breast cancer tissues (accession number GSE124647 and GSE56493) from the Gene Expression Omnibus database (https://www.ncbi.nlm.nih.gov/geo/) (21, 22). The normalized gene expression data and the clinical and demographic characteristics of both cohorts were downloaded using the “GEOquery” package in R (23). Subjects with metastatic lesions in “lung/pleura” lesions were excluded due to uncertainty whether metastatic lesions located in a lung or pleura. Subjects with metastatic lesions in “others” were excluded due to a lack of data for a specific location. Subjects with ovary or soft tissue were excluded, because they were included in only one dataset (GSE124647), not in another dataset (GSE56493). A total of 102 and 117 subjects were included for the analysis in GSE124647 and GSE56493, respectively. The specific sites of breast cancer metastasis are summarized in Supplementary Table 1. For single-cell analysis, we obtained an scRNA-seq dataset of breast cancer patients from the Gene Expression Omnibus database (accession number GSE75688) (24). It included 11 breast cancer patients. Additionally, an scRNA-seq dataset of head and neck cancer (accession number GSE103322) was analyzed to validate the results of the breast cancer dataset (25). It included 17 head and neck cancer patients.
Immune Cell Enrichment Analysis
To evaluate the cellular landscape of TIME, cell type enrichment scores and proportions of cellular populations were estimated. We used three available bioinformatics analytical tools that provide an estimation of the cell type enrichment scores or relative levels of distinct cell types from gene expression data: xCell, CIBERSORT, and TIMER (26–28). We estimated the cell type enrichment scores of 64 immune and stromal cell types of tumor microenvironment using “xCellAnalysis” function included in the “xCell” package in R. Subsequently, composite scores of 35 immune cells were selected for the analysis. Using CIBERSORT, we obtained the relative proportion of 22 immune cell types in TIME. The CIBERSORT analysis process was performed with web-based resource (https://cibersort.stanford.edu/). TIMER provided the enrichment score reflecting the relative abundance of 6 immune cell types. The TIMER analysis process was performed with web-based resource (http://timer.cistrome.org/).
scRNA-Seq Analysis
A read count matrix was generated from an scRNA-seq dataset of breast cancer patients (GSE75688). The scRNA-seq data were scaled by log-normalization after the read counts were divided by the total number of transcripts and multiplied by 10,000. Two thousand highly variable genes were selected using the “FindVariableFeatures” function of Seurat (version 3.0) (29, 30). Data were then scaled to z-scores with regression analysis from total cellular read counts and mitochondrial read counts. Cell types were determined using the graph-based clustering approach implemented by the “FindClusters” function. Before clustering, dimension reduction was performed by principal component analysis, and 10 dimensions were used for the clustering. The conservative resolution was set to 1.0. To identify the marker genes of the clusters, the “FindAllMarkers” function of Seurat was used, and 10 high-ranked marker genes were identified according to the fold-change. The scRNA-seq data were embedded by two-dimensional projection using t-distributed stochastic neighborhood embedding. To identify cell types, the expression levels of 8 markers (CD68, CD79A, CD3D, CD8A, CD4, FOXP3, COL1A1, and CDH1) were assessed. CD3D, CD8A, CD4, and FOXP3 are well-known markers for tumor-infiltrating lymphocytes (31). CD68 and CD79A are the most representative markers of macrophages and B-cells, respectively (32, 33). CDH1 and COL1A1 are markers for epithelium and stroma (34, 35). Based on the expression levels, each cluster was classified into 4 cell types: cancer cell, T-cell, B-cell, and macrophage. Cell numbers were compared between the primary tumor and metastatic lymph nodes in the entire dataset and in the paired samples. Subsequently, differentially expressed genes (DEGs) were extracted by Wilcoxon rank sum test using the “FindMarkers” function. The threshold of adjusted p-value was 0.05 to determine DEGs. Genes with adjusted p-value below 0.1 were employed in gene ontology (GO) analysis. It was performed to assess biological processes related to the DEGs with the enrichGO function. The cutoff value for p-value and q-value was 0.05.
Statistics
The cell type enrichment scores estimated from the primary tumor and other metastatic sites were compared using the Wilcoxon rank test. Log2 fold change (log2FC) was calculated from mean of enrichment scores of each cell type. A threshold of log2FC was 0.5 to determine differential expression, and a p-value of 0.05 was considered significant in volcano plots showing differences in cell type enrichment scores or relative levels between the metastatic site and breast lesion. Hierarchical clustering was performed in Euclidean distances between each sample with Ward's minimum variance method. Chi-square analyses were applied to identify association between clusters and metastatic sites. All statistical analyses were performed using R (v 3.6.1). The overall scheme of analysis methods is described in Supplementary Figure 1.
Results
Differences in Cell Types in TIME as Per the Tumor Site
We first evaluated the enrichment of cell types in TIME using a gene expression dataset from breast cancer patients (GSE124647). As shown in the heatmap of the cell enrichment analysis, immune cells were highly different in different samples (Figure 1A). We also evaluated the effect of metastatic tumor sites on the cell enrichment scores/proportions. To identify specific cell types enriched in metastatic tumors, cell enrichment scores/proportions of each immune cell from each metastatic site were compared with those from breast lesions (Figures 1B–D). In CIBERSORT analysis, fewer macrophages were observed in metastatic lymph nodes than in breast lesions. In xCell and CIBERSORT analysis, more neutrophils were found in metastatic bone lesions than in breast lesions. Log2FC and p-value were listed in Supplementary Table 2.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/649004/fonc-11-649004-HTML/image_m/fonc-11-649004-g001.jpg)
Figure 1. Enrichment of immune cell populations in different metastatic sites. (A) Heatmaps depicting the distribution of immune cell enrichment scores according to metastatic site estimated using xCell, CIBERSORT, and TIMER. The distribution of immune cells varied across all samples. (B–D) Volcano plots showing enriched immune cells in specific metastatic sites compared to breast lesions. Wilcoxon rank test was applied to compare the enrichment of immune cells. A threshold of log2FC was 0.5 to determine differential expression, and a p-value of 0.05 was considered significant. Values outside the range of x-axis were displayed at the margin of plots. Results from the xCell, CIBERSORT, and TIMER analyses are shown in (B–D), respectively. In CIBERSORT analysis, there were fewer macrophages in metastatic lymph nodes than in breast lesions. In xCell and CIBERSORT analyses, there were more neutrophils in metastatic bone lesions than in breast lesions.
Different Myeloid Cell Compositions as Per the Tumor Site
The previous analysis showed that macrophages were less enriched in some metastatic tumors compared to the primary site, especially in xCell and CIBERSORT analyses. Therefore, we focused on the enrichment level of macrophages in metastatic lymph nodes and bone lesions. Metastatic lymph nodes had significantly fewer macrophages than breast lesions in the CIBERSORT analysis (p = 0.004, W = 608). There was no significant difference in xCell and TIMER analyses (Figure 2A). A lower macrophage enrichment score was also observed in metastatic bone lesions. Metastatic bone lesions had fewer macrophages than breast lesions (Figure 2B). The xCell analysis showed a significantly lower macrophage enrichment score in metastatic bone lesions than in breast tumors (p = 0.012, W = 46). There was no significant difference in CIBERSORT and TIMER analyses. Conversely, there were more neutrophils, another myeloid cell type, in metastatic bone lesions than in breast lesions in all 3 analyses (Figure 2C, p = 0.016, W = 161 for xCell; p = 0.050, W = 150.5 for CIBERSORT; and p = 0.005, W = 169 for TIMER).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/649004/fonc-11-649004-HTML/image_m/fonc-11-649004-g002.jpg)
Figure 2. Different compositions of myeloid cells according to metastatic tumor site. (A) In the CIBERSORT analysis, macrophages showed significantly lower enrichment scores in metastatic lymph nodes than in breast lesions (p = 0.0039, CI = 0.04~0.19). (B) In the xCell analysis, macrophages showed significantly lower enrichment scores in metastatic bone lesions than in breast lesions (p = 0.012, CI = −8.38*10−2~-6.06*10−5). (C) In the xCell and TIMER analyses, neutrophils showed significantly higher enrichment scores in metastatic bone lesions than in breast lesions (p = 0.016, CI = 0.03~0.23 and p = 0.0045, CI = 4.83*10−3~0.02, respectively).
Tumors were then clustered according to TIME cell types to identify any association between TIME-based clusters and metastatic sites. Three clusters were identified via hierarchical clustering of TIME cell enrichment scores. In the xCell analysis, cluster 1 was associated with enriched macrophages, plasma cells, and dendrocytes. Cluster 2 was associated with enriched monocytes, neutrophils, and eosinophils. Cluster 3 was associated with enriched B-cells and CD4+ T-cells. The clusters based on TIME were significantly associated with metastatic tumor sites, particularly bone and lymph node lesions, in all 3 analyses (Figures 3A–C; p < 0.001, chi-square = 30.818 for xCell; p < 0.001, chi-square = 30.937 for CIBERSORT; and p < 0.001, chi-square = 31.338 for TIMER, Supplementary Table 3). Most of metastatic bone lesions were included in cluster 2. Although metastatic lymph nodes were included in all clusters, cluster 3 contained only metastatic lymph nodes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/649004/fonc-11-649004-HTML/image_m/fonc-11-649004-g003.jpg)
Figure 3. Association between TIME clusters and tumor sites. Mosaic plots showing the association of specific clusters with specific metastatic sites. Heatmaps showing top 3 immune cells enriched in each cluster. For example, macrophages showed high enrichment in cluster 1, neutrophils in cluster 2, and B-cells in cluster 3 in the xCell analysis. (A) Results from xCell. (B) Results from CIBERSORT. (C) Results from TIMER. TIME, tumor immune microenvironment.
We subsequently examined whether this finding was recapitulated in another dataset (GSE56493). As in the previous analyses, heatmaps and volcano plots were created (Supplementary Figure 2). Metastatic lymph nodes showed significantly fewer macrophages than breast lesions in the xCell analysis (p = 0.034, W = 560). This was also observed in the TIMER analysis (Supplementary Figure 3A, p < 0.001, W = 663). Metastatic bone lesions contained significantly fewer macrophages than breast lesions in both the xCell and TIMER analyses (Supplementary Figure 3B; p = 0.003, W = 6 for xCell and p = 0.024, W = 16 for TIMER). Neutrophils were significantly enriched in bone metastases compared to breast lesions in the CIBERSORT analysis (Supplementary Figure 3C, p = 0.015, W = 81). The clusters based on TIME were significantly associated with lymph node lesions in all 3 analyses (Supplementary Figures 3D–F; p < 0.039, chi-square = 21.847 for xCell; p < 0.001, chi-square 27.693 for CIBERSORT; and p < 0.001, chi-square = 32.06 for TIMER, Supplementary Table 3).
Immune Cell Types in Metastatic Lymph Nodes in scRNA-Seq
To elucidate the differences in cell types in TIME according to the metastatic site, we analyzed single-cell-level data of metastatic lymph nodes from breast cancer (GSE75688). A total of 515 cells from 11 tumor samples were clustered into 10 clusters and 4 cell types (Figure 4A). The markers of each cell cluster and their expression are shown in Supplementary Figure 4. The TIMEs of metastatic lymph nodes had fewer macrophages than those of primary tumors (Figures 4B,C). There were more B-cells and T-cells in metastatic lymph nodes than in primary breast tumors, which was consistent with our results derived from bulk gene expression data. In addition, we analyzed whether cancer cells in different tumor sites had different gene expression, which could be affected by different TIMEs. When cancer cells from primary and metastatic tumors were compared, seven DEGs were identified: small EDRK-rich factor 2 (SERF2), calmodulin-like protein 5 (CALML5), fatty acid binding protein 7 (FABP7), homeodomain-only protein (HOPX), androgen receptor (AR), CD9, and phospholipase A2 group IIA (PLA2G2A) (Figures 4D,E). Functional enrichment analysis using GO analysis demonstrated that macrophage-related terms were significantly enriched in this transcriptional profile (Figure 4F). There were two sets of paired primary tumors and metastatic lymph nodes from the same patient, and we performed additional analyses in these paired samples. In these samples, the TIMEs of the metastatic lymph nodes had fewer macrophages than those of the primary tumors, as observed in the total dataset (Supplementary Figures 5A,B). We identified nine DEGs in these paired samples: SERF2, CALML5, FABP7, HOPX, PLA2G2A, sperm protein associated with the nucleus on the X chromosome D (SPANXD), endoplasmic reticulum oxidoreductase 1 alpha (ERO1L), AR, and CD9 (Supplementary Figures 5C,D). These genes included seven of the DEGs identified in all primary and metastatic tumors. In GO analyses, macrophage-related terms were also selected as significantly enriched (Supplementary Figure 5E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/649004/fonc-11-649004-HTML/image_m/fonc-11-649004-g004.jpg)
Figure 4. TIME of metastatic lymph nodes analyzed by single-cell RNA-sequencing. (A) A total of 515 cells from an scRNA-seq dataset clustered into 10 clusters and 4 cell types. (B) There were fewer macrophages in metastatic lymph nodes than in primary tumors in the t-SNE plot. (C) Bar plot showing macrophages in metastatic lymph nodes and primary tumors. (D) Volcano plots representing seven differentially expressed genes: SERF2, CALML5, FABP7, HOPX, AR, CD9, and PLA2G2A. All the seven DEGs showed negative log2 fold change; lower expression in primary tumors than in metastatic lymph nodes. (E) Heatmap demonstrating the expression of each differentially expressed gene. (F) In gene ontology analyses, macrophage-related terms were selected as significant. (G) In a head and neck cancer scRNA-seq dataset, there were fewer macrophages in metastatic lymph nodes than in primary tumors in the t-SNE plot. (H) Bar plot showing macrophages in metastatic lymph nodes and primary tumors. AR, androgen receptor; CALML5, calmodulin-like protein 5; FABP7, fatty acid binding protein 7; HOPX, homeodomain-only protein; PLA2G2A, phospholipase A2 group IIA; scRNA-seq, single-cell RNA sequencing; SERF2, small EDRK-rich factor 2; TIME, tumor immune microenvironment; t-SNE, t-distributed stochastic neighborhood embedding.
To validate our finding of the difference in the enrichment of macrophages in TIME, we performed additional analyses using single-cell-level transcriptome data derived from head and neck cancer (GSE103322). Even though the cancer type was different from that in the other datasets, head and neck cancer also commonly metastasizes to the lymph nodes. Therefore, it was expected to show similar results, i.e., fewer macrophages and more lymphocytes in metastatic lymph nodes than in the primary tumor. A total of 5782 cells from 22 tumors were clustered into 11 clusters and 4 cell types (Supplementary Figure 6A). The markers of each cell cluster and their expression are shown in Supplementary Figure 6B. The TIMEs of metastatic lymph nodes had fewer macrophages than those of primary tumors (Figures 4G,H). In six paired samples, the TIMEs of metastatic lymph nodes had fewer macrophages than those of primary tumors, as in the total dataset (Supplementary Figures 6C,D).
Discussion
The present study focused on differences in TIME based on metastatic sites of breast cancer. We explored the TIMEs of different metastatic sites using bulk gene expression datasets. This was followed by single-cell analysis to support differences in the TIMEs of metastatic lymph nodes, which was revealed in the analysis of the bulk datasets. The TIME showed a high variance according to the metastatic site. Myeloid cells were significantly enriched in metastatic lymph nodes and metastatic bone lesions compared to breast lesions. This was further shown by the cluster analysis, as metastatic lymph nodes and bone lesions were classified into specific clusters. Myeloid cells, including macrophages and neutrophils, are closely associated with tumorigenesis and metastasis, although they have antitumor roles. However, their protumor roles are more common, such as promoting tissue remodeling and suppressing innate immunity to tumor cells (36–40). Additionally, they may interact with each other via cytokines such as interlukin-8 or tumor necrosis factor alpha (41, 42). Although macrophages and neutrophils cooperate and have functional similarity, our study showed that the distribution of myeloid cells can vary in different metastatic tissues. This suggests that the pre-metastatic niche and TIME formation may employ different mechanisms in different metastatic sites.
It is unclear whether macrophages have protumoral or antitumoral effects. One of the difficulties in identifying the role of macrophages in metastasis is that there are various subtypes of macrophages and myeloid cells affecting TIME. Previous studies found that the abundance or immaturity of dendritic cells in lymph nodes is related to a lower antitumor effect of T-cells (43, 44). The effects of fewer macrophages in metastatic lymph nodes on metastasis may be inferred from these previous reports. It can be hypothesized that a smaller number of TAMs causes less antigen processing/presenting and decreased T-cell responses in the lymph nodes. Conversely, another subtype of macrophages induces protumoral effects by suppressing the T-cell response via the cytotoxic T-lymphocyte-associated protein 4 and programmed death-ligand 1 (PD-L1) pathways (45). Therefore, there is controversy surrounding the reason for reduced macrophages in metastatic lymph nodes. Scholars revealed that numbers of macrophages were smaller in normal lymph nodes than in primary breast tumor (46). However, another study demonstrated that macrophages in metastatic lymph nodes were lesser in normal lymph nodes, by immunohistochemistry analysis in breast cancer patients (47). Considering these previous studies, the finding that there are fewer macrophages in metastatic lymph nodes than in primary lesions may be deemed not only a result of normal immune cell distribution in lymph nodes but also feature of disease progression. This can provide an insight into the underlying biology of metastasis. For cancer cells to survive in the metastatic site, they may change the enrichment of specific immune cell types in TIME. The decrease in macrophages may result from this phenomenon.
Results from the bulk datasets were supported by scRNA-seq data of breast cancer and head and neck cancer. In a part of samples, metastatic lymph nodes showed lesser macrophages than primary tumors. Our results imply that the immune populations in TIME may be, at least partly, dependent on the anatomical site. Our results highlight that the TIMEs of specific metastatic sites may be similar, even if the primary site differs. This is supported by the concept of the pre-metastatic niche, wherein the environment of the metastatic site is altered to proper “soil” for the tumor “seed” (48).
Interestingly, there were fewer macrophages in bone metastases than in breast lesions in our datasets. However, neutrophils were enriched in metastatic bone lesions compared to breast lesions. Therefore, antigen presentation is not well-activated in bone lesions as macrophages and dendritic cells which showed lower enrichment in the bone lesions compared with breast lesions. These types of cells were associated with antigen presentation of tumors; a key process of the tumor immunity mediated by lymphocytes (49, 50). Instead, it seems that neutrophils may be a key factor in bone metastasis. Neutrophils release cytokines, including chemokine receptor type 4 and matrix metallopeptidase 9, to foster bone invasion (5). Additionally, previous studies showed that the degree of neutrophil recruitment was associated with prognosis in patients with bone metastasis (51, 52). Although the lung is another representative metastatic site of breast cancer, characteristics of TIME in metastatic lung lesions were not shown in this study due to a lack of data. However, previous animal studies have shown that TAMs and neutrophils promote metastasis of breast cancer (53, 54). These knowledges also support that TIME may vary according to the metastatic sites.
It is notable that FABP7 and AR were selected as differentially expressed genes which were upregulated in cancer cells in metastatic lymph nodes. FABP7 is known as a molecular factor affecting brain metastasis and survival/proliferation of breast cancer (55, 56). Also, AR is known to promote proliferation and migration of breast cancer (57). As previous knowledges, the present study shows potential relevance of these genes toward promoting lymph node metastasis of breast cancer, by analyzing scRNA-seq data.
Precision oncology is based on the unique genomic or molecular characteristics of each patient and tumor (58, 59). Immuno-oncology therapy can enhance the host immunity or inhibit the protumor effect of tumor-infiltrated immune cells (19, 60). Immune checkpoint inhibitor (ICI) therapies, including PD-L1 checkpoint inhibitors, are the most commonly applied immunotherapies in the clinic. ICIs have caused unprecedented prognostic improvement in lung cancer, lymphoma, and melanoma patients (61–63). However, several recent studies have reported that the response of ICI therapy differs based on the metastatic site. It has been postulated that the different immunologic environment of each tissue may affect the response to ICIs (64–66). This would indicate that the concept of precision medicine needs to move from an individual-based approach to a lesion-based approach. In other words, it is important to investigate the TIME of each metastatic lesion or site to select an appropriate treatment option, thereby improving prognosis. The presence of TAMs in TIME might be of particular interest, as targeting TAMs is thought to affect the response to ICI therapy (45). There are experimental lines of evidence that suppressing TAMs enhances the efficacy of ICIs and inhibits cancer progression (67–69).
Our findings provide insight into the clinical translation of treatment strategies that consider differences in TIME based on the metastatic site. There have been few studies documenting the diverse microenvironments of different metastatic tissues (70). This is the first study to demonstrate differences in TIME based on the metastatic site in breast cancer. We clearly showed differences in the TIME components in metastatic lymph nodes and bone lesions compared to breast lesions. Further, this finding was validated using a single-cell dataset and data from another cancer. Our data also highlight the necessity of different therapies according to the metastatic site, even in the same individual. However, we were unable to include patients with metastatic lesions in more than one anatomical site. Nevertheless, our results can be used to select treatment options based on metastatic site.
In this study, we applied multiple analysis tools to conduct quantitative analysis of immune populations in TIME. There were a few dissimilarities between three methods. In addition, some results were not completely replicated in the same analysis methods. It is caused by different deconvolutional methods and marker genes in each tool. Also, different demographics and clinical settings of each dataset may cause partly heterogeneous results. Nonetheless, the key findings were supported by similar tendency despite a degree of statistical insignificance between the three methods. It is a limitation of this study focusing on analyses for different pre-existing datasets. Further study is warranted to corroborate these findings with experimental methods.
In conclusion, we analyzed 2 gene expression datasets and one scRNA-seq dataset to analyze the TIMEs in different metastatic breast cancer sites. The enrichment of cell types in TIME differed in different metastatic sites. Specifically, there were fewer macrophages in metastatic lymph nodes and metastatic bone lesions than in breast lesions. Our findings suggest that the TIME of each metastatic lesion must be considered when selecting the optimal immuno-oncology treatment option.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found at: https://www.ncbi.nlm.nih.gov/geo/, GSE124647, GSE56493, GSE75688, and GSE103322.
Author Contributions
HL, KN, and HC designed the study, performed analysis and interpretation, and wrote the manuscript. All authors read and approved the final manuscript.
Funding
This research was supported by the National Research Foundation of Korea (NRF) and funded by the Korean government (MSIT) (Nos. 2020R1C1C1007105 and 2020M3A9B6038086), as well as the Korea Health Technology R&D Project through the Korea Health Industry Development Institute (KHIDI), and by the Ministry of Health & Welfare, Republic of Korea (No. HI19C0339).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.649004/full#supplementary-material
References
1. Yuan Y, Jiang Y-C, Sun C-K, Chen Q-M. Role of the tumor microenvironment in tumor progression and the clinical applications. Oncol Rep. (2016) 35:2499–515. doi: 10.3892/or.2016.4660
2. Gonzalez H, Hagerling C, Werb Z. Roles of the immune system in cancer: from tumor initiation to metastatic progression. Genes Dev. (2018) 32:1267–84. doi: 10.1101/gad.314617.118
3. El-Kenawi A, Hänggi K, Ruffell B. The immune microenvironment and cancer metastasis. Cold Spring Harb Perspect Med. (2020) 10:a037424. doi: 10.1101/cshperspect.a037424
4. Quail DF, Joyce JA. Microenvironmental regulation of tumor progression and metastasis. Nat Med. (2013) 19:1423–37. doi: 10.1038/nm.3394
5. Leach J, Morton JP, Sansom OJ. Neutrophils: homing in on the myeloid mechanisms of metastasis. Mol Immunol. (2019) 110:69–76. doi: 10.1016/j.molimm.2017.12.013
6. Wu L, Saxena S, Awaji M, Singh RK. Tumor-associated neutrophils in cancer: going pro. Cancers. (2019) 11:564. doi: 10.3390/cancers11040564
7. Lin Y, Xu J, Lan H. Tumor-associated macrophages in tumor metastasis: biological roles and clinical therapeutic applications. J Hematol Oncol Rep. (2019) 12:76. doi: 10.1186/s13045-019-0760-3
8. Sanchez LR, Borriello L, Entenberg D, Condeelis JS, Oktay MH, Karagiannis GS. The emerging roles of macrophages in cancer metastasis and response to chemotherapy. J Leukoc Biol. (2019) 106:259–74. doi: 10.1002/JLB.MR0218-056RR
9. Vitale I, Manic G, Coussens LM, Kroemer G, Galluzzi L. Macrophages and metabolism in the tumor microenvironment. Cell Metab. (2019) 30:36–50. doi: 10.1016/j.cmet.2019.06.001
10. Shanker A, Verdeil G, Buferne M, Inderberg-Suso E-M, Puthier D, Joly F, et al. CD8 T cell help for innate antitumor immunity. J Immunol. (2007) 179:6651–62. doi: 10.4049/jimmunol.179.10.6651
11. Woan KV, Miller JS. Harnessing natural killer cell antitumor immunity: from the bench to bedside. Cancer Immunol Res. (2019) 7:1742–7. doi: 10.1158/2326-6066.CIR-19-0404
12. Peinado H, Zhang H, Matei IR, Costa-Silva B, Hoshino A, Rodrigues G, et al. Pre-metastatic niches: organ-specific homes for metastases. Nat Rev Cancer. (2017) 17:302–17. doi: 10.1038/nrc.2017.6
13. Pereira ER, Jones D, Jung K, Padera TP. The lymph node microenvironment and its role in the progression of metastatic cancer. Semin Cell Dev Biol. (2015) 38:98–105. doi: 10.1016/j.semcdb.2015.01.008
14. Jones D, Pereira ER, Padera TP. Growth and immune evasion of lymph node metastasis. Front Oncol. (2018) 8:36. doi: 10.3389/fonc.2018.00036
15. You H, Baluszek S, Kaminska B. Immune microenvironment of brain metastases-are microglia little helpers? Front Immunol. (2019) 10:1941. doi: 10.3389/fimmu.2019.01941
16. Cacho-Díaz B, García-Botello DR, Wegman-Ostrosky T, Reyes-Soto G, Ortiz-Sánchez E, Herrera-Montalvo LA. Tumor microenvironment differences between primary tumor and brain metastases. J Transl Med. (2020) 18:1. doi: 10.1186/s12967-019-02189-8
17. Roato I, Ferracini R. Cancer stem cells, bone and tumor microenvironment: key players in bone metastases. Cancers. (2018) 10:56. doi: 10.3390/cancers10020056
18. Finotello F, Eduati F. Multi-omics profiling of the tumor microenvironment: paving the way to precision immuno-oncology. Front Oncol. (2018) 8:430. doi: 10.3389/fonc.2018.00430
19. Marshall HT, Djamgoz M. Immuno-oncology: emerging targets and combination therapies. Front Oncol. (2018) 8:315. doi: 10.3389/fonc.2018.00315
20. Reuben A, Spencer CN, Prieto PA, Gopalakrishnan V, Reddy SM, Miller JP, et al. Genomic and immune heterogeneity are associated with differential responses to therapy in melanoma. NPJ Genom Med. (2017) 2:10.
21. Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, et al. NCBI GEO: archive for functional genomics data sets—update. Nucleic Acids Res. (2012) 41:D991–5. doi: 10.1093/nar/gks1193
22. Sinn BV, Fu C, Lau R, Litton J, Tsai T-H, Murthy R, et al. SET ER/PR: a robust 18-gene predictor for sensitivity to endocrine therapy for metastatic breast cancer. NPJ Breast Cancer. (2019) 5:16. doi: 10.1038/s41523-019-0111-0
23. Davis S, Meltzer PS. GEOquery: a bridge between the gene expression omnibus (GEO) and bioconductor. Bioinformatics. (2007) 23:1846–7. doi: 10.1093/bioinformatics/btm254
24. Chung W, Eum HH, Lee HO, Lee KM, Lee HB, Kim KT, et al. Single-cell RNA-seq enables comprehensive tumour and immune cell profiling in primary breast cancer. Nat Commun. (2017) 8:15081. doi: 10.1038/ncomms15081
25. Puram SV, Tirosh I, Parikh AS, Patel AP, Yizhak K, Gillespie S, et al. Single-cell transcriptomic analysis of primary and metastatic tumor ecosystems in head and neck cancer. Cell. (2017) 171:1611–24. doi: 10.1016/j.cell.2017.10.044
26. Li B, Severson E, Pignon J-C, Zhao H, Li T, Novak J, et al. Comprehensive analyses of tumor immunity: implications for cancer immunotherapy. Genome Biol. (2016) 17:174. doi: 10.1186/s13059-016-1028-7
27. Aran D, Hu Z, Butte AJ. xCell: digitally portraying the tissue cellular heterogeneity landscape. Genome Biol. (2017) 18:220. doi: 10.1186/s13059-017-1349-1
28. Chen B, Khodadoust MS, Liu CL, Newman AM, Alizadeh AA. Profiling tumor infiltrating immune cells with CIBERSORT. Methods Mol Biol. (2018) 1711:243–59. doi: 10.1007/978-1-4939-7493-1_12
29. Butler A, Hoffman P, Smibert P, Papalexi E, Satija R. Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat Biotechnol. (2018) 36:411–20. doi: 10.1038/nbt.4096
30. Stuart T, Butler A, Hoffman P, Hafemeister C, Papalexi E, Mauck III, et al. Comprehensive integration of single-cell data. Cell. (2019) 177:1888–902. doi: 10.1016/j.cell.2019.05.031
31. Gooden MJM, Bock GH, Leffers N, Daemen T, Nijman HW. The prognostic influence of tumour-infiltrating lymphocytes in cancer: a systematic review with meta-analysis. Br J Cancer. (2011) 105:93–103. doi: 10.1038/bjc.2011.189
32. Holness CL, Simmons DL. Molecular cloning of CD68, a human macrophage marker related to lysosomal glycoproteins. Blood. (1993) 81:1607–13. doi: 10.1182/blood.V81.6.1607.1607
33. Chu PG, Arber DA. CD79: a review. Appl Immunohistochem Mol Morphol. (2001) 9:97–106. doi: 10.1097/00129039-200106000-00001
34. Gall TMH, Frampton AE. Gene of the month: E-cadherin (CDH1). J Clin Pathol. (2013) 66:928–32. doi: 10.1136/jclinpath-2013-201768
35. Wang Y, Xu H, Zhu B, Qui Z, Lin Z. Systematic identification of the key candidate genes in breast cancer stroma. Cell Mol Biol Lett. (2018) 23:44. doi: 10.1186/s11658-018-0110-4
36. Mantovani A, Schioppa T, Porta C, Allavena P, Sica A. Role of tumor-associated macrophages in tumor progression and invasion. Cancer Metastasis Rev. (2006) 25:315–22. doi: 10.1007/s10555-006-9001-7
37. Fridlender ZG, Albelda SM. Tumor-associated neutrophils: friend or foe? Carcinogenesis. (2012) 33:949–55. doi: 10.1093/carcin/bgs123
38. Galdiero MR, Bonavita E, Barajon I, Garlanda C, Mantovani A, Jaillon S. Tumor associated macrophages and neutrophils in cancer. Immunobiology. (2013) 218:1402–10. doi: 10.1016/j.imbio.2013.06.003
39. Galdiero MR, Garlanda C, Jaillon S, Marone G, Mantovani A. Tumor associated macrophages and neutrophils in tumor progression. J Cell Physiol. (2013) 228:1404–12. doi: 10.1002/jcp.24260
40. Masucci MT, Minopoli M, Carriero MV. Tumor associated neutrophils. Their role in tumorigenesis, metastasis, prognosis and therapy. Front Oncol. (2019) 9:1146. doi: 10.3389/fonc.2019.01146
41. Lefkowitz DL, Lefkowitz SS. Macrophage–neutrophil interaction: a paradigm for chronic inflammation revisited. Immunol Cell Biol. (2001) 79:502–6. doi: 10.1046/j.1440-1711.2001.01020.x
42. Kumar V, Sharma A. Neutrophils: cinderella of innate immune system. Int Immunopharmacol. (2010) 10:1325–34. doi: 10.1016/j.intimp.2010.08.012
43. Huang RR, Wen DR, Guo J, Giuliano AE, Nguyen M, Offodile R, et al. Selective modulation of paracortical dendritic cells and T-lymphocytes in breast cancer sentinel lymph nodes. Breast J. (2000) 6:225–32. doi: 10.1046/j.1524-4741.2000.98114.x
44. Ito M, Minamiya Y, Kawai H, Saito S, Saito H, Nakagawa T, et al. Tumor-derived TGFβ-1 induces dendritic cell apoptosis in the sentinel lymph node. J Immunol. (2006) 176:5637–43. doi: 10.4049/jimmunol.176.9.5637
45. Cassetta L, Kitamura T. Targeting tumor-associated macrophages as a potential strategy to enhance the response to immune checkpoint inhibitors. Front Cell Dev Biol. (2018) 6:38. doi: 10.3389/fcell.2018.00038
46. Azizi E, Carr AJ, Plitas G, Cornish AE, Konopacki C, Prabhakaran S, et al. Single-cell map of diverse immune phenotypes in the breast tumor microenvironment. Cell. (2018) 174:1293–308. doi: 10.1016/j.cell.2018.05.060
47. Yang J, Li X, Liu X, Liu Y. The role of tumor-associated macrophages in breast carcinoma invasion and metastasis. Int J Clin Exp Pathol. (2015) 8:6656–64.
48. Paget S. The distribution of secondary growths in cancer of the breast. Cancer Metastasis Rev. (1989) 8:98–101.
49. Boczkowski D, Nair SK, Nam J-H, Lyerly HK, Gilboa E. Induction of tumor immunity and cytotoxic T lymphocyte responses using dendritic cells transfected with messenger RNA amplified from tumor cells. Cancer Res. (2000) 60:1028–34.
50. DeNardo DG, Andreu P, Coussens LM. Interactions between lymphocytes and myeloid cells regulate pro-versus anti-tumor immunity. Cancer Metastasis Rev. (2010) 29:309–16. doi: 10.1007/s10555-010-9223-6
51. Wang S, Zhang Z, Fang F, Gao X, Sun W, Liu H. The neutrophil/lymphocyte ratio is an independent prognostic indicator in patients with bone metastasis. Oncol Lett. (2011) 2:735–40. doi: 10.3892/ol.2011.304
52. Thio QC, Goudriaan WA, Janssen SJ, Pereira NRP, Sciubba DM, Rosovksy RP, et al. Prognostic role of neutrophil-to-lymphocyte ratio and platelet-to-lymphocyte ratio in patients with bone metastases. Br J Cancer. (2018) 119:737–43. doi: 10.1038/s41416-018-0231-6
53. DeNardo DG, Barreto JB, Andreu P, Vasquez L, Tawfik D, Kolhatkar N, et al. CD4+ T cells regulate pulmonary metastasis of mammary carcinomas by enhancing protumor properties of macrophages. Cancer Cell. (2009) 16:91–102. doi: 10.1016/j.ccr.2009.06.018
54. Wculek SK, Malanchi I. Neutrophils support lung colonization of metastasis-initiating breast cancer cells. Nature. (2015) 528:413–7. doi: 10.1038/nature16140
55. Liu RZ, Graham K, Glubrecht DD, Lai R, Mackey JR, Godbout R. A fatty acid-binding protein 7/RXRβ pathway enhances survival and proliferation in triple-negative breast cancer. J Pathol. (2012) 228:310–21. doi: 10.1002/path.4001
56. Cordero A, Kanojia D, Miska J, Panek WK, Xiao A, Han Y, et al. FABP7 is a key metabolic regulator in HER2+ breast cancer brain metastasis. Oncogene. (2019) 38:6445–60. doi: 10.1038/s41388-019-0893-4
57. Giovannelli P, Di Donato M, Galasso G, Di Zazzo E, Bilancio A, Migliaccio A. The androgen receptor in breast cancer. Front Endocrinol. (2018) 9:492. doi: 10.3389/fendo.2018.00492
58. Garraway LA, Verweij J, Ballman KV. Precision oncology: an overview. J Clin Oncol. (2013) 31:1803–5. doi: 10.1200/JCO.2013.49.4799
59. Collins FS, Varmus H. A new initiative on precision medicine. N Eng J Med. (2015) 372:793–5. doi: 10.1056/NEJMp1500523
60. Antonia SJ, Larkin J, Ascierto PA. Immuno-oncology combinations: a review of clinical experience and future prospects. Clin Cancer Res. (2014) 20:6258–68. doi: 10.1158/1078-0432.CCR-14-1457
61. Wolchok JD, Kluger H, Callahan MK, Postow MA, Rizvi NA, Lesokhin AM, et al. Nivolumab plus ipilimumab in advanced melanoma. N Eng J Med. (2013) 369:122–33. doi: 10.1056/NEJMoa1302369
62. Ansell SM, Lesokhin AM, Borrello I, Halwani A, Scott EC, Gutierrez M, et al. PD-1 blockade with nivolumab in relapsed or refractory Hodgkin's lymphoma. N Eng J Med. (2015) 372:311–9. doi: 10.1056/NEJMoa1411087
63. Reck M, Rodríguez-Abreu D, Robinson AG, Hui R, Csoszi T, Fülöp A, et al. Pembrolizumab versus chemotherapy for PD-L1–positive non–small-cell lung cancer. N Eng J Med. (2016) 375:1823–33. doi: 10.1056/NEJMoa1606774
64. Nishino M, Ramaiya NH, Chambers ES, Adeni AE, Hatabu H, Jänne PA, et al. Immune-related response assessment during PD-1 inhibitor therapy in advanced non-small-cell lung cancer patients. J Immunother Cancer. (2016) 4:84. doi: 10.1186/s40425-016-0193-2
65. Schmid S, Diem S, Li Q, Krapf M, Flatz L, Leschka S, et al. Organ-specific response to nivolumab in patients with non-small cell lung cancer (NSCLC). Cancer Immunol Immunother. (2018) 67:1825–32. doi: 10.1007/s00262-018-2239-4
66. Osorio JC, Arbour KC, Le DT, Durham JN, Plodkowski AJ, Halpenny DF, et al. Lesion-level response dynamics to programmed cell death protein (PD-1) blockade. J Clin Oncol. (2019) 37:3546–55. doi: 10.1200/JCO.19.00709
67. Zhu Y, Knolhoff BL, Meyer MA, Nywening TM, West BL, Luo J, et al. CSF1/CSF1R blockade reprograms tumor-infiltrating macrophages and improves response to T-cell checkpoint immunotherapy in pancreatic cancer models. Cancer Res. (2014) 74:5057–69. doi: 10.1158/0008-5472.CAN-13-3723
68. Georgoudaki AM, Prokopec KE, Boura VF, Hellqvist E, Sohn S, Östling J, et al. Reprogramming tumor-associated macrophages by antibody targeting inhibits cancer progression and metastasis. Cell Rep. (2016) 15:2000–11. doi: 10.1016/j.celrep.2016.04.084
69. Steggerda SM, Bennett MK, Chen J, Emberley E, Huang T, Janes JR, et al. Inhibition of arginase by CB-1158 blocks myeloid cell-mediated immune suppression in the tumor microenvironment. J Immunother Cancer. (2017) 5:101. doi: 10.1186/s40425-017-0308-4
Keywords: cancer, metastasis, tumor immune microenvironment, macrophage, neutrophil
Citation: Lee H, Na KJ and Choi H (2021) Differences in Tumor Immune Microenvironment in Metastatic Sites of Breast Cancer. Front. Oncol. 11:649004. doi: 10.3389/fonc.2021.649004
Received: 03 January 2021; Accepted: 22 February 2021;
Published: 18 March 2021.
Edited by:
Benjamin Frey, University Hospital Erlangen, GermanyReviewed by:
Stephen J. Bush, University of Oxford, United KingdomBernardo Cacho-Díaz, National Institute of Cancerology (INCAN), Mexico
Copyright © 2021 Lee, Na and Choi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hongyoon Choi, chy1000@snu.ac.kr; Kwon Joong Na, kjna85@gmail.com
†These authors have contributed equally to this work