Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Oncogene. Author manuscript; available in PMC 2009 May 4.
Published in final edited form as:
PMCID: PMC2676849

Androgen-induced programs for prostate epithelial growth and invasion arise in embryogenesis and are reactivated in cancer


Cancer cells differentiate along specific lineages that largely determine their clinical and biologic behavior. Distinct cancer phenotypes from different cells and organs likely result from unique gene expression repertoires established in the embryo and maintained after malignant transformation. We used comprehensive gene expression analysis to examine this concept in the prostate, an organ with a tractable developmental program and a high propensity for cancer. We focused on gene expression in the murine prostate rudiment at three time points during the first 48 h of exposure to androgen, which initiates proliferation and invasion of prostate epithelial buds into surrounding urogenital sinus mesenchyme. Here, we show that androgen exposure regulates genes previously implicated in prostate carcinogenesis comprising pathways for the phosphatase and tensin homolog (PTEN), fibroblast growth factor (FGF)/mitogen-activated protein kinase (MAPK), and Wnt signaling along with cellular programs regulating such ‘hallmarks’ of cancer as angiogenesis, apoptosis, migration and proliferation. We found statistically significant evidence for novel androgeninduced gene regulation events that establish and/or maintain prostate cell fate. These include modulation of gene expression through microRNAs, expression of specific transcription factors, and regulation of their predicted targets. By querying public gene expression databases from other tissues, we found that rather than generally characterizing androgen exposure or epithelial budding, the early prostate development program more closely resembles the program for human prostate cancer. Most importantly, early androgen-regulated genes and functional themes associated with prostate development were highly enriched in contrasts between increasingly lethal forms of prostate cancer, confirming a ‘reactivation’ of embryonic pathways for proliferation and invasion in prostate cancer progression. Among the genes with the most significant links to the development and cancer, we highlight coordinate induction of the transcription factor Sox9 and suppression of the proapoptotic phospholipid-binding protein Annexin A1 that link early prostate development to early prostate carcinogenesis. These results credential early prostate development as a reliable and valid model system for the investigation of genes and pathways that drive prostate cancer.

Keywords: prostate organogenesis, androgen signaling, prostate cancer, microarray


The discovery of proto-oncogenes in the 1970’s (Stehelin et al., 1976) and tumor suppressor genes a decade later (Friend et al., 1986) launched an effort to identify a common genetic basis for all cancers. This approach has yielded significant insights into common molecular machinery regulating tumor initiation and growth (Hanahan and Weinberg, 2000), but fails to account for dramatic differences in the behaviors of tumors arising in different sites. As revealed by recent genomic approaches, the site of origin is the dominant influence on gene expression in cancers (Ramaswamy et al., 2001) and specifies particular oncogenic mutations (Garraway and Sellers, 2006). Tissue-specific behaviors of cancers almost certainly reflect lineage-specific epigenetic programs that operate during embryogenesis in cells and tissues from which cancers arise. Indeed, reawakening of embryonic programs has long been posited for cancer (Bailey and Cushing, 1925), and could underlie tissue specific modes of regulating critical aspects of the malignant phenotype such as survival, angiogenesis, invasion and migration. The advent of comprehensive genomic profiling techniques now permits this hypothesis to be tested, and to link cancer to embryogenesis through cellular pathways that define specific lineages and organs (Garraway and Sellers, 2006). Identification of these pathways and their functions should identify targets for more specific and less toxic cancer therapies.

The prostate gland represents a prime target for such an analysis as the genetic causes of prostate cancer are poorly understood and because there are excellent models of prostate development and homeostasis. While other cancers have canonical cancer initiating mutations (for example, K-ras for pancreatic carcinoma and Adenomatosis Polyposis Coli for colon carcinoma), the driving force for prostate cancer initiation in humans remains uncertain. Answers may lie in cellular programs activated by androgen receptor (AR) signaling, which strictly controls prostate epithelial cell fate in embryogenesis and in cancer.

Upon binding androgen, AR signals through genomic and nongenomic modes (Manolagas et al., 2002). Genomic responses entail nuclear translocation of liganded AR and activation of transcription at regulatory regions containing AR-binding sites. Nongenomic signaling, in contrast, comprises protein-protein interactions in the cytoplasm and is exceedingly rapid, with measurable responses within minutes of androgen exposure (Kousteni et al., 2001). AR signaling, through genomic and/or nongenomic routes, is necessary and sufficient for prostate organogenesis. In response to circulating testosterone and its local conversion to the more potent androgen, dihydrotestosterone, AR signaling induces epithelial-mesenchymal interactions in the urogenital sinus (UGS), an embryonic rudiment present in both sexes. These interactions lead to outgrowth of buds from the UGS epithelium (UGE) that proliferate and invade surrounding urogenital sinus mesenchyme. AR expression has been described in both UGE and urogenital sinus mesenchyme (Drews et al., 2001; B Simons, EM Schaeffer and DM Berman, unpublished observations), although epithelial AR is dispensable for induction of prostate development (Cunha et al., 2004).


Hormonal manipulation of mouse prostate development

In mouse, epithelial expression of the androgen-sensitive homeobox family member Nkx3. 1 begins by embryonic day 16 (e16), and is the earliest reported marker of prostate development (Bhatia-Gaur et al., 1999). We comprehensively profiled gene expression in the UGS between e16 and e17.5, when prostate buds first emerge. The developmental fate of the UGS is bipotential in both sexes and depends soley on androgens to drive prostate formation in males (Jost, 1968). The absence of androgens leads to vaginal/urethral formation in females (Wilson et al., 1980). Taking advantage of this dynamic, we induced prostate development in the androgen-naïve yet androgen-sensitive female UGS with precisely timed 6 and 12 h intrauterine exposure to pharmacological levels of dihydrotestosterone (Figure 1a). For the latest time point, approximating 48 h of physiologic androgen exposure, we compared UGS gene expression in unmanipulated male and female e17.5 littermates. In all cases, we isolated high quality RNA, prepared labeled cRNA pools and competitively hybridized androgen-exposed pools against androgennaïve pools on Agilent 44 000 probe whole genome mouse microarrays. High quality expression data was confirmed as described in Supplementary (S) methods.

