Dynamic Competition of Polycomb and Trithorax in Transcriptional Programming.

Predicting regulatory potential from primary DNA sequences or transcription factor binding patterns is not possible. However, the annotation of the genome by chromatin proteins, histone modifications, and differential compaction is largely sufficient to reveal the locations of genes and their differential activity states. The Polycomb Group (PcG) and Trithorax Group (TrxG) proteins are the central players in this cell type-specific chromatin organization. PcG function was originally viewed as being solely repressive and irreversible, as observed at the homeotic loci in flies and mammals. However, it is now clear that modular and reversible PcG function is essential at most developmental genes. Focusing mainly on recent advances, we review evidence for how PcG and TrxG patterns change dynamically during cell type transitions. The ability to implement cell type-specific transcriptional programming with exquisite fidelity is essential for normal development.

Integrated loss- and gain-of-function screens define a core network governing human embryonic stem cell behavior.

Understanding the genetic control of human embryonic stem cell function is foundational for developmental biology and regenerative medicine. Here we describe an integrated genome-scale loss- and gain-of-function screening approach to identify genetic networks governing embryonic stem cell proliferation and differentiation into the three germ layers. We identified a deep link between pluripotency maintenance and survival by showing that genetic alterations that cause pluripotency dissolution simultaneously increase apoptosis resistance. We discovered that the chromatin-modifying complex SAGA and in particular its subunit TADA2B are central regulators of pluripotency, survival, growth, and lineage specification. Joint analysis of all screens revealed that genetic alterations that broadly inhibit differentiation across multiple germ layers drive proliferation and survival under pluripotency-maintaining conditions and coincide with known cancer drivers. Our results show the power of integrated multilayer genetic screening for the robust mapping of complex genetic networks.

A MOZ-TIF2 leukemia mouse model displays KAT6-dependent H3K23 propionylation and overexpression of a set of active developmental genes.

Aberrant regulation of chromatin modifiers is a common occurrence across many cancer types, and a key priority is to determine how specific alterations of these proteins, often enzymes, can be targeted therapeutically. MOZ, a histone acyltransferase, is recurrently fused to coactivators CBP, p300, and TIF2 in cases of acute myeloid leukemia (AML). Using either pharmacological inhibition or targeted protein degradation in a mouse model for MOZ-TIF2-driven leukemia, we show that KAT6 (MOZ/MORF) enzymatic activity and the MOZ-TIF2 protein are necessary for indefinite proliferation in cell culture. MOZ-TIF2 directly regulates a small subset of genes encoding developmental transcription factors, augmenting their high expression. Furthermore, transcription levels in MOZ-TIF2 cells positively correlate with enrichment of histone H3 propionylation at lysine 23 (H3K23pr), a recently appreciated histone acylation associated with gene activation. Unexpectedly, we also show that MOZ-TIF2 and MLL-AF9 regulate transcription of unique gene sets, and their cellular models exhibit distinct sensitivities to multiple small-molecule inhibitors directed against AML pathways. This is despite the shared genetic pathways of wild-type MOZ and MLL. Overall, our data provide insight into how aberrant regulation of MOZ contributes to leukemogenesis. We anticipate that these experiments will inform future work identifying targeted therapies in the treatment of AML and other diseases involving MOZ-induced transcriptional dysregulation.

Uncoupling histone modification crosstalk by engineering lysine demethylase LSD1.

Biochemical crosstalk between two or more histone modifications is often observed in epigenetic enzyme regulation, but its functional significance in cells has been difficult to discern. Previous enzymatic studies revealed that Lys14 acetylation of histone H3 can inhibit Lys4 demethylation by lysine-specific demethylase 1 (LSD1). In the present study, we engineered a mutant form of LSD1, Y391K, which renders the nucleosome demethylase activity of LSD1 insensitive to Lys14 acetylation. K562 cells with the Y391K LSD1 CRISPR knockin show decreased expression of a set of genes associated with cellular adhesion and myeloid leukocyte activation. Chromatin profiling revealed that the cis-regulatory regions of these silenced genes display a higher level of H3 Lys14 acetylation, and edited K562 cells show diminished H3 mono-methyl Lys4 near these silenced genes, consistent with a role for enhanced LSD1 demethylase activity. These findings illuminate the functional consequences of disconnecting histone modification crosstalk for a key epigenetic enzyme.

Bivalent complexes of PRC1 with orthologs of BRD4 and MOZ/MORF target developmental genes in Drosophila.

Regulatory decisions in Drosophila require Polycomb group (PcG) proteins to maintain the silent state and Trithorax group (TrxG) proteins to oppose silencing. Since PcG and TrxG are ubiquitous and lack apparent sequence specificity, a long-standing model is that targeting occurs via protein interactions; for instance, between repressors and PcG proteins. Instead, we found that Pc-repressive complex 1 (PRC1) purifies with coactivators Fs(1)h [female sterile (1) homeotic] and Enok/Br140 during embryogenesis. Fs(1)h is a TrxG member and the ortholog of BRD4, a bromodomain protein that binds to acetylated histones and is a key transcriptional coactivator in mammals. Enok and Br140, another bromodomain protein, are orthologous to subunits of a mammalian MOZ/MORF acetyltransferase complex. Here we confirm PRC1-Br140 and PRC1-Fs(1)h interactions and identify their genomic binding sites. PRC1-Br140 bind developmental genes in fly embryos, with analogous co-occupancy of PRC1 and a Br140 ortholog, BRD1, at bivalent loci in human embryonic stem (ES) cells. We propose that identification of PRC1-Br140 "bivalent complexes" in fly embryos supports and extends the bivalency model posited in mammalian cells, in which the coexistence of H3K4me3 and H3K27me3 at developmental promoters represents a poised transcriptional state. We further speculate that local competition between acetylation and deacetylation may play a critical role in the resolution of bivalent protein complexes during development.

A sequence motif within chromatin entry sites directs MSL establishment on the Drosophila X chromosome.

The Drosophila MSL complex associates with active genes specifically on the male X chromosome to acetylate histone H4 at lysine 16 and increase expression approximately 2-fold. To date, no DNA sequence has been discovered to explain the specificity of MSL binding. We hypothesized that sequence-specific targeting occurs at "chromatin entry sites," but the majority of sites are sequence independent. Here we characterize 150 potential entry sites by ChIP-chip and ChIP-seq and discover a GA-rich MSL recognition element (MRE). The motif is only slightly enriched on the X chromosome ( approximately 2-fold), but this is doubled when considering its preferential location within or 3' to active genes (>4-fold enrichment). When inserted on an autosome, a newly identified site can direct local MSL spreading to flanking active genes. These results provide strong evidence for both sequence-dependent and -independent steps in MSL targeting of dosage compensation to the male X chromosome.

The oncogenic BRD4-NUT chromatin regulator drives aberrant transcription within large topological domains.

NUT midline carcinoma (NMC), a subtype of squamous cell cancer, is one of the most aggressive human solid malignancies known. NMC is driven by the creation of a translocation oncoprotein, BRD4-NUT, which blocks differentiation and drives growth of NMC cells. BRD4-NUT forms distinctive nuclear foci in patient tumors, which we found correlate with ∼100 unprecedented, hyperacetylated expanses of chromatin that reach up to 2 Mb in size. These "megadomains" appear to be the result of aberrant, feed-forward loops of acetylation and binding of acetylated histones that drive transcription of underlying DNA in NMC patient cells and naïve cells induced to express BRD4-NUT. Megadomain locations are typically cell lineage-specific; however, the cMYC and TP63 regions are targeted in all NMCs tested and play functional roles in tumor growth. Megadomains appear to originate from select pre-existing enhancers that progressively broaden but are ultimately delimited by topologically associating domain (TAD) boundaries. Therefore, our findings establish a basis for understanding the powerful role played by large-scale chromatin organization in normal and aberrant lineage-specific gene transcription.

Sex comb on midleg (Scm) is a functional link between PcG-repressive complexes in Drosophila.

The Polycomb group (PcG) proteins are key regulators of development in Drosophila and are strongly implicated in human health and disease. How PcG complexes form repressive chromatin domains remains unclear. Using cross-linked affinity purifications of BioTAP-Polycomb (Pc) or BioTAP-Enhancer of zeste [E(z)], we captured all PcG-repressive complex 1 (PRC1) or PRC2 core components and Sex comb on midleg (Scm) as the only protein strongly enriched with both complexes. Although previously not linked to PRC2, we confirmed direct binding of Scm and PRC2 using recombinant protein expression and colocalization of Scm with PRC1, PRC2, and H3K27me3 in embryos and cultured cells using ChIP-seq (chromatin immunoprecipitation [ChIP] combined with deep sequencing). Furthermore, we found that RNAi knockdown of Scm and overexpression of the dominant-negative Scm-SAM (sterile α motif) domain both affected the binding pattern of E(z) on polytene chromosomes. Aberrant localization of the Scm-SAM domain in long contiguous regions on polytene chromosomes revealed its independent ability to spread on chromatin, consistent with its previously described ability to oligomerize in vitro. Pull-downs of BioTAP-Scm captured PRC1 and PRC2 and additional repressive complexes, including PhoRC, LINT, and CtBP. We propose that Scm is a key mediator connecting PRC1, PRC2, and transcriptional silencing. Combined with previous structural and genetic analyses, our results strongly suggest that Scm coordinates PcG complexes and polymerizes to produce broad domains of PcG silencing.

Heterochromatin-associated interactions of Drosophila HP1a with dADD1, HIPP1, and repetitive RNAs.

Heterochromatin protein 1 (HP1a) has conserved roles in gene silencing and heterochromatin and is also implicated in transcription, DNA replication, and repair. Here we identify chromatin-associated protein and RNA interactions of HP1a by BioTAP-XL mass spectrometry and sequencing from Drosophila S2 cells, embryos, larvae, and adults. Our results reveal an extensive list of known and novel HP1a-interacting proteins, of which we selected three for validation. A strong novel interactor, dADD1 (Drosophila ADD1) (CG8290), is highly enriched in heterochromatin, harbors an ADD domain similar to human ATRX, displays selective binding to H3K9me2 and H3K9me3, and is a classic genetic suppressor of position-effect variegation. Unexpectedly, a second hit, HIPP1 (HP1 and insulator partner protein-1) (CG3680), is strongly connected to CP190-related complexes localized at putative insulator sequences throughout the genome in addition to its colocalization with HP1a in heterochromatin. A third interactor, the histone methyltransferase MES-4, is also enriched in heterochromatin. In addition to these protein-protein interactions, we found that HP1a selectively associated with a broad set of RNAs transcribed from repetitive regions. We propose that this rich network of previously undiscovered interactions will define how HP1a complexes perform their diverse functions in cells and developing organisms.

Sex-specific phenotypes of histone H4 point mutants establish dosage compensation as the critical function of H4K16 acetylation in Drosophila.

Acetylation of histone H4 at lysine 16 (H4K16) modulates nucleosome-nucleosome interactions and directly affects nucleosome binding by certain proteins. In Drosophila, H4K16 acetylation by the dosage compensation complex subunit Mof is linked to increased transcription of genes on the single X chromosome in males. Here, we analyzed Drosophila containing different H4K16 mutations or lacking Mof protein. An H4K16A mutation causes embryonic lethality in both sexes, whereas an H4K16R mutation permits females to develop into adults but causes lethality in males. The acetyl-mimic mutation H4K16Q permits both females and males to develop into adults. Complementary analyses reveal that males lacking maternally deposited and zygotically expressed Mof protein arrest development during gastrulation, whereas females of the same genotype develop into adults. Together, this demonstrates the causative role of H4K16 acetylation by Mof for dosage compensation in Drosophila and uncovers a previously unrecognized requirement for this process already during the onset of zygotic gene transcription.

Conservation and de novo acquisition of dosage compensation on newly evolved sex chromosomes in Drosophila.

Dosage compensation has arisen in response to the evolution of distinct male (XY) and female (XX) karyotypes. In Drosophila melanogaster, the MSL complex increases male X transcription approximately twofold. X-specific targeting is thought to occur through sequence-dependent binding to chromatin entry sites (CESs), followed by spreading in cis to active genes. We tested this model by asking how newly evolving sex chromosome arms in Drosophila miranda acquired dosage compensation. We found evidence for the creation of new CESs, with the analogous sequence and spacing as in D. melanogaster, providing strong support for the spreading model in the establishment of dosage compensation.

Ectopic protein interactions within BRD4-chromatin complexes drive oncogenic megadomain formation in NUT midline carcinoma.

To investigate the mechanism that drives dramatic mistargeting of active chromatin in NUT midline carcinoma (NMC), we have identified protein interactions unique to the BRD4-NUT fusion oncoprotein compared with wild-type BRD4. Using cross-linking, affinity purification, and mass spectrometry, we identified the EP300 acetyltransferase as uniquely associated with BRD4 through the NUT fusion in both NMC and non-NMC cell types. We also discovered ZNF532 associated with BRD4-NUT in NMC patient cells but not detectable in 293T cells. EP300 and ZNF532 are both implicated in feed-forward regulatory loops leading to propagation of the oncogenic chromatin complex in BRD4-NUT patient cells. Adding key functional significance to our biochemical findings, we independently discovered a ZNF532-NUT translocation fusion in a newly diagnosed NMC patient. ChIP sequencing of the major players NUT, ZNF532, BRD4, EP300, and H3K27ac revealed the formation of ZNF532-NUT-associated hyperacetylated megadomains, distinctly localized but otherwise analogous to those found in BRD4-NUT patient cells. Our results support a model in which NMC is dependent on ectopic NUT-mediated interactions between EP300 and components of BRD4 regulatory complexes, leading to a cascade of misregulation.

Combination Targeting of the Bromodomain and Acetyltransferase Active Site of p300/CBP.

p300 and CBP are highly related histone acetyltransferase (HAT) enzymes that regulate gene expression, and their dysregulation has been linked to cancer and other diseases. p300/CBP is composed of a number of domains including a HAT domain, which is inhibited by the small molecule A-485, and an acetyl-lysine binding bromodomain, which was recently found to be selectively antagonized by the small molecule I-CBP112. Here we show that the combination of I-CBP112 and A-485 can synergize to inhibit prostate cancer cell proliferation. We find that the combination confers a dramatic reduction in p300 chromatin occupancy compared to the individual effects of blocking either domain alone. Accompanying this loss of p300 on chromatin, combination treatment leads to the reduction of specific mRNAs including androgen-dependent and pro-oncogenic prostate genes such as KLK3 (PSA) and c-Myc. Consistent with p300 directly affecting gene expression, mRNAs that are significantly reduced by combination treatment also exhibit a strong reduction in p300 chromatin occupancy at their gene promoters. The relatively few mRNAs that are up-regulated upon combination treatment show no correlation with p300 occupancy. These studies provide support for the pharmacologic advantage of concurrent targeting of two domains within one key epigenetic modification enzyme.

Streamlined discovery of cross-linked chromatin complexes and associated histone modifications by mass spectrometry.

Posttranslational modifications (PTMs) are key contributors to chromatin function. The ability to comprehensively link specific histone PTMs with specific chromatin factors would be an important advance in understanding the functions and genomic targeting mechanisms of those factors. We recently introduced a cross-linked affinity technique, BioTAP-XL, to identify chromatin-bound protein interactions that can be difficult to capture with native affinity techniques. However, BioTAP-XL was not strictly compatible with similarly comprehensive analyses of associated histone PTMs. Here we advance BioTAP-XL by demonstrating the ability to quantify histone PTMs linked to specific chromatin factors in parallel with the ability to identify nonhistone binding partners. Furthermore we demonstrate that the initially published quantity of starting material can be scaled down orders of magnitude without loss in proteomic sensitivity. We also integrate hydrophilic interaction chromatography to mitigate detergent carryover and improve liquid chromatography-mass spectrometric performance. In summary, we greatly extend the practicality of BioTAP-XL to enable comprehensive identification of protein complexes and their local chromatin environment.

X chromosome dosage compensation via enhanced transcriptional elongation in Drosophila.

The evolution of sex chromosomes has resulted in numerous species in which females inherit two X chromosomes but males have a single X, thus requiring dosage compensation. MSL (Male-specific lethal) complex increases transcription on the single X chromosome of Drosophila males to equalize expression of X-linked genes between the sexes. The biochemical mechanisms used for dosage compensation must function over a wide dynamic range of transcription levels and differential expression patterns. It has been proposed that the MSL complex regulates transcriptional elongation to control dosage compensation, a model subsequently supported by mapping of the MSL complex and MSL-dependent histone 4 lysine 16 acetylation to the bodies of X-linked genes in males, with a bias towards 3' ends. However, experimental analysis of MSL function at the mechanistic level has been challenging owing to the small magnitude of the chromosome-wide effect and the lack of an in vitro system for biochemical analysis. Here we use global run-on sequencing (GRO-seq) to examine the specific effect of the MSL complex on RNA Polymerase II (RNAP II) on a genome-wide level. Results indicate that the MSL complex enhances transcription by facilitating the progression of RNAP II across the bodies of active X-linked genes. Improving transcriptional output downstream of typical gene-specific controls may explain how dosage compensation can be imposed on the diverse set of genes along an entire chromosome.

Comprehensive analysis of the chromatin landscape in Drosophila melanogaster.

Chromatin is composed of DNA and a variety of modified histones and non-histone proteins, which have an impact on cell differentiation, gene regulation and other key cellular processes. Here we present a genome-wide chromatin landscape for Drosophila melanogaster based on eighteen histone modifications, summarized by nine prevalent combinatorial patterns. Integrative analysis with other data (non-histone chromatin proteins, DNase I hypersensitivity, GRO-Seq reads produced by engaged polymerase, short/long RNA products) reveals discrete characteristics of chromosomes, genes, regulatory elements and other functional domains. We find that active genes display distinct chromatin signatures that are correlated with disparate gene lengths, exon patterns, regulatory functions and genomic contexts. We also demonstrate a diversity of signatures among Polycomb targets that include a subset with paused polymerase. This systematic profiling and integrative analysis of chromatin signatures provides insights into how genomic elements are regulated, and will serve as a resource for future experimental investigations of genome structure and function.

Reciprocal interactions of human C10orf12 and C17orf96 with PRC2 revealed by BioTAP-XL cross-linking and affinity purification.

Understanding the composition of epigenetic regulators remains an important challenge in chromatin biology. Traditional biochemical analysis of chromatin-associated complexes requires their release from DNA under conditions that can also disrupt key interactions. Here we develop a complementary approach (BioTAP-XL), in which cross-linking (XL) enhances the preservation of protein interactions and also allows the analysis of DNA targets under the same tandem affinity purification (BioTAP) regimen. We demonstrate the power of BioTAP-XL through analysis of human EZH2, a core subunit of polycomb repressive complex 2 (PRC2). We identify and validate two strong interactors, C10orf12 and C17orf96, which display enrichment with EZH2-BioTAP at levels similar to canonical PRC2 components (SUZ12, EED, MTF2, JARID2, PHF1, and AEBP2). ChIP-seq analysis of BioTAP-tagged C10orf12 or C17orf96 revealed the similarity of each binding pattern with the location of EZH2 and the H3K27me3-silencing mark, validating their physical interaction with PRC2 components. Interestingly, analysis by mass spectrometry of C10orf12 and C17orf96 interactions revealed that these proteins may be mutually exclusive PRC2 subunits that fail to interact with each other or with JARID2 and AEBP2. C10orf12, in addition, shows a strong and unexpected association with components of the EHMT1/2 complex, thus potentially connecting PRC2 to another histone methyltransferase. Similarly, results from CBX4-BioTAP protein pulldowns are consistent with reports of a diversity of PRC1 complexes. Our results highlight the importance of reciprocal analyses of multiple subunits and suggest that iterative use of BioTAP-XL has strong potential to reveal networks of chromatin-based interactions in higher organisms.

Long-range spreading of dosage compensation in Drosophila captures transcribed autosomal genes inserted on X.

Dosage compensation in Drosophila melanogaster males is achieved via targeting of male-specific lethal (MSL) complex to X-linked genes. This is proposed to involve sequence-specific recognition of the X at approximately 150-300 chromatin entry sites, and subsequent spreading to active genes. Here we ask whether the spreading step requires transcription and is sequence-independent. We find that MSL complex binds, acetylates, and up-regulates autosomal genes inserted on X, but only if transcriptionally active. We conclude that a long-sought specific DNA sequence within X-linked genes is not obligatory for MSL binding. Instead, linkage and transcription play the pivotal roles in MSL targeting irrespective of gene origin and DNA sequence.

Impact of sequencing depth in ChIP-seq experiments.

In a chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) experiment, an important consideration in experimental design is the minimum number of sequenced reads required to obtain statistically significant results. We present an extensive evaluation of the impact of sequencing depth on identification of enriched regions for key histone modifications (H3K4me3, H3K36me3, H3K27me3 and H3K9me2/me3) using deep-sequenced datasets in human and fly. We propose to define sufficient sequencing depth as the number of reads at which detected enrichment regions increase <1% for an additional million reads. Although the required depth depends on the nature of the mark and the state of the cell in each experiment, we observe that sufficient depth is often reached at <20 million reads for fly. For human, there are no clear saturation points for the examined datasets, but our analysis suggests 40-50 million reads as a practical minimum for most marks. We also devise a mathematical model to estimate the sufficient depth and total genomic coverage of a mark. Lastly, we find that the five algorithms tested do not agree well for broad enrichment profiles, especially at lower depths. Our findings suggest that sufficient sequencing depth and an appropriate peak-calling algorithm are essential for ensuring robustness of conclusions derived from ChIP-seq data.

Nature and function of insulator protein binding sites in the Drosophila genome.

Chromatin insulator elements and associated proteins have been proposed to partition eukaryotic genomes into sets of independently regulated domains. Here we test this hypothesis by quantitative genome-wide analysis of insulator protein binding to Drosophila chromatin. We find distinct combinatorial binding of insulator proteins to different classes of sites and uncover a novel type of insulator element that binds CP190 but not any other known insulator proteins. Functional characterization of different classes of binding sites indicates that only a small fraction act as robust insulators in standard enhancer-blocking assays. We show that insulators restrict the spreading of the H3K27me3 mark but only at a small number of Polycomb target regions and only to prevent repressive histone methylation within adjacent genes that are already transcriptionally inactive. RNAi knockdown of insulator proteins in cultured cells does not lead to major alterations in genome expression. Taken together, these observations argue against the concept of a genome partitioned by specialized boundary elements and suggest that insulators are reserved for specific regulation of selected genes.