Identification of a 9-gene prognostic signature for breast cancer

Txnrd1 as a prognosticator for recurrence, metastasis and response to neoadjuvant chemotherapy and radiotherapy in breast cancer patients

Thioredoxin reductase 1 (Txnrd1) is known to have prognostic significance in a subset of breast cancer patients. Despite the pivotal role of Txnrd1 in regulating several cellular and physiological processes in cancer progression and metastasis, its clinical significance is largely unrecognized. Here, we undertook a retrospective comprehensive meta-analysis of 13,322 breast cancer patients from 43 independent cohorts to assess prognostic and predictive roles of Txnrd1. We observed that Txnrd1 has a positive correlation with tumor grade and size and it is over-expressed in higher-grade and larger tumors. Further, hormone receptor-negative and HER2-positive tumors exhibit elevated Txnrd1 gene expression. Patients with elevated Txnrd1 expression exhibit significant hazards for shorter disease-specific and overall survival. While Txnrd1 has a positive correlation with tumor recurrence and metastasis, it has a negative correlation with time to recurrence and metastasis. Txnrd1High patients exhibit 2.5 years early recurrence and 1.3 years early metastasis as compared to Txnrd1Low cohort. Interestingly, patients with high Txnrd1 gene expression exhibit a pathologic complete response (pCR) to neoadjuvant chemotherapy, but they experience early recurrence after radiotherapy. Txnrd1High MDA-MB-231 cells exhibit significant ROS generation and reduced viability after doxorubicin treatment compared to Txnrd1Low MCF7 cells. Corroborating with findings from meta-analysis, Txnrd1 depletion leads to decreased survival, enhanced sensitivity to radiation induced killing, poor scratch-wound healing, and reduced invasion potential in MDA-MB-231 cells. Thus, Txnrd1 appears to be a potential predictor of recurrence, metastasis and therapy response in breast cancer patients.

Prognostic implication of novel immune-related signature in breast cancer

Checkpoint inhibitor therapy has become increasingly important and has been endorsed as a treatment regimen in breast cancer. But benefits were limited to a small proportion of patients. We aimed to develop an improved signature on the basis of immune genes for detection of potential benefit from immunotherapy. Gene expression data of patients with breast cancer initially extracted from The Cancer Genome Atlas were analyzed. Ten genes were selected from the interaction of differentially expressed genes as well as immune-related genes to develop a survival signature. We compared the high-risk and low-risk groups by gene set enrichment analysis, immune infiltration, checkpoint molecule expression and immunophenoscore. Ten genes were extracted from interactions of differentially expressed and immune-related genes. The immune risk score was determined on the basis of the Cox regression coefficient of hub genes and validated with the GSE96058 dataset. Immune cell infiltrates, including CD8 + T cells, plasma cells, follicular helper T cells, CD4 + memory T cells, M1 macrophages, regulatory T cells and resting NK cells, were more highly infiltrated in the high-risk group as compared to the low-risk group. Checkpoint molecules, including CTLA-4, PD-L1, TIM-3, VISTA, ICOS, PD-1, and PD-L2, were expressed at markedly lower levels in the high-risk group as compared to the low-risk group. Immunophenoscores, as a surrogate of response to immune checkpoint therapy, was observed significant lower in the high-risk group. The 10-gene prognostic signature could identify patients' survival and was correlated with the biomarkers of immune checkpoint inhibitor therapy, which may guide precise therapeutic decisions in clinical practice.

miR‑302d‑3p regulates the viability, migration and apoptosis of breast cancer cells through regulating the TMBIM6‑mediated ERK signaling pathway

MicroRNAs (miRs/miRNAs) play important roles in the occurrence, metastasis and prognosis of multiple types of cancers. However, the specific role of miR-302d-3p and its underlying mechanism in breast cancer (BC) have not yet been reported. The present study aimed to identify the role of miR-302D-3p in BC and its potential mechanism using BC cell lines MCF7 and MDA-MB-231 and normal breast epithelial cell MCF-10A. Cancer and paracancerous tissue from patients with BC were also used. Reverse transcription-quantitative PCR was performed to detect the expression of miR-302d-3p and transmembrane Bax inhibitor motif containing 6 (TMBIM6). Dual-luciferase reporter assays verified the binding sites of miR-302d-3p and TMBIM6. Immunohistochemistry was used to measure the expression of TMBIM6. Cell transfection techniques were used to overexpress or interfere with miR-302d-3p and TMBIM6. A Cell Counting Kit-8 assay was performed to detect cell viability, and migration was measured using a wound healing assay. Apoptosis was detected by flow cytometry. The expression levels of apoptosis-related proteins and pathway-related proteins were detected by western blotting. The expression of miR-302d-3p in BC cell lines was found to be downregulated. It was also demonstrated that miR-302d-3p could inhibit cell viability and migration and promote apoptosis. The expression of TMBIM6 in BC cell lines and tissues was upregulated. Upregulated miR-302d-3p was shown to inhibit viability and migration, and promote apoptosis by targeting TMBIM6, during which extracellular signal-regulated kinase (ERK) and its phosphorylation were inhibited in the ERK signaling pathway in cells. Overall, the present study demonstrated that miR-302d-3p could regulate the viability, migration and apoptosis of BC cells through regulating TMBIM6-mediated ERK signaling pathway.

