Clinical Implications of Tumor-Infiltrating Immune Cells in Breast Cancer.

The immune infiltration of tumors is closely related to clinical outcomes. The composition of tumor-infiltrating immune cells (TIICs) can serve as biomarkers for predicting response to treatment and survival in different patient subgroups in terms of chemotherapy and immunotherapy. This study is focused on investigating the clinical implications of TIICs in breast cancer patients. We performed several in silico analyses of gene expression profiles in 2976 nonmetastatic tumor samples. CIBERSORT was used to estimate the proportion of 22 immune cell types to analyze their correlation with overall survival (OS) and disease-free survival (DFS) in different breast cancer subtypes and stages. Our results showed that a higher fraction of plasma cells in estrogen receptor (ER)-positive breast cancer patients indicated an increase in DFS (hazard ratio [HR]=0.66, 95% confidence interval [CI] 0.54~0.82, p<0.01), while a decreased OS was correlated with a greater number of M0 macrophages (HR=2.02, 95% CI 1.27~3.30, p=0.01) and regulatory T cells (HR=1.90, 95% CI 1.20~3.02, p=0.02). In ER-negative or progesterone receptor (PR)-negative subtypes or in a combined subtype, the increase in activated memory CD4+ T cells was correlated with increased DFS (HR=0.46, 95% CI 0.33~0.63, p<0.01). In all breast cancer patients, a higher proportion of M0 macrophages indicated a decreased DFS (HR=1.67, 95% CI 1.22~2.27, p<0.01), while increased OS was associated with relatively larger fractions of resting memory CD4+ T cells (HR=0.70, 95% CI 0.55~0.90, p=0.02) and γδ T cells (HR=0.66, 95% CI 0.51~0.85, p<0.01). Therefore, this study revealed that the composition of TIICs is different in patients with various subtypes of breast cancer and is directly related to prognosis, suggesting that TIICs are important participants in tumor progression and may, potentially be used for future diagnosis and treatment. © The author(s).

Background

Breast cancer is the most frequent and common malignant tumor and the leading cause of cancer death among females 1. Moreover, it is a heterogeneous disease with high clinical diversity, and thus, the standard treatment regimens are not suitable for all patients 2. In recent years, great advances have been achieved in early diagnosis, surgical techniques, neoadjuvant and adjuvant treatments, but the clinical outcomes are still unsatisfactory 3. Therefore, it is essential to predict treatment responses, select more personalized treatment regimens, and explore new drug targets for different patient subgroups 4.

Tumor-infiltrating immune cells (TIICs) are indispensable components of the tumor microenvironment and have been used for the prediction of prognosis and treatment in cancer patients 4-6. Tumor-infiltrating lymphocytes (TILs) are the most widely studied populations of TIICs and their correlation with favorable outcomes in breast cancer has been reported previously 7-9. In patients with human epidermal growth factor receptor 2 (HER2)-positive tumors, a high number of TILs is associated with increased overall survival (OS) and may predict an improved response to anthracyclines and trastuzumab 9-11. Moreover, the presence of TILs can decrease distal recurrence and increase metastasis-free survival in naive triple-negative breast cancer 11-14. However, the number of TILs is not always proportional to the treatment response, suggesting that the immune phenotypes of TILs are also an important prognostic factor 7. In addition to TILs, tumor-associated macrophages (TAMs) are another immune population that can regulate the interaction between cancer and the immune system. Clinicopathological studies have suggested that TAM accumulation in tumors usually correlates with a poor clinical outcome 15. Furthermore, there have been some inconsistent opinions regarding the effects of TIICs on cancer patients, which may be due to differences in tumor molecular subtypes between tumor serial progression stages 5-7.

Immunotherapies based on immune checkpoints have significantly improved clinical outcomes in a subset of patients with solid tumors, including breast cancer. Recently, cytotoxic T-lymphocyte-associated protein 4 (CTLA4), programmed death-1 (PD-1) and programmed death-ligand 1 (PD-L1) have been shown to be highly promising targets for modulating the interaction of immune cells and tumor cells. However, successful cases have been limited to specific types of cancers, which has been attributed to the insufficient and heterogeneous expression of checkpoint molecules in the tumor microenvironment 16, 17. Moreover, the distribution of TIICs may an important factor in determining whether patients show a good response to anticancer therapies 18, 19. Previous studies concerning the composition of TIICs have relied mainly on flow cytometry and immunohistochemistry, but these methods can only evaluate a few immune cell types at once and are limited by phenotypic markers 20. Interestingly, the algorithmic analysis of a large amount of gene expression data provides a fast and effective method for estimating the composition of TIICs.

