Breast cancer is the second leading cause of cancer-related death in women in the United States. During the advanced stages of disease, many breast cancer patients suffer from bone metastasis. These metastases are predominantly osteolytic and develop when tumor cells interact with bone. In vivo models that mimic the breast cancer-specific osteolytic bone microenvironment are limited. Previously, we developed a mouse model of tumor-bone interaction in which three mouse breast cancer cell lines were implanted onto the calvaria. Analysis of tumors from this model revealed that they exhibited strong bone resorption, induction of osteoclasts and intracranial penetration at the tumor bone (TB)-interface.
Methods
In this study, we identified and used a TB microenvironment-specific gene expression signature from this model to extend our understanding of the metastatic bone microenvironment in human disease and to predict potential therapeutic targets.
Results
We identified a TB signature consisting of 934 genes that were commonly (among our 3 cell lines) and specifically (as compared to tumor-alone area within the bone microenvironment) up- and down-regulated >2-fold at the TB interface in our mouse osteolytic model. By comparing the TB signature with gene expression profiles from human breast metastases and an in vitro osteoclast model, we demonstrate that our model mimics both the human breast cancer bone microenvironment and osteoclastogenesis. Furthermore, we observed enrichment in various signaling pathways specific to the TB interface; that is, TGF-β and myeloid self-renewal pathways were activated and the Wnt pathway was inactivated. Lastly, we used the TB-signature to predict cyclopenthiazide as a potential inhibitor of the TB interface.
Conclusion
Our mouse breast cancer model morphologically and genetically resembles the osteoclastic bone microenvironment observed in human disease. Characterization of the gene expression signature specific to the TB interface in our model revealed signaling mechanisms operative in human breast cancer metastases and predicted a therapeutic inhibitor of cancer-mediated osteolysis.
Bone is one of the most common sites for metastasis in human breast cancer. Bone metastasis results in cancer-related pain, pathological fracture, hypercalcemia, neurological defects, and immobility; all of which increase the risk of mortality and decrease the quality of life for breast cancer patients [1–4]. While a number of strategies exist to treat breast cancer bone metastases (e.g., surgery, radiation and/or chemotherapy), none are curative. Furthermore, these treatment methods have limited efficacy due in part to the fact that they do not effectively target the interaction between tumor cells and bone [5]. Even though the bisphosphonate class of drugs (which target the tumor-bone interface) have been shown to improve the quality of life and disease-free survival in some patients, more therapeutic targets and agents are desirable [6].
Within the osteolytic lesions of bone metastases, tumor cells interact with osteoclasts (bone resorbing cells) and osteoblasts (bone forming cells), thereby inhibiting normal bone development and ultimately leading to bone destruction [1–4]. As for osteoclasts, their interaction with tumor cells is reciprocal: tumor cells produce factors (e.g., parathyroid hormone-related peptide; interleukin-6; tumor necrosis factor; and macrophage colony stimulating factor, M-CSF) that directly or indirectly induce the formation of osteoclasts, and activated osteoclasts produce factors (e.g., transforming growth factor, TGF-β; insulin growth factor, IGF; and bone morphogenetic proteins, BMPs) that stimulate tumor growth and bone destruction [1]. Despite a general comprehension of this process, we are still far from a complete mechanistic understanding and lack well defined targets for therapeutic intervention.
Several animal models have been developed to study the mechanisms governing cancer-mediated osteolysis. However, there is no single animal model that ideally replicates the entire metastatic process from primary breast tumor to bone metastasis. Nevertheless, several models that represent various aspects of bone metastasis have been used successfully to study specific features of the disease. For example, Arguello, et al. developed a model in which melanoma cells injected into the left ventricle of the heart ultimately form bone metastases [7]. This model was later used to study various mechanisms behind breast cancer-specific osteoclast formation and bone metastasis [8–10]. Our group has also developed a rat model to study bone metastatic microenvironment in which prostate tumors were directly transplanted onto the calvariae of syngeneic animals. These tumors exhibited pathological osteoblastic and osteoclastic changes [11]. More recently, we used this approach with mouse breast cancer cell lines and found that the tumor cells induce osteolytic changes in the bone microenvironment [12–15]. With this model, we found that cathepsin G cleaves the receptor activator of nuclear factor-B ligand (Rankl) leading to enhanced activation of osteoclasts in the breast cancer bone microenvironment [15]. Furthermore, we also demonstrated the importance of TGF-β signaling and osteoclast activation in the breast cancer bone microenvironment [12, 14]. While this series of observations has furthered our understanding of the mechanisms underlying osteolysis, their relevance to human breast cancer remained unknown.
To address this question, we reanalyzed gene expression profiles generated from our previous studies using the syngeneic mouse model of breast cancer specific osteolysis that was developed by implanting 3 different cell lines - 4T1 (highly metastatic), Cl66 (moderately metastatic) and Cl66-M2 (low metastatic) - onto the calvariae bone of BALB/C mice [12–15]. The gene expression profiles were generated from microdissected tumors in which the tumor-bone (TB) interface and the tumor alone (TA) area were isolated independently. Then we identified a TB signature involved in bone destruction by comparing the gene expression profiles of the TA area and TB interface from the dissected tumors. Lastly, using our TB signature, open-access gene expression data, and pathway analytics, we demonstrated that our model mimics human disease and predicted key pathways and a potential therapeutic agent for breast cancer osteolysis.
Methods
Mouse osteolytic model and microarray
Mouse breast cancer cell lines - 4T1, Cl66 (66Cl4) and Cl66M2 - with different metastatic potential [16–18] were maintained in culture and were implanted under the dorsal skin flap onto the calvaria of female BALB/c mice, as described [12]. Mice were euthanized and necropsied to examine osteolytic lesions at 4 weeks post implantation. The tissues for histological examination were prepared as described [12]. All studies were carried out in accordance with the Institutional Animal Use and Care Committee (IAUCC) of the University of the Nebraska Medical Center. Calcified frozen tissues were serially sectioned into 10 μm slices and then microdissected to separate the TB interface from the TA area. RNA isolation and gene expression profiling of the TB interface and TA area were performed using Affymetrix GeneChip® Mouse Genome 430A 2.0 Array, as described [14].
Analysis of gene arrays and public microarray datasets
The CEL files for all the samples from Affymetrix GeneChip® were processed and MAS 5.0-normalized using the SimpleAffy [19] program and robust multiarray (RMA)-normalized using BRB Array tools [20]. The log2 MAS 5.0-normalized data was used for subsequent analyses. Fold-change at the TB interface with respect to the TA area for tissues, standard deviation (SD) across TA samples, and median-centered analysis within the TA area were calculated for each of the cell lines (4T1, Cl66 and Cl66-M2) to identify genes up- and down-regulated in the respective samples. The genes were ranked from highest to lowest expression based on the values from fold-change or median-centered analysis.
The following publicly available Affymetrix microarray data were obtained from Gene Expression Omnibus (GEO) [21]: GSE13563 for normal bone from mouse calvaria (n = 2), mandible (n = 2) and ulna (n = 2); GSE14017 (n = 29) and GSE14018 (n = 36) for metastases from breast cancer [22]; GSE11259 (n = 9) for 4T1 primary tumor data [23]; and GSE17563 (n = 3) for osteoclast precursors treated with human RANKL at different time points [24]. All the GEO data were processed and normalized as described above. Affymetrix microarray data for breast tumors (n = 118) [25] and cancer cell lines (n = 54) [26] were also compared with the TA area gene expression profile.
The NearestTemplatePrediction algorithm (NTP) [27, 28] was used to predict the class of a given sample with statistical significance (false discovery rate, FDR < 0.2) using a predefined set of markers that are specific to multiple classes (TB interface vs. TA area). Microarray data from different studies and platforms were sample- and gene-normalized and then pooled using the Distance Weighted Discrimination (DWD) algorithm, as described [29, 30]. The significance of expression between the mouse model and human bone metastases was estimated using SubMap [31]. Hierarchical clustering of genes and samples were performed using the Cluster software [32]. Visualization was performed with TreeView [32] and Hierarchical Clustering Viewer from GenePattern [33].
Gene ontology (GO) and pathway analysis
The association of gene signature with known pathways was determined using gene ontology (GO) [34], pathways from Kyoto Encyclopedia of Genes and Genomes (KEGG) [35], and Broad Institute based Molecular Signature Databases (MSigDB) [36]. The enrichment analysis was performed using the TB signature and the GlobalTest package (Version 4.2.0) [37, 38].
Connectivity Map analysis
Gene symbols were mapped to HG-U133A array probes. They were then used to query the Connectivity Map database [39].
Results
The TA area resembles the primary tumor
Previously, we transplanted three breast cancer cell lines - 4T1 (highly metastatic), Cl66 (moderately metastatic) and Cl66-M2 (low metastatic) - onto the calvarial bone of BALB/C mice [12, 14]. Irrespective of the cell lines used, histochemical analysis of these tumors demonstrated that they exhibited tumor-induced osteolysis and osteoclast activation similar to that observed in breast cancer bone metastasis [12, 14]. Metastatic lesions from the osteolytic tumors were microdissected into two cohorts - TB interface and TA area - and gene expression profile analyses were performed [12, 14]. Herein, we reanalyzed these gene expression data sets in search of a breast cancer osteolysis-specific gene signature.
Our reanalysis illustrates that there is little similarity in gene expression in the TA area samples among the three cell lines (Figure 1A). This is altogether not too surprising given that these cell lines were originally derived from different mouse tumors [16–18]. Consistently, the sublines Cl66 and Cl66-M2 [16–18], share the most similarity in gene expression (Figure 1A).
Figure 1
The TA area resembles matched primary tumors, human breast tumors and cancer cell lines. A. Hierarchical clustering of median centered gene expression from the TA area of the three cell lines (4T1, Cl66 and Cl66M2). B. The gene expression profile of primary syngeneic tumors derived using the 4T1, CI66 and 67NR mouse cell lines [23] resembled the TA area in our osteolytic breast cancer model as predicted by NTP algorithm (FDR < 0.2). Primary tumor samples and prediction significance (FDR < 0.2) are color coded and displayed at the top of the hierarchical cluster. C. The gene expression profile of primary human breast tumors [25] and D. human breast cancer cell lines resembled the TA area of either the 4T1 (up-regulated genes from 4T1 TA area, top bar in pink) or Cl66 (down-regulated genes from 4T1 TA area, top bar in green) of the osteolytic breast cancer model as predicted by NTP algorithm. Prediction significance (FDR < 0.2) for each sample is shown in the colored bar on the top of the hierarchical cluster. Samples predicted with FDR p > 0.2 do not resemble either the TA area or the TB interface.
The TA area was grown in a non-canonical tumor microenvironment (i.e., bone instead of breast) and as such can be considered a metastatic tumor. Nevertheless, we still expect that the gene expression profile from the TA area (grown near bone) will resemble previously reported profiles for the cell lines (grown in the breast) used in this study, especially given the fact that the primary tumor and its metastatic tumor have been reported to have similar gene expression profiles [40]. To confirm that the TA area expression signature of each cell line resembles that of primary tumors, we used a public gene expression profile of tumors grown in the breast from the 4T1 and Cl66 (66Cl4) cell lines [23]. Reassuringly, the up-regulated genes from the TA area of 4T1 cells significantly (FDR p < 0.2) predicted primary tumors from 4T1 cells and the down-regulated genes predicted tumors from Cl66 using the NTP algorithm (Figure 1B). Since the gene signature from the TA area of 4T1 cells are reported relative to Cl66 and Cl66-M2, most of the down-regulated genes represent those up-regulated in Cl66 and Cl66-M2. These results demonstrate that the gene expression profile from our microdissected TA area samples represents that of primary tumors.
In an effort to translate our findings from our mouse breast tumor model to human disease, we compared the gene expression profile from the TA area of our mouse model to that of primary human breast tumors and cancer cell lines using the NTP algorithm. Specifically, we compared microarray data from 118 primary breast tumor samples [25] to the gene expression profile from the 4T1 and Cl66 TA areas. Interestingly, 37 breast tumor samples (top bar in pink) were significantly associated with 4T1 TA area and 34 breast tumor samples were significantly associated with Cl66 TA area (top bar in green) with an FDR p < 0.2 (Figure 1C). Our analysis also predicted that 16 (top bar in pink) and 3 (top bar in green) out of 54 human breast cancer cell lines [26] resemble 4T1 and Cl66 tumors, respectively (Figure 1D). Again, the down-regulated TA area genes represent the TA area of Cl66 and Cl66-M2. This analysis predicts that it is possible to use these 19 human breast cancer cell lines in our mouse model and that similar results could be obtained.
TB interface-specific gene expression signature
In order to identify genes that are important for the interaction of breast cancer cells with the tumor microenvironment, we reanalyzed the gene expression at the TB interface and compared that profile to the gene expression profile at the TA area for each of the cell lines. Despite the expected heterogeneity in gene expression from cell line to cell line, we were able to identify 934 genes (TB signature) that were consistently different between the TB interface and the TA area. Among these, 359 were up-regulated and 575 were down-regulated with at least a 2-fold change at the TB interface across all the three cell lines. Figure 2A illustrates the top 50 known up- and down-regulated genes. The top differentially expressed genes are detailed in Tables 1 and 2.
Figure 2
The TB interface gene signature does not resemble normal bone. A. Top genes differentially expressed at the TB interface common to all the three cell lines. Data are expressed in fold-change relative to TA area: (up-regulated genes (≥2-fold) are shown in red, and down-regulated genes (≤-2-fold) are shown in green. B. The NTP algorithm was used to predict whether normal calvarial bone, ulnar bone and mandiblular bone resemble the TB interface of our mouse osteoclastic model. Samples predicted with FDR p > 0.2 do not resemble either the TA area or the TB interface. Samples and prediction significance are shown on the top of the hierarchical cluster.
Table 1
Genes up-regulated in the TB interface and their fold-change relative to the TA area.
Gene
Affymetrix probe ID
Description
Fold change
p- value
FDR p- value
Ibsp
1417484_at or
1417485_at
Integrin binding sialoprotein
7.2
0.05
0.2
Tnfsf11
(Rankl)
1419083_at or 1451944_a_at
Tumor necrosis factor (ligand) superfamily, member 11
5.3
0.006
0.1
Aftph
1426861_at
Afitiphilin
4.8
0.05
0.2
Mmp13
1417256_at
Matrix metalloproteinase 13
4.7
0.05
0.2
Anks1b
1447464_at
Ankyrin repeat and sterile alpha motif domain containing 1B
4.4
0.02
0.1
Zic1
1423477_at or 1439627_at
Zinc finger protein of the cerebellum 1
3.7
0.09
0.2
Drd1a
1455629_at or 1456051_at
Dopamine receptor 1A
3.6
0.005
0.1
Alpk2
1452478_at
Alpha-kinase 2
3.2
0.005
0.1
Smoc2
1415935_at or 1431362_a_at
SPARC related modular calcium binding 2
3.2
0.01
0.1
Table 2
Genes down-regulated in the TB interface and their fold-change relative to the TA area.
The gene expression profile of the TB interface was identified relative to the TA area, and, as such, should be enriched for transcriptional processes associated with the TB microenvironment. Indeed, three of the top four genes up-regulated at the TB interface (i.e., receptor activator of NF-kB ligand, Rankl; integrin binding sialoprotein, Ibsp; and matrix metalloproteinase 13, Mmp13) are well established as mediators of bone metastasis [7, 11–13, 15, 41, 42]. Table 1 highlights the fold-change of these genes at the TB interface as compared to the TA area (from the Affymetrix microarray profiling). Furthermore, we have previously validated the expression and function of several of these genes in our mouse model [12–15]. Collectively, these data strongly suggest that our analysis identified genes uniquely enriched in and important for the metastatic bone microenvironment.
The TB microenvironment is different than normal bone
Next, we compared the specificity of our TB specific gene set against that from the normal bone microenvironment. To this end, we used a public gene expression profile containing data for normal mouse calvarial bone, normal mouse ulnar bone and normal mouse mandibular bone (GSE13563). Our TB signature was compared against this data set using the NTP algorithm. As shown in Figure 2B, none of the calvarial or ulnar samples are enriched for the TB-signature (FDR p > 0.2), though one of the mandibular bone samples is predicted to be similar to TB microenvironment. This data demonstrates that the TB interface is genetically different from the microenvironment of normal bone.
The TB interface resembles the metastatic bone microenvironment of human breast cancer
A primary concern with any animal model is whether it accurately represents human disease. To address this, we applied NTP using the TB signature and publicly available gene expression profiles of human breast metastases (i.e., brain, lung and bone) [22]. As shown in Figure 3A, 60% of the samples from bone metastases were significantly predicted (FDR p < 0.2) to belong to the TB interface of our model. Importantly, the gene expression profiles of metastases from both brain and lung did not correlate with the TB interface data.
Figure 3
The TB interface gene signature resembles human bone metastases and osteoclastogenesis. A. The NTP algorithm was used to predict whether human breast cancer metastases (lung, bone and brain) significantly (FDR < 0.2) resemble the TB interface of our mouse osteoclastic model. B. SubMap based de novo analysis was used to compare the TB interface and TA area gene expression profiles to human breast cancer metastases samples. p-values are color coded in the bar. C. The NTP algorithm was used to predict whether the gene expression profile of the TB interface resembles osteoclasts (72 h RANKL treatment; blue) or OCPs (0 and 24 h RANKL treatment; yellow and red, respectively). Samples predicted with FDR p > 0.2 do not resemble either the TA area or the TB interface. Samples and prediction significance are shown in the colored bar on the top of the hierarchical cluster.
In addition, we also performed the Gene Set Enrichment Analysis (GSEA) [36] based SubMap algorithm [31] to predict if the TB interface gene expression profile resembles bone metastases from humans. Here, SubMap analysis with the TB signature was used to compare different human metastases samples (brain, bone and lung) to the sample sets from our mouse model (TB interface and TA area). Interestingly, de novo analysis showed that TB interface samples significantly (FDR < 0.2) resemble bone metastases samples but not lung or brain samples. TA area samples also do not resemble any of the metastases (Figure 3B). Furthermore, the Rankl and Mmp13 genes, which are up-regulated at the TB interface, are also up-regulated in the human bone metastases samples. Collectively, these data demonstrate that the osteolytic bone microenvironment in our mouse model mimics the bone microenvironment in human breast cancer but not that of other metastatic microenvironments (i.e., lung and brain metastases).
The TB interface resembles osteoclastogenesis in culture
The Rankl-mediated differentiation of osteoclast precursors (OCPs) to mature osteoclasts is a key step in breast cancer-specific bone metastasis [43]. Since Rankl is among the most highly up-regulated genes at the TB interface, we suspected that osteoclastogenesis may be occurring at the TB interface in our mouse model. To address this possibility, we performed NTP analysis using our TB signature and a publicly available gene expression profile from OCPs that have been differentiated into osteoclasts in vitro [24]. The osteoclasts used in the aforementioned data set were generated following a two stage differentiation protocol: OCPs were pretreated with macrophage colony stimulating factor (M-CSF) and then treated with human RANKL for 0, 24 or 72 h. Terminal osteoclast differentiation requires at least 72 h of RANKL treatment [24]. NTP analysis of our TB signature predicted that it was similar to OCPs treated with RANKL for 72 h with a FDR of p = 0.2. Interestingly, our TB signature did not correlate with either RANKL-untreated OCPs or those only treated for 24 h (Figure 3C). This analysis suggests that osteoclastogenesis is occurring at the TB interface in our model.
Pathways associated with the TB interface
To assess whether mechanisms that govern bone metastasis in humans are also present in our osteolytic model, we performed Gene Ontology (GO) [34]; pathway Kyoto Encyclopedia of Genes and Genomes, KEGG [44]; and Broad Institute based Molecular Signature Databases, MSigDB [45] canonical pathway enrichment analysis. The enrichment analysis was performed using the TB signature and the GlobalTest package [37, 38]. Table 3 shows GO terms significantly (FDR p < 0.05) associated with our osteolytic model. Among the GO terms significantly associated with the TB signature is TGF-β signaling (Figure 4A). Indeed, the TGF-β superfamily ligand Bmp10 is up-regulated at the TB interface in all three cell lines (greater than 2-fold in Cl66 and Cl66M2; data not shown). This would suggest that TGF-β superfamily signaling is mediated in part by the Bmp10 ligand in our model. Consistently, negative regulators (Sostdc1 and Cer1) of the TGF-β pathway are down-regulated at the TB interface and up-regulated in TA area (Figure 4A). These data suggest that Bmp-10 mediated TGF-β superfamily signaling is active at the TB interface but not in the TA area. Future studies specifically over-expressing and knocking-down members of the TGF-β signaling pathway will be required to specifically determine the role of TGF-β signaling at the TB interface.
Table 3
GO terms significantly associated with TB signature.
GO ID
Alias
BH
p-value
Statistic
GO:0007219
Notch signaling pathway
0.001
2.65e-05
78.7
GO:0007155
Cell adhesion
0.001
3.39e-06
56.3
GO:0022610
Biological adhesion
0.001
3.39e-06
56.3
GO:0007178
Transmembrane receptor protein serine/threonine kinase signaling pathway
Regulation of transforming growth factor beta receptor signaling pathway
0.008
0.0014
93.5
Figure 4
Genes involved in pathways related to the TB interface and the TA area and prediction of a therapeutic agent that targets the TB interface. p-values of pathway-specific genes enriched in the TB signature and in the TA area are plotted in green and red, respectively. Differential expression of genes in: A. the TGF-β pathway as determined using the GO database; B. the Wnt-signaling pathway using the KEGG pathway database; and C. the myeloid proliferation and self-renewal pathway [46] using MSigDB. All the enrichment analyses were done using GlobalTest package software. D. Connectivity Map analysis predicted cyclopenthiazide as a candidate drug against the TB interface gene signature with 4 instances. All four instances fall in or near the red area (represented by black lines in the bar on the bottom), which suggests that cyclopenthiazide reverses the TB interface gene expression signature.
Pathways identified using KEGG analysis that were significantly (FDR p < 0.05) associated with our osteolytic model are shown in Table 4. Interestingly, the Wnt-signaling pathway is significantly associated with the TB signature (Figure 4B), and it appears to be inhibited. Indeed, two Wnt pathway antagonists (Wif1 and Sfrp4) are expressed greater than 2-fold at the TB interface for all the mouse cell lines (data not shown). Among the four most down-regulated genes at the TB interface, relative to the TA area, two are Wnt pathway agonists (Wnt8b and Wnt2). These data suggest that the Wnt signaling pathway is active in the TA area but inhibited in the TB interface. Again, future studies specifically over-expressing and knocking-down members of the Wnt signaling pathway may be performed to further elucidate the role of Wnt signaling at the TB interface and in the TA area.
Table 4
KEGG pathways significantly associated with TB-signature.
KEGG ID
Alias
p-value
BH
Statistic
04060
Cytokine-cytokine receptor interaction
7.60e-05
0.003
75.3
01100
Metabolic pathways
8.88e-05
0.003
52.8
05016
Huntington's disease
9.74e-05
0.003
57.9
04310
Wnt signaling pathway
1.07e-04
0.003
72.0
00061
Fatty acid biosynthesis
1.26e-04
0.003
98.2
04912
GnRH signaling pathway
1.26e-04
0.003
57.6
04020
Calcium signaling pathway
3.59e-04
0.007
67.1
05010
Alzheimer's disease
4.15e-04
0.007
66.5
04650
Natural killer cell mediated cytotoxicity
5.35e-04
0.008
60.9
04012
ErbB signaling pathway
6.28e-04
0.009
59.9
We also performed enrichment analysis of the TB signature using MSigDB canonical pathway database and GlobalTest package [37, 38]. Among the pathways significantly associated with the TB interface (Table 5) were myeloid proliferation and self-renewal [46]. Consistently, two genes (Rankl and Myb) highly expressed at the TB interface were significantly associated with this pathway (Figure 4C) [45]. This data further corroborates the NTP analysis comparing osteoclasts to our TB signature (Figure 3C) and provides additional evidence for a role of osteoclastogenesis at the TB interface.
Table 5
MSigDB pathway signatures significantly associated with TB signature.
Alias
p-value
BH
Statistic
INTRINSICPATHWAY
0.0004
0.0007
96.9
BLOOD_CLOTTING_CASCADE
0.0004
0.0007
96.9
HSIAO_LIVER_SPECIES_GENES
0.0004
0.0007
96.9
TPA_SENS_MIDDLE_DN
0.0004
0.0007
96.9
TPA__SENS_LATE_DN
0.0004
0.0007
96.9
TPA__SENS_EARLY_DN
0.0004
0.0007
96.9
HSA04610_COMPLEMENT_AND_COAGULATION_CASCADES
0.0004
0.0007
96.9
BROWN_MYELOID_PROLIF_AND_SELF_RENEWAL
0.0062
0.010
60.2
BROWN_GRAN_MONO_DIFFERENTIATION
0.0121
0.017
38.9
KAMMINGA_EZH2_TARGETS
0.0521
0.068
65.2
Prediction and validation of therapeutic targets using the TB signature
To predict a therapeutic agent that specifically targets the TB interface, we queried Connectivity Map database [39] using the TB gene signature. Probeset identifiers from the Affymetrix Mouse Genome 430A 2.0 array were mapped to Affymetrix Human Genome U133A array. This was then used to query the Connectivity Map database. Of the 6,100 potential therapeutic candidates, cyclopenthiazide had the most highly significant negative mean connectivity scores. In other words, cyclopenthiazide was predicted to reverse the gene expression signature of the TB interface (Figure 4D). This analysis suggests that cyclopenthiazide may be a potential agent against human osteoclastic bone metastasis. Future studies aim to address this possibility by therapeutically dosing our mouse model with cyclopenthiazide and monitoring for changes in the TB microenvironment.
Discussion
Mouse Model of the Osteolytic Microenvironment in Breast Cancer
Animal models that faithfully recapitulate aspects of human breast cancer-specific bone metastasis provide powerful tools to study the complex molecular mechanism(s) by which breast cancer cells metastasize to and interact with the bone microenvironment [47, 48]. Previously, we developed mouse models of bone osteolysis for prostate and breast cancer by implanting syngeneic tumor cells onto the calvaria of animals using a simple surgical technique. These models produced osteolytic lesions at the TB interface of the implant region, thereby allowing us to explore the cellular and molecular interactions between malignant cells and skeletal tissue [11, 12, 15]. Because the tumor cells are implanted directly into the bone microenvironment (albeit, at an atypical location for breast cancer bone metastasis), it was important to confirm that the interactions observed in our model reflect those observed between metastatic human breast cells and the bone microenvironment. Building on our previous work, we now demonstrate that the TB microenvironment in our model appears very similar to that of human breast cancer bone metastases based on gene expression data. As such, this mouse model can be readily used to study the cellular and molecular mechanisms driving human breast cancer-metastasis and osteolysis. Furthermore, this model also provides a powerful preclinical setting to test cyclopenthiazide and other therapeutic agents that specifically target breast cancer osteolysis.
Gene Expression Profile Analysis
There has been tremendous growth in both the development of high-throughput microarray technology to measure gene expression in tissue and cells and in high-dimensional methods to analyze such data [49, 50]. Together with this, many of the gene expression microarray data sets generated from different labs are now available in open-access databases [21, 51], which enables the comparison and integration of data acquired from different batches, laboratories and experimental platforms [52]. Importantly, this has opened up opportunities to perform cross-species comparisons of mouse models and human disease [30].
In the current study, we applied microarray technology to generate a signature specific to the TB interface of our mouse model. The robustness of our TB-signature is supported by the fact that it was derived from a common set of genes regulated at the TB interface across a heterogeneous set of three mouse breast cancer cell lines. Combining gene expression profiling and molecular pathology, we demonstrated that the TB interface of our model truly represents the tumor microenvironment and not the normal bone microenvironment. Subsequent cross-species comparative transcriptomic analysis demonstrated that many human bone metastases samples are associated with the TB interface in a statistically significant manner. Importantly, there was no association between our breast TB interface and human brain or lung metastases. Together, these data demonstrate that our model specifically mimics human breast (and not lung or brain) cancer bone metastases. Furthermore, analysis of a panel of human breast cancer cell lines predicted 16 that have similar gene expression characteristics to those of the 4T1 tumors. This suggests that our osteolytic model may be adapted to study human breast cancer bone metastasis directly using any of these 16 human cell lines.
Pathways involved in the Breast Cancer Osteolytic Microenvironment
The TGF-β pathway has a well established role in bone metastasis [53], and previously we demonstrated the importance of TGF-β signaling in the TB interface using our model [12]. Here, we demonstrate that the TGF-β receptor I is expressed and that the TGF-β pathway is active in tumor cells and osteoclasts at the TB interface. On the other hand, TGF-β signaling is not active in the TA area [12]. Interestingly, the TGF-β signaling ligand Bmp10 [54] is highly expressed at the TB interface and TGF-β pathway inhibitors (i.e., Sostdc1 [55] and Cer1 [56]) are suppressed at the TB interface. These data suggest that Bmp-10 is responsible for mediating TGF-β pathway activation at the TB interface.
The canonical and noncanonical Wnt signaling pathways are involved in the formation, growth and development of normal bone [57] and bone metastasis [58]. Activation of canonical Wnt signaling through β-catenin both promotes osteoblast differentiation and inhibits osteoclast formation and bone resorption [59]. Our KEGG pathway enrichment analysis showed a significant association of the Wnt signaling pathway at the TB interface. Indeed, we observed that Wnt pathway antagonists - Wif1, which is associated with decreased bone mineral density [60], and Sfrp4, which is associated with the suppression of osteoblast proliferation [61] - were over-expressed at the TB interface. Furthermore, we observed a down-regulation of the Wnt pathway ligands Wnt2 [62] and Wnt8b [63] at the TB interface relative to the TA area. Together these data suggest that our mouse model exhibits (i) Wnt pathway activation in the TA area and (ii) increased bone resorption and suppressed bone formation (at least in part through Wnt pathway activation) at the TB interface.
Osteoclasts are derived from hematopoietic precursor cells of the myeloid lineage upon CSF-1 stimulation followed by RANKL-mediated maturation [43]. In our current study, we used a publicly available microarray dataset from RANKL-differentiated OCPs. Interestingly, we found that the gene expression profile of in vitro differentiated osteoclasts (72 h RANKL treatment) was similar to that of the TB interface. In addition, pathway analysis using the MSigDB showed an enrichment of the TB-signature in a myeloid cell line model [46] (Figure 4C). Overall, these results suggest that osteolysis is operative at the TB interface of our mouse model.
Prediction of a Therapeutic Agent that Targets the TB interface
The identification of new therapeutic agents that inhibit the establishment of tumor cells in the TB microenvironment will benefit patients with breast cancer bone metastases [5]. This will require a thorough understanding of the mechanisms governing breast-to-bone metastasis to determine appropriate biological targets for intervention. In one example, we previously demonstrated that TGF-β signaling activity may provide such a target as pathway attenuation in our mouse model led to a reduction in breast tumor-induced osteolysis [12]. Herein, we used gene expression profiles from our mouse model and Connectivity Map database to find therapeutic agents that target the TB interface, rather than a given pathway.
The advantage of Connectivity Map database is that it can predict potential therapeutic agents based solely on gene signatures [39]. In the current study, our query of Connectivity Map database with the TB signature flagged cyclopenthiazide in the MCF7 cell line (Figure 4D). This analysis suggests that cyclopenthiazide has the potential to inhibit the establishment of breast cancer cells at TB interface.
Thiazides comprise a class of diuretic agents (of which cyclopenthiazide is a member) that are traditionally used to treat hypertension and edema [64]. Although thiazides have not been widely viewed as therapeutic agents for bone metastasis, reports abound noting that treatment of hypertension using thiazides has the beneficial side effect of strengthening bone [65–69]. Furthermore, Devorak et al. have demonstrated that the bone strengthening activity of thiazides results from their direct action on OCPs, where thiazide analogs are able to directly induce osteoblast differentiation [70]. These data suggest that cyclopenthiazide may be a useful agent against osteoclastic bone metastasis. Future efforts are aimed at validating this prediction in the osteolytic mouse model. This study serves as an example of how mouse breast cancer-specific osteolytic models and gene expression analysis can be used to identify treatment strategies for human disease.
Conclusions
In summary, we have demonstrated that the TB microenvironment in our mouse model of osteolytic breast cancer metastasis is highly similar to that of human breast cancer-to-bone metastases. Furthermore, gene expression profile analysis of tumors from this model: (i) identified a TB interface specific gene signature; (ii) revealed signaling pathways that were differentially activated at the TB interface and TA area; (iii) demonstrated a role for osteoclasts in metastatic osteolysis; and (iv) predicted a novel therapeutic agent that specifically targets the TB interface. These data clearly demonstrate that this mouse model can be used to study the cellular and molecular mechanisms driving human breast cancer-to-bone metastasis and osteolysis. Moreover, this model also provides a powerful preclinical setting to test thiazides and other therapeutic agents that specifically target breast cancer osteolysis.
Declarations
Acknowledgements
We thank Dr. Ke Cheng for editing the manuscript. This work was supported in part by grants CA71781 (R.K.S.), and Cancer Center Support Grant (P30CA036727) from National Cancer Institute, National Institute of Health and Nebraska Research Initiative Molecular Therapeutics Program (R.K.S.). This work was supported in part by Susan G. Komen for the Cure grant KG090860 (R.K.S).
Authors’ Affiliations
(1)
Department of Pathology and Microbiology, University of Nebraska Medical Center
(2)
Department of Molecular Toxicology, Graduate School of Medical Sciences, Nagoya City University
(3)
Department of Systems Biology, Harvard Medical School
(4)
Division of Signal Transduction, Beth Israel Deaconess Medical Center
(5)
Department of Life Sciences, Lawrence Berkeley National Lab
(6)
Swiss Institute for Experimental Cancer Research (ISREC), Ecole Polytechnique Federale de Lausanne (EPFL)
References
Roodman GD: Mechanisms of bone metastasis. The New England journal of medicine. 2004, 350 (16): 1655-1664.View ArticlePubMedGoogle Scholar
Arguello F, Baggs RB, Frantz CN: A murine model of experimental metastasis to bone and bone marrow. Cancer research. 1988, 48 (23): 6876-6881.PubMedGoogle Scholar
Harms JF, Welch DR: MDA-MB-435 human breast carcinoma metastasis to bone. Clinical & experimental metastasis. 2003, 20 (4): 327-334.View ArticleGoogle Scholar
Thomas RJ, Guise TA, Yin JJ, Elliott J, Horwood NJ, Martin TJ, Gillespie MT: Breast cancer cells interact with osteoblasts to support osteoclast formation. Endocrinology. 1999, 140 (10): 4451-4458.PubMedGoogle Scholar
Yoneda T, Williams PJ, Hiraga T, Niewolna M, Nishimura R: A bone-seeking clone exhibits different biological properties from the MDA-MB-231 parental human breast cancer cells and a brain-seeking clone in vivo and in vitro. J Bone Miner Res. 2001, 16 (8): 1486-1495.View ArticlePubMedGoogle Scholar
Lynch CC, Hikosaka A, Acuff HB, Martin MD, Kawai N, Singh RK, Vargo-Gogola TC, Begtrup JL, Peterson TE, Fingleton B, Shirai T, Matrisian LM, Futakuchi M: MMP-7 promotes prostate cancer-induced osteolysis via the solubilization of RANKL. Cancer cell. 2005, 7 (5): 485-496.View ArticlePubMedGoogle Scholar
Futakuchi M, Nannuru KC, Varney ML, Sadanandam A, Nakao K, Asai K, Shirai T, Sato SY, Singh RK: Transforming growth factor-beta signaling at the tumor-bone interface promotes mammary tumor growth and osteoclast activation. Cancer science. 2009, 100 (1): 71-81.View ArticlePubMedGoogle Scholar
Nannuru KC, Futakuchi M, Sadanandam A, Wilson TJ, Varney ML, Myers KJ, Li X, Marcusson EG, Singh RK: Enhanced expression and shedding of receptor activator of NF-kappaB ligand during tumor-bone interaction potentiates mammary tumor-induced osteolysis. Clinical & experimental metastasis. 2009, 26 (7): 797-808.View ArticleGoogle Scholar
Nannuru KC, Futakuchi M, Varney ML, Vincent TM, Marcusson EG, Singh RK: Matrix metalloproteinase (MMP)-13 regulates mammary tumor-induced osteolysis by activating MMP9 and transforming growth factor-beta signaling at the tumor-bone interface. Cancer research. 2010, 70 (9): 3494-3504.View ArticlePubMedPubMed CentralGoogle Scholar
Wilson TJ, Nannuru KC, Futakuchi M, Sadanandam A, Singh RK: Cathepsin G enhances mammary tumor-induced osteolysis by generating soluble receptor activator of nuclear factor-kappaB ligand. Cancer research. 2008, 68 (14): 5803-5811.View ArticlePubMedGoogle Scholar
Aslakson CJ, Miller FR: Selective events in the metastatic process defined by analysis of the sequential dissemination of subpopulations of a mouse mammary tumor. Cancer research. 1992, 52 (6): 1399-1405.PubMedGoogle Scholar
Murphy BO, Joshi S, Kessinger A, Reed E, Sharp JG: A murine model of bone marrow micrometastasis in breast cancer. Clinical & experimental metastasis. 2002, 19 (7): 561-569.View ArticleGoogle Scholar
Varney ML, Singh S, Backora M, Chen Z, Singh RK: Lymphangiogenesis and anti-tumor immune responses. Current molecular medicine. 2009, 9 (6): 694-701.View ArticlePubMedGoogle Scholar
Wilson CL, Miller CJ: Simpleaffy: a BioConductor package for Affymetrix Quality Control and data analysis. Bioinformatics (Oxford, England). 2005, 21 (18): 3683-3685.View ArticleGoogle Scholar
Simon R, Lam A, Li MC, Ngan M, Menenzes S, Zhao Y: Analysis of gene expression data using BRB-ArrayTools. Cancer informatics. 2007, 3: 11-17.PubMedPubMed CentralGoogle Scholar
Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, Evangelista C, Kim IF, Soboleva A, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Muertter RN, Edgar R: NCBI GEO: archive for high-throughput functional genomic data. Nucleic acids research. 2009, D885-890. 37 DatabaseGoogle Scholar
Zhang XH, Wang Q, Gerald W, Hudis CA, Norton L, Smid M, Foekens JA, Massague J: Latent bone metastasis in breast cancer tied to Src-dependent survival signals. Cancer cell. 2009, 16 (1): 67-78.View ArticlePubMedPubMed CentralGoogle Scholar
Lou Y, Preobrazhenska O, auf dem Keller U, Sutcliffe M, Barclay L, McDonald PC, Roskelley C, Overall CM, Dedhar S: Epithelial-mesenchymal transition (EMT) is not sufficient for spontaneous murine breast cancer metastasis. Dev Dyn. 2008, 237 (10): 2755-2768.View ArticlePubMedGoogle Scholar
Zhao B, Takami M, Yamada A, Wang X, Koga T, Hu X, Tamura T, Ozato K, Choi Y, Ivashkiv LB, Takayanagi H, Kamijo R: Interferon regulatory factor-8 regulates bone metabolism by suppressing osteoclastogenesis. Nature medicine. 2009, 15 (9): 1066-1071.View ArticlePubMedPubMed CentralGoogle Scholar
Chin K, DeVries S, Fridlyand J, Spellman PT, Roydasgupta R, Kuo WL, Lapuk A, Neve RM, Qian Z, Ryder T, Chen F, Feiler H, Tokuyasu T, Kingsley C, Dairkee S, Meng Z, Chew K, Pinkel D, Jain A, Ljung BM, Esserman L, Albertson DG, Waldman FM, Gray JW: Genomic and transcriptional aberrations linked to breast cancer pathophysiologies. Cancer cell. 2006, 10 (6): 529-541.View ArticlePubMedGoogle Scholar
Neve RM, Chin K, Fridlyand J, Yeh J, Baehner FL, Fevr T, Clark L, Bayani N, Coppe JP, Tong F, Speed T, Spellman PT, DeVries S, Lapuk A, Wang NJ, Kuo WL, Stilwell JL, Pinkel D, Albertson DG, Waldman FM, McCormick F, Dickson RB, Johnson MD, Lippman M, Ethier S, Gazdar A, Gray JW: A collection of breast cancer cell lines for the study of functionally distinct cancer subtypes. Cancer cell. 2006, 10 (6): 515-527.View ArticlePubMedPubMed CentralGoogle Scholar
Hoshida Y, Villanueva A, Kobayashi M, Peix J, Chiang DY, Camargo A, Gupta S, Moore J, Wrobel MJ, Lerner J, Reich M, Chan JA, Glickman JN, Ikeda K, Hashimoto M, Watanabe G, Daidone MG, Roayaie S, Schwartz M, Thung S, Salvesen HB, Gabriel S, Mazzaferro V, Bruix J, Friedman SL, Kumada H, Llovet JM, Golub TR: Gene expression in fixed tissues and outcome in hepatocellular carcinoma. The New England journal of medicine. 2008, 359 (19): 1995-2004.View ArticlePubMedPubMed CentralGoogle Scholar
Hoshida Y: Nearest template prediction: a single-sample-based flexible class prediction with confidence assessment. PloS one. 2010, 5 (11): e15543-View ArticlePubMedPubMed CentralGoogle Scholar
Benito M, Parker J, Du Q, Wu J, Xiang D, Perou CM, Marron JS: Adjustment of systematic microarray data biases. Bioinformatics (Oxford, England). 2004, 20 (1): 105-114.View ArticleGoogle Scholar
Herschkowitz JI, Simin K, Weigman VJ, Mikaelian I, Usary J, Hu Z, Rasmussen KE, Jones LP, Assefnia S, Chandrasekharan S, Backlund MG, Yin Y, Khramtsov AI, Bastein R, Quackenbush J, Glazer RI, Brown PH, Green JE, Kopelovich L, Furth PA, Palazzo JP, Olopade OI, Bernard PS, Churchill GA, Van Dyke T, Perou CM: Identification of conserved gene expression features between murine mammary carcinoma models and human breast tumors. Genome biology. 2007, 8 (5): R76-View ArticlePubMedPubMed CentralGoogle Scholar
Hoshida Y, Brunet JP, Tamayo P, Golub TR, Mesirov JP: Subclass mapping: identifying common subtypes in independent disease data sets. PloS one. 2007, 2 (11): e1195-View ArticlePubMedPubMed CentralGoogle Scholar
Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proceedings of the National Academy of Sciences of the United States of America. 1998, 95 (25): 14863-14868.View ArticlePubMedPubMed CentralGoogle Scholar
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nature genetics. 2000, 25 (1): 25-29.View ArticlePubMedPubMed CentralGoogle Scholar
Kanehisa M, Goto S, Furumichi M, Tanabe M, Hirakawa M: KEGG for representation and analysis of molecular networks involving diseases and drugs. Nucleic acids research. D355-360. 38 DatabaseGoogle Scholar
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proceedings of the National Academy of Sciences of the United States of America. 2005, 102 (43): 15545-15550.View ArticlePubMedPubMed CentralGoogle Scholar
Goeman JJ, Oosting J, Cleton-Jansen AM, Anninga JK, van Houwelingen HC: Testing association of a pathway with survival using gene expression data. Bioinformatics (Oxford, England). 2005, 21 (9): 1950-1957.View ArticleGoogle Scholar
Goeman JJ, van de Geer SA, de Kort F, van Houwelingen HC: A global test for groups of genes: testing association with a clinical outcome. Bioinformatics (Oxford, England). 2004, 20 (1): 93-99.View ArticleGoogle Scholar
Lamb J, Crawford ED, Peck D, Modell JW, Blat IC, Wrobel MJ, Lerner J, Brunet JP, Subramanian A, Ross KN, Reich M, Hieronymus H, Wei G, Armstrong SA, Haggarty SJ, Clemons PA, Wei R, Carr SA, Lander ES, Golub TR: The Connectivity Map: using gene-expression signatures to connect small molecules, genes, and disease. Science (New York, NY. 2006, 313 (5795): 1929-1935.View ArticleGoogle Scholar
Ramaswamy S, Ross KN, Lander ES, Golub TR: A molecular signature of metastasis in primary solid tumors. Nature genetics. 2003, 33 (1): 49-54.View ArticlePubMedGoogle Scholar
Kong YY, Yoshida H, Sarosi I, Tan HL, Timms E, Capparelli C, Morony S, Oliveira-dos-Santos AJ, Van G, Itie A, Khoo W, Wakeham A, Dunstan CR, Lacey DL, Mak TW, Boyle WJ, Penninger JM: OPGL is a key regulator of osteoclastogenesis, lymphocyte development and lymph-node organogenesis. Nature. 1999, 397 (6717): 315-323.View ArticlePubMedGoogle Scholar
Lacey DL, Timms E, Tan HL, Kelley MJ, Dunstan CR, Burgess T, Elliott R, Colombero A, Elliott G, Scully S, Hsu H, Sullivan J, Hawkins N, Davy E, Capparelli C, Eli A, Qian YX, Kaufman S, Sarosi I, Shalhoub V, Senaldi G, Guo J, Delaney J, Boyle WJ: Osteoprotegerin ligand is a cytokine that regulates osteoclast differentiation and activation. Cell. 1998, 93 (2): 165-176.View ArticlePubMedGoogle Scholar
Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M: KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic acids research. 1999, 27 (1): 29-34.View ArticlePubMedPubMed CentralGoogle Scholar
Subramanian A, Kuehn H, Gould J, Tamayo P, Mesirov JP: GSEA-P: a desktop application for Gene Set Enrichment Analysis. Bioinformatics (Oxford, England). 2007, 23 (23): 3251-3253.View ArticleGoogle Scholar
Brown AL, Wilkinson CR, Waterman SR, Kok CH, Salerno DG, Diakiw SM, Reynolds B, Scott HS, Tsykin A, Glonek GF, Goodall GJ, Solomon PJ, Gonda TJ, D'Andrea RJ: Genetic regulators of myelopoiesis and leukemic signaling identified by gene profiling and linear modeling. Journal of leukocyte biology. 2006, 80 (2): 433-447.View ArticlePubMedGoogle Scholar
Kang Y, Siegel PM, Shu W, Drobnjak M, Kakonen SM, Cordon-Cardo C, Guise TA, Massague J: A multigenic program mediating breast cancer metastasis to bone. Cancer cell. 2003, 3 (6): 537-549.View ArticlePubMedGoogle Scholar
Reddi AH, Roodman D, Freeman C, Mohla S: Mechanisms of tumor metastasis to the bone: challenges and opportunities. J Bone Miner Res. 2003, 18 (2): 190-194.View ArticlePubMedGoogle Scholar
Stoughton RB: Applications of DNA microarrays in biology. Annual review of biochemistry. 2005, 74: 53-82.View ArticlePubMedGoogle Scholar
Wang Y, Miller DJ, Clarke R: Approaches to working in high-dimensional data spaces: gene expression microarrays. British journal of cancer. 2008, 98 (6): 1023-1028.View ArticlePubMedPubMed CentralGoogle Scholar
Sarkans U, Parkinson H, Lara GG, Oezcimen A, Sharma A, Abeygunawardena N, Contrino S, Holloway E, Rocca-Serra P, Mukherjee G, Shojatalab M, Kapushesky M, Sansone SA, Farne A, Rayner T, Brazma A: The ArrayExpress gene expression database: a software engineering and implementation perspective. Bioinformatics (Oxford, England). 2005, 21 (8): 1495-1501.View ArticleGoogle Scholar
Allison DB, Cui X, Page GP, Sabripour M: Microarray data analysis: from disarray to consolidation and consensus. Nat Rev Genet. 2006, 7 (1): 55-65.View ArticlePubMedGoogle Scholar
Neuhaus H, Rosen V, Thies RS: Heart specific expression of mouse BMP-10 a novel member of the TGF-beta superfamily. Mechanisms of development. 1999, 80 (2): 181-184.View ArticlePubMedGoogle Scholar
Yanagita M, Oka M, Watabe T, Iguchi H, Niida A, Takahashi S, Akiyama T, Miyazono K, Yanagisawa M, Sakurai T: USAG-1: a bone morphogenetic protein antagonist abundantly expressed in the kidney. Biochemical and biophysical research communications. 2004, 316 (2): 490-500.View ArticlePubMedGoogle Scholar
Yanagita M: BMP modulators regulate the function of BMP during body patterning and disease progression. BioFactors (Oxford, England). 2009, 35 (2): 113-119.View ArticleGoogle Scholar
Krishnan V, Bryant HU, Macdougald OA: Regulation of bone mass by Wnt signaling. The Journal of clinical investigation. 2006, 116 (5): 1202-1209.View ArticlePubMedPubMed CentralGoogle Scholar
Voorzanger-Rousselot N, Goehrig D, Journe F, Doriath V, Body JJ, Clezardin P, Garnero P: Increased Dickkopf-1 expression in breast cancer bone metastases. British journal of cancer. 2007, 97 (7): 964-970.PubMedPubMed CentralGoogle Scholar
Glass DA, Bialek P, Ahn JD, Starbuck M, Patel MS, Clevers H, Taketo MM, Long F, McMahon AP, Lang RA, Karsenty G: Canonical Wnt signaling in differentiated osteoblasts controls osteoclast differentiation. Developmental cell. 2005, 8 (5): 751-764.View ArticlePubMedGoogle Scholar
Baron R, Rawadi G: Targeting the Wnt/beta-catenin pathway to regulate bone formation in the adult skeleton. Endocrinology. 2007, 148 (6): 2635-2643.View ArticlePubMedGoogle Scholar
Nakanishi R, Akiyama H, Kimura H, Otsuki B, Shimizu M, Tsuboyama T, Nakamura T: Osteoblast-targeted expression of Sfrp4 in mice results in low bone mass. J Bone Miner Res. 2008, 23 (2): 271-277.View ArticlePubMedGoogle Scholar
Katoh M: Frequent up-regulation of WNT2 in primary gastric cancer and colorectal cancer. International journal of oncology. 2001, 19 (5): 1003-1007.PubMedGoogle Scholar
Lako M, Strachan T, Curtis AR, Lindsay S: Isolation and characterization of WNT8B, a novel human Wnt gene that maps to 10q24. Genomics. 1996, 35 (2): 386-388.View ArticlePubMedGoogle Scholar
Ernst ME, Moser M: Use of diuretics in patients with hypertension. The New England journal of medicine. 2009, 361 (22): 2153-2164.View ArticlePubMedGoogle Scholar
Brickman AS, Massry SG, Coburn JW: changes in serum and urinary calcium during treatment with hydrochlorothiazide: studies on mechanisms. The Journal of clinical investigation. 1972, 51 (4): 945-954.View ArticlePubMedPubMed CentralGoogle Scholar
Gamba G: Molecular physiology and pathophysiology of electroneutral cation-chloride cotransporters. Physiological reviews. 2005, 85 (2): 423-493.View ArticlePubMedGoogle Scholar
Jones G, Nguyen T, Sambrook PN, Eisman JA: Thiazide diuretics and fractures: can meta-analysis help?. J Bone Miner Res. 1995, 10 (1): 106-111.PubMedGoogle Scholar
LaCroix AZ, Wienpahl J, White LR, Wallace RB, Scherr PA, George LK, Cornoni-Huntley J, Ostfeld AM: Thiazide diuretic agents and the incidence of hip fracture. The New England journal of medicine. 1990, 322 (5): 286-290.View ArticlePubMedGoogle Scholar
Wasnich RD, Benfante RJ, Yano K, Heilbrun L, Vogel JM: Thiazide effect on the mineral content of bone. The New England journal of medicine. 1983, 309 (6): 344-347.View ArticlePubMedGoogle Scholar
Dvorak MM, De Joussineau C, Carter DH, Pisitkun T, Knepper MA, Gamba G, Kemp PJ, Riccardi D: Thiazide diuretics directly induce osteoblast differentiation and mineralized nodule formation by interacting with a sodium chloride co-transporter in bone. J Am Soc Nephrol. 2007, 18 (9): 2509-2516.View ArticlePubMedPubMed CentralGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