Figure 1
Flowchart of data acquisition and analysis. (a) Schematic of early prostate development. The embryonic prostate rudiment, the urogenital sinus (UGS). Mesenchyme (light blue) surrounds epithelium (darker green). In the mouse, prostate-specific ...

Distinct and prostate-specific phases of gene expression in early prostate development

The genomic response to androgen varied dramatically with the duration of exposure to the hormone. After 6, 12 and 48 h of androgen exposure, respectively, 693, 177 and 829 genes were differentially expressed (adj. P<0.03) (Figure 2a). Complete lists are provided in the Supplementary material (Supplementary Tables 1-3). Pairwise comparisons between any two time points revealed an overlap between 10-42% of differentially expressed genes. Despite some diversity in the most highly differentially expressed genes across different time points, we successfully validated the androgen-mediated developmental program using a variety of methods. Starting with the latest time point, reverse transcriptase-PCR analysis of 20 transcripts confirmed differential expression that was concordant with the array results (Supplementary Figure 1). As expected, Nkx3.1 was highly induced (adj. P<10-4) among the transcripts confirmed by reverse transcriptase-PCR. We further validated these results with separate biologic replicates and with RNA from embryos treated with pharmacologic AR blockade (Supplementary Figure 1). Using the 48 h time point as a standard, we performed ‘correspondence at the top (CAT)’ analysis (Irizarry et al., 2005) to compare identities of top ranking (most differentially expressed) genes in lists from different time points. CAT analysis showed that both pharmacologic and physiologic androgen exposure induced highly related gene expression programs in the UGS, (Figure 2d). This was particularly true when comparing gene expression at 12 h and 48 h; gene lists from these two time points showed 40% identity for top ranking genes. By connecting physiologic and pharmacologic androgen responses in the UGS, these results validated our experimental approach.

Figure 2
Androgen-induced gene expression in early prostate development is dynamic and organ specific. (a) Distinct and overlapping genomic responses to androgen at 6, 12 or 48 h of exposure (see text). Values represent differential expression at adj. P<0.05. ...

Several genes were consistently androgen regulated; either induced (Aspn, Klf9, Synpr, Gadd45 g, Sox9, Adamdec1 and Tle1) or suppressed (Prrx1 and Inhba) at all three time points (Figure 2b). This cohort may represent particularly important genes in establishing and maintaining prostate lineage.

The genomic response to androgen is dynamic

The overall pattern of genomic response to androgen was suppression, mild induction and robust induction of gene expression at 6, 12, and 48 h, respectively (Figure 2c). To better understand this pattern of regulation, we analyzed functional annotations of differentially expressed genes (Supplementary Methods). At 6 h, an overwhelming majority (76%) were suppressed (Figure 2c), and most (58%) of the suppressed genes contained one or more predicted binding site for miRNAs, far exceeding the number expected by chance (P<1 × 10-5; χ2-test). Considering previously described roles for miRNAs in developmental biology, and specific functions in targeting complementary transcripts for destruction (Stefani and Slack, 2008), our results suggest a mechanism for rapid gene suppression in response to androgen. After longer exposure to androgen, the majority of genes were induced, with 71 (69%) and 403 (54%) genes showing increased expression at 12 and 48 h, respectively.

At equal statistical stringency (log odd>1.1, adjusted P<0.03), we found more differentially expressed genes at the 6 and 48 h time points than at 12 h. A portion of this effect stems from additional biologic replicates available for 6 and 48 h (Supplementary Methods). However, biological factors likely explain most of the variability in numbers of differentially expressed genes. These include expression levels passing through neutral as they change from repressed at 6 h to induced at 48 h, and an amplification effect of transcription factor expression. Accordingly, 50% of genes with enrichment at both 6 and 48 h were repressed at 6 h and induced at 48 h. Several transcription factors were induced at 12 h and their predicted target genes were differentially expressed at 48 h (Supplementary Tables 2, 12). Examples include induction of Sox9 followed by differential expression of 33 of its predicted targets (Supplementary Table 12), including fibroblast growth factor R3 (FgfR3), Hoxb9 and Tle1 (see below). Altogether, these data indicate that androgen exposure reorganizes the genomic repertoire of the UGS toward prostatic growth and differentiation-first by suppressing a bipotential gene expression program, and then by inducing transcription factors that initiate and carry out prostate-specific gene expression.

Prostate-specific responses to androgen in the embryo

To further define the specificity of expression programs in embryonic prostate, we used publicly available gene expression data to compare early prostate development to embryonic and androgen responsive gene expression in other tissues. As the prostate utilizes programs of branching morphogenesis during development, we employed CAT analysis to compare top ranking genes in prostate development with top ranking genes associated with embryonic branching in the lung (Lu et al., 2004). Similarly, we compared gene expression in the developing prostate and another androgen-responsive organ, the adult salivary gland (Treister et al., 2005). We observed distinct clustering of top ranking genes by organ (Figures 2d and e). Rather than generic programs for branching morphogenesis and androgen response, these analyses highlight a distinct developmentally and hormonally restricted lineage program governing prostate development.

Discrete and overlapping features of prostate development and prostate homeostasis

