Nuclear receptors exert profound effects on mammary gland physiology and have complex roles in the etiology of breast cancer. In addition to receptors for classic steroid hormones such as estrogen and progesterone, the nuclear vitamin D receptor (VDR) interacts with its ligand 1α,25(OH)2D3 to modulate the normal mammary epithelial cell genome and subsequent phenotype. Observational studies suggest that vitamin D deficiency is common in breast cancer patients and that low vitamin D status enhances the risk for disease development or progression. Genomic profiling has characterized many 1α,25(OH)2D3 responsive targets in normal mammary cells and in breast cancers, providing insight into the molecular actions of 1α,25(OH)2D3 and the VDR in regulation of cell cycle, apoptosis, and differentiation. New areas of emphasis include regulation of tumor metabolism and innate immune responses. However, the role of VDR in individual cell types (i.e., epithelial, adipose, fibroblast, endothelial, immune) of normal and tumor tissues remains to be clarified. Furthermore, the mechanisms by which VDR integrates signaling between diverse cell types and controls soluble signals and paracrine pathways in the tissue/tumor microenvironment remain to be defined. Model systems of carcinogenesis have provided evidence that both VDR expression and 1α,25(OH)2D3 actions change with transformation but clinical data regarding vitamin D responsiveness of established tumors is limited and inconclusive. Because breast cancer is heterogeneous, analysis of VDR actions in specific molecular subtypes of the disease may help to clarify the conflicting data. The expanded use of genomic, proteomic and metabolomic approaches on a diverse array of in vitro and in vivo model systems is clearly warranted to comprehensively understand the network of vitamin D regulated pathways in the context of breast cancer.

To elucidate potential mediators of vitamin D receptor (VDR) action in breast cancer, we profiled the genomic effects of its ligand 1,25-dihydroxyvitamin D3 (1,25D) in cells derived from normal mammary tissue and breast cancer. In non-transformed hTERT-HME cells, 483 1,25D responsive entities in 42 pathways were identified, whereas in MCF7 breast cancer cells, 249 1,25D responsive entities in 31 pathways were identified. Only 21 annotated genes were commonly altered by 1,25D in both MCF7 and hTERT-HME cells. Gene set enrichment analysis highlighted eight pathways (including senescence/autophagy, TGFβ signaling, endochondral ossification, and adipogenesis) commonly altered by 1,25D in hTERT-HME and MCF7 cells. Regulation of a subset of immune (CD14, IL1RL1, MALL, CAMP, SEMA6D, TREM1, CSF1, IL33, TLR4) and metabolic (ITGB3, SLC1A1, G6PD, GLUL, HIF1A, KDR, BIRC3) genes by 1,25D was confirmed in hTERT-HME cells and similar changes were observed in another comparable non-transformed mammary cell line (HME cells). The effects of 1,25D on these genes were retained in HME cells expressing SV40 large T antigen but were selectively abrogated in HME cells expressing SV40 + RAS and in MCF7 cells. Integration of the datasets from hTERT-HME and MCF7 cells with publically available RNA-SEQ data from 1,25D treated SKBR3 breast cancer cells enabled identification of an 11-gene signature representative of 1,25D exposure in all three breast-derived cell lines. Four of these 11 genes (CYP24A1, CLMN, EFTUD1, and SERPINB1) were also identified as 1,25D responsive in human breast tumor explants, suggesting that this gene signature may prove useful as a biomarker of vitamin D exposure in breast tissue.

Previous genomic profiling of immortalized, non-tumorigenic human breast epithelial cells identified a set of 1,25-dihydroxyvitamin D3 (1,25D) regulated genes with potential relevance to breast cancer prevention. In this report, we characterized the effect of 1,25D on a subset of these genes in six cell lines derived from mammary tissue and breast cancers. Non-tumorigenic cell lines included hTERT-HME1, HME and MCF10A cells which are often used to model normal breast epithelial cells. Breast cancer cell lines included MCF7 cells (a model of early stage, estrogen-dependent disease), DCIS.com cells (a derivative of MCF10A cells that models in situ breast cancer) and Hs578T cells (a model of metastatic disease). All of these cell lines express the vitamin D receptor (VDR) and exhibit anti-cancer responses to 1,25D such as changes in proliferation, apoptosis, metabolism, or invasion. Our comparative data demonstrate highly variable responses to 1,25D (100nM, 24h) between the cell lines. In both hTERT-HME1 and HME cell lines, CYP24A1, SLC1A1 and ITGB3 were up-regulated whereas KDR, GLUL and BIRC3 were down-regulated in response to 1,25D. In contrast, no changes in SLC1A1, ITGB3 or GLUL expression were detected in 1,25D treated MCF10A cells although KDR and BIRC3 were down-regulated by 1,25D. The effects of 1,25D on these genes in the breast cancer cell lines were blunted, with the DCIS.com cells exhibiting the most similar responses to the immortalized hTERT-HME1 and HME cells. The differences in cellular responses were not due to general impairment in VDR function as robust CYP24A1 induction was observed in all cell lines. Thus, our data indicate that the genomic changes induced by 1,25D are highly cell-type specific even in model cell lines derived from the same tissue. The implication of these findings is that genomic responses to changes in vitamin D status in vivo are likely to be distinct from individual to individual, particularly in neoplastic tissue.

Genomic profiling has identified a subset of metabolic genes that are altered by 1,25-dihydroxyvitamin D (1,25D) in breast cells, including GLUL, the gene that encodes glutamine synthetase (GS). In this study, we explored the relevance of vitamin D modulation of GLUL and other metabolic genes in the context of glutamine utilization and dependence. We show that exposure of breast epithelial cells to glutamine deprivation or a GS inhibitor reduced growth and these effects were exacerbated by cotreatment with 1,25D. 1,25D downregulation of GLUL was sufficient to reduce abundance and activity of GS. Flow cytometry demonstrated that glutamine deprivation induced S phase arrest, likely due to reduced availability of glutamine for DNA synthesis. In contrast, 1,25D induced G0/G1 arrest, indicating that its effects are not solely due to reduced glutamine synthesis. Indeed, 1,25D also reduced expression of GLS1 and GLS2 genes, which code for glutaminases that shunt glutamine into the tricarboxylic acid (TCA) cycle. Consistent with reduced entry of glutamine into the TCA cycle, 1,25D inhibited glutamine oxidation and the metabolic response to exogenous glutamine as analyzed by Seahorse Bioscience extracellular flux assays. Effects of 1,25D on GLUL/GS expression and glutamine oxidation were retained in human mammary epithelial (HME) cells that express SV-40 (HME-LT cells) but not in those that express SV-40 and oncogenic H-Ras (HME-PR cells). Furthermore, HME-PR cells exhibited glutamine independence and expressed constitutively high levels of GLUL/GS, which were unaffected by 1,25D. Collectively, these data suggest that 1,25D alters glutamine availability, dependence, and metabolism in nontransformed and preneoplastic mammary epithelial cells in association with cell cycle arrest.