- Research article
- Open Access
- Open Peer Review
Wnt signaling in triple negative breast cancer is associated with metastasis
BMC Cancer volume 13, Article number: 537 (2013)
Triple Negative subset of (TN) Breast Cancers (BC), a close associate of the basal-like subtype (with limited discordance) is an aggressive form of the disease which convey unpredictable, and poor prognosis due to limited treatment options and lack of proven effective targeted therapies.
We conducted an expression study of 240 formalin-fixed, paraffin-embedded (FFPE) primary biopsies from two cohorts, including 130 TN tumors, to identify molecular mechanisms of TN disease.
The annotation of differentially expressed genes in TN tumors contained an overrepresentation of canonical Wnt signaling components in our cohort and others. These observations were supported by upregulation of experimentally induced oncogenic Wnt/β-catenin genes in TN tumors, recapitulated using targets induced by Wnt3A. A functional blockade of Wnt/β-catenin pathway by either a pharmacological Wnt-antagonist, WntC59, sulidac sulfide, or β-catenin (functional read out of Wnt/β-catenin pathway) SiRNA mediated genetic manipulation demonstrated that a functional perturbation of the pathway is causal to the metastasis- associated phenotypes including fibronectin-directed migration, F-actin organization, and invasion in TNBC cells. A classifier, trained on microarray data from β-catenin transfected mammary cells, identified a disproportionate number of TNBC breast tumors as compared to other breast cancer subtypes in a meta-analysis of 11 studies and 1,878 breast cancer patients, including the two cohorts published here. Patients identified by the Wnt/β-catenin classifier had a greater risk of lung and brain, but not bone metastases.
These data implicate transcriptional Wnt signaling as a hallmark of TNBC disease associated with specific metastatic pathways.
Stratification of breast cancer (BC) into distinct histological or molecular subtypes has clinical utility for prognosis of outcome and prediction of treatment . Breast cancer subtypes based on clinical or molecular characteristics are typically referred to as hormone-receptor positive (HR+) or luminal, HER2-amplified, and triple negative (TN) or basal-like. TNBC are defined by negative of expression of ER, PR, and HER2 amplification and are associated with a higher grade, undifferentiated metaplastic histology, stem cell-like characteristics, invasiveness, higher metastatic potential, and inconsistently effective therapies. Triple negative and basal-like subtypes have a significant overlap, and lack standardized clinical markers that differentiate the two subtypes, something that underlines the inherent heterogeneous nature of these subtypes [2–4]. Triple negative and basal-like subtypes portend some of the worst prognoses in BC [5, 6], and have the most challenging diagnosis among BC patients due to the potential aggressive nature of the disease and limited number of therapeutic options available. Studies of Shah et al. indicated that understanding the biology and therapeutic responses of patients with TNBC will require the determination of individual tumor clonal genotypes . Despite of making considerable progress in cancer research, the mortality rate of TNBC has remained unchanged in the last decade primarily due to the lack of specific target identification. The allure of the emerging genomic technologies in cancer in their ability to generate new biomarkers that predict how individual patients will respond to various treatments has not been completely successful in TNBC. Only recently, the first comprehensive genomic analysis of a basal-like breast cancer was performed by using massively parallel sequencing technology . Despite a few recent reports that indicated the involvement of certain genes/signaling molecules related to tumorigenic pathways [4, 9–11] in these subsets, there remains an unmet need for an in-depth study to identify driver pathways in these closely associated subtypes of BC.
We analyzed mRNA expression from archival formalin-fixed paraffin-embedded (FFPE) tumor specimens from two breast cancer cohorts. Our data, together with meta-analysis of previous breast cancer microarray studies, indicate Wnt pathway activation in TNBC subtypes and provides evidence for an increased Wnt/β-catenin signaling associated with high grade, poor prognosis, and metastatic disease. In light of reports from Reis-Filho’s team , who established a preferential increase in β-catenin protein (immunohistochemistry) in TNBC patients, our study not only identifies clinical markers associated with Wnt signaling such as histological grade 3 tumors and triple negative pathological subtype, but also indicated an upregulated state of Wnt signaling increasing the risk for brain and lung metastases, thus classifying Wnt signaling as a rational target in TNBC. The inhibition of metastasis-associated phenotypes, integrin-directed migration and invasion following Wnt-antagogist, WntC59 as well as β-catenin SiRNA provided mechanistic proof of concept for the involvement of this pathway in the progression of the disease and its clinical outcome.
Archived FFPE tumor specimens were obtained from St. Mary’s Hospital, Montreal, Quebec, Canada (Quebec cohort) and Grady Memorial Hospital, Atlanta, GA (Georgia cohort) according to institutional guidelines (Emory University). Tumor content (> 50% inclusion criteria) was assessed by a board certified pathologist. Cohort sizes of 107 and 166 patients from St. Mary’s and Grady Hospitals respectively were analyzed at the Emory Biomarker Service Center (Winship Cancer Institute, Atlanta, GA). Both Georgia and Quebec samples (FFPE) were acquired following the acceptance of our protocols by the ethics committees (Emory University, USA for Georgia study and Canada for Quebec study) of the respective institutions. Archived FFPE samples were used for the study. All archived FFPE tumour specimens were obtained from St Mary’s Hospital (Montreal, Quebec, Canada) according to institutional guidelines.
Data Deposition: Microarray data curetted under GEO series accession numbers [GSE17650 & GSE18539].
RNA preparation, quality control, and DASL assay
FFPE samples were analyzed using DASL (cDNA-mediated, Annealing, Selection, Extension, and Ligation) expression chemistry (Illumina, Inc., San Diego, CA), on the Illumina human cancer panel and a custom panel with breast cancer related genes. Depending upon tissue availability, specimens were obtained via either three 5 μm sections or 0.5 mm cores. Deparaffinization, RNA extraction and RNA purification were carried out using commercially available RNA High Pure Kit (Roche, Mannheim, Germany) modified as previously described . Prior to DASL analysis, RNA quality was assessed via RPL13a TaqMan assay using a threshold Ct of less than 29.5. For samples that passed QC, 200 ng of total RNA was used as input for the DASL assay according to the manufacturer’s protocol (Illumina, San Diego, CA). When ample RNA was available, RNA replicates were run to test for reproducibility of the DASL assay.
To determine whether the DASL assay yields comparable data to IHC data, the DASL assay gene intensity (expression) data were compared with the available IHC protein expression data for ER, PR and HER2 on the set of tumor samples. Standard clinical IHC testing was conducted for ER, PR and HER2 according to guidelines based on ASCO/CAP. This is now stated in the Methods section (p. 7, last sentence). “In total, we obtained 87 FFPE breast carcinomas that had previously been scored for the breast cancer markers, ER, PR and HER2 by immunohistochemistry (IHC) according to guidelines based on the ASCO/CAP recommendations for ER, PR and HER2 testing (ER/PR testing [http://www.cap.org/apps/docs/laboratory_accreditation/summary_of_recommendations.pdf]; HER2 testing [http://www.cap.org/apps/docs/committees/immunohistochemistry/summary_of_recommendations.pdf])”. Once standardized, the similar protocol was followed in both the cohorts. The comparison of IHC data with DASL data for ER, PR and HER2 had been carried out as mentioned in our earlier publication. Our data show that the concordance of DASL data with IHC data for all three receptors is very high, which is consistent with our previous published work relating mRNA and IHC protein levels . Once standardized, the similar protocol was followed in both the cohorts. A detailed list of the gene-composition of the Illumina human cancer panel, and the breast cancer related genes has been published earlier by our group . We have added a supplementary table with the gene-list (Additional file 1: Table S1, Additional file 2: Table S2) to the revised-MS for the convenience of the reader. The detailed procedure to test the reproducibility of the DASL assay has been already published by our group  (The specific details are noted in the Additional file 3 of the above mentioned published article). We used 80 replicates for the standard panel, and 100 replicates for breast cancer custom panel.
Western blots were performed by solubilizing cell lines with lysis buffer (150 mM NaCl, 6 mM Na2HPO4, 4 mM NaH2PO4, 2 mM EDTA, 1% sodium deoxycholate, 1% NP-40, 0.1% SDS, 1% aprotinin, 0.2 M sodium orthovanadate, and 0.1 M phenylarsineoxide), and lysates were assayed for total protein (Bio-Rad protein assay kit) using BSA as a standard. The normalized lysates (20–40 μg protein) were resolved by 12.5% SDS-PAGE and transferred to nitrocellulose membranes . Membranes were probed with anti-β-catenin (Abcam Inc, Cambridge, MA) antibody, and visualized with enhanced chemiluminescence reagent combined with peroxidase-conjugated IgG .
In vitroknockdown of β-catenin protein by SiRNA
Breast cancer cell line (MDA-MB231) was seeded onto 6-well tissue culture dishes, and allowed to attach in culture medium supplemented with 10% FBS. A cell density of 60% to 70% was used for the transient transfection (Lipofectamine 2000) of β-catenin-specific SiRNA (Invitrogen, NY; CTNNB1 VHS50819) into MDA-MB231 cells according to the manufacturer’s instructions. Transfected cells were collected after 24, 48, and 72 hours for analyses .
TCF/LEF promoter activity assay
A luciferase-based reporter gene was used to measure promoter activity of the TCF/LEF transcription factor . For SiRNA based study, cells were transiently transfected with beta-catenin SiRNA . After beta-catenin siRNA transfection for 24 hours, the cells were transiently transfected with the reporter construct TOPflash or FOPflash. In brief, cells were co-transfected with 2.5 μg TOP flash, a synthetic luciferase-based promoter plasmid (sensitive to the activity of the β-catenin/ TCF-4 complex, containing three copies of the TCF-4 binding site upstream of a firefly luciferase reporter gene) using the Lipofectamine 2000. In the other set of cells, an equal amount of the mutant form of the above promoter (FOP flash) was co-transfected using the same transfection reagent. FOPflash has mutated copies of Tcf/Lef sites and is used as a control for measuring nonspecific activation of the reporter. Twenty hours after TOPflash or FOPflash transfection, luciferase assay was performed. Relative luciferase activity (in arbitrary units) was reported. In a separate set of experiments, cells were co-transfected either with TOP flash or FOPflash using lipofectamine. After 12 hour incubation, each set was treated with sulindac sulfide for 24 hours. The relative luciferase activity (TOP flash/FOP flash) was calculated from triplicate experiments.
Cell line based phenotypic assays
Fibronectin directed migration assay was performed on Wnt-antagonist, WntC59 (Cellagen Technology, LLC, San Diego, CA) treated or β-catenin SiRNA transfected MDA-MB231 cells by transwell assay and scratch assay. Invasion assay was performed by transwell assay. Haptotaxis assays were carried out using transwell migration chambers (Costar Corp.) as previously described . Cells were added into the upper chamber of the transwell containing the through which they were allowed to migrate over time to the fibronectin-coated side. Control experiments involved coating both sides of the membrane with fibronectin. In vitro wound healing assays were performed as previously described . In brief, after coating plates with fibronectin, wounds were created by scratching the confluent monolayer of cells. The width of the “scratched” area was measured from randomly chosen fields using either Olympus DP72 system or Axiovert 200 M, Zeiss system. Student’s t test was used to determine the statistical significance.
Confocal microscopy and real-time video microscopy of live cells
To study the cytoskeletal arrangement, HCC38 and MDA-MB231TNBC cells were fixed, and permeabilized with PHEMO buffer. Phalloidin 555 was used for staining the cytoskeleton filamentous-actin and DAPI as a counter stain. Cells were imaged using a Zeiss (Thornwood, NY) LSM 510 Meta confocal microscope with a Plan-Apochromat oil objective. Images were acquired using Zeiss LSM 510 software and processed using Adobe Photoshop CS3. To study the involvement of Wnt-pathway in integrin-directed migration in real time, video microscopy was performed. A scratch-would healing assay was performed on the confluent layer of cells (grown on fibronectin-coated glass-cover slip culture- dishes; Mattek, Ashland, MA). Time-lapse images were acquired with a Perkin Elmer Ultraview ERS (Norwalk, CT) disk-spinning confocal system, mounted on a Zeiss Axiovert 200 M inverted microscope equipped with a 37°C stage warmer, incubator, and humidified CO2 perfusion system. Bright-field images were acquired with a Hamamatsu Orca-ER camera with a Plan-Neoflour 10× objective (NA 0.75; 1×1 binning) at 10 minutes intervals for each image set.
Data and statistical analysis
A full description can be found in the supplementary methods. In brief, DASL transcript intensities were quantile normalized in GenomeStudio and replicates were mean combined. Differential transcripts were determined using permutation testing  with a false discovery rate (FDR) less than 1% and a 1.5 fold-change. Hierarchical clustering was performed using the “heatmap.2” function of the R/Bioconductor package “gplots” . KEGG signal transduction pathways  were analyzed for overrepresentation in the triple negative subtype using Fisher’s exact test and pathway expression was determined using the mean of normalized pathway components. Differential pathway regulation was assessed by a t-test with Bonferroni’s correction applied and permutation testing. HMEC oncogenic data  was downloaded from Gene Expression Omnibus (GEO)  and pathway expressions were calculated as mean of differentially expressed genes with induced genes weighted positively and inhibited genes negatively. Analysis of Wnt regulation induced by Wnt3A from Nguyen et al.  employed the same methods for the LWS-81 genes. Meta-analysis data [25–32] were downloaded from GEO  and Affymetrix CEL files were MAS 5.0 normalized with a target intensity of 600. Agilent normalized series matrix files were downloaded and duplicate samples between studies were removed. A nearest shrunken centroid classifier implemented in the R/Bioconductor package “pamr”  was trained on β-catenin induced data from Bild et al.  and applied to 11 studies and 1,878 patients. Overrepresented pathological subtype, intrinsic subtype, histological grade, and lymph node status in Wnt + (Wnt classifier signature) patients were analyzed using Fisher’s exact test and Kaplan-Meier survival curves were created in R/Bioconductor using the “survival” package with significant differences in risk calculated using the log-rank test. The details of the analyses and statistical methods are presented in the supplementary section as “Additional file 4”.
Gene expression was reproducible and concordant with clinical pathological subtype
Messenger RNA expression from FFPE samples were characterized using cDNA-mediated, Annealing, Selection, Extension, and Ligation (DASL)  for two populations of breast cancer patients: one from Quebec, and the other from Georgia. Quebec samples came from an unselected patient population from a community hospital in Canada whereas the Georgia cohort was preferentially selected for TNBC. The Quebec cohort was representative of the breast cancer population at the local hospital, whereas the Georgia cohort was preferentially selected for TNBC (Additional file 3: Figure S1H, K and L). Other inclusion criteria required specimens with more than 50% tumor content, matching pathology records for ER, PR, and HER2, and at least 200 ng of RNA. RNA was extracted and further quality controlled as previously described . Samples were run on two DASL panels; one a commercially available panel targeted at cancer related genes, the other a panel targeting genes relevant to breast cancer. RNA was run on both DASL panels with technical replicates when sufficient RNA was available, and had an average Pearson correlation coefficient (R2) of 0.96. Common genes (152) between the two DASL panels measured a correlation (R2) of 0.88 (Additional file 3: Figure S1A-F). Expression of ESR1 (transcript for ER), PGR (transcript for PR), and ERBB2 (transcript for HER2) corresponded well to pathology immunohistochemistry (IHC) records. DASL expression of ESR1 and PGR were greater in HR + as compared to TNBC and HER2+ subtypes, and likewise ERBB2 expression was higher in HER2+ samples compared to TNBC and HR + subtypes in all cohorts (Additional file 3: Figure S1G-L). These results are reproducible and consistent with pathology records.
Wnt signaling is upregulated in TNBC
Differentially expressed transcripts in the TNBC subtype of the Quebec and Georgia cohorts were annotated for canonical signal transduction pathways using the KEGG database . After Bonferroni correction for multiple hypotheses testing, the only significant result was the Wnt signaling pathway in the Quebec cohort (p = 0.048, Fisher’s exact test; (Additional file 3: Figure S2A-C). Further analysis of three published microarray studies, including a 99 patient cohort with pathology subtype from Memorial Sloan-Kettering Cancer Center (hereafter MSKCC-99) , a 186 patient cohort with intrinsic subtype from University of North Carolina Chapel Hill (hereafter UNCCH-186) , and a 159 patient cohort also with intrinsic subtype from University Hospital, Stockholm, Sweden (hereafter Stockholm-159) , the only significant results were Wnt (p = 0.001) and TGF-β (p = 0.035) in the Stockholm-159 basal-like BC subtype (Additional file 3: Figure S2D-F). Canonical pathway expression, assessed as the mean of normalized gene components, was measured for differential regulation between breast cancer subtypes. After correction for multiple hypotheses, the most notable pathways included Wnt signaling upregulated in the TNBC/basal-like BC subtypes of the Quebec, UNCCH-186, and Stockholm-159 cohorts; TGF-β in the TNBC subtypes of the Quebec and UNCCH-186 cohort; and both ErbB and VEGF signaling upregulated in the MSKCC-99 triple negative and UNCCH-186 basal-like subtypes (p < 0.05, Additional file 3: Figure S3). Of the 10 KEGG canonical signal transduction pathways investigated, Wnt was the most commonly overexpressed pathway in the TNBC subtypes indicating pathway perturbation in these types of breast cancers (Additional file 3: Figure S4). Based on these data we proceeded to examine experimental data sets to confirm these findings. Oncogenic signaling was investigated using gene sets derived from adenoviral vector transfected human mammary epithelial cells (HMECs) for β-catenin, E2F3, Myc, Ras, and Src genes . Transcripts uniquely induced and inhibited from each HMEC model were used to define oncogenic pathway expression in the same manner as that applied to the canonical pathways, except in this case both those genes expressed and inhibited were used to measure pathway expression (see supplementary methods – HMEC Oncogenic Pathways). Analysis of pathway regulation by subtype indicated Wnt/β-catenin as the most commonly upregulated pathway in TNBC subtypes (Additional file 3: Figure S5). Indeed, the Wnt pathway was upregulated relative to the HR + or luminal subtypes in each cohort examined (Additional file 3: Figure S6). Thus, these data indicate that both oncogenic and canonical Wnt signaling pathways are uniquely upregulated in the TNBC subtypes. Expression of Wnt/β-catenin components in patients classified as Wnt-compared to Wnt+ (Wnt classifier signature) in each of the 11 studies analyzed in the meta-analysis shows greater expression of tumors classified as Wnt+ (Wnt classifier signature) (Additional file 3: Figure S7). To further validate Wnt transcriptional activation in TNBC subtypes an independent experimental Wnt gene set was used to analyze breast cancer subtypes. This Wnt classifier signature (mentioned in figure as Wnt+) is composed of 81 genes identified by treating lung cancer cell lines with the Wnt3A ligand (hereafter LWS-81) . Despite a small overlap of the LWS-81 genes and probes available on the breast cancer DASL panel (n =15), this gene set clustered TNBC patients together in both Quebec and Georgia cohorts (Figure 1). Moreover, Wnt signaling defined by the LWS-81 genes was significantly upregulated in TNBC as compared HR + or luminal subtypes in the Quebec, Georgia, MSKCC-99, UNCCH-186, and Stockholm-159 studies (Figure 2). These data, in agreement with the β-catenin induced Wnt signaling from Bild et al.  (Additional file 3: Figure S5), uniformly indicate elevated oncogenic Wnt signaling in TNBC subtypes.
Functional involvement of Wnt signaling in metastasis–associated tumor cell phenotypes
TNBC is a highly metastatic disease. A transcriptionally active β-catenin (unphosphorylated on Serines 33 and 37 as well as Threonine 41)  is a direct functional readout of Wnt/β-catenin signaling. As an independent line of evidence for Wnt pathway activation in TNBC, and its functional association with the metastatic disease, the role of Wnt signaling was examined in vitro using TNBC cell lines. A functional blockade of Wnt/β-catenin pathway by either a pharmacological Wnt-antagonist, WntC59, or β-catenin SiRNA mediated genetic manipulation demonstrated that a functional perturbation of the pathway is causal to the metastasis- associated phenotypes including fibronectin-mediated and invasion in TNBC cells. Similarly, Wnt/β-catenin signaling attenuator, sulindac sulfide also inhibited migration of TNBC cells on fibronectin (data not shown). We observed that the treatment with a pharmacological Wnt-antagonist, WntC59 blocked fibronectin-mediated migration and invasion in MDA-MB231 cells (Figure 3A). This is consistent with our observation that the SiRNA-mediated decrease in levels of β-catenin protein at different time points following tranfection of β-catenin SiRNA in MDA-MB231 cells (Figure 3B) caused significant inhibition of fibronectin-mediated migration and invasion. Since F-actin organization is one of the key effectors of cell movement, we tested the effect of inhibition of Wnt-signals on the F-actin organization in TNBC cells. Out data show that sulindac sulfide substantially abrogated the cellular organization of F-actin (Z-section) on fibronectin in MDA-MB231 and HCC38 TNBC cells (Figure 3C i, ii and iii). Finally, we carried out a real time video microscopy of the cell movement following the inhibition of Wnt-signals in HCC38 TNBC cells (Additional files 5 and 6; two AVI files).
To address the issue of whether or not the modulation of Wnt-beta-catenin pathway causes the attenuation of its downstream signals, we tested the expression of active beta-catenin following beta-catenin SiRNA or sulindac sulfide treatment in MDA-MB-231 cells (Figure 3D left panel). Since active beta-catenin levels were decreased in both the conditions, we also conducted the experiment to study the transcriptional activity of beta-catenin. Our data show a significant decrease in the relative luciferase activity following beta-catenin SiRNA or sulindac sulfide treatment in MDA-MB-231 cells (Figure 3D right panel).
Wnt/β-catenin signaling is associated with metastatic disease
We have performed Kaplan-Meier survival curves of Wnt/β-catenin positive (Wnt+: Wnt classifier signature) and Wnt/β- catenin negative (Wnt-) patients with respect to overall survival, recurrence-free survival, metastasis-free survival, lung metastasis-free survival, brain metastasis-free survival, and lung metastasis-free survival in TNBC patients (Figure 4). To investigate the implications of Wnt signaling, we used a nearest shrunken centroid classifier  to stratify patients by Wnt/β-catenin transcriptional activity. This Wnt/β-catenin classifier was trained to identify β-catenin transformed HMECs as opposed to normal and other oncogenic (E2F3, Myc, Ras, and Src) models  and subsequently applied to a meta-analysis of 11 studies and 1,878 expression profiles from primary breast cancers [25–32, 37] including the Quebec and Georgia cohorts. Well known components of the Wnt signaling pathway is represented in Figure 5. Wnt/β-catenin positive (Wnt+) tumors accounted for 188 of the 1,878 patients. Patients were subsequently analyzed for intrinsic and pathology determined subtype, lymph node status, and grade, as well as metastasis-free (MFS), recurrence-free (RFS), and overall survival (OS). Complete pathology records for ER, PR, and HER2 were available for 310 patients, and of these 56 were categorized as Wnt+, 52 of which were TNBC. This strongly supports earlier observations that Wnt/β-catenin is preferentially activated in the TNBC subtype (p = 6.3 × 10-14, Fisher’s exact test). Analysis of 465 patients with intrinsic subtype found 53 of 71 Wnt + tumors corresponded to the basal-like subtype, likewise implicating elevated Wnt signaling in this subtype (p = 2.2 × 10-16). Data regarding spread of disease to lymph nodes were available for 1,202 patients where Wnt + patients composed 7.1% and 10.2% of lymph node negative and positive cases, respectively, marginally associating Wnt signaling with positive lymph node status (p = 0.042). In contrast, analysis of 912 patients with histological grades 1, 2, and 3 consisted of 1.1%, 5.7%, and 15.2% of Wnt + patients, correlating Wnt signaling with grade 3 carcinomas (p = 4.0 × 10-11, Additional file 3: Table S3). These data support our earlier observations of elevated Wnt signaling in TNBC subtypes, and associate Wnt signaling with high grade carcinomas. Survival analyses of breast cancer patients with elevated Wnt signaling distinguish these cancers as having greater metastatic potential and overall worse prognoses. Kaplan-Meier survival curves for OS, RFS, and MFS, including lung, brain, and bone specific metastases where analyzed for Wnt + and Wnt- patients stratified by pathology and intrinsic subtypes, grade, and lymph node status (Figure 6). Significantly increased risk was found for Wnt + patients with respect to OS, RFS, and MFS (p < 0.05, Figure 4). However, these differences were most significant in the metastatic setting (p = 9.6 × 10-7) and specifically in lung and brain but not bone metastases (Figure 6A vi). Stratifying patients by pathological determined subtype limited this analysis to a much smaller cohort of patients with both outcome and pathology records; however, Wnt + TNBC patients had greater risk of lung metastases (p = 0.0393, Figure 4F). Other significant differences include a worse prognosis for Grade 2 Wnt + patients with respect to OS as well as RFS (Figure 6D), and increased risk for lung and brain metastases for lymph node negative and positive patients, respectively (Figure 6E). These data cumulatively suggest that increased Wnt/β-catenin signaling is associated with metastatic pathways to the brain and lung.
Our study demonstrates two major findings that are consistent with current knowledge and advance our understanding of subset specific breast cancer etiology. First, using novel technologies to characterize FFPE materials we ascertained data consistent with clinical pathological subtype that also identified canonical and oncogenic Wnt signaling as an activated pathway in TNBC. These findings were confirmed using multiple published studies in diverse cohorts of patients across different microarray platforms, as well as independent experiments that identified Wnt induced targets which were consistently upregulated in TNBC subtypes. In the context of reports indicating associations of Wnt signaling with the basal-like subtype [38, 39], these data strongly suggest that the Wnt pathway is preferentially activated in TNBC subtypes, and may represent a possible therapeutic target in the treatment of these cancers. Reis-Filho’s group reported that β-catenin pathway activation in BC is associated with the TNBC phenotype but not with CTNNB1 mutation . We have observed (by immuno-fluorescence) a higher level of β-catenin (active) in the nucleus of MDA-MB231 TNBC cells as compared to non-TNBC (MCF7 and BT474) cell lines (data not shown).
Khramtsov et al., and others reported the association of Wnt signaling in TNBC with higher metastasis and poor prognosis [40, 41]. This observation can be explained by the fact that Wnt-β-catenin pathway plays a critical role in the regulation of metastasis-associated phenotypes in tumor cells. Breast cancer metastases are osteolytic in nature, and osteolytic bone lesions are formed due to tumor-induced bone resorption and destruction . Regulatory mechanisms underlying osteolytic metastasis to bone is a vicious cycle reflecting complex interplay of molecules which is propagated by four contributors: tumor cells, osteoblasts, osteoclasts and factors within bone matrix [43, 44]. Wnt-pathway has emerged as a crucial regulator of bone formation, and regeneration as Wnt signaling stimulates bone formation, and is also reported as a therapeutic target for bone diseases [45, 46]. Wnt signaling in osteoblasts regulates expression of the receptor activator of NFkappaB ligand, and inhibits osteoclastogenesis in vitro (while Dickkopf1, a secreted Wnt/beta-catenin antagonist, produced by breast tumor cells is an important mechanistic link between primary breast tumors, and secondary osteolytic bone metastases . In case of hormone-receptor positive tumors, it has been reported that in contrast to its role in breast cancer initiation, estrogen signaling has a protective effect in later stages, where estrogen receptor (ER) loss associates with aggressive metastatic disease . In the context of the above mentioned reports, it appears that Wnt-positivity in hormone-receptor positive patients (HR+/Wnt+) may have a negative regulatory influence on bone metastases as observed by the lower rate of bone metastatic events as compared to hormone receptors-positive Wnt-negative patients (HR+/Wnt-).
In our data, both TNBC and HER2+ have somewhat similar/identical average levels of wnt-expression, and expression of wnt/beta-catenin genes (Figure 2 and Additional file 3: Figure S5). This similarity in the expression pattern(s) between the TNBC and the HER2+ group can be due to an intrinsic heterogeneity within HER2-enriched/amplified subtype. Indeed, it has been shown that, when only HER2-amplified breast cancers are taken into account, approximately 50% are ER-positive . Hence the other 50% are ER-negative HER2-amplified breast cancers. Previous comparative genomic hybridization (CGH) studies have demonstrated that ER- negative disease differs significantly from ER-positive cancers in terms of the pattern, type, and complexity of genetic aberrations [50–53]. Furthermore, data at the European Society for Medical Oncology (ESMO) Vienna 2012 congress on the duration of adjuvant trastuzumab therapy hint at a difference between HER2 + /ER + and HER2 + /ER − disease, in keeping with the concept that HER2 + /luminal is biologically distinct from HER2 + /HER2-enriched disease, which is predominantly ER − [30, 54]. Sircoulomb et al. have shown that ER + and ER- ERBB2-amplified BCs are different and the WNT/b catenin signaling pathway was involved in ERERBB2-amplified BCs . Thus it is highly possible that ER-negative HER2-amplified tumors present within the HER2-amplified group in our study can influence the levels of wnt-expression and expression of wnt/β-catenin genes.
Wnt-C59 is a potent and selective Wnt signaling modulator with IC50 <0.11 nM in Wnt-Luc reporter assay for Wnt pathway inhibition, and with chemical/physical properties, suitable for in vitro/in vivo studies. Wnt-C59 prevents palmitylation of Wnt proteins by Porcupine, thereby blocking Wnt secretion and activity, similar to Wnt inhibitors IWP-2, IWP-3 and IWP-4. The observed inhibition of integrin-directed migration and invasion of MDA-MB231 cells following Wnt-C59 treatment in our results provides mechanistic explanation to our observation that, (1) Wnt signaling is upregulated in TNBC, and (2) Wnt/β-catenin signaling is associated with metastatic disease. Recently, Craig et al., have reported genome and transcriptome sequencing in prospective metastatic triple negative breast cancer . To further ascertain the functional significance of the pathway in metastatic disease, we also genetically manipulated the cellular levels of β-catenin, the functional readout of Wnt/β-catenin pathway (Figure 5). Together, our functional data demonstrate a direct involvement of Wnt-β-catenin pathway in the metastasis-associated phenotypes in tumor cell (Figure 3 and Additional file 5). We have studied the cause-effect relationship of Wnt-β-catenin pathway with metastasis in TNBC cell line models using three tools, genetic (beta-catenin SiRNA), pharmacological (Wnt-β-catenin pathway modulator; Wnt-C59), and functional (sulindac sulfide). Our data show that perturbation of the Wnt-β-catenin pathway abrogated metastasis-associated phenotypes in TNBC cells following attenuation of beta-catenin transcriptional activity, proving a direct mechanism based relationship between Wnt-β-catenin pathway and metastasis in TNBC. The results of the study may have implications for therapeutic target identification in future. The functional data would benefit from validation in other in vitro models.
Furthermore, other observations in this analysis include upregulation of Myc regulated genes (Additional file 3: Figure S5) that is consistent with recent reports identifying this Wnt transcriptional target  as upregulated in the basal-like subtype . Thus, these data indicate that a significant subset of TNBC is characterized by Wnt activation. Notable Wnt transcriptional targets upregulated in TNBC (Figure 5) included matrix metallopeptidase 7 (MMP7) , inteleukin 8 (IL8) , MYC , VEGF , frizzled 7 (FZD7) , survivin (BIRC5) , CD44 , MET [38, 64], peroxisome proliferator-activated receptor gamma (PPARD) , uPAR (PLAUR) , and snail (SNAI1)  (Figure 5). Furthermore, several Wnt antagonists were downregulated in the TNBC subtypes such as the androgen receptor (AR) [68, 69], FOXA1 , and MYB . These data highlight some of the Wnt components differentially regulated in TNBC.
In addition to the observed upregulation of Wnt signaling in TNBC, we also found association of Wnt signaling with metastatic disease. The Wnt/β-catenin classifier trained to identify oncogenic β-catenin signaling identified a disproportionate number of TNBC patients, supporting earlier observations of preferential Wnt activation in this subtype. Importantly, this classifier identified patients that were more likely to experience lung and brain metastases. These two metastatic routes have been associated with the basal-like subtype where Wnt signaling was noted as upregulated , and more recently, Wnt has been causally implicated in lung metastases . Our analyses further establish these findings and suggest that Wnt signaling confers a greater risk of lung metastases within the TNBC subtype (Figure 4F). There are certain limitations of the study. Sample size of HER2 enriched/amplified group in both Quebec and Georgia cohorts is lower than HR + and TNBC groups. We have restricted our study only to the 3 major clinical treatment categories of breast cancer, HR+, HER2+ and TNBC based on comprehensive gene expression profiling. However, clinically, each of three major categories of breast cancer is also a heterogeneous group by themselves. Several recent studies have described that even the relatively small class of breast tumors like TNBC can be further divided into five or six subclasses, each with its own molecular features, and unique sensitivity to therapeutic agents. A number of hypotheses have been proposed to explain the origin of inter-tumor heterogeneity in breast cancer, including subtype-specific tumor cell–of–origin and transforming events [72, 73]. We also have not stratified our study in the node-negative and node-positive settings.
Targeting the Wnt pathway has traditionally been difficult, but emerging modalities provide potential opportunities. Examples of small molecule inhibitors include sulindac, XAV939, ICG-001, and thiazolidinediones (TZDs) (Figure 5). Sulindac is a non-steroidal anti-inflammatory drug (NSAID) which also inhibits Wnt signaling by binding the PDZ domain of disheveled (DVL1) [74, 75] and like other NSAIDS, inhibits cyclooxygenase-2 (COX2), a gene recently implicated in breast cancer metastasis to the brain [25, 76]. XAV939 is a small molecule inhibitor that has recently been identified as targeting the poly-ADP ribose polymerase (PARP) gene tankyrase, which degrades Axin and allows β-catenin to avoid phospohorylation, subsequent poly-ubiquitination, and proteasomal degradation. We have observed that both sulindac and XAV939 blocks metastasis associated phenotypes (e.g. integrin-dependent migration, matrigel invasion, vascular mimicry) and clonogenic growth in TNBC cells lines (data not shown).
Our study contributes not only by identifying clinical markers associated with Wnt signaling such as histological grade 3 tumors and TNBC pathological subtype, but also increased risk for brain and lung metastases, thus recognizing Wnt signaling as a rational target in TNBC. The results of the study have implications for therapeutic target identification and the design of future clinical trials for this aggressive group of breast cancer. More genomic studies like this, however, are needed to create a genetic landscape of TNBC which will be utilized to differentiate “driver mutations” from “carrier mutations” and will guide therapeutics development. Individualized treatment will be possible only once we fully appreciate the biology of these genetic abnormalities.
Introduction to Additional file 3: Details of Figure S1: Assay reproducibility in the Quebec cohort measured by a distribution of Pearson R2 coefficients between RNA replicates on the A) human, and B) breast cancer (BC) DASL panels. C) Overlap of gene content on the human (dashed line) and BC (solid line) DASL panels, and D) replicated samples on both platforms. E) Distribution of inter-platform Pearson R2 correlations for samples run on both platforms, and F) a Venn diagram of differentially regulated genes found using the overlap of patients and genes on the human and BC DASL panels. Expression of ESR1, PGR, and ERBB2 by pathology determined subtype for HR + (white), HER2+ (grey), and TN (black) subtypes in the Quebec cohort on the G) human and I) BC DASL panels, as well as the K) Georgia cohort on the BC DASL panel correspond with expected clinical pathological subtype. Cohort sizes by pathology subtype for the Quebec cohort on the H) human and J) BC DASL panels, as well as the L) Georgia cohort on the BC DASL panel are depicted by pie charts. Details of Figure S2: Analysis of upregulated probes in TNBC as compared to other subtypes in context with KEGG  signal transduction pathways. Overrepresented probes in each pathway were compared to the number of pathway probes available on the specific platform and the total number of probes upregulated using Fisher’s exact test. A) Quebec cohort on the human cancer DASL panel (QC-HCP), B) Quebec on the BC DASL panel (QC-BCP), C) Georgia cohort on the BC DASL panel (GA-BCP), D) MSKCC-99 , E) UNCCH-186 , and F) Stockholm-159. P-values were Bonferroni corrected and significant pathways (p <0.05) are in bold. Details of Figure S3: Pathway expression of canonical KEGG signal transduction pathways  were measured as the normalized mean of the pathway components and subsequently used to calculate pathway perturbation between BC subtypes (see Additional file 4). Pathways that were differentially expressed relative to the TNBC subtype have significance lines and a corresponding p-value with Bonferroni’s correction for multiple hypothesis testing applied. Analyses included cohorts from Quebec profiled on the A) human (QCHCP) and B) BC DASL panels (QC-BCP), C) Georgia on the BC DASL panel (GA-BCP), D) MSKCC-99 , E) UNCCH-186  and F) Stockholm-159 . Details of Figure S4: Canonical Wnt expression of KEGG  pathway components. Pathway regulation was higher in TNBC as compared to other subtypes in A & B) Quebec (QC-HCP & QC-BCP), C) Georgia (GABCP), D) MSKCC-99 , E) UNCCH-186 , and F) Stockholm-159 cohorts. Significant p-values after Bonferroni’s correction are shown relative to the TNBC subtypes (see Additional file 4). Details of Figure S5: Pathway expression of experimentally derived oncogenic signaling pathways from Bild et al.  measured between BC subtypes (see Additional file 4). Pathways that were differentially expressed after Bonferroni’s correction are in bold. Analysis include cohorts from Quebec profiled on the A) human (QC-HCP) and B) BC DASL panels (QCBCP), C) Georgia on the BC DASL panel (GA-BCP), D) MSKCC-99 , E) UNCCH-186 , and F) Stockholm-159 . Hierarchical clustering of pathway expression depicts patterns in pathway regulation (rows) in context with BC subtype (columns). Details of Figure S6: Experimentally induced Wnt/β-catenin pathway expression from Bild et al. . Pathway regulation was assessed in the Quebec cohort on the A) human (QC-HCP) and B) BC DASL panels (QC-BCP), C) Georgia (GA-BCP), D) MSKCC-99 , E) UNCCH-186 , and F) Stockholm-159 cohorts. Significant p-values after Bonferroni’s correction are shown relative to the TNBC subtypes.
Human epidermal growth factor receptor 2
cDNA mediated annealing, selection, and ligation.
Perou CM, Sorlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, Pollack JR, Ross DT, Johnsen H, Akslen LA, et al: Molecular portraits of human breast tumours. Nature. 2000, 406 (6797): 747-752. 10.1038/35021093.
Bertucci F, Finetti P, Cervera N, Esterni B, Hermitte F, Viens P, Birnbaum D: How basal are triple-negative breast cancers?. Int J Cancer. 2008, 123 (1): 236-240. 10.1002/ijc.23518.
Hennessy BT, Gonzalez-Angulo AM, Stemke-Hale K, Gilcrease MZ, Krishnamurthy S, Lee JS, Fridlyand J, Sahin A, Agarwal R, Joy C, et al: Characterization of a naturally occurring breast cancer subset enriched in epithelial-to-mesenchymal transition and stem cell characteristics. Cancer Res. 2009, 69 (10): 4116-4124. 10.1158/0008-5472.CAN-08-3441.
Peddi PF, Ellis MJ, Ma C: Molecular basis of triple negative breast cancer and implications for therapy. Int J Breast Cancer. 2012, 2012: 217185-
Ben-Porath I, Thomson MW, Carey VJ, Ge R, Bell GW, Regev A, Weinberg RA: An embryonic stem cell-like gene expression signature in poorly differentiated aggressive human tumors. Nat Genet. 2008, 40 (5): 499-507. 10.1038/ng.127.
Creighton CJ, Li X, Landis M, Dixon JM, Neumeister VM, Sjolund A, Rimm DL, Wong H, Rodriguez A, Herschkowitz JI, et al: Residual breast cancers after conventional therapy display mesenchymal as well as tumor-initiating features. Proc Natl Acad Sci USA. 2009, 106 (33): 13820-13825. 10.1073/pnas.0905718106.
Shah SP, Roth A, Goya R, Oloumi A, Ha G, Zhao Y, Turashvili G, Ding J, Tse K, Haffari G, et al: The clonal and mutational evolution spectrum of primary triple-negative breast cancers. Nature. 2012, 486 (7403): 395-399.
Ding L, Ellis MJ, Li S, Larson DE, Chen K, Wallis JW, Harris CC, McLellan MD, Fulton RS, Fulton LL, et al: Genome remodelling in a basal-like breast cancer metastasis and xenograft. Nature. 2010, 464 (7291): 999-1005. 10.1038/nature08989.
Carey LA: Directed therapy of subtypes of triple-negative breast cancer. Oncologist. 2011, 16 (Suppl 1): 71-78.
Moulder SL: Does the PI3K pathway play a role in basal breast cancer?. Clin Breast Cancer. 2010, 10 (Suppl 3): S66-S71.
Hoeflich KP, O’Brien C, Boyd Z, Cavet G, Guerrero S, Jung K, Januario T, Savage H, Punnoose E, Truong T, et al: In vivo antitumor activity of MEK and phosphatidylinositol 3-kinase inhibitors in basal-like breast cancer models. Clin Cancer Res. 2009, 15 (14): 4649-4664. 10.1158/1078-0432.CCR-09-0317.
Geyer FC, Lacroix-Triki M, Savage K, Arnedos M, Lambros MB, MacKay A, Natrajan R, Reis-Filho JS: Beta-catenin pathway activation in breast cancer is associated with triple-negative phenotype but not with CTNNB1 mutation. Mod Pathol. 2011, 24 (2): 209-231. 10.1038/modpathol.2010.205.
Abramovitz M, Ordanic-Kodani M, Wang Y, Li Z, Catzavelos C, Bouzyk M, Sledge GW, Moreno CS, Leyland-Jones B: Optimization of RNA extraction from FFPE tissues for expression profiling in the DASL assay. Biotechniques. 2008, 44 (3): 417-423. 10.2144/000112703.
Abramovitz M, Barwick BG, Willis S, Young B, Catzavelos C, Li Z, Kodani M, Tang W, Bouzyk M, Moreno CS, et al: Molecular characterisation of formalin-fixed paraffin-embedded (FFPE) breast tumour specimens using a custom 512-gene breast cancer bead array-based platform. Br J Cancer. 2011, 105 (10): 1574-1581. 10.1038/bjc.2011.355.
Dey N, De PK, Wang M, Zhang H, Dobrota EA, Robertson KA, Durden DL: CSK controls retinoic acid receptor (RAR) signaling: a RAR-c-SRC signaling axis is required for neuritogenic differentiation. Mol Cell Biol. 2007, 27 (11): 4179-4197. 10.1128/MCB.01352-06.
Dey N, Crosswell HE, De P, Parsons R, Peng Q, Su JD, Durden DL: The protein phosphatase activity of PTEN regulates SRC family kinases and controls glioma migration. Cancer Res. 2008, 68 (6): 1862-1871. 10.1158/0008-5472.CAN-07-1182.
Dey N, Young B, Abramovitz M, Bouzyk M, Barwick B, De P, Leyland-Jones B: Differential activation of Wnt-b-catenin pathway in triple negative breast cancer increases MMP7 in a PTEN dependent manner. PLoS One. 2013, In Press
Zeng G, Apte U, Cieply B, Singh S, Monga SP: siRNA-mediated beta-catenin knockdown in human hepatoma cells results in decreased growth and survival. Neoplasia. 2007, 9 (11): 951-959. 10.1593/neo.07469.
Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001, 98 (9): 5116-5121. 10.1073/pnas.091062498.
Reimers M, Carey VJ: Bioconductor: an open source framework for bioinformatics and computational biology. Methods Enzymol. 2006, 411: 119-134.
Kanehisa M: The KEGG database. Novartis Found Symp. 2002, 247: 91-101. discussion 101–103, 119–128, 244–152
Bild AH, Yao G, Chang JT, Wang Q, Potti A, Chasse D, Joshi MB, Harpole D, Lancaster JM, Berchuck A, et al: Oncogenic pathway signatures in human cancers as a guide to targeted therapies. Nature. 2006, 439 (7074): 353-357. 10.1038/nature04296.
Edgar R, Domrachev M, Lash AE: Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30 (1): 207-210. 10.1093/nar/30.1.207.
Nguyen DX, Chiang AC, Zhang XH, Kim JY, Kris MG, Ladanyi M, Gerald WL, Massague J: WNT/TCF signaling through LEF1 and HOXB9 mediates lung adenocarcinoma metastasis. Cell. 2009, 138 (1): 51-62. 10.1016/j.cell.2009.04.030.
Bos PD, Zhang XH, Nadal C, Shu W, Gomis RR, Nguyen DX, Minn AJ, van de Vijver MJ, Gerald WL, Foekens JA, et al: Genes that mediate breast cancer metastasis to the brain. Nature. 2009, 459 (7249): 1005-1009. 10.1038/nature08021.
Herschkowitz JI, Simin K, Weigman VJ, Mikaelian I, Usary J, Hu Z, Rasmussen KE, Jones LP, Assefnia S, Chandrasekharan S, et al: Identification of conserved gene expression features between murine mammary carcinoma models and human breast tumors. Genome Biol. 2007, 8 (5): R76-10.1186/gb-2007-8-5-r76.
Hoadley KA, Weigman VJ, Fan C, Sawyer LR, He X, Troester MA, Sartor CI, Rieger-House T, Bernard PS, Carey LA, et al: EGFR associated expression profiles vary with breast tumor subtype. BMC Genomics. 2007, 8: 258-10.1186/1471-2164-8-258.
Loi S, Haibe-Kains B, Desmedt C, Lallemand F, Tutt AM, Gillet C, Ellis P, Harris A, Bergh J, Foekens JA, et al: Definition of clinically distinct molecular subtypes in estrogen receptor-positive breast carcinomas through genomic grade. J Clin Oncol. 2007, 25 (10): 1239-1246. 10.1200/JCO.2006.07.1522.
Minn AJ, Gupta GP, Siegel PM, Bos PD, Shu W, Giri DD, Viale A, Olshen AB, Gerald WL, Massague J: Genes that mediate breast cancer metastasis to lung. Nature. 2005, 436 (7050): 518-524. 10.1038/nature03799.
Parker JS, Mullins M, Cheang MC, Leung S, Voduc D, Vickery T, Davies S, Fauron C, He X, Hu Z, et al: Supervised risk predictor of breast cancer based on intrinsic subtypes. J Clin Oncol. 2009, 27 (8): 1160-1167. 10.1200/JCO.2008.18.1370.
Pawitan Y, Bjohle J, Amler L, Borg AL, Egyhazi S, Hall P, Han X, Holmberg L, Huang F, Klaar S, et al: Gene expression profiling spares early breast cancer patients from adjuvant therapy: derived and validated in two population-based cohorts. Breast Cancer Res. 2005, 7 (6): R953-R964. 10.1186/bcr1325.
Wang Y, Klijn JG, Zhang Y, Sieuwerts AM, Look MP, Yang F, Talantov D, Timmermans M, Meijer-van Gelder ME, Yu J, et al: Gene-expression profiles to predict distant metastasis of lymph-node-negative primary breast cancer. Lancet. 2005, 365 (9460): 671-679. 10.1016/S0140-6736(05)17947-1.
Tibshirani R, Hastie T, Narasimhan B, Chu G: Diagnosis of multiple cancer types by shrunken centroids of gene expression. Proc Natl Acad Sci USA. 2002, 99 (10): 6567-6572. 10.1073/pnas.082099299.
Fan JB, Yeakley JM, Bibikova M, Chudin E, Wickham E, Chen J, Doucet D, Rigault P, Zhang B, Shen R, et al: A versatile assay for high-throughput gene expression profiling on universal array matrices. Genome Res. 2004, 14 (5): 878-885. 10.1101/gr.2167504.
Minn AJ, Gupta GP, Padua D, Bos P, Nguyen DX, Nuyten D, Kreike B, Zhang Y, Wang Y, Ishwaran H, et al: Lung metastasis genes couple breast tumor size and metastatic spread. Proc Natl Acad Sci USA. 2007, 104 (16): 6740-6745. 10.1073/pnas.0701138104.
Wang PS, Chou FS, Bloomston M, Vonau MS, Saji M, Espinosa A, Pinzone JJ: Thiazolidinediones downregulate Wnt/beta-catenin signaling via multiple mechanisms in breast cancer cells. J Surg Res. 2009, 153 (2): 210-216. 10.1016/j.jss.2008.05.032.
Ivshina AV, George J, Senko O, Mow B, Putti TC, Smeds J, Lindahl T, Pawitan Y, Hall P, Nordgren H, et al: Genetic reclassification of histologic grade delineates new clinical subtypes of breast cancer. Cancer Res. 2006, 66 (21): 10292-10301. 10.1158/0008-5472.CAN-05-4414.
Ponzo MG, Lesurf R, Petkiewicz S, O’Malley FP, Pinnaduwage D, Andrulis IL, Bull SB, Chughtai N, Zuo D, Souleimanova M, et al: Met induces mammary tumors with diverse histologies and is associated with poor outcome and human basal breast cancer. Proc Natl Acad Sci USA. 2009, 106 (31): 12903-12908. 10.1073/pnas.0810402106.
Smid M, Wang Y, Zhang Y, Sieuwerts AM, Yu J, Klijn JG, Foekens JA, Martens JW: Subtypes of breast cancer show preferential site of relapse. Cancer Res. 2008, 68 (9): 3108-3114. 10.1158/0008-5472.CAN-07-5644.
Khramtsov AI, Khramtsova GF, Tretiakova M, Huo D, Olopade OI, Goss KH: Wnt/beta-catenin pathway activation is enriched in basal-like breast cancers and predicts poor outcome. Am J Pathol. 2010, 176 (6): 2911-2920. 10.2353/ajpath.2010.091125.
DiMeo TA, Anderson K, Phadke P, Fan C, Perou CM, Naber S, Kuperwasser C: A novel lung metastasis signature links Wnt signaling with cancer cell self-renewal and epithelial-mesenchymal transition in basal-like breast cancer. Cancer Res. 2009, 69 (13): 5364-5373. 10.1158/0008-5472.CAN-08-4135.
Bu G, Lu W, Liu CC, Selander K, Yoneda T, Hall C, Keller ET, Li Y: Breast cancer-derived Dickkopf1 inhibits osteoblast differentiation and osteoprotegerin expression: implication for breast cancer osteolytic bone metastases. Int J Cancer. 2008, 123 (5): 1034-1042. 10.1002/ijc.23625.
Vashisht S, Bagler G: An approach for the identification of targets specific to bone metastasis using cancer genes interactome and gene ontology analysis. PLoS One. 2012, 7 (11): e49401-10.1371/journal.pone.0049401.
Kozlow W, Guise TA: Breast cancer metastasis to bone: mechanisms of osteolysis and implications for therapy. J Mammary Gland Biol Neoplasia. 2005, 10 (2): 169-180. 10.1007/s10911-005-5399-8.
Hoeppner LH, Secreto FJ, Westendorf JJ: Wnt signaling as a therapeutic target for bone diseases. Expert Opin Ther Targets. 2009, 13 (4): 485-496. 10.1517/14728220902841961.
Krishnan V, Bryant HU, Macdougald OA: Regulation of bone mass by Wnt signaling. J Clin Invest. 2006, 116 (5): 1202-1209. 10.1172/JCI28551.
Spencer GJ, Utting JC, Etheridge SL, Arnett TR, Genever PG: Wnt signalling in osteoblasts regulates expression of the receptor activator of NFkappaB ligand and inhibits osteoclastogenesis in vitro. Journal of cell science. 2006, 119 (Pt 7): 1283-1296.
Chimge NO, Baniwal SK, Little GH, Chen YB, Kahn M, Tripathy D, Borok Z, Frenkel B: Regulation of breast cancer metastasis by Runx2 and estrogen signaling: the role of SNAI2. Breast Cancer Res. 2011, 13 (6): R127-10.1186/bcr3073.
Marchio C, Natrajan R, Shiu KK, Lambros MB, Rodriguez-Pinilla SM, Tan DS, Lord CJ, Hungermann D, Fenwick K, Tamber N, et al: The genomic profile of HER2-amplified breast cancers: the influence of ER status. J Pathol. 2008, 216 (4): 399-407. 10.1002/path.2423.
Loo LW, Grove DI, Williams EM, Neal CL, Cousens LA, Schubert EL, Holcomb IN, Massa HF, Glogovac J, Li CI, et al: Array comparative genomic hybridization analysis of genomic alterations in breast cancer subtypes. Cancer Res. 2004, 64 (23): 8541-8549. 10.1158/0008-5472.CAN-04-1992.
Simpson PT, Reis-Filho JS, Gale T, Lakhani SR: Molecular evolution of breast cancer. J Pathol. 2005, 205 (2): 248-254. 10.1002/path.1691.
Melchor L, Honrado E, Huang J, Alvarez S, Naylor TL, Garcia MJ, Osorio A, Blesa D, Stratton MR, Weber BL, et al: Estrogen receptor status could modulate the genomic pattern in familial and sporadic breast cancer. Clin Cancer Res. 2007, 13 (24): 7305-7313. 10.1158/1078-0432.CCR-07-0711.
Korsching E, Packeisen J, Helms MW, Kersting C, Voss R, van Diest PJ, Brandt B, van der Wall E, Boecker W, Burger H: Deciphering a subgroup of breast carcinomas with putative progression of grade during carcinogenesis revealed by comparative genomic hybridisation (CGH) and immunohistochemistry. Br J Cancer. 2004, 90 (7): 1422-1428. 10.1038/sj.bjc.6601658.
Ellis MJ, Perou CM: The genomic landscape of breast cancer as a therapeutic roadmap. Cancer discovery. 2013, 3 (1): 27-34. 10.1158/2159-8290.CD-12-0462.
Sircoulomb F, Bekhouche I, Finetti P, Adelaide J, Ben Hamida A, Bonansea J, Raynaud S, Innocenti C, Charafe-Jauffret E, Tarpin C, et al: Genome profiling of ERBB2-amplified breast cancers. BMC cancer. 2010, 10: 539-10.1186/1471-2407-10-539.
Craig DW, O’Shaughnessy JA, Kiefer JA, Aldrich J, Sinari S, Moses TM, Wong S, Dinh J, Christoforides A, Blum JL, et al: Genome and transcriptome sequencing in prospective metastatic triple-negative breast cancer uncovers therapeutic vulnerabilities. Mol Cancer Ther. 2013, 12 (1): 104-116. 10.1158/1535-7163.MCT-12-0781.
He TC, Sparks AB, Rago C, Hermeking H, Zawel L, da Costa LT, Morin PJ, Vogelstein B, Kinzler KW: Identification of c-MYC as a target of the APC pathway. Science. 1998, 281 (5382): 1509-1512.
Chandriani S, Frengen E, Cowling VH, Pendergrass SA, Perou CM, Whitfield ML, Cole MD: A core MYC gene expression signature is prominent in basal-like breast cancer but only partially overlaps the core serum response. PLoS One. 2009, 4 (8): e6693-10.1371/journal.pone.0006693.
Brabletz T, Jung A, Dag S, Hlubek F, Kirchner T: Beta-catenin regulates the expression of the matrix metalloproteinase-7 in human colorectal cancer. Am J Pathol. 1999, 155 (4): 1033-1038. 10.1016/S0002-9440(10)65204-2.
Masckauchan TN, Shawber CJ, Funahashi Y, Li CM, Kitajewski J: Wnt/beta-catenin signaling induces proliferation, survival and interleukin-8 in human endothelial cells. Angiogenesis. 2005, 8 (1): 43-51. 10.1007/s10456-005-5612-9.
Zhang T, Otevrel T, Gao Z, Ehrlich SM, Fields JZ, Boman BM: Evidence that APC regulates survivin expression: a possible mechanism contributing to the stem cell origin of colon cancer. Cancer Res. 2001, 61 (24): 8664-8667.
Willert J, Epping M, Pollack JR, Brown PO, Nusse R: A transcriptional response to Wnt protein in human embryonic carcinoma cells. BMC Dev Biol. 2002, 2: 8-10.1186/1471-213X-2-8.
Wielenga VJ, Smits R, Korinek V, Smit L, Kielman M, Fodde R, Clevers H, Pals ST: Expression of CD44 in Apc and Tcf mutant mice implies regulation by the WNT pathway. Am J Pathol. 1999, 154 (2): 515-523. 10.1016/S0002-9440(10)65297-2.
Boon EM, van der Neut R, van de Wetering M, Clevers H, Pals ST: Wnt signaling regulates expression of the receptor tyrosine kinase met in colorectal cancer. Cancer Res. 2002, 62 (18): 5126-5128.
He TC, Chan TA, Vogelstein B, Kinzler KW: PPARdelta is an APC-regulated target of nonsteroidal anti-inflammatory drugs. Cell. 1999, 99 (3): 335-345. 10.1016/S0092-8674(00)81664-5.
Mann B, Gelos M, Siedow A, Hanski ML, Gratchev A, Ilyas M, Bodmer WF, Moyer MP, Riecken EO, Buhr HJ, et al: Target genes of beta-catenin-T cell-factor/lymphoid-enhancer-factor signaling in human colorectal carcinomas. Proc Natl Acad Sci USA. 1999, 96 (4): 1603-1608. 10.1073/pnas.96.4.1603.
ten Berge D, Koole W, Fuerer C, Fish M, Eroglu E, Nusse R: Wnt signaling mediates self-organization and axis formation in embryoid bodies. Cell Stem Cell. 2008, 3 (5): 508-518. 10.1016/j.stem.2008.09.013.
Aarnisalo P, Palvimo JJ, Janne OA: CREB-binding protein in androgen receptor-mediated signaling. Proc Natl Acad Sci USA. 1998, 95 (5): 2122-2127. 10.1073/pnas.95.5.2122.
Mulholland DJ, Read JT, Rennie PS, Cox ME, Nelson CC: Functional localization and competition between the androgen receptor and T-cell factor for nuclear beta-catenin: a means for inhibition of the Tcf signaling axis. Oncogene. 2003, 22 (36): 5602-5613. 10.1038/sj.onc.1206802.
Gerin I, Bommer GT, Lidell ME, Cederberg A, Enerback S, Macdougald OA: On the role of FOX transcription factors in adipocyte differentiation and insulin-stimulated glucose uptake. J Biol Chem. 2009, 284 (16): 10755-10763. 10.1074/jbc.M809115200.
Kurahashi T, Nomura T, Kanei-Ishii C, Shinkai Y, Ishii S: The Wnt-NLK signaling pathway inhibits A-Myb activity by inhibiting the association with coactivator CBP and methylating histone H3. Mol Biol Cell. 2005, 16 (10): 4705-4713. 10.1091/mbc.E05-05-0470.
Polyak K: Heterogeneity in breast cancer. J Clin Invest. 2011, 121 (10): 3786-3788. 10.1172/JCI60534.
Bertos NR, Park M: Breast cancer - one term, many entities?. J Clin Invest. 2011, 121 (10): 3789-3796. 10.1172/JCI57100.
Han A, Song Z, Tong C, Hu D, Bi X, Augenlicht LH, Yang W: Sulindac suppresses beta-catenin expression in human cancer cells. Eur J Pharmacol. 2008, 583 (1): 26-31. 10.1016/j.ejphar.2007.12.034.
Lee HJ, Wang NX, Shi DL, Zheng JJ: Sulindac inhibits canonical Wnt signaling by blocking the PDZ domain of the protein dishevelled. Angew Chem Int Ed Engl. 2009, 48 (35): 6448-6452. 10.1002/anie.200902981.
Lu W, Tinsley HN, Keeton A, Qu Z, Piazza GA, Li Y: Suppression of Wnt/beta-catenin signaling inhibits prostate cancer cell proliferation. Eur J Pharmacol. 2009, 602 (1): 8-14. 10.1016/j.ejphar.2008.10.053.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/13/537/prepub
The authors would like to acknowledge the Emory Biomarker Service Center for microarray analysis, and the Cell Imaging and Microscopy Core of Winship Cancer Institute of Emory University. We thank Dr. Yajun Mei of Georgia Institute of Technology for statistical advice. We thank Dr. Kevin Ward for providing SEER registry information regarding pathological subtype for the Georgia cohort. We further acknowledge funding from the Georgia Research Alliance and Winship Cancer Institute. CSM was supported by NIH R01 CA106826. This work was supported by the Department of Defense Center of Excellence: “Center of Excellence for Individualization of Therapy for Breast Cancer” grant # W81XWH-04-1-0468.
The authors declare that they have no competing interests.
ND was involved in designing / conducting the cell-line based experiments, and actively contributed to writing the manuscript. BB is the key contributor in the data analyses, and one of the key authors on writing the manuscript together with ND. CM helped design the DASL component of the study, performed the DASL analysis, and contributed to writing the manuscript. MO prepared RNA. ZC performed statistical analysis. PD contributed in designing experiments and contributed to writing the manuscript. WT prepared RNA and performed some DASL experiments. CC conducted pathological analysis. KK was involved in coordinating the study. GS conceived of the study, and participated in its design. MA participated in the design of the study and helped to draft the manuscript. MB helped conceive and design the molecular component of the study and contributed to writing the manuscript. BLJ participated in the overall conception and the design of the study. All authors read and approved the final manuscript.
Nandini Dey, Benjamin G Barwick contributed equally to this work.
Electronic supplementary material
Additional file 3: Figure S1-Table S3: Assay reproducibility in the Quebec cohort measured by a distribution of Pearson R2 coefficients between RNA replicates. Figure S2: Analysis of upregulated probes in TNBC as compared to other subtypes in context with KEGG  signal transduction pathways. P-values were Bonferroni corrected and significant pathways (p <0.05) are in bold. Figure S3: Pathway expression of canonical KEGG signal transduction pathways  were measured as the normalized mean of the pathway components and subsequently used to calculate pathway perturbation between BC subtypes (see Additional file 4). Figure S4: Canonical Wnt expression of KEGG  pathway components. Significant p-values after Bonferroni’s correction are shown relative to the TNBC subtypes (see Additional file 4). Figure S5: Pathway expression of experimentally derived oncogenic signaling pathways from Bild et al.  measured between BC subtypes (see Additional file 4). Hierarchical clustering of pathway expression depicts patterns in pathway regulation (rows) in context with BC subtype (columns). Figure S6: Experimentally induced Wnt/β-catenin pathway expression from Bild et al. . Significant p-values after Bonferroni’s correction are shown relative to the TNBC subtypes. Figure S7: Expression of Wnt/β-catenin components in patients classified as Wnt compared to Wnt + (Wnt classifier signature) in each of the 11 studies analyzed in the meta-analysis shows greater expression of tumors classified as Wnt + (Wnt classifier signature). Table S3: Table identifying the number of patients from each cohort in a meta-analysis of 11 studies and 1,878 patients with pathological or intrinsic determined subtype, grade, and lymph node status. Each category is broken down by total number of Wnt + (Wnt classifier signature) and Wnt- patients identified by the Wnt/β-catenin classifier and then analyzed for overrepresentation using Fisher’s exact test. (PDF 149 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.