Prostate homeostasis in adults can be studied by hormonal manipulation of the mature gland and consists of three phases: active cell death upon androgen withdrawal (regression), castrate steady state and androgen induced growth (regeneration; English et al., 1987). Cross platform analysis of the 1536 and 2197 genes differentially expressed (adj. P<0.001) in development or regeneration (Wang et al., 2007b), respectively, revealed only a limited number of differentially expressed genes enriched in both processes (333 genes, 9%) (Supplementary Figure 2). Annotation analysis reported previously for prostate regeneration noted significant enrichment of adhesion, proliferation and metabolic themes (Wang et al., 2007b). We hypothesized that prostate organogenesis would involve a wider variety of themes that govern additional processes, such as lineage differentiation, extracellular matrix remodeling and migration.

Androgen-activated embryonic gene expression invokes the ‘hallmarks’ of cancer

By analysis of functional annotation, we identified significant themes in early prostate development (Table 1, Supplementary Tables 4-12). A self renewing population of stem cells has recently been described in the adult murine prostate (Lawson et al., 2007), and stem cell-associated genes were enriched (P<0.05) among androgen response genes (Supplementary Table 19). Sca1 expression, which in part characterize stem cells in the adult prostate (Xin et al., 2005) was highly induced at the 48 h time point and may signal the emergence of an androgen-regulated progenitor population at this stage. All three time points exhibited enrichment (adj. P<0.001) of AR targets as predicted by transcription factor-binding site analysis. We also identified significant (adj. P<0.05) enrichment of genes associated with ‘hallmark’ characteristic activities of cancer (Hanahan and Weinberg, 2000) including angiogenesis, apoptosis, migration, motility and proliferation (Table 1). Similarly, and despite the dominant role of post-translational modifications in their regulation, several extracellular and intracellular signaling pathways showed significant androgen-induced responses. These included oncogenic pathways, such as FGF/mitogen-activated protein kinase (MAPK), phosphatase and tensin (PTEN)-PI3K-mTOR (mammalian target of rapamaycin) and Wnt (Table 1). These findings indicate that cancers may reactivate cellular programs initiated in development, a concept we test more formally below.

Table 1
Functional annotation enrichment during early androgen induced prostate development

We also captured evidence of transcription factor networks acting downstream of the AR. Transcription factors including AR, Foxf2 and Sox9 were induced within 12 h of androgen exposure and associated with coordinate expression of their predicted target genes (Table 2, Supplementary Tables 11,12). This observation not only identifies new androgen-regulated activities of these transcription factors, but also cross validates the effects of androgens across different time points in our analysis. Altogether, this work, including detailed annotated gene lists (see Supplementary material), provides the first comprehensive map of the androgen response in prostate development and highlights androgen-mediated programs for regulated invasion, proliferation, survival and other cellular behaviors that become aberrantly regulated in cancer.

Table 2
Embryonic genes enriched in cancer transitions

Signaling in prostate development: FGF/MAPK

