Identification of Breast Cancer Immune Subtypes by Analyzing Bulk Tumor and Single Cell Transcriptomes

Background: The histological and molecular classification of breast cancer (BC) is being used in the clinical management of this disease. However, subtyping of BC based on the tumor immune microenvironment (TIME) remains insufficiently explored, although such investigation may provide new insights into intratumor heterogeneity in BC and potential clinical implications for BC immunotherapy. Methods: Based on the enrichment scores of 28 immune cell types, we performed clustering analysis of transcriptomic data to identify immune-specific subtypes of BC using six different datasets, including five bulk tumor datasets and one single-cell dataset. We further analyzed the molecular and clinical features of these subtypes. Results: Consistently in the six datasets, we identified three BC subtypes: BC-ImH, BC-ImM, and BC-ImL, which had high, medium, and low immune signature scores, respectively. BC-ImH displayed a significantly better survival prognosis than BC-ImL. Triple-negative BC (TNBC) and human epidermal growth factor receptor-2-positive (HER2+) BC were likely to have the highest proportion in BC-ImH and the lowest proportion in BC-ImL. In contrast, hormone receptor-positive (HR+) BC had the highest proportion in BC-ImL and the lowest proportion in BC-ImH. Furthermore, BC-ImH had the highest tumor mutation burden (TMB) and predicted neoantigens, while BC-ImL had the highest somatic copy number alteration (SCNA) scores. It is consistent with that TMB and SCNA correlate positively and negatively with anti-tumor immune response, respectively. TP53 had the highest mutation rate in BC-ImH and the lowest mutation rate in BC-ImL, supporting that TP53 mutations promote anti-tumor immune response in BC. In contrast, PIK3CA displayed the highest mutation rate in BC-ImM, while GATA3 had the highest mutation rate in BC-ImL. Besides immune pathways, many oncogenic pathways were upregulated in BC-ImH, including ErbB, MAPK, VEGF, and Wnt signaling pathways; the activities of these pathways correlated positively with immune signature scores in BC. Conclusions: The tumors with the strong immune response (“hot” tumors) have better clinical outcomes than the tumors with the weak immune response (“cold” tumors) in BC. TNBC and HER2+ BC are more immunogenic, while HR + BC is less immunogenic. Certain HER2+ or HR + BC patients could be propitious to immunotherapy in addition to TNBC.