Transcriptomic profiling of Indian breast cancer patients revealed subtype-specific mRNA and lncRNA signatures

Abstract

Breast cancer (BC) is one of the leading causes of cancer-associated death in women. Despite the progress in therapeutic regimen, resistance and recurrence of breast cancer have affected the overall survival of patients. The present signatures, such as PAM50 and Oncotype DX, do not segregate the Indian breast samples based on molecular subtypes. This study aims at finding signatures of long noncoding RNA (lncRNA) and mRNA in Indian breast cancer patients using RNA-seq. We have analyzed the survival based on the menopausal and hormone status of 380 Indian breast cancer patients, and of these, we have sequenced and analyzed matched tumor–normal transcriptome of 17 (pre- and postmenopausal) Indian breast cancer patients representing six different subtypes, namely, four patients in triple-positive, three patients in estrogen receptor–positive (ER+ve), three patients in estrogen and progesterone receptors–positive (ER+ve, PR+ve), two patients in human epidermal growth factor receptor (Her2+ve), three patients in triple-negative, and one patient in ER+ve and Her2+ve subtypes. We have identified a 25 mRNA–27 lncRNA gene set, which segregated the subtypes in our data. A pathway analysis of the differentially expressed genes revealed downregulated ECM interaction and upregulated immune regulation, cell cycle, DNA damage response and repair, and telomere elongation in premenopausal women. Postmenopausal women showed downregulated metabolism, innate immune system, upregulated translation, sumoylation, and AKT2 activation. A Kaplan–Meier survival analysis revealed that menopausal status, grade of the tumor, and hormonal status displayed statistically significant effects (p < 0.05) on the risk of mortality due to breast cancer. Her2+ve patients showed low overall survival. One of the unique lncRNA-mRNA pairs specific to the EP-subtype, SNHG12 and EPB41, showed interaction, which correlates with their expression level; SNHG12 is downregulated and EPB41 is upregulated in EP samples.