Relaxin Modulates the Genomic Actions and Biological Effects of Estrogen in the Myometrium by Reducing Estrogen Receptor Alpha Phosphorylation.

Estradiol (E2) and relaxin (Rln) are steroid and polypeptide hormones, respectively, with important roles in the female reproductive tract, including myometrium. Some actions of Rln, which are mediated by its membrane receptor RXFP1, require or are augmented by E2 signaling through its cognate nuclear steroid receptor, estrogen receptor alpha (ERα). In contrast, other actions of Rln act in opposition to the effects of E2. Here we explore the molecular and genomic mechanisms that underlie the functional interplay between E2 and Rln in the myometrium. We used both ovariectomized female mice and immortalized human myometrial cells expressing wild type or mutant ERα (hTERT-HM-ERα cells). Our results indicate that Rln attenuates the genomic actions and biological effects of estrogen in the myometrium and myometrial cells by reducing phosphorylation ERα on serine 118 (S118). Interestingly, we observed a potent inhibitory effect of Rln on the E2-dependent binding of ERα across the genome. The reduction in ERα binding was associated with changes in the hormone-regulated transcriptome, including a decrease in the E2-dependent expression of neighboring genes. The inhibitory effects of Rln cotreatment on the E2-dependent phosphorylation of ERα required the nuclear dual-specificity phosphatases DUSP1 and DUSP5. Moreover, the inhibitory effects of Rln were reflected in a concomitant inhibition of the E2-dependent contraction of myometrial cells. Collectively, our results identify a pathway that integrates Rln/RXFP1 and E2/ERα signaling, resulting in a convergence of membrane and nuclear signaling pathways to control genomic and biological outcomes.

Genome-wide identification of transcriptional enhancers during human placental development and association with function, differentiation, and disease†.

The placenta is a dynamic organ that must perform a remarkable variety of functions during its relatively short existence in order to support a developing fetus. These functions include nutrient delivery, gas exchange, waste removal, hormone production, and immune barrier protection. Proper placenta development and function are critical for healthy pregnancy outcomes, but the underlying genomic regulatory events that control this process remain largely unknown. We hypothesized that mapping sites of transcriptional enhancer activity and associated changes in gene expression across gestation in human placenta tissue would identify genomic loci and predicted transcription factor activity related to critical placenta functions. We used a suite of genomic assays [i.e., RNA-sequencing (RNA-seq), Precision run-on-sequencing (PRO-seq), and Chromatin immunoprecipitation-sequencing (ChIP-seq)] and computational pipelines to identify a set of >20 000 enhancers that are active at various time points in gestation. Changes in the activity of these enhancers correlate with changes in gene expression. In addition, some of these enhancers encode risk for adverse pregnancy outcomes. We further show that integrating enhancer activity, transcription factor motif analysis, and transcription factor expression can identify distinct sets of transcription factors predicted to be more active either in early pregnancy or at term. Knockdown of selected identified transcription factors in a trophoblast stem cell culture model altered the expression of key placental marker genes. These observations provide a framework for future mechanistic studies of individual enhancer-transcription factor-target gene interactions and have the potential to inform genetic risk prediction for adverse pregnancy outcomes.

Oncohistone Mutations Occur at Functional Sites of Regulatory ADP-Ribosylation.

Recent studies have identified cancer-associated mutations in histone genes that lead to the expression of mutant versions of core histones called oncohistones. Many oncohistone mutations occur at Asp and Glu residues, two amino acids known to be ADP-ribosylated (ADPRylated) by PARP1. We screened 25 Glu or Asp oncohistone mutants for their effects on cell growth in breast and ovarian cancer cells. Ectopic expression of six mutants of three different core histones (H2B, H3, and H4) altered cell growth in at least two different cell lines. Two of these sites, H2B-D51 and H4-D68, were indeed sites of ADPRylation in wild-type (unmutated) histones, and mutation of these sites inhibited ADPRylation. Mutation of H2B-D51 dramatically altered chromatin accessibility at enhancers and promoters, as well as gene expression outcomes, whereas mutation of H4-D68 did not. Additional biochemical, cellular, proteomic, and genomic analyses demonstrated that ADPRylation of H2B-D51 inhibits p300-mediated acetylation of H2B at many Lys residues. In breast cancer cell xenografts in mice, H2B-D51A promoted tumor growth, but did not confer resistance to the cytotoxic effects of PARP inhibition. Collectively, these results demonstrate that functional Asp and Glu ADPRylation sites on histones are mutated in cancers, allowing cancer cells to escape the growth-regulating effects of post-translational modifications via distinct mechanisms. SIGNIFICANCE: This study identifies cancer-driving mutations in histones as sites of PARP1-mediated ADP-ribosylation in breast and ovarian cancers, providing a molecular pathway by which cancers may subvert the growth-regulating effects of PARP1.

PARP-1 Regulates Estrogen-Dependent Gene Expression in Estrogen Receptor α-Positive Breast Cancer Cells.

