Homeobox transcription factors are developmentally regulated genes that play crucial roles in tissue patterning. Homeobox C6 (HOXC6) is overexpressed in prostate cancers and correlated with cancer progression, but the downstream targets of HOXC6 are largely unknown. We have performed genome-wide localization analysis to identify promoters bound by HOXC6 in prostate cancer cells. This analysis identified 468 reproducibly bound promoters whose associated genes are involved in functions such as cell proliferation and apoptosis. We have complemented these data with expression profiling of prostates from mice with homozygous disruption of the Hoxc6 gene to identify 31 direct regulatory target genes of HOXC6. We show that HOXC6 directly regulates expression of bone morphogenic protein 7, fibroblast growth factor receptor 2, insulin-like growth factor binding protein 3, and platelet-derived growth factor receptor α (PDGFRA) in prostate cells and indirectly influences the Notch and Wnt signaling pathways in vivo. We further show that inhibition of PDGFRA reduces proliferation of prostate cancer cells, and that overexpression of HOXC6 can overcome the effects of PDGFRA inhibition. HOXC6 regulates genes with both oncogenic and tumor suppressor activities as well as several genes such as CD44 that are important for prostate branching morphogenesis and metastasis to the bone microenvironment. [Cancer Res 2008;68(6):1988–96]

Homeobox (HOX) transcription factors are developmentally regulated genes that play crucial roles in tissue patterning (1). However, HOX gene expression has also been observed in many adult tissues with possible roles in cellular differentiation. Dysregulation of HOX genes has been observed in many human cancers including leukemias, and solid tumors of the breast, colon, lung, kidney, ovary, and prostate (2). HOX family member HOXC6 is expressed in osteosarcomas, medulloblastomas, as well as carcinomas of the breast, lung, and prostate, and is overexpressed in the LNCaP prostate cancer cell line (3). In a recent prostate cancer study, HOXC6 was identified as the gene most strongly correlated with Gleason grade out of a newly identified 16 gene signatures (4).

Previously, we have shown that siRNA knockdown of HOXC6 expression induces apoptosis, and overexpression of HOXC6 results in increased proliferation and decreased apoptosis in LNCaP cells (5). This previous analysis identified T-cell receptor alternate reading frame protein (TARP), insulin-like growth factor binding protein 3 (IGFBP3), and neutral endopeptidase/membrane metalloendopeptidase as biologically relevant HOXC6 target genes due to their change in expression levels in prostate tumor normal tissue as well as when HOXC6 levels were perturbed in LNCaP cells. However, to date, few direct targets of HOXC6 are known and none have been identified in prostate cancer cells. Indeed, there is evidence that HOX genes function differently in different tissues, so targets identified in one cell type might not reflect the targets in another (6).

Here, we have developed a comprehensive transcriptional network of genes under the control of HOXC6 in prostate cancer. We developed a stable LNCaP cell line that expresses epitope-tagged HOXC6-HA protein and used these cells in chromatin immunoprecipitation (ChIP) followed by microarray analysis (ChIP-chip). This analysis identified 468 reproducibly bound promoters whose associated genes are involved in functions such as cell proliferation and apoptosis. We correlated HOXC6-induced gene expression changes with gene promoters bound by HOXC6 and identified several key developmental ligands and receptors under direct transcriptional regulation by HOXC6. We further showed the biological relevance of many identified direct target genes by performing gene expression profiling on mature prostates of Hoxc6−/− and Hoxc6+/+ mice. Bone morphogenic protein 7 (BMP7), fibroblast growth factor receptor 2 (FGFR2), and platelet-derived growth factor receptor α (PDGFRA), among others, were identified in these studies as in vivo direct regulatory targets of HOXC6 in prostate tissue. Along with elucidation of direct targets of HOXC6, we have also identified several biological pathways, including the Notch and Wnt signaling pathways that are indirectly affected by HOXC6 expression and are important in cancer progression. Additionally, several of the targets and pathways regulated by HOXC6 effect prostate duct branching morphogenesis and bone development, and may facilitate metastasis of prostate cancers to the bone microenvironment.

Lentiviral vectors and cell culture. A hemagglutinin (HA) tag was inserted by restriction enzyme digestion onto the 3′ end of a HOXC6 cDNA clone. The resulting HOXC6-HA open reading frame was subcloned into pUBQ-IRES-eYFP-ΔU3 (a gift from Dr. Hihn Ly, Emory University, Atlanta, GA) producing pUBQ-HOXC6-IRES-eYFP-ΔU3. The virus-producing cell line, HEK293FT, was transfected with either empty vector or the HOXC6-HA–containing construct along with viral-packaging vectors, pVSVG (Clontech) and pCMV_R8.2. At 48 h posttransfection, virus-containing supernatant was collected, filtered through a 45-μm filter, and used to infect either LNCaP or RWPE-1 cells. Infected cells were isolated from noninfected cells by flow sorting for yellow fluorescent protein (YFP) expression at the Emory University flow core.4

Immunoprecipitation and Western blot. Cells and tissues were lysed by mechanical disruption in immunoprecipitation Lysis Buffer [20 mmol/L Tris (pH 8.0), 10% glycerol, 137 mmol/L NaCl, 10 mmol/L EDTA, and 1% NP40] and quantitated using Bradford Protein Assay (Bio-Rad). Immunoprecipitations were performed using a mouse monoclonal antibody (mAb) to the 12CA5 HA tag or normal mouse IgG and protein-G Sepharose beads. Immunoblots were probed with polyclonal antibodies to HOXC6 (5).

ChIP. Two 150-mm dishes each of HOXC6-HA-LNCaP or HOXC6-HA-RWPE-1 stable and HOXC6-HA-YFP-LNCaP or HOXC6-HA-RWPE-1 control cell lines were grown to 80% confluency. The cells were fixed with formaldehyde, lysed, and sonicated as previously described (7). Sonicated chromatin was immunoprecipitated with 12CA5 ascites or mouse IgG overnight and collected with Dynal M280 sheep anti-mouse IgG magnetic beads for 2 h at 4°C. The beads were washed, complexes eluted, and the DNA was purified as previously described (7). Purified DNA was subjected to PCR using the primers in Supplementary Table S1.

ChIP-chip data analysis. DNA isolated from ChIP with 12CA5 mAb from the HOXC6-HA-LNCaP and YFP-LNCaP stable cell lines, along with input control, was amplified through linker-mediated PCR as previously described (8). Amplified DNA was labeled and hybridized to the Nimblegen 25K human promoter array, which tiles ∼5 kb upstream of 25,000 human transcribed sequences. The data from each hybridization was Z-score normalized, and ratios of immunoprecipitation to input signal were determined. ChIPOTle analysis, which identifies peaks using a sliding window approach, was performed for each array (9). The window was set to the average shear length of the DNA (500 bp) and the step was set to the probe size (50 bp). The peaks of signal over the input DNA background identified anti-HA 12CA5–enriched DNA. Each cell line was analyzed in duplicate. Chip-chip data have been submitted to the Gene Expression Omnibus (GEO) repository (accession number, GSE9772).

