We have integrated and analyzed a large number of data sets from a variety of genomic assays using a novel computational pipeline to provide a global view of estrogen receptor 1 (ESR1; a.k.a. ERα) enhancers in MCF-7 human breast cancer cells. Using this approach, we have defined a class of primary transcripts (eRNAs) that are transcribed uni- or bidirectionally from estrogen receptor binding sites (ERBSs) with an average transcription unit length of ∼3–5 kb. The majority are up-regulated by short treatments with estradiol (i.e., 10, 25, or 40 min) with kinetics that precede or match the induction of the target genes. The production of eRNAs at ERBSs is strongly correlated with the enrichment of a number of genomic features that are associated with enhancers (e.g., H3K4me1, H3K27ac, EP300/CREBBP, RNA polymerase II, open chromatin architecture), as well as enhancer looping to target gene promoters. In the absence of eRNA production, strong enrichment of these features is not observed, even though ESR1 binding is evident. We find that flavopiridol, a CDK9 inhibitor that blocks transcription elongation, inhibits eRNA production but does not affect other molecular indicators of enhancer activity, suggesting that eRNA production occurs after the assembly of active enhancers. Finally, we show that an enhancer transcription “signature” based on GRO-seq data can be used for de novo enhancer prediction across cell types. Together, our studies shed new light on the activity of ESR1 at its enhancer sites and provide new insights about enhancer function.

PARP inhibitors (PARPi) prevent cancer cell growth by inducing synthetic lethality with DNA repair defects (e.g., in BRCA1/2 mutant cells). We have identified an alternative pathway for PARPi-mediated growth control in BRCA1/2-intact breast cancer cells involving rDNA transcription and ribosome biogenesis. PARP-1 binds to snoRNAs, which stimulate PARP-1 catalytic activity in the nucleolus independent of DNA damage. Activated PARP-1 ADP-ribosylates DDX21, an RNA helicase that localizes to nucleoli and promotes rDNA transcription when ADP-ribosylated. Treatment with PARPi or mutation of the ADP-ribosylation sites reduces DDX21 nucleolar localization, rDNA transcription, ribosome biogenesis, protein translation, and cell growth. The salient features of this pathway are evident in xenografts in mice and human breast cancer patient samples. Elevated levels of PARP-1 and nucleolar DDX21 are associated with cancer-related outcomes. Our studies provide a mechanistic rationale for efficacy of PARPi in cancer cells lacking defects in DNA repair whose growth is inhibited by PARPi.

SUMMARY The interplay between mitogenic and proinflammatory signaling pathways play key roles in determining the phenotypes and clinical outcomes of breast cancers. Using GRO-seq in MCF-7 cells, we defined the immediate transcriptional effects of crosstalk between estradiol (E2) and TNFα, identifying a large set of target genes whose expression is rapidly altered with combined E2+TNFα treatment, but not with either agent alone. The pleiotropic effects on gene transcription in response to E2+TNFα are orchestrated by extensive remodeling of the ERα enhancer landscape in an NF-κB- and FoxA1-dependent manner. In addition, expression of the de novo and synergistically regulated genes is strongly associated with clinical outcomes in breast cancers. Together, our genomic and molecular analyses indicate that TNFα signaling, acting in pathways culminating in the redistribution of NF-κB and FoxA1 binding sites across the genome, creates latent ERα binding sites that underlie altered patterns of gene expression and clinically relevant cellular responses.

Noncoding transcription is a defining feature of active enhancers, linking transcription factor (TF) binding to the molecular mechanisms controlling gene expression. To determine the relationship between enhancer activity and biological outcomes in breast cancers, we profiled the transcriptomes (using GRO-seq and RNA-seq) and epigenomes (using ChIP-seq) of 11 different human breast cancer cell lines representing five major molecular subtypes of breast cancer, as well as two immortalized ("normal") human breast cell lines. In addition, we developed a robust and unbiased computational pipeline that simultaneously identifies putative subtype-specific enhancers and their cognate TFs by integrating the magnitude of enhancer transcription, TF mRNA expression levels, TF motif -values, and enrichment of H3K4me1 and H3K27ac. When applied across the 13 different cell lines noted above, the Total Functional Score of Enhancer Elements (TFSEE) identified key breast cancer subtype-specific TFs that act at transcribed enhancers to dictate gene expression patterns determining growth outcomes, including Forkhead TFs, FOSL1, and PLAG1. FOSL1, a Fos family TF, (1) is highly enriched at the enhancers of triple negative breast cancer (TNBC) cells, (2) acts as a key regulator of the proliferation and viability of TNBC cells, but not Luminal A cells, and (3) is associated with a poor prognosis in TNBC breast cancer patients. Taken together, our results validate our enhancer identification pipeline and reveal that enhancers transcribed in breast cancer cells direct critical gene regulatory networks that promote pathogenesis.

