Alternative splicing decouples local from global PRC2 activity.

The Polycomb repressive complex 2 (PRC2) mediates epigenetic maintenance of gene silencing in eukaryotes via methylation of histone H3 at lysine 27 (H3K27). Accessory factors define two distinct subtypes, PRC2.1 and PRC2.2, with different actions and chromatin-targeting mechanisms. The mechanisms orchestrating PRC2 assembly are not fully understood. Here, we report that alternative splicing (AS) of PRC2 core component SUZ12 generates an uncharacterized isoform SUZ12-S, which co-exists with the canonical SUZ12-L isoform in virtually all tissues and developmental stages. SUZ12-S drives PRC2.1 formation and favors PRC2 dimerization. While SUZ12-S is necessary and sufficient for the repression of target genes via promoter-proximal H3K27me3 deposition, SUZ12-L maintains global H3K27 methylation levels. Mouse embryonic stem cells (ESCs) lacking either isoform exit pluripotency more slowly and fail to acquire neuronal cell identity. Our findings reveal a physiological mechanism regulating PRC2 assembly and higher-order interactions in eutherians, with impacts on H3K27 methylation and gene repression.

The Dynamic Regulatory Genome of Capsaspora and the Origin of Animal Multicellularity.

The unicellular ancestor of animals had a complex repertoire of genes linked to multicellular processes. This suggests that changes in the regulatory genome, rather than in gene innovation, were key to the origin of animals. Here, we carry out multiple functional genomic assays in Capsaspora owczarzaki, the unicellular relative of animals with the largest known gene repertoire for transcriptional regulation. We show that changing chromatin states, differential lincRNA expression, and dynamic cis-regulatory sites are associated with life cycle transitions in Capsaspora. Moreover, we demonstrate conservation of animal developmental transcription-factor networks and extensive network interconnection in this premetazoan organism. In contrast, however, Capsaspora lacks animal promoter types, and its regulatory sites are small, proximal, and lack signatures of animal enhancers. Overall, our results indicate that the emergence of animal multicellularity was linked to a major shift in genome cis-regulatory complexity, most notably the appearance of distal enhancer regulation.

Thymine DNA glycosylase regulates cell-cycle-driven p53 transcriptional control in pluripotent cells.

Cell cycle progression is linked to transcriptome dynamics and variations in the response of pluripotent cells to differentiation cues, mostly through unknown determinants. Here, we characterized the cell-cycle-associated transcriptome and proteome of mouse embryonic stem cells (mESCs) in naive ground state. We found that the thymine DNA glycosylase (TDG) is a cell-cycle-regulated co-factor of the tumor suppressor p53. Furthermore, TDG and p53 co-bind ESC-specific cis-regulatory elements and thereby control transcription of p53-dependent genes during self-renewal. We determined that the dynamic expression of TDG is required to promote the cell-cycle-associated transcriptional heterogeneity. Moreover, we demonstrated that transient depletion of TDG influences cell fate decisions during the early differentiation of mESCs. Our findings reveal an unanticipated role of TDG in promoting molecular heterogeneity during the cell cycle and highlight the central role of protein dynamics for the temporal control of cell fate during development.

In vivo temporal resolution of acute promyelocytic leukemia progression reveals a role of Klf4 in suppressing early leukemic transformation.

Genome organization plays a pivotal role in transcription, but how transcription factors (TFs) rewire the structure of the genome to initiate and maintain the programs that lead to oncogenic transformation remains poorly understood. Acute promyelocytic leukemia (APL) is a fatal subtype of leukemia driven by a chromosomal translocation between the promyelocytic leukemia (PML) and retinoic acid receptor α (RARα) genes. We used primary hematopoietic stem and progenitor cells (HSPCs) and leukemic blasts that express the fusion protein PML-RARα as a paradigm to temporally dissect the dynamic changes in the epigenome, transcriptome, and genome architecture induced during oncogenic transformation. We found that PML-RARα initiates a continuum of topologic alterations, including switches from A to B compartments, transcriptional repression, loss of active histone marks, and gain of repressive histone marks. Our multiomics-integrated analysis identifies Klf4 as an early down-regulated gene in PML-RARα-driven leukemogenesis. Furthermore, we characterized the dynamic alterations in the Klf4 cis-regulatory network during APL progression and demonstrated that ectopic Klf4 overexpression can suppress self-renewal and reverse the differentiation block induced by PML-RARα. Our study provides a comprehensive in vivo temporal dissection of the epigenomic and topological reprogramming induced by an oncogenic TF and illustrates how topological architecture can be used to identify new drivers of malignant transformation.