Luciferase assays. Promoters were amplified from human genomic DNA and directionally cloned into pGL3-basic (Promega) using restrictions sites incorporated into the primers. LNCaP cells were plated at 70% confluence and transfected using DMRIE-C (Invitrogen) with luciferase expression vector, Renilla expression control vector, and either pCDNA3.1 or pCDNA3.1-HOXC6 expression vectors. Forty-eight hours posttransfection, cells were rinsed with 1× PBS followed lysis in Passive Lysis buffer (Promega). Fifty micrograms of each lysate was used in a Dual-Luciferase assay according to the manufacturer (Promega) and normalized for transfection efficiency by Renilla light units. Fold change upon HOXC6 expression was computed relative to the pcDNA control and normalized to the fold change observed on the empty pGL3-basic vector.

Mouse tissue collection and genotyping. All animal procedures were performed according to protocols approved by the Institutional Animal Care and Use Committee at Emory University. C57Bl/6 Hoxc6 knockout mice are null for HOXC6 expression due to disruption of the Hoxc6 gene homeodomain (10). Hoxc6−/− and Hoxc6+/+ male mice were generated from intercrosses of heterozygous parents. Mice were genotyped using PCR amplification from tail DNA as previously described (10). Prostates were dissected from Hoxc6−/− and Hoxc6+/+ male mice at postnatal week 32 and immediately snap frozen in liquid nitrogen and stored at −80°C before RNA and protein isolation.

Microarray analysis. Total RNA was extracted from three Hoxc6−/− and three Hoxc6+/+ mouse prostates using RNeasy kit (Qiagen) according to the manufacturer's instructions. All RNA samples were DNase treated (Rnase-free; Invitrogen) to eliminate residual genomic DNA. RNA was submitted to the Winship Cancer Institute Microarray Core Facility at Emory University, subjected to quality control analysis, labeled using the TotalPrep RNA labeling kit (Ambion), and hybridized to the Illumina Mouse-6 BeadChip for analysis of 47,000 transcripts. After data normalization, genes with uniformly low expression (detection P value of >0.05) were removed from consideration, leaving 23,386 probe sets for analysis using Significance Analysis of Microarrays (SAM 2.23) software (11). Relevant variables for the SAM analysis were as follows: two-class, unpaired data (wild-type and knockout prostates), imputation engine, 10-nearest neighbor; number of permutations, 200; RNG seed, 1234567; Δ, 1.345; fold change, 1.7; and false discovery rate, <1%. Microarray data are available at the GEO repository with accession numbers GSE9182 (Illumina) and GSE9207 (Affymetrix).

Bioinformatic analysis of transcription factor binding sites. For the in vivo position weight matrix (PWM) for HOXC6, we used THEME software (12) using the homeobox (PF00046) matrix as an initial hypothesis and the 21 confirmed peaks. Co-occurring transcription factor binding sites (TFBS) were identified using CONFAC software,5

which identifies TFBS that are conserved between mouse and human sequences (13). A total of 515 promoter fragments that were bound by HOXC6 in LNCaP cells were analyzed together with 10 sets of 515 promoter sequences that were not bound by HOXC6 in LNCaP cells. Mann-Whitney tests were used to identify TFBS that were significantly enriched in the HOXC6-bound fragments compared with the 10 random sets of fragments.

Quantitative real-time PCR. Total RNA was extracted from and treated as above. cDNA was synthesized by reverse transcription (Superscript III; Invitrogen), and transcript levels were quantitated by quantitative real-time PCR (QRT-PCR) using SYBR Green (Molecular Probes). Cycle threshold differences were determined using an I-cycler (Bio-Rad) relative to 18s to control for total starting RNA levels in each prostate sample. All primers used in this analysis are listed in Supplementary Table S2.

Apoptosis assay. LNCaP cells (1 × 104) were plated in a 96-well dish and allowed to adhere overnight. Cells were transiently transfected with pcDNA-HOXC6 expression vector or pcDNA empty vector using DMRI-C (Invitrogen) per the manufacturer's directions. Twenty-four hours posttransfection, cells were treated with 25 ng/mL recombinant BMP7 (Leinco Technologies) or PBS. At 72 h posttreatment, the caspase-3/7 activity was measured using the Caspase 3/7 kit (Promega) and luminescence was read using a Bio-Tek 96-well plate reader luminometer.

Proliferation assay. LNCaP cells (1 × 104) were plated in a 96-well dish and allowed to adhere overnight. Cells were transiently transfected with pcDNA-HOXC6 expression vector or pcDNA empty vector using DMRI-C (Invitrogen) per the manufacturer's directions. Twenty-four hours posttransfection, cells were treated with 25 ng/mL recombinant BMP7 (Leinco Technologies), 5 μmol/L imatinib (a generous gift from Dr. Jing Chen, Emory University, Atlanta, GA) or vehicle. At 72 h posttreatment, 20 μL of 5 μg/mL 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide (MTT; Sigma) was added to each well, incubated for 2 h at 37°C, supernatant removed, and precipitate was dissolved in DMSO and absorbance read at 570 nm.

HOXC6-HA associates with the IGFBP3 and TARP promoters in LNCaP cells. Although HOXC6 expression is up-regulated in human prostate cancers, little is known about the role it plays in disease progression (35). To understand the molecular role that HOXC6 plays in prostate cancer progression, it is essential to identify its direct downstream targets. Because antibodies able to specifically immunoprecipitate HOXC6 are not available, we created HOXC6-HA expressing stable LNCaP cells that could be used in ChIP experiments. A bicistronic HOXC6-HA-IRES-eYFP construct under an ubiquitin promoter was used in a lentiviral infection system to create LNCaP cells that stably express HOXC6-HA (Fig. 1A). The ubiquitin promoter allows for low, physiologic levels of constitutive HOXC6-HA expression. The infected cells were isolated from the noninfected cells using flow cytometry for YFP expression (Fig. 1B). In addition, a control LNCaP cell line that stably expressed YFP alone was created (data not shown). The resulting stable cell lines were tested for expression by immunoprecipitation with monoclonal anti-HA 12CA5 antibodies followed by immunoblot with polyclonal anti-HOXC6 (Fig. 1C).

Figure 1.

Creation of stable cells expressing HOXC6-HA/YFP or YFP only. A, transient virus-producing cell lines were created using a lentiviral vector modified with a bicistronic HOXC6-HA-IRES-eYFP expression cassette or IRES-eYFP alone. B, YFP-based fluorescence-activated cell sorting isolated near-pure populations of cells expressing the HOXC6-HA/YFP transgene. C, anti-HA antibody specifically immunoprecipitates HOXC6-HA. Lysates from LNCaP cells stably expressing HOXC6-HA/eYFP or eYFP alone were immunoprecipitated with anti-HA or anti-IgG and immunoblotted with HOXC6 antisera. D, ChIP using anti-IgG or anti-HA were performed with stable LNCaP cell lines expressing HOXC6-HA/eYFP (HOXC6-HA) or eYFP alone (eYFP) followed by PCR using primers to the candidate target promoter regions of TARP (two sets) and IGFBP3.

Figure 1.

Creation of stable cells expressing HOXC6-HA/YFP or YFP only. A, transient virus-producing cell lines were created using a lentiviral vector modified with a bicistronic HOXC6-HA-IRES-eYFP expression cassette or IRES-eYFP alone. B, YFP-based fluorescence-activated cell sorting isolated near-pure populations of cells expressing the HOXC6-HA/YFP transgene. C, anti-HA antibody specifically immunoprecipitates HOXC6-HA. Lysates from LNCaP cells stably expressing HOXC6-HA/eYFP or eYFP alone were immunoprecipitated with anti-HA or anti-IgG and immunoblotted with HOXC6 antisera. D, ChIP using anti-IgG or anti-HA were performed with stable LNCaP cell lines expressing HOXC6-HA/eYFP (HOXC6-HA) or eYFP alone (eYFP) followed by PCR using primers to the candidate target promoter regions of TARP (two sets) and IGFBP3.