Developing ZNF Gene Signatures Predicting Radiosensitivity of Patients with Breast Cancer

Adjuvant radiotherapy is one of the main treatment methods for breast cancer, but its clinical benefit depends largely on the characteristics of the patient. This study aimed to explore the relationship between the expression of zinc finger (ZNF) gene family proteins and the radiosensitivity of breast cancer patients. Clinical and gene expression data on a total of 976 breast cancer samples were obtained from The Cancer Genome Atlas (TCGA) database. ZNF gene expression was dichotomized into groups with a higher or lower level than the median level of expression. Univariate and multivariate Cox regression analyses were used to evaluate the relationship between ZNF gene expression levels and radiosensitivity. The Molecular Taxonomy Data of the International Federation of Breast Cancer (METABRIC) database was used for validation. The results revealed that 4 ZNF genes were possible radiosensitivity markers. High expression of ZNF644 and low expression levels of the other 3 genes (ZNF341, ZNF541, and ZNF653) were related to the radiosensitivity of breast cancer. Hierarchical cluster, Cox, and CoxBoost analysis based on these 4 ZNF genes indicated that patients with a favorable 4-gene signature had better overall survival on radiotherapy. Thus, this 4-gene signature may have value for selecting those patients most likely to benefit from radiotherapy. ZNF gene clusters could act as radiosensitivity signatures for breast cancer patients and may be involved in determining the radiosensitivity of cancer.

A Novel Risk Model Based on Lipid Metabolism-Associated Genes Predicts Prognosis and Indicates Immune Microenvironment in Breast Cancer

Background

Breast cancer (BRCA) is the most common tumor in women, and lipid metabolism involvement has been demonstrated in its tumorigenesis and development. However, the role of lipid metabolism-associated genes (LMAGs) in the immune microenvironment and prognosis of BRCA remains unclear.

Methods

A total of 1076 patients with BRCA were extracted from The Cancer Genome Atlas database and randomly assigned to the training cohort (n = 760) or validation cohort (n = 316). Kaplan–Meier analysis was used to assess differences in survival. Consensus clustering was performed to categorize the patients with BRCA into subtypes. Using multivariate Cox regression analysis, an LMAG-based prognostic risk model was constructed from the training cohort and validated using the validation cohort. The immune microenvironment was evaluated using the ESTIMATE and tumor immune estimation resource algorithms, CIBERSORT, and single sample gene set enrichment analyses.

Results

Consensus clustering classified the patients with BRCA into two subgroups with significantly different overall survival rates and immune microenvironments. Better prognosis was associated with high immune infiltration. The prognostic risk model, based on four LMAGs (MED10, PLA2G2D, CYP4F11, and GPS2), successfully stratified the patients into high- and low-risk groups in both the training and validation sets. High risk scores predicted poor prognosis and indicated low immune status. Subgroup analysis suggested that the risk model was an independent predictor of prognosis in BRCA.

Conclusion

This study demonstrated, for the first time, that LMAG expression plays a crucial role in BRCA. The LMAG-based risk model successfully predicted the prognosis and indicated the immune microenvironment of patients with BRCA. Our study may provide inspiration for further research on BRCA pathomechanisms.

Identification of a 9-gene prognostic signature for breast cancer

Breast cancer (BRCA) is the most common cancer among women and is the second leading cause of cancer death in women. In this study, we developed a 9‐gene prognostic signature to predict the prognosis of patients with BRCA. GSE20685, GSE42568, GSE20711, and GSE88770 were used as training sets. The Kaplan–Meier plot was constructed to assess survival differences and log‐rank test was performed to evaluate the statistical significance. The overall survival (OS) of patients in the low‐risk group was significantly higher than that in the high‐risk group. ROC analysis indicated that this 9‐gene signature shows good diagnostic efficiency both in OS and disease‐free survival (DFS). The 9‐gene signature was further validated through GSE16446, GSE7390, and TCGA‐BRCA datasets. We also established a nomogram that integrates clinicopathological features and 9‐gene signature. The analysis of the calibration plot showed that the nomogram has good prognostic performance. More convincingly, real‐time reverse transcription‐polymerase chain reaction (RT‐PCR) results indicated that the protective prognostic factors in BRCA patients were downregulated, whereas the dangerous prognostic factors were upregulated. The innovation of this article is not only constructing a prognostic gene signature, but also combining with clinical information to further establish a nomogram to better predict the survival probability of patients. It is worth mentioning that this signature also does not depend on other clinical factors or variables.