Integrative genomics of gene and metabolic regulation by estrogen receptors α and ß, and their coregulators.

The closely related transcription factors (TFs), estrogen receptors ERα and ERß, regulate divergent gene expression programs and proliferative outcomes in breast cancer. Utilizing breast cancer cells with ERα, ERß, or both receptors as a model system to define the basis for differing response specification by related TFs, we show that these TFs and their key coregulators, SRC3 and RIP140, generate overlapping as well as unique chromatin-binding and transcription-regulating modules. Cistrome and transcriptome analyses and the use of clustering algorithms delineated 11 clusters representing different chromatin-bound receptor and coregulator assemblies that could be functionally associated through enrichment analysis with distinct patterns of gene regulation and preferential coregulator usage, RIP140 with ERß and SRC3 with ERα. The receptors modified each other's transcriptional effect, and ERß countered the proliferative drive of ERα through several novel mechanisms associated with specific binding-site clusters. Our findings delineate distinct TF-coregulator assemblies that function as control nodes, specifying precise patterns of gene regulation, proliferation, and metabolism, as exemplified by two of the most important nuclear hormone receptors in human breast cancer.

Genetic variation of ESR1 and its co-activator PPARGC1B is synergistic in augmenting the risk of estrogen receptor-positive breast cancer.

INTRODUCTION: Given the role of estrogen in breast carcinogenesis and the modification of estrogen receptor (ER) activity by its biochemical cofactors, we hypothesize that genetic variation within ER cofactor genes alters cellular response to estrogen exposure and consequently modifies the risk for ER-positive breast cancer. METHODS: We genotyped 790 tagging SNPs within 60 ER cofactor genes in 1,257 cases and 1,464 controls from Sweden and in 2,215 cases and 1,265 controls from Finland, and tested their associations with either ER-positive or ER-negative breast cancer. RESULTS: Seven SNPs showed consistent association with ER-positive breast cancer in the two independent samples, and six of them were located within PPARGC1B, encoding an ER co-activator, with the strongest association at rs741581 (odds ratio = 1.41, P = 4.84 × 10⁻5) that survived Bonferroni correction for multiple testing in the combined ER-positive breast cancer sample (Pcorrected = 0.03). Moreover, we also observed significant synergistic interaction (Pinteraction = 0.008) between the genetic polymorphisms within PPARGC1B and ESR1 in ER-positive breast cancer. By contrast, no consistent association was observed in ER-negative breast cancer. Furthermore, we found that administration of estrogen in the MCF-7 cell line induced PPARGC1B expression and enhanced occupancies of ER and RNA polymerase II within the region of SNP association, suggesting the upregulation of PPARGC1B expression by ESR1 activation. CONCLUSIONS: Our study revealed that DNA polymorphisms of PPARGC1B, coding a bona fide ER co-activator, are associated with ER-positive breast cancer risk. The feed-forward transcriptional regulatory loop between PPARGC1B and ESR1 further augments their protein interaction, which provides a plausible mechanistic explanation for the synergistic genetic interaction between PPARGC1B and ESR1 in ER-positive breast cancer. Our study also highlights that biochemically and genomically informed candidate gene studies can enhance the discovery of interactive disease susceptibility genes.

Whole-genome cartography of estrogen receptor alpha binding sites.

Using a chromatin immunoprecipitation-paired end diTag cloning and sequencing strategy, we mapped estrogen receptor alpha (ERalpha) binding sites in MCF-7 breast cancer cells. We identified 1,234 high confidence binding clusters of which 94% are projected to be bona fide ERalpha binding regions. Only 5% of the mapped estrogen receptor binding sites are located within 5 kb upstream of the transcriptional start sites of adjacent genes, regions containing the proximal promoters, whereas vast majority of the sites are mapped to intronic or distal locations (>5 kb from 5' and 3' ends of adjacent transcript), suggesting transcriptional regulatory mechanisms over significant physical distances. Of all the identified sites, 71% harbored putative full estrogen response elements (EREs), 25% bore ERE half sites, and only 4% had no recognizable ERE sequences. Genes in the vicinity of ERalpha binding sites were enriched for regulation by estradiol in MCF-7 cells, and their expression profiles in patient samples segregate ERalpha-positive from ERalpha-negative breast tumors. The expression dynamics of the genes adjacent to ERalpha binding sites suggest a direct induction of gene expression through binding to ERE-like sequences, whereas transcriptional repression by ERalpha appears to be through indirect mechanisms. Our analysis also indicates a number of candidate transcription factor binding sites adjacent to occupied EREs at frequencies much greater than by chance, including the previously reported FOXA1 sites, and demonstrate the potential involvement of one such putative adjacent factor, Sp1, in the global regulation of ERalpha target genes. Unexpectedly, we found that only 22%-24% of the bona fide human ERalpha binding sites were overlapping conserved regions in whole genome vertebrate alignments, which suggest limited conservation of functional binding sites. Taken together, this genome-scale analysis suggests complex but definable rules governing ERalpha binding and gene regulation.

