- Open Access
Classification of triple-negative breast cancers based on Immunogenomic profiling
Journal of Experimental & Clinical Cancer Researchvolume 37, Article number: 327 (2018)
Abundant evidence shows that triple-negative breast cancer (TNBC) is heterogeneous, and many efforts have been devoted to identifying TNBC subtypes on the basis of genomic profiling. However, few studies have explored the classification of TNBC specifically based on immune signatures that may facilitate the optimal stratification of TNBC patients responsive to immunotherapy.
Using four publicly available TNBC genomics datasets, we classified TNBC on the basis of the immunogenomic profiling of 29 immune signatures. Unsupervised and supervised machine learning methods were used to perform the classification.
We identified three TNBC subtypes that we named Immunity High (Immunity_H), Immunity Medium (Immunity_M), and Immunity Low (Immunity_L) and demonstrated that this classification was reliable and predictable by analyzing multiple different datasets. Immunity_H was characterized by greater immune cell infiltration and anti-tumor immune activities, as well as better survival prognosis compared to the other subtypes. Besides the immune signatures, some cancer-associated pathways were hyperactivated in Immunity_H, including apoptosis, calcium signaling, MAPK signaling, PI3K–Akt signaling, and RAS signaling. In contrast, Immunity_L presented depressed immune signatures and increased activation of cell cycle, Hippo signaling, DNA replication, mismatch repair, cell adhesion molecule binding, spliceosome, adherens junction function, pyrimidine metabolism, glycosylphosphatidylinositol (GPI)-anchor biosynthesis, and RNA polymerase pathways. Furthermore, we identified a gene co-expression subnetwork centered around five transcription factor (TF) genes (CORO1A, STAT4, BCL11B, ZNF831, and EOMES) specifically significant in the Immunity_H subtype and a subnetwork centered around two TF genes (IRF8 and SPI1) characteristic of the Immunity_L subtype.
The identification of TNBC subtypes based on immune signatures has potential clinical implications for TNBC treatment.
Triple-negative breast cancer (TNBC) is a breast cancer subtype that lacks the expression of hormone receptors (estrogen receptor (ER) and progesterone receptor (PR)) and human epidermal growth factor receptor 2 (HER2). TNBC is associated with a high risk of mortality for its aggressiveness and the lack of effective targeted therapies. Moreover, abundant evidence shows that TNBC is very heterogeneous [1,2,3,4]. Lehmann et al. identified six gene expression profile-based TNBC subtypes, including an immunomodulatory (IM) subtype that was enriched in immune cell processes . Bonsang-Kitzis et al. identified six TNBC subgroups based on a biological network-driven approach, which included two immunity clusters whose stromal immune module gene signatures exhibited a strong prognostic value . Burstein et al. identified four stable TNBC subgroups based on mRNA expression and DNA genomic profiling, which included Luminal/Androgen Receptor, Mesenchymal, Basal-Like Immune Suppressed, and Basal-Like Immune Activated (BLIA); furthermore, the authors identified potential therapeutic targets for these specific subtypes . These efforts to classify TNBC might lay the foundation for developing targeted therapies for TNBC.
Recently, cancer immunotherapy has been successful in treating many refractory malignancies . Thus, it is worth considering immunotherapy for TNBC, since the therapeutic options for this disease are significantly limited. Indeed, many experimental and clinical studies have explored the possibility of treating TNBC patients with immunotherapy [6,7,8,9,10,11]. Moreover, numerous studies have demonstrated that TNBC is more immunogenic than other breast cancer (BC) subtypes, which may warrant an immunotherapeutic approach for TNBC [12, 13]. However, currently, immunotherapeutic strategies exhibit beneficial effects in less than 20% of cancer patients. This suggests that not all TNBC patients could respond to immunotherapy. In fact, certain genetic or genomic features, such as tumor mutation burden (TMB), neoantigen load, PD-L1 expression, and deficient DNA mismatch repair, have been associated with cancer immunotherapeutic responsiveness [14,15,16,17,18].
In this study, we classified TNBC into three distinct subtypes by immunogenomic profiling: Immunity High (Immunity_H), Immunity Medium (Immunity_M), and Immunity Low (Immunity_L). We demonstrated the stability and reproducibility of this classification in four independent datasets by a machine learning approach. Furthermore, we identified the subtype-specific molecular features, including genes, gene ontology, pathways, and networks. The identification of immune signature-associated TNBC subtypes may facilitate the optimal selection of TNBC patients responsive to immunotherapy.
For each TNBC dataset, we first quantified the enrichment levels of the 29 immune signatures in each TNBC sample by the single-sample gene-set enrichment analysis (ssGSEA) score [19, 20]. Based on the enrichment levels (ssGSEA scores) of the 29 immune signatures, we performed hierarchical clustering of TNBC.
Evaluation of immune cell infiltration level, tumor purity, and stromal content in TNBC
ESTIMATE  was used to evaluate the immune cell infiltration level (immune score), tumor purity, and stromal content (stromal score) for each TNBC sample.
Gene-set enrichment analysis
We performed gene-set enrichment analysis of the METABRIC and TCGA datasets by GSEA (R implementation) [22,23,24]. This analysis identified the KEGG  pathways that were upregulated in Immunity_H and Immunity_L (FDR < 0.05), respectively. The common pathways identified in both datasets were selected.
Correlation of pathway activities with immune cell infiltration levels in TNBC
We quantified the activity of a pathway with the ssGSEA score of the set of genes included in the pathway, and the immune cell infiltration level with the immune score. The Spearman correlation of the ssGSEA score and the immune score were used to evaluate the correlation of pathway activities with immune cell infiltration levels in TNBC.
Identification of TNBC subtype-specific gene ontology and networks
We used WGCNA  to identify the gene modules (gene ontology) that were significantly associated with the genes highly correlated with immune cell infiltration based on gene co-expression analysis. The gene modules specifically amplified in different TNBC subtypes were identified. On the basis of the expression correlations between the hub genes in the gene modules, we built gene–gene interaction networks. A hub gene was defined as a gene that was connected to no less than 10 other genes, with a connectedness weight greater than 0.25.
We compared the survival prognosis (overall survival (OS), disease-free survival (DFS), and metastasis-free survival (MFS) of TNBC patients considering tumor subtype and the expression level of the identified genes, i.e., higher expression level (expression levels > median) versus lower expression level (expression levels < median). The log-rank test was used to calculate the significance of survival time differences using a threshold of P-value < 0.05. Kaplan–Meier curves were plotted to show the survival time differences. We performed the survival analyses using the METABRIC, TCGA, and GSE103091 datasets, where the survival data were available.
We transformed each attribute (immune signature or gene set) value (ssGSEA score) xi into xi′ by the equation xi′ = (xi − xmin)/(xmax − xmin), where xmin and xmax represent the minimum and maximum of the ssGSEA scores for the gene set across all TNBC samples, respectively. The Random Forest (RF) classifier was used to classify the TNBC subtypes. We set the number of trees to 100 and all 29 immune signatures as features for the RF classifier. The classification performance was evaluated by the accuracy and the weighted F-score. We carried out the classification in Weka .
Comparison of the proportions of immune cell subsets between TNBC subtypes
CIBERSORT  was used to calculate the proportions of 22 human immune cell subsets. We set 1000 permutations and P < 0.05 as the criteria for the successful deconvolution of a sample. We compared the proportions of the immune cell subsets between TNBC subtypes using the Mann–Whitney U test.
Comparison of clonal heterogeneity between the TNBC subtypes
We used the ABSOLUTE algorithm  to assess the ploidy score, representing clonal heterogeneity, for each TNBC sample. We compared the ploidy scores between the TNBC subtypes using the Kruskal–Wallis test.
Comparison of biological processes between the TNBC subtypes
We compared the activities (ssGSEA scores) of stem cell-associated (marker genes ABCA8 and ALDH1A1), proliferation (MKI67), and epithelial-to-mesenchymal transition (EMT) (ZEB1, ZEB2, SNAIL, CDH2 and TGFB1) biological processes between the TNBC subtypes. The Kruskal–Wallis test was used to determine the statistical significance of the results.
Comparison of somatic copy number alteration (SCNA) levels between the TNBC subtypes
We applied GISTIC2  to the SNP6 file of the SCNA data for TNBC in TCGA. We obtained arm-level SCNA frequencies for Immunity_H and Immunity_L TNBC samples and compared them. Moreover, we calculated focal SCNA levels for each TNBC samples and compared them between Immunity_H and Immunity_L.
Immunogenomic profiling identifies three TNBC subtypes
We analyzed 29 immune-associated gene sets which represented diverse immune cell types, functions, and pathways (Additional file 1: Table S1). We used the ssGSEA score [19, 20] to quantify the activity or enrichment levels of immune cells, functions, or pathways in the cancer samples. On the basis of the ssGSEA scores of the 29 gene sets, we hierarchically clustered TNBC in four BC datasets (METABRIC , TCGA , GSE75688 , and GSE103091 ). Interestingly, all four datasets showed similar clustering results, with three clusters being clearly separated (Fig. 1). We defined the three clusters as: Immunity High (Immunity_H), Immunity Medium (Immunity_M), and Immunity Low (Immunity_L). We found that the immune scores were significantly higher in Immunity_H and significantly lower in Immunity_L in all four datasets (Kruskal–Wallis test, P < 0.001) (Fig. 2a). Moreover, we found that the percentage of lymphocyte infiltration was significantly higher in Immunity_H and significantly lower in Immunity_L in TCGA based on the pathological slides data (Kruskal–Wallis test, P = 0.05). These features directed the classification. In addition, when comparing the tumor purity and stromal score of the three TNBC subtypes, we obtained opposite trends, with tumor purity increasing from Immunity_H to Immunity_L (Immunity_H < Immunity_M < Immunity_L) and stromal score decreasing from Immunity_H to Immunity_L (Immunity_H > Immunity_M > Immunity_L) (Kruskal–Wallis test, P < 0.001) (Additional file 2: Figure S1). Collectively, these results suggest that Immunity_H contains the highest number of immune cells and stromal cells, while Immunity_L contains the highest number of tumor cells.
Notably, most HLA genes showed significantly higher expression levels in Immunity_H and significantly lower expression levels in Immunity_L (ANOVA test, P < 0.05) (Fig. 2b, Additional file 3: Figure S2A). Moreover, the expression levels of various immune cell subpopulation marker genes  were the highest in Immunity_H and the lowest in Immunity_L, such as CD8A (cytotoxic T cell), CD45RO (memory T cell), CD20 (B cell), CXCR5 (Tfh cell), FOXP3 (Treg), IL-17 (Th17 cell), CD1A (iDC), and IL3RA (pDC) (Additional file 3: Figure S2B).
We examined the expression of PD-L1 (programmed cell death 1 ligand) in the three TNBC subtypes and found that Immunity_H had the highest PD-L1 expression levels and Immunity_L had the lowest PD-L1 expression levels (ANOVA test, P < 0.05) (Fig. 2c). This suggest that the TNBC subtype Immunity_H might better respond to anti-PD-L1 immunotherapy than the other TNBC subtypes, since PD-L1 expression tends to be positively associated with immunotherapeutic responsiveness .
Survival analyses showed that these TNBC subtypes had distinct clinical outcomes. The Immunity_H subtype likely had a better survival prognosis than the Immunity_M and Immunity_L subtypes, but there was no significant survival difference between the Immunity_M and the Immunity_L subtypes (Fig. 2d). This is consistent with previous studies showing that TNBC with elevated immune activity were associated with more favorable clinical outcomes [4, 12, 34].
Comparisons of the immunogenomic profiling-based TNBC classification with other TNBC classification methods
We used the TNBCtype method [2, 37] to classify the four TNBC datasets. We found that the immunomodulatory (IM) subtype of TNBCs was most frequently associated with Immunity_H and least frequently associated with Immunity_L (Fisher’s exact test, P < 0.001) (Fig. 3). This is consistent with the enrichment of immune cell processes in the IM subtype . In contrast, the mesenchymal (M) subtype of TNBCs was mostly detected in Immunity_L and least frequently detected in Immunity_H (Fisher’s exact test, P < 0.001) (Fig. 3). The M subtype is mainly characterized by pathways involved in cell motility, ECM receptor interaction, and cell differentiation, such as Wnt, ALK, and TGF-β signaling . Our results suggest that the activities of these pathways may be associated with reduced tumor immunity in TNBC.
Identification of TNBC subtype-specific pathways, gene ontology, and networks
Identification of TNBC subtype-specific pathways
GSEA identified a number of KEGG  pathways enriched in Immunity_H and Immunity_L (Fig. 4a, Additional file 4: Figure S3A). Typically, the immune-associated pathways were highly active in Immunity_H and included antigen processing and presentation pathways, B and T cell receptor signaling, chemokine signaling, cytokine–cytokine receptor interactions, IL-17 signaling, Jak–STAT signaling, natural killer cell-mediated cytotoxicity, NF-kappa B signaling, NOD-like receptor signaling, TNF signaling, and Toll-like receptor signaling (Fig. 4a, Additional file 4: Figure S3A). This result confirmed the elevated immune activity in Immunity_H. Besides, we identified various cancer-associated pathways that were hyperactivated in Immunity_H, including apoptosis, calcium signaling, MAPK signaling, PI3K–Akt signaling, and RAS signaling (Fig. 4a, Additional file 4: Figure S3A). This suggests that the activities of these cancer-associated pathways are positively associated with TNBC immunity. In contrast, the TNBC subtype Immunity_L was enriched in pathways related to Hippo signaling, DNA replication, mismatch repair, spliceosome, adherens junctions, pyrimidine metabolism, glycosylphosphatidylinositol (GPI)-anchor biosynthesis, and RNA polymerase (Fig. 4a, Additional file 4: Figure S3A). This indicates that the activities of these pathways could be negatively associated with TNBC immunity. In fact, a previous study has shown that the activities of MAPK and PI3K–Akt cascades positively correlated with the activation of various immune pathways, while the activity of the mismatch repair pathway showed a negative correlation with immune activation in TNBC . Furthermore, we confirmed that all the cancer-associated pathways hyperactivated in Immunity_H were positively associated with the immune scores, whereas the pathways hyperactivated in Immunity_L likely showed a negative correlation (Spearman’s correlation test, P < 0.05) (Fig. 4b).
Identification of TNBC subtype-specific gene ontology
We performed a weighted gene co-expression network analysis of the METRABRIC dataset by WGCNA  and identified a set of gene modules (gene ontology) associated with the highly expressed genes previously determined. We found several gene modules that significantly differentiated TNBC by subtype, survival time, or survival status (Fig. 4c). As expected, the immune response was significantly elevated in Immunity_H (P = 4.0*10− 54), while was depressed in Immunity_L (P = 1.0*10− 32). Moreover, a high immune response was associated with a better survival prognosis in TNBC patients (P = 5.0*10− 4). This finding is in line with the previous observation that the subtype Immunity_H is associated with better clinical outcomes than the other subtypes. Similar results were observed for the TCGA dataset (Additional file 4: Figure S3B). The other two immune-associated gene modules, i.e., myeloid leukocyte activation and response to type I interferon, were also enriched in Immunity_H (P = 3.0*10− 14 and 2.0*10− 13, respectively), and were reduced in Immunity_L (P = 2.0*10− 16 and 2.0*10− 11, respectively). In contrast, cell adhesion molecule (CAMD) binding activity was significantly increased in Immunity_L (P = 1.0*10− 30) and decreased in Immunity_H (P = 2.0*10− 35). This suggests that CAMD activity has a strong inverse correlation with tumor immunity in TNBC. Interestingly, CAMD activity correlated with reduced survival (P = 0.001 for OS, and P = 0.002 for DFS). Cell cycle process was also increased in Immunity_L (P = 0.04), suggesting that the cell cycle signature correlates with reduced tumor immunity. This finding is consistent with results from previous studies [38, 39].
Identification of TNBC subtype-specific networks
WGCNA generated a gene module (green color, Fig. 4c) that was specifically significant in Immunity_H. We identified 98 hub genes from the gene module, including five transcription factor (TF) genes, i.e., CORO1A, STAT4, BCL11B, ZNF831, and EOMES. The five TFs interact with each other and form a subnetwork with diverse immune and cancer-related genes that they regulate (Fig. 4d). Typically, CD247 (the marker gene for a T cell subpopulation) was regulated by all these TFs, and the cytotoxic T cell marker gene CD8A was co-regulated by CORO1A, STAT4, and EOMES. MAP4K1 (Mitogen-Activated Protein Kinase Kinase Kinase Kinase 1), which is involved in multiple immune and cancer-related pathways including B cell receptor signaling, JNK, EGF/EGFR, TGF-β, and MAPK signaling, was also regulated by the five TFs. CORO1A encodes a member of the WD repeat protein family which is involved in diverse cellular processes including cell cycle, apoptosis, signal transduction, and gene regulation. The main pathways related to CORO1A include cytoskeletal signaling and phagosome function, and its relatedness with immune regulation has been revealed [40, 41]. The association of the other TFs STAT4 , BCL11B , and EOMES  with immunity has been examined, whereas the role of ZNF831 in immune regulation remains unexplored.
WGCNA also generated a gene module (turquoise color, Fig. 4c) that was more enriched in Immunity_L. This module included 112 hub genes, two of which encode the TFs IRF8 and SPI1. A subnetwork of the hub genes centered on IRF8 and SPI1is shown in Fig. 4e. IRF8 (interferon regulatory factor 8) has been shown to play a negative role in immune cell regulation . Thus, the IRF8-centered regulatory network may be responsible for the depressed immunity of the TNBC subtype Immunity_L. SPI1 (Spi-1 proto-oncogene) encodes a transcription factor that activates gene expression during immune cell development. As a result, the deregulation of SPI1 may affect immunity. In fact, SPI1 showed significantly lower expression levels in Immunity_L than in Immunity_H (Student’s t test, P = 9.1*10− 28, fold change > 2). Therefore, the down regulation of SPI1 may contribute to the decreased immunity of the Immunity_L subtype. The contribution of the IRF8- and SPI1-centered regulatory network to the depressed immunity of Immunity_L is evidenced by a previous study showing that IRF8 and SPI1 together negatively regulated immune cell differentiation .
Interestingly, survival analyses showed that elevated expression levels of these TF genes (except SPI11) were consistently associated with better survival prognosis in TNBC (Fig. 4f), suggesting the pivotal role of these TFs in TNBC immunity and prognosis.
Class prediction of TNBC subtypes based on immunogenomic profiling
We first used 10-fold cross validation (CV) to evaluate the classification performance in METABRIC and then predicted the TNBC subtypes in the other three datasets using the METABRIC dataset as the training set. The 10-fold CV accuracy was 89% in classifying the METABRIC dataset. The classification accuracies were 70, 84, and 63% in TCGA, GSE75688, and GSE103091, respectively. The weighted F-scores in these classifications were 89, 71, 83, and 63% for METABRIC, TCGA, GSE75688, and GSE103091, respectively (Fig. 5). These results demonstrate that the immunogenomic profiling-based classification of TNBC is stable and predictable.
A number of prior studies have identified TNBC subtypes on the basis of genomic profiling [2,3,4, 34]. However, very few studies have investigated the classification of TNBC specifically based on immune signatures. To fill this knowledge gap, we focused on identifying immune-related TNBC subtypes using immunogenomic profiling. Our results show that TNBC could be classified into three stable subtypes: Immunity High, Immunity Medium, and Immunity Low. Furthermore, we demonstrated that this classification was reproducible and predictable. The Immunity High TNBC subtype was enriched not only in immune signatures, but also in many cancer-associated pathways including apoptosis, calcium signaling, MAPK signaling, PI3K–Akt signaling, and RAS signaling (Fig. 4a). This is in line with our previous study showing that diverse immune signatures positively correlated with the MAPK and PI3K–Akt signaling pathways in TNBC . In contrast, the Immunity Low TNBC subtype was impoverished in immune signatures but enriched in Hippo signaling, DNA replication, mismatch repair, spliceosome, adherens junction, pyrimidine metabolism, glycosylphosphatidylinositol (GPI)-anchor biosynthesis, and RNA polymerase pathways (Fig. 4a). It is rational that the mismatch repair pathway activity was significantly negatively correlated with immune signatures in cancer, since deficient mismatch repair often results in elevated tumor immunity . Interestingly, we found that the Hippo signaling pathway had a significantly negative correlation with immune signatures in TNBC. This observation is in agreement with findings from previous studies showing that the Hippo signaling pathway plays a key role in regulating tumor immunity [46,47,48]. Deficiency of Hippo pathway components such as kinases LATS1/2 (large tumor suppressor 1 and 2) , effector YAP (Yes-associated protein) , and transcriptional co-activator TAZ (WW domain-containing transcription regulator 1)  could promote anti-tumor immunity. Overall, these results revealed potential positive or negative associations between pathway activities and immune activities in TNBC.
Immunity_H had stronger immune cell infiltration and anti-tumor immune activities, e.g., high levels of cytotoxic T cells and B cells infiltration (Fig. 1). When we used CIBERSORT  to calculate the proportions of 22 immune cell subsets in TNBC, we found that CD8 T cells, CD4 T cells, NK cells, and M1 macrophages tended to be present in significantly higher numbers in Immunity_H than in Immunity_L (Mann–Whitney U test, P < 0.05) (Fig. 6). This analysis again confirmed elevated anti-tumor immune activity in Immunity_H. The high anti-tumor immune activation could explain why Immunity_H had more favorable clinical outcomes compared to the other subtypes (Fig. 2d). In fact, numerous studies have demonstrated that the density of tumor-infiltrating lymphocytes (TILs) is positively associated with survival prognosis in various cancers [35, 49, 50]. Immunity_H more highly expressed most HLA genes, which is indicative of stronger immunogenicity compared to the other subtypes. However, Immunity_H did not show higher TMB or neoantigen load than the other subtypes. This suggests that the differential immunogenicity between the TNBC subtypes cannot be attributed to TMB and neoantigens. In addition, we did not find lower clonal heterogeneity in Immunity_H, as estimated by ABSOLUTE , than in the other subtypes, although, in some studies, clonal heterogeneity was shown to have a significant negative correlation with tumor immunity [13, 51]. Interestingly, Immunity_H exhibited more active stem cell-associated biological processes than the other subtypes (Additional file 5: Figure S4), while it showed no significant differences in proliferation and epithelial-to-mesenchymal transition (EMT) processes.
Furthermore, we compared SCNA levels between Immunity_H and Immunity_L subtypes. We found that Immunity_H had significantly lower arm-level SCNAs than Immunity_L (Wilcoxon signed-rank test, P = 0.04, 0.001, 0.0006 for comparisons of amplification, deletion, and total alteration frequencies, respectively) (Fig. 7a). Moreover, Immunity_H had significantly lower focal SCNA levels than Immunity_L (Mann–Whitney U test, P = 0.01, 0.02, 0.01 for comparisons of amplification, deletion, and total alteration levels, respectively) (Fig. 7b). These findings demonstrated that Immunity_H had lower levels of SCNAs compared to Immunity_L, supporting the notion that high tumor aneuploidy correlates with reduced tumor immune infiltration .
Currently, immunotherapy for TNBC is an active field of investigation , and the stronger immunogenicity exhibited by TNBC compared to other breast cancer subtypes suggests that immunotherapy could be a viable option for TNBC patients . However, some preliminary TNBC immunotherapy clinical trials have not shown significant patients’ improvement (personal communication). Thus, the immune signature-based classification of TNBC may aid the stratification of TNBC patients to identify those responsive to immunotherapy. It is imaginable that patients with an Immunity_H subtype of TNBC would be more likely to respond to anti-PD-1/PD-L1 therapy than patients with other TNBC subtypes, since PD-L1 is more highly expressed in Immunity_H TNBC, and PD-L1 expression is a predictive biomarker for the response to PD-1/PD-L1-directed immunotherapy [36, 54].
The identification of TNBC subtypes based on immune signatures has potential clinical implications for TNBC treatment.
Activated dendritic cells
Basal-like immune activated
Cell adhesion molecule binding
Cytokine and cytokine receptor
False discovery rate
Gene-set enrichment analysis
Human leukocyte antigen
Human epidermal growth factor receptor 2
Immature dendritic cells
Major histocompatibility complex
- NK cells:
Natural killer cells
Plasmacytoid dendritic cells
Single-sample gene-set enrichment analysis
Somatic copy number alteration
The Cancer Genome Atlas
- Tfh cells:
Follicular helper T cells
- Th17 cells:
T helper 17 cells
Triple-negative breast cancer
Regulatory T cells
Weighted gene co-expression network analysis
Lehmann BD, et al. Refinement of triple-negative breast Cancer molecular subtypes: implications for neoadjuvant chemotherapy selection. PLoS One. 2016;11(6):e0157368.
Lehmann BD, et al. Identification of human triple-negative breast cancer subtypes and preclinical models for selection of targeted therapies. J Clin Invest. 2011;121(7):2750–67.
Bonsang-Kitzis H, et al. Biological network-driven gene selection identifies a stromal immune module as a key determinant of triple-negative breast carcinoma prognosis. Oncoimmunology. 2016;5(1):e1061176.
Burstein MD, et al. Comprehensive genomic analysis identifies novel subtypes and targets of triple-negative breast cancer. Clin Cancer Res. 2015;21(7):1688–98.
Del Paggio JC. Immunotherapy: Cancer immunotherapy and the value of cure. Nat Rev Clin Oncol. 2018;15(5):268–270.
Salatino M, Girotti MR, Rabinovich GA. Glycans pave the way for immunotherapy in triple-negative breast Cancer. Cancer Cell. 2018;33(2):155–7.
Liu L, et al. Combination immunotherapy of MUC1 mRNA Nano-vaccine and CTLA-4 blockade effectively inhibits growth of triple negative breast Cancer. Mol Ther. 2018;26(1):45–55.
Hu Z, et al. Targeting tissue factor for immunotherapy of triple-negative breast Cancer using a second-generation ICON. Cancer Immunol Res. 2018;6(6):671–84.
Choi AH, et al. Novel chimeric parapoxvirus CF189 as an oncolytic immunotherapy in triple-negative breast cancer. Surgery. 2018;163(2):336–42.
Bernier C, et al. DZ-2384 has a superior preclinical profile to taxanes for the treatment of triple-negative breast cancer and is synergistic with anti-CTLA-4 immunotherapy. Anti-Cancer Drugs. 2018.
Nanda R, et al. Pembrolizumab in patients with advanced triple-negative breast Cancer: phase Ib KEYNOTE-012 study. J Clin Oncol. 2016;34(21):2460–7.
Liu Z, et al. A comprehensive immunologic portrait of triple-negative breast Cancer. Transl Oncol. 2018;11(2):311–29.
Safonov A, et al. Immune gene expression is associated with genomic aberrations in breast Cancer. Cancer Res. 2017;77(12):3317–24.
Snyder A, et al. Genetic basis for clinical response to CTLA-4 blockade in melanoma. N Engl J Med. 2014;371(23):2189–99.
Rizvi NA, et al. Cancer immunology. Mutational landscape determines sensitivity to PD-1 blockade in non-small cell lung cancer. Science. 2015;348(6230):124–8.
Hugo W, et al. Genomic and transcriptomic features of response to anti-PD-1 therapy in metastatic melanoma. Cell. 2016;165(1):35–44.
Van Allen EM, et al. Genomic correlates of response to CTLA-4 blockade in metastatic melanoma. Science. 2015;350(6257):207–11.
Le DT, et al. PD-1 blockade in tumors with mismatch-repair deficiency. N Engl J Med. 2015;372(26):2509–20.
Barbie DA, et al. Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1. Nature. 2009;462(7269):108–12.
Hanzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics. 2013;14:7.
Yoshihara K, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612.
Subramanian A, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50.
Mootha VK, et al. PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003;34(3):267–73.
Yu G, et al. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.
Kanehisa M, et al. KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017;45(D1):D353–61.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.
Witten IH, et al. Data Mining: Practical Machine Learning Tools and Techniques. 4th ed: Morgan Kaufmann, San Francisco; 2016.
Newman AM, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7.
Carter SL, et al. Absolute quantification of somatic DNA alterations in human cancer. Nat Biotechnol. 2012;30(5):413–21.
Mermel CH, et al. GISTIC2.0 facilitates sensitive and confident localization of the targets of focal somatic copy-number alteration in human cancers. Genome Biol. 2011;12(4):R41.
Curtis C, et al. The genomic and transcriptomic architecture of 2,000 breast tumours reveals novel subgroups. Nature. 2012;486(7403):346–52.
Cancer Genome Atlas, N. Comprehensive molecular portraits of human breast tumours. Nature. 2012;490(7418):61–70.
Chung W, et al. Single-cell RNA-seq enables comprehensive tumour and immune cell profiling in primary breast cancer. Nat Commun. 2017;8:15081.
Jezequel P, et al. Gene-expression molecular subtyping of triple-negative breast cancer tumours: importance of immune response. Breast Cancer Res. 2015;17:43.
Bindea G, et al. Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer. Immunity. 2013;39(4):782–95.
Patel SP, Kurzrock R. PD-L1 expression as a predictive biomarker in Cancer immunotherapy. Mol Cancer Ther. 2015;14(4):847–56.
Chen X, et al. TNBCtype: a subtyping tool for triple-negative breast Cancer. Cancer Inform. 2012;11:147–56.
Goel S, et al. CDK4/6 inhibition triggers anti-tumour immunity. Nature. 2017;548(7668):471–5.
Jiang Z, et al. Immunogenomics analysis reveals that TP53 mutations inhibit tumor immunity in gastric Cancer. Transl Oncol. 2018;11(5):1171–87.
Pick R, et al. Coronin 1A, a novel player in integrin biology, controls neutrophil trafficking in innate immunity. Blood. 2017;130(7):847–58.
Pieters J, Muller P, Jayachandran R. On guard: coronin proteins in innate and adaptive immunity. Nat Rev Immunol. 2013;13(7):510–8.
Zhang S, et al. CD4+ T-cell-mediated anti-tumor immunity can be uncoupled from autoimmunity via the STAT4/STAT6 signaling axis. Eur J Immunol. 2009;39(5):1252–9.
Avram D, Califano D. The multifaceted roles of Bcl11b in thymic and peripheral T cells: impact on immune diseases. J Immunol. 2014;193(5):2059–65.
Knox JJ, et al. Characterization of T-bet and eomes in peripheral human immune cells. Front Immunol. 2014;5:217.
Carotta S, et al. The transcription factors IRF8 and PU.1 negatively regulate plasma cell differentiation. J Exp Med. 2014;211(11):2169–81.
Moroishi, T., et al., The hippo pathway kinases LATS1/2 suppress Cancer immunity. Cell, 2016. 167(6): p. 1525–1539 e17.
Guo X, et al. Single tumor-initiating cells evade immune clearance by recruiting type II macrophages. Genes Dev. 2017;31(3):247–59.
Janse van Rensburg HJ, et al. The hippo pathway component TAZ promotes immune evasion in human Cancer through PD-L1. Cancer Res. 2018;78(6):1457–70.
Ali HR, et al. Association between CD8+ T-cell infiltration and breast cancer survival in 12,439 patients. Ann Oncol. 2014;25(8):1536–43.
Garnelo M, et al. Interaction between tumour-infiltrating B cells and T cells controls the progression of hepatocellular carcinoma. Gut. 2017;66(2):342–51.
Thorsson V, et al. The immune landscape of Cancer. Immunity. 2018;48(4):812–830 e14.
Davoli T, et al. Tumor aneuploidy correlates with markers of immune evasion and with reduced response to immunotherapy. Science. 2017;355(6322).
Wang X, et al. Immunological therapy: A novel thriving area for triple-negative breast cancer treatment. Cancer Lett. 2018;442:409–428.
Aguiar PN Jr, et al. PD-L1 expression as a predictive biomarker in advanced non-small-cell lung cancer: updated survival data. Immunotherapy. 2017;9(6):499–506.
We thank Dr. Agnese Mariotti for editing the manuscript.
This study was funded by the China Pharmaceutical University (grant number 2632018YX01, 3150120001).
Availability of data and materials
We used four publicly available TNBC genomic datasets: METABRIC , TCGA , GSE75688  and GSE103091 . The METABRIC dataset was obtained from cBioPortal (http://www.cbioportal.org/study?id=brca_metabric#summary). The TCGA dataset was obtained from the TCGA data portal (https://portal.gdc.cancer.gov/). The GSE75688 and GSE103091 datasets were downloaded from the NCBI gene expression omnibus (https://www.ncbi.nlm.nih.gov/geo/). We obtained 29 immune signatures (represented by 29 different gene sets, respectively) from the publications [12, 35] (Additional file 1: Table S1).
Ethics approval and consent to participate
Ethical approval and consent to participate was waived since we used only publicly available data and materials in this study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. The 29 immune signatures represented by 29 different gene sets. (XLSX 41 kb)
Figure S1. Comparisons of the stromal content and tumor purity between TNBC subtypes (Mann–Whitney U test). (PDF 368 kb)
Figure S2. Comparisons of the expression levels of immune-related genes between TNBC subtypes. A. Comparisons of the expression levels of HLA genes between TNBC subtypes. B. Comparisons of the expression levels of immune cell subpopulation marker genes between TNBC subtypes. ANOVA test. *P < 0.05, **P < 0.01, ***P < 0.001. (PDF 168 kb)
Figure S3. Identification of TNBC subtype-specific pathways and gene ontology. A. KEGG pathways enriched in Immunity_H and Immunity_L. B. Gene modules significantly differentiating TNBC by subtype, survival time, or survival status. (PDF 164 kb)
Figure S4. Stem cell-associated activity is higher in Immunity_H than in the other subtypes. (PDF 110 kb)