- Research article
- Open Access
Evolutionary conservation of zinc finger transcription factor binding sites in promoters of genes co-expressed with WT1 in prostate cancer
BMC Genomics volume 9, Article number: 337 (2008)
Gene expression analyses have led to a better understanding of growth control of prostate cancer cells. We and others have identified the presence of several zinc finger transcription factors in the neoplastic prostate, suggesting a potential role for these genes in the regulation of the prostate cancer transcriptome. One of the transcription factors (TFs) identified in the prostate cancer epithelial cells was the Wilms tumor gene (WT1). To rapidly identify coordinately expressed prostate cancer growth control genes that may be regulated by WT1, we used an in silico approach.
Evolutionary conserved transcription factor binding sites (TFBS) recognized by WT1, EGR1, SP1, SP2, AP2 and GATA1 were identified in the promoters of 24 differentially expressed prostate cancer genes from eight mammalian species. To test the relationship between sequence conservation and function, chromatin of LNCaP prostate cancer and kidney 293 cells were tested for TF binding using chromatin immunoprecipitation (ChIP). Multiple putative TFBS in gene promoters of placental mammals were found to be shared with those in human gene promoters and some were conserved between genomes that diverged about 170 million years ago (i.e., primates and marsupials), therefore implicating these sites as candidate binding sites. Among those genes coordinately expressed with WT1 was the kallikrein-related peptidase 3 (KLK3) gene commonly known as the prostate specific antigen (PSA) gene. This analysis located several potential WT1 TFBS in the PSA gene promoter and led to the rapid identification of a novel putative binding site confirmed in vivo by ChIP. Conversely for two prostate growth control genes, androgen receptor (AR) and vascular endothelial growth factor (VEGF), known to be transcriptionally regulated by WT1, regulatory sequence conservation was observed and TF binding in vivo was confirmed by ChIP.
Overall, this targeted approach rapidly identified important candidate WT1-binding elements in genes coordinately expressed with WT1 in prostate cancer cells, thus enabling a more focused functional analysis of the most likely target genes in prostate cancer progression. Identifying these genes will help to better understand how gene regulation is altered in these tumor cells.
In the United States, prostate cancer is the most common form of cancer in men and is the second most deadly cancer in men killing more than 27,000 annually . Nearly one in six men will develop prostate cancer at some point in their life, with the majority of incidences occurring after the age of 50. The major biomarker for prostate cancer diagnosis is prostate specific antigen (PSA), however, the sensitivity and specificity of the PSA assay is limited . Improved biomarkers will result from a better understanding of molecular mechanisms that regulate this disease.
Global gene expression analyses have led to a better understanding of growth control of prostate cancer cells [3–5]. Ongoing studies identified more than 200 genes predominantly expressed in prostate cancer epithelial cells  and included genes likely to influence growth of prostate cancer cells, such as growth factors, growth factor receptors and TFs (as identified by Gene Ontology and KEGG pathway analyses). Two of the TFs identified in the prostate cancer epithelial cells were the Wilms tumor gene (WT1) and the early growth response gene (EGR1), zinc finger transcription factors that bind at G-rich promoters of genes that regulate growth. In fact, the WT1 TF binds at several G-rich sites (GNGNGGGNG), including the EGR1 consensus binding site GCGGGGGCG [7–9]. Both WT1 and EGR1 have been identified in prostate cancer cells, although their function in prostate epithelium is unknown [10–12]. WT1 has an essential role in the normal development of the urogenital system and has been shown to suppress transcription of the promoters of many important growth factors .
While identifying prostate growth control pathways potentially regulated by WT1, we have focused our studies on candidate genes belonging to known growth regulatory pathways. We have previously described WT1 regulation of the androgen receptor (AR) and vascular endothelial growth factor (VEGF) gene promoters [14, 15]. To go beyond the candidate genes approach and identify novel gene targets coordinately expressed with WT1 in tumor epithelial cells, a more systematic and unbiased high-throughput computational approach was used. These in silico analyses were based on 24 genes expressed in prostate cancer epithelium that were likely to influence growth of prostate cancer cells. Putative TFBS were computationally predicted; however, the identification of functional TFBS is a challenge and requires an alternative approach. Availability of complete genomic sequence from multiple species allows identification of evolutionary conserved elements, e.g. cis-regulatory elements. Functionally important elements are likely to experience purifying selection pressure [16–20], thus, we can utilize the degree of evolutionary conservation to identify TFBS that are likely to be functional. Our approach was to identify regions (and TFBS) evolutionary conserved across multiple mammalian genomes, including those separated by 170 million years (human and opossum) . Overall, this targeted approach identified important candidate binding elements in genes coordinately expressed with WT1 in prostate cancer epithelial cells. Identifying genes regulated by zinc finger TFs expressed in prostate cancer cells will enhance understanding of the altered pathways in these tumor cells and provide useful biomarkers for prostate cancer progression.
Evolutionary conservation analysis: TFBS conserved in prostate cancer growth genes
Genomic sequences of proximal promoter regions of 24 genes expressed in prostate cancer epithelial cells (Additional file 1) were analyzed to determine the degree of evolutionary conservation and to identify potentially important regulatory regions. Binding sites for six TFs (WT1, EGR1, SP1, SP2, AP2, and GATA1) were investigated for evolutionary conservation over a range of eight different mammalian species (human, chimpanzee, macaque, cow, dog, mouse, rat and opossum) (Table 1). Tables 2 and 3 highlight 11 of these genes whose promoter sequences could be aligned in at least five mammalian species (human, chimpanzee, macaque, rat and mouse) and were found to have at least one evolutionary conserved TFBS.
Among the TFBS investigated, WT1, EGR1 and SP1 sites showed the highest frequency of evolutionary conservation in the gene promoters surveyed. For example, the promoters of EGR1, GATA2 and WT1 were found to have multiple WT1, EGR1 and SP1 candidate binding sites that were conserved through multiple species (Table 3). In the EGR1 promoter, 50% of WT1 sites are conserved between human and primates. Additionally, in the GATA2 gene promoter, 94% of WT1 sites, 70% of SP1 sites, and 100% of EGR1 sites are conserved between human and other primates (Table 3). Similarly, in the WT1 gene promoter 50% of SP1, 43% of WT1 and 100% of EGR1 sites are conserved between human and other primates (Table 3). WT1, EGR1, and SP1 TFBS within the promoters of IGFBP2, KLK3, NPY, SOX4, SOX9, and TFAP2C are also conserved between human and other primates (Table 3).
Importantly, for the WT1 and EGR1 gene promoters this conservation extended into the marsupials (Table 4). The EGR1 gene promoter is relatively conserved between human and opossum with 20% of predicted EGR1, 12% of predicted WT1 and 14% of predicted SP1 sites conserved between human and opossum. Similarly, the WT1 gene promoter exhibited conservation between human and opossum, with 33% of predicted SP1 and 14% of predicted WT1 sites shared between human and opossum. In the GATA2 promoter only 12% of predicted WT1 sites are shared between human and opossum (Table 4). Overall, TFBS for the three TF (WT1, EGR1, and SP1) were evolutionary conserved between human and the distantly related species, opossum, in seven different promoters (WT1, EGR1, GATA2, IGFBP2, SOX4, SOX9, and TFAP2C).
Tables 3 and 4 show that there were fewer SP2, AP2 and GATA1 than WT1, EGR1 and SP1 TFBS in the 11 gene promoters analyzed. While evolutionary conservation between primates was similar for all six TFBS, conservation between human and rodents diminished for SP2 and AP2 TFBS. AP2 sites in the promoters of the GATA2, WT1, and NPY genes showed 25% to 100% conservation between human and other primates. Conservation of AP2 sites was the strongest in the NPY gene promoter as these sites are also conserved between human and opossum (Table 4). In addition to conservation of GC-rich TFBS, the AT-rich GATA1 binding sites were shown to be highly conserved in several gene promoters including SOX4, EGR1, IGFBP2 and NPY (Table 3). All of the GATA1 sites in these four promoters are conserved between human and chimpanzee, and for the SOX4 gene promoter this strong conservation extends to rodents as well.
The overall evolutionary conservation of predicted TFBS of these 11 different genes expressed in prostate cancer cells was analyzed. As would be expected, conservation of TFBS decreased as species became more evolutionarily divergent (Table 5). TFBS were found to be the most conserved among primates, followed by rodents, and the least amount of conservation was found between human and opossum. Of the 47 predicted WT1 sites in the 11 genes analyzed, primates had 68% of these sites conserved between human and primate genomes, while rodent genomes had only 15% of these sites being conserved, and opossum only 6% of these sites conserved, clearly showing a drastic drop in conservation as species diverge. This same pattern was shown for the other TFBS that were analyzed, including EGR1 and SP1. In particular, 85% of the EGR1 sites were conserved between human and other primates, 26% between human and rodents, and 19% between human and opossum. Similarly, there were 50 predicted SP1 binding sites, of which 62%, 22% and 12% were conserved between human and primates, rodents, and opossum genomes, respectively, therefore, exhibiting decreasing conservation of these sites with evolutionary divergence. Thus, with this approach of identifying evolutionary conserved sequences we were able to pinpoint specific candidate binding sites that could be tested for functional relevance.
Conservation of overlapping WT1, EGR1, and SP1 TFBS
Several of the genes investigated have multiple overlapping WT1, EGR1, and SP1 binding sites in their proximal promoter regions. For example, the promoter of the human EGR1 gene has evolutionary conserved overlapping WT1/SP1 binding sites (one of which is shown in Figure 1A). Both the overlapping WT1 (human 565–581) and SP1 (human 563–577) sites are conserved between seven of eight species compared, and the SP1 site is also conserved between human and opossum. A second WT1 site (human 614–630) located 33 bp downstream overlaps an EGR1 site (human 608–624) and both sites are conserved among all eight species, including opossum (Figure 1A). The promoter of the GATA2 gene also contained overlapping SP1 and WT1 TFBS (located in human positions 1125–1139 and 1127–1145, respectively) that are conserved among several mammalian genomes (Figure 1B). The WT1 gene promoter also has overlapping WT1/SP1 binding sites and when aligned with multiple species, one 3' WT1 site (human 1444–1468) was conserved between all primates, rodents, and opossum, thus, depicting millions of years of conservation of this particular site (Figure 1C). The SP1 site (human 1420–1434) is conserved between all primates and rodents tested, and overlaps with a WT1 site (human 1409–1425) that is conserved between human and chimpanzee (Figure 1C). Interestingly, the sequence similarity is so great between human and chimpanzee for this WT1 promoter region that no insertions or deletions were observed in either genomic sequence; thus, these TFBS were located in exactly the same positions relative to the start ATG codon.
Identification of overlapping TFBS in the gene promoters indicated that WT1 and EGR1 may compete for binding. Analyses of the promoter regions of 11 genes expressed in prostate cancer epithelial cells showed that WT1 TFBS overlapped SP1 and EGR1 TFBS, either separately or together. Overall, it was found that there were 25 overlapping sites in the promoter regions of these genes. There were 12 WT1/SP1, seven SP1/EGR1, three WT1/EGR1, and three WT1/SP1/EGR1 overlapping sites (Table 6). These overlapping sites were found in 10 of the 11 gene promoters analyzed. Seven overlapping sites were identified in the promoter region of the EGR1 gene, and three of these seven overlapping sites are conserved between human and other species. Three other gene promoters, GATA2, IGFBP2, and TFAP2C, have three overlapping sites each, with one SP1/EGR1 site conserved between human and opossum for both the TFAP2C and IGFBP2 promoters. The WT1 and KLK3 promoters have overlapping WT1/SP1 and SP1/EGR1 sites, respectively. All of these overlapping TFBS are excellent candidates for functional testing to determine whether competition for TF binding at these sites results in activation or suppression of the genes they are regulating.
Sequence conservation of TFBS indicates a potentially functional WT1 binding site in the KLK3 (PSA) promoter
One of the 24 genes differentially expressed in prostate cancer epithelial cells was KLK3 (PSA), an important diagnostic marker. Sequence alignment of the KLK3 promoter revealed three WT1 sites and two SP1 sites, with two-thirds of the WT1 and one-half of the SP1 sites conserved between human and other primates (Table 3). Given the premise that evolutionary conserved sites are more likely to be functionally relevant, we tested these conserved sites for their ability to bind TF in vivo. PCR primers were designed to flank the region where adjacent conserved WT1 (human 1332–1352) and the SP1 sites (human 1404–1418) were identified (Figure 2A). Both of these binding sites in the PSA promoter were tested by chromatin immunoprecipitation (ChIP) in hormone responsive LNCaP prostate cancer cells (Figure 2B). Since LNCaP cells express little WT1 , they were transfected with a green fluorescent protein (GFP)-tagged WT1 expression construct 48 hours prior to the ChIP assay. After crosslinking, the chromatin and TF complexes were immunoprecipitated by both WT1 and SP1antibodies, as demonstrated by PCR amplification of the promoter region. WT1 and SP1 may bind at adjacent sites within the PSA promoter or at overlapping sites, since the SP1 site overlaps the EGR1 site, to which WT1 may also bind [7–9]. The importance of these WT1 and SP1 TFBS as candidate binding sites was confirmed by the in vivo ChIP assay.
Functional WT1 and SP1 binding sites in the VEGF promoter are conserved between human and other primates
Having tested the significance of identified evolutionary conserved sites, we then asked whether TFBS known to mediate transcriptional regulation would also be conserved. Two genes that regulate prostate cancer progression by enhancing growth and blood supply, AR and VEGF, have multiple WT1 and SP1 binding sites in their proximal promoter regions [14, 15, 23–25]. We have previously identified an EGR1 site in the VEGF promoter that binds both WT1 and SP1 protein in vitro , and here demonstrate by ChIP assay that this promoter region binds WT1 and SP1 in vivo (Figure 3). Chromatin from both embryonic kidney 293 cells and LNCaP cells expressing a GFP-tagged WT1 expression construct was immunoprecipitated by WT1 and SP1 antibodies and amplified by PCR. Using primers specific for the VEGF proximal promoter region, products ~140 bp in size were amplified from chromatin of both 293 and LNCaP cells (Figure 3A and 3B). These ChIP assays also demonstrated selective WT1 binding, since an adjacent site 190 nucleotides downstream failed to bind WT1 in the same assay (data not shown). These sites were validated as being transcriptionally regulated in several different assays, including luciferase reporter assays , so we asked whether they were evolutionary conserved in different species. In silico analyses predicted that an overlapping EGR1 (human 1717–1733) and SP1 (human 1721–1735) site and a WT1 site (human 1755–1771) were conserved between primates and dogs, but not in rodents (Figure 3C). Furthermore, as seen with the PSA promoter region, WT1 and SP1 may bind at adjacent sites or potentially at overlapping sites since WT1 also binds at EGR1 sites [7–9]. Both PSA and VEGF promoter regions contain evolutionary conserved WT1 sites adjacent to overlapping EGR1/SP1 TFBS, to which WT1 is also likely to bind, thus facilitating either cooperation or competition between TFs.
Similarly, WT1 binding sites previously identified in the AR proximal promoter region by EMSA analysis and verified to mediate transcriptional regulation in luciferase reporter assays [14, 23] were confirmed by ChIP using PCR primers flanking the WT1 and SP1 TFBS (Figure 4A). Since these binding sites were tested in vivo, evidence of sequence conservation was sought, as described. As shown in Figure 4B, both a WT1 site (human 1434–1450) and an EGR1 site (human 1524–1537) were identified within the region amplified by ChIP. This less common pyrimidine-rich EGR1 TFBS, consisting of TCC repeats, has been shown to bind both WT1 and SP1 [7, 14, 26], thus all three zinc finger TFs could compete for binding at this site. Evidence for evolutionary conservation between human and other primates was limited by the lack of genomic sequence information available for chimpanzee (and lack of conservation between human and macaque).
Identification of evolutionary conserved sequences derived from comparisons of multiple genomes (so-called "phylogenetic footprints") has been successful in identifying functionally important regions, including those regions that regulate gene expression [19, 27–34]. However, some regulatory genomic sequences do not appear to be conserved or the level of evolutionary conservation varies between different genomic comparisons [35, 36]. Importantly, some functional regions have been reported to experience a relatively fast rate of turnover, where the functional significance of the element is retained despite changes at the nucleotide sequence level (e.g., transcription start sites, ). Thus, it is likely that gene expression in mammalian genomes is controlled by both types of regulatory elements, i.e., those elements that exhibit evolutionary and functional conservation and those that exhibit functional conservation only. Moreover, while numerous algorithms are available to computationally predict potential regulatory elements, it is often challenging to narrow down the list of those that are likely to be functional, particularly for relatively short elements such as TFBS. One of the approaches that utilizes evolutionary conservation as a predictor of TFBS functionality is the rVISTA tool that uses pairwise sequence alignments to identify the most highly conserved TFBS between the pair of genomic sequences . Another set of tools, the Mulan, takes advantage of evolutionary conservation information obtained from multi-sequence alignments of several genomes . However, the latter requires the TFBS to be shared among all genomes present in the alignment  and may potentially miss the lineage-specific regulatory elements that are absent from some subsets of genomes. Therefore, in this work we used TFBS elements shared between some but not necessarily all of the available genomes.
We used evolutionary sequence conservation, as determined by both the multi-species sequence alignments and the in silico TFBS predictions, to identify those sites most likely to regulate expression of target genes that influence growth of prostate cancer cells. Regulatory regions with functional importance can be expected to exhibit sequence conservation due to selection. Thus, predicted TFBS that are located in the orthologous positions in multiple genomes are likely to be truly functional. Our identification of evolutionary conserved WT1 and SP1 binding sites in the PSA promoter indeed supports this notion (Figure 2). As expected, conservation of TFBS decreased as species became more evolutionarily divergent , so those TFBS that were conserved between multiple species including opossum are more likely to be functionally important in the regulation of gene expression.
The abundance of overlapping zinc finger TFBS also supported the functional importance of these regulatory regions. Thus, we identified many TFBS in potential target genes that were co-expressed with WT1 in prostate cancer epithelial cells. Evolutionary conserved WT1 and SP1 sites in the PSA promoter were confirmed by ChIP to bind both WT1 and SP1 in LNCaP prostate cancer cells chromatin. Although it is a novel discovery that both SP1 and WT1 bind the PSA promoter and may play a role in its regulation, reporter assays are needed to confirm their contribution to transcription. In addition, a WT1 binding site known to transcriptionally regulate the VEGF promoter  was confirmed by ChIP and found to be in an evolutionary conserved region. Interestingly, transcriptionally active WT1 and EGR1 binding sites in the AR promoter  were not conserved between human and macaque, although adjacent genomic regions could be aligned between multiple species (Figure 4). This suggests that the AR promoter may have experienced faster turn-over than the VEGF promoter, yet remained functionally conserved despite sequence changes at the nucleotide level.
Many of the genes expressed in prostate cancer epithelial cells have previously been reported to interact and regulate each other, suggesting multiple potential targets for altered pathways that may lead to prostate cancer progression. We and others have identified gene interactions [8, 14, 15, 23, 41–47] that are consistent with WT1 regulating the progression and/or growth of tumors in the prostate. However, PSA was a candidate gene target identified by our in silico evolutionary conservation approach and confirmed by in vivo chromatin binding assays. PSA is a member of the kallikrein family of serine proteases and is a marker of epithelial differentiation in the prostate . It is up-regulated in prostate cancer cells when compared to normal adjacent tissue  and its expression is regulated by the ligand bound androgen receptor (AR) . Since WT1 activates the AR promoter in prostate cancer cells , this suggests that WT1 may directly or indirectly regulate PSA gene expression.
In addition to PSA, genes that were co-expressed with WT1 in prostate cancer epithelial cells and that could potentially interact with, or be regulated by, WT1 included GATA2, ECAD, EGR1, and NDRG1 . GATA binding proteins are zinc finger transcription factors that bind the WGATAR consensus motif and are expressed in multiple tissues, including endocrine glands [50–52]. Interestingly, GATA TFs regulate WT1 expression, as multiple GATA TFBS are found within the WT1 promoter and enhancer regions [53–55]. GATA binding protein 2 (GATA2) has been shown to be one of the main GATA family members expressed in the prostate of human and mouse . It has been suggested that GATA2 plays a role in androgen mediated regulation of PSA expression, possibly through interaction with AR, as GATA sites are adjacent to AR TFBS in the PSA promoter . WT1 could contribute to GATA2 mediated regulation of target genes in prostate cancer cells, if WT1 also physically interacts with GATA2. This notion is consistent with the observation that WT1 interacts with GATA4 to regulate SRY gene expression . This complex pattern of zinc finger-protein interaction between WT1 and GATA, along with regulation of WT1 expression by GATA TF, suggests a potential for WT1 feedback control of GATA activity.
The WT1 promoter is itself a target of autoregulation by WT1 . WT1 is a multifunctional transcription factor; its four major isoforms are formed by alternative splicing at two sites resulting in the inclusion or exclusion of (1) exon V and/or (2) a tripeptide (KTS) in exon 9 that alters the zinc finger DNA binding structure . While the functions of the various isoforms of WT1 are still being discovered, the -KTS isoform is a transcriptional regulator with G-rich recognition sequence . The +KTS isoform is also likely to be present in prostate cancer cells but would contribute to gene regulation via splicing and post-transcriptional gene regulation [59, 60]. Here we have identified potential target genes with well-described DNA binding sites recognized by the -KTS isoform and have not assessed the less well understood RNA binding sites recognized by the +KTS isoform .
The early growth response 1 gene (EGR1) is a homolog of WT1 . Although it has only three zinc-fingers, it shares some TFBS with WT1. EGR1 has been implicated as a cancer suppressor gene and activates genes required for differentiation . In human prostate cancer, EGR1 is over-expressed [11, 12] and in a mouse model of prostate cancer, EGR1 regulates genes essential for progression of tumor growth . Since WT1 regulates the EGR1 promoter in vitro  it may indirectly regulate other EGR1 target genes, such as the N-myc downstream regulated gene 1 (NDRG1), an α/β hydrolase. In many cancer cell lines it has been shown to be up-regulated by both hypoxia and hormone treatment suggesting that it could be linked to androgen induced differentiation and signaling in the prostate [62, 63]. Since EGR1 regulates NDRG1, WT1 could either directly or indirectly regulate NDRG1.
While analyzing the homologous sequences of the different gene promoters, numerous overlapping TFBS were found, suggesting competition for binding and differential regulation of these gene promoters. Several studies have shown that EGR1 and SP1 TFBS often overlap [7, 64, 65]. When EGR1 binds to a site also bound by SP1, it displaces the SP1 "activator" from the binding site and represses transcription of these genes . For example, the promoter of NDRG1 was shown to be regulated by an overlapping EGR1/SP1 binding site  (located outside of the surveyed region of our study). It was shown that this evolutionary conserved site was vital in positively regulating expression of NDRG1 . Similarly, our results showed evolutionary conserved overlapping EGR1/SP1 sites in several other gene promoters, including VEGF and PSA. In the latter, overlapping EGR1/SP1 sites were found to be conserved between human and two other primate species (chimpanzee and macaque).
Additionally, WT1 and EGR1 compete for binding at shared TFBS. WT1 recognizes and binds to EGR1 sites on the promoters of many different genes [7, 9, 66–68]. WT1 generally functions as a transcriptional repressor when bound to EGR1 TFBS in the transforming growth factor-beta 1 (TGFβ1) and EGR1 promoters, while EGR1 functions as an activator [8, 9]. Many gene promoters with overlapping WT1, EGR1, and SP1 binding sites have been identified (reviewed in ). For example, three-way competition occurs between EGR1, SP1 and WT1 for binding and regulation of superoxide dismutase expression . However, the mechanisms of gene regulation at overlapping sites, including TF competition, are not well understood.
Combinations of adjacent and overlapping EGR1, WT1 and SP1 TFBS conserved between multiple species were found in multiple gene promoters. Adjacent sites were found in the PSA promoter where an overlapping EGR1/SP1 site is 50 base pairs downstream of a WT1 site and in the VEGF promoter where an EGR1/SP1 overlapping site is 20 base pairs away from a WT1 site. Such sites can facilitate synergistic interactions or may be required for inducible expression, as described for AR and GATA2 interactions in the PSA promoter . Additionally, in the VEGF promoter an SP1 site adjacent to a non-canonical estrogen receptor (ER) TFBS contributes to hormone induction of VEGF expression . Similarly, WT1 appears to interact with ER at neighboring sites in the insulin like growth factor 1 receptor (IGF1R) promoter . These complex arrangements of EGR1, WT1 and SP1 TFBS could facilitate cooperative or competitive binding by these TFs and would have pleiotropic effects on the regulation of these genes. Genes with evolutionary conserved overlapping TFBS could be part of a prostate epithelial cell transcriptome regulated by WT1.
Genes coordinately expressed in prostate cancer epithelial cells have conserved regulatory elements and an abundance of overlapping zinc finger TFBS. Potential WT1 gene targets were identified based on TFBS sequence conservation, and the significance of the WT1 TFBS in the PSA promoter was verified in vivo by ChIP assays. Similarly, a transcriptionally active WT1 binding site in the VEGF promoter was confirmed by ChIP and found to be in a region conserved amongst primates. Thus, these genes could be part of a novel network of regulatory pathways initiated by WT1 and have important implications in the progression of prostate cancer.
Promoter sequence compilation
For each of the 24 prostate cancer growth regulatory genes differentially expressed, the complete or draft genomes of eight different mammalian species were downloaded from the Ensembl Genome Browser [71, 72]. The following genome assemblies were used: the NCBI 36 assembly of human (Homo sapiens) genome, the NCBI m36 assembly of mouse (Mus musculus) genome, the Pan Tro 2.1 assembly of chimp (Pan troglodytes) genome, a whole genome shotgun (WGS) preliminary assembly Btau_3.1 of cow (Bos Taurus) genome, a WGS assembly Can Fam2.0 of dog (Canis familiaris) genome, a WGS preliminary assembly Mmul_1 of rhesus monkey (Macaca mulatta) genome, the Mon Dom5 assembly of opossum (Monodelphis domestica) genome, and the RGSC3.4 assembly of rat (Rattus noregicus) genome. Since major regulatory elements are located within several hundred base pairs of transcription start sites , 1.5 kb of human nucleotide sequence 5' of the translational start site (that is, 5' of the first exon as defined in Ensembl ) was collected. Orthologous sequences from other mammalian genomes were obtained from respective genome assemblies. In the case of the EGR1 promoter this extended beyond 1.5 kb, so was assigned a negative number. The genome viewer and annotation program Artemis was used to ensure the correct context of genomic sequences . In each sequence, the nucleotide positions were numbered sequentially, with the targeted promoter region occupying positions 1 through 1500 (5' to 3' direction) of the forward strand, and ATG start codon located at positions 1501–1503 of the genomic sequence analyzed.
AR and VEGF promoter sequences containing the functional WT1 TFBS for the human AR and VEGF promoters were obtained from Ensembl (ENSG00000169083 and ENSG00000112715, respectively). For alignment analyses of known functional sites [14, 15], an orthologous promoter region (3 kb) was then collected from eight mammalian genomes as described above.
TFBS predictions, evolutionary conservation and multiple sequence alignments
TFBS of WT1, EGR1, SP1, SP2, AP2 and GATA1 were predicted for each gene by the program MatInspector  that utilizes the TRANSFAC libraries of TF binding motifs [75, 76]. The default parameters of similarity thresholds were used for all examined genes, and they were as follows: core similarity > 0.75 and optimized matrix similarity thresholds (i.e., those that minimize false positives for each individual matrix as available in the library) . In MatInspector, core similarity is one of the built-in program parameters that determines whether the observed sequence match will be analyzed further. It refers to the four most conserved consecutive nucleotides of the matrix, usually the most critical sites for protein binding, and reaches 1.0 only when there is a perfect match [75, 77]. Sequence matches with low core similarity (less than 0.75) are not, by default, reported to the user. Vertebrate matrices of the Matrix Family Library Version 6.2 (October 2006) that included 464 matrices were used . Multiple sequence alignments of the promoter sequences were reconstructed with the program blastZ using MultiPipMaker , and predicted human TFBS were mapped onto the alignments.
Regions that are conserved in multiple genomes are often found to correspond to functionally important ones . However, because of the species-specific differences in gene regulation due to underlying differences in morphogenesis and development, such as those between different segments of human and rodent prostate , it can be expected that some functionally important regions will be conserved only in a limited set of genomes where they play a critical role. Thus, we used a flexible definition of "evolutionary conservation" to accommodate such potential differences between genes and/or TFBS: here a TFBS was considered evolutionary conserved if it was predicted as a respective TFBS in orthologous position in at least three of eight surveyed genomes. In other words, the same genomic region was predicted to function as a candidate binding site for a particular TF in at least 3 surveyed genomes. Further, because differences in presence/absence of particular TFB sites between genomes may also be attributed to differences in the role of respective genes in each of the organisms, we examined evolutionary conserved sites at different levels of resolution: Human-Primates, Human-Rodents, and Human-Opossum, thereby, allowing us to identify genes and TFB sites that are functionally relevant to each of these comparisons.
Cell culture and chromatin immunoprecipitation
LNCaP prostate cancer cells (ATCC-CRL 1740) and human embryonic kidney 293 cells (ATCC-CRL 1573) were cultured in RPMI or DEM/F12 (HyClone Laboratories, Utah) media, respectively, as described . The cytomegalovirus (CMV) promoter-driven pGFP-WT1 (A) expression construct encoding the murine Wt1 gene (lacking both KTS insertion and exon 5) fused to GFP coding region were obtained from Dr. A. Ward . All DNA was purified by the Qiagen plasmid Maxi Kit (Qiagen, Carlsbad CA) and transfections were performed using lipofectamine 2000 (Invitrogen; Carlsbad CA) in serum- and antibiotic-free media as described . Green fluorescing cells were visualized by epifluorescence microscopy (Olympic) at 100–400× magnification at 24 and 48 hrs after transfection prior to cell harvest for chromatin isolation.
The Farnham ChIP protocol  was used with some modifications. Two million cells were treated with formaldehyde to crosslink proteins to DNA and lysed in PBS-PI as recommended for the EZ ChIP Assay (Upstate Biotechnology Inc). Lysates were centrifuged and DNA sheared by sonication (Biosonik III, Bronwill Scientific, Rochester, NY) to fragments of 100–400 bp in length. The supernatant was pre-cleared by incubation with Protein G Agarose and incubated overnight at 4°C with either SP1 antibody (Upstate Biotechnology Inc) or WT1 antibody (a mixture of C19 and N18 polyclonal Abs, Santa Cruz Biotechnology) or non-immune IgG. The antibody/protein/DNA complex was collected by incubation with Protein G Agarose and washed in increasing salt buffers, then rinsed in TE as recommended (Upstate Biotechnology Inc). The complexes were recovered from agarose beads with an elution buffer, crosslinks were reversed and DNA was purified using G-50 spin columns. Four percent of both immunoprecipitated and input chromatin were amplified by PCR using Taq polymerase (Applied Biosystems by Roche Molecular System, Inc) and the following set of primers: VEGF primers (F) 5'TTCCTAGCAAAGAGGGAACG3' and (R) 5'ACCAAGGTTCACAGCCTGAA3'; AR primers (F) 5'TATCTGCTGGCTTGGTCATGGCTTG3' and (R) 5'CTGCTTCCTGAATAGCTCCTGCTT3'; and PSA primers (F) 5'TCTGCCTTTGTCCCCTAGAT3' and (R) 5'AACCTTCATTCCCCAGGACT3'. Following an initial 10 min denaturation at 95°C, DNA was amplified by 32 cycles of: 1) 20 sec denaturation at 95°C, 2) 30 sec annealing at either 58°C (for VEGF primers) or 59°C (for AR and PSA primers) and 3) 30 sec extension at 72°C; amplification was completed with a 2 min final extension at 72°C. PCR products were fractionated on 1% agarose gel, and ethidium bromide stained DNA was visualized by a gel doc system (BIORAD, CA). Specificity controls are shown in Additional file 2.
- WT1 :
Wilms tumor 1
transcription factor binding sites
- AR :
- VEGF :
vascular endothelial growth factor
- PSA :
prostate specific antigen
- KLK3 :
kallikrein-related peptidase 3
polymerase chain reaction
green fluorescent protein.
Jemal A, Siegel R, Ward E, Murray T, Xu J, Thun MJ: Cancer statistics, 2007. CA Cancer J Clin. 2007, 57: 43-66.
Stamey TA, Caldwell M, McNeal JE, Nolley R, Hemenez M, Downs J: The prostate specific antigen era in the United States is over for prostate cancer: what happened in the last 20 years?. J Urol. 2004, 172: 1297-1301. 10.1097/01.ju.0000139993.51181.5d.
Luo JH, Yu YP, Cieply K, Lin F, Deflavia P, Dhir R, Finkelstein S, Michalopoulos G, Becich M: Gene expression analysis of prostate cancers. Mol Carcinog. 2002, 33: 25-35. 10.1002/mc.10018.
Singh D, Febbo PG, Ross K, Jackson DG, Manola J, Ladd C, Tamayo P, Renshaw AA, D'Amico AV, Richie JP, Lander ES, Loda M, Kantoff PW, Golub TR, Sellers WR: Gene expression correlates of clinical prostate cancer behavior. Cancer Cell. 2002, 1: 203-209. 10.1016/S1535-6108(02)00030-2.
Chandran UR, Dhir R, Ma C, Michalopoulos G, Becich M, Gilbertson J: Differences in gene expression in prostate cancer, normal appearing prostate tissue adjacent to cancer and prostate tissue from cancer free organ donors. BMC Cancer. 2005, 5: 45-10.1186/1471-2407-5-45.
Brown K: Differential Gene Expression Patterns in Prostate Cancer Epithelial and Interstitial Stromal Cells. 2006, MS Thesis, Kent State University
Liu C, Calogero A, Ragona G, Adamson E, Mercola D: EGR-1, the reluctant suppression factor: EGR-1 is known to function in the regulation of growth, differentiation, and also has significant tumor suppressor activity and a mechanism involving the induction of TGF-beta1 is postulated to account for this suppressor activity. Crit Rev Oncog. 1996, 7 (1-2): 101-125.
Madden SL, Cook DM, Morris JF, Gashler A, Sukhatme VP, Rauscher FJ: Transcriptional repression mediated by the WT1 Wilms tumor gene product. Science. 1991, 253: 1550-1553. 10.1126/science.1654597.
Dey BR, Sukhatme VP, Roberts AB, Sporn MB, Rauscher FJ, Kim SJ: Repression of the transforming growth factor-beta 1 gene by the Wilms' tumor suppressor WT1 gene product. Mol Endocrinol. 1994, 8: 595-602. 10.1210/me.8.5.595.
Devilard E, Bladou F, Ramuz O, Karsenty G, Dales JP, Gravis G, Nguyen C, Bertucci F, Xerri L, Birnbaum D: FGFR1 and WT1 are markers of human prostate cancer progression. BMC Cancer. 2006, 6: 272-10.1186/1471-2407-6-272.
Ogishima T, Shiina H, Breault JE, Tabatabai L, Bassett WW, Enokida H, Li LC, Kawakami T, Urakami S, Ribeiro-Filho LA, Terashima M, Fujime M, Igawa M, Dahiya R: Increased heparanase expression is caused by promoter hypomethylation and up-regulation of transcriptional factor early growth response-1 in human prostate cancer. Clin Cancer Res. 2005, 11 (3): 1028-1036.
Eid MA, Kumar MV, Iczkowski KA, Bostwick DG, Tindall DJ: Expression of early growth response genes in human prostate cancer. Cancer Res. 1998, 58: 2461-2468.
Rivera MN, Haber DA: Wilms' tumour: connecting tumorigenesis and organ development in the kidney. Nat Rev Cancer. 2005, 5: 699-712. 10.1038/nrc1696.
Shimamura R, Fraizer GC, Trapman J, Lau Y, Saunders GF: The Wilms' tumor gene WT1 can regulate genes involved in sex determination and differentiation: SRY, Mullerian-inhibiting substance, and the androgen receptor. Clin Cancer Res. 1997, 3 (12 Pt 2): 2571-2580.
Hanson J, Gorman J, Reese J, Fraizer G: Regulation of vascular endothelial growth factor, VEGF, gene promoter by the tumor suppressor, WT1. Front Biosci. 2007, 12: 2279-2290. 10.2741/2230.
Thomas JW, Touchman JW: Vertebrate genome sequencing: building a backbone for comparative genomics. Trends in Genetics. 2002, 18: 104-108. 10.1016/S0168-9525(02)02599-4.
Dubchak I, Brudno M, Loots GG, Pachter L, Mayor C, Rubin EM, Frazer KA: Active conservation of noncoding sequences revealed by three-way species comparisons. Genome Res. 2000, 10: 1304-1306. 10.1101/gr.142200.
Stojanovic N, Florea L, Riemer C, Gumucio D, Slightom J, Goodman M, Miller W, Hardison R, Journals O: Comparison of five methods for finding conserved sequences in multiple alignments of gene regulatory regions. Nucleic Acids Res. 1999, 27: 3899-3910. 10.1093/nar/27.19.3899.
Wang H, Zhang Y, Cheng Y, Zhou Y, King DC, Taylor J, Chiaromonte F, Kasturi J, Petrykowska H, Gibb B, Dorman C, Miller W, Dore LC, Welch J, Weiss MJ, Hardison RC: Experimental validation of predicted mammalian erythroid cis-regulatory modules. Genome Res. 2006, 16: 1480-1492. 10.1101/gr.5353806.
Johnson DS, Davidson B, Brown CD, Smith WC, Sidow A: Noncoding regulatory sequences of Ciona exhibit strong correspondence between evolutionary constraint and functional importance. Genome Res. 2004, 14: 2448-2456. 10.1101/gr.2964504.
Hedges S, Kumar S: Genomic clocks and evolutionary timescales. Trends in Genetics. 2003, 19: 200-206. 10.1016/S0168-9525(03)00053-2.
Fraizer G, Leahy R, Priyadarshini S, Graham K, Delacerda J, Diaz M: Suppression of prostate tumor cell growth in vivo by WT1, the Wilms' tumor suppressor gene. Int J Oncol. 2004, 24: 461-471.
Kohler B, Delezoide AL, Boizet-Bonhoure B, McPhaul MJ, Sultan C, Lumbroso S: Coexpression of Wilms' tumor suppressor 1 (WT1) and androgen receptor (AR) in the genital tract of human male embryos and regulation of AR promoter activity by WT1. J Mol Endocrinol. 2007, 38: 547-554. 10.1677/JME-06-0020.
Husbeck B, Bhattacharyya RS, Feldman D, Knox SJ: Inhibition of androgen receptor signaling by selenite and methylseleninic acid in prostate cancer cells: two distinct mechanisms of action. Mol Cancer Ther. 2006, 5: 2078-2085. 10.1158/1535-7163.MCT-06-0056.
Pore N, Gupta AK, Cerniglia GJ, Jiang Z, Bernhard EJ, Evans SM, Koch CJ, Hahn SM, Maity A: Nelfinavir down-regulates hypoxia-inducible factor 1alpha and VEGF expression and increases tumor oxygenation: implications for radiotherapy. Cancer Res. 2006, 66: 9252-9259. 10.1158/0008-5472.CAN-06-1239.
Englert C, Hou X, Maheswaran S, Bennett P, Ngwu C, Re GG, Garvin AJ, Rosner MR, Haber DA: WT1 suppresses synthesis of the epidermal growth factor receptor and induces apoptosis. EMBO J. 1995, 14 (19): 4662-4675.
Tagle DA, Koop BF, Goodman M, Slightom JL, Hess DL, Jones RT: Embryonic and globin genes of a prosimian primate (Galago crassicaudatus): Nucleotide and amino acid sequences, developmental regulation and phylogenetic footprints. J Mol Biol. 1988, 203: 439-455. 10.1016/0022-2836(88)90011-3.
Elnitski L, Miller W, Hardison R: Conserved E boxes function as part of the enhancer in hypersensitive site 2 of the beta-globin locus control region. Role of basic helix-loop-helix proteins. J Biol Chem. 1997, 272: 369-378. 10.1074/jbc.272.1.369.
Hardison RC, Oeltjen J, Miller W: Long human-mouse sequence alignments reveal novel regulatory elements: a reason to sequence the mouse genome. Genome Res. 1997, 7 (10): 959-966.
Vuillaumier S, Dixmeras I, Messai H, Lapoumeroulie C, Lallemand D, Gekas J, Chehab FF, Perret C, Elion J, Denamur E: Cross-species characterization of the promoter region of the cystic fibrosis transmembrane conductance regulator gene reveals multiple levels of regulation. Biochem J. 1997, 327 (Pt 3): 651-662.
Loots GG, Locksley RM, Blankespoor CM, Wang ZE, Miller W, Rubin EM, Frazer KA: Identification of a coordinate regulator of interleukins 4, 13, and 5 by cross-species sequence comparisons. Science. 2000, 288: 136-140. 10.1126/science.288.5463.136.
Cliften P, Sudarsanam P, Desikan A, Fulton L, Fulton B, Majors J, Waterston R, Cohen BA, Johnston M: Finding functional features in saccharomyces genomes by phylogenetic footprinting. Science. 2003, 301: 71-76. 10.1126/science.1084337.
Kellis M, Patterson N, Endrizzi M, Birren B, Lander ES: Sequencing and comparison of yeast species to identify genes and regulatory elements. Nature. 2003, 423: 241-254. 10.1038/nature01644.
Xie X, Lu J, Kulbokas EJ, Golub TR, Mootha V, Lindblad-Toh K, Lander ES, Kellis M: Systematic discovery of regulatory motifs in human promoters and 3' UTRs by comparison of several mammals. Nature. 2005, 434: 338-345. 10.1038/nature03441.
Hughes JR, Cheng JF, Ventress N, Prabhakar S, Clark K, Anguita E, De Gobbi M, de Jong P, Rubin E, Higgs DR: Annotation of cis-regulatory elements by identification, subclassification, and functional assessment of multispecies conserved sequences. Proc Natl Acad Sci USA. 2005, 102: 9830-10.1073/pnas.0503401102.
King DC, Taylor J, Elnitski L, Chiaromonte F, Miller W, Hardison RC: Evaluation of regulatory potential and conservation scores for detecting cis-regulatory modules in aligned mammalian genome sequences. Genome Res. 2005, 15: 1051-1060. 10.1101/gr.3642605.
Frith MC, Ponjavic J, Fredman D, Kai C, Kawai J, Carninci P, Hayshizaki Y, Sandelin A: Evolutionary turnover of mammalian transcription start sites. Genome Res. 2006, 16: 713-722. 10.1101/gr.5031006.
Loots GG, Ovcharenko I: rVISTA 2.0: evolutionary analysis of transcription factor binding sites. Nucleic Acids Res. 2004, 32: W217-221. 10.1093/nar/gkh383.
Ovcharenko I, Loots GG, Giardine BM, Hou M, Ma J, Hardison RC, Stubbs L, Miller W: Mulan: Multiple-sequence local alignment and visualization for studying function and evolution. Genome Res. 2005, 15: 184-194. 10.1101/gr.3007205.
Dermitzakis ET, Clark AG: Evolution of transcription factor binding sites in mammalian gene regulatory regions: conservation and turnover. Mol Biol Evol. 2002, 19: 1114-1121.
Hewitt SM, Hamada S, McDonnell TJ, Rauscher FJ, Saunders GF: Regulation of the proto-oncogenes bcl-2 and c-myc by the Wilms' tumor suppressor gene WT1. Cancer Res. 1995, 55 (22): 5386-5389.
Mayo MW, Wang CY, Drouin SS, Madrid LV, Marshall AF, Reed JC, Weissman BE, Baldwin AS: WT1 modulates apoptosis by transcriptionally upregulating the bcl-2 proto-oncogene. EMBO J. 1999, 18: 3990-4003. 10.1093/emboj/18.14.3990.
Werner H, Rauscher FJ, Sukhatme VP, Drummond IA, Roberts CT, LeRoith D: Transcriptional repression of the insulin-like growth factor I receptor (IGF-I-R) gene by the tumor suppressor WT1 involves binding to sequences both upstream and downstream of the IGF-I-R gene transcription start site. J Biol Chem. 1994, 269 (17): 12577-12582.
Werner H, Re GG, Drummond IA, Sukhatme VP, Rauscher FJ, Sens DA, Garvin AJ, LeRoith D, Roberts CT: Increased expression of the insulin-like growth factor I receptor gene, IGF1R, in Wilms tumor is correlated with modulation of IGF1R promoter activity by the WT1 Wilms tumor gene product. Proc Natl Acad Sci U S A. 1993, 90: 5828-5832. 10.1073/pnas.90.12.5828.
Drummond IA, Madden SL, Rohwer-Nutter P, Bell GI, Sukhatme VP, Rauscher FJ: Repression of the insulin-like growth factor II gene by the Wilms tumor suppressor WT1. Science. 1992, 257: 674-678. 10.1126/science.1323141.
Hosono S, Gross I, English MA, Hajra KM, Fearon ER, Licht JD: E-cadherin is a WT1 target gene. J Biol Chem. 2000, 275: 10943-10953. 10.1074/jbc.275.15.10943.
Hewitt SM, Fraizer GC, Wu YJ, Rauscher FJ, Saunders GF: Differential function of Wilms' tumor gene WT1 splice isoforms in transcriptional regulation. J Biol Chem. 1996, 271: 8588-8592. 10.1074/jbc.271.15.8588.
Yin H, Radomska HS, Tenen DG, Glass J: Down regulation of PSA by C/EBPalpha is associated with loss of AR expression and inhibition of PSA promoter activity in the LNCaP cell line. BMC Cancer. 2006, 6: 158-10.1186/1471-2407-6-158.
Borgono CA, Diamandis EP: The emerging roles of human tissue kallikreins in cancer. Nat Rev Cancer. 2004, 4: 876-890. 10.1038/nrc1474.
Merika M, Orkin SH: DNA-binding specificity of GATA family transcription factors. Mol Cell Biol. 1993, 13: 3999-4010.
Tremblay JJ, Viger RS: Novel roles for GATA transcription factors in the regulation of steroidogenesis. J Steroid Biochem Mol Biol. 2003, 85: 291-298. 10.1016/S0960-0760(03)00211-5.
Viger RS, Guittot SM, Anttonen M, Wilson DB, Heikinheimo M: Role of the GATA family of transcription factors in endocrine development, function, and disease. Mol Endocrinol. 2008, 22: 781-798. 10.1210/me.2007-0513.
Fraizer GC, Wu YJ, Hewitt SM, Maity T, Ton CC, Huff V, Saunders GF: Transcriptional regulation of the human Wilms' tumor gene (WT1). cell type-specific enhancer and promiscuous promoter. J Biol Chem. 1994, 269 (12): 8892-8900.
Wu Y, Fraizer GC, Saunders GF: GATA-1 Transactivates the WT1 Hematopoietic Specific Enhancer. J Biol Chem. 1995, 270: 5944-5949. 10.1074/jbc.270.11.5944.
Zhang X, Xing G, Fraizer GC, Saunders GF: Transactivation of an Intronic Hematopoietic-specific Enhancer of the Human Wilms' Tumor 1 Gene by GATA-1 and c-Myb. J Biol Chem. 1997, 272: 29272-29280. 10.1074/jbc.272.46.29272.
Perez-Stable CM, Pozas A, Roos BA: A role for GATA transcription factors in the androgen regulation of the prostate-specific antigen gene enhancer. Mol Cell Endocrinol. 2000, 167: 43-53. 10.1016/S0303-7207(00)00300-2.
Miyamoto Y, Silversides D, Viger R: GATA4 Enhances SRY Gene Transcription Through a Direct Interaction with Wilms Tumor 1 (WT1). Biology of Reproduction. 2005, 97-[http://abstracts.co.allenpress.com/pweb/ssr2005/document/47755]
Rauscher FJ, Morris JF, Tournay OE, Cook DM, Curran T: Binding of the Wilms' tumor locus zinc finger protein to the EGR-1 consensus sequence. Science. 1990, 250: 1259-1262. 10.1126/science.2244209.
Larsson SH, Charlieu JP, Miyagawa K, Engelkamp D, Rassoulzadegan M, Ross A, Cuzin F, van Heyningen V, Hastie ND: Subnuclear localization of WT1 in splicing or transcription factor domains is regulated by alternative splicing. Cell. 1995, 81: 391-401. 10.1016/0092-8674(95)90392-5.
Ladomery M, Sommerville J, Woolner S, Slight J, Hastie N: Expression in Xenopus oocytes shows that WT1 binds transcripts in vivo, with a central role for zinc finger one. J Cell Sci. 2003, 116: 1539-1549. 10.1242/jcs.00324.
Abdulkadir SA, Qu Z, Garabedian E, Song SK, Peters TJ, Svaren J, Carbone JM, Naughton CK, Catalona WJ, Ackerman JJ, Gordon JI, Humphrey PA, Milbrandt J: Impaired prostate tumorigenesis in Egr1-deficient mice. Nat Med. 2001, 7: 101-107. 10.1038/83231.
Ellen T, Ke Q, Zhang P, Costa M: NDRG1, a Growth and Cancer Related Gene: Regulation of Gene Expression and Function in Normal and Disease States. Carcinogenesis. 2008, 29: 2-8. 10.1093/carcin/bgm200.
Caruso RP, Levinson B, Melamed J, Wieczorek R, Taneja S, Polsky D, Chang C, Zeleniuch-Jacquotte A, Salnikow K, Yee H, Costa M, Osman I: Altered N-myc downstream-regulated gene 1 protein expression in African-American compared with caucasian prostate cancer patients. Clin Cancer Res. 2004, 10: 222-227. 10.1158/1078-0432.CCR-0604-3.
Rong Y, Hu F, Huang R, Mackman N, Horowitz JM, Jensen RL, Durden DL, Van Meir EG, Brat DJ: Early growth response gene-1 regulates hypoxia-induced expression of tissue factor in glioblastoma multiforme through hypoxia-inducible factor-1-independent mechanisms. Cancer Res. 2006, 66: 7067-7074. 10.1158/0008-5472.CAN-06-0346.
Zhang P, Tchou-Wong KM, Costa M: Egr-1 Mediates Hypoxia-Inducible Transcription of the NDRG1 Gene through an Overlapping Egr-1/Sp1 Binding Site in the Promoter. Cancer Res. 2007, 67: 9125-10.1158/0008-5472.CAN-07-1525.
Minc E, de Coppet P, Masson P, Thiery L, Dutertre S, Amor-Gueret M, Jaulin C: The human copper-zinc superoxide dismutase gene (SOD1) proximal promoter is regulated by Sp1, Egr-1, and WT1 via non-canonical binding sites. J Biol Chem. 1999, 274: 503-509. 10.1074/jbc.274.1.503.
Harrington MA, Konicek B, Song A, Xia XL, Fredericks WJ, Rauscher FJ: Inhibition of colony-stimulating factor-1 promoter activity by the product of the Wilms' tumor locus. J Biol Chem. 1993, 268 (28): 21271-21275.
Wang ZY, Madden SL, Deuel TF, Rauscher FJ: The Wilms' tumor gene product, WT1, represses transcription of the platelet-derived growth factor A-chain gene. J Biol Chem. 1992, 267: 21999-22002.
Stoner M, Wormke M, Saville B, Samudio I, Qin C, Abdelrahim M, Safe S: Estrogen regulation of vascular endothelial growth factor gene expression in ZR-75 breast cancer cells through interaction of estrogen receptor alpha and SP proteins. Oncogene. 2004, 23: 1052-1063. 10.1038/sj.onc.1207201.
Reizner N, Maor S, Sarfstein R, Abramovitch S, Welshons WV, Curran EM, Lee AV, Werner H: The WT1 Wilms' tumor suppressor gene product interacts with estrogen receptor-alpha and regulates IGF-I receptor gene transcription in breast cancer cells. J Mol Endocrinol. 2005, 35: 135-144. 10.1677/jme.1.01761.
Ensembl Genome Browser. [http://www.ensembl.org/]
Hubbard TJ, Aken BL, Beal K, Ballester B, Caccamo M, Chen Y, Clarke L, Coates G, Cunningham F, Cutts T, Down T, Dyer SC, Fitzgerald S, Fernandez-Banet J, Graf S, Haider S, Hammond M, Herrero J, Holland R, Howe K, Howe K, Johnson N, Kahari A, Keefe D, Kokocinski F, Kulesha E, Lawson D, Longden I, Melsopp C, Megy K, Meidl P, Ouverdin B, Parker A, Prlic A, Rice S, Rios D, Schuster M, Sealy I, Severin J, Slater G, Smedley D, Spudich G, Trevanion S, Vilella A, Vogel J, White S, Wood M, Cox T, Curwen V, Durbin R, Fernandez-Suarez XM, Flicek P, Kasprzyk A, Proctor G, Searle S, Smith J, Ureta-Vidal A, Birney E: Ensembl 2007. Nucleic Acids Res. 2007, 35: D610-7. 10.1093/nar/gkl996.
Dieterich C, Cusack B, Wang H, Rateitschak K, Krause A, Vingron M: Annotating regulatory DNA based on man-mouse genomic comparison. Bioinformatics. 2002, 18 (Suppl 2): S84-90.
Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, Rajandream MA, Barrell B: Artemis: sequence visualization and annotation. Bioinformatics. 2000, 16: 944-945. 10.1093/bioinformatics/16.10.944.
Cartharius K, Frech K, Grote K, Klocke B, Haltmeier M, Klingenhoff A, Frisch M, Bayerlein M, Werner T: MatInspector and beyond: promoter analysis based on transcription factor binding sites. Bioinformatics. 2005, 21: 2933-2942. 10.1093/bioinformatics/bti473.
Wingender E, Chen X, Hehl R, Karas H, Liebich I, Matys V, Meinhardt T, Prüß M, Reuter I, Schacherer F: TRANSFAC: an integrated system for gene expression regulation. Nucleic Acids Res. 2000, 28: 316-319. 10.1093/nar/28.1.316.
Werner T: Computer-assisted analysis of transcription control regions. Methods Mol Biol. 2000, 132: 337-349.
Genomatix-understanding gene regulation. [http://www.genomatix.de]
Schwartz S, Elnitski L, Li M, Weirauch M, Riemer C, Smit A, NISC Comparative Sequencing Program, Green ED, Hardison RC, Miller W: MultiPipMaker and supporting tools: alignments and analysis of multiple genomic DNA sequences. Nucleic Acids Res. 2003, 31: 3518-3524. 10.1093/nar/gkg579.
Frazer KA, Tao H, Osoegawa K, de Jong PJ, Chen X, Doherty MF, Cox DR: Noncoding Sequences Conserved in a Limited Number of Mammals in the SIM2 Interval are Frequently Functional. Genome Res. 2004, 14: 367-372. 10.1101/gr.1961204.
Berquin IM, Min Y, Wu R, Wu H, Chen YQ: Expression signature of the mouse prostate. J Biol Chem. 2005, 280: 36442-36451. 10.1074/jbc.M504945200.
Dutton JR, Lahiri D, Ward A: Different isoforms of the Wilms' tumour protein WT1 have distinct patterns of distribution and trafficking within the nucleus. Cell Prolif. 2006, 39: 519-535. 10.1111/j.1365-2184.2006.00409.x.
Weinmann AS, Farnham PJ: Identification of unknown target genes of human transcription factors using chromatin immunoprecipitation. Methods. 2002, 26: 37-47. 10.1016/S1046-2023(02)00006-3.
Hanson J, Brown K, Reese J, Gorman J, Cash J, Graham K, Fraizer GC: Wilms' tumor suppressor gene, WT1, modulates VEGF expression and is differentially expressed in prostate tumor cells. Proc Amer Assoc Cancer Res. 2006, 47: 2492-[http://www.aacrmeetingabstracts.org/cgi/content/abstract/2006/1/590]
Sroka IC, Nagle RB, Bowden GT: Membrane-type 1 matrix metalloproteinase is regulated by sp1 through the differential activation of AKT, JNK, and ERK pathways in human prostate tumor cells. Neoplasia. 2007, 9: 406-417. 10.1593/neo.07193.
Phan D, Cheng CJ, Galfione M, Vakar-Lopez F, Tunstead J, Thompson NE, Burgess RR, Najjar SM, Yu-Lee LY, Lin SH: Identification of Sp2 as a Transcriptional Repressor of Carcinoembryonic Antigen-Related Cell Adhesion Molecule 1 in Tumorigenesis. Cancer Res. 2004, 64: 3072-3078. 10.1158/0008-5472.CAN-03-3730.
Ho SM, Leung YK, Chung I: Estrogens and antiestrogens as etiological factors and therapeutics for prostate cancer. Ann N Y Acad Sci. 2006, 1089: 177-193. 10.1196/annals.1386.005.
Ruiz M, Pettaway C, Song R, Stoeltzing O, Ellis L, Bar-Eli M: Activator protein 2alpha inhibits tumorigenicity and represses vascular endothelial growth factor transcription in prostate cancer cells. Cancer Res. 2004, 64: 631-638. 10.1158/0008-5472.CAN-03-2751.
Zhang X, Leung YK, Ho SM: AP-2 regulates the transcription of estrogen receptor (ER)-beta by acting through a methylation hotspot of the oN promoter in prostate cancer cells. Oncogene. 2007, 26: 7346-7354. 10.1038/sj.onc.1210537.
Rhodes DR, Sanda MG, Otte AP, Chinnaiyan AM, Rubin MA: Multiplex biomarker approach for determining risk of prostate-specific antigen-defined recurrence of prostate cancer. J Natl Cancer Inst. 2003, 95: 661-668.
Mehrian-Shai R, Chen CD, Shi T, Horvath S, Nelson SF, Reichardt JK, Sawyers CL: Insulin growth factor-binding protein 2 is a candidate biomarker for PTEN status and PI3K/Akt pathway activation in glioblastoma and prostate cancer. Proc Natl Acad Sci USA. 2007, 104: 5563-5568. 10.1073/pnas.0609139104.
Chatterjee S, Park ES, Soloff MS: Proliferation of DU145 prostate cancer cells is inhibited by suppressing insulin-like growth factor binding protein-2. Int J Urol. 2004, 11: 876-884. 10.1111/j.1442-2042.2004.00898.x.
Le H, Arnold JT, McFann KK, Blackman MR: DHT and testosterone, but not DHEA or E2, differentially modulate IGF-I, IGFBP-2, and IGFBP-3 in human prostatic stromal cells. Am J Physiol Endocrinol Metab. 2006, 290: E952-E960. 10.1152/ajpendo.00451.2005.
Nagpal ML, Davis J, Lin T: Overexpression of CXCL10 in human prostate LNCaP cells activates its receptor (CXCR3) expression and inhibits cell proliferation. Biochim Biophys Acta. 2006, 1762: 811-818.
Tu LC, Yan X, Hood L, Lin B: Proteomics analysis of the interactome of N-myc downstream regulated gene 1 and its interactions with the androgen response program in prostate cancer cells. Mol Cell Proteomics. 2007, 6: 575-588. 10.1074/mcp.M600249-MCP200.
Kim KH, Dobi A, Shaheduzzaman S, Gao CL, Masuda K, Li H, Drukier A, Gu Y, Srikantan V, Rhim JS, Srivastava S: Characterization of the androgen receptor in a benign prostate tissue-derived human prostate epithelial cell line: RC-165N/human telomerase reverse transcriptase. Prostate Cancer Prostatic Dis. 2007, 10: 30-38. 10.1038/sj.pcan.4500915.
Ruscica M, Dozio E, Boghossian S, Bovo G, Martos Riano V, Motta M, Magni P: Activation of the Y1 receptor by neuropeptide Y regulates the growth of prostate cancer cells. Endocrinology. 2006, 147: 1466-1473. 10.1210/en.2005-0925.
Rasiah KK, Kench JG, Gardiner-Garden M, Biankin AV, Golovsky D, Brenner PC, Kooner R, O'neill GF, Turner JJ, Delprado W, Lee CS, Brown DA, Breit SN, Grygiel JJ, Horvath LG, Stricker PD, Sutherland RL, Henshall SM: Aberrant neuropeptide Y and macrophage inhibitory cytokine-1 expression are early events in prostate cancer development and are associated with poor prognosis. Cancer Epidemiol Biomarkers Prev. 2006, 15: 711-716. 10.1158/1055-9965.EPI-05-0752.
Liu P, Ramachandran S, Ali Seyed M, Scharer CD, Laycock N, Dalton WB, Williams H, Karanam S, Datta MW, Jaye DL, Moreno CS: Sex-determining region Y box 4 is a transforming oncogene in human prostate cancer cells. Cancer Res. 2006, 66: 4011-4019. 10.1158/0008-5472.CAN-05-3055.
Vanaja DK, Ballman KV, Morlan BW, Cheville JC, Neumann RM, Lieber MM, Tindall DJ, Young CY: PDLIM4 repression by hypermethylation as a potential biomarker for prostate cancer. Clin Cancer Res. 2006, 12: 1128-1136. 10.1158/1078-0432.CCR-05-2072.
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. 2007, 67: 528-536. 10.1158/0008-5472.CAN-06-1672.
Drivdahl R, Haugk KH, Sprenger CC, Nelson PS, Tennant MK, Plymate SR: Suppression of growth and tumorigenicity in the prostate tumor cell line M12 by overexpression of the transcription factor SOX9. Oncogene. 2004, 23: 4584-4593. 10.1038/sj.onc.1207603.
ExPASy-UniProt Knowledgebase: Swiss-Prot and TrEMBL. [http://ca.expasy.org/sprot/]
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556.
Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4: P3-10.1186/gb-2003-4-5-p3.
This research was supported by NIH-1CA331160 (GF). We gratefully acknowledge Dr. A. Ward for the GFP-WT1 fusion construct  and Kimberly Bean for transfected cells.
KE performed ChIP analysis, in silico analysis, and drafted the text. ST performed in silico analysis. ABazarov performed in silico analysis. ABrett performed ChIP analysis and contributed to the text. GF planned functional assays, guided student research, and drafted the text. HP planned in silico analyses, guided student research, and drafted the text.
Electronic supplementary material
Additional file 1: Evolutionary conserved TFBS in promoters of 24 genes expressed in prostate cancer epithelium. This table lists evolutionary conserved transcription factor binding sites in promoters of 24 genes expressed in prostate cancer. (PDF 106 KB)
Additional file 2: WT1 bound the proximal, but not the distal, region of the amphiregulin (AREG) gene promoter in chromatin of LNCaP cells. This specificity control illustrates WT1 binding to the proximal region (known to bind WT1), but not the distal region of the AREG promoter in LNCaP chromatin. (PDF 74 KB)
About this article
Cite this article
Eisermann, K., Tandon, S., Bazarov, A. et al. Evolutionary conservation of zinc finger transcription factor binding sites in promoters of genes co-expressed with WT1 in prostate cancer. BMC Genomics 9, 337 (2008). https://doi.org/10.1186/1471-2164-9-337
- Androgen Receptor
- Prostate Specific Antigen
- Transcription Factor Binding Site
- Vascular Endothelial Growth Factor Promoter
- Transcription Factor Binding Site Prediction