Differential contribution to gene expression prediction of histone modifications at enhancers or promoters.

The ChIP-seq signal of histone modifications at promoters is a good predictor of gene expression in different cellular contexts, but whether this is also true at enhancers is not clear. To address this issue, we develop quantitative models to characterize the relationship of gene expression with histone modifications at enhancers or promoters. We use embryonic stem cells (ESCs), which contain a full spectrum of active and repressed (poised) enhancers, to train predictive models. As many poised enhancers in ESCs switch towards an active state during differentiation, predictive models can also be trained on poised enhancers throughout differentiation and in development. Remarkably, we determine that histone modifications at enhancers, as well as promoters, are predictive of gene expression in ESCs and throughout differentiation and development. Importantly, we demonstrate that their contribution to the predictive models varies depending on their location in enhancers or promoters. Moreover, we use a local regression (LOESS) to normalize sequencing data from different sources, which allows us to apply predictive models trained in a specific cellular context to a different one. We conclude that the relationship between gene expression and histone modifications at enhancers is universal and different from promoters. Our study provides new insight into how histone modifications relate to gene expression based on their location in enhancers or promoters.

Distinct structural transitions of chromatin topological domains correlate with coordinated hormone-induced gene regulation.

The human genome is segmented into topologically associating domains (TADs), but the role of this conserved organization during transient changes in gene expression is not known. Here we describe the distribution of progestin-induced chromatin modifications and changes in transcriptional activity over TADs in T47D breast cancer cells. Using ChIP-seq (chromatin immunoprecipitation combined with high-throughput sequencing), Hi-C (chromosome capture followed by high-throughput sequencing), and three-dimensional (3D) modeling techniques, we found that the borders of the ∼ 2000 TADs in these cells are largely maintained after hormone treatment and that up to 20% of the TADs could be considered as discrete regulatory units where the majority of the genes are either transcriptionally activated or repressed in a coordinated fashion. The epigenetic signatures of the TADs are homogeneously modified by hormones in correlation with the transcriptional changes. Hormone-induced changes in gene activity and chromatin remodeling are accompanied by differential structural changes for activated and repressed TADs, as reflected by specific and opposite changes in the strength of intra-TAD interactions within responsive TADs. Indeed, 3D modeling of the Hi-C data suggested that the structure of TADs was modified upon treatment. The differential responses of TADs to progestins and estrogens suggest that TADs could function as "regulons" to enable spatially proximal genes to be coordinately transcribed in response to hormones.

Nucleosome-driven transcription factor binding and gene regulation.

Elucidating the global function of a transcription factor implies the identification of its target genes and genomic binding sites. The role of chromatin in this context is unclear, but the dominant view is that factors bind preferentially to nucleosome-depleted regions identified as DNaseI-hypersensitive sites (DHS). Here we show by ChIP, MNase, and DNaseI assays followed by deep sequencing that the progesterone receptor (PR) requires nucleosomes for optimal binding and function. In breast cancer cells treated with progestins, we identified 25,000 PR binding sites (PRbs). The majority of these sites encompassed several copies of the hexanucleotide TGTYCY, which is highly abundant in the genome. We found that functional PRbs accumulate around progesterone-induced genes, mainly in enhancers. Most of these sites overlap with DHS but exhibit high nucleosome occupancy. Progestin stimulation results in remodeling of these nucleosomes with displacement of histones H1 and H2A/H2B dimers. Our results strongly suggest that nucleosomes are crucial for PR binding and hormonal gene regulation.

CDK2-dependent activation of PARP-1 is required for hormonal gene regulation in breast cancer cells.

