Sensitive ChIP-DSL technology reveals an extensive estrogen receptor alpha-binding program on human gene promoters

Leveraging sex-genetic interactions to understand brain disorders: recent advances and current gaps

It is established that there are sex differences in terms of prevalence, age of onset, clinical manifestations, and response to treatment for a variety of brain disorders, including neurodevelopmental, psychiatric, and neurodegenerative disorders. Cohorts of increasing sample sizes with diverse data types collected, including genetic, transcriptomic and/or phenotypic data, are providing the building blocks to permit analytical designs to test for sex-biased genetic variant-trait associations, and for sex-biased transcriptional regulation. Such molecular assessments can contribute to our understanding of the manifested phenotypic differences between the sexes for brain disorders, offering the future possibility of delivering personalized therapy for females and males. With the intention of raising the profile of this field as a research priority, this review aims to shed light on the importance of investigating sex-genetic interactions for brain disorders, focusing on two areas: (i) variant-trait associations and (ii) transcriptomics (i.e. gene expression, transcript usage and regulation). We specifically discuss recent advances in the field, current gaps and provide considerations for future studies.

Predictive models of eukaryotic transcriptional regulation reveals changes in transcription factor roles and promoter usage between metabolic conditions

Transcription factors (TF) are central to transcriptional regulation, but they are often studied in relative isolation and without close control of the metabolic state of the cell. Here, we describe genome-wide binding (by ChIP-exo) of 15 yeast TFs in four chemostat conditions that cover a range of metabolic states. We integrate this data with transcriptomics and six additional recently mapped TFs to identify predictive models describing how TFs control gene expression in different metabolic conditions. Contributions by TFs to gene regulation are predicted to be mostly activating, additive and well approximated by assuming linear effects from TF binding signal. Notably, using TF binding peaks from peak finding algorithms gave distinctly worse predictions than simply summing the low-noise and high-resolution TF ChIP-exo reads on promoters. Finally, we discover indications of a novel functional role for three TFs; Gcn4, Ert1 and Sut1 during nitrogen limited aerobic fermentation. In only this condition, the three TFs have correlated binding to a large number of genes (enriched for glycolytic and translation processes) and a negative correlation to target gene transcript levels.

Genome-Wide Investigation of Genes Regulated by ERα in Breast Cancer Cells

Estrogen receptor alpha (ERα), which has been detected in over 70% of breast cancer cases, is a driving factor for breast cancer growth. For investigating the underlying genes and networks regulated by ERα in breast cancer, RNA-seq was performed between ERα transgenic MDA-MB-231 cells and wild type MDA-MB-231 cells. A total of 267 differentially expressed genes (DEGs) were identified. Then bioinformatics analyses were performed to illustrate the mechanism of ERα. Besides, by comparison of RNA-seq data obtained from MDA-MB-231 cells and microarray dataset obtained from estrogen (E2) stimulated MCF-7 cells, an overlap of 126 DEGs was screened. The expression level of ERα was negatively associated with metastasis and EMT in breast cancer. We further verified that ERα might inhibit metastasis by regulating of VCL and TNFRSF12A, and suppress EMT by the regulating of JUNB and ID3. And the relationship between ERα and these genes were validated by RT-PCR and correlation analysis based on TCGA database. By PPI network analysis, we identified TOP5 hub genes, FOS, SP1, CDKN1A, CALCR and JUNB, which were involved in cell proliferation and invasion. Taken together, the whole-genome insights carried in this work can help fully understanding biological roles of ERα in breast cancer.

Nuclear receptors connect progenitor transcription factors to cell cycle control

The specification and growth of organs is controlled simultaneously by networks of transcription factors. While the connection between these transcription factors with fate determinants is increasingly clear, how they establish the link with the cell cycle is far less understood. Here we investigate this link in the developing Drosophila eye, where two transcription factors, the MEIS1 homologue hth and the Zn-finger tsh, synergize to stimulate the proliferation of naïve eye progenitors. Experiments combining transcriptomics, open-chromatin profiling, motif analysis and functional assays indicate that these progenitor transcription factors exert a global regulation of the proliferation program. Rather than directly regulating cell cycle genes, they control proliferation through an intermediary layer of nuclear receptors of the ecdysone/estrogen-signaling pathway. This regulatory subnetwork between hth, tsh and nuclear receptors might be conserved from Drosophila to mammals, as we find a significant co-overexpression of their human homologues in specific cancer types.

Systematic identification of Ctr9 regulome in ERα-positive breast cancer

Background

We had previously identified Ctr9, the key scaffold subunit of the human RNA polymerase II (RNAPII) associated factor complex (PAFc), as a key factor regulating a massive ERα target gene expression and ERα-positive breast cancer growth. Furthermore, we have shown that knockdown of Ctr9 reduces ERα protein stability and decreases the occupancy of ERα and RNAPII at a few ERα-target loci. However, it remains to be determined whether Ctr9 controls ERα-target gene expression by regulating the global chromatin occupancy of ERα and RNAPII in the presence of estrogen.

Results

In this study, we determined the genome-wide ERα and RNAPII occupancy in response to estrogen treatment and/or Ctr9 knockdown by performing chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-seq). We found that loss of Ctr9 dramatically decreases the global occupancy of ERα and RNAPII, highlighting the significance of Ctr9 in regulating estrogen signaling in ERα-positive breast cancer cells. Combining this resource with previously published genomic data sets, we identified a unique subset of ERα and Ctr9 target genes, and further delineated the independent function of Ctr9 from other subunits in PAFc when regulating transcription.

Conclusions

Our data demonstrated that Ctr9, independent of other PAFc subunits, controls ERα-target gene expression by regulating global chromatin occupancies of ERα and RNAPII.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3248-3) contains supplementary material, which is available to authorized users.

SERPINA1 is a direct estrogen receptor target gene and a predictor of survival in breast cancer patients

Of all breast cancer patients, about 70% are ER+ and 10% are ER+/HER2+. The ER+/HER2+ patients have a worse outcome compared to ER+/HER2- patients. Currently there is a lack of effective prognosis biomarkers for the prediction of outcome in ER+/HER2+ patients. Genome-wide differences in ER binding between the endocrine-responsive and endocrine-resistant cells were discovered using ChIP-seq, and combined with gene expression microarray data to identify direct ER target genes. These genes were correlated to survival outcome using publicly available breast cancer patient cohorts. We found the expression of the gene SERPINA1 to have a significant predictive value for the overall survival (OS) of ER+ patients in the TCGA cohort, and validated this finding in the Curtis cohort. SERPINA1 also has a significant predictive value for the OS of ER+/HER2+ patients in the TCGA cohort, with validation in the Bild cohort. The expression of SERPINA1 can be suppressed by fulvestrant and HER2 siRNA. Our results indicate that ER is constitutively activated, resulting in an E2-independent ER binding to the SERPINA1 gene and upregulation of SERPINA1 expression. Importantly, results of survival correlation suggests that high expression of SERPINA1 could be predictive for a better clinical outcome of ER+ and ER+/HER2+ patients.

Epigenetics of the failing heart

With the impressive advancement in high-throughput 'omics' technologies over the past two decades, epigenetic mechanisms have emerged as the regulatory interface between the genome and environmental factors. These mechanisms include DNA methylation, histone modifications, ATP-dependent chromatin remodeling and RNA-based mechanisms. Their highly interdependent and coordinated action modulates the chromatin structure controlling access of the transcription machinery and thereby regulating expression of target genes. Given the rather limited proliferative capability of human cardiomyocytes, epigenetic regulation appears to play a particularly important role in the myocardium. The highly dynamic nature of the epigenome allows the heart to adapt to environmental challenges and to respond quickly and properly to cardiac stress. It is now becoming evident that histone-modifying and chromatin-remodeling enzymes as well as numerous non-coding RNAs play critical roles in cardiac development and function, while their dysregulation contributes to the onset and development of pathological cardiac remodeling culminating in HF. This review focuses on up-to-date knowledge about the epigenetic mechanisms and highlights their emerging role in the healthy and failing heart. Uncovering the determinants of epigenetic regulation holds great promise to accelerate the development of successful new diagnostic and therapeutic strategies in human cardiac disease.