In this study, we performed several in silico analyses on the gene expression profiles of 2976 unrelated tumor samples from nonmetastatic breast cancer patients with known clinical follow-up. CIBERSORT 21, 22, an established computational approach, was applied to estimate the relative proportions of 22 types of immune cells in these tumor samples. This is conducive to confirming the cellular composition in terms of the immune response in tumor patients and disclosing the associations among molecular subtypes, clinicopathological variables and survival. In addition, immune checkpoint modulators in breast cancer were analyzed in an attempt to reveal their potential prognostic value.

Materials and methods

Data acquisition

Data from the public domain were used in this study. Raw microarray expression data from primary tumors from patients with nonmetastasized breast cancer were downloaded from the Gene Expression Omnibus (GEO). GEO accession numbers of these data were GSE16391, GSE19615, GSE20685, GSE20713, GSE21653, GSE22035, GSE22646, GSE26639, GSE28583, GSE28694, GSE28826, GSE29044, GSE31192, GSE31448, GSE32646, GSE36771, GSE40837, GSE41568, GSE43502, GSE43568, GSE47389, GSE48390, GSE48391, GSE58984, GSE71258, GSE76124, GSE87007, GSE87377, GSE88770. Clinicopathological data were also collected, including gender, age, TNM stage, grade, Ki67 status, tumor size, the use of adjuvant chemotherapy, OS and DFS, and estrogen receptor (ER), progesterone receptor (PR) and HER2 status. To avoid the influence of differences among platforms, our analyses were confined to samples that hybridized to the Affymetrix HG-U133 Plus 2.0 platform (accession number GPL570).

Estimation of tumor-infiltrating immune cells

CIBERSORT, a deconvolution algorithm, was applied to estimate the abundances of specific cell types in complex tissues based on gene expression profiles, which has been proven to be a method resulting in a high degree of consistency during actual assessment in tumor studies 4, 21, 22. Normalized gene expression data were used to estimate the relative proportions of 22 types of TIICs using CIBERSORT. Simultaneously, the LM22 signature matrix with 1000 permutations was used for the performance of the algorithm. The samples for which the CIBERSORT analysis resulted in p<0.05 were used for further analysis. These TIICs included seven types of T cells (resting memory CD4+ T cells, activated memory CD4+ T cells, CD8+ T cells, naive CD4+ T cells, γδ T cells, follicular helper T cells, and regulatory T cells), memory B cells, naive B cells, activated mast cells, resting mast cells, activated dendritic cells (DCs), resting DCs, macrophages (M0, M1 and M2), eosinophils, monocytes, activated natural killer (NK) cells, resting NK cells, plasma cells and neutrophils. For each sample, the sum of all estimated immune cell-type fractions equalled 1.

Statistical analysis