Eukaryotic gene regulation implies that transcription factors gain access to genomic information via poorly understood processes involving activation and targeting of kinases, histone-modifying enzymes, and chromatin remodelers to chromatin. Here we report that progestin gene regulation in breast cancer cells requires a rapid and transient increase in poly-(ADP)-ribose (PAR), accompanied by a dramatic decrease of cellular NAD that could have broad implications in cell physiology. This rapid increase in nuclear PARylation is mediated by activation of PAR polymerase PARP-1 as a result of phosphorylation by cyclin-dependent kinase CDK2. Hormone-dependent phosphorylation of PARP-1 by CDK2, within the catalytic domain, enhances its enzymatic capabilities. Activated PARP-1 contributes to the displacement of histone H1 and is essential for regulation of the majority of hormone-responsive genes and for the effect of progestins on cell cycle progression. Both global chromatin immunoprecipitation (ChIP) coupled with deep sequencing (ChIP-seq) and gene expression analysis show a strong overlap between PARP-1 and CDK2. Thus, progestin gene regulation involves a novel signaling pathway that connects CDK2-dependent activation of PARP-1 with histone H1 displacement. Given the multiplicity of PARP targets, this new pathway could be used for the pharmacological management of breast cancer.

Four enzymes cooperate to displace histone H1 during the first minute of hormonal gene activation.

Gene regulation by external signals requires access of transcription factors to DNA sequences of target genes, which is limited by the compaction of DNA in chromatin. Although we have gained insight into how core histones and their modifications influence this process, the role of linker histones remains unclear. Here we show that, within the first minute of progesterone action, a complex cooperation between different enzymes acting on chromatin mediates histone H1 displacement as a requisite for gene induction and cell proliferation. First, activated progesterone receptor (PR) recruits the chromatin remodeling complexes NURF and ASCOM (ASC-2 [activating signal cointegrator-2] complex) to hormone target genes. The trimethylation of histone H3 at Lys 4 by the MLL2/MLL3 subunits of ASCOM, enhanced by the hormone-induced displacement of the H3K4 demethylase KDM5B, stabilizes NURF binding. NURF facilitates the PR-mediated recruitment of Cdk2/CyclinA, which is required for histone H1 displacement. Cooperation of ATP-dependent remodeling, histone methylation, and kinase activation, followed by H1 displacement, is a prerequisite for the subsequent displacement of histone H2A/H2B catalyzed by PCAF and BAF. Chromatin immunoprecipitation (ChIP) and sequencing (ChIP-seq) and expression arrays show that H1 displacement is required for hormone induction of most hormone target genes, some of which are involved in cell proliferation.

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.

Progesterone receptor induces bcl-x expression through intragenic binding sites favoring RNA polymerase II elongation.

Steroid receptors were classically described for regulating transcription by binding to target gene promoters. However, genome-wide studies reveal that steroid receptors-binding sites are mainly located at intragenic regions. To determine the role of these sites, we examined the effect of progestins on the transcription of the bcl-x gene, where only intragenic progesterone receptor-binding sites (PRbs) were identified. We found that in response to hormone treatment, the PR is recruited to these sites along with two histone acetyltransferases CREB-binding protein (CBP) and GCN5, leading to an increase in histone H3 and H4 acetylation and to the binding of the SWI/SNF complex. Concomitant, a more relaxed chromatin was detected along bcl-x gene mainly in the regions surrounding the intragenic PRbs. PR also mediated the recruitment of the positive elongation factor pTEFb, favoring RNA polymerase II (Pol II) elongation activity. Together these events promoted the re-distribution of the active Pol II toward the 3'-end of the gene and a decrease in the ratio between proximal and distal transcription. These results suggest a novel mechanism by which PR regulates gene expression by facilitating the proper passage of the polymerase along hormone-dependent genes.

Quantitative Comparison of Multiple Chromatin Immunoprecipitation-Sequencing (ChIP-seq) Experiments with spikChIP.

The chromatin immunoprecipitation coupled with the next-generation sequencing (ChIP-seq) is a powerful technique that enables to characterize the genomic distribution of chromatin-associated proteins, histone posttranslational modifications, and histone variants. However, in the absence of a reference control for monitoring experimental and biological variations, the standard ChIP-seq scheme is unable to accurately assess changes in the abundance of chromatin targets across different experimental samples. To overcome this limitation, the combination of external spike-in material with the experimental chromatin is offered as an effective solution for quantitative comparison of ChIP-seq data across different conditions. Here, we detail (i) the experimental protocol for preparing quality control spike-in chromatin from Drosophila melanogaster cells and (ii) the computational protocol to compare ChIP-seq samples with spike-in based on the use of the spikChIP software.

MAF amplification licenses ERα through epigenetic remodelling to drive breast cancer metastasis.