BackgroundEpigenetic regulators are frequently mutated or aberrantly expressed in a variety of cancers, leading to altered transcription states that result in changes in cell identity, behavior, and response to therapy.ResultsTo define alterations in epigenetic landscapes in breast cancers, we profiled the distributions of 8 key histone modifications by ChIP-Seq, as well as primary (GRO-seq) and steady state (RNA-Seq) transcriptomes, across 13 distinct cell lines that represent 5 molecular subtypes of breast cancer and immortalized human mammary epithelial cells.DiscussionUsing combinatorial patterns of distinct histone modification signals, we defined subtype-specific chromatin signatures to nominate potential biomarkers. This approach identified AFAP1-AS1 as a triple negative breast cancer-specific gene associated with cell proliferation and epithelial-mesenchymal-transition. In addition, our chromatin mapping data in basal TNBC cell lines are consistent with gene expression patterns in TCGA that indicate decreased activity of the androgen receptor pathway but increased activity of the vitamin D biosynthesis pathway.ConclusionsTogether, these datasets provide a comprehensive resource for histone modification profiles that define epigenetic landscapes and reveal key chromatin signatures in breast cancer cell line subtypes with potential to identify novel and actionable targets for treatment.Electronic supplementary materialThe online version of this article (10.1186/s12864-018-4533-0) contains supplementary material, which is available to authorized users.

Successfully fighting infection requires a properly tuned immune system. Recent epidemiological studies link exposure to pollutants that bind the aryl hydrocarbon receptor (AHR) during development with poorer immune responses later in life. Yet, how developmental triggering of AHR durably alters immune cell function remains unknown. Using a mouse model, we show that developmental activation of AHR leads to long-lasting reduction in the response of CD8+ T cells during influenza virus infection, cells critical for resolving primary infection. Combining genome-wide approaches, we demonstrate that developmental activation alters DNA methylation and gene expression patterns in isolated CD8+ T cells prior to and during infection. Altered transcriptional profiles in CD8+ T cells from developmentally exposed mice reflect changes in pathways involved in proliferation and immunoregulation, with an overall pattern that bears hallmarks of T cell exhaustion. Developmental exposure also changed DNA methylation across the genome, but differences were most pronounced following infection, where we observed inverse correlation between promoter methylation and gene expression. This points to altered regulation of DNA methylation as one mechanism by which AHR causes durable changes in T cell function. Discovering that distinct gene sets and pathways were differentially changed in developmentally exposed mice prior to and after infection further reveals that the process of CD8+ T cell activation is rendered fundamentally different by early life AHR signaling. These findings reveal a novel role for AHR in the developing immune system: regulating DNA methylation and gene expression as T cells respond to infection later in life.

Although many features of active transcriptional enhancers have been defined by genomic assays, we lack a clear understanding of the order of events leading to enhancer formation and activation as well as the dynamics of coregulator interactions within the enhancer complex. Here, we used selective loss-or gain-of-function mutants of estrogen receptor α (ERα) to define two distinct phases of ligand-dependent enhancer formation. In the first phase (0-20 min), p300 is recruited to ERα by Mediator as well as p300's acetylhistone-binding bromodomain to promote initial enhancer formation, which is not competent for sustained activation. In the second phase (20-45 min), p300 is recruited to ERα by steroid receptor coregulators (SRCs) for enhancer maturation and maintenance. Successful transition between these two phases ("coregulator switching") is required for proper enhancer function. Failure to recruit p300 during either phase leads to abortive enhancer formation and a lack of target gene expression. Our results reveal an ordered and cooperative assembly of ERα enhancers requiring functional interplay among p300, Mediator, and SRCs, which has implications for hormone-dependent gene regulation in breast cancers. More broadly, our results demonstrate the unexpectedly dynamic nature of coregulator interactions within enhancer complexes, which are likely to be a defining feature of all enhancers.

The bromodomain family member proteins (BRD; BET proteins) are key coregulators for estrogen receptor alpha (ERa)-mediated transcriptional enhancers. The use of BRDselective inhibitors has gained much attention as a potential treatment for various solid tumors, including ER-positive breast cancers. However, the roles of individual BET family members have largely remained unexplored. Here, we describe the role of BRDs in estrogen (E2)-dependent gene expression in ERa-positive breast cancer cells. We observed that chemical inhibition of BET family proteins with JQ1 impairs E2regulated gene expression and growth in breast cancer cells. In addition, RNAi-mediated depletion of each BET family member (BRDs 2, 3, and 4) revealed partially redundant roles at ERa enhancers and for target gene transcription. Furthermore , we found a unique role of BRD3 as a molecular sensor of total BET family protein levels and activity through compensatory control of its own protein levels. Finally, we observed that BRD3 is recruited to a subset of ERa-binding sites (ERBS) that are enriched for active enhancer features, located in clusters of ERBSs likely functioning as "super enhancers," and associated with highly E2-responsive genes. Collectively, our results illustrate a critical and specific role for BET family members in ERa-dependent gene transcription. Implications: BRD3 is recruited to and controls the activity of a subset ERa transcriptional enhancers, providing a therapeutic opportunity to target BRD3 with BET inhibitors in ERapositive breast cancers.