Close modal

To ensure the ability of the stable cell lines to identify HOXC6-bound promoters, ChIP experiments were performed for promoter regions of candidate HOXC6 target genes. We previously identified IGFBP3 as a repression target of HOXC6 and TARP as an activation target of HOXC6 (5). Primers were designed to amplify regions in the promoters of IGFBP3 and TARP that contain consensus HOX family binding sites (Supplementary Table S1). Enrichment was seen in the anti-HA sample from the HOXC6-HA stable LNCaP cells compared with IgG control immunoprecipitations from these same cells. No enrichment was observed in anti-HA immunoprecipitations from the control YFP LNCaP cell line (Fig. 1D). Thus, IGFBP3 and TARP, whose transcript levels are altered upon perturbation of HOXC6 expression, are direct downstream targets of HOXC6 in LNCaP cells. Additionally, these results justify the use of stable LNCaP HOXC6-HA and YFP cell lines to globally identify HOXC6-bound promoters.

Global identification of novel gene promoters bound by HOXC6. To globally identify promoters bound by HOXC6 in LNCaP cells, genome-wide localization (ChIP-chip) was performed using the HOXC6-HA-LNCaP and YFP-LNCaP stable cell lines. Each ChIP sample was collected and analyzed in two independent replicates. Isolated DNA from each sample and total input DNA were amplified by ligation-mediated PCR. To identify the promoters bound by HOXC6, the anti-HA specific DNA and total input control DNA samples from the HOXC6-HA stable cells were hybridized to the Nimblegen 25k human promoter microarray set. The same analysis was performed using the anti-HA immunoprecipitation DNA sample and its total input control from the YFP alone stable cells to identify promoters that were enriched due to nonspecific DNA that coprecipitated with the 12CA5 anti-HA antibody, thus eliminating false positive results.

We used ChIPOTle analysis software (9) to identify 514 reproducible peaks, representing 468 different probe sets that were specific to the HOXC6-HA ChIP samples (Supplementary Table S3). To determine which annotated functional pathways were overrepresented by these 468 HOXC6 direct target genes, we performed an Ingenuity Pathways Analysis (IPA) to identify the canonical pathways and gene functions that are significantly enriched relative to random sets of genes. Interestingly, the 468 HOXC6 direct target genes were enriched for gene functions involved in apoptosis (P = 5.17E-5) and cellular development (P = 1.71E-4; Table 1). This finding supports a possible role in prostate cancer progression for these newly identified transcriptional targets of HOXC6.

Table 1.

GO biological processes significantly enriched in HOXC6 direct targets

CategoryProcess annotationSignificanceNo. of genes
Cancer Apoptosis of tumor cell lines 5.71E-05 31 
Cellular development Development of cells 1.71E-04 28 
Cell morphology Morphogenesis of cells 4.76E-04 16 
Cellular movement Cell movement 5.04E-04 36 
Tissue development Developmental process of tissue 3.34E-03 21 
CategoryProcess annotationSignificanceNo. of genes
Cancer Apoptosis of tumor cell lines 5.71E-05 31 
Cellular development Development of cells 1.71E-04 28 
Cell morphology Morphogenesis of cells 4.76E-04 16 
Cellular movement Cell movement 5.04E-04 36 
Tissue development Developmental process of tissue 3.34E-03 21 

To validate the ChIP-chip data set, PCR primers (Supplementary Table S1) were designed to 24 genomic regions representing HOXC6-HA–bound peaks and conventional ChIP-PCR was performed. The 24 peaks were chosen based on the function of their associated genes being relevant to cancer progression. Twenty-one of twenty-four peaks tested (87.5%) were reproducibly enriched by anti-HA 12CA5 ChIP in the HOXC6-HA cells compared with the YFP control (Fig. 2A). Thirteen of twenty-four genes, including BMP7, FGFR2, HOXC6, caspase recruitment domain family, member 8 (CARD8), caspase 1 (CASP1), and Runt-Related Transcription Factor 1 (RUNX1) were also confirmed as direct targets in RWPE-1 cells that stably express HOXC6-HA (Fig. 2A). To quantitate the enrichment of HOXC6-HA anti-HA DNA over YFP anti-HA DNA in LNCaP cells, QRT-PCR was performed for seven peaks (Fig. 2A). All seven promoters showed 2- to 7-fold enrichment over QRT-PCR from control YFP cells in LNCaP cells, and five of seven promoters were confirmed in RWPE-1 cells. These results verify the validity of the ChIP-chip analysis of HOXC6-HA bound promoters in LNCaP and RWPE-1 cells. Furthermore, they show that HOXC6 binds to its own promoter.

Figure 2.

A, ChIP-CHIP experiments were ∼86% (21 of 24) accurate as determined by conventional ChIP. PCR verification was performed with input DNA, water control, anti-HA ChIP DNA from LNCaP cells, and RWPE-1 cells stably expressing HOXC6-HA/YFP (HOXC6-HA) or YFP alone. *, genes whose expression was changed by HOXC6 expression in transient expression assays. QRT-PCR shows fold enrichment of HOXC6 promoter targets in HOXC6-HA–expressing cells versus YFP-expressing cells. The results are consistent with what is observed using traditional PCR. B, PDGFRA, CARD8, and CASP1 promoters were cloned into the pGL3-basic reporter vector and cotransfected in LNCaP cells with pDNA3.1 or pCDNA3.1-HOXC6. HOXC6 addition resulted in significant increase in luciferase reporter activity from the PDGFRA and CARD8 promoters and decrease from the CASP1 promoter compared with pCDNA3.1 vector alone. All fold changes were normalized to the HOXC6 effect on empty pGL3-basic vector. C, sequence for the 21 verified HOXC6-bound promoters were input for THEME analysis using the homeobox family PWM as an initial hypothesis. The most significant motif is pictured along with the cross-validation, significance, and β, similarity, scores. The HOXC6-predicted motif varies from the initial hypothesis by showing nonspecificity at the third position. D, CONFAC analysis of the 515 HOXC6-bound sequences compared with 10 random sets of 515 unbound sequences of similar size resulted in 18 enriched PWM. The significant PWM include homeodomain sites and forkhead protein sites.

Figure 2.

A, ChIP-CHIP experiments were ∼86% (21 of 24) accurate as determined by conventional ChIP. PCR verification was performed with input DNA, water control, anti-HA ChIP DNA from LNCaP cells, and RWPE-1 cells stably expressing HOXC6-HA/YFP (HOXC6-HA) or YFP alone. *, genes whose expression was changed by HOXC6 expression in transient expression assays. QRT-PCR shows fold enrichment of HOXC6 promoter targets in HOXC6-HA–expressing cells versus YFP-expressing cells. The results are consistent with what is observed using traditional PCR. B, PDGFRA, CARD8, and CASP1 promoters were cloned into the pGL3-basic reporter vector and cotransfected in LNCaP cells with pDNA3.1 or pCDNA3.1-HOXC6. HOXC6 addition resulted in significant increase in luciferase reporter activity from the PDGFRA and CARD8 promoters and decrease from the CASP1 promoter compared with pCDNA3.1 vector alone. All fold changes were normalized to the HOXC6 effect on empty pGL3-basic vector. C, sequence for the 21 verified HOXC6-bound promoters were input for THEME analysis using the homeobox family PWM as an initial hypothesis. The most significant motif is pictured along with the cross-validation, significance, and β, similarity, scores. The HOXC6-predicted motif varies from the initial hypothesis by showing nonspecificity at the third position. D, CONFAC analysis of the 515 HOXC6-bound sequences compared with 10 random sets of 515 unbound sequences of similar size resulted in 18 enriched PWM. The significant PWM include homeodomain sites and forkhead protein sites.