MAF amplification increases the risk of breast cancer (BCa) metastasis through mechanisms that are still poorly understood yet have important clinical implications. Oestrogen-receptor-positive (ER+) BCa requires oestrogen for both growth and metastasis, albeit by ill-known mechanisms. Here we integrate proteomics, transcriptomics, epigenomics, chromatin accessibility and functional assays from human and syngeneic mouse BCa models to show that MAF directly interacts with oestrogen receptor alpha (ERα), thereby promoting a unique chromatin landscape that favours metastatic spread. We identify metastasis-promoting genes that are de novo licensed following oestrogen exposure in a MAF-dependent manner. The histone demethylase KDM1A is key to the epigenomic remodelling that facilitates the expression of the pro-metastatic MAF/oestrogen-driven gene expression program, and loss of KDM1A activity prevents this metastasis. We have thus determined that the molecular basis underlying MAF/oestrogen-mediated metastasis requires genetic, epigenetic and hormone signals from the systemic environment, which influence the ability of BCa cells to metastasize.

Pyicos: a versatile toolkit for the analysis of high-throughput sequencing data.

MOTIVATION: High-throughput sequencing (HTS) has revolutionized gene regulation studies and is now fundamental for the detection of protein-DNA and protein-RNA binding, as well as for measuring RNA expression. With increasing variety and sequencing depth of HTS datasets, the need for more flexible and memory-efficient tools to analyse them is growing. RESULTS: We describe Pyicos, a powerful toolkit for the analysis of mapped reads from diverse HTS experiments: ChIP-Seq, either punctuated or broad signals, CLIP-Seq and RNA-Seq. We prove the effectiveness of Pyicos to select for significant signals and show that its accuracy is comparable and sometimes superior to that of methods specifically designed for each particular type of experiment. Pyicos facilitates the analysis of a variety of HTS datatypes through its flexibility and memory efficiency, providing a useful framework for data integration into models of regulatory genomics. AVAILABILITY: Open-source software, with tutorials and protocol files, is available at http://regulatorygenomics.upf.edu/pyicos or as a Galaxy server at http://regulatorygenomics.upf.edu/galaxy CONTACT: eduardo.eyras@upf.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Polycomb Factor PHF19 Controls Cell Growth and Differentiation Toward Erythroid Pathway in Chronic Myeloid Leukemia Cells.

Polycomb group (PcG) of proteins are a group of highly conserved epigenetic regulators involved in many biological functions, such as embryonic development, cell proliferation, and adult stem cell determination. PHD finger protein 19 (PHF19) is an associated factor of Polycomb repressor complex 2 (PRC2), often upregulated in human cancers. In particular, myeloid leukemia cell lines show increased levels of PHF19, yet little is known about its function. Here, we have characterized the role of PHF19 in myeloid leukemia cells. We demonstrated that PHF19 depletion decreases cell proliferation and promotes chronic myeloid leukemia (CML) differentiation. Mechanistically, we have shown how PHF19 regulates the proliferation of CML through a direct regulation of the cell cycle inhibitor p21. Furthermore, we observed that MTF2, a PHF19 homolog, partially compensates for PHF19 depletion in a subset of target genes, instructing specific erythroid differentiation. Taken together, our results show that PHF19 is a key transcriptional regulator for cell fate determination and could be a potential therapeutic target for myeloid leukemia treatment.

PHF19 mediated regulation of proliferation and invasiveness in prostate cancer cells.

The Polycomb-like protein PHF19/PCL3 associates with PRC2 and mediates its recruitment to chromatin in embryonic stem cells. PHF19 is also overexpressed in many cancers. However, neither PHF19 targets nor misregulated pathways involving PHF19 are known. Here, we investigate the role of PHF19 in prostate cancer cells. We find that PHF19 interacts with PRC2 and binds to PRC2 targets on chromatin. PHF19 target genes are involved in proliferation, differentiation, angiogenesis, and extracellular matrix organization. Depletion of PHF19 triggers an increase in MTF2/PCL2 chromatin recruitment, with a genome-wide gain in PRC2 occupancy and H3K27me3 deposition. Transcriptome analysis shows that PHF19 loss promotes deregulation of key genes involved in growth, metastasis, invasion, and of factors that stimulate blood vessels formation. Consistent with this, PHF19 silencing reduces cell proliferation, while promotes invasive growth and angiogenesis. Our findings reveal a role for PHF19 in controlling the balance between cell proliferation and invasiveness in prostate cancer.

RING1B recruits EWSR1-FLI1 and cooperates in the remodeling of chromatin necessary for Ewing sarcoma tumorigenesis.

