NRF1 motif sequence-enriched genes involved in ER/PR -ve HER2 +ve breast cancer signaling pathways.

Nuclear respiratory factor 1 (NRF1) transcription factor has recently been shown to control breast cancer progression. However, mechanistic aspects by which NRF1 may contribute to susceptibility to different breast tumor subtypes are still not fully understood. Since transcriptional control of NRF1 seems to be dependent on epidermal growth factor receptor signaling, herein, we investigated the role of NRF1 in estrogen receptor/progesterone receptor negative, but human epidermal growth factor receptor 2-positive (ER/PR -ve HER2 +ve) breast cancer. We found that both mRNA and protein levels of NRF1 and its transcriptional activity were significantly higher in ER/PR -ve HER2 +ve breast cancer samples compared to normal breast tissues. This was consistent with our observation of higher NRF1 protein expression in the experimental model of HER2+ breast cancer brain metastasis. To identify network-based pathways involved in the susceptibility to the ER/PR -ve HER2 +ve breast cancer subtype, the NRF1 transcriptional regulatory genome-wide landscape was analyzed using the approach consisting of a systematic integration of ChIP DNA-seq, RNA-Microarray, NRF1 protein-DNA motif binding, signal pathway analysis, and Bayesian machine learning. Our findings show that a high percentage of known HER2+ breast cancer susceptibility genes, including EGFR, IGFR, and E2F1, are under transcriptional control of NRF1. Promoters of several genes from the KEGG HER2+ breast cancer pathway and 11 signaling pathways linked to 6 hallmarks of cancer contain the NRF1 motif. By pathway analysis, key breast cancer hallmark genes of epithelial-mesenchymal transition, stemness, cell apoptosis, cell cycle regulation, chromosomal integrity, and DNA damage/repair were highly enriched with NRF1 motifs. In addition, we found using Bayesian network-based machine learning that 30 NRF1 motif-enriched genes including growth factor receptors-FGFR1, IGF1R; E2Fs transcription factor family-E2F1, E2F3; MAPK pathway-SHC2, GRB2, MAPK1; PI3K-AKT-mTOR signaling pathway-PIK3CD, PIK3R1, PIK3R3, RPS6KB2; WNT signaling pathway-WNT7B, DLV1, DLV2, GSK3B, NRF1, and DDB2, known for its role in DNA repair and involvement in early events associated with metastatic progression of breast cancer cells, were associated with HER2-amplified breast cancer. Machine learning search further revealed that the likelihood of HER2-positive breast cancer is almost 100% in a patient with the high NRF1 expression combined with expression patterns of high E2F3, GSK3B, and MAPK1, low or no change in E2F1 and FGFR1, and high or no change in PIK3R3. In summary, our findings suggest novel roles of NRF1 and its regulatory networks in susceptibility to the ER/PR -ve HER2 +ve aggressive breast cancer subtype. Clinical confirmation of our machine learned Bayesian networks will have significant impact on our understanding of the role of NRF1 as a valuable biomarker for breast cancer diagnosis and prognosis as well as provide strong rationale for future studies to develop NRF1 signaling-based therapeutics to target HER2+ breast cancer.