Close modal

HOXC6 activates expression from the PDGFRA and CARD8 promoters and represses expression from the CASP1 promoter. To explore the ability of HOXC6 to regulate transcription from target gene promoters identified by ChIP-chip, we cloned the promoters of three genes whose HOXC6-bound peak was relatively close to the start of transcription, PDGFRA (1,225 bp), CARD8 (736 bp), and CASP1 (1,298 bp) into a luciferase expression vector. The cloned promoters, at a minimum, contained the region bound by HOXC6 in LNCaP cells as well as the transcription start site of the target gene. Co-expression of HOXC6 with PDGFRA and CARD8 reporter plasmids increased luciferase activity, whereas coexpression of HOXC6 with CASP1 reporter plasmid decreased luciferase activity when compared with control transfections lacking HOXC6 (Fig. 2B). The direct interaction of HOXC6 protein with these promoters as well as the ability of HOXC6 to regulate reporter expression from these promoters indicate that HOXC6 is functionally capable of regulating gene expression through direct promoter interaction in living cells.

HOXC6 target promoters identify a novel in vivo PWM. The HOX family of transcription factors has little sequence specificity in vitro (14), which results in a PWM in the TRANSFAC database that is general to all HOX family members. To develop an in vivo PWM for HOXC6, we used THEME software (12). THEME uses a principled statistical method to test hypotheses about the binding site specificity of a transcription factor given a set of input sequences. We used the sequences for the 21 verified HOXC6-bound peaks as input for THEME analysis and the homeobox family PWM as an initial hypothesis.

To determine the significance of the THEME-discovered motif, an enrichment score, -log10 (P value), and a cross-validation error (CV), a measure of how often the discovered motif of THEME correctly predicts binding sites in a subset of input sequences, were calculated. CV values of <0.4 are considered significant. To determine the degree of similarity between the THEME-discovered motif and the initial hypothesis, a β score was calculated. The β score can range from 0 to 1, with higher values representing a motif more similar to the initial hypothesis. The most significant HOXC6-predicted PWM has an enrichment score of 1.66 and a CV value of 0.32. This motif is highly different from the initial hypothesis with a β score of 0.1 by showing nonspecificity at the third position (Fig. 2C). This analysis identified an in vivo derived PWM for HOXC6 and the first in vivo PWM available for any HOX family protein. Sequence analysis of the 468 promoters bound by HOXC6 determined that 263 (56%) contained matches to this motif.

PDGFRA, FGFR2, and RUNX1 expression are significantly increased by HOXC6 overexpression in LNCaP cells. To begin to address the functional role that HOXC6 plays when bound to its target promoters, we cross-referenced our previously generated HOXC6 gene expression microarray data (5) with our ChIP-chip data (Supplementary Table S3). We identified 50 genes with bound promoters that were significantly altered in our previously reported gene expression data, five of which were verified in ChIP-PCR experiments (Fig. 2A; Supplementary Table S4). Expression of FGFR2, PDGFRA, and RUNX1 were up-regulated, whereas p21-activated kinase 2 and RAB9A were down-regulated by HOXC6 overexpression in LNCaP cells. Up-regulation of PDGFRA upon overexpression of HOXC6 is consistent with our luciferase reporter assays demonstrating HOXC6 activation of the PDGFRA promoter region (Fig. 2B). Although CARD8 and CASP1 did not show gene expression changes upon overexpression of HOXC6, reporter assays using their promoters showed < 2-fold changes, suggesting a level of regulation by HOXC6 that may be below the level of detection of the microarray assay. Alternatively, in LNCaP cells, the genomic promoters may be constrained from regulation by HOXC6 due to chromatin structure or cell-specific co-factor differences.

To identify potential DNA-binding co-factors that could be interacting with HOXC6 on its target gene promoters, we performed CONFAC analysis. CONFAC identifies TFBS that are conserved between mouse and human sequences (13). We input the 515 promoter fragments that were bound by HOXC6 in LNCaP cells. We also generated 10 sets of 515 promoter sequences that were not bound by HOXC6 in LNCaP cells and performed a similar analysis. Mann-Whitney tests were used to identify TFBS that were significantly enriched in the HOXC6-bound fragments relative to at least 9 of the 10 random sets of sequences. This analysis identified 18 enriched PWMs (Fig. 2D). The PWM identified several homeodomain sites and forkhead sites, suggesting that HOXC6 may interact with one or more forkhead transcription factors.

To determine if these targets are perturbed in human prostate cancers, publicly available microarray expression data were downloaded6

(15) and integrated with our own public microarray expression data (16) to identify genes significantly changed in localized and metastatic prostate cancers (q value < 0.01). Of the 50 identified direct target genes of HOXC6, 33 have previously been shown to be altered in prostate cancer relative to benign prostate epithelium (Supplementary Table S4), including FGFR2, PDGFRA, Presenilin-1, and RUNX1. These data suggest that HOXC6 directly regulates the expression of these genes in LNCaP cells and prostate tumors.

Gene expression changes upon knockout of HOXC6 in mouse prostate overlap with HOXC6 direct target genes identified in LNCaP cells. To identify biologically relevant targets of Hoxc6 in vivo, we used mice homozygous for a targeted disruption of Hoxc6 DNA binding domain (10). Prostates were isolated from three Hoxc6+/+ and three Hoxc6−/− mice at ages 32 weeks, assayed for Hoxc6 RNA and protein expression (Fig. 3A and B), and subsequently used for whole-genome expression profiling. We identified 1,213 genes with differential gene expression between these groups using SAM analysis (Fig. 3C; Supplementary Table S3). Eight of these target genes were confirmed by QRT-PCR (Fig. 3D). Comparison of microarray data from Hoxc6+/+ relative to Hoxc6−/− prostates revealed that expression of Pdgfra, Fgfr2, Cd44, Hey1, Wif1, and Dkk3 was up-regulated, whereas Bmp7 expression was down-regulated.

Figure 3.

Three prostates from 32-wk-old Hoxc6−/− and Hoxc6+/+ mice were isolated and used for gene expression profiling. A, RT-PCR for Hoxc6 shows expression of Hoxc6 in each wild-type (WT) prostate and absence of expression in each knockout (MUT) prostate. B, Western blot analysis shows protein expression of HOXC6 in the wild-type mouse prostates and absence in the MUT mouse prostates. C, microarray analysis of each mouse prostate sample resulted in 1,213 genes that were significantly differentially expressed between wild-type and MUT mice. D, QRT-PCR-verified gene expression changes for eight genes in the wild-type versus knockout mouse prostates. Bmp7 shows decreased expression in the wild-type versus knockout mice prostates, whereas Pdgfra, Fgfr2, Gas2l3, Dkk3, Cd44, Hey1, and Wif1 show increased expression in wild-type versus knockout mouse prostates.