Estrogen Receptors alpha and beta as determinants of gene expression: influence of ligand, dose, and chromatin binding.

Estrogen receptors alpha and beta (ERalpha and ERbeta) mediate the actions of estrogens in a variety of normal and cancer target cells. Estrogens differ in their preference for these ERs, and many phytoestrogens bind preferentially to ERbeta. To investigate how phytoestrogens such as genistein impact ER-regulated gene expression, we used adenoviral gene delivery of ERbeta coupled with ERalpha depletion with small interfering RNA to generate human breast cancer (MCF-7) cells expressing four complements of ERalpha and ERbeta. We examined the dose-dependent effects of genistein on genome-wide gene expression by DNA microarrays and monitored the recruitment of ERs and coregulators to responsive regions of estrogen-regulated genes. At a low (6 nm) concentration, genistein regulated gene expression much more effectively in cells coexpressing ERalpha and ERbeta than in cells expressing ERalpha alone, whereas at high concentration (300 nm), genistein induced transcriptome changes very similar to that of 17beta-estradiol. We demonstrate that ERbeta is preferentially activated by genistein and is recruited to estrogen-responsive genomic sites and that differential occupancy of ERalpha and ERbeta by genistein and 17beta-estradiol in turn influences the recruitment patterns of coregulators such as steroid receptor coactivator 3 (SRC3) and receptor-interacting protein 140 (RIP140). Our observations indicate that genistein is a potency-selective ligand for gene expression regulation by ERalpha and ERbeta and that the ability of ERalpha and ERbeta to serve as determinants of gene expression is greatly influenced by the nature of the ligand, by ligand dose, and by the differential abilities of ligand-ER complexes to recruit different coregulators at ER binding sites of hormone-regulated genes.

The estrogen-regulated transcription factor PITX1 coordinates gene-specific regulation by estrogen receptor-alpha in breast cancer cells.

The estrogen receptor α (ERα) is a master regulator of gene expression and works along with cooperating transcription factors in mediating the actions of the hormone estradiol (E2) in ER-positive tissues and breast tumors. Here, we report that expression of paired-like homeodomain transcription factor (PITX1), a tumor suppressor and member of the homeobox family of transcription factors, is robustly up-regulated by E2 in several ERα-positive breast cancer cell lines via ERα-dependent interaction between the proximal promoter and an enhancer region 5' upstream of the PITX1 gene. Overexpression of PITX1 selectively inhibited the transcriptional activity of ERα and ERß, while enhancing the activities of the glucocorticoid receptor and progesterone receptor. Reduction of PITX1 by small interfering RNA enhanced ERα-dependent transcriptional regulation of a subset of ERα target genes. The consensus PITX1 binding motif was found to be present in 28% of genome-wide ERα binding sites and was in close proximity to estrogen response elements in a subset of ERα binding sites, and E2 treatment enhanced PITX1 as well as ERα recruitment to these binding sites. These studies identify PITX1 as a new ERα transcriptional target that acts as a repressor to coordinate and fine tune target-specific, ERα-mediated transcriptional activity in human breast cancer cells.

Genome-wide dynamics of chromatin binding of estrogen receptors alpha and beta: mutual restriction and competitive site selection.

