Transcriptional effects of 1,25 dihydroxyvitamin D3physiological and supra-physiological concentrations in breast cancer organotypic culture
- Cintia Milani1,
- Maria Lucia Hirata Katayama1,
- Eduardo Carneiro de Lyra1, 2,
- JoEllen Welsh3,
- Laura Tojeiro Campos1,
- M Mitzi Brentani1,
- Maria do Socorro Maciel4,
- Rosimeire Aparecida Roela1,
- Paulo Roberto del Valle1,
- João Carlos Guedes Sampaio Góes2,
- Suely Nonogaki5,
- Rodrigo Esaki Tamura6 and
- Maria Aparecida Azevedo Koike Folgueira1Email author
© Milani et al.; licensee BioMed Central Ltd. 2013
Received: 11 May 2012
Accepted: 8 March 2013
Published: 15 March 2013
Vitamin D transcriptional effects were linked to tumor growth control, however, the hormone targets were determined in cell cultures exposed to supra physiological concentrations of 1,25(OH)2D3 (50-100nM). Our aim was to evaluate the transcriptional effects of 1,25(OH)2D3 in a more physiological model of breast cancer, consisting of fresh tumor slices exposed to 1,25(OH)2D3 at concentrations that can be attained in vivo.
Tumor samples from post-menopausal breast cancer patients were sliced and cultured for 24 hours with or without 1,25(OH)2D3 0.5nM or 100nM. Gene expression was analyzed by microarray (SAM paired analysis, FDR≤0.1) or RT-qPCR (p≤0.05, Friedman/Wilcoxon test). Expression of candidate genes was then evaluated in mammary epithelial/breast cancer lineages and cancer associated fibroblasts (CAFs), exposed or not to 1,25(OH)2D3 0.5nM, using RT-qPCR, western blot or immunocytochemistry.
1,25(OH)2D3 0.5nM or 100nM effects were evaluated in five tumor samples by microarray and seven and 136 genes, respectively, were up-regulated. There was an enrichment of genes containing transcription factor binding sites for the vitamin D receptor (VDR) in samples exposed to 1,25(OH)2D3 near physiological concentration. Genes up-modulated by both 1,25(OH)2D3 concentrations were CYP24A1, DPP4, CA2, EFTUD1, TKTL1, KCNK3. Expression of candidate genes was subsequently evaluated in another 16 samples by RT-qPCR and up-regulation of CYP24A1, DPP4 and CA2 by 1,25(OH)2D3 was confirmed. To evaluate whether the transcripitonal targets of 1,25(OH)2D3 0.5nM were restricted to the epithelial or stromal compartments, gene expression was examined in HB4A, C5.4, SKBR3, MDA-MB231, MCF-7 lineages and CAFs, using RT-qPCR. In epithelial cells, there was a clear induction of CYP24A1, CA2, CD14 and IL1RL1. In fibroblasts, in addition to CYP24A1 induction, there was a trend towards up-regulation of CA2, IL1RL1, and DPP4. A higher protein expression of CD14 in epithelial cells and CA2 and DPP4 in CAFs exposed to 1,25(OH)2D3 0.5nM was detected.
In breast cancer specimens a short period of 1,25(OH)2D3 exposure at near physiological concentration modestly activates the hormone transcriptional pathway. Induction of CYP24A1, CA2, DPP4, IL1RL1 expression appears to reflect 1,25(OH)2D3 effects in epithelial as well as stromal cells, however, induction of CD14 expression is likely restricted to the epithelial compartment.
KeywordsBreast cancer Calcitriol Gene expression Organotypic culture
Epidemiological data indicates higher incidence and mortality rates from breast cancer in low latitude regions. Among the mechanisms suggested for a relationship between sunlight and cancer is the genesis of vitamin D in the skin, resulting from the UV light action. In accordance with this hypothesis, there is evidence that lower 25(OH)D3[1–5] and 1,25(OH)2D3[6, 7] serum concentrations are encountered in patients with breast cancer, as compared with women without cancer, as well as in patients with advanced or metastatic disease in comparison with those with early-stage disease [8, 9]. In addition, 25(OH)D3 deficiency at diagnosis was related with poor prognosis, evaluated as metastasis-free and overall survival .
In human breast xenografts established in immunossupressed mice 1,25(OH)2D3 exerts growth inhibitory effects, and in mouse mammary organ culture exposed to chemical carcinogens, both 25(OH)D3 and 1,25(OH)2D3 mediate preventive effects [11–13]. However, the chemopreventive effect of vitamin D is still controversial, as supplementation trials on vitamin D3 and colon or breast cancer incidence have been inconsistent [14, 15]. One critical issue is that the appropriate supplementation dose for cancer prevention trials was not well established . On the other hand, clinical studies point to a clinical benefit for 1,25(OH)2D3 (or analogues) alone or in combination with chemotherapy in the treatment of hormone refractory prostate cancer and breast cancer skin lesions [17, 18]. However, concerns about hypercalcemic side effects limit the dose of 1,25(OH)2D3 (or analogues) that can be safely administered in vivo.
Phase I clinical studies indicate that subcutaneous doses of calcitriol given every other day result in peak 1,25(OH)2D3 serum concentration of 0.25-0.75 nM  while weekly pulses of oral calcitriol allow higher dose administration and peak serum concentrations of 1–15 nM . Although these vitamin D concentrations represent about 1.3-83 times the upper limit of physiologic serum levels, they are well below the concentrations (10-100nM) typically used to investigate hormone actions in cell culture studies. At these concentrations, 1,25(OH)2D3 exerts antiproliferative and pro apoptotic effects  and modulates angiogenesis [22, 23], invasion and metastasis [24, 25]. Among the downstream targets of the hormone are cyclin dependent kinase inhibitors as p21WAF1/CIP1 and p27KIP1; growth factors, receptors and associated proteins as TGFβ, TGFβ receptors and insulin-like growth factor binding protein-3 (IGFBP-3) [26–31]. In addition, gene expression profiling of breast cancer cell lines MCF7 and MDA-MB-231 have identified many potential 1,25(OH)2D3 target genes,  but again, these studies were conducted with supra physiological concentrations of calcitriol (50-100nM). Furthermore, experiments in cell lines do not reflect the complex array of interactions among malignant and stromal cells, secreted factors and extracellular matrix proteins taking place in the tumor microenvironment, which also modulate the hormone actions.
Although the majority of human breast cancers express vitamin D receptors (VDR) [7, 32, 33], there have been no demonstrations that 1,25(OH)2D3 modulates gene expression in human breast cancer samples. To address this research gap, a physiologically relevant in vitro model to study 1,25(OH)2D3 actions, represented by short term culture of fresh breast cancer tissue slices, which maintain the epithelial mesenchymal relationship and preserve tissue morphology and proliferation rate, was established [25, 34, 35]. With this organotypic culture system the transcriptional effects of 1,25(OH)2D3 at 0.5nM, a concentration that can be safely attained in vivo, and 100nM, the concentration typically used in cell culture studies, was compared. In addition, mammary cell lines and fibroblasts obtained from breast cancer samples were used to validate transcriptional targets of 1,25(OH)2D3 in epithelial and stromal cell types. Cancer associated fibroblasts (CAF) are interactive cells that infiltrate tumor specimens, influencing their behavior [36–38], which are also potential targets of the hormone. Although VDRs have been detected in fibroblasts obtained from prostate and breast tumors, few studies have compared 1,25(OH)2D3 mediated genomic effects in epithelial and stromal cells [39, 40]. The present study indicates that physiologically relevant concentrations of 1,25(OH)2D3 may influence gene expression in breast tumor slices cultured ex vivo, and that regulation of target genes likely occurs in both epithelial and stromal compartments of the tumor.
Characteristics of patients
Training group (n=5)
Validation group (n=16)
Tissue slice preparation and treatment
Tumor fragments were obtained immediately after tumor resection by the pathologist, who selected an involved area for this study. Fragments were placed into culture medium (RPMI 1640 with antibiotics and fungicide) and tissue slices were prepared using the Krumdieck tissue slicing system (Alabama Research and Development Corporation, Birmingham, AL, USA). Fragment thickness varied between 400–800 μm. Slices were cultured for 24 hours in 6-well plates (1 slice/well; 1–3 slices per treatment) containing 2 mL of culture media, RPMI supplemented with 10% v/v FBS, antibiotics and 0.001% ethanol (vehicle) or 1,25(OH)2D3 (Calbiochem, Darmstadt, Germany) 0.5nM or 100nM (from now on called physiological and supra-physiological concentrations, respectively). One slice of each sample was processed by FFPE and hematoxilin-eosin stained slides revealed that tumor samples contained > 50% malignant cells.
Fibroblasts primary culture
Primary fibroblast culture was established from tumor samples obtained from another five post-menopausal patients, diagnosed with invasive ductal carcinoma (histological grades II or III, three of them hormone receptor positive). Tumor samples were cut into small pieces and fibroblast primary culture was established through the explant methodology. After three cell passages, mesenchymal origin of the cells was confirmed by their spindle cell morphology and positive expression of vimentin [mouse anti human vimentin monoclonal antibody, clone Vim 3b4 (1:200); DAKO Corporation, Carpinteria, CA, USA] and alpha smooth muscle actin [(mouse monoclonal antibody anti human alpha smooth muscle actin, clone 1A4 (1:50); R&D Systems] and negative expression of cytokeratin [mouse monoclonal antibody anti human cytokeratin clone AE1/AE3 (1:100); DAKO] by immunocytochemistry (data not shown). Fibroblasts were then exposed to 1,25(OH)2D3 (Calbiochem, Darmstadt, Germany) 0.5nM or vehicle for 24 hours and after RNA extraction, RT-qPCR was performed to evaluate expression of candidate genes.
Culture of mammary epithelial cell lines
HB4A (normal mammary epithelial cell line) and C5.2a (HB4A transfected with HER2), both donated by Drs. Mike O’Hare and Alan Mackay, Ludwig Institute for Cancer Research, London, UK; SKBR3: breast cancer cell line overexpressing HER2; MDA MB-231: breast cancer cell line triple negative; and MCF-7: breast cancer cell line ER(+), acquired from American Type Culture Colection (Manassas, Virginia, USA), were cultured in RPMI-1640 supplemented with 10% fetal calf serum (FCS). After 24 hours, medium was replaced and 1,25(OH)2D3 0.5 nM (treated cells) or ethanol (control cells) was added. After 24 hs of treatment, total RNA was isolated using Trizol reagent and used in RT-qPCR.
RNA extraction and microarray hybridization
Tumor specimens were pulverized (Bio-Pulverizer™ BioSpec Products Inc., Oklahoma, USA) under liquid nitrogen and total RNA was isolated using RNeasy kit (Qiagen, Valencia, CA, USA), according to the manufacturer’s protocol. RNA integrity was verified in a Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA) and samples with RNA integrity number ≥ 6.6 were analyzed. Beginning with 100 ng total RNA, a two-round linear amplification was carried out, according to Affymetrix protocol (Two Cycle Target Labeling Kit, Affymetrix, Santa Clara, CA, USA). Afterwards, biotin-labeled cRNA was synthesized from double strand cDNA, using IVT labeling kit (Affymetrix) and 20 μg of biotinylated fragmented aRNA was hybridized onto Human Genome U133 Plus 2.0 GeneChip (Affymetrix
Hybridized arrays were scanned using Affymetrix GeneChip Scanner 3000 and after visual inspection, images were subjected to Affymetrix GeneChip Operating Software (GCOS) analysis to generate report files for quality control. Data normalization was performed using the Robust Multi-Array Average (RMA). Samples were categorized according to treatment in three groups: 1,25(OH)2D3 0.5nM, 1,25(OH)2D3 100nM and control. To establish a differential gene expression profile between vitamin D treated and untreated samples, SAM two class paired, provided on MEV (MultiExperiment Viewer – Boston, MA, USA) was used, after selecting 50% of the genes with the highest standard deviation. False discovery ratio (FDR) ≤0.10 was considered significant. In addition, results obtained with FDR≤0.01 are presented. Unsupervised hierarchical clustering based on Euclidean distance and average linkage was used to verify association patterns. The reliability of the clustering was assessed by the Bootstrap technique. Raw data complying with MIAME format was deposited at the Gene Expression Omnibus (GEO) data repository (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE27220) accession number GSE27220.
To explore functional enrichment associated with calcitriol treatment based on Ontologies (GO, Pathway), Regulome (TFBS, transcription factor binding site) Pharmacome (Drug-gene associations) among other features, differentially expressed genes were subject to subsequent analysis using ToppFun, available on ToppGene Suite (http://toppgene.cchmc.org/enrichment.jsp) and were considered significant if P < 0.05 .
Gene set enrichment analysis (GSEA) method was used to identify whether predefined gene sets might associate with gene expression differences between phenotypes. In this pairwise comparison, all genes are ranked based on signal-to-noise ratio and the alternative hypothesis that rank ordering of distinct pathway members is associated with a specific phenotype is tested . This methodology makes it possible to detect situations where all genes in a predefined set change in a small but coordinated way. FDR<0.10 was considered significant.
Real time RT-PCR
Reverse transcription was performed with random primers and Superscript III (Invitrogen Corporation, Carlsbad, CA, USA). Quantitative PCR (qPCR) was carried out using specific primers (Additional file 1: Table S1) and SYBR-green I (Sigma, St. Louis, MO, USA) in a Rotor-gene system (Corbett Research, Mortlake, Australia). Relative expression of target genes was calculated as 2-ΔΔCT, using GAPDH or ACTB as internal control (as indicated) and the average value of the target gene in control samples, as reference level.
Protein lysates from cell lines were made using RIPA buffer (1% NP-40, 0.1% SDS, 0.5% Sodium Deoxycholate in 1 × PBS) supplemented with complete mini protease inhibitor cocktail tablets (Roche; cat 04693124001). Afterwards, 50 μg of protein was subjected to SDS-PAGE and transferred to Hybond ECL membrane (GE Lifesciences), which was probed with the following primary antibodies overnight at 4°C: CD26 (DPP4, clone H-270, rabbit polyclonal antibody, 1:500; Santa Cruz Biotechonology Inc. Santa Cruz, CA, USA); CD14 (clone M305, sc9150, rabbit polyclonal antibody 1:500; Santa Cruz); β-actin (monoclonal Anti-β-Actin antibody produced in mouse, clone AC-15, ascites fluid, A5441, 1:2000, Sigma-Aldrich) and then with appropriate secondary antibodies (170–6515 Goat Anti-Rabbit Ig-G (H+L) HRP conjugate; 170–6516 Goat Anti-Mouse Ig-G (H+L) HRP conjugate; Bio-Rad.). Protein expression was detected with ECL Plus Western Blotting Detection Reagents (GE Lifesciences) in a ImageQuant LAS 4000 (GE Healthcare).
Fibroblasts were grown on coverslips in the absence or presence of 1,25(OH)2D3 0.5nM for 24 hours. Samples were fixed in 4% paraformaldehyde and permeabilized with 0.5% Triton X-100/PBS, in case of intracellular targets. Blocking of unspecific binding was performed with 2% BSA/PBS. Afterwards, cells were incubated with the primary antibody (CD26, clone H-270 rabbit polyclonal antibody, anti-DPP4, 1:200, Santa Cruz; CA II, clone G-2, sc-48351 mouse monoclonal antibody, 1:100, Santa Cruz) overnight in humid chamber at 4°C and then with the secondary antibody conjugated with Alexa Fluor 488 (1:700, species specific: goat anti-rabbit IgG, n° A11008; goat anti-mouse IgG n° A11001; Molecular Probes) for 1 h at room temperature in the dark. DAPI was added for nuclear staining. Images were acquired in a Olympus fluorescence microscope DX-5AI, using an Image Pro-PLUS 6,0 software.
Breast cancer slices from seven patients (six samples cultured in the absence (control) or presence of 1,25(OH)2D3 100nM and one sample cultured in the presence of 1,25(OH)2D3 0.5nM) were available for analysis (from the patients described in Table 1). Sections of 3 μm thickness were cut from paraffin blocks and antigen retrieval was carried out in 10 mM citrate buffer at pH 6.0 in humid heat under pressure cooker. Staining with the following specific antibody took place overnight at 4°C: CD14, clone M-305, (sc-9150Santa Cruz Biotechnology) rabbit polyclonal IgG, 1:800. Reaction was revealed with Novolink Polymer Detection Systems (Leica Biosystems, Newcastle, UK, cat: RE7280-k), followed by analysis in a Olympus fluorescence microscope DX-5AI (40x objective) and acquisition with an Image Pro-PLUS 6,0 software.
Detection of soluble CD14 (sCD14) in culture medium of tumor samples
Tumor slices from another four post-menopausal patients (median age 56 years) diagnosed with invasive ductal carcinoma clinical stages I-II, HER2 negative and hormone receptor positive (except for one tumor triple negative) were cultured with or without 1,25(OH)2D3 0.5nM or 100nM for 24 hours and 100 μL of the conditioned medium was used for soluble CD14 (sCD14) quantitative determination, through an enzyme-linked immunosorbent assay (Quantikine ELISA Human sCD14 Immunoassay, R&D Systems, Minneapolis, MN, USA). For every sample, two analyses on the same plate were carried out and the mean value was used.
Kolmogorov-Smirnov test was applied to check for normality of the data, followed by parametric or non-parametric tests, as appropriate. To detect an association between variables, Pearson chi-square or Fisher exact tests were used. A two-tailed p value ≤ 0.05 was considered significant. Analysis was undertaken using Instat (GraphPad Software, Inc., La Jolla, CA, USA) or SPSS (Chicago, IL, USA).
Twenty one post-menopausal patients with breast cancer clinical stages I-III were included in this study. Samples from five patients were analyzed in a training group, using microarray, and from another 16 patients were analyzed in a validation group, using RT-qPCR. There were no differences between groups concerning age, clinical stage, lymph node involvement; ductal histology; ER, PR and HER2 immunoexpression (Table 1).
Vitamin D transcriptional effects in breast cancer slices
Genes differentially modulated in breast tumor slices incubated with 0.5nM 1,25(OH) 2 D 3
cytochrome P450, family 24, subfamily A, polypeptide 1
dipeptidyl-peptidase 4 (CD26, adenosine deaminase complexing protein 2)
cytochrome P450, family 26, subfamily B, polypeptide 1
potassium channel, subfamily K, member 3
spindlin family, member 3
elongation factor Tu GTP binding domain containing 1
Fc fragment of IgG, low affinity IIc, receptor for (CD32)
SAM domain, SH3 domain and nuclear localization signals 1
Using GSEA (motif, transcription factors), to compare samples treated with 0.5nM 1,25(OH)2D3 treated and untreated, only one gene set was enriched at FDR ≤ 0.1, namely DR3, comprising genes containing a motif for vitamin D receptor (VDR) around the transcription start site (Additional file 2: Table S2).
Functional categories of genes up-regulated in breast tumor slices incubated in 100 nM 1,25(OH) 2 D 3
Genes up-regulated in VD3 100nM treated samples
Response to external stimulus
TGFBR2, LXN, THBD,PTEN, PTGER3, HBEGF, CYP24A1, ACVRL1, CCL19, FOXF1, FYN, OSM, CD28, CD14, BMP6, BMP2, PLAT, CD1D, PKD2, SERPINA1, PROCR, ALDH1A2
Response to wounding
TGFBR2, THBD, PTGER3, HBEGF, ACVRL1, CCL19, FOXF1, OSM, CD28, CD14, BMP6, BMP2, PLAT, SERPINA1, PROCR
Regulation of leukocyte mediated immunity
CD226, FOXF1, DPP4, CD28, CD1D, GIMAP1
Pathway-restricted SMAD protein phosphorylation
TGFBR2, ACVR1B, BMP6, BMP2
Urogenital system development
CRLF1, PTEN, CYP19A1, BMP6, BMP2, CA2, PKD2, ALDH1A2
Vitamin metabolic process
TKTL1, SLC22A4, CYP24A1, DHRS9, MTAP, PSAT1, ALDH1A2
Positive regulation of alpha-beta T cell activation
TGFBR2, CD28, CD1D, GIMAP1
TGF-beta signaling pathway
TGFBR2, ACVRL1, ACVR1B, IDA, BMP6, BMP2
EFTUD1, EHBP1, TRIM56, HBEGF, CYP24A1, CYP19A1, ACVRL1, ACVR1B, SEMAD6D, ARRDC4, CD14, BMP6, PRKD1, BMP2, CLMN, PLAT, IL1RL1, PRKCH, SLC1A1, CA2, FAM20C, SHE
To determine overlapping genes up-regulated by both calcitriol concentrations (at FDR ≤ 0.1), a Venn diagram was assembled. This approach identified five commonly up-modulated genes: CYP24A1, DPP4, EFTUD1, TKTL1 and KCNK3.
A subset of seven genes was selected for further analysis in samples from another group of patients, using qPCR. Candidates were chosen from microarray analysis and included two genes modulated by both calcitriol concentrations: CYP24A1 and DPP4; and five genes regulated by 100nM calcitriol at a fold change > 2, compared to control samples: IL1RL1, SHE, CD14, CA2 and BMP6. Initially, significant correlations between gene expression values obtained from the microarray dataset and those obtained by subsequent qPCR analysis in the first group of five patients were evaluated, as a technical validation procedure. In these 15 samples (control; 0.5nM 1,25(OH)2D3 and 100nM 1,25(OH)2D3) significant direct correlations were demonstrated for all genes, except for BMP6 (Additional file 5: Table S4).
Vitamin D transcriptional effects in epithelial and stromal cells
The effects of 1,25(OH)2D3 0.5nM on the expression of CYP24A1, DPP4, IL1RL1, CD14, CA2 and BMP6, were further explored in breast tumor derived cells, representing the epithelial and stromal compartments, using RT-qPCR. For this analysis, normal and cancerous breast cell lines (HB4A, C5.4, SKBR3, MDA-MB231, MCF-7) and cancer associated fibroblasts (primary cultures obtained from fresh tumor samples) were used. In the breast-derived epithelial cell lines, robust expression of CYP24A1 was observed in all lineages, indicating functional VDR expression. Breast cell lines (HB4a, C5.4, SKBR3) that exhibited low baseline CYP24A1 expression showed larger fold-induction of this gene than cell lines (MCF-7, MD-MBA-231) presenting high baseline CYP24A1. Expression of CA2, CD14 and IL1RL1, was significantly induced by 1,25(OH)2D3 0.5nM, but considerable variability in the response of individual lineages was observed, and cells displaying the most robust up-regulation of CYP24A1 in response to 1,25(OH)2D3 did not necessarily exhibit the highest induction of the other target genes. Three of the breast cancer cell lines demonstrated up-regulation of BMP6 in response to 1,25(OH)2D3 0.5nM however, the group response was not statistically significant.
Vitamin D effects on protein expression
Vitamin D effects in protein expression were analyzed in tumor slices and culture medium, as well as in epithelial cell lines and fibroblasts.
The primary goal of this work was to evaluate the transcriptional responses of breast cancer samples to physiologically relevant concentrations of 1,25(OH)2D3, using a culture model that retains features of intact tumors, such as stromal-epithelial interactions. Microarray analysis identified nine genes that were significantly altered within 24 h of exposure to 1,25(OH)2D3 0.5nM, a concentration that is physiologically achievable in patients. Of these, the vitamin D target gene CYP24A1(which codes a cytochrome P450 enzyme, that hydroxylates 25(OH)D3 and 1,25(OH)2D3 to less active forms 24,25(OH)2D3 and 1,24,25(OH)3D3) was induced over 7-fold in microarray analysis and was validated in another set of tumor samples, clearly indicating activation of VDR signaling. Additional evidence for activation of the VDR pathway in this dataset was obtained by GSEA, which indicated a trend towards the enrichment of genes sharing DR3 binding sites, a consensus motif for VDR.
Comparison of microarray data from tumor slices cultured with 0.5nM vs. 100nM 1,25(OH)2D3 indicated a clear concentration effect, as the number of differentially expressed transcripts increased from nine at 0.5nM to 186 at 100nM (20 fold increment). Induction of CYP24A1 increased from 7-fold (at 0.5 nM) to 70-fold (at 100nM) - a 10 fold enhancement. In both datasets, the majority of genes (approximately 75%) were up-regulated rather than down-regulated by 1,25(OH)2D3, consistent with other array data from established cell lines cultured with high dose 1,25(OH)2D3in vitro[43–45].
In addition to CYP24A1, five other genes were commonly up-regulated in tumor slices exposed to both low and high concentrations of 1,25(OH)2D3: DPP4, KCKN3, EFTUD1, TKTL1 and CA2. All, except TKTL1 (transketolase-like 1) have been previously identified as VDR target genes in various model systems. DPP4 (dipeptidyl-peptidase 4, also called CD26) was up-regulated in artery smooth muscle cells exposed to 1,25(OH)2D3 and its overexpression in distinct cell types (melanocytes, non-small cell lung, prostate and neuroblastoma cells) triggered anti-tumorigenic effects including cell growth arrest, inhibition of cell migration and increased apoptosis . KCNK3 (potassium channel, subfamily K, member) was induced by 1,25(OH)2D3 in artery smooth muscle cells, and EFTUD1 (elongation factor Tu GTP binding domain containing 1) in oral squamous carcinoma, breast cancer associated fibroblasts, immortalized prostate cells and lymphoblastoid cell lines [40, 43–46]. CA2 (carbonic anhydrase II) mRNA appeared to be directly induced by 1,25(OH)2D3 in myelomonocytic cell lines but indirectly regulated in osteoclast progenitors, where the physical communication with stromal cells seems to be required [48, 49]. CYP26B1 (cytochrome P450, family 26, subfamily b, polypeptide 1) which was up-regulated in samples treated with 1,25(OH)2D3 0.5nM, was previously identified as a vitamin D induced gene in immortalized non-transformed prostate epithelial and oral squamous carcinoma cell lines, and in silico analysis has tentatively identified a VDR binding site at this genomic region [43, 44].
Other authors have analyzed physiological concentration effects of vitamin D using animal models. Vitamin D supplemented diet as well as calcitriol injections were shown to stimulate the VDR pathway, mildly increasing CYP24A1 expression (x2) in MCF-7 xenografts in immunocompromised mice . Interestingly, vitamin D transcriptional effects may not overlap in tumor specimens and non-transformed mammary glands in the MMTV-neu transgenic mouse model of breast cancer, fed a high vitamin D diet . Comparison between cancer and normal cells is an interesting issue, as vitamin D potential effects in cancer prevention have also been claimed. In accordance with the previous work , differences in transcriptional targets were also described for breast cancer associated fibroblasts (CAF) and normal adjacent fibroblasts (NAF) exposed to 1,25(OH)2D3 in a supra-physiological concentration. Among up-regulated genes 45.7% were commonly modulated in CAFs and NAFs, however, 36.4% were exclusively up-regulated in NAFs and 17.4% exclusively up-regulated in CAFs . In addition, looking at overlapping genes in the Venn diagram of vitamin D up-regulated transcripts in six works [40, 43–46], only seven intersections were found in non-cancer cells: AKR1B1, CRIP1, FZD8, MREG (in immortalized prostate cells and NAF), BCAT1, GCLC (in coronary artery smooth muscle cells and NAFs) and PRR6 (in immortalized prostate cells and coronary artery smooth muscle cells). Furthermore, it was reported that vitamin D response is blunted in transformed HME normal mammary cells as compared with parental normal cells . The last works evaluating vitamin D effects in normal cells however, were performed using supra-physiological concentrations of 1,25(OH)2D3 (10-100nM) or analogs and the role of physiological concentrations of the hormone in normal cells is not fully established.
At 100nM, 1,25(OH)2D3 exerted more extensive transcriptional effects, and at least 40 of the induced genes in breast cancer organotypic culture have already been reported as up regulated by the hormone, such as ALCAM, ARRDC4, BMP2, BMP6, CA2, CD14, CLIC6, CILP, CLMN, CYP19A1, DCLDB1, EFTUD1, EHBP1, FAM20C, FOXF1, FRAS1, GOS2, GRK5, HBGEF, HSMPP8, IL1RL1, KCNK3, KIAA0500, PKD2, RGNEF, SEMA6D, SERPINB1, SLC1A1, THBD, TIMP1, TRIM56[40, 43–46]. However, co-aggregation of paired samples (treated and untreated) upon cluster analysis suggests that an individual dominant transcriptional profile was maintained, regardless of treatment. These results were not unexpected, as a high degree of transcriptional similarity was also demonstrated for matched pre and post-neoadjuvant chemotherapy, even though the chemotherapy exerts a more pronounced acute cellular effect than hormonal treatments [53–55].
Some of the genes induced by 100nM 1,25(OH)2D3 concentration are involved in TGF beta signaling pathway, in accordance with other authors [56, 57]. Other genes are involved in regulation of leukocyte mediated immunity and positive regulation of alpha-beta T cell activation, including CD14, which encodes a receptor to bacterial lipopolysaccharide, as previously reported in a variety of cells as mononuclear phagocytes, normal human epidermal keratinocytes, oral squamous carcinoma, immortalized non-transformed prostate epithelial cell lines and malignant breast cells [43, 56, 58].
The present tumor slice model represents a heterogeneous combination of epithelial and stromal cells, in which the complex array of reciprocal interactions taking place in the tumor microenvironment, including cell-cell contacts and a variety of secreted factors, might modulate the overall response to 1,25(OH)2D3. Hence, after evaluating the hormone effects in tumor slices, the effects of 1,25(OH)2D3 0.5nM in defined populations of cancer associated fibroblasts and epithelial cells were compared. This data indicated that even though CYP24A1 was induced in both fibroblasts and epithelial cells, CD14, CA2, and IL1RL1 were primarily induced in epithelial cells. There was also a trend towards up-regulation of CA2, DPP4 and IL1RL1 in cancer associated fibroblasts.
One major strengthen of this work was the comparison of achievable versus supra-physiological concentrations of 1,25(OH)2D3 in breast cancer slices, a model that preserves the epithelial-mesenchimal interactionss, indicating that effects are much less intense in near physiological concentrations. However, a weakness of this work was the small number of samples used in microarray experiments. These effects however, were later confirmed in a larger number of tumor samples and cell lines, using RT-PCR, even though they were more difficult to detect at the protein level, in face of the discrete changes induced by 0.5nM 1,25(OH)2D3.
Our main conclusion is that a very modest transcriptional response may be observed after exposure to 1,25(OH)2D3, within the physiological concentration range. Gene targets in breast cancer samples, including CYP24A1, DPP4 and CA2, seem to be shared by both fibroblasts and epithelial cells. A higher number of genes may be induced by a supra-physiological concentration of the hormone. Further studies employing physiological and supra-physiological concentrations may help to elucidate the hormone’s potential effects in breast cancer prevention and treatment, including calcitriol supplementation effects in post-menopausal women and calcitriol intra-tumoral effects in breast cancer xenografts.
The authors would like to acknowledge the assistance of Dr. Sridar Chittur for helpful discussions on microarray data analysis and Dr. Igor Moysés Longo Snitcovsky for critical review and important suggestions to the manuscript.
This work was supported by the Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP grants 07/04799-2; 09/10088-7; 08/51750-1) and Coordenação de Aperfeiçoamento de Pessoal de nível Superior (CAPES).
- Lowe LC, Guy M, Mansi JL, Peckitt C, Bliss J, Wilson RG, Colston KW: Plasma 25-hydroxy vitamin D concentrations, vitamin D receptor genotype and breast cancer risk in a UK Caucasian population. Eur J Cancer. 2005, 41 (8): 1164-1169. 10.1016/j.ejca.2005.01.017.View ArticlePubMed
- Bertone-Johnson ER, Chen WY, Holick MF, Hollis BW, Colditz GA, Willett WC, Hankinson SE: Plasma 25-hydroxyvitamin D and 1,25-dihydroxyvitamin D and risk of breast cancer. Cancer Epidemiol Biomarkers Prev. 2005, 14 (8): 1991-1997. 10.1158/1055-9965.EPI-04-0722.View ArticlePubMed
- Abbas S, Linseisen J, Slanger T, Kropp S, Mutschelknauss EJ, Flesch-Janys D, Chang-Claude J: Serum 25-hydroxyvitamin D and risk of post-menopausal breast cancer–results of a large case–control study. Carcinogenesis. 2008, 29 (1): 93-99.View ArticlePubMed
- Crew KD, Gammon MD, Steck SE, Hershman DL, Cremers S, Dworakowski E, Shane E, Terry MB, Desai M, Teitelbaum SL: Association between plasma 25-hydroxyvitamin D and breast cancer risk. Cancer Prev Res (Phila). 2009, 2 (6): 598-604. 10.1158/1940-6207.CAPR-08-0138.View Article
- Yao S, Sucheston LE, Millen AE, Johnson CS, Trump DL, Nesline MK, Davis W, Hong CC, McCann SE, Hwang H: Pretreatment serum concentrations of 25-hydroxyvitamin D and breast cancer prognostic characteristics: a case–control and a case-series study. PLoS One. 2011, 6 (2): e17251-10.1371/journal.pone.0017251.PubMed CentralView ArticlePubMed
- Janowsky EC, Lester GE, Weinberg CR, Millikan RC, Schildkraut JM, Garrett PA, Hulka BS: Association between low levels of 1,25-dihydroxyvitamin D and breast cancer risk. Public Health Nutr. 1999, 2 (1): 283-291.PubMed
- de Lyra EC, da Silva IA, Katayama ML, Brentani MM, Nonogaki S, Goes JC, Folgueira MA: 25(OH)D3 and 1,25(OH)2D3 serum concentration and breast tissue expression of 1alpha-hydroxylase, 24-hydroxylase and Vitamin D receptor in women with and without breast cancer. J Steroid Biochem Mol Biol. 2006, 100 (4–5): 184-192.View ArticlePubMed
- Mawer EB, Walls J, Howell A, Davies M, Ratcliffe WA, Bundred NJ: Serum 1,25-dihydroxyvitamin D May Be related inversely to disease activity in breast cancer patients with bone metastases. J Clin Endocrinol Metab. 1997, 82 (1): 118-122. 10.1210/jc.82.1.118.PubMed
- Palmieri C, MacGregor T, Girgis S, Vigushin D: Serum 25-hydroxyvitamin D levels in early and advanced breast cancer. J Clin Pathol. 2006, 59 (12): 1334-1336. 10.1136/jcp.2006.042747.PubMed CentralView ArticlePubMed
- Goodwin PJ, Ennis M, Pritchard KI, Koo J, Hood N: Prognostic effects of 25-hydroxyvitamin D levels in early breast cancer. J Clin Oncol. 2009, 27 (23): 3757-3763. 10.1200/JCO.2008.20.0725.View ArticlePubMed
- Eisman JA, Barkla DH, Tutton PJM: Suppression of in vivo growth of human cancer solid tumor xenografts by 1,25-dihydroxyvitamin D 3. Cancer Res. 1987, 47: 21-25.PubMed
- Pourgholami MH, Akhter J, Lu Y, Morris DL: In vitro and in vivo inhibition of liver cancer cells by 1,25-dihydroxyvitamin D3. Cancer Lett. 2000, 151: 97-102. 10.1016/S0304-3835(99)00416-4.View ArticlePubMed
- Peng X, Hawthorne M, Vaishnav A, St-Arnaud R, Mehta RG: 25-Hydroxyvitamin D3 is a natural chemopreventive agent against carcinogen induced precancerous lesions in mouse mammary gland organ culture. Breast Cancer Res Treat. 2009, 113 (1): 31-41. 10.1007/s10549-008-9900-0.PubMed CentralView ArticlePubMed
- Wactawski-Wende J, Kotchen JM, Anderson GL, Assaf AR, Brunner RL, O’Sullivan MJ, Margolis KL, Ockene JK, Phillips L, Pottern L: Calcium plus vitamin D supplementation and the risk of colorectal cancer. N Engl J Med. 2006, 354 (7): 684-696. 10.1056/NEJMoa055222.View ArticlePubMed
- Chlebowski RT, Johnson KC, Kooperberg C, Pettinger M, Wactawski-Wende J, Rohan T, Rossouw J, Lane D, O’Sullivan MJ, Yasmeen S: Calcium plus vitamin D supplementation and the risk of breast cancer. J Natl Cancer Inst. 2008, 100 (22): 1581-1591. 10.1093/jnci/djn360.PubMed CentralView ArticlePubMed
- Garland CF, French CB, Baggerly LL, Heaney R: Vitamin D supplement doses and serum 25-hydroxyvitamin D in the range associated with cancer prevention. Anticancer Res. 2011, 31: 607-611.PubMed
- Bower M, Colston KW, Stein RC, Hedley A, Gazet JC, Ford HT, Coombes RC: Topical calcipotriol treatment in advanced breast cancer. Lancet. 1991, 337: 701-702. 10.1016/0140-6736(91)90280-3.View ArticlePubMed
- Beer TM, Ryan CW, Venner PM, Petrylak DP, Chatta GS, Ruether JD, Redfern CH, Fehrenbacher L, Saleh MN, Waterhouse DM: Double-blinded randomized study of high-dose calcitriol plus docetaxel compared with placebo plus docetaxel in androgen-independent prostate cancer: a report from the ASCENT Investigators. J Clin Oncol. 2007, 25 (6): 669-674. 10.1200/JCO.2006.06.8197.View ArticlePubMed
- Smith DC, Johnson CS, Freeman CC, Muindi J, Wilson JW, Trump DL: A Phase I Trial of Calcitriol (1,25-Dihydroxycholecalciferol) in Patients with Advanced Malignancy. Clin Cancer Res. 1999, 5: 2431-2439.
- Beer TM, Munar M, Henner WD: A phase I trial of pulse calcitriol in patients with refractory malignancies: pulse dosing permits substantial dose escalation. Cancer. 2001, 91 (12): 2431-2439. 10.1002/1097-0142(20010615)91:12<2431::AID-CNCR1278>3.0.CO;2-3.View ArticlePubMed
- G-JCMvd B, Pols HAP, Leeuwen JPTM: Anti-Tumor Effects of 1,25-Dihydroxyvitamin D3 and Vitamin D Analogs. Curr Pharm Des. 2000, 6: 717-732. 10.2174/1381612003400498.View Article
- Nakagawa K, Sasaki Y, Kato S, Kubodera N, Okano T: 22-Oxa-1alpha,25-dihydroxyvitamin D3 inhibits metastasis and angiogenesis in lung cancer. Carcinogenesis. 2005, 26 (6): 1044-1054. 10.1093/carcin/bgi049.View ArticlePubMed
- Mantell DJ, Owens PE, Bundred NJ, Mawer EB, Canfield AE: 1α,25-Dihydroxyvitamin D3 Inhibits Angiogenesis In Vitro and In Vivo. Circ Res. 2000, 87: 214-220. 10.1161/01.RES.87.3.214.View ArticlePubMed
- Swami S, Raghavachari N, Muller UR, Bao YP, Feldman D: Vitamin D growth inhibition of breast cancer cells: gene expression patterns assessed by cDNAmicroarray. Breast Cancer Res Treat. 2003, 80: 49-62. 10.1023/A:1024487118457.View ArticlePubMed
- Barbosa EM, Nonogaki S, Katayama ML, Folgueira MA, Alves VF, Brentani MM: Vitamin D3 modulation of plasminogen activator inhibitor type-1 in human breast carcinomas under organ culture. Virchows Arch. 2004, 444 (2): 175-182. 10.1007/s00428-003-0929-5.View ArticlePubMed
- Wu G, Fan RS, Li W, Ko TC, Brattain MG: Modulation of cell Cycle control by vitamin D3 and its analogue, EB1089, in human breast cancer cells. Oncogene. 1997, 15: 155-1563.View Article
- Verlinden L, Verstuyf A, Convents R, Marcelis S, Camp MV, Bouillon R: Action of 1,25(OH)2D3 on the cell cycle genes, cyclin D1, p21 and p27 in MCF-7 cells. Mol Cell Endo. 1998, 142: 57-65. 10.1016/S0303-7207(98)00117-8.View Article
- Jensen SS, Mandsen MW, Lukas J, Binderup L, Bartek J: Inhibitory Effects of 1,25-Dihydroxyvitamin D3 on the G1–S Phase-Controlling Machinery. Mol Endocrinol. 2001, 15 (8): 1370-1380. 10.1210/me.15.8.1370.PubMed
- Koli K, Keski-Oja J: 1,25-Dihydroxyvitamin D3 enhances the expression of transforming growth factor β1 and its latent form binding protein in cultured breast carcinoma cells. Cancer Res. 1995, 55: 1540-1546.PubMed
- Wu G, Fan RS, Li W, Srinivas V, Brattain MG: Regulation of transforming growth factor-BetaType II receptor expression in human breast cancer MCF-7 cells by vitamin D3 and its analogues. J Biol Chem. 1998, 273 (13): 7749-7756. 10.1074/jbc.273.13.7749.View ArticlePubMed
- Katayama MLH, Pasini FS, Folgueira MAAK, Snitcovsky IML, Brentani MM: Molecular targets of 1,25(OH)2D3 in HC11 normal mouse mammary cell line. J Steroid Biochem Mol Biol. 2003, 84 (1): 57-69. 10.1016/S0960-0760(03)00004-9.View ArticlePubMed
- Bortman P, Folgueira MAAK, Katayama ML, Snitcovsky IML, Brentani MM: Antiproliferative effects of 1,25-dihydroxyvitamin D3 on breast cells - A mini review. Braz J Med Biol Res. 2002, 35 (1): 1-9.View ArticlePubMed
- Lopes N, Sousa B, Martins D, Gomes M, Vieira D, Veronese LA, Milanezi F, Paredes J, Costa JL, Schmitt F: Alterations in Vitamin D signalling and metabolic pathways in breast cancer progression: a study of VDR, CYP27B1 and CYP24A1 expression in benign and malignant breast lesions. BMC Cancer. 2010, 10: 483-10.1186/1471-2407-10-483.PubMed CentralView ArticlePubMed
- Sobral RA, Honda ST, Katayama ML, Brentani H, Brentani MM, Patrao DF, Folgueira MA: Tumor slices as a model to evaluate doxorubicin in vitro treatment and expression of trios of genes PRSS11, MTSS1, CLPTM1 and PRSS11, MTSS1, SMYD2 in canine mammary gland cancer. Acta Vet Scand. 2008, 50: 27-10.1186/1751-0147-50-27.PubMed CentralView ArticlePubMed
- Milani C, Welsh J, Katayama ML, Lyra EC, Maciel MS, Brentani MM, Folgueira MA: Human breast tumor slices: a model for identification of vitamin D regulated genes in the tumor microenvironment. J Steroid Biochem Mol Biol. 2010, 121 (1–2): 151-155.View ArticlePubMed
- Gache C, Berthois Y, Cvitkovic E, Martin P, Saez S: Differential regulation of normal and tumoral breast epithelial cell growth by fibroblasts and 1,25-dihydroxyvitaminD3. Breast Cancer Res Treat. 1999, 55: 29-39. 10.1023/A:1006163418479.View ArticlePubMed
- Rozenchan PB, Carraro DM, Brentani H, de Carvalho Mota LD, Bastos EP, e Ferreira EN, Torres CH, Katayama ML, Roela RA, Lyra EC: Reciprocal changes in gene expression profiles of cocultured breast epithelial cells and primary fibroblasts. Int J Cancer. 2009, 125 (12): 2767-2777. 10.1002/ijc.24646.View ArticlePubMed
- Santos RP, Benvenuti TT, Honda ST, Del Valle PR, Katayama ML, Brentani HP, Carraro DM, Rozenchan PB, Brentani MM, de Lyra EC: Influence of the interaction between nodal fibroblast and breast cancer cells on gene expression. Tumour Biol. 2011, 32 (1): 145-157. 10.1007/s13277-010-0108-7.PubMed CentralView ArticlePubMed
- Hidalgo AA, Montecinos VP, Paredes R, Godoy AS, McNerney EM, Tovar H, Pantoja D, Johnson C, Trump D, Onate SA: Biochemical characterization of nuclear receptors for vitamin D3 and glucocorticoids in prostate stroma cell microenvironment. Biochem Biophys Res Commun. 2011, 412 (1): 13-19. 10.1016/j.bbrc.2011.06.181.View ArticlePubMed
- Campos LT, Brentani H, Roela RA, Katayama ML, Lima L, Rolim CF, Milani C, Folgueira MA, Brentani MM: Differences in transcriptional effects of 1α,25 dihydroxyvitamin D3 on fibroblasts associated to breast carcinomas and from paired normal breast tissues. J Steroid Biochem Mol Biol. 2013, 133: 12-24.View ArticlePubMed
- Chen J, Bardes EE, Aronow BJ, Jegga AG: ToppGene Suite for gene list enrichment analysis and candidate gene prioritization. Nucleic Acids Res. 2009, 37: W305-311. 10.1093/nar/gkp427.PubMed CentralView ArticlePubMed
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102 (43): 15545-15550. 10.1073/pnas.0506580102.PubMed CentralView ArticlePubMed
- Wang TT, Tavera-Mendoza LE, Laperriere D, Libby E, MacLeod NB, Nagai Y, Bourdeau V, Konstorum A, Lallemant B, Zhang R: Large-scale in silico and microarray-based identification of direct 1,25-dihydroxyvitamin D3 target genes. Mol Endocrinol. 2005, 19 (11): 2685-2695. 10.1210/me.2005-0106.View ArticlePubMed
- Kovalenko PL, Zhang Z, Cui M, Clinton SK, Fleet JC: 1,25 dihydroxyvitamin D-mediated orchestration of anticancer, transcript-level effects in the immortalized, non-transformed prostate epithelial cell line, RWPE1. BMC Genomics. 2010, 13: 11-26.
- Ramagopalan SV, Heger A, Berlanga AJ, Maugeri NJ, Lincoln MR, Burrell A, Handunnetthi L, Handel AE, Disanto G, Orton SM: A ChIP-seq defined genome-wide map of vitamin D receptor binding: associations with disease and evolution. Genome Res. 2010, 20 (10): 1352-1360. 10.1101/gr.107920.110.PubMed CentralView ArticlePubMed
- Shalhoub V, Shatzen EM, Ward SC, Young JI, Boedigheimer M, Twehues L, McNinch J, Scully S, Twomey B, Baker D: Chondro/osteoblastic and cardiovascular gene modulation in human artery smooth muscle cells that calcify in the presence of phosphate and calcitriol or paricalcitol. J Cell Biochem. 2010, 111 (4): 911-921. 10.1002/jcb.22779.PubMed CentralView ArticlePubMed
- Yu DM, Yao TW, Chowdhury S, Nadvi NA, Osborne B, Church WB, McCaughan GW, Gorrell MD: The dipeptidyl peptidase IV family in cancer and cell biology. FEBS J. 2010, 277 (5): 1126-1144. 10.1111/j.1742-4658.2009.07526.x.View ArticlePubMed
- Biskobing DM, Fan D, Fan X, Rubin J: Induction of carbonic anhydrase II expression in osteoclast progenitors requires physical contact with stromal cells. Endocrinology. 1997, 138 (11): 4852-4857. 10.1210/en.138.11.4852.PubMed
- Quelo I, Machuca I, Jurdic P: Identification of a vitamin D response element in the proximal promoter of the chicken carbonic anhydrase II gene. J Biol Chem. 1998, 273 (17): 10638-10646. 10.1074/jbc.273.17.10638.View ArticlePubMed
- Swami S, Krishnan AV, Wang JY, Jensen K, Horst R, Albertelli MA, Feldman D: Dietary vitamin D3 and 1,25-dihydroxyvitamin D3 (calcitriol) exhibit equivalent anticancer activity in mouse xenograft models of breast and prostate cancer. Endocrinology. 2012, 153 (6): 2576-87. 10.1210/en.2011-1600.PubMed CentralView ArticlePubMed
- Matthews D, LaPorta E, Zinser GM, Narvaez CJ, Welsh J: Genomic vitamin D signaling in breast cancer: Insights from animal models and human cells. J Steroid Biochem Mol Biol. 2010, 121 (1–2): 362-367.PubMed CentralView ArticlePubMed
- Kemmis CM, Welsh JE: Mammary epithelial cell transformation is associated with deregulation of the vitamin D pathway. J Cell Biochem. 2008, 105: 980-988. 10.1002/jcb.21896.PubMed CentralView ArticlePubMed
- Modlich O, Prisack HB, Munnes M, Audretsch W, Bojar H: Immediate gene expression changes after the first course of neoadjuvant chemotherapy in patients with primary breast cancer disease. Clin Cancer Res. 2004, 10 (19): 6418-6431. 10.1158/1078-0432.CCR-04-1031.View ArticlePubMed
- Hannemann J, Oosterkamp HM, Bosch CA, Velds A, Wessels LF, Loo C, Rutgers EJ, Rodenhuis S, van de Vijver MJ: Changes in gene expression associated with response to neoadjuvant chemotherapy in breast cancer. J Clin Oncol. 2005, 23 (15): 3331-3342. 10.1200/JCO.2005.09.077.View ArticlePubMed
- Folgueira MAAK, Brentani H, Carraro DM, Barros MD, Katayama MLH, De Abreu APS, Barbosa EM, De Oliveira CT, Patrao DFC, Mota LD: Gene expression profile of residual breast cancer after doxorubicin and cyclophosphamide neoadjuvant chemotherapy. Oncol Rep. 2009, 22 (4): 805-813.
- Lee HJ, Liu H, Goodman C, Ji Y, Maehr H, Uskokovic M, Notterman D, Reiss M, Suh N: Gene expression profiling changes induced by a novel Gemini Vitamin D derivative during the progression of breast cancer. Biochem Pharmacol. 2006, 72 (3): 332-343. 10.1016/j.bcp.2006.04.030.View ArticlePubMed
- Towsend K, Trevino V, Falciani F, Stewart PM, Hewison M, Campbell MJ: Identification of VDR-responsive gene signatures in breast cancer cells. Oncology. 2006, 71 (1–2): 111-123.View ArticlePubMed
- Schauber J, Oda Y, Buchau AS, Yun QC, Steinmeyer A, Zugel U, Bikle DD, Gallo RL: Histone acetylation in keratinocytes enables control of the expression of cathelicidin and CD14 by 1,25-dihydroxyvitamin D3. J Invest Dermatol. 2008, 128 (4): 816-824. 10.1038/sj.jid.5701102.View ArticlePubMed
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/13/119/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.