Figure 3.

Three prostates from 32-wk-old Hoxc6−/− and Hoxc6+/+ mice were isolated and used for gene expression profiling. A, RT-PCR for Hoxc6 shows expression of Hoxc6 in each wild-type (WT) prostate and absence of expression in each knockout (MUT) prostate. B, Western blot analysis shows protein expression of HOXC6 in the wild-type mouse prostates and absence in the MUT mouse prostates. C, microarray analysis of each mouse prostate sample resulted in 1,213 genes that were significantly differentially expressed between wild-type and MUT mice. D, QRT-PCR-verified gene expression changes for eight genes in the wild-type versus knockout mouse prostates. Bmp7 shows decreased expression in the wild-type versus knockout mice prostates, whereas Pdgfra, Fgfr2, Gas2l3, Dkk3, Cd44, Hey1, and Wif1 show increased expression in wild-type versus knockout mouse prostates.

Close modal

Thirty-one genes that showed mRNA expression changes in Hoxc6−/− compared with Hoxc6+/+ mouse prostates were also bound by HOXC6 in LNCaP cells (Table 2). Of these 31 genes, 21 were positively correlated with Hoxc6 expression and 10 were negatively correlated. We believe that this list of genes represent the most biologically relevant direct targets of HOXC6 in normal prostate tissue. Additionally, of those 31 direct targets of Hoxc6 in mouse prostate, 23 have previously been shown to be significantly altered in prostate cancers, including FGFR2 and PDGFRA (Table 2), making them strong candidates for regulation by HOXC6 in prostate cancer progression. Furthermore, 7 of these 31 genes, including FGFR2 and PDGFRA, were also affected by overexpression of HOXC6 in LNCaP cells (Table 2).

Table 2.

HOXC6 direct targets affected by loss of HOXC6 expression in normal mouse prostates

Hs symbolHs entrez IDMm entrez IDIllumina probe IDHOXC6−/− FCAffymetrix HOXC6 FCUp in PCaUp in MetsDown in PCaDown in Mets
CHRNA2 1135 110902 5050315 6.58      
AMELY 266 11832 2940722 3.49     
PDE4DIP 9659 83679 4730131 2.43 2.74 LV DTY   
BCDO2 83875 170752 4560279 2.25 4.57     
BMP7 655 12162 100360181 2.23   
SLC1A3 6507 20512 1500070 2.01      
SEPT7 989 235072 105570082 1.93 2.23   LW 
SIPA1L1 26037 217692 50292 1.92      
IRAK4 51135 266632 6130594 1.83      
RFXDC2 64864 319758 100110095 1.75      
ERRFI1 54206 74155 3450273 −1.73 −6.65    
FAM49A 81553 76820 102190458 −1.76    LV LY 
SENP8 123228 71599 2350184 −1.82      
DCN 1634 13179 106110022 −1.82     
SH3D19 152503 27059 6020520 −1.83 2.49    
ZBTB7C 201501 207259 2970609 −1.86      
PDGFRA 5156 18595 104670735 −1.87 2.32  MW 
PAK1 5058 18479 4540315 −1.88    LV  
FGFR2 2263 14183 6940377 −1.9 4.16   LMVW DLY 
SLC44A1 23446 100434 2510164 −1.96   
APBB1 322 11785 104070044 −2.05      
COL23A1 91522 237759 2030358 −2.05      
TNFAIP8 25816 106869 2760551 −2.08   
VCAM1 7412 22329 101580309 −2.11     
PLXDC1 57125 72324 540132 −2.2      
GAS2L3 283431 237436 105700397 −2.47  VY DY   
ZNF536 9745 243937 4920008 −2.69     TY 
C14orf37 145407 67419 5860195 −2.71      
FAM101A 144347 73121 5570040 −2.94     
CNTN1 1272 12805 3800204 −3.15    LMW DLY 
MEGF10 84466 70417 101990463 −3.71     
Hs symbolHs entrez IDMm entrez IDIllumina probe IDHOXC6−/− FCAffymetrix HOXC6 FCUp in PCaUp in MetsDown in PCaDown in Mets
CHRNA2 1135 110902 5050315 6.58      
AMELY 266 11832 2940722 3.49     
PDE4DIP 9659 83679 4730131 2.43 2.74 LV DTY   
BCDO2 83875 170752 4560279 2.25 4.57     
BMP7 655 12162 100360181 2.23   
SLC1A3 6507 20512 1500070 2.01      
SEPT7 989 235072 105570082 1.93 2.23   LW 
SIPA1L1 26037 217692 50292 1.92      
IRAK4 51135 266632 6130594 1.83      
RFXDC2 64864 319758 100110095 1.75      
ERRFI1 54206 74155 3450273 −1.73 −6.65    
FAM49A 81553 76820 102190458 −1.76    LV LY 
SENP8 123228 71599 2350184 −1.82      
DCN 1634 13179 106110022 −1.82     
SH3D19 152503 27059 6020520 −1.83 2.49    
ZBTB7C 201501 207259 2970609 −1.86      
PDGFRA 5156 18595 104670735 −1.87 2.32  MW 
PAK1 5058 18479 4540315 −1.88    LV  
FGFR2 2263 14183 6940377 −1.9 4.16   LMVW DLY 
SLC44A1 23446 100434 2510164 −1.96   
APBB1 322 11785 104070044 −2.05      
COL23A1 91522 237759 2030358 −2.05      
TNFAIP8 25816 106869 2760551 −2.08   
VCAM1 7412 22329 101580309 −2.11     
PLXDC1 57125 72324 540132 −2.2      
GAS2L3 283431 237436 105700397 −2.47  VY DY   
ZNF536 9745 243937 4920008 −2.69     TY 
C14orf37 145407 67419 5860195 −2.71      
FAM101A 144347 73121 5570040 −2.94     
CNTN1 1272 12805 3800204 −3.15    LMW DLY 
MEGF10 84466 70417 101990463 −3.71     

NOTE: Also shown are genes affected by overexpression of HOXC6 in LNCaP cells. HOXC6−/− FC represents the fold change of normal mouse prostates lacking HOXC6 relative to wild-type prostates, and thus, HOXC6 activation targets show negative fold changes for that column. Affymetrix HOXC6 FC indicates fold change from overexpression of HOXC6 in LNCaP cells, and thus, HOXC6 activation targets have positive fold changes in that column. Those genes with concordant effects in mouse prostate and LNCaP cells are shown in bold. Also indicated are microarray expression data from oncomine and our own previous data (16) indicating studies that showed changes in these genes in localized and metastatic prostate cancers.

Abbreviations: D, Dhanasekaran (41); L, LaPointe (43); M, Liu (16); S, Singh (46); T, Tomlins (42); V, Varambally (44); W, Welsh (45); Y, Yu (47).

Recombinant BMP7 decreases proliferation and increases apoptosis in LNCaP cells. Previously, we reported that overexpression of HOXC6 leads to increased proliferation and decreased apoptosis in LNCaP cells (5). Because BMP7 is a downstream repression target of HOXC6, we predicted that treatment of LNCaP cells with exogenous recombinant BMP7 could prevent the proliferative and antiapoptotic effects of HOXC6 overexpression. To test this prediction, we transiently transfected LNCaP cells with a HOXC6 expression vector or empty vector and added 25 ng/mL BMP7 or PBS vehicle to the medium. At 72 h posttreatment, we observed that recombinant BMP7 eliminates the reduction in apoptosis and increased proliferation caused by exogenous HOXC6 expression, as expected (Fig. 4A). These data suggest that one mechanism by which HOXC6 exerts effects on proliferation and apoptosis is through repression of BMP7.