Secreted fibroblast growth factor (FGF) ligands signal through their receptors to activate intracellular pathways controlling proliferation. FGFs have multiple known roles in development and constitutive activation of the pathway in mouse prostate leads to cancer (Acevedo et al., 2007; Memarzadeh et al., 2007). In development, Fgf7 and Fgf10 have been intensively studied as candidate paracrine signals that facilitate epithelial budding in response to androgen (Donjacour et al., 2003). However, despite dramatic effects of these peptides on prostate ductal outgrowth, Fgf7 and Fgf10 demonstrate largely equivalent expression patterns in male and female UGS tissue. Thus, other factors must be required to induce prostate development, possibly by modulating responsiveness to Fgfs (Thomson, 2001; Donjacour et al., 2003). In the developing prostate, we identified sexually dimorphic expression of Fgf receptors and ligands as part of the MAPK pathway through which they signal (Table 1). Fgf18, which binds to FgfR3, (reviewed in (Mohammadi et al., 2005) is induced at 6 h and then is highly suppressed later in prostate development. Fgf12 was induced in males, although it reportedly does not bind FgfRs (Olsen et al., 2003). Two Fgf receptors showed increased expression in androgen-induced tissue; FgfR3, which can bind a wide range of Fgf ligands, and Fgfr2 which encodes the preferred receptor for Fgfs 7 and 10 (Mohammadi et al., 2005). These results implicate Fgf12 and Fgf18 in UGS differentiation and suggest that previously observed sexually dimorphic responses to Fgf7 and Fgf10 could stem from male-specific expression of FgfR2.

Signaling in prostate development: PTEN/PI3K/mTOR signaling

Aberrant signaling through PTEN and downstream effectors is highly implicated in prostate cancer, but its expression in early prostate development is uncharacterized. PTEN function suppresses a downstream signaling cascade whose effectors PI3 K, AKT (a proto-oncogene), mTOR and ribosomal S6 kinase (S6K) promote a wide variety of cellular functions, including growth, survival, migration and angiogenesis. PTEN loss is common in advanced prostate cancer (McMenamin et al., 1999; Wang et al., 2003) and PTEN loss is oncogenic in the mouse prostate, perhaps due to its promotion of epithelial cell growth and survival in the face of androgen withdrawal (Shen and Abate-Shen, 2007). Intriguingly, androgen rapidly induced several components of this pathway, including Pik3r, a PI3K regulatory subunit Rheb, a small GTPase Ras homolog required for mTOR activity (Long et al., 2005), and Rps6kb1, a subunit of the mTOR target, ribosomal protein S6 kinase. By the 48-h time point, PTEN transcripts were suppressed, along with the Pik3r2 regulatory subunit of PI3K. These results indicate dynamic androgen-responsive regulation of the Pten signaling pathway. Future studies in prostate development may therefore elucidate normal Pten/PI3K/mTOR regulation, its functional roles and more effective strategies for pathway measurement and manipulation.

Signaling in prostate development: the Wnt pathway

Extending an observation made previously by Serial Analysis of Gene Expression in the urogenital sinus (Zhang et al., 2006), analyses of functional annotation revealed Wnt signaling as one of the most androgen regulated pathways in early prostate development (adj. P<0.001), with regulation of multiple Wnt ligands and soluble inhibitors (Table 1). Wnt ligands are highly conserved secreted molecules that play critical but pleiotropic roles in cell-cell signaling during embryogenesis (Nusse, 2005). Studies of this pathway in prostate cancer show up- and downregulation of Wnt ligands, both stimulatory and inhibitory, with little or no agreement as to their net effect (Verras and Sun, 2006). Likewise, nearly all Wnts are expressed in the developing prostate (not shown) along with Wnt inhibitory factor (Wif-1), and secreted frizzled proteins, which can modulate UGS development in culture (Joesting et al., 2005; Table 1). Androgen-responsive expression was also observed for Wnt ligand coactivator, R-spondin1 (Rspo1; Nam et al., 2007), and several intracellular components of the pathway, including disheveled 1 (Dvl1), lymphoid enhancer-binding factor 1 (lef1) and transducin-like enhancer of split 1 (Tle1; Table 1). Of special interest, the Wnt target gene matrix metalloproteinase 7 (Mmp7; Crawford et al., 1999) has a well-established role in cancer invasion (Ii et al., 2006), and is induced during prostate budding (Table 1). These results suggest a general role for Wnt signal modulation in prostate development that may underlie the mechanics of epithelial invasion in prostate development and prostate cancer.

Prostate carcinogenesis and progression reactivates androgen-induced embryonic gene expression programs

Having demonstrated that androgen-responsive embryonic gene expression programs regulate processes that are fundamental to cancer, we examined the expression of embryonic prostate genes in human prostate cancer tissues. Like prostate organogenesis, initiation and progression of prostate cancer depends on AR signaling (Scher and Sawyers, 2005). Neoplastic transformation in the prostate is first recognized as high grade prostatic intraepithelial neoplasia (PIN), a noninvasive lesion that precedes invasive cancer (Bostwick and Qian, 2004). Invasive prostate cancer can be indolent, particularly in lower grade tumors (Gleason 6 or less). Higher grade tumors (Gleason 8-10) more often metastasize and cause death. Using public domain gene expression data from a large and well-characterized series (Lapointe et al., 2004), we assembled lists of genes that characterize transitions to more invasive lesions (that is, normal vs locally invasive cancers; lower grade vs higher grade locally invasive cancers and localized vs metastatic cancers). Comparing these lists to lists of differentially expressed genes at each of the developmental time points revealed significant (adj P<0.001) enrichment of developmental genes in invasive transitions, especially in more aggressive prostate cancers (Figure 3a, top). In contrast, gene expression in adult prostate regeneration showed little enrichment in cancer (Figure 3a, bottom). The latter result was somewhat surprising, given the previously reported enrichment of cell cycle, cell adhesion and metabolic genes in the prostate regeneration study (Wang et al., 2007b).

Figure 3
Embryonic gene expression in human prostate cancer. (a) Genes differentially expressed at each time point (6, 12, or 48 h of androgen exposure; Y-axis) in early prostate development (top half of heat map) are identical to genes differentially ...

Using laser capture microdissection to isolate and profile pure epithelial cell populations from the prostate, Tomlins et al. (2007) have assigned much of the variation in previously published prostate cancer gene expression studies to varying ratios of cancer epithelia to benign stroma. We therefore performed a second analysis, looking for enrichment of developmental genes in the cancer contrasts from Tomlins et al. (2007). Surprisingly, these analyses (Figure 3b) showed increased enrichment of embryonic genes compared with macrodissected bulk cancer samples (Figure 3a). Enrichment was significant for all time points (adj P<0.001), and included transitions for categories indicating malignant transformation (normal vs PIN), invasion (PIN vs cancer) and aggressiveness (Gleason grade; Figure 3b). The genes driving these processes are listed in Supplementary Tables 13-18. Thus, invasive prostate cancer adopts embryonic programs particularly during neoplastic transformation and invasion (PIN vs cancer) as indicated by enrichment across all three developmental time points (boxes 2, 5 and 6 in Figure 3b). In addition, genes differentially expressed in prostate cancer progression overlap most significantly with genes differentially expressed at the earliest (6 h) time point in prostate development, when prostate identity is first becoming established. Although most of these genes are repressed at the 6-h time point in development (Figure 2c), the same genes show a more equal balance between induction and repression (Supplementary Tables 13-18) in cancer. These observations identify androgen-responsive gene networks that establish prostate identity and operate in prostate cancer. Comparing and contrasting normal and malignant operation of these networks will likely yield useful insights into the molecular basis of prostate cancer.

Embryonic genes drive neoplastic transformation and invasion

Several shared functional themes from early prostate development also operate in human prostate cancer tissues, although at distinct stages of disease progression (Figure 4). Selected differentially expressed genes (P<0.05) and themes enriched in development and cancer are listed in Table 2. Complete lists are provided in Supplementary Tables 13-18. The transition between normal prostate and PIN is characterized by enrichment in pathways regulating transcription and apoptosis (Figure 4a). The latter observation is in keeping with striking suppression of programmed cell death previously observed in PIN (Zeng and Kyprianou, 2005). In contrast, the transition from noninvasive PIN to invasive cancer is enriched for extracellular matrix remodeling, cell motility and angiogenesis, as well as signaling pathways including Wnt and bone morphogenetic proteins (Figure 4a). These observations outline a program for prostate invasion that serves prostate development and cancer.

Figure 4
Transitions to increasingly invasive cancers are characterized by the activation of distinct pathways, transcription factors and microRNA target genes. Analysis of functional annotation in development and cancer transition reveals distinct (a) Gene ontology ...

Interestingly, other investigators have reported that metastatic potential can pre-exist in localized tumors rather than being induced by the external environment (Bernards and Weinberg, 2002; Kang et al., 2003). This concept may explain why the transition from localized to metastatic lesions did not represent themes expected to underlie metastasis, such as cell migration and adhesion. Instead, enriched processes included microtubule dynamics, mitotic spindle regulation and enhanced intracellular signaling (tyrosine kinase activity and adaptor proteins). Proteins in the latter class make excellent drug targets (Zwick et al., 2002), and to the extent that these genes might facilitate prostate cancer metastasis, they warrant further investigation.

Prostate budding recapitulates the aspects of neoplastic transformation

As a complement to statistical analyses and reverse transcriptase-PCR confirmation, we performed immunohistochemical analysis of early prostate development. Here, we show examples of two genes, Anxa1 and Sox9, which are coordinately suppressed and activated, respectively in prostate epithelial buds, mirroring their regulation in prostate carcinogenesis.

Annexin A (Anxa1) is a potent inducer of cell death in prostate cancer cell lines (Ornstein and Tyson, 2006). Reduced Anxa1 expression characterizes 91% percent of PIN lesions/cases, and is maintained in 94% of prostate cancers (Kang et al., 2002). Anxa1 expression has not been described in prostate development. In contrast to its loss in cancer, Anxa1 transcript expression was strongly induced at the 48-h time point (Supplementary Table 3). Immunohistochemical staining of developing prostates confirmed increased expression of Anxa1 in male tissue; however, there was a striking reduction of expression in emerging epithelial buds (Figures 5a-c). This pattern of expression is consistent with a role for Anxa1 repression in growth, survival, and/or invasion of prostate epithelium. Future studies involving genetic and/or pharmacologic manipulation of the developing prostate should further elucidate this role.

Figure 5
Annexin and Sox9 in epithelial invasion. (a and b) Immunohistochemical localization of Annexin A1 (AnxA1) in UGS tissue from female (a) and male (b) e17.5 littermates. (c) Immunohistochemical localization of Annexin A1 in male e18 UGS showing decreased. ...

Sox9 is a member of the SOX (Sry-related high mobility group) family of transcription factors that binds a DNA consensus sequence and enhances steroid hormone receptor binding. Functionally, it has been linked to a number of themes shared by development and cancer (Table 2). In developing prostate, Sox9 and its predicted transcriptional targets are differentially expressed during each of the three time points after androgen exposure (Supplementary Tables 1-3, 10-12). Immunohistochemical analysis revealed that Sox9 is expressed in early prostate epithelium and is particularly concentrated in the tips of invading prostate buds (Figure 5d). This expression pattern contrasts to that of Anxa1 with its reduced staining in bud tips and differs from p63, which is diffusely expressed in UGE at this stage of development. Altogether, this suggests Sox9 is critically located to facilitate prostate outgrowth, a hypothesis supported by recent work demonstrating its essential role in prostate development (Thomsen et al., 2008; Z Huang, B Simons, EM Schaeffer and DM Berman, unpublished observations). We noted that Sox9 transcripts and predicted targets were differentially expressed across prostate cancer progression (Supplementary Tables 14-18), and others have linked Sox9 to invasion, cell growth and metastasis in prostate cancer models.

We further characterized SOX9 expression in a panel of 219 human prostate cases, using immunohistochemistry and Immunofluorescence. In benign prostate tissue, Sox9 is expressed predominantly in basal epithelial cells, the stem cell compartment for prostate epithelia. (Xu et al., 2005; Wang et al., 2007a). We discovered aberrant Sox9 expression in neoplastic luminal epithelial cells in 51 of 53 (96%) PIN lesions (Supplementary Table 20), indicating an association with neoplastic transformation (Figures 5e and f). In prostate carcinomas (n = 105), basal cell gene expression is lost; however Sox9 expression is maintained in cancer cells (n = 59, 56%) with a trend toward increased expression with increasing grade (P = 0.1; Supplementary Table 20). Nuclear Sox9 was also noted in similar numbers (17/31, 55%) of advanced (lymph node metastasis) prostate cancer cases. In sum, alterations in the geographic and temporal location of Sox9 expression parallel the earliest events in the neoplastic transformation of prostate epithelial cells. These results suggest that Sox9 may play a critical role in the early, initiating phase of prostate carcinogenesis and contribute aspects of the basal/stem cell phenotype to prostate cancer. This pattern differs from recent work, which identified increased Sox9 in metastatic rather than localized lesions in a murine prostate cancer model (Acevedo et al., 2007). It is likely that Sox9 has several context-dependent functions in prostate epithelium and prostate cancer. The list of predicted Sox9 target genes that participate in prostate development (Supplementary Table 21) represents a useful starting point for further investigation.


Previous studies have linked primitive embryonic gene expression profiles to aggressive subsets of brain and lung cancers (Kho et al., 2004), thereby supporting balances between differentiation state and growth potential that operate similarly in organogenesis and tumorigenesis. Here, we identified prostate-specific programs for growth, survival, angiogenesis, and invasion that originate in organogenesis and are reactivated at specific steps in cancer progression. This curated list of androgen-regulated programs in prostate development and carcinogenesis provides a roadmap for understanding prostate growth and invasion. Prostate development in the mouse becomes a tractable experimental system in which to investigate the specific functions of these genes. This model bypasses the difficulty of probing gene regulation in human prostate cancer cell lines, which represent rare and possibly skewed exceptions to the rule that human prostate cancers do not adapt to growth in the laboratory. Unlike the extant models of prostate cancer involving transgenic mice, prostate development is unbiased by a preselected genetic lesion. Development is reproducible, genetically (Xin et al., 2003) and pharmacologically (Berman et al., 2004) tractable, and shown here to be reflected the entire spectrum of human prostate cancer progression.

A unique feature of these studies was the ability to induce prostatic lineage commitment and growth by controlled induction of signaling by AR, a gene with lineage-specific oncogenic properties (Garraway and Sellers, 2006). Evidence indicates the operation of other lineage-specific oncogenes in melanocytes (Garraway and Sellers, 2006), lung (NKX2-1; Weir et al., 2007) and elsewhere in the body (Garraway and Sellers, 2006), suggesting that gene expression programs relevant to other types of cancers can be identified, manipulated and modeled in the embryo.

Materials and methods

Mice and tissues

C57/Bl6J (The Jackson Laboratories, Bar Harbor, ME, USA) pregnancies were timed according to scheduled 4 h pairings. Paired pregnant females were injected intraperitoneally with dihydrotestosterone (Sigma-Aldrich, St Louis, MO, USA) at 50 mg/kg or triolein vehicle at e16.0. 6 or 12 h later, embryos were sexed and female UGS tissue dissected. For the 48-h time point tissue from unmanipulated male and female littermates was harvested at e17.5.

RNA isolation

Frozen tissue was homogenized and total RNA purified using the RNeasy system (Qiagen, Hilden, Germany) and analyzed using a Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA). Each time point yielded pools of ‘androgen exposed’ and ‘androgen naïve’ RNA. There were a total of 8,6 and 10 pools of RNA at the 6-, 12- and 48-h time points, respectively. Each pool was obtained from 3-4 individual UGS.