Estrogen receptors ERalpha and ERbeta, members of the nuclear receptor superfamily, exert profound effects on the gene expression and biological response programs of their target cells. Herein, we explore the dynamic interplay between these two receptors in their selection of chromatin binding sites when present separately or together in MCF-7 breast cancer cells. Treatment of cells (containing ERalpha only, ERbeta only, or ERalpha and ERbeta) with estradiol or ER subtype-selective ligands was followed by chromatin immunoprecipitation analysis with a custom-designed tiling array for ER binding sites across the genome to examine the effects of ligand-occupied and unoccupied ERalpha and ERbeta on chromatin binding. There was substantial overlap in binding sites for these estradiol-liganded nuclear receptors when present alone, but many fewer sites were shared when both ERs were present. Each ER restricted the binding site occupancy of the other, with ERalpha generally being dominant. Binding sites of both receptors were highly enriched in estrogen response element motifs, but when both ERs were present, ERalpha displaced ERbeta, shifting it into new sites less enriched in estrogen response elements. Binding regions of the two ERs also showed differences in their enrichments for other transcription factor binding motifs. Studies with ER subtype-specific ligands revealed that it was the liganded subtype that principally determined the spectrum of chromatin binding. These findings highlight the dynamic interplay between the two ERs in their selection of chromatin binding sites, with competition, restriction, and site shifting having important implications for the regulation of gene expression by these two nuclear receptors.

Genome-wide analysis of estrogen receptor alpha DNA binding and tethering mechanisms identifies Runx1 as a novel tethering factor in receptor-mediated transcriptional activation.

Nuclear receptor estrogen receptor alpha (ER alpha) controls the expression of hundreds of genes responsible for target cell phenotypic properties, but the relative importance of direct versus tethering mechanisms of DNA binding has not been established. In this first report, we examine the genome-wide chromatin localization of an altered-specificity mutant ER with a DNA binding domain deficient in binding to estrogen response element (ERE)-containing DNA (DBDmut ER) versus wild-type ER alpha. Using high-throughput sequencing of ER chromatin immunoprecipitations (ChIP-Seq) and mRNA transcriptional profiling, we show that direct ERE binding is required for most of (75%) estrogen-dependent gene regulation and 90% of hormone-dependent recruitment of ER to genomic binding sites. De novo motif analysis of the chromatin binding regions in MDA-MB-231 human breast cancer cells defined unique transcription factor profiles responsible for genes regulated through tethering versus direct ERE binding, with Runx motifs enriched in ER-tethered sites. We confirmed a role for Runx1 in mediating ER alpha genomic recruitment and regulation of tethering genes. Our findings delineate the contributions of direct receptor ERE binding versus binding through response elements for other transcription factors in chromatin localization and ER-dependent gene regulation, paradigms likely to underlie the gene regulatory actions of other nuclear receptors as well.

Estrogen receptor regulation of carbonic anhydrase XII through a distal enhancer in breast cancer.

The expression of carbonic anhydrase XII (CA12), a gene that encodes a zinc metalloenzyme responsible for acidification of the microenvironment of cancer cells, is highly correlated with estrogen receptor alpha (ER alpha) in human breast tumors. Here, we show that CA12 is robustly regulated by estrogen via ER alpha in breast cancer cells, and that this regulation involves a distal estrogen-responsive enhancer region. Upon the addition of estradiol, ER alpha binds directly to this distal enhancer in vivo, resulting in the recruitment of RNA polymerase II and steroid receptor coactivators SRC-2 and SRC-3, and changes in histone acetylation. Mutagenesis of an imperfect estrogen-responsive element within this enhancer region abolishes estrogen-dependent activity, and chromosome conformation capture and chromatin immunoprecipitation assays show that this distal enhancer communicates with the transcriptional start site of the CA12 gene via intrachromosomal looping upon hormone treatment. This distal enhancer element is observed in the homologous mouse genomic sequence, and the expression of the mouse homologue, Car12, is rapidly and robustly stimulated by estradiol in the mouse uterus in vivo, suggesting that the ER regulation of CA12 is mechanistically and evolutionarily conserved. Our findings highlight the crucial role of ER in the regulation of the CA12 gene, and provide insight into the transcriptional regulatory mechanism that accounts for the strong association of CA12 and ER in human breast cancers.