In silico screening using bulk and single-cell RNA-seq data identifies RIMS2 as a prognostic marker in basal-like breast cancer

Abstract Single-cell RNA-seq has become a powerful tool to understand tumor cell heterogenicity. This study tried to screen prognosis-related genes in basal-like breast tumors and evaluate their correlations with cellular states at the single-cell level. Bulk RNA-seq data of basal-like tumor cases from The Cancer Genome Atlas-Breast Cancer (TCGA-BRCA) and single-cell RNA-seq from GSE75688 were retrospectively reviewed. Kaplan–Meier survival curves, univariate and multivariate analysis based on Cox regression model were conducted for survival analysis. Gene set enrichment analysis (GSEA) and single-cell cellular functional state analysis were performed. Twenty thousand five hundred thirty genes with bulk RNA-seq data in TCGA were subjected to screening. Preliminary screening identified 10 candidate progression-related genes, including CDH19, AQP5, SDR16C5, NCAN, TTYH1, XAGE2, RIMS2, GZMB, LY6D, and FAM3B . By checking their profiles using single-cell RNA-seq data, only CDH19, SDR16C5, TTYH1, and RIMS2 had expression in primary triple-negative breast cancer (TNBC) cells. Prognostic analysis only confirmed that RIMS2 expression was an independent prognostic indicator of favorable progression free survival (PFS) (HR: 0.78, 95%: 0.64–0.95, P = .015). GSEA analysis showed that low RIMS2 group expression had genes significantly enriched in DNA Repair, and MYC Targets V2. Among the 89 basal-like cells, RIMS2 expression was negatively correlated with DNA repair and epithelial-to-mesenchymal transition (EMT). RIMS2 expression was negatively associated with DNA repair capability of basal-like breast tumor cells and might serve as an independent indicator of favorable PFS.