Dynamics of Protein Expression Reveals Primary Targets and Secondary Messengers of Estrogen Receptor Alpha Signaling in MCF-7 Breast Cancer Cells

Estrogen receptor alpha (ERα)-mediated proliferation of breast cancer cells is facilitated through expression of multiple primary target genes, products of which induce a secondary response to stimulation. To differentiate between the primary and secondary target proteins of ERα signaling, we measured dynamics of protein expression induced by 17β-estradiol in MCF-7 breast cancer cells. Measurement of the global proteomic effects of estradiol by stable isotope labeling by amino acids in cell culture (SILAC) resulted in identification of 103 estrogen-regulated proteins, with only 40 of the corresponding genes having estrogen response elements. Selected reaction monitoring (SRM) assays were used to validate the differential expression of 19 proteins and measure the dynamics of their expression within 72 h after estradiol stimulation, and in the absence or presence of 4-hydroxytamoxifen, to confirm ERα-mediated signaling. Dynamics of protein expression unambiguously revealed early and delayed response proteins and well correlated with presence or absence of estrogen response elements in the corresponding genes. Finally, we quantified dynamics of protein expression in a rarely studied network of transcription factors with a negative feedback loop (ERα-EGR3-NAB2). Because NAB2 protein is a repressor of EGR3-induced transcription, siRNA-mediated silencing of NAB2 resulted in the enhanced expression of the EGR3-induced protein ITGA2. To conclude, we provided a high-quality proteomic resource to supplement genomic and transcriptomic studies of ERα signaling.

Wbp2 is required for normal glutamatergic synapses in the cochlea and is crucial for hearing

WBP2 encodes the WW domain-binding protein 2 that acts as a transcriptional coactivator for estrogen receptor α (ESR1) and progesterone receptor (PGR). We reported that the loss of Wbp2 expression leads to progressive high-frequency hearing loss in mouse, as well as in two deaf children, each carrying two different variants in the WBP2 gene. The earliest abnormality we detect in Wbp2-deficient mice is a primary defect at inner hair cell afferent synapses. This study defines a new gene involved in the molecular pathway linking hearing impairment to hormonal signalling and provides new therapeutic targets.

miRNAs regulated by estrogens, tamoxifen, and endocrine disruptors and their downstream gene targets

MicroRNAs (miRNAs) are short (22 nucleotides), single-stranded, non-coding RNAs that form complimentary base-pairs with the 3' untranslated region of target mRNAs within the RNA-induced silencing complex (RISC) and block translation and/or stimulate mRNA transcript degradation. The non-coding miRBase (release 21, June 2014) reports that human genome contains ∼ 2588 mature miRNAs which regulate ∼ 60% of human protein-coding mRNAs. Dysregulation of miRNA expression has been implicated in estrogen-related diseases including breast cancer and endometrial cancer. The mechanism for estrogen regulation of miRNA expression and the role of estrogen-regulated miRNAs in normal homeostasis, reproduction, lactation, and in cancer is an area of great research and clinical interest. Estrogens regulate miRNA transcription through estrogen receptors α and β in a tissue-specific and cell-dependent manner. This review focuses primarily on the regulation of miRNA expression by ligand-activated ERs and their bona fide gene targets and includes miRNA regulation by tamoxifen and endocrine disrupting chemicals (EDCs) in breast cancer and cell lines.

Enhancer-bound LDB1 regulates a corticotrope promoter-pausing repression program

Substantial evidence supports the hypothesis that enhancers are critical regulators of cell-type determination, orchestrating both positive and negative transcriptional programs; however, the basic mechanisms by which enhancers orchestrate interactions with cognate promoters during activation and repression events remain incompletely understood. Here we report the required actions of LIM domain-binding protein 1 (LDB1)/cofactor of LIM homeodomain protein 2/nuclear LIM interactor, interacting with the enhancer-binding protein achaete-scute complex homolog 1, to mediate looping to target gene promoters and target gene regulation in corticotrope cells. LDB1-mediated enhancer:promoter looping appears to be required for both activation and repression of these target genes. Although LDB1-dependent activated genes are regulated at the level of transcriptional initiation, the LDB1-dependent repressed transcription units appear to be regulated primarily at the level of promoter pausing, with LDB1 regulating recruitment of metastasis-associated 1 family, member 2, a component of the nucleosome remodeling deacetylase complex, on these negative enhancers, required for the repressive enhancer function. These results indicate that LDB1-dependent looping events can deliver repressive cargo to cognate promoters to mediate promoter pausing events in a pituitary cell type.

SOX2, a predictor of survival in gastric cancer, inhibits cell proliferation and metastasis by regulating PTEN

Inconsistent results of SOX2 expression have been reported in gastric cancer (GC). Here, we demonstrated that SOX2 was progressively downregulated during GC development via immunochemistry in 755 human gastric specimens. Low SOX2 levels were associated with pathological stage and clinical outcome. Multivariate analysis indicated that SOX2 protein expression served as an independent prognostic marker for GC. Gain-and loss-of function studies showed the anti-proliferative, anti-metastatic, and pro-apoptotic effects of SOX2 in GC. PTEN was selected as SOX2 targets by cDNA microarray and ChIP-DSL, further identified by luciferase assays, EMSA and ChIP-PCR. PTEN upregulation in response to SOX2-enforced expression suppressed GC malignancy via regulating Akt dephosphorylation. PTEN inhibition reversed SOX2-induced anticancer effects. Moreover, concordant positivity of SOX2 and PTEN proteins in nontumorous tissues but lost in matched GC specimens predicted a worse patient prognosis. Thus, SOX2 proved to be a new marker for evaluating GC outcome.

Enriched trimethylation of lysine 4 of histone H3 of WDR63 enhanced osteogenic differentiation potentials of stem cells from apical papilla

INTRODUCTION

Dental tissue-derived mesenchymal stem cells (MSCs) are a reliable cell source for dental tissue regeneration. However, the molecular mechanisms underlying their directed differentiation remain unclear, thus limiting their use. Trimethylation of lysine 4 of histone H3 (H3K4Me3) correlates with gene activation and osteogenic differentiation. We used stem cells from apical papilla (SCAPs) to investigate the effects of genomic changes in H3K4Me3 modification at gene promoter regions on MSC osteogenic differentiation.

METHODS

ChIP-on-chip assays were applied to compare the H3K4Me3 profiles at gene promoter regions of undifferentiated and differentiated SCAPs. Alkaline phosphatase activity assay, alizarin red staining, quantitative analysis of calcium, the expressions of osteogenesis-related genes, and transplantation in nude mice were used to investigate the osteogenic differentiation potentials of SCAPs.

RESULTS

In differentiated SCAPs, 119 gene promoters exhibited >2-fold increases of H3K4Me3; in contrast, the promoter regions of 21 genes exhibited >2-fold decreases of H3K4Me3. On the basis of enriched H3K4Me3 and up-regulated gene expression on the osteogenic differentiation of SCAPs, WDR63 may be a potential regulator for mediating SCAP osteogenic differentiation. Through gain-of-function and loss-of-function studies, we discovered that WDR63 enhances alkaline phosphatase activity, mineralization, and the expression of BSP, OSX, and RUNX2 in vitro. In addition, transplant experiments in nude mice confirmed that SCAP osteogenesis is triggered by activated WDR63.