Poly(ADP-ribose) polymerase-1 (PARP-1) has gained considerable attention as a target for therapeutic inhibitors in breast cancers. Previously we showed that PARP-1 localizes to active gene promoters to regulate histone methylation and RNA polymerase II activity (Pol II), altering the expression of various tumor-related genes. Here we report a role for PARP-1 in estrogen-dependent transcription in estrogen receptor alpha (ERα)-positive (ER+) breast cancers. Global nuclear run-on and sequencing analyses functionally linked PARP-1 to the direct control of estrogen-regulated gene expression in ER+ MCF-7 breast cancer cells by promoting transcriptional elongation by Pol II. Furthermore, chromatin immunoprecipitation sequencing analyses revealed that PARP-1 regulates the estrogen-dependent binding of ERα and FoxA1 to a subset of genomic ERα binding sites, promoting active enhancer formation. Moreover, we found that the expression levels of the PARP-1- and estrogen-coregulated gene set are enriched in the luminal subtype of breast cancer, and high PARP-1 expression in ER+ cases correlates with poor survival. Finally, treatment with a PARP inhibitor or a transcriptional elongation inhibitor attenuated estrogen-dependent growth of multiple ER+ breast cancer cell lines. Taken together, our results show that PARP-1 regulates critical molecular pathways that control the estrogen-dependent gene expression program underlying the proliferation of ER+ breast cancer cells. IMPLICATIONS: PARP-1 regulates the estrogen-dependent genomic binding of ERα and FoxA1 to regulate critical gene expression programs by RNA Pol II that underlie the proliferation of ER+ breast cancers, providing a potential therapeutic opportunity for PARP inhibitors in estrogen-responsive breast cancers.

Total Functional Score of Enhancer Elements Identifies Lineage-Specific Enhancers That Drive Differentiation of Pancreatic Cells.

The differentiation of embryonic stem cells into various lineages is highly dependent on the chromatin state of the genome and patterns of gene expression. To identify lineage-specific enhancers driving the differentiation of progenitors into pancreatic cells, we used a previously described computational framework called Total Functional Score of Enhancer Elements (TFSEE), which integrates multiple genomic assays that probe both transcriptional and epigenomic states. First, we evaluated and compared TFSEE as an enhancer-calling algorithm with enhancers called using GRO-seq-defined enhancer transcripts (method 1) versus enhancers called using histone modification ChIP-seq data (method 2). Second, we used TFSEE to define the enhancer landscape and identify transcription factors (TFs) that maintain the multipotency of a subpopulation of endodermal stem cells during differentiation into pancreatic lineages. Collectively, our results demonstrate that TFSEE is a robust enhancer-calling algorithm that can be used to perform multilayer genomic data integration to uncover cell type-specific TFs that control lineage-specific enhancers.

Distinct Roles for BET Family Members in Estrogen Receptor α Enhancer Function and Gene Regulation in Breast Cancer Cells.

The bromodomain family member proteins (BRD; BET proteins) are key coregulators for estrogen receptor alpha (ERα)-mediated transcriptional enhancers. The use of BRD-selective 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 ERα-positive breast cancer cells. We observed that chemical inhibition of BET family proteins with JQ1 impairs E2-regulated 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 ERα 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 ERα-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 ERα-dependent gene transcription. IMPLICATIONS: BRD3 is recruited to and controls the activity of a subset ERα transcriptional enhancers, providing a therapeutic opportunity to target BRD3 with BET inhibitors in ERα-positive breast cancers.

Activation of PARP-1 by snoRNAs Controls Ribosome Biogenesis and Cell Growth via the RNA Helicase DDX21.

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.

Histone modification profiling in breast cancer cell lines highlights commonalities and differences among subtypes.

BACKGROUND: Epigenetic 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. RESULTS: To 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. DISCUSSION: Using 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. CONCLUSIONS: Together, 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.

Enhancer transcription reveals subtype-specific gene expression programs controlling breast cancer pathogenesis.

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 P-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.

Dynamic assembly and activation of estrogen receptor α enhancers through coregulator switching.

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.

Computational Approaches for Mining GRO-Seq Data to Identify and Characterize Active Enhancers.

Transcriptional enhancers are DNA regulatory elements that are bound by transcription factors and act to positively regulate the expression of nearby or distally located target genes. Enhancers have many features that have been discovered using genomic analyses. Recent studies have shown that active enhancers recruit RNA polymerase II (Pol II) and are transcribed, producing enhancer RNAs (eRNAs). GRO-seq, a method for identifying the location and orientation of all actively transcribing RNA polymerases across the genome, is a powerful approach for monitoring nascent enhancer transcription. Furthermore, the unique pattern of enhancer transcription can be used to identify enhancers in the absence of any information about the underlying transcription factors. Here, we describe the computational approaches required to identify and analyze active enhancers using GRO-seq data, including data pre-processing, alignment, and transcript calling. In addition, we describe protocols and computational pipelines for mining GRO-seq data to identify active enhancers, as well as known transcription factor binding sites that are transcribed. Furthermore, we discuss approaches for integrating GRO-seq-based enhancer data with other genomic data, including target gene expression and function. Finally, we describe molecular biology assays that can be used to confirm and explore further the function of enhancers that have been identified using genomic assays. Together, these approaches should allow the user to identify and explore the features and biological functions of new cell type-specific enhancers.

Interplay between inflammatory and estrogen signaling in breast cancer.

Inflammation is known to have a paradoxical effect in cancers, in some cases promoting pathogenesis while in others inhibiting pathogenesis, depending on the cellular context. In an effort to answer a number of fundamental questions about two of the major signaling cascades that affect breast tumorigenesis and impact clinical outcome, we examined the genome-wide consequences of treating ERα-positive breast cancer cells with both estrogen and TNFα. Below, we highlight our observations, their biological significance, and how they provide a framework for understanding the molecular basis for integration of proinflammatory and estrogen signaling in breast cancer.

Linking the aryl hydrocarbon receptor with altered DNA methylation patterns and developmentally induced aberrant antiviral CD8+ T cell responses.

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.

TNFα signaling exposes latent estrogen receptor binding sites to alter the breast cancer cell transcriptome.

The interplay between mitogenic and proinflammatory signaling pathways plays 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.

Enhancer transcripts mark active estrogen receptor binding sites.

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.