Single-cell and bulk RNA sequencing analysis of B cell marker genes in TNBC TME landscape and immunotherapy

This study amied to investigate the prognostic characteristics of triple negative breast cancer (TNBC) patients by analyzing B cell marker genes based on single-cell and bulk RNA sequencing.

Utilizing single-cell sequencing data from TNBC patients, we examined tumor-associated B cell marker genes. Transcriptomic data from The Cancer Genome Atlas (TCGA) database were used as the foundation for predictive modeling. Independent validation set was conducted using the GSE58812 dataset. Immune cell infiltration into the tumor was assessed through various, including XCELL, TIMER, QUANTISEQ, CIBERSORT, CIBERSORT-ABS, and ssGSEA. The TIDE score was utilized to predict immunotherapy outcomes. Additional investigations were conducted on the immune checkpoint blockade gene, tumor mutational load, and the GSEA enrichment analysis.

Our analysis encompassed 22,106 cells and 20,556 genes in cancerous tissue samples from four TNBC patients, resulting in the identification of 116 B cell marker genes. A B cell marker gene score (BCMG score) involving nine B cell marker genes (ZBP1, SEL1L3, CCND2, TNFRSF13C, HSPA6, PLPP5, CXCR4, GZMB, and CCDC50) was developed using TCGA transcriptomic data, revealing statistically significant differences in survival analysis (P<0.05). Functional analysis demonstrated that marker genes were predominantly associated with immune-related pathways. Notably, substantial differences between the higher and lower- BCMG score groups were observed in terms of immune cell infiltration, immune cell activity, tumor mutational burden, TIDE score, and the expression of immune checkpoint blockade genes.

This study has established a robust model based on B-cell marker genes in TNBC, which holds significant potential for predicting prognosis and response to immunotherapy in TNBC patients.