CONCLUSIONS

These results indicate that WDR63 is a positive enhancer for SCAP osteogenic differentiation and suggest that activation of WDR63 signaling might improve tissue regeneration mediated by MSCs of dental origin.

Biased, non-equivalent gene-proximal and -distal binding motifs of orphan nuclear receptor TR4 in primary human erythroid cells

We previously reported that TR2 and TR4 orphan nuclear receptors bind to direct repeat (DR) elements in the ε- and γ-globin promoters, and act as molecular anchors for the recruitment of epigenetic corepressors of the multifaceted DRED complex, thereby leading to ε- and γ-globin transcriptional repression during definitive erythropoiesis. Other than the ε- and γ-globin and the GATA1 genes, TR4-regulated target genes in human erythroid cells remain unknown. Here, we identified TR4 binding sites genome-wide using chromatin immunoprecipitation followed by massively parallel sequencing (ChIP-seq) as human primary CD34⁺ hematopoietic progenitors differentiated progressively to late erythroid precursors. We also performed whole transcriptome analyses by RNA-seq to identify TR4 downstream targets after lentiviral-mediated TR4 shRNA knockdown in erythroid cells. Analyses from combined ChIP-seq and RNA-seq datasets indicate that DR1 motifs are more prevalent in the proximal promoters of TR4 direct target genes, which are involved in basic biological functions (e.g., mRNA processing, ribosomal assembly, RNA splicing and primary metabolic processes). In contrast, other non-DR1 repeat motifs (DR4, ER6 and IR1) are more prevalent at gene-distal TR4 binding sites. Of these, approximately 50% are also marked with epigenetic chromatin signatures (such as P300, H3K27ac, H3K4me1 and H3K27me3) associated with enhancer function. Thus, we hypothesize that TR4 regulates gene transcription via gene-proximal DR1 sites as TR4/TR2 heterodimers, while it can associate with novel nuclear receptor partners (such as RXR) to bind to distant non-DR1 consensus sites. In summary, this study reveals that the TR4 regulatory network is far more complex than previously appreciated and that TR4 regulates basic, essential biological processes during the terminal differentiation of human erythroid cells.

Sequential genome-wide binding studies investigated by deep sequencing (ChIP-seq) represent a powerful tool for investigating the temporal sequence of gene activation and repression events that take place as cells differentiate. Here, we report the binding of an “orphan” nuclear receptor (one for which no ligand has been identified) to its cognate genomic regulatory sites and perform the functional analysis to validate its downstream targets as precursor cells differentiate from very early human hematopoietic progenitors into red blood cells. We discovered that when this receptor is bound at gene proximal promoters, it recognizes a different DNA sequence than when it binds to more distant regulatory sites (enhancers and silencers). Since this receptor can either activate or repress specific target genes, the data suggest the intriguing possibility that the two different modes of DNA recognition may reflect association of the receptor with different partner molecules when regulating gene expression from proximal or distal sequences.

An integrated approach for the identification of USF1-centered transcriptional regulatory networks during liver regeneration

Liver regeneration after partial hepatectomy (PH) is a synchronized process that is precisely controlled by system-wide transcriptional regulatory networks. To clarify the transcriptional changes and regulatory networks that involve transcription factors (TFs) and their target genes during the priming phase, an advanced mouse oligonucleotide array-based transcription factor assay (MOUSE OATFA), mRNA microarray analysis, bioinformatic analysis and ChIP-on-chip experiments were used. A total of 774 genes were upregulated or downregulated in PH liver samples compared with the sham operation (SH) group. Seventeen TFs showed significant changes in activity in the regenerating livers, some of which have not been extensively studied in previous reports, including upstream stimulatory transcription factor 1 (USF1). The TF signatures from MOUSE OATFA were combined with mRNA expression profiles and ChIP-on-chip analyses to construct experimental transcriptional regulatory networks in regenerating livers. USF1-centered regulatory networks were further confirmed by ChIP assays, revealing some of its target genes and novel coregulatory networks. The combination of MOUSE OATFA with transcriptome profiling and bioinformatic analysis represents a novel paradigm for the comprehensive prediction of transcriptional coregulatory networks during the early phase of liver regeneration.

Genome-wide activity of unliganded estrogen receptor-α in breast cancer cells

Estrogen receptor-α (ERα) has central role in hormone-dependent breast cancer and its ligand-induced functions have been extensively characterized. However, evidence exists that ERα has functions that are independent of ligands. In the present work, we investigated the binding of ERα to chromatin in the absence of ligands and its functions on gene regulation. We demonstrated that in MCF7 breast cancer cells unliganded ERα binds to more than 4,000 chromatin sites. Unexpectedly, although almost entirely comprised in the larger group of estrogen-induced binding sites, we found that unliganded-ERα binding is specifically linked to genes with developmental functions, compared with estrogen-induced binding. Moreover, we found that siRNA-mediated down-regulation of ERα in absence of estrogen is accompanied by changes in the expression levels of hundreds of coding and noncoding RNAs. Down-regulated mRNAs showed enrichment in genes related to epithelial cell growth and development. Stable ERα down-regulation using shRNA, which caused cell growth arrest, was accompanied by increased H3K27me3 at ERα binding sites. Finally, we found that FOXA1 and AP2γ binding to several sites is decreased upon ERα silencing, suggesting that unliganded ERα participates, together with other factors, in the maintenance of the luminal-specific cistrome in breast cancer cells.

Inhibition of KIF22 suppresses cancer cell proliferation by delaying mitotic exit through upregulating CDC25C expression

KIF22 is a microtubule-dependent molecular motor protein with DNA-binding capacity. It is well known that KIF22 plays a critical role in cell mitosis as a motor protein; however, the role of altered KIF22 expression and its transcriptional regulatory function in cancer development have not yet been defined. This study showed that KIF22 was overexpressed in human cancer tissues, and inhibition of KIF22 significantly led to accumulation of cells in the G2/M phases, resulting in suppression of cancer cell proliferation. The investigation of the molecular mechanisms demonstrated that cell division cycle 25C (CDC25C) is a direct transcriptional target of KIF22, and inhibition of KIF22 increased CDC25C expression and cyclin-dependent kinase 1 (CDK1) activity, resulting in delayed mitotic exit. Phosphorylation of KIF22 was required for its transcriptional regulatory function and the reduction of CDK1 activity. Thus, we conclude that inhibition of KIF22 suppresses cancer cell proliferation by delaying mitotic exit through the transcriptional upregulation of CDC25C.

Estrogen receptor-mediated long-range chromatin interactions and transcription in breast cancer

Estrogen induces the binding of ERα to thousands of locations in the breast cancer genome, preferring intergenic and distal regions rather than near the promoters of estrogen-regulated genes. With recent technological innovations in mapping and characterization of global chromatin organization, evidence now indicates ERα mediates long-range chromatin interactions to control gene transcription. The principles that govern how ERα communicates with their putative target genes via chromosomal interactions are also beginning to unravel. Herein, we summarize our current knowledge on the functional significance of chromatin looping in estrogen-mediated transcription. ERα collaborative factors and other players that contribute to define the genomic interactions in breast cancer cells will also be discussed. Defects in chromatin organization are emerging key players in diseases such as cancer, thus understanding how ERα-mediated chromatin looping affects genome organization will clarify the receptor's role in estrogen responsive pathways sensitive to defects in chromatin organization.