Probe synthesis, hybridization to oligonucleotide arrays and detection

Amplification, labeling, hybridization, and detection of 250 ng samples were carried out according to the manufacturer’s directions (Agilent).

Statistical analysis

Data were processed without background subtraction with packages from R/Bioconductor ( (Ihaka and Gentleman, 1996; Gentleman et al., 2004). Within-array ‘loess’ and between-arrays ‘scale’ normalization methods were applied to log2 expression ratios. Moderated t-statistics (by empirical Bayes shrinkage of standard errors (Smyth, 2004)) log odds ratios of differential expression and adjusted P-values (Benjamini et al., 2001) were obtained from a linear model accounting for biological replicates’ correlation, and dye/group effects.

Functional themes were obtained from Gene ontology (Ashburner et al., 2003), KEGG (Kanehisa et al., 2004), MsigDb(Subramanian et al., 2005). Enrichment analysis was performed by one-sided Wilcoxon test and multiple testing correction performed separately (Benjamini et al., 2001). Detailed statistical methods, public microarray data set identification and MIAME (Minimal Information about Microarray Experiments) information are provided in Supplementary materials. Data will be hosted in the Gene Expression Omnibus database.


Detection was performed as described in Berman et al.(2004) using Sox9 (Chemicon, Billerica, MA, USA) and Annexin A1 (Invitrogen, Carlsbad, CA, USA) antibodies.

Double immunofluorescence

Detection was performed with anti-Sox9 and anti-p63 (Santacruz, Santa Cruz, CA, USA) antibodies as described in supplement (Bivalacqua et al., 2007).

Tissue microarrays

Human tissues (Summarized in Supplementary Table 20) were stained with Sox9 antibody and scored as described in Supplementary Methods.

Supplementary Material

Supplementary Methods

Supplementary Figures 1-2

Supplementary Tables 1-3

Supplementary Tables 4-6

Supplementary Tables 7-9

Supplementary Tables 10-12

Supplementary Tables 13-18

Supplementary Tables 19-21


We thank W Matsui and N Watkins for comments on the manuscript, A DeMarzo and J Epstein for TMAs and YQ Chen for prostate regeneration data. These studies were funded by the Evensen Family, Passano and Patrick C Walsh Prostate Cancer Foundations, and NIH5K08DK059375 (DB) NIHK08 DK081019 (EMS), NIH5P30CA06973-39 (GP) and NSF034211 (GP and LM).


Supplementary Information accompanies the paper on the Oncogene website (


  • Acevedo VD, Gangula RD, Freeman KW, Li R, Zhang Y, Wang F, et al. Inducible FGFR-1 activation leads to irreversible prostate adenocarcinoma and an epithelial-to-mesenchymal transition. Cancer Cell. 2007;12:559–571. [PubMed]
  • Ashburner M, Mungall CJ, Lewis SE. Ontologies for biologists: a community model for the annotation of genomic data. Cold Spring Harb Symp Quant Biol. 2003;68:227–235. [PubMed]
  • Bailey P, Cushing H. Medulloblastoma cerbelli: a common type of midcerebellar glioma of childhood. Arch Neurol Psychiatry. 1925;14:192–223.
  • Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I. Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 2001;125:279–284. [PubMed]
  • Berman DM, Desai N, Wang X, Karhadkar SS, Reynon M, Abate-Shen C, et al. Roles for Hedgehog signaling in androgen production and prostate ductal morphogenesis. Dev Biol. 2004;267:387–398. [PubMed]
  • Bernards R, Weinberg RA. A progression puzzle. Nature. 2002;418:823. [PubMed]
  • Bhatia-Gaur R, Donjacour AA, Sciavolino PJ, Kim M, Desai N, Young P, et al. Roles for Nkx3.1 in prostate development and cancer. Genes Dev. 1999;13:966–977. [PubMed]
  • Bivalacqua TJ, Deng W, Kendirci M, Usta MF, Robinson C, Taylor BK, et al. Mesenchymal stem cells alone or ex vivo gene modified with endothelial nitric oxide synthase reverse age-associated erectile dysfunction. Am J Physiol Heart Circ Physiol. 2007;292:H1278–H1290. [PubMed]
  • Bostwick DG, Qian J. High-grade prostatic intraepithelial neoplasia. Mod Pathol. 2004;17:360–379. [PubMed]
  • Crawford HC, Fingleton BM, Rudolph-Owen LA, Goss KJ, Rubinfeld B, Polakis P, et al. The metalloproteinase matrilysin is a target of beta-catenin transactivation in intestinal tumors. Oncogene. 1999;18:2883–2891. [PubMed]
  • Cunha GR, Cooke PS, Kurita T. Role of stromal-epithelial interactions in hormonal responses. Arch Histol Cytol. 2004;67:417–434. [PubMed]
  • Donjacour AA, Thomson AA, Cunha GR. FGF-10 plays an essential role in the growth of the fetal prostate. Dev Biol. 2003;261:39–54. [PubMed]
  • Drews U, Sulak O, Oppitz M. Immunohistochemical localisation of androgen receptor during sex-specific morphogenesis in the fetal mouse. Histochem Cell Biol. 2001;116:427–439. [PubMed]
  • English HF, Santen RJ, Isaacs JT. Response of glandular versus basal rat ventral prostatic epithelial cells to androgen withdrawal and replacement. Prostate. 1987;11:229–242. [PubMed]
  • Friend SH, Bernards R, Rogelj S, Weinberg RA, Rapaport JM, Albert DM, et al. A human DNA segment with properties of the gene that predisposes to retinoblastoma and osteosarcoma. Nature. 1986;323:643–646. [PubMed]
  • Garraway LA, Sellers WR. Lineage dependency and lineage-survival oncogenes in human cancer. Nat Rev Cancer. 2006;6:593–602. [PubMed]
  • Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5:R80. [PMC free article] [PubMed]
  • Hanahan D, Weinberg RA. The hallmarks of cancer. Cell. 2000;100:57–70. [PubMed]
  • Ihaka R, Gentleman R. R: A language for data analysis and graphics. J Comput Graph Stat. 1996;5:299–314.
  • Ii M, Yamamoto H, Adachi Y, Maruyama Y, Shinomura Y. Role of matrix metalloproteinase-7 (matrilysin) in human cancer invasion, apoptosis, growth, and angiogenesis. Exp Biol Med (Maywood) 2006;231:20–27. [PubMed]
  • Irizarry RA, Warren D, Spencer F, Kim IF, Biswal S, Frank BC, et al. Multiple-laboratory comparison of microarray platforms. Nat Methods. 2005;2:345–350. [PubMed]
  • Joesting MS, Perrin S, Elenbaas B, Fawell SE, Rubin JS, Franco OE, et al. Identification of SFRP1 as a candidate mediator of stromal-to-epithelial signaling in prostate cancer. Cancer Res. 2005;65:10423–10430. [PubMed]
  • Jost A. Full or partial maturation of fetal endocrine systems under pituitary control. Perspect Biol Med. 1968;11:371. [PubMed]
  • Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M. The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004;32:D277–D280. [PMC free article] [PubMed]
  • Kang JS, Calvo BF, Maygarden SJ, Caskey LS, Mohler JL, Ornstein DK. Dysregulation of annexin I protein expression in highgrade prostatic intraepithelial neoplasia and prostate cancer. Clin Cancer Res. 2002;8:117–123. [PubMed]
  • Kang Y, Siegel PM, Shu W, Drobnjak M, Kakonen SM, Cordon-Cardo C, et al. A multigenic program mediating breast cancer metastasis to bone. Cancer Cell. 2003;3:537–549. [PubMed]
  • Kho AT, Zhao Q, Cai Z, Butte AJ, Kim JY, Pomeroy SL, et al. Conserved mechanisms across development and tumorigenesis revealed by a mouse development perspective of human cancers. Genes Dev. 2004;18:629–640. [PubMed]
  • Kousteni S, Bellido T, Plotkin LI, O’Brien CA, Bodenner DL, Han L, et al. Nongenotropic, sex-nonspecific signaling through the estrogen or androgen receptors: dissociation from transcriptional activity. Cell. 2001;104:719–730. [PubMed]
  • Lapointe J, Li C, Higgins JP, van de Rijn M, Bair E, Montgomery K, et al. Gene expression profiling identifies clinically relevant subtypes of prostate cancer. Proc Natl Acad Sci USA. 2004;101:811–816. [PubMed]
  • Lawson DA, Xin L, Lukacs RU, Cheng D, Witte ON. Isolation and functional characterization of murine prostate stem cells. Proc Natl Acad Sci USA. 2007;104:181–186. [PubMed]
  • Long X, Lin Y, Ortiz-Vega S, Yonezawa K, Avruch J. Rheb binds and regulates the mTOR kinase. Curr Biol. 2005;15:702–713. [PubMed]
  • Lu J, Qian J, Izvolsky KI, Cardoso WV. Global analysis of genes differentially expressed in branching and non-branching regions of the mouse embryonic lung. Dev Biol. 2004;273:418–435. [PubMed]
  • Manolagas SC, Kousteni S, Jilka RL. Sex steroids and bone. Recent Prog Horm Res. 2002;57:385–409. [PubMed]
  • McMenamin ME, Soung P, Perera S, Kaplan I, Loda M, Sellers WR. Loss of PTEN expression in paraffin-embedded primary prostate cancer correlates with high Gleason score and advanced stage. Cancer Res. 1999;59:4291–4296. [PubMed]
  • Memarzadeh S, Xin L, Mulholland DJ, Mansukhani A, Wu H, Teitell MA, et al. Enhanced paracrine FGF10 expression promotes formation of multifocal prostate adenocarcinoma and an increase in epithelial androgen receptor. Cancer Cell. 2007;12:572–585. [PMC free article] [PubMed]
  • Mohammadi M, Olsen SK, Ibrahimi OA. Structural basis for fibroblast growth factor receptor activation. Cytokine Growth Factor Rev. 2005;16:107–137. [PubMed]
  • Nam JS, Park E, Turcotte TJ, Palencia S, Zhan X, Lee J, et al. Mouse R-spondin2 is required for apical ectodermal ridge maintenance in the hindlimb. Dev Biol. 2007;311:124–135. [PMC free article] [PubMed]
  • Nusse R. Wnt signaling in disease and in development. Cell Res. 2005;15:28–32. [PubMed]
  • Olsen SK, Garbi M, Zampieri N, Eliseenkova AV, Ornitz DM, Goldfarb M, et al. Fibroblast growth factor (FGF) homologous factors share structural but not functional homology with FGFs. J Biol Chem. 2003;278:34226–34236. [PubMed]
  • Ornstein DK, Tyson DR. Proteomics for the identification of new prostate cancer biomarkers. Urol Oncol. 2006;24:231–236. [PubMed]
  • Ramaswamy S, Tamayo P, Rifkin R, Mukherjee S, Yeang CH, Angelo M, et al. Multiclass cancer diagnosis using tumor gene expression signatures. Proc Natl Acad Sci USA. 2001;98:15149–15154. [PubMed]
  • Scher HI, Sawyers CL. Biology of progressive, castration-resistant prostate cancer: directed therapies targeting the androgen-receptor signaling axis. J Clin Oncol. 2005;23:8253–8261. [PubMed]
  • Shen MM, Abate-Shen C. Pten inactivation and the emergence of androgen-independent prostate cancer. Cancer Res. 2007;67:6535–6538. [PubMed]
  • Smyth GK. Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3 article3 (e-pub 2004 Feb 12) [PubMed]
  • Stefani G, Slack FJ. Small non-coding RNAs in animal development. Nat Rev Mol Cell Biol. 2008;9:219. [PubMed]
  • Stehelin D, Varmus HE, Bishop JM, Vogt PK. DNA related to the transforming gene(s) of avian sarcoma viruses is present in normal avian DNA. Nature. 1976;260:170–173. [PubMed]
  • Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005;102:15545–15550. [PubMed]
  • Thomsen MK, Butler CM, Shen MM, Swain A. Sox9 is required for prostate development. Dev Biol. 2008;316:302–311. [PubMed]
  • Thomson AA. Role of androgens and fibroblast growth factors in prostatic development. Reproduction. 2001;121:187–195. [PubMed]
  • Tomlins SA, Mehra R, Rhodes DR, Cao X, Wang L, Dhanasekaran SM, et al. Integrative molecular concept modeling of prostate cancer progression. Nat Genet. 2007;39:41–51. [PubMed]
  • Treister NS, Richards SM, Lombardi MJ, Rowley P, Jensen RV, Sullivan DA. Sex-related differences in gene expression in salivary glands of BALB/c mice. J Dent Res. 2005;84:160–165. [PubMed]
  • Verras M, Sun Z. Roles and regulation of Wnt signaling and beta-catenin in prostate cancer. Cancer Lett. 2006;237:22–32. [PubMed]
  • Wang H, McKnight NC, Zhang T, Lu ML, Balk SP, Yuan X. SOX9 is expressed in normal prostate basal cells and regulates androgen receptor expression in prostate cancer cells. Cancer Res. 2007a;67:528–536. [PubMed]
  • Wang S, Gao J, Lei Q, Rozengurt N, Pritchard C, Jiao J, et al. Prostate-specific deletion of the murine Pten tumor suppressor gene leads to metastatic prostate cancer. Cancer Cell. 2003;4:209–221. [PubMed]
  • Wang XD, Wang BE, Soriano R, Zha J, Zhang Z, Modrusan Z, et al. Expression profiling of the mouse prostate after castration and hormone replacement: implication of H-cadherin in prostate tumorigenesis. Differentiation. 2007b;75:219–234. [PubMed]
  • Weir BA, Woo MS, Getz G, Perner S, Ding L, Beroukhim R, et al. Characterizing the cancer genome in lung adenocarcinoma. Nature. 2007;450:893–898. [PMC free article] [PubMed]
  • Wilson JD, Griffin JE, George FW. Sexual differentiation: early hormone synthesis and action. Biol Reprod. 1980;22:9. [PubMed]
  • Xin L, Ide H, Kim Y, Dubey P, Witte ON. In vivo regeneration of murine prostate from dissociated cell populations of postnatal epithelia and urogenital sinus mesenchyme. Proc Natl Acad Sci USA. 2003;100(Suppl 1):11896–11903. [PubMed]
  • Xin L, Lawson DA, Witte ON. The Sca-1 cell surface marker enriches for a prostate-regenerating cell subpopulation that can initiate prostate tumorigenesis. Proc Natl Acad Sci USA. 2005;102:6942–6947. [PubMed]
  • Xu J, Lowey J, Wiklund F, Sun J, Lindmark F, Hsu FC, et al. The interaction of four genes in the inflammation pathway significantly predicts prostate cancer risk. Cancer Epidemiol Biomarkers Prev. 2005;14:2563–2568. [PubMed]
  • Yeager N, Brewer C, Cai KQ, Xu XX, Di Cristofano A. Mammalian target of rapamycin is the key effector of phosphatidylinositol-3-OH-initiated proliferative signals in the thyroid follicular epithelium. Cancer Res. 2008;68:444. [PubMed]
  • Zeng L, Kyprianou N. Apoptotic regulators in prostatic intraepithelial neoplasia (PIN): value in prostate cancer detection and prevention. Prostate Cancer Prostatic Dis. 2005;8:7–13. [PubMed]
  • Zhang TJ, Hoffman BG, Ruiz de Algara T, Helgason CD. SAGE reveals expression of Wnt signalling pathway members during mouse prostate development. Gene Expr Patterns. 2006;6:310–324. [PubMed]
  • Zwick E, Bange J, Ullrich A. Receptor tyrosine kinases as targets for anticancer drugs. Trends Mol Med. 2002;8:17–23. [PubMed]