Abstract
Purpose:
Breast cancer is a heterogeneous disease, and although advances in molecular subtyping have been achieved in recent years, most subtyping strategies target individual genes independent of one another and primarily concentrate on proliferative markers. The contributions of biological processes and immune patterns have been neglected in breast cancer subtype stratification.
Methods:
We performed a gene set variation analysis to simplify the information on biological processes using hallmark terms and to decompose immune cell data using the immune cell gene terms on 985 breast invasive ductal/lobular carcinoma RNAseq samples in the TCGA database.
Results:
The samples were gathered into three clusters following implementation of the t-SNE and DBSCAN algorithms and were categorized as 'hallmark-tsne' subtypes. Here, we identified a high-risk luminal A dominant breast cancer subtype (C3) that displayed increased motility, cancer stem cell-like features, a higher expression of hormone/luminal-related genes, a lower expression of proliferation-related genes and immune dysfunction. With regard to immune dysfunction, we observed that the motility-increased C3 subtype exhibited high granulocyte colony stimulating factor (G-CSF) expression accompanied by neutrophil aggregation. Cancer cells that produce high levels of G-CSF can stimulate neutrophils to form neutrophil extracellular traps, which promote cancer cell migration. This finding sheds light on one potential explanation for why the C3 subtype correlates with poor prognosis.
Conclusions:
The hallmark-tsne subtypes confirmed again that even the luminal A subtype is heterogeneous and can be further subdivided. The biological processes and immune heterogeneity of breast cancer must be understood to facilitate the improvement of clinical treatments.