Decoding a signature-based model of transcription cofactor recruitment dictated by cardinal cis-regulatory elements in proximal promoter regions

Genome-wide maps of DNase I hypersensitive sites (DHSs) reveal that most human promoters contain perpetually active cis-regulatory elements between −150 bp and +50 bp (−150/+50 bp) relative to the transcription start site (TSS). Transcription factors (TFs) recruit cofactors (chromatin remodelers, histone/protein-modifying enzymes, and scaffold proteins) to these elements in order to organize the local chromatin structure and coordinate the balance of post-translational modifications nearby, contributing to the overall regulation of transcription. However, the rules of TF-mediated cofactor recruitment to the −150/+50 bp promoter regions remain poorly understood. Here, we provide evidence for a general model in which a series of cis-regulatory elements (here termed ‘cardinal’ motifs) prefer acting individually, rather than in fixed combinations, within the −150/+50 bp regions to recruit TFs that dictate cofactor signatures distinctive of specific promoter subsets. Subsequently, human promoters can be subclassified based on the presence of cardinal elements and their associated cofactor signatures. In this study, furthermore, we have focused on promoters containing the nuclear respiratory factor 1 (NRF1) motif as the cardinal cis-regulatory element and have identified the pervasive association of NRF1 with the cofactor lysine-specific demethylase 1 (LSD1/KDM1A). This signature might be distinctive of promoters regulating nuclear-encoded mitochondrial and other particular genes in at least some cells. Together, we propose that decoding a signature-based, expanded model of control at proximal promoter regions should lead to a better understanding of coordinated regulation of gene transcription.

Human cells exploit different mechanisms to coordinate the expression of both protein-coding and non-coding RNAs. Elucidating these mechanisms is essential to understanding normal physiology and disease. In our attempt to identify new regulatory layers acting particularly at proximal promoters, we have computationally analyzed the genomic sequences located from −150 bp to +50 bp relative to the transcriptional start site (TSS), which are often at the center of ‘open’ chromatin regions in human promoters. We have confirmed the presence of a series of cis-regulatory elements (here referred to as ‘cardinal’ motifs) that show a strong preference for these short regions. Interestingly, these elements tend to act independently rather than in fixed combinations. Therefore, we propose that they confer unique regulatory features to the human promoter subsets that contain each of these particular elements. In agreement with this model, we have identified a large repertoire of preferential partnerships between transcription factors recognizing cardinal motifs and their associated proteins (cofactors), thus decoding a signature-based model that distinguishes distinctive regulatory types of promoters based on cardinal motifs. These signatures may underlie a new layer of transcriptional regulation to orchestrate coordinated gene expression in human promoters.

Proteomic analysis of coregulators bound to ERα on DNA and nucleosomes reveals coregulator dynamics

Chromatin immunoprecipitation studies have mapped protein occupancies at many genomic loci. However, a detailed picture of the complexity of coregulators (CoRs) bound to a defined enhancer along with a transcription factor is missing. To address this, we used biotin-DNA pull-down assays coupled with mass spectrometry-immunoblotting to identify at least 17 CoRs from nuclear extracts bound to 17β-estradiol (E2)-liganded estrogen receptor-α on estrogen response elements (EREs). Unexpectedly, these complexes initially are biochemically stable and contain certain atypical corepressors. Addition of ATP dynamically converts these complexes to an "activated" state by phosphorylation events, primarily mediated by DNA-dependent protein kinase. Importantly, a "natural" ERE-containing enhancer and nucleosomal EREs recruit similar complexes. We further discovered the mechanism whereby H3K4me3 stimulates ERα-mediated transcription as compared with unmodified nucleosomes. H3K4me3 templates promote specific CoR dynamics in the presence of ATP and AcCoA, as manifested by CBP/p300 and SRC-3 dismissal and SAGA and TFIID stabilization/recruitment.

Interaction of glucocorticoid receptor (GR) with estrogen receptor (ER) α and activator protein 1 (AP1) in dexamethasone-mediated interference of ERα activity

The role of glucocorticoids in the inhibition of estrogen (17-β-estradiol (E2))-regulated estrogen receptor (ER)-positive breast cancer cell proliferation is well established. We and others have seen that synthetic glucocorticoid dexamethasone (Dex) antagonizes E2-stimulated endogenous ERα target gene expression. However, how glucocorticoids negatively regulate the ERα signaling pathway is still poorly understood. ChIP studies using ERα- and glucocorticoid receptor (GR)-positive MCF-7 cells revealed that GR occupies several ERα-binding regions (EBRs) in cells treated with E2 and Dex simultaneously. Interestingly, there was little or no GR loading to these regions when cells were treated with E2 or Dex alone. The E2+Dex-dependent GR recruitment is associated with the displacement of ERα and steroid receptor coactivator-3 from the target EBRs leading to the repression of ERα-mediated transcriptional activation. The recruitment of GR to EBRs requires assistance from ERα and FOXA1 and is facilitated by AP1 binding within the EBRs. The GR binding to EBRs is mediated via direct protein-protein interaction between the GR DNA-binding domain and ERα. Limited mutational analyses indicate that arginine 488 located within the C-terminal zinc finger domain of the GR DNA-binding domain plays a critical role in stabilizing this interaction. Together, the results of this study unravel a novel mechanism involved in glucocorticoid inhibition of ERα transcriptional activity and E2-mediated cell proliferation and thus establish a foundation for future exploitation of the GR signaling pathway in the treatment of ER-positive breast cancer.

Functional roles of enhancer RNAs for oestrogen-dependent transcriptional activation

The functional importance of gene enhancers in regulated gene expression is well established. In addition to widespread transcription of long non-coding RNAs (lncRNAs) in mammalian cells, bidirectional ncRNAs are transcribed on enhancers, and are thus referred to as enhancer RNAs (eRNAs). However, it has remained unclear whether these eRNAs are functional or merely a reflection of enhancer activation. Here we report that in human breast cancer cells 17β-oestradiol (E2)-bound oestrogen receptor α (ER-α) causes a global increase in eRNA transcription on enhancers adjacent to E2-upregulated coding genes. These induced eRNAs, as functional transcripts, seem to exert important roles for the observed ligand-dependent induction of target coding genes, increasing the strength of specific enhancer-promoter looping initiated by ER-α binding. Cohesin, present on many ER-α-regulated enhancers even before ligand treatment, apparently contributes to E2-dependent gene activation, at least in part by stabilizing E2/ER-α/eRNA-induced enhancer-promoter looping. Our data indicate that eRNAs are likely to have important functions in many regulated programs of gene transcription.

β-Estradiol-dependent activation of the JAK/STAT pathway requires p/CIP and CARM1

The steroid receptor coactivator p/CIP, also known as SRC-3, is an oncogene commonly amplified in breast and ovarian cancers. p/CIP is known to associate with coactivator arginine methyltransferase 1 (CARM1) on select estrogen responsive genes. We have shown, using a ChIP-on-chip approach, that in response to stimulation with 17β-estradiol (E2), the p/CIP/CARM1 complex is recruited to 204 proximal promoters in MCF-7 cells. Many of the complex target genes have been previously implicated in signaling pathways related to oncogenesis. Jak2, a member of the Jak/Stat signaling cascade, is one of the direct E2-dependent targets of the p/CIP/CARM1 complex. Following E2-treatment, histone modifications at the Jak2 promoter are reflective of a transcriptionally permissive gene, and modest changes in RNA and protein expression lead us to suggest that an additional factor(s) may be required for a more notable transcriptional and functional response. Bioinformatic examination of the 204 proximal promoter sequences of p/CIP/CARM1 targets supports the idea that transcription factor crosstalk is likely the favored mechanism of E2-dependent p/CIP/CARM1 complex recruitment. This data may have implications towards understanding the oncogenic role of p/CIP in breast cancer and ultimately allow for the identification of new prognostic indicators and/or viable therapeutic targets.