Statistical analyses were conducted using R v3.3.2 and Bioconductor (https://www.bioconductor.org/). The aggregation and preparation of the raw data were performed using the “affy” package. Each dataset was processed by a weighted average method to compare the differences in the compositions of TIICs in different subtypes and grades of breast cancer. We used a univariate linear regression model to assess the association between the percentage of immune cells and clinicopathological variables. Cox regression analysis was performed to evaluate the prognostic values of TIICs and immune checkpoint molecules. The hazard ratio (HR) and 95% confidence interval (CI) were determined, and a p-value less than 0.05 was considered statistically significant for all statistical analyses.

Results

Data set of breast cancer samples

A data set that included 2976 breast cancer samples was chosen for the analyses in this study. A summary of the primary tumor characteristics of these samples is presented in Table 1. Samples were classified according to their ER, PR and HER2 statuses (Figure 1).

Table 1

Primary tumor characteristics

Variable	No. of sample	Valid%
Age at diagnosis, y
≤50	727	43.59%
>50	941	56.41%
Gender
Female	2166	100%
Male	0	0
Tumor grade
1	142	12.24%
2	473	40.78%
3	545	46.98%
ER status
Positive	938	52.64%
Negative	844	47.36%
HER2 status
Positive	251	22.55%
Negative	862	77.45%
PR status
Positive	741	42.83%
Negative	989	57.17%

Figure 1

Overview of breast cancer subtypes based on receptor status.

Composition of TIICs

The composition of immune cells in breast cancers of different subtypes and grades are shown in Figure 2. Breast tumors generally contained an abundance of TAMs (36.3%), plasma cells (17.9%) and follicular helper T (Tfh) cells (6.9%), whereas eosinophils (0.1%), monocytes (0.5%) and resting NK cells (0.6%) were scarce. The composition of TIICs was greatly different in various subtypes of breast cancer. Compared with HER2-positive tumors, there was a higher proportion of activated mast cells in HER2-negative samples. The number of activated memory CD4+ T cells, activated DCs and naive CD4+ T cells was increased in ER-negative cancers, while more resting mast cells and resting memory CD4+ T cells existed in ER-positive cancers (Figures 2A and 3A). As the tumor grade increased, the proportions of activated DCs, activated memory CD4+ T cells, M0 and M1 macrophages were gradually increased, whereas activated NK cells, monocytes, and resting DCs were decreased in cancer patients with higher grade tumors (Figures 2B and 3B).

Figure 2

Distribution of immune cell-type frequencies in breast cancer of different subtypes (A) and grades (B). The fractions of each immune cell type in different breast cancer subtypes and stages were compared. The size of the bubble represents the fraction of the immune cell-type.

Figure 3

Differences in the composition of TIICs in various subtypes (A) and grades (B) of breast cancer.

Immune cell-type fractions and patient prognosis

To assess the prognostic values of immune cells, a Cox regression analysis was used to analyze the correlations between survival and the proportions of TIICs. As shown in Figure 4, the sizes of the bubbles in the heatmap represent the level of statistical significance, while blue and red colors indicate a negative and positive correlation, respectively, between TIICs and OS/DFS. In breast cancer patients, a higher proportion of M0 macrophages indicated decreased DFS (HR=1.67, 95% CI 1.22~2.27, p<0.01), while increased OS was associated with relatively higher fractions of resting memory CD4+ T cells (HR=0.70, 95% CI 0.55~0.90, p=0.02) and γδ T cells (HR=0.66, 95% CI 0.51~0.85, p<0.01). In ER-positive breast cancer patients, an increased fraction of plasma cells indicated an increase in DFS (HR=0.66, 95% CI=0.54~0.82, p<0.01), while a decreased OS was correlated with an increased number of M0 macrophages (HR=2.02, 95% CI 1.27~3.30, p=0.01) and regulatory T cells (HR=1.90, 95% CI 1.20~3.02, p=0.02). Conversely, an increase in activated memory CD4+ T cells in ER-negative or PR-negative subtypes or in tumors with a combined subtype was related to an increase in DFS (HR=0.46, 95% CI 0.33~0.63, p<0.01). In addition, larger fraction of activated mast cells in HER2-negative tumors was associated with an improved OS (HR=0.30, 95% CI 0.12~0.71, p=0.02), while a higher fraction of activated NK cells in HER2- positive patients were associated with decreased OS (HR=6.25, 95% Cl 1.74~22.49, p<0.02) (Figure 4A).

Figure 4

Bubble heat map showing the predictive and prognostic value of immune cell-type frequencies in breast cancer of different subtypes (A) and grades (B). The associations between the fractions and OS and DFS were assessed by Cox regression analysis. A blue bubble indicates that a higher fraction of immune cells is associated with decreased DFS or OS; a red bubble indicates that a higher fraction of immune cells is associated with increased DFS or OS. The size of the bubble indicates the statistical significance level.

In patients with different stages of cancer, the correlation between TIICs and OS/DFS exhibited large differences (Figure 4B). In grade 2 tumors, increased OS was associated with a relatively higher proportion of resting CD4+ T cells (HR=0.41, 95% CI 0.21~0.79, p=0.03) and a smaller fraction of regulatory T cells (HR=2.40, 95% CI 1.28~4.49, p=0.02), while larger fractions of activated NK cells (HR=2.43, 95% CI 1.39~4.23, p<0.01), eosinophils (HR=3.75, 95% CI 1.38~10.18, p=0.03), and M2 macrophages (HR=1.91, 95% CI 1.13~3.25, p=0.04) were correlated with decreased DFS. In grade 3 tumors, increased OS was related to a smaller fraction of activated NK cells (HR=2.16, 95% CI 1.29~3.61, p=0.01), while increased DFS was related to higher proportions of resting NK cells (HR=0.60, 95% CI 0.42~0.85, p=0.01), naive B cells (HR=0.64, 95% CI 0.45~0.89, p=0.03) and a relatively smaller fraction of Tfh cells (HR=1.52, 95% CI 1.08~2.11, p=0.04).

Association between TIICs and clinicopathological variables

Based on the above analyses, we found that plasma cells, resting memory CD4+ T cells, regulatory T cells, γδ T cells, activated NK cells, and M0 macrophages played relatively important roles in breast cancer. Therefore, a univariate analysis was performed to further reveal the potential correlations between these TIICs and clinicopathological variables including age, TNM stage, Ki67 status, tumor size, and the use of adjuvant chemotherapy. As shown in Figure 5, resting memory CD4+ T cells were associated with a smaller tumor size. M0 macrophages and regulatory T cells were associated with increased proliferation (≥20% Ki67). In patients receiving adjuvant chemotherapy, high numbers of plasma cells and γδ T cells and a smaller fraction of regulatory T cells were observed.

Figure 5

Univariate associations between the composition of TIICs and clinicopathological variables in a breast cancer cohort.

Immunomodulators and tumor subtypes

To evaluate the potential prognostic value of immune checkpoint molecules participating in tumor immune escape, a total of ten immunomodulators were analyzed, including two co-stimulatory molecules (CD27 and ICOS) and eight co-inhibitory molecules (PD-1, TIM-3, LAG3, PD-L1, CTLA4, PD-L2, TIGIT and IDO1) (Figure 6). A poor prognosis was associated with low expression levels of TIM-3, CTLA4, TIGIT and IDO1 in the HER2-positive subtype and with low expression levels of ICOS, TIM-3, PD-L1, CTLA4 and PD-L2 in ER-positive breast cancer patients. Thus, these results indicated that immunotherapies targeting to these immune checkpoint molecules may result in limited or no added value for HER2-negative or ER-negative breast cancer patients.

Figure 6

Cox regression analyses of the prognostic value of immunomodulators in ER-positive, ER-negative, HER2-positive and HER2-negative patients.

Discussion

In this study, we revealed distinct patterns of immune infiltration in various breast cancer patient subgroups based on the deconvolution of bulk gene expression data from a large set of samples. Classically, breast cancer is classified into four subtypes ([ER+|PR+]HER2+, [ER+|PR+]HER2-, ER-PR-HER2+, and ER-PR-HER2-) according to the dichotomized immunohistochemical evaluation of three receptors 23, 24. In our analyses, to obtain much more detailed information concerning the clinical implications of TIICs in breast cancer patients by analyzing ER- and HER2-predominant subtypes, a total of 2976 nonmetastatic breast tumor samples were classified into nine subtypes based on receptor status, including ER+PR+, ER+PR-, ER-PR+, ER-PR-, ER+HER2+, ER+HER2-, ER-HER2+, ER-HER2-, and ER-PR-HER2-. Among them, the proportions of the ER-PR- and ER+PR+ subtypes were greatest, whereas the ER-PR+ subgroup was the smallest, which is consistent with the fact that ER and PR status are strongly associated with each other, since PR-positive status depends on ER-positive status 2, 23. Furthermore, the prognostic significance of TIICs was investigated and the associations between TIICs and clinicopathological parameters were revealed to provide guidance for the development of immune-modulating therapies. We found that the proportion of regulatory T cells was higher in the ER-positive than in the ER-negative subtype, and that an increased fraction of regulatory T cells was related to decreased OS in ER-positive tumors. Regulatory T cells can express inhibitory receptors to maintain self-tolerance and suppress the anti-tumor activity of effector T cells by secreting soluble immunosuppressive cytokines 25, 26. Hence, the anti-CTLA4 antibody (ipilimumab), which can kill effector regulatory T cells or attenuate their suppressive activity, is considered a promising treatment drug 26. Our observations are in line with previous findings that indicate that the regulatory T cells, compared to other immune cell populations, predict statistically significantly worse prognosis for many malignancies 27, 28. In addition, regulatory T cells are decreased after adjuvant chemotherapy, which may profoundly drive the development of tumor-specific CD4+ T effector cells. This provides a valuable insight and may aid in the development of a potential intervention strategy for ER-positive subtypes.

Macrophages can be divided into classical M1 and alternative M2 macrophages based on their function 29. It is generally believed that M1 macrophages participate in inflammatory reactions and anti-tumor immunity, while M2 macrophages have pro-tumorigenic properties. However, the “polarized states” of macrophages are controversial. M1 and M2 phenotypes may represent the extremes of a range of functional states rather than truly different cell types 30. In this study, we observed that a high fraction of M0 macrophages were related to decreased OS and DFS in ER-positive tumors. In terms of the development of interventions that can affect macrophage polarization, our results deserve further attention. Moreover, an increased fraction of M0 macrophages was observed in breast tumors with higher grades, which implied that circulating macrophages can be gradually recruited into tumors to alter the tumor microenvironment and promote tumor progression. Thus, our results further suggested that TAMs could be considered biomarkers for breast cancer.

Compared with the ER-positive subtype, ER-negative tumors had a higher fraction of activated memory CD4+ T cells, which was related to increased DFS. Memory T cells have a long lifespan, which contributes to their key role in tumor development 31. The subsets of memory T cells are CD4+ and CD8+ memory T cells. Memory CD8+ T cells can kill tumor cells upon secondary recognition of the tumor-associated antigens 32. Meanwhile, memory CD4+ T cells suppress the outgrowth of tumor cells by promoting the proliferation of CD8+ cells, which migrate to tumor-associated tissues and differentiate into effector cells 33, 34. Our findings showed that increased DFS was directly correlated with an increase in resting and activated memory CD4+ T cells in breast cancer, indicating the anti-tumor properties of memory CD4+ T cells. Unfortunately, due to the reduction in resting memory CD4+ T cells after adjuvant chemotherapy, this clinical strategy may show little efficacy against ER-positive, PR-positive or combined subtypes.

Our results showed that the presence of an increased fraction of activated NK cells suggested a poor prognosis in HER2-positive tumors, especially for breast cancer patients with advanced stage tumors. Occasionally, NK cells exhibit immunoregulatory functions and can limit T cell responses 35, which may contribute to metabolic reprogramming-induced aberrant effector functions 36, 37. Therefore, the interactions of immune cells in the tumor microenvironment are complicated, and chemotherapy-induced reduction of activated NK cells may be favorable for strengthening the overall antitumor immunity in host. B cells can differentiate into antibody-secreting plasma cells and memory B cells upon antigen activation, contributing to humoral immunity 38. Our study found that a higher fraction of plasma cells was associated with an increased DFS in breast cancer patients regardless of receptor status. Consistently, a relationship between tumor-infiltrating plasma cells and improved survival has been reported in non-small cell lung cancer 39.

Mast cells are considered crucial effector cells of the immune system. However, their role in tumorigenesis is controversial, which may be a result of local stromal conditions 40, 41. It is generally believed that mast cells could stimulate growth, neo-angiogenesis and metastasis in tumors, and therefore they are linked to worse outcome. Conversely, some studies have confirmed that mast cells can mobilize T cells and DCs for the regulation of adaptive T-cell mediated immunity 41. In this study, a higher proportion of activated mast cells was associated with an improved OS in HER2-negative tumors. Thus, the function of mast cells in breast cancer patients may be different depending on the tumor subtype. This indicates that mast cells could serve as a potential therapeutic target in HER2-negative tumors.

Finally, the prognostic value of some important immune checkpoint modulators was revealed. The interactions of immune checkpoint modulators can regulate immune responses. Recently, an attractive strategy for immunotherapy has been developed that involves the suppression of immune tolerance by inhibiting the interaction between tumor cells and immune cells 19, 42. Our observations in this study showed that the expression levels of TIM-3, TIGIT, IDO1, ICOS, CTLA4, PD-L1 and PD-L2 were associated with a favorable prognosis in HER2-positive or ER-positive breast cancer. A similar situation has also been reported in breast cancer and other carcinomas, and this may contribute to the insufficient and heterogeneous expression of checkpoint molecules 17 and may be influenced by age, tumor size, histologic grade, nodal metastasis, hormone receptor status, HER2 status, and the extent of TILs 43, 44. Clinical trials have shown that anti-tumor immunotherapies targeting immune checkpoints have achieved considerable success in both positive and negative cancer patients, and this is an argument for the use of immunomodulator inhibiting drugs in breast cancer patients and also indicates that it may be worthwhile to study the relationship between immune checkpoint modulators in TIICs and cancer prognosis.

Conclusion

In this study, we revealed the complex relationships between infiltrating immune cells and tumor subtypes, clinicopathological parameters and survival in breast cancer. The results showed that a multitude of immune cells were associated with therapeutic responses and prognosis. Moreover, there were large changes in TIIC composition in various subtypes, suggesting the complex involvement of the immune system in breast cancer patients. Based on our analyses, it is worth determining the proportions of M0 macrophages, resting memory CD4+ T cells and γδ T cells in all breast cancer patients. Plasma cells, M0 macrophages and regulatory T cells are valuable in ER-positive subtypes, while activated memory CD4+ T cells are favorable in terms of DFS in ER- or PR-negative patients or those with a combined subtype. Activated mast cells play a critical role in HER2-negative tumors, but activated NK cells have a significant adverse effect on HER2-positive tumors, as well as patients in the middle and terminal stage of breast cancer. Overall, this study can promote the understanding of immune phenotypes and provide valuable guidance to optimize immunotherapeutic regimens or combined chemoimmunotherapy strategies to achieve improved outcomes for breast cancer patients that are not responsive to existing treatments.