Identification of Breast Cancer Immune-related Prognostic Characteristics in Tumor Microenvironment

Accumulated evidence suggest that tumor microenvironment (TME) plays a crucial role in breast cancer (BRCA) progression and therapeutic effects.

This study aimed to characterize immune-related BRCA subtypes in TME, and identify genes with prognostic value.

RNA sequencing profiles with corresponding clinical data from The Cancer Genome Atlas (TCGA) database of BRCA patients were downloaded to evaluate immune infiltration using the single-sample gene set enrichment (ssGAEA) algorithm. Further, BRCA was clustered according to immune infiltration status by consensus clustering analysis. Using Venn analysis, differentially expressed genes (DEGs) were overlapped to obtain candidate genes. Kaplan-Meier (K-M) analysis was performed to identify prognostic genes, and the results were verified in the GEO and METABRIC datasets. RT-qPCR was conducted to detect the mRNA expression of prognostic genes.

In the TCGA database, 3 immune-related BRCA subtypes were identified [cluster1 (C1), cluster2 (C2), and cluster3 (C2)]. The C2 subtype had better overall survival (OS) compared to the C1 subtype. Higher levels of immune markers and checkpoint protein were found in the C2 subtype than in others. By combining DEGs between BRCA and normal tissues, with the C1 and C2 subtypes associated with different OS, 25 BRCA candidate genes were identified. Among these, 8 genes were identified as prognostic genes for BRCA. RT-qPCR showed that the expressions of 2 genes were significantly elevated in BRCA tissues, while that of other genes were decreased.

Three BRCA subtypes were identified with the immune index, which may help design advanced treatment of BRCA. The data code used for the analysis in this article was available on GitHub (https://github.com/tangzhn/BRCA1.git).