Figure 4.

A, LNCaP cells were transiently transfected with empty vector of HOXC6 expression vector and treated with 25 ng/mL recombinant BMP7 or PBS vehicle. At 72 h posttreatment, cells were tested for proliferation by MTT assay or apoptosis by Caspase 3/7 activation assay. Transient expression of HOXC6 results in a significant increase in proliferation and decrease in apoptosis (*, P < 0.05 compared with vector + PBS) in LNCaP cells. Treatment of cells with recombinant BMP7 results in a significant increase in apoptosis and decrease in proliferation in LNCaP cells (*, P < 0.05 compared with Vector +PBS; **, P < 0.05 compared with HOXC6 + PBS). BMP7 blocks the antiapoptotic and proproliferative effect of increased HOXC6 levels. B, LNCaP cells were transiently transfected with HOXC6 expression vector or vector alone and treated with 5 μmol/L imatinib or DMSO. At 72 h, cells were tested for proliferation by MTT assay. Imatinib largely inhibited the proliferative effects of HOXC6 transient expression in LNCaP cells (*, P < 0.05 compared with vector + DMSO; **, P < 0.05 compared with vector + imatinib). C, IPA for genes demonstrating at least two lines of evidence for direct or indirect regulation by HOXC6. A subset of the direct and indirect target genes of HOXC6 creates a network that affects proliferation and survival through Akt, PI3K, and Caspase pathways. *, direct transcriptional regulation from this study. D, an overlapping subset of HOXC6 targets affects bone morphogenesis and branching morphogenesis. Red, up-regulated genes; green, down-regulated genes. The intensity of color reflects the level of expression change.

Figure 4.

A, LNCaP cells were transiently transfected with empty vector of HOXC6 expression vector and treated with 25 ng/mL recombinant BMP7 or PBS vehicle. At 72 h posttreatment, cells were tested for proliferation by MTT assay or apoptosis by Caspase 3/7 activation assay. Transient expression of HOXC6 results in a significant increase in proliferation and decrease in apoptosis (*, P < 0.05 compared with vector + PBS) in LNCaP cells. Treatment of cells with recombinant BMP7 results in a significant increase in apoptosis and decrease in proliferation in LNCaP cells (*, P < 0.05 compared with Vector +PBS; **, P < 0.05 compared with HOXC6 + PBS). BMP7 blocks the antiapoptotic and proproliferative effect of increased HOXC6 levels. B, LNCaP cells were transiently transfected with HOXC6 expression vector or vector alone and treated with 5 μmol/L imatinib or DMSO. At 72 h, cells were tested for proliferation by MTT assay. Imatinib largely inhibited the proliferative effects of HOXC6 transient expression in LNCaP cells (*, P < 0.05 compared with vector + DMSO; **, P < 0.05 compared with vector + imatinib). C, IPA for genes demonstrating at least two lines of evidence for direct or indirect regulation by HOXC6. A subset of the direct and indirect target genes of HOXC6 creates a network that affects proliferation and survival through Akt, PI3K, and Caspase pathways. *, direct transcriptional regulation from this study. D, an overlapping subset of HOXC6 targets affects bone morphogenesis and branching morphogenesis. Red, up-regulated genes; green, down-regulated genes. The intensity of color reflects the level of expression change.

Close modal

PDGFRA inhibition blocks the HOXC6-induced increase in prostate cancer cell proliferation. Imatinib is a potent inhibitor of cell proliferation through inhibition of PDGFRA activity (17). Because PDGFRA is a downstream activation target of HOXC6, we predicted that overexpression of HOXC6 might overcome the antiproliferative effects of imatinib by increasing PDGFRA levels. To test this prediction, we treated LNCaP cells with 5 μmol/L imatinib with or without HOXC6 overexpression. After 72 h of imatinib treatment, LNCaP cells transfected with vector alone showed a significant decrease in proliferation (Fig. 4B). In contrast, LNCaP cells transiently transfected with a HOXC6 expression vector and treated for 72 h with 5 μmol/L imatinib showed no reduction in proliferation by imatinib treatment. Interestingly, imatinib treatment did block the proliferative effect of HOXC6 in LNCaP cells, suggesting that HOXC6 exerts its proliferative effect partially through PDGFRA (Fig. 4B).

HOXC6 directly and indirectly regulates multiple genes involved in cell proliferation, cell death, and tissue development. To better understand the role that HOXC6 plays in prostate cancer progression, we performed IPA on those genes with at least two lines of evidence supporting their transcriptional regulation by HOXC6 in prostate cells (Supplementary Table S3). Those genes with promoters bound by ChIP experiments and whose expression is perturbed in in vitro or in vivo expression profiling experiments were considered direct targets of HOXC6. Those genes whose promoters were not bound in the ChIP experiments but whose expression was changed in both in vitro and in vivo expression profiling experiments were considered indirect targets of HOXC6. The 115 direct and indirect targets of HOXC6 were significantly enriched for cellular growth (P = 5.3E-9), apoptosis (P = 1.7E-6), development (P = 3.5E-6), and cancer (P = 8.3E-6) biological functions (Supplementary Table S5).

Some HOXC6 activation targets such as FGFR2 and CD44 have crucial roles in promoting prostate ductal branching morphogenesis (18, 19). To determine if prostate ductal branching was also reduced in Hoxc6−/− mice, we microdissected mature prostates from Hoxc6−/− and Hoxc6+/+ mice at 12 weeks. Although we did not observe a statistically significant quantitative difference in the number of prostate branches between these two groups of mice, we did observe that prostates from Hoxc6−/− exhibited a variable phenotype of fewer and thicker ducts with some aborted branches (Supplementary Fig. S1). The reduced penetrance of this phenotype may be due to compensatory activity from other HOX family members such as Hoxb6 (Supplementary Table S3), which is common among HOX transcription factors.

Here we have developed a comprehensive transcriptional network of genes under control of the developmental transcription factor HOXC6 in the mature prostate gland. Using a variety of complementary approaches, we have identified several direct targets of HOXC6, as well as a number of indirect targets that are affected by loss of HOXC6 in vivo. A major finding of this study is that several key developmental ligands and receptors are under direct regulation of HOXC6, including IGFBP3, BMP7, RUNX1, FGFR2, and PDGFRA. Increases in HOXC6 expression results in repression of IGFBP3 and BMP7 and activation of RUNX1, FGFR2, and PDGFRA.

IPA analysis of the direct and indirect targets of HOXC6 identified several pathways affected by changes in HOXC6 target genes expression level. The PI3K/Akt proproliferative pathway is regulated by the HOXC6 direct targets BMP7 (20), IGFBP3 (21), and PDGFRA (Fig. 4C; ref. 22). This suggests that one way HOXC6 may exert its prosurvival function is through the PI3K/Akt pathway. Moreover, our functional analysis of BMP7 and PDGFRA indicate that these genes are key effectors that mediate part of the proliferative and antiapoptotic effects of HOXC6 in LNCaP cells.

