Open Access Research article

Transcriptional effects of 1,25 dihydroxyvitamin D3 physiological and supra-physiological concentrations in breast cancer organotypic culture

Cintia Milani1, Maria Lucia Hirata Katayama1, Eduardo Carneiro de Lyra12, 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 Folgueira1*

Author Affiliations

1 Disciplina de Oncologia, LIM24, Departamento de Radiologia e Oncologia, Faculdade de Medicina da Universidade de São Paulo, Av. Dr. Arnaldo 455,Sala 4124, São Paulo, SP, 01246-903, Brazil

2 Instituto Brasileiro de Controle do Câncer, IBCC, São Paulo, Brazil

3 Department of Environmental Health Sciences, School of Public Health, University at Albany, Rensselaer, NY, 12144, USA

4 Hospital do Câncer A. C. Camargo, São Paulo, Brazil

5 Instituto Adolfo Lutz, Centro de Patologia, Núcleo de Patologia Quantitativa, São Paulo, Brazil

6 Post-doctoral fellow, Viral Vectors Laboratory, Instituto do Câncer do Estado de São Paulo, ICESP, São Paulo, Brazil

For all author emails, please log on.

BMC Cancer 2013, 13:119  doi:10.1186/1471-2407-13-119


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2407/13/119


Received:11 May 2012
Accepted:8 March 2013
Published:15 March 2013

© 2013 Milani et al.; licensee 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.

Abstract

Background

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.

Methods

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.

Results

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.

Conclusions

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.

Keywords:
Breast cancer; Calcitriol; Gene expression; Organotypic culture

Background

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 [10].

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 [16]. 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 [19] while weekly pulses of oral calcitriol allow higher dose administration and peak serum concentrations of 1–15 nM [20]. 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 [21] 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, [24] 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.

Methods

Patients

Post-menopausal breast cancer patients clinical stages I-III were invited to take part in the study. This protocol was carried out in compliance with the Helsinki Declaration and was approved by the Institutional Ethics Committee (Comitê de Ética do Instituto Brasileiro de Controle do Câncer, protocol number 108/2006/7; Comitê de Ética em Pesquisa do Hospital das Clínicas da Faculdade de Medicina da Universidade de São Paulo, protocol number 626/06; Comitê de Ética do Hospital do Câncer A. C. Camargo, protocol number 1131/08). A written informed consent was signed by all participants. Twenty one patients were prospectively accrued at Instituto Brasileiro de Controle do Câncer and Hospital do Câncer A. C. Camargo, São Paulo, from August 2007 to September 2009. Characteristics of these patients are described on Table 1.

Table 1. Characteristics of patients

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 webcite) 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 webcite) and were considered significant if P < 0.05 [41].

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 [42]. 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.

Additional file 1: Table S1. Primers.

Format: DOC Size: 51KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Western blot

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).

Immunocytochemistry

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.

Immunohystochemistry

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.

Statistics

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).

Results

Patients characteristics

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

At first, the transcriptional effects of 0.5 nM 1,25(OH)2D3 vs control in breast cancer slices were compared, using SAM paired analysis (FDR ≤ 0.1). As shown in Table 2, seven genes were up-regulated and two genes were down-regulated in tumor slices exposed to 0.5nM 1,25(OH)2D3 for 24 h. Enrichment of genes involved in vitamin metabolic process (TKTL1, CYP24A1, CYP26B1) was observed. Unsupervised clustering of the differentially expressed genes identified two branches, however there was no aggregation of samples according to 1,25(OH)2D3. treatment (Figure 1). At a more stringent FDR level (≤ 0.01), only five (CYP24A1, DPP4, EFTUD1, FCGR2C, SAMSN1) genes were differentially expressed.

thumbnailFigure 1. Unsupervised hierarchical clustering of breast cancer slices exposed to vehicle or 0.5nM 1,25(OH)2D3. Five tumor samples (6, 11, 12, 13, 15) were cultured with or without 0.5nM 1,25(OH)2D3 for 24 h. Gene expression was evaluated by microarray and analysis performed using SAM paired analysis (FDR ≤ 0.1). The upper colored bar and values indicate gene expression in target samples (i.e., red, more expressed and green, less expressed). Colored lines of the dendrogram and numbers stand for the support for each clustering, black and gray, more reliable; yellow and red, less reliable. Nine genes were differentially expressed (plus one probe without a gene name, NA). Two branches are identified, one of them including three control samples.

Table 2. Genes differentially modulated in breast tumor slices incubated with 0.5nM 1,25(OH)2D3

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).

Additional file 2: Table S2. Gene sets enriched in breast tumor slices incubated in 0.5nM 1,25(OH)2D3.

Format: DOC Size: 39KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Next, the effects of a 24 h exposure to 1,25(OH)2D3 100nM was evaluated in this model system. Using FDR cut offs of ≤ 0.1 and ≤ 0.01, 196 (136 more and 60 less expressed in treated samples) and 30 (28 more and 2 less expressed in treated samples) candidate target genes were found, respectively (Additional file 3: Table S3). Up-regulated genes were involved in vitamin metabolic process, regulation of leukocyte mediated immunity and positive regulation of alpha-beta T cell activation. In addition, some genes were associated with calcitriol and TGF-beta signaling pathway (Table 3). There was no separation of treated and untreated samples, upon unsupervised hierarchical cluster analysis, and paired tumors co-aggregated in the same branch (Additional file 4: Figure S1).

Additional file 3: Table S3. Genes differentially modulated in breast tumor slices incubated in 100 nM 1,25(OH)2D3.

Format: DOC Size: 185KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Additional file 4. Unsupervised hierarchical clustering of breast cancer tissue slices exposed to vehicle or 1,25(OH)2D3 100 nM.

Format: PDF Size: 621KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Table 3. Functional categories of genes up-regulated in breast tumor slices incubated in 100 nM 1,25(OH)2D3

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.

The reproducibility of the present gene list was further tested against gene lists determined in other cell lines. To this end, vitamin D (100 nM) up-regulated genes were cross checked in breast cancer slices and derived fibroblasts [40], squamous carcinoma [43] immortalized prostate [44], and lymphoblastoid cell lines [45] as well as in carotid artery smooth muscle cells [46] (Figure 2). These cell lines were treated with supra-physiological concentrations, ranging from 10-100nM, of 1,25(OH)2D3 or EB1089 (vitamin D analog) for 12–36 hours. This analysis revealed CYP24A1 as the universal vitamin D target gene in all cell types. Expression of CLMN, EFTUD1 and SERPINB1 was up-regulated in five of the six studies and BMP6, CD14, FAM20C, and THBD in four studies. CA2, CILP, CYP19A1, DCBLD1, DPP4, FOXF1, G0S2, GRK5, IL1RL1, KCNK3, SEMA6D and SLC1A1 were up-regulated in another two studies, in addition to the present one. Many of these genes were also regulated by vitamin D in this organotypic culture.

thumbnailFigure 2. Overlapping genes up-regulated by 1,25(OH)2D3 10-100nM in cell lines. Venn diagram of 1,25(OH)2D3 up-regulated genes in breast cancer slices (A) cancer associated fibroblasts and normal adjacent fibroblasts (B) [40]; non-transformed prostate epithelial cell line (C) [44]; human coronary artery smooth muscle cells (D) [46]; SCC25 cells (floor of the mouth/base of the tongue squamous tumor, (E) [43] and lymphoblastoid cell lines (F) [45]. B, C, and D include non-transformed cells. Number of genes regulated in two cell lines appears on the diagram; number of genes regulated in more than two cell lines appears on the text.

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).

Additional file 5: Table 4. Correlation of gene expression values (microarray vs qPCR) evaluated in breast cancer slices.

Format: DOC Size: 44KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Subsequently, the expression of these seven genes was determined in samples from an additional group of 16 patients (validation subset). In these samples, CYP24A1, DPP4 and CA2 were up-regulated by both 1,25(OH)2D3 0.5 and 100nM whereas CD14 expression was induced only by 1,25(OH)2D3 100nM (Figure 3). Median expression of IL1RL1, SHE, and BMP6 was not significantly up-regulated by either dose of 1,25(OH)2D3 in these additional samples, even though elevated mRNA levels were detected in a subset of tumors after treatment.

thumbnailFigure 3. Effect of 1,25(OH)2D3 on expression of CYP24A1, CA2, CD14, DPP4, SHE, BMP6 and ILRL1 in breast cancer slices. Tumor specimens from 12–16 patients were treated or not (control) with 1,25(OH)2D3 (0.5nM or 100nM) for 24 hours and analyzed by qPCR. Relative expression was calculated as 2-ΔΔCT, using GAPDH as internal control and the average value of the target gene in control samples, as reference. Horizontal line represents the median value of gene expression. Friedman test, p value inside the box; * (over the horizontal bracket): p ≤ 0.05, Wilcoxon signed ranks test.

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.

In five independently-derived primary cultures of cancer-associated fibroblasts, CYP24A1 expression was consistently induced in response to 1,25(OH)2D3 0.5nM indicating active VDR signaling in the tumor stroma. However, none of the other target genes, identified in the microarray analysis, were significantly up-regulated in tumor fibroblasts cultured with 0.5nM 1,25(OH)2D3ex vivo , even though there was a trend towards up-regulation of CA2, IL1RL1 and DPP4 (Figure 4).

thumbnailFigure 4. Effect of 0.5nM 1,25(OH)2D3 on expression of CYP24A1, CD14, CA2, DPP4, ILRL1 and BMP6 in cancer infiltrating fibroblasts and mammary epithelial cells. Cancer infiltrating fibroblasts (n=5) and breast derived cell lines (left to right: HB4A, C5.4, SKBR3, MCF7, MDA-MB231) were cultured in the presence or absence of 0.5nM 1,25(OH)2D3 for 24 h. Three independent assays for each cell line were performed and the mean relative gene expression value was calculated. Gene expression was determined by RT-qPCR and relative expression is shown on the Y axis (which is in log scale for CYP24A1, CA2, CD14, DPP4 and IL1RL1). p values: Wilcoxon signed ranks test. White bar: control samples; Gray bars: 1,25(OH)2D3 treated samples.

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.

Tumor slices from seven patients (in FFPE, six samples cultured in the absence or presence of 1,25(OH)2D3 100nM and one sample in the presence of 1,25(OH)2D3 0.5nM) were available for immunohistochemistry. CD14 moderate cytoplasmic staining was observed in at least 50% of tumor cells and weak staining of 10% of the fibroblasts. No differences could be detected between 1,25(OH)2D3 treated and untreated tumor samples (Figure 5).

thumbnailFigure 5. CD14 expression in breast cancer samples. Tumor slices were cultured in the absence (control: B, E) or presence of 1,25(OH)2D3 100nM (C) or 0.5nM (F) for 24 h. Protein expression was examined by immunohistochemistry. CD14 expression was detected as a cytoplasmic staining. Negative reaction with secondary antibodies exclusively, appears on the left (A, D). scale bar: 25μm.

CD14 may be either soluble (sCD14) or membrane-bound (mCD14). There is evidence that sCD14 may be detected in plasma samples from breast cancer patients, hence we determined whether sCD14 concentration might be regulated in the culture medium of breast cancer slices. Although no significant statistical differences were found between control and calcitriol treated samples (0.5nM and 100nM), there was a trend towards higher values of sCD14 in 3/4 samples exposed to 1,25(OH)2D3 100nM (Figure 6).

thumbnailFigure 6. Soluble CD14 in the conditioned medium of breast cancer slices. Tumor specimens from four patients with invasive ductal carcinoma (not included in previous analysis of microarray or PCR) were cultured with or without (control, white bars) 1,25(OH)2D3 0.5 (hatched bars) or 100 nM (gray bars) for 24 h. sCD14 concentration was evaluated in the conditioned medium by ELISA (p=0.42; Friedman’s test).

In mammary epithelial cell lines and cancer associated fibroblasts protein expression was analyzed through western blot and immunocytochemistry. CD14 was more expressed in MCF7 and C5.2 treated with 0.5nM 1,25(OH)2D3, as compared with untreated cells, however no differences in CD14 expression were observed in fibroblasts exposed to 1,25(OH)2D3. In MCF7 cells, DPP4 (CD26) expression was also induced after 1,25(OH)2D3 0.5nM exposure (Figure 7). In addition, in fibroblasts, CA2 as well as DPP4 cytoplasmic immunoexpression was more intense in 1,25(OH)2D3 treated than in control cells (Figure 8).

thumbnailFigure 7. CD14 and DPP4 (CD26) expression in mammary epithelial cell lines and fibroblasts. MCF7 and C5.2 cells as well as fibroblasts were cultured in the absence (control, C) or presence of 1,25(OH)2D3 0.5 nM (D 0.5). Protein expression was analyzed through western blot. Two assays were performed with similar results.

thumbnailFigure 8. CA2 and DPP4 (CD26) expression in fibroblasts. Cancer associated fibroblasts from two patients with IDC were cultured in the presence or absence (control-1st row) of 1,25(OH)2D3 0.5 nM (2nd row) or 100 nM (3rd row) and processed for immunostaining with specific antibodies and DAPI for nuclear staining (on the left). (A) CA2 cytoplasmic granular staining and (B) DPP4 cytoplasmic staining appears in green. scale bar: 25μm.

Discussion

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[46] 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 [47]. 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 [50]. 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 [51]. 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 [51], 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 [40]. 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 [52]. 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.

Conclusion

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.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

CM designed the study, participated in data acquisition, participated in data analysis and interpretation, drafted and revised the manuscript; MLHK designed the study, participated in data analysis and interpretation, drafted and revised the manuscript; ECL, designed the study and participated in data acquisition; JW participated in data analysis and interpretation, drafted and revised the manuscript; LTC participated in data acquisition; MMB participated in data analysis and interpretation, drafted and revised the manuscript; MSM was responsible for the study in Hosp. A. C. Camargo and participated in data acquisition; RAR participated in data acquisition; PRV participated in data acquisition; JCGSG was responsible for the study in IBCC and participated in data acquisition; SN participated in data acquisition; RET participated in data acquisition; MAAKF designed the study, participated in data analysis and interpretation, drafted and revised the manuscript. All authors read and approved the final manuscript.

Acknowledgments

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).

References

  1. 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. PubMed Abstract | Publisher Full Text OpenURL

  2. 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. PubMed Abstract | Publisher Full Text OpenURL

  3. 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. PubMed Abstract | Publisher Full Text OpenURL

  4. 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. Publisher Full Text OpenURL

  5. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. 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 Abstract | Publisher Full Text OpenURL

  7. 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. PubMed Abstract | Publisher Full Text OpenURL

  8. 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. PubMed Abstract | Publisher Full Text OpenURL

  9. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. 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. PubMed Abstract | Publisher Full Text OpenURL

  11. 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 Abstract | Publisher Full Text OpenURL

  12. 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. PubMed Abstract | Publisher Full Text OpenURL

  13. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. 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. PubMed Abstract | Publisher Full Text OpenURL

  15. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. 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 Abstract | Publisher Full Text OpenURL

  17. 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. PubMed Abstract | Publisher Full Text OpenURL

  18. 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. PubMed Abstract | Publisher Full Text OpenURL

  19. 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. PubMed Abstract | Publisher Full Text OpenURL

  20. 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. PubMed Abstract | Publisher Full Text OpenURL

  21. 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. PubMed Abstract | Publisher Full Text OpenURL

  22. 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. PubMed Abstract | Publisher Full Text OpenURL

  23. 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. PubMed Abstract | Publisher Full Text OpenURL

  24. 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. PubMed Abstract | Publisher Full Text OpenURL

  25. 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. PubMed Abstract | Publisher Full Text OpenURL

  26. 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. OpenURL

  27. 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. Publisher Full Text OpenURL

  28. 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. PubMed Abstract | Publisher Full Text OpenURL

  29. 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 Abstract | Publisher Full Text OpenURL

  30. 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. PubMed Abstract | Publisher Full Text OpenURL

  31. 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. PubMed Abstract | Publisher Full Text OpenURL

  32. 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. PubMed Abstract | Publisher Full Text OpenURL

  33. 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. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  34. 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. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  35. 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. PubMed Abstract | Publisher Full Text OpenURL

  36. 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. PubMed Abstract | Publisher Full Text OpenURL

  37. 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. PubMed Abstract | Publisher Full Text OpenURL

  38. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. 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. PubMed Abstract | Publisher Full Text OpenURL

  40. 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. PubMed Abstract | Publisher Full Text OpenURL

  41. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  42. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  43. 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. PubMed Abstract | Publisher Full Text OpenURL

  44. 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. OpenURL

  45. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  46. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  47. 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. PubMed Abstract | Publisher Full Text OpenURL

  48. 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. PubMed Abstract | Publisher Full Text OpenURL

  49. 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. PubMed Abstract | Publisher Full Text OpenURL

  50. Swami S, Krishnan AV, Wang JY, Jensen K, Horst R, Albertelli MA, Feldman D: Dietary vitamin D₃ and 1,25-dihydroxyvitamin D₃ (calcitriol) exhibit equivalent anticancer activity in mouse xenograft models of breast and prostate cancer.

    Endocrinology 2012, 153(6):2576-87. PubMed Abstract | Publisher Full Text OpenURL

  51. 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 Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  52. Kemmis CM, Welsh JE: Mammary epithelial cell transformation is associated with deregulation of the vitamin D pathway.

    J Cell Biochem 2008, 105:980-988. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  53. 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. PubMed Abstract | Publisher Full Text OpenURL

  54. 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. PubMed Abstract | Publisher Full Text OpenURL

  55. 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. PubMed Abstract | Publisher Full Text OpenURL

  56. 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. PubMed Abstract | Publisher Full Text OpenURL

  57. 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. PubMed Abstract | Publisher Full Text OpenURL

  58. 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. PubMed Abstract | Publisher Full Text OpenURL

Pre-publication history

The pre-publication history for this paper can be accessed here:

http://www.biomedcentral.com/1471-2407/13/119/prepub