Chromatin immunoprecipitation: advancing analysis of nuclear hormone signaling

Recent decades have been filled with groundbreaking research in the field of endocrine hormone signaling. Pivotal events like the isolation and purification of the estrogen receptor, the cloning of glucocorticoid receptor cDNA, or dissemination of nuclear hormone receptor (NHR) DNA binding sequences are well recognized for their contributions. However, the novel genome-wide and gene-specific information obtained over the last decade describing NHR association with chromatin, cofactors, and epigenetic modifications, as well as their role in gene regulation, has been largely facilitated by the adaptation of the chromatin immunoprecipitation (ChIP) technique. Use of ChIP-based technologies has taken the field of hormone signaling from speculating about the transcription-enabling properties of acetylated chromatin and putative transcription (co-)factor genomic occupancy to demonstrating the detailed, stepwise mechanisms of factor binding and transcriptional initiation; from treating hormone-induced transcription as a steady-state event to understanding its dynamic and cyclic nature; from looking at the DNA sequences recognized by various DNA-binding domains in vitro to analyzing the cell-specific genome-wide pattern of nuclear receptor binding and interpreting its physiological implications. Not only have these events propelled hormone research, but, as some of the pioneering studies, have also contributed tremendously to the field of molecular endocrinology as a whole. In this review, we give a brief summary of some of the most important discoveries in hormone signaling using ChIP and other derivative techniques and speculate on what the future may hold.

Identification of direct forkhead box O1 targets involved in palmitate-induced apoptosis in clonal insulin-secreting cells using chromatin immunoprecipitation coupled to DNA selection and ligation

AIMS/HYPOTHESIS

The transcription factor, forkhead box (FOX)O1, is involved in fatty acid-induced apoptosis in pancreatic beta cells, but the precise mechanism is poorly understood. We aimed to identify which direct downstream targets of FOXO1 are involved in palmitate-induced apoptosis in the pancreatic beta cell line MIN6.

METHODS

Chromatin immunoprecipitation (ChIP) coupled to a DNA selection and ligation technique (ChIP-DSL) was used to identify the direct targets of FOXO1. The mRNA level was examined by real-time PCR assay. The ChIP-DSL results were verified using ChIP-PCR and luciferase assay, respectively. The cell apoptosis rate was determined by TUNEL assay and by scoring cells with pycnotic nuclei.

RESULTS

We identified 189 target genes and selected 106 targets for expression analysis in MIN6 cells treated with palmitate. The results showed that six genes were significantly upregulated and four were downregulated. Binding of FOXO1 to the promoters was determined by ChIP-PCR and confirmed by luciferase assay. Among the ten up- and downregulated genes, mRNA expression of A930038C07Rik was significantly decreased and that of Ppa1 was increased in 8-week-old db/db mice. The apoptosis assay showed that overproduction of the protein 'RIKEN cDNA A930038C07' (A930038C07Rik) drastically enhanced palmitate-induced apoptosis, while pyrophosphatase (inorganic) 1 (PPA1) partially protected the cells from apoptosis. Knockdown of PPA1, moreover, significantly increased apoptosis.

CONCLUSIONS/INTERPRETATION

We identified for the first time FOXO1 targets in MIN6 cells treated with palmitate, thus revealing the important roles of A930038C07Rik and PPA1 in palmitate-induced cell apoptosis. These results shed light on the mechanisms of palmitate-induced apoptosis in pancreatic beta cells.

[DNA-binding profiles of mammalian transcription factors]

The differential gene expression is the molecular base of development and responses to stimuli of organisms. Transcription factors (TFs) play important regulatory roles in this kind of differential gene expression. Therefore, to elucidate how these TFs regulate the complex differential gene expression, it is necessary to identify all target genes of them and construct the gene transcription regulatory network controlled by them. DNA binding is a key step for TFs regulating gene transcription. Therefore, in order to identify their target genes, it is indispensable to identify all possible DNA sequences that can be recognized and bound by TFs at the molecular level of their interactions with DNA, i.e., construction of the DNA-binding profiles of TFs. In recent years, along with the development of DNA microarray and high-throughput DNA sequencing techniques, there appeared some revolutionary new techniques for constructing DNA-binding profiles of TFs, which greatly promotes studies in this field. These techniques include ChIP-chip and ChIP-Seq for constructing in vivo DNA-binding profiles of TFs, dsDNA microarray, SELEX-SAGE, Bind-n-Seq, MMP-SELEX, EMSA-Seq, and HiTS-FLIP for constructing in vitro DNA-binding profiles of TFs. This paper reviewed these techniques.

Dual transcriptional activator and repressor roles of TBX20 regulate adult cardiac structure and function