Because HOXC6 expression increases with tumor aggressiveness and metastasis in human prostate cancers, one would expect that the downstream targets of HOXC6 would favor proliferation and inhibit differentiation. Targets such as IGFBP3, PDGFRA, RUNX1, and BMP7 do support this aspect of the role of HOXC6. For example, overexpression of RUNX1 accelerates cell cycle progression in a myeloid progenitor cell line (23). Our identification of these ligands as direct targets of HOXC6 suggests that their expression changes may work together to regulate the proliferative and antiapoptotic effects of HOXC6.

Interestingly, there were also a number of tumor suppressor targets of HOXC6 activation that are believed to inhibit proliferation and promote differentiation. The most definitive of these targets is FGFR2 (18, 24), which is directly activated by HOXC6 and whose expression is reported most often to be lost in prostate cancers (25). Other indirect tumor suppressive genes that were activated by HOXC6 included CD44 and four inhibitors of Wnt signaling: WIF, DKK3, SFRP1, and SFRP2. Importantly, although HOXC6 activates expression of these genes in normal mouse prostates, all six of these genes have been previously shown to be silenced by hypermethylation in tumors and cancer cell lines (2631). Thus, epigenetic hypermethylation of these tumor suppressor genes via alternative mechanisms likely prevents their activation when HOXC6 is overexpressed in prostate cancers.

Many of the direct targets of HOXC6 have also been shown to play roles in bone morphogenesis (Fig. 4D), and metastatic prostate cancer is characterized by a strong predisposition for metastasis to the bone. BMP7 is a direct repression target of HOXC6 that has osteogenic properties that can produce ectopic bone formation (32). Importantly, BMP7 has recently been shown to inhibit growth of prostate cancer cells in mouse bone, suggesting it plays a role in inhibiting prostate cancer bone metastases (33). Furthermore, BMP7 represses Notch signaling and Hey1 during prostate development (20). Thus, our results support a model in which HOXC6 directly represses BMP7, indirectly activating the Notch pathway and Hey1 expression.

HOXC6 direct target PDGFRA and indirect target CD44 have also been shown to regulate bone morphogenesis. PDGFRA is a direct HOXC6 activation target that is essential for normal neural crest, vertebrae, and rib development (34) and has been shown to be a significant predictor of tumor progression in osteosarcomas (35). PDGFRA expression is increased in prostate cancer bone metastases (36) and correlates with propensity of prostate cancer cell lines to metastasize to the bone (37). The HOXC6 indirect activation target CD44 is critical for adhesion of prostate cancer cells to bone marrow endothelial cells (38). Moreover, CD44+ prostate cancer cells are enriched for metastatic and tumorigenic progenitor cells (39). Thus, HOXC6 seems to affect multiple genes that influence metastasis of prostate cancers to the bone microenvironment. HOXC6 may promote prostate tumor progression and metastasis to the bone through repression of BMP7 and activation of PDGFRA, CD44, and Hey1.

Another method by which HOXC6 may be exerting its prosurvival phenotype is through prostate ductal branching morphogenesis. HOXC6 activation targets FGFR2 (18) and CD44 (19) have shown crucial roles in promoting prostate ductal branching morphogenesis (Fig. 4D). The HOXC6 repression target BMP7 (Fig. 4D) has been shown to inhibit branching morphogenesis and Notch signaling in the prostate (20). Notch signaling is crucial for branching morphogenesis during prostate development (40). Hoxc6−/− mice have reduced thoracic mammary gland branching and ductal development in newborn pups (10). Thus, HOXC6 may promote branching morphogenesis by direct repression of BMP7, direct induction of FGFR2, and indirect induction of CD44 and Notch signaling.

Taken together, these data suggest a model in which HOXC6 is a critical regulator of normal prostate development, directly controlling expression of BMP7, FGFR2, and PDGFRA, indirectly inhibiting Wnt signals and activating Notch signals. However, in prostate cancers, epigenetic silencing of Wnt-suppressing target genes enables HOXC6 to activate Notch signals without interfering with Wnt signaling. Our data show that the transcriptional network regulated by HOXC6 plays crucial roles in proliferation, survival, and metastasis of prostate cancer cells. Further studies will be required to determine whether HOXC6 represents a potential therapeutic target in prostate cancer.