Ewing sarcoma (EwS) is an aggressive tumor that affects adolescents and young adults. EwS is defined by a chromosomal translocation, EWSR1-FLI1 being the most common, that causes genome reprogramming through remodeling of enhancers. Here, we describe an unexpected function of RING1B, which is highly expressed in EwS. While retaining its repressive activity at Polycomb developmental regulated genes, RING1B colocalizes with EWSR1-FLI1 at active enhancers. We demonstrate that RING1B is necessary for the expression of key EWSR1-FLI1 targets by facilitating oncogene recruitment to their enhancers. Knockdown of RING1B impairs growth of tumor xenografts and expression of genes regulated by EWSR1-FLI1 bound enhancers. Pharmacological inhibition of AURKB with AZD1152 increases H2Aub levels causing down-regulation of RING1B/EWSR1-FLI1 common targets. Our findings demonstrate that RING1B is a critical modulator of EWSR1-FLI1-induced chromatin remodeling, and its inhibition is a potential therapeutic strategy for the treatment of these tumors.

Convergence on chromatin of non-genomic and genomic pathways of hormone signaling.

Gene regulation by steroid hormones involves genomic and non-genomic signaling pathways and the relationship between these two pathways is unknown. Genomic actions are often mediated by binding of the ligand-activated hormone receptors to hormone responsive elements (HREs) followed by recruitment of co-regulators, remodeling of chromatin and formation of the transcription initiation complex. The non-genomic effects of steroid hormones involve the rapid and transient activation of several kinase cascades often mediated by a subpopulation of "nuclear" receptors located in the cytoplasmic side of the cell membrane. The progesterone effect on breast cancer cell proliferation involves activation of the Src/Ras/Erk cascade mediated by a specific interaction between two domains of the N-terminal half of PR and the ligand-binding domain of ERalpha. Unexpectedly, selective inhibition of Erk, or its target kinase Msk1, interferes with chromatin remodeling and blocks MMTV transcriptional activation. A complex of activated PR, Erk and Msk1 is recruited to promoter already 5 min after hormone treatment and phosphorylates histone H3 at serine 10, leading to displacement of HP1gamma, as a requisite for recruitment of Src1, chromatin remodeling complexes (hSnf2h and Brg1) and RNA polymerase II. Thus, activation of signaling cascades in the cytoplasm is essential for chromatin remodeling and transcriptional activation of a subset of steroid hormone target genes.

Promoter bivalency favors an open chromatin architecture in embryonic stem cells.

In embryonic stem cells (ESCs), developmental gene promoters are characterized by their bivalent chromatin state, with simultaneous modification by MLL2 and Polycomb complexes. Although essential for embryogenesis, bivalency is functionally not well understood. Here, we show that MLL2 plays a central role in ESC genome organization. We generate a catalog of bona fide bivalent genes in ESCs and demonstrate that loss of MLL2 leads to increased Polycomb occupancy. Consequently, promoters lose accessibility, long-range interactions are redistributed, and ESCs fail to differentiate. We pose that bivalency balances accessibility and long-range connectivity of promoters, allowing developmental gene expression to be properly modulated.

Two domains of the progesterone receptor interact with the estrogen receptor and are required for progesterone activation of the c-Src/Erk pathway in mammalian cells.

In breast cancer cells, estrogens activate the Src/Erk pathway through an interaction of the estrogen receptor alpha (ERalpha) with the SH2 domain of c-Src. Progestins have been reported to activate also this pathway either via an interaction of the progesterone receptor isoform B (PRB) with ERalpha, which itself activates c-Src, or by direct interaction of PRB with the SH3 domain of c-Src. Here we identify two domains of PRB, ERID-I and -II, mediating a direct interaction with the ligand-binding domain of ERalpha. ERID-I and ERID-II flank a proline cluster responsible for binding of PRB to c-Src. In mammalian cells, the interaction of PRB with ERalpha and the progestin activation of the Src/Erk cascade are abolished by deletion of either ERID-I or ERID-II. These regions are not required for transactivation of a progesterone-responsive reporter gene. Mutations in the proline cluster of PRB that prevent a direct interaction with c-Src do not affect the strong activation of c-Src by progestins in the presence of ERalpha. Thus, in cells with ERalpha, ERID-I and ERID-II are necessary and sufficient for progestin activation of the endogenous Src/Erk pathway.