The ongoing requirement in adult heart for transcription factors with key roles in cardiac development is not well understood. We recently demonstrated that TBX20, a transcriptional regulator required for cardiac development, has key roles in the maintenance of functional and structural phenotypes in adult mouse heart. Conditional ablation of Tbx20 in adult cardiomyocytes leads to a rapid onset and progression of heart failure, with prominent conduction and contractility phenotypes that lead to death. Here we describe a more comprehensive molecular characterization of the functions of TBX20 in adult mouse heart. Coupling genome-wide chromatin immunoprecipitation and transcriptome analyses (RNA-Seq), we identified a subset of genes that change expression in Tbx20 adult cardiomyocyte-specific knockout hearts which are direct downstream targets of TBX20. This analysis revealed a dual role for TBX20 as both a transcriptional activator and a repressor, and that each of these functions regulates genes with very specialized and distinct molecular roles. We also show how TBX20 binds to its targets genome-wide in a context-dependent manner, using various cohorts of co-factors to either promote or repress distinct genetic programs within adult heart. Our integrative approach has uncovered several novel aspects of TBX20 and T-box protein function within adult heart. Sequencing data accession number (http://www.ncbi.nlm.nih.gov/geo): GSE30943.

A modulated empirical Bayes model for identifying topological and temporal estrogen receptor α regulatory networks in breast cancer

BACKGROUND

Estrogens regulate diverse physiological processes in various tissues through genomic and non-genomic mechanisms that result in activation or repression of gene expression. Transcription regulation upon estrogen stimulation is a critical biological process underlying the onset and progress of the majority of breast cancer. Dynamic gene expression changes have been shown to characterize the breast cancer cell response to estrogens, the every molecular mechanism of which is still not well understood.

RESULTS

We developed a modulated empirical Bayes model, and constructed a novel topological and temporal transcription factor (TF) regulatory network in MCF7 breast cancer cell line upon stimulation by 17β-estradiol stimulation. In the network, significant TF genomic hubs were identified including ER-alpha and AP-1; significant non-genomic hubs include ZFP161, TFDP1, NRF1, TFAP2A, EGR1, E2F1, and PITX2. Although the early and late networks were distinct (<5% overlap of ERα target genes between the 4 and 24 h time points), all nine hubs were significantly represented in both networks. In MCF7 cells with acquired resistance to tamoxifen, the ERα regulatory network was unresponsive to 17β-estradiol stimulation. The significant loss of hormone responsiveness was associated with marked epigenomic changes, including hyper- or hypo-methylation of promoter CpG islands and repressive histone methylations.

CONCLUSIONS

We identified a number of estrogen regulated target genes and established estrogen-regulated network that distinguishes the genomic and non-genomic actions of estrogen receptor. Many gene targets of this network were not active anymore in anti-estrogen resistant cell lines, possibly because their DNA methylation and histone acetylation patterns have changed.

Pre-mRNA splicing during transcription in the mammalian system

Splicing of RNA polymerase II transcripts is a crucial step in gene expression and a key generator of mRNA diversity. Splicing and transcription have generally been studied in isolation, although in vivo pre-mRNA splicing occurs in concert with transcription. The two processes appear to be functionally connected because a number of variables that regulate transcription have been identified as also influencing splicing. However, the mechanisms that couple the two processes are largely unknown. This review highlights the observations that implicate splicing as occurring during transcription and describes the evidence supporting functional interactions between the two processes. I discuss postulated models of how splicing couples to transcription and consider the potential impact that such coupling might have on exon recognition. WIREs RNA 2011 2 700-717 DOI: 10.1002/wrna.86 For further resources related to this article, please visit the WIREs website.

ER-/ER+ breast cancer cell lines exhibited different resistance to paclitaxel through pulse selection

In this study, we established two PTX-resistant breast cancer cell lines, 231 TIM10 and MCF-7 TIM10, from ER-negative MDA-MB-231 cells and ER-positive MCF-7 cells by pulse selection, respectively. We found that 231 TIM10 variants acquired higher drug resistance than MCF-7 TIM10 variants by the pulse selection, although ER-positive MCF-7 cells were not as sensitive as ER-negative MDA-MB-231 to the initial pulses with PTX. 231 TIM10 had 11.9-fold greater resistance (RI = 11.9) than the parental MDA-MB-231 cells, while MCF-7 TIM10 got 5.5-fold resistance (RI = 5.5) when compared with the parental MCF-7 cells. In the presence of 5nM PTX, 231 TIM10 cells formed colonies, but no colony formed when MCF-7 TIM10 cells were cultured in the same condition. These data have two implications. First, the ER expression state might be an important determinant for the response of breast cancer cells to paclitaxel treatment. Second, ER-negative and ER-positive breast cancer cells develop drug-resistance phenotype with distinctive mechanisms. Our work not only established useful models for studying the paclitaxel resistance but also provides interesting clues to understand the mechanisms underlying the drug resistance of ER-negative and ER-positive breast cancer cells.

9p21 DNA variants associated with coronary artery disease impair interferon-γ signalling response

Genome wide association studies (GWAS) have identified SNPs in the 9p21 gene desert associated with coronary artery disease (CAD)1–4 and Type 2 diabetes (T2D)5–7. Despite evidence for a role of the associated interval in neighboring gene regulation8–10, the biological underpinnings of these genetic associations to CAD or T2D have not yet been explained. Here we identify 33 enhancers in 9p21; the interval is the second densest gene-desert for predicted enhancers and 6 times denser than the whole genome (p<6.55 10⁻³³). The CAD risk alleles of SNPs rs10811656/rs10757278 are located in one of these enhancers and disrupt a binding site for STAT1. Lymphoblastoid cell lines (LCL) homozygous for the CAD risk haplotype exhibit no binding of STAT1, and in LCL homozygous for the CAD non-risk haplotype binding of STAT1 inhibits CDKN2BAS expression, which is reversed by siRNA knock-down of STAT1. Using a new, open-ended approach to detect long-distance interactions (3D-DSL), we find that in human vascular endothelium cells (HUVEC) the enhancer interval containing the CAD locus physically interacts with the CDKN2A/B locus, the MTAP gene and an interval downstream of INFA21. In HUVEC, IFNγ activation strongly affects the structure of the chromatin and the transcriptional regulation in the 9p21 locus, including STAT1 binding, long-range enhancer interactions and altered expression of neighboring genes. Our findings establish a link between CAD genetic susceptibility and the response to inflammatory signaling in a vascular cell type and thus demonstrate the utility of GWAS findings to direct studies to novel genomic loci and biological processes important for disease etiology.

The cardiac transcription network modulated by Gata4, Mef2a, Nkx2.5, Srf, histone modifications, and microRNAs

The transcriptome, as the pool of all transcribed elements in a given cell, is regulated by the interaction between different molecular levels, involving epigenetic, transcriptional, and post-transcriptional mechanisms. However, many previous studies investigated each of these levels individually, and little is known about their interdependency. We present a systems biology study integrating mRNA profiles with DNA–binding events of key cardiac transcription factors (Gata4, Mef2a, Nkx2.5, and Srf), activating histone modifications (H3ac, H4ac, H3K4me2, and H3K4me3), and microRNA profiles obtained in wild-type and RNAi–mediated knockdown. Finally, we confirmed conclusions primarily obtained in cardiomyocyte cell culture in a time-course of cardiac maturation in mouse around birth. We provide insights into the combinatorial regulation by cardiac transcription factors and show that they can partially compensate each other's function. Genes regulated by multiple transcription factors are less likely differentially expressed in RNAi knockdown of one respective factor. In addition to the analysis of the individual transcription factors, we found that histone 3 acetylation correlates with Srf- and Gata4-dependent gene expression and is complementarily reduced in cardiac Srf knockdown. Further, we found that altered microRNA expression in Srf knockdown potentially explains up to 45% of indirect mRNA targets. Considering all three levels of regulation, we present an Srf-centered transcription network providing on a single-gene level insights into the regulatory circuits establishing respective mRNA profiles. In summary, we show the combinatorial contribution of four DNA–binding transcription factors in regulating the cardiac transcriptome and provide evidence that histone modifications and microRNAs modulate their functional consequence. This opens a new perspective to understand heart development and the complexity cardiovascular disorders.

An evolutionary conserved orchestra of transcription factors controls cardiac development and function. More recently the contributions of epigenetic and post-transcriptional mechanisms like histone modifications and microRNAs have been identified. The interplay between these regulatory mechanisms is still an open question. However, perturbations of the cardiac transcriptome, triggered by all three levels of regulation, are underlying cardiovascular disease such as congenital heart malformations. Here, we show the impact of the interdependencies of four key transcription factors (Gata4, Mef2a, Nkx2.5, and Srf) and the contribution of activating histone modifications and microRNAs on the cardiac transcriptome. We found that even these non-paralogous transcription factors can partially compensate each other's function. Our data show that histone 3 acetylation correlates with Srf- and Gata4- dependent gene activation. Moreover, we predict a large proportion of indirect Srf targets to be regulated by Srf-dependent microRNAs, which thus might represent an important intermediate layer of regulation. Taken together, we suggest that the different levels regulating cardiac mRNA profiles have a high degree of interdependency and the potential to buffer each other, which presents a starting point to understand the phenotypic variability typically seen in complex cardiovascular disorders.

Epigenetic regulation by nuclear receptors

Nuclear receptors (NRs) represent a vital class of ligand-activated transcription factors responsible for coordinately regulating the expression of genes involved in numerous biological processes. Transcriptional regulation by NRs is conducted through interactions with multiple coactivator or corepressor complexes that modify the chromatin environment to facilitate or inhibit RNA polymerase II binding and transcription initiation. In recent years, studies have identified specific biological roles for cofactors mediating NR signaling through epigenetic modifications such as acetylation and methylation of histones. Intriguingly, genome-wide analysis of NR and cofactor localization has both confirmed findings from single-gene studies and revealed new insights into the relationships between NRs, cofactors and target genes in determining gene expression. Here, we review recent developments in the understanding of epigenetic regulation by NRs across the genome within the context of the well-established background of cofactor complexes and their roles in histone modification.

Aberrant estrogen regulation of PEMT results in choline deficiency-associated liver dysfunction

When dietary choline is restricted, most men and postmenopausal women develop multiorgan dysfunction marked by hepatic steatosis (choline deficiency syndrome (CDS)). However, a significant subset of premenopausal women is protected from CDS. Because hepatic PEMT (phosphatidylethanolamine N-methyltransferase) catalyzes de novo biosynthesis of choline and this gene is under estrogenic control, we hypothesized that there are SNPs in PEMT that disrupt the hormonal regulation of PEMT and thereby put women at risk for CDS. In this study, we performed transcript-specific gene expression analysis, which revealed that estrogen regulates PEMT in an isoform-specific fashion. Locus-wide SNP analysis identified a risk-associated haplotype that was selectively associated with loss of hormonal activation. Chromatin immunoprecipitation, analyzed by locus-wide microarray studies, comprehensively identified regions of estrogen receptor binding in PEMT. The polymorphism (rs12325817) most highly linked with the development of CDS (p < 0.00006) was located within 1 kb of the critical estrogen response element. The risk allele failed to bind either the estrogen receptor or the pioneer factor FOXA1. These data demonstrate that allele-specific ablation of estrogen receptor-DNA interaction in the PEMT locus prevents hormone-inducible PEMT expression, conferring risk of CDS in women.

SLUG-induced elevation of D1 cyclin in breast cancer cells through the inhibition of its ubiquitination

UbcH5c, a member of the UbcH5 family of protein ubiquitin conjugase E2 enzymes, is a critical component of biological processes in human cells, being the initial ubiquitinating enzyme of substrates like IκB, TP53, and cyclin D1. We report here that the metastasis regulator protein SLUG inhibits the expression of UbcH5c directly through chromatin remodeling and thus, among other downstream effects, elevates the level of cyclin D1, thus enhancing the growth rates of breast cancer cells. Overexpression of SLUG in the SLUG-deficient breast cancer cells significantly decreased the levels of mRNA and protein of UbcH5c but only elevated the protein levels of cyclin D1. On the contrary, knockdown of SLUG in SLUG-high breast cancer cells elevated the levels of UbcH5c while decreasing the level of cyclin D1 protein. SLUG is recruited at the E2-box sequence at the UbcH5c gene promoter along with the corepressor CtBP1 and the effector HDAC1 to silence the expression of this gene. Knockdown of UbcH5c in the SLUG-deficient human breast cells elevated the level of cyclin D1 as well as the rates of proliferation and invasiveness of these cells. Whereas the growth rates of the cells are enhanced due to overexpression of SLUG or knockdown of UbcH5c in the breast cancer cells tested, ER(+) cells also acquire resistance to the anti-estrogen 4-hydroxytamoxifen due to the rise of cyclin D1 levels in these cells. This study thus implicates high levels of SLUG and low levels of UbcH5c as a determinant in the progression of metastatic breast cancer.

The current state of chromatin immunoprecipitation

The biological significance of interactions of nuclear proteins with DNA in the context of gene expression, cell differentiation, or disease has immensely been enhanced by the advent of chromatin immunoprecipitation (ChIP). ChIP is a technique whereby a protein of interest is selectively immunoprecipitated from a chromatin preparation to determine the DNA sequences associated with it. ChIP has been widely used to map the localization of post-translationally modified histones, histone variants, transcription factors, or chromatin modifying enzymes on the genome or on a given locus. In spite of its power, ChIP has for a long time remained a cumbersome procedure requiring large numbers of cells. These limitations have sparked the development of modifications to shorten the procedure, simplify sample handling and make ChIP amenable to small numbers of cells. In addition, the combination of ChIP with DNA microarray and high-throughput sequencing technologies has in recent years enabled the profiling of histone modification, histone variants, and transcription factor occupancy on a genome-wide scale. This review highlights the variations on the theme of the ChIP assay, the various detection methods applied downstream of ChIP, and examples of their application.

Microarray analysis of equine endometrium at days 8 and 12 of pregnancy

Establishment and maintenance of pregnancy in equids is only partially understood. To provide new insights into early events of this process, we performed a systematic analysis of transcriptome changes in the endometrium at Days 8 and 12 of pregnancy. Endometrial biopsy samples from pregnant and nonpregnant stages were taken from the same mares. Composition of the collected biopsy samples was analyzed using quantitative stereological techniques to determine proportions of surface and glandular epithelium and blood vessels. Microarray analysis did not reveal detectable changes in gene expression at Day 8, whereas at Day 12 of pregnancy 374 differentially expressed genes were identified, 332 with higher and 42 with lower transcript levels in pregnant endometrium. Expression of selected genes was validated by quantitative real-time RT-PCR. Gene set enrichment analysis, functional annotation clustering, and cocitation analysis were performed to characterize the genes differentially expressed in Day 12 pregnant endometrium. Many known estrogen-induced genes and genes involved in regulation of estrogen signaling were found, but also genes known to be regulated by progesterone and prostaglandin E2. Additionally, differential expression of a number of genes related to angiogenesis and vascular remodeling suggests an important role of this process. Furthermore, genes that probably have conserved functions across species, such as CRYAB, ERRFI1, FGF9, IGFBP2, NR2F2, STC1, and TNFSF10, were identified. This study revealed the potential target genes and pathways of conceptus-derived estrogens, progesterone, and prostaglandin E2 in the equine endometrium probably involved in the early events of establishment and maintenance of pregnancy in the mare.

The tri-nucleotide spacer sequence between estrogen response element half-sites is conserved and modulates ERalpha-mediated transcriptional responses

The estrogen response element (ERE) consensus sequence is AGGTCAnnnTGACCT, where nnn is known as the tri-nucleotide spacer sequence. Studying 1017 high-confidence ERalpha-bound loci, we found that genomic EREs are enriched for spacers composed of C(A/T)G, suggesting that the spacer may influence receptor binding and transcriptional responses. We designed consensus EREs containing variable spacer sequences and compared ERalpha binding in gel shift assays and enhancer function in reporter assays. We found that ERalpha-ERE binding affinity is modulated by the tri-nucleotide spacer sequence and is favored by spacer sequences of CTG>GCC>TTT. Similarly, luciferase reporter assays indicated that the estrogen-stimulated transcriptional response is modulated by the spacer and parallels the gel shift data: CTG>GCC>TTT. Reporter assays demonstrated that the spacer sequence also modulates the sensitivity of EREs to repression engendered by the receptor antagonist hydroxytamoxifen. These experiments indicate that the sequence of the tri-nucleotide spacer is non-random at receptor-bound genomic loci, influences ERalpha-DNA-binding affinity, and modulates transactivation potential of the receptor-ligand-DNA complex. This work has implications for understanding which genomic EREs are targeted by ERalpha, should improve computational prediction of functional EREs within genomic sequences, and describes novel sequence determinants of the estrogen response.

Genomic analyses of hormone signaling and gene regulation

Many cellular signaling pathways ultimately control specific patterns of gene expression in the nucleus through a variety of signal-regulated transcription factors (TFs), including nuclear hormone receptors (NRs). The advent of genomic technologies for examining signal-regulated transcriptional responses and TF binding on a genomic scale has dramatically increased our understanding of the cellular programs that control hormonal signaling and gene regulation. Studies of TFs, especially NRs, using genomic approaches have revealed novel and unexpected features of hormone-regulated transcription, and a global view is beginning to emerge. In this review, we discuss the genomic methodologies that have been applied to the study of hormone-regulated gene expression, the results that have been obtained from using them, and the future prospects for these approaches. Given the wealth of information about hormone-dependent gene regulation by NRs, we have focused this review on the knowledge gained from genomic studies of their function.

Activation of estrogen-responsive genes does not require their nuclear co-localization

The spatial organization of the genome in the nucleus plays a role in the regulation of gene expression. Whether co-regulated genes are subject to coordinated repositioning to a shared nuclear space is a matter of considerable interest and debate. We investigated the nuclear organization of estrogen receptor alpha (ERalpha) target genes in human breast epithelial and cancer cell lines, before and after transcriptional activation induced with estradiol. We find that, contrary to another report, the ERalpha target genes TFF1 and GREB1 are distributed in the nucleoplasm with no particular relationship to each other. The nuclear separation between these genes, as well as between the ERalpha target genes PGR and CTSD, was unchanged by hormone addition and transcriptional activation with no evidence for co-localization between alleles. Similarly, while the volume occupied by the chromosomes increased, the relative nuclear position of the respective chromosome territories was unaffected by hormone addition. Our results demonstrate that estradiol-induced ERalpha target genes are not required to co-localize in the nucleus.

Estrogen receptor alpha controls a gene network in luminal-like breast cancer cells comprising multiple transcription factors and microRNAs

Luminal-like breast tumor cells express estrogen receptor alpha (ERalpha), a member of the nuclear receptor family of ligand-activated transcription factors that controls their proliferation, survival, and functional status. To identify the molecular determinants of this hormone-responsive tumor phenotype, a comprehensive genome-wide analysis was performed in estrogen stimulated MCF-7 and ZR-75.1 cells by integrating time-course mRNA expression profiling with global mapping of genomic ERalpha binding sites by chromatin immunoprecipitation coupled to massively parallel sequencing, microRNA expression profiling, and in silico analysis of transcription units and receptor binding regions identified. All 1270 genes that were found to respond to 17beta-estradiol in both cell lines cluster in 33 highly concordant groups, each of which showed defined kinetics of RNA changes. This hormone-responsive gene set includes several direct targets of ERalpha and is organized in a gene regulation cascade, stemming from ligand-activated receptor and reaching a large number of downstream targets via AP-2gamma, B-cell activating transcription factor, E2F1 and 2, E74-like factor 3, GTF2IRD1, hairy and enhancer of split homologue-1, MYB, SMAD3, RARalpha, and RXRalpha transcription factors. MicroRNAs are also integral components of this gene regulation network because miR-107, miR-424, miR-570, miR-618, and miR-760 are regulated by 17beta-estradiol along with other microRNAs that can target a significant number of transcripts belonging to one or more estrogen-responsive gene clusters.

Global characterization of transcriptional impact of the SRC-3 coregulator

The nuclear receptor and bona fide oncogene, steroid receptor coactivator-3 (SRC-3, AIB1), acts as a master transcriptional regulator of breast cancer by transducing growth signals via the estrogen receptor alpha (ER). In this resource paper, we present the genome-wide localization analysis of SRC-3 chromatin affinity sites in MCF-7 human breast cancer chromatin and compare the cis binding sites to global cartographies for ER and FoxA1. By correlating their gene proximal binding sites to integrated gene expression signatures, and in combination with gene ontology analyses, we provide a functional classification of estradiol-induced gene regulation that further highlights an intricate transcriptional control of interdependent cellular pathways by SRC-3. Furthermore, by presenting proteomics analyses of in vivo SRC-3- and ER-associated proteins, we give strong evidence to support the idea that the interpretative power of SRC-3 in estrogen signaling is mediated through the formation of distinct, cell state-dependent protein complexes. Altogether, we present the first approach in complementary comparative analyses that converges results obtained by three discovery-driven methods (cistromics, transcriptomics, and proteomics) into testable hypotheses, thus providing a valuable resource for follow-up studies that further our understanding of estrogen signaling in human diseases in general and breast cancer in particular.

Estrogen receptor-dependent regulation of CYP2B6 in human breast cancer cells

Estrogen receptor alpha (ERalpha) mediates the biological actions of estrogens and also contributes to the development and progression of breast cancer. To gain a more comprehensive understanding of ERalpha-mediated transcription, we used chromatin immunoprecipitation and promoter focused microarrays (ChIP-chip) to identify ERalpha binding sites in T-47D human breast cancer cells. Transcription factor binding site analysis revealed that the estrogen response element (ERE) was significantly over-represented and was found in 50% of the 243 ERalpha-bound regions identified. Interestingly, multiple ERalpha-bound regions were detected in the upstream regulatory sequences of the CYP2B gene cluster. Because ERalpha has been reported to regulate the expression of other cytochrome P450 enzymes and CYP2B6 is highly expressed in ERalpha-positive breast tumors, we focused on characterizing the ERalpha-dependent regulation of CYP2B6. Reporter gene assays revealed that ERalpha and ERbeta increased CYP2B6-regulated gene expression through a functional ERE located at -1669 to -1657 in the upstream regulatory region of CYP2B6. E2 increased ERalpha and nuclear receptor coactivator 3 (NCoA3) recruitment to the 5'-flanking region of CYP2B6, and increased CYP2B6 mRNA levels in T-47D but not in MCF-7 human breast cancer cells. RNAi-mediated knockdown of ERalpha in the T-47D cells resulted in a significant decrease in CYP2B6 mRNA levels. Taken together, our study provides evidence for cell-type specific transcriptional regulation of the CYP2B6 gene by ERs.

Location analysis for the estrogen receptor-alpha reveals binding to diverse ERE sequences and widespread binding within repetitive DNA elements

Location analysis for estrogen receptor-α (ERα)-bound cis-regulatory elements was determined in MCF7 cells using chromatin immunoprecipitation (ChIP)-on-chip. Here, we present the estrogen response element (ERE) sequences that were identified at ERα-bound loci and quantify the incidence of ERE sequences under two stringencies of detection: <10% and 10–20% nucleotide deviation from the canonical ERE sequence. We demonstrate that ∼50% of all ERα-bound loci do not have a discernable ERE and show that most ERα-bound EREs are not perfect consensus EREs. Approximately one-third of all ERα-bound ERE sequences reside within repetitive DNA sequences, most commonly of the AluS family. In addition, the 3-bp spacer between the inverted ERE half-sites, rather than being random nucleotides, is C(A/T)G-enriched at bona fide receptor targets. Diverse ERα-bound loci were validated using electrophoretic mobility shift assay and ChIP-polymerase chain reaction (PCR). The functional significance of receptor-bound loci was demonstrated using luciferase reporter assays which proved that repetitive element ERE sequences contribute to enhancer function. ChIP-PCR demonstrated estrogen-dependent recruitment of the coactivator SRC3 to these loci in vivo. Our data demonstrate that ERα binds to widely variant EREs with less sequence specificity than had previously been suspected and that binding at repetitive and nonrepetitive genomic targets is favored by specific trinucleotide spacers.

Robustness of transcriptional regulatory program influences gene expression variability

Background

Most genes are not affected when any transcription factor (TF) is knocked out, indicating that they have robust transcriptional regulatory program. Yet the mechanism underlying robust transcriptional regulatory program is less clear.

Results

Here, we studied the cause and effect of robust transcriptional regulatory program. We found that cooperative TFs in the robust transcriptional regulatory program regulate their common target genes in an activity-redundant fashion, and they are able to compensate for each other's loss. As a result, their target genes are insensitive to their single perturbation. We next revealed that the degree of robustness of transcriptional regulatory program influences gene expression variability. Genes with fragile (unrobust) transcriptional regulatory program under normal growth condition could be readily reprogrammed to significantly modulate gene expression upon changing conditions. They also have high evolutionary rates of gene expression. We further showed that the fragile transcriptional regulatory program is a major source of expression variability.

Conclusion

We showed that activity-redundant TFs guarantee the robustness of transcriptional regulatory programs, and the fragility of transcriptional regulatory program plays a major role in gene expression variability. These findings reveal the mechanisms underlying robust transcription and expression variability.