Note: Supplementary data for this article are available at Cancer Research Online (http://cancerres.aacrjournals.org/).

Grant support: NCI R01 CA106826 (C.S. Moreno) and Department of Defense Congressionally Directed Medical Research Programs Prostate Cancer Postdoctoral Training Fellowship PC060114 (C.D. McCabe).

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

We thank Dr. Maja Kodani in the Winship Cancer Institute Microarray Core Facility, the Robert Karaffa at the Emory University FlowCore, Dr. Hinh Ly for the IRES-eYFP lentiviral vector, Dr. Jing Chen for imatinib, and the Winship Cancer Institute Transgenic Core Facility for help with prostate microdissections.

1
Veraksa A, Del Campo M, McGinnis W. Developmental patterning genes and their conserved functions: from model organisms to humans.
Mol Genet Metab
2000
;
69
:
85
–100.
2
Abate-Shen C. Deregulated homeobox gene expression in cancer: cause or consequence?
Nat Rev Cancer
2002
;
2
:
777
–85.
3
Miller GJ, Miller HL, van Bokhoven A, et al. Aberrant HOXC expression accompanies the malignant phenotype in human prostate.
Cancer Res
2003
;
63
:
5879
–88.
4
Bibikova M, Chudin E, Arsanjani A, et al. Expression signatures that correlated with Gleason score and relapse in prostate cancer.
Genomics
2007
;
89
:
666
–72.
5
Ramachandran S, Liu P, Young AN, et al. Loss of HOXC6 expression induces apoptosis in prostate cancer cells.
Oncogene
2005
;
24
:
188
–98.
6
Zhao Y, Potter SS. Functional comparison of the Hoxa 4, Hoxa 10, and Hoxa 11 homeoboxes.
Dev Biol
2002
;
244
:
21
–36.
7
Odom DT, Zizlsperger N, Gordon DB, et al. Control of pancreas and liver gene expression by HNF transcription factors.
Science
2004
;
303
:
1378
–81.
8
Ren B, Robert F, Wyrick JJ, et al. Genome-wide location and function of DNA binding proteins.
Science
2000
;
290
:
2306
–9.
9
Buck MJ, Nobel AB, Lieb JD. ChIPOTle: a user-friendly tool for the analysis of ChIP-chip data.
Genome Biol
2005
;
6
:
R97
.
10
Garcia-Gasca A, Spyropoulos DD. Differential mammary morphogenesis along the anteroposterior axis in Hoxc6 gene targeted mice.
Dev Dyn
2000
;
219
:
261
–76.
11
Tusher VG, Tibshirani R, Chu G. Significance analysis of microarrays applied to the ionizing radiation response.
Proc Natl Acad Sci U S A
2001
;
98
:
5116
–21.
12
Macisaac KD, Gordon DB, Nekludova L, et al. A hypothesis-based approach for identifying the binding specificity of regulatory proteins from chromatin immunoprecipitation data.
Bioinformatics
2006
;
22
:
423
–9.
13
Karanam S, Moreno CS. CONFAC: automated application of comparative genomic promoter analysis to DNA microarray data sets.
Nucleic Acids Res
2004
;
32
:
W475
–84.
14
Shang Z, Ebright YW, Iler N, et al. DNA affinity cleaving analysis of homeodomain-DNA interaction: identification of homeodomain consensus sites in genomic DNA.
Proc Natl Acad Sci U S A
1994
;
91
:
118
–22.
15
Rhodes DR, Yu J, Shanker K, et al. ONCOMINE: a cancer microarray database and integrated data-mining platform.
Neoplasia
2004
;
6
:
1
–6.
16
Liu P, Ramachandran S, Ali-Seyed M, et al. Sex-determining region Y box 4 is a transforming oncogene in human prostate cancer cells.
Cancer Res
2006
;
46
:
4011
–9.
17
Matei D, Chang DD, Jeng MH. Imatinib mesylate (Gleevec) inhibits ovarian cancer cell growth through a mechanism dependent on platelet-derived growth factor receptor α and Akt inactivation.
Clin Cancer Res
2004
;
10
:
681
–90.
18
Lin Y, Liu G, Zhang Y, et al. Fibroblast growth factor receptor 2 tyrosine kinase is required for prostatic morphogenesis and the acquisition of strict androgen dependency for adult tissue homeostasis.
Development
2007
;
134
:
723
–34.
19
Gakunga P, Frost G, Shuster S, Cunha G, Formby B, Stern R. Hyaluronan is a prerequisite for ductal branching morphogenesis.
Development
1997
;
124
:
3987
–97.
20
Grishina IB, Kim SY, Ferrara C, Makarenkova HP, Walden PD. BMP7 inhibits branching morphogenesis in the prostate gland and interferes with Notch signaling.
Dev Biol
2005
;
288
:
334
–47.
21
Ricort JM, Binoux M. Insulin-like growth factor-binding protein-3 activates a phosphotyrosine phosphatase. Effects on the insulin-like growth factor signaling pathway.
J Biol Chem
2002
;
277
:
19448
–54.
22
van der Geer P, Hunter T, Lindberg RA. Receptor protein-tyrosine kinases and their signal transduction pathways.
Annu Rev Cell Biol
1994
;
10
:
251
–337.
23
Strom DK, Nip J, Westendorf JJ, et al. Expression of the AML-1 oncogene shortens the G(1) phase of the cell cycle.
J Biol Chem
2000
;
275
:
3438
–45.
24
Matsubara A, Kan M, Feng S, McKeehan WL. Inhibition of growth of malignant rat prostate tumor cells by restoration of fibroblast growth factor receptor 2.
Cancer Res
1998
;
58
:
1509
–14.
25
Naimi B, Latil A, Fournier G, Mangin P, Cussenot O, Berthon P. Down-regulation of (IIIb) and (IIIc) isoforms of fibroblast growth factor receptor 2 (FGFR2) is associated with malignant progression in human prostate.
Prostate
2002
;
52
:
245
–52.
26
Lodygin D, Epanchintsev A, Menssen A, Diebold J, Hermeking H. Functional epigenomics identifies genes frequently silenced in prostate cancer.
Cancer Res
2005
;
65
:
4218
–27.
27
Wissmann C, Wild PJ, Kaiser S, et al. WIF1, a component of the Wnt pathway, is down-regulated in prostate, breast, lung, and bladder cancer.
J Pathol
2003
;
201
:
204
–12.
28
Park S, Kim JH, Jang JH. Aberrant hypermethylation of the FGFR2 gene in human gastric cancer cell lines.
Biochem Biophys Res Commun
2007
;
357
:
1011
–5.
29
Zhu X, Lee K, Asa SL, Ezzat S. Epigenetic silencing through DNA and histone methylation of fibroblast growth factor receptor 2 in neoplastic pituitary cells.
Am J Pathol
2007
;
170
:
1618
–28.
30
Lou W, Krill D, Dhir R, et al. Methylation of the CD44 metastasis suppressor gene in human prostate cancer.
Cancer Res
1999
;
59
:
2329
–31.
31
Zou H, Molina JR, Harrington JJ, et al. Aberrant methylation of secreted frizzled-related protein genes in esophageal adenocarcinoma and Barrett's esophagus.
Int J Cancer
2005
;
116
:
584
–91.
32
Franceschi RT, Wang D, Krebsbach PH, Rutherford RB. Gene therapy for bone formation: in vitro and in vivo osteogenic activity of an adenovirus expressing BMP7.
J Cell Biochem
2000
;
78
:
476
–86.
33
Buijs JT, Rentsch CA, van der Horst G, et al. BMP7, a putative regulator of epithelial homeostasis in the human prostate, is a potent inhibitor of prostate cancer bone metastasis in vivo.
Am J Pathol
2007
;
171
:
1047
–57.
34
Soriano P. The PDGF α receptor is required for neural crest cell development and for normal patterning of the somites.
Development
1997
;
124
:
2691
–700.
35
Sulzbacher I, Birner P, Trieb K, Traxler M, Lang S, Chott A. Expression of platelet-derived growth factor-AA is associated with tumor progression in osteosarcoma.
Mod Pathol
2003
;
16
:
66
–71.
36
Chott A, Sun Z, Morganstern D, et al. Tyrosine kinases expressed in vivo by human prostate cancer bone marrow metastases and loss of the type 1 insulin-like growth factor receptor.
Am J Pathol
1999
;
155
:
1271
–9.
37
Dolloff NG, Shulby SS, Nelson AV, et al. Bone-metastatic potential of human prostate cancer cells correlates with Akt/PKB activation by α platelet-derived growth factor receptor.
Oncogene
2005
;
24
:
6848
–54.
38
Draffin JE, McFarlane S, Hill A, Johnston PG, Waugh DJ. CD44 potentiates the adherence of metastatic prostate and breast cancer cells to bone marrow endothelial cells.
Cancer Res
2004
;
64
:
5702
–11.
39
Patrawala L, Calhoun T, Schneider-Broussard R, et al. Highly purified CD44+ prostate cancer cells from xenograft human tumors are enriched in tumorigenic and metastatic progenitor cells.
Oncogene
2006
;
25
:
1696
–708.
40
Wang XD, Leow CC, Zha J, et al. Notch signaling is required for normal prostatic epithelial cell proliferation and differentiation.
Dev Biol
2006
;
290
:
66
–80.
41
Dhanasekaran SM, Barrette TR, Ghosh D, et al. Delineation of prognostic biomarkers in prostate cancer.
Nature
2001
;
412
:
822
–6.
42
Tomlins SA, Mehra R, Rhodes DR, et al. Integrative molecular concept modeling of prostate cancer progression.
Nat Genet
2007
;
39
:
41
–51.
43
Lapointe J, Li C, Higgins JP, et al. Gene expression profiling identifies clinically relevant subtypes of prostate cancer.
Proc Natl Acad Sci U S A
2004
;
101
:
811
–6.
44
Varambally S, Yu J, Laxman B, et al. Integrative genomic and proteomic analysis of prostate cancer reveals signatures of metastatic progression.
Cancer Cell
2005
;
8
:
393
–406.
45
Welsh JB, Sapinoso LM, Su AI, et al. Analysis of gene expression identifies candidate markers and pharmacological targets in prostate cancer.
Cancer Res
2001
;
61
:
5974
–8.
46
Singh D, Febbo PG, Ross K, et al. Gene expression correlates of clinical prostate cancer behavior.
Cancer Cell
2002
;
1
:
203
–9.
47
Yu YP, Landsittel D, Jing L, et al. Gene expression alterations in prostate cancer predicting tumor aggression and preceding development of malignancy.
J Clin Oncol
2004
;
22
:
2790
–9.