Abstract
Retinoic acid receptor β (RARβ) has been proposed to act as a tumor suppressor in breast cancer. In contrast, recent data have shown that RARβ promotes ERBB2-induced mammary gland tumorigenesis through remodeling of the stromal compartment and activation of cancer-associated fibroblasts. However, it is currently unknown whether RARβ oncogenic activity is specific to ERBB2-induced tumors, or whether it influences the initiation and progression of other breast cancer subtypes. Accordingly, we set out to investigate the involvement of RARβ in basal-like breast cancer using mouse mammary tumor virus (MMTV)-wingless-related integration site 1 (Wnt1)-induced mammary gland tumorigenesis as a model system. We found that compared with wild type mice, inactivation of
Introduction
Retinoic acid (RA), which is derived from dietary vitamin A, plays an essential role in the maintenance of the differentiated state of epithelial cells and tissue remodeling [1–3]. The family of retinoic acid receptors (RARs) that includes the RARα, β and γ subtypes and associated alternatively spliced isoforms are the main mediators for the biologic effects of RA [4, 5]. RA and other natural and synthetic RAR agonists are well known to inhibit the growth of human breast tumor cell lines in vitro and in mouse xenograft models, a process that generally involves the induction of both apoptosis and cell-cycle arrest [6–12]. In addition, the absence of RARβ-dependent retinoid signaling has been associated with the development of carcinomas, including breast cancer, and resistance to the anti-cancer activity of RA [13]. Despite the potent growth-suppressing activity of retinoids on human breast cancer cell lines, this class of compounds has not
shown relevant anti-cancer activity in the treatment of breast cancer [14, 15].
Recently, RAR-null mice have been used to address the role of the RARs in mammary gland development and the response of the gland to specific oncogenic insults [16, 17]. Ablation of RARα1 led to an increase in the density of the mammary epithelial tree and the content of luminal mammary progenitors, in addition to reducing the size of the mammary stem cell-containing compartment [17]. In contrast, loss of RARβ resulted in a delay in mammary tree development in the pubertal stage as a significant shortening of the distance between the lymph node and the terminal end buds were observed in the RARβ-null gland [16]. The developmental abnormalities observed in both RAR-deficient mouse models were transient since, as previously reported [18–20], the glands were fully functional at later stages of development. While ablation of RARα1 and β resulted in distinct developmental mammary gland abnormalities, the absence of RA signaling transduced by either receptor resulted in a delay in mammary gland tumorigenesis. In the RARα1-null model, the reduction observed in the mammary stem cell-containing compartment might be a factor in the delay in tumor formation by the Wnt1 oncogene, but the exact mechanism responsible for the prooncogenic property of RARα1 in this context remains to be determined. On the other hand, absence of RARβ led to a remodeling of the stroma during tumor progression induced by the Erbb2 oncogene that was marked by a decrease in angiogenesis, in the recruitment of inflammatory cells and in the number of myofibroblasts [16]. In agreement with this observation, tissue recombination experiments demonstrated that the presence of RARβ in the stromal compartment is essential for the growth of mammary carcinoma, and that RA-dependent activity of the Cxcl12/Cxcr4/Erbb2 signaling axis played a significant role in this process. Together, these studies suggest that retinoid-based approaches for the prevention and treatment of breast cancer should be re-evaluated and that a better understanding of the mechanisms of retinoid action in mammary gland tumorigenesis is required to achieve this goal.
Estrogen receptor β (ERβ) is an attractive target for drug development. While orginal models of estrogen action suggested that only a single ER gene (which encodes a protein that is now termed ERα) was responsible for transducing signals of estradiol and other ligands [11], the discovery in 1996 of a second ER gene, encoding ERβ, prompted a reevaluation of the estrogen signaling system. It is now known that ERα and ERβ play different roles in gene regulation [12] and that ERα and ERβ have overlapping but distinctive tissue distributions and non-redundant roles [13]. These considerations have led to the suggestion that ERβ could be an attractive therapeutic target for the development of selective agonists to treat and prevent neurodegenerative diseases [14] and other diseases, including autoimmune diseases, endometriosis, depression, hypertension, and colon, breast, prostate, lung, and skin cancer [15]. It is important that such ligands should not cross-react with ERα, which triggers classical estrogenic side effects, such as breast or uterine stimulation, thereby increasing a woman's chance of developing breast or uterine cancer, and gynecomastia and decreased libido in men [16]. Presently, natural and synthetic estrogens for ERβ are being studied in colon cancer, breast cancer, lung cancer, schizophrenia, and metabolic syndrome [17]. While the agonist ERB-041 failed to demonstrate efficacy in a Phase II double-blind clincial trial for rheumatoid arthritis, further studies are warranted to examine ERβ-selective efficacy in other inflammatory disorders [18]. With many other clinical trials still in progress awaiting completion of the study followed by publication of the findings, it is still too early to make any definitive conclusions about these drugs and their effects. The discovery of potentially beneficial effects of selective ERβ ligands on prostate cancer proliferation and apoptosis in the absence of full ERα or ERβ agonism has raised hopes that applications for new safe selective ERβ modulator ligands could emerge in the context of this disease [19–22].
In the current study, we further tested the specificity of the RAR/oncogene interaction in the development of mammary gland carcinoma using bigenic mice generated by crossing the RARβ-null mice with the MMTV-Wnt1 transgenic mice. The Wnt1 oncogene targets mammary stem cells or early progenitor cells as transformation cells [21]. Accordingly, Wnt1-induced tumors show a co-existence of both luminal cells and basal/myoepithelial cells, which is reminiscent of human basal-like breast cancer that is associated with poor prognosis and lack of effective treatment [22]. Wnt1-induced tumors have also been reported to have abundant activated stroma that correlated with more aggressive tumor development [23]. Therefore, we sought to investigate the influence of RARβ on Wnt1-induced mammary carcinoma. We found that inactivation of RARβ resulted in a protective effect against Wnt1-induced mammary gland tumorigenesis, abrogation of Wnt signaling in both the epithelial and stromal compartments, and suppression of epithelial-mesenchymal transition (EMT) in the tumors. Our work indicates therefore that RARβ is essential for the bi-directional interaction between tumor and stromal cells and that RARβ specific antagonists may represent a novel therapeutic strategy to treat poor outcome breast cancer.
Materials and Methods
Mice
All animals were FVB/NJ background and housed in a pathogen-free facility at McGill University and were given a standard rodent diet and water ad libitum. MMTV-Wnt1 mice were purchased from the Jackson Laboratory [24]. To determine the initiation of mammary tumor formation, mice were palpated weekly (sensitive to 1-mm tumor masses). All mouse manipulations were performed in accordance with the McGill Facility Animal Care Committee and the Canadian Council on Animal Care.
Cell culture
Mouse mammary fibroblast cell lines were generated from the MMTV-Wnt1 animals as described before [16]. To generate tumor cells, tumors were dissected, mechanically dissociated, and forced through 40 μm mesh. Viable cells were plated, grown in DMEM supplemented with 10% FBS and penicillin/streptomycin. The cells had virtually identical epithelial morphology under microscope and had similar cell surface marker expression. All cultured cells were maintained in DMEM supplemented with 10% FBS and incubated in 5% CO2 at 37oC. The fibroblasts were treated with all-trans RA (Sigma #R2625) at a concentration of 0.5 μM for 2 days to generate conditioned media. Conditioned media from
Reactive stroma index
Hematoxylin and eosin staining was carried out in the histology service core of the Goodman Cancer Research Centre. Stained sections were examined in a blinded manner and scored for the reactive stroma index based on the percentage contribution of a stroma area in tumor mass (i.e., 0–10% stroma area = 0, 11–20% stroma area = 1, and so on, >50% stroma area = 5). The average of five areas was used as the reactive stroma index for each sample.
Histological analyses
Immunohistochemical analyses were performed on 4-μm formalin-fixed paraffin-embedded sections of tumors from mice killed one month after tumor initiation. Immunostaining was done using the VECTASTAIN avidin-biotin complex kit (Vector Laboratories) as previously described [46]. Immunofluorescence staining was conducted on formalin-fixed paraffin-embedded tumor sections. Five random areas per section from 3 separate sections obtained from different mice for each group were analyzed. The following dilutions of primary antibodies were used: rabbit polyclonal anti–human-RARβ (Abcam), 1:100; mouse monoclonal anti-human Ki67 (clone B56; BD Pharmingen), 1:100; rabbit polyclonal anti-mouse Wnt1 (Upstate), 1:50; rabbit polyclonal anti–human-caspase-3 (Cell Signaling), 1:50; rat monoclonal anti–mouse-CD31 (clone Mec13.3; Biocare), 1:50; mouse monoclonal anti-CD45, 1:100, and rabbit polyclonal anti-mouse collagen I, 1:100 (Abcam Inc.), horse polyclonal anti-human smooth muscle actin-α (α-SMA) (Abcam), 1:100,; rabbit polyclonal anti-mouse IGF1 (Abcam), 1:100. For the negative control, primary antibody was replaced with non-immunized rabbit or mouse IgG (Vector Laboratories). All images were taken with 10x objectives (100x magnification). For multiple antigen labeling the VECTASTAIN System (Vector ABCDAB kit, Vector ABC-AP kit, and Vector ImmPRESSNovaRED and VIP kit) (Vector Laboratories) was used. Peroxidase and alkaline phosphatase substrates were used to develop color of the antigens with different localization. To compare positive cell numbers, 10 random fields per section were documented in the Aperio Image System and were analyzed using the Spectrum software (Aperio Technologies). Mean values shown beside the corresponding images were determined from results from at least six different mice for each genotype.
Apoptosis analysis
The number of TUNEL-positive or caspase 3-postive tumor cell nuclei was calculated relative to the total number of 1,000 cancer cell nuclei, which was analyzed with the Spectrum software. Mean values were determined from results from at least six different mice.
Western blotting
Lysates were prepared from the mammary tumors of three
Laser capture microdissection, RNA extraction and linear amplification
All tissues included in this study were re-examined by an animal pathologist dedicated to the project. Tissue specimens were micro-dissected into epithelium and stroma using a PixCell IIe laser micro-dissection (LCM) system (Arcturus). All micro-dissections were performed within two hours following tissue staining. Total RNA was extracted from each population of micro-dissected cells using a guanidinium isothiocyanate (GITC) extraction protocol. Briefly, LCM caps were incubated for 5 minutes (room temperature) in 200 μl GITC extraction buffer (4 M GITC, 25 mM sodium citrate pH 7.0, 0.1 M β-mercaptoethanol, 0.5% N-lauroylsarcosine) supplemented with 1.6 μl β-mercaptoethanol. Subsequently, 20 μl of 2 M NaOAc, pH 4.0, 220 μl of water-saturated phenol and 60 μl of chloroformisoamyl alcohol (23:1) were added to the extraction buffer. Following a 15 minute incubation on ice and centrifugation (12,000 rpm, 15 minutes) the aqueous phase was removed and RNA was precipitated with 2 μl glycogen (GenHunter, Nashville, Tennessee, USA) and 200 μl isopropanol. Samples were placed at −80°C for 30 minutes and centrifuged at 4 °C (12,000 rpm) for 30 minutes to pellet RNA. Pellets were washed with 70% ethanol, air-dried and subjected to DNase I treatment (Roche, Basel, Switzerland). DNAse treatment was performed in the presence of an RNase inhibitor (Invitrogen, Carlsbad, California, USA). Subsequently, samples were re-extracted as described above and re-suspended in 10 μl of diethylpyrocarbonate-treated water. RNA was quantified using a RiboGreen assay (Molecular Probes, Carlsbad, California, USA). Subsequently, 2 to 4 ng of total RNA was subjected to two rounds of T7 linear amplification using Ambion Amino Allyl MessageAmp kit (Ambion, Austin, Texas, USA) and labeled with Cy3 and Cy5 dyes according to the manufacturer's procedure. Prior to microarray hybridizations, amplified products were quantified using a spectrophotometer (Nanodrop, Wilmington, Delaware, USA) and analyzed using a BioAnalyzer to assay for quality (Agilent Technologies, Santa Clara, California, USA).
RNA microarrays
SurePrint G3 mouse Genome 8×60 K arrays (Agilent Technologies, product G4852A) were used for all experiments. RNA samples (500 ng) were subjected to fragmentation followed by 18 hours hybridization, washing, and scanning (Agilent Technologies, model G2505B) according to the manufacturer's protocol (manual ID #G4140-90030). Samples were hybridized against Universal mouse Reference RNA (Stratagene, ID #750600, La Jolla, California, USA). Duplicate hybridizations were performed for all samples using reverse-dye labeling. Microarray data were feature extracted using Feature Extraction Software (v. 7.11) from Agilent with the default parameters. Raw data were uploaded to the NCBI Gene Expression Omnibus database (GSE56391). Outlier features on arrays were flagged by the software. Arrays were required to have an average raw signal intensity of 1,000 in each channel, and a signal to noise ratio above 16 per channel. MvA plots were examined for signs of hybridization or labeling problems. Replicate arrays were required to have a concordance above 0.944. This level was established empirically using sets of known good replicate arrays in our database.
Real-time PCR for Igf1 expression
Total RNA from LCM samples was amplified using Amino Allyl MessageAmp II aRNA Kit (Ambion #1753). cDNA was made from 2 μg of RNA by reverse transcription with Oligo(dT) primer, dNTPs, 5X 1st strand buffer, DTT, RNase inhibitor, and Superscript II RNase H Reverse Transcriptase (Invitrogen). cDNA was amplified by qRT-PCR using the specific forward 5'-TTCTACCTGGCGCTCTGCTTGC and reverse 5'-CCCTCCGAATGCTGGAGCCATA primers, a QuantiTect SYBR Green PCR Kit (Qiagen) and a LightCycler instrument (Roche) following the Qiagen software protocol.
ELISA
Submandibular venous blood was collected into microcentrifuge tubes using a sterile lancet (Medipoint, Mineola, NY) according to the manufacturer's instructions. Whole blood was centrifuged at 2500xg for 20 minutes at 4 °C to coll ect serum, which was then stored at −80°C. Serum was processed for ELISA of IGF-1 (R&D Systems, Minneapolis, MN) according to the manufacturer's protocol. Serum IGF-1 levels are reported in ng/ml. An IGF-1 control sample, with aliquots stored at −80°C, was included on each plate, and all data are reported using simple ratio normalization to the initial reading of the control sample.
Statistical analysis
Two-tailed paired Student t test was calculated by Excel software. Prism software was used for one or 2-way ANOVA and log-rank (Mantel-Cox) test.
Results
Global ablation of Rarb in mice delays Wnt1-induced mammary gland tumorigenesis and alters the stromal compartment
To further study the potential role of RARβ as a tumor promoter in mammary gland tumorigenesis, we used the well-characterized mouse model of human breast cancer expressing the oncogene Wnt1 in combination with the Rarb-null allele [18, 24]. While both Rarb-/-and wild-type mice displayed similar high penetrance for tumor formation (78% and 83%, respectively), the Rarb-/- mice showed a significant delay (∼12 weeks, P<0.01) in tumor formation when strong expression of Wnt1 was driven by the MMTV promoter (Figure 1A). In addition, the tumors in mice with the Rarb-/- genetic background grew at a significantly slower rate than tumors induced in their wild-type siblings (6.3 vs 4.0 weeks, P<0.05) (Figure 1B). Consistent with these results, examination of mammary gland sections taken from MMTV-Wnt1 mice by stroma index than tumors developing in

Ablation of Rarb suppresses Wnt1-induced mammary tumorigenesis and affects the composition of the stroma.
Loss of Rarb results in remodeling of the stromal compartment in Wnt1-induced mammary gland tumors
Stroma may initiate a desmoplastic reaction that includes activation of CAFs and trans-differentiation of tumor or epithelial cells into myofibroblasts, infiltration of immune cells, increased secretion of growth factors and cytokines as well as elevated matrix synthesis and remodeling that manifests as matrix stiffening. Analysis of differential gene expression in the stroma and epithelial compartments using laser capture microdissection (LCM) showed that several genes associated with activation of fibroblasts such as

Loss of Rarb alters the gene expression profile of the stromal compartment and its composition.
Inhibition of epithelial-mesenchymal transition in Wnt1-induced Rarb-null tumors
Besides activation from local fibroblasts, emerging evidence indicates that matrix-producing myofibroblasts also arise through epithelial-mesenchymal transition (EMT) [25]. Interestingly, LCM-mRNA microarray data showed that several genes involved in EMT, such as Ecm2, Mmp3, Fgf2, and Igf1, were significantly down regulated in MMTV-Wnt1/

Inhibition of epithelial-mesenchymal transition in Wnt1-induced Rarb-null tumors.
Down regulation of the IGF1/AKT axis in Wnt1-induced Rarb-null tumors
IGF-1 plays a central role in cell growth, differentiation, survival, transformation and metastasis, and deregulation of the IGF signaling pathway is well recognized as a key contributor to the progression of multiple cancers including mammary gland tumors. The IGF-1 signaling axis can also affect EMT. IGF-1 binds IGF1R to induce PI3K and activates AKT to repress GSK3β. GSK3β is a repressor of β-catenin and Snail able of inhibiting E-cadherin expression and subsequently promotes EMT and tumor progression. While we observed a reduction in the expression of Igf1 in the tumor compartment, IGF-1 is also produced by stromal cells. Western blotting using protein lysates obtained from whole tumors showed a decrease in endogenous IGF-1 levels and an increase in the expression of IGFBP5, an antagonist of IGF signaling, in the

Down regulation of the IGF1/AKT axis in Wnt1-induced Rarb-null tumors.
Discussion
The antitumor activities of retinoids and their receptors have been investigated extensively and one of the receptors, RARβ, has been suggested to possess many of the functional characteristics of a tumor suppressor [13]. However, previous work from our laboratory showed that in the context of the whole organism,
In the Wnt1 model, oncogenic transformation has been proposed to target mammary gland stem cells or early progenitor cells. Therefore, Wnt1-induced tumors are morphologically composed of a variety of cell types, including distinct tumor luminal epithelial, basal/myoepithelial cells and host-derived stroma cells [26–28]. Wnt1 ligand produced by tumor cells can directly activate the surrounding stromal cells in a paracrine manner [29]. In return, stromal cells can affect tumor cells by expressing IGF-1, which activates IGF1R signaling in tumor cells [30]. This important molecular cross talk between tumor and stromal cells is severely weakened in the absence of RARβ. Deletion of
Besides affecting stromal cells, Wnt1 ligand secreted by luminal epithelia can also target myoepithelial cells [29]. A recent report indicates that the Wnt1-driven tumor growth can be inhibited by a soluble Wnt receptor inhibitor [30]. The mechanism through which this effect is mediated involves the induction of IGFBP5, an antagonist of IGF signaling that has been shown to facilitate tumor growth by paracrine communication. The regulation of
Conclusions
The failure of retinoid-based therapies to treat breast cancer led us and others to hypothesize that retinoic acid receptor isoform might have distinct and perhaps antagonistic functions in mammary gland development and in oncogenesis. Indeed, while activation or inhibition of the RARα isoform in mice are consistent with the anti-proliferative action of this receptor [17, 45], the activity of the RARβ and RARγ suggests that these isoforms act as proto-oncogenes in the context of the whole gland [16, 45]. In this work, we extended these findings by demonstrating that

Public Datasets
Gene Expression Omnibus
GSE56391
Footnotes
Acknowledgements
This work was supported in part by grant DK043093 to MRS from the National Institutes of Health. Gene expression microarray data were generated by the Genomics Core Facility and analyzed with assistance from the Biostatistics Core Facility of the USC Norris Comprehensive Cancer Center, supported by Cancer Center Support Grant P30CA014089 from the National Cancer Institute. RC and D-YW were supported in part by National Institutes of Health-funded Training Grants T32 CA009320 and T32 GM067587, respectively.
