Cell state-dependent chromatin targeting in NUT carcinoma

Aberrant transcriptional programming and chromatin dysregulation are common to most cancers. Whether by deranged cell signaling or environmental insult, the resulting oncogenic phenotype is typically manifested in transcriptional changes characteristic of undifferentiated cell growth. Here we analyze targeting of an oncogenic fusion protein, BRD4-NUT, composed of 2 normally independent chromatin regulators. The fusion causes the formation of large hyperacetylated genomic regions or megadomains, mis-regulation of c-MYC, and an aggressive carcinoma of squamous cell origin. Our previous work revealed largely distinct megadomain locations in different NUT carcinoma patient cell lines. To assess whether this was due to variations in individual genome sequences or epigenetic cell state, we expressed BRD4-NUT in a human stem cell model and found that megadomains formed in dissimilar patterns when comparing cells in the pluripotent state with the same cell line following induction along a mesodermal lineage. Thus, our work implicates initial cell state as the critical factor in the locations of BRD4-NUT megadomains. These results, together with our analysis of c-MYC protein-protein interactions in a patient cell line, are consistent with a cascade of chromatin misregulation underlying NUT carcinoma.

Cell state-dependent chromatin targeting in NUT carcinoma

Aberrant transcriptional programming and chromatin dysregulation are common to most cancers. Whether by deranged cell signaling or environmental insult, the resulting oncogenic phenotype is typically manifested in transcriptional changes characteristic of undifferentiated cell growth. Here we analyze targeting of an oncogenic fusion protein, BRD4-NUT, composed of two normally independent chromatin regulators. The fusion causes the formation of large hyperacetylated genomic regions or megadomains, mis-regulation of c-MYC , and an aggressive carcinoma of squamous cell origin. Our previous work revealed largely distinct megadomain locations in different NUT carcinoma patient cell lines. To assess whether this was due to variations in individual genome sequences or epigenetic cell state, we expressed BRD4-NUT in a human stem cell model and found that megadomains formed in dissimilar patterns when comparing cells in the pluripotent state with the same cell line following induction along a mesodermal lineage. Thus, our work implicates initial cell state as the critical factor in the locations of BRD4-NUT megadomains. These results, together with our analysis of c-MYC protein-protein interactions in a patient cell line, are consistent with a cascade of chromatin misregulation underlying NUT carcinoma.

Chemical Screen Identifies Diverse and Novel Histone Deacetylase Inhibitors as Repressors of NUT Function: Implications for NUT Carcinoma Pathogenesis and Treatment

NUT carcinoma (NC), characterized most commonly by the BRD4-NUTM1 fusion, is a rare, aggressive variant of squamous carcinoma with no effective treatment. BRD4-NUT drives growth and maintains the poorly differentiated state of NC by activating pro-growth genes such as MYC, through the formation of massive, hyperacetylated, superenhancer-like domains termed megadomains. BRD4-NUT-mediated hyperacetylation of chromatin is facilitated by the chromatin-targeting tandem bromodomains of BRD4, combined with NUT, which recruits the histone acetyltransferase, p300. Here, we developed a high-throughput small-molecule screen to identify inhibitors of transcriptional activation by NUT. In this dCAS9-based GFP-reporter assay, the strongest hits were diverse histone deacetylase (HDAC) inhibitors. Two structurally unrelated HDAC inhibitors, panobinostat and the novel compound, IRBM6, both repressed growth and induced differentiation of NC cells in proportion to their inhibition of NUT transcriptional activity. These two compounds repressed transcription of megadomain-associated oncogenic genes, such as MYC and SOX2, while upregulating pro-differentiation, non-megadomain-associated genes, including JUN, FOS, and key cell-cycle regulators, such as CDKN1A. The transcriptional changes correlate with depletion of BRD4-NUT from megadomains, and redistribution of the p300/CBP-associated chromatin acetylation mark, H3K27ac, away from megadomains toward regular enhancer regions previously populated by H3K27ac. In NC xenograft models, we demonstrated that suppression of tumor growth by panobinostat was comparable with that of bromodomain inhibition, and when combined they improved both survival and growth suppression. IMPLICATIONS: The findings provide mechanistic and preclinical rationale for the use of HDAC inhibitors, alone or combined with other agents, in the treatment of NUT carcinoma.

"Z4" Complex Member Fusions in NUT Carcinoma: Implications for a Novel Oncogenic Mechanism

Nuclear protein in testis (NUT) carcinoma (NC) is a rare, distinctly aggressive subtype of squamous carcinoma defined by the presence of NUT-fusion oncogenes resulting from chromosomal translocation. In most cases, the NUT gene (NUTM1) is fused to bromodomain containing 4 (BRD4) forming the BRD4-NUT oncogene. Here, a novel fusion partner to NUT was discovered using next-generation sequencing and FISH from a young patient with an undifferentiated malignant round cell tumor. Interestingly, the NUT fusion identified involved ZNF592, a zinc finger containing protein, which was previously identified as a component of the BRD4-NUT complex. In BRD4-NUT-expressing NC cells, wild-type ZNF592 and other associated "Z4" complex proteins, including ZNF532 and ZMYND8, colocalize with BRD4-NUT in characteristic nuclear foci. Furthermore, ectopic expression of BRD4-NUT in a non-NC cell line induces sequestration of Z4 factors to BRD4-NUT foci. Finally, the data demonstrate the specific dependency of NC cells on Z4 modules, ZNF532 and ZNF592. IMPLICATIONS: This study establishes the oncogenic role of Z4 factors in NC, offering potential new targeted therapeutic strategies in this incurable cancer.Visual Overview: http://mcr.aacrjournals.org/content/molcanres/16/12/1826/F1.large.jpg.

Abstract 2487: Novel ZNF-fusions in NUT carcinoma

INTRODUCTION. NUT Carcinoma (NC) is a rare, aggressive squamous carcinoma characterized by chromosome translocations involving rearrangement of the NUT gene (also known as NUTM1) on chromosome 15. While the BRD4-NUT fusion oncoprotein drives the majority of NCs, approximately 25% of NCs harbor non-BRD4 fusion partner genes, including BRD3-, NSD3-, and ZNF532-NUT. All of these partner genes encode proteins that are key components of the BRD4-NUT oncogenic complex and have led to greater understanding of the oncogenic mechanism of NC. Despite these advances, several NCs harbor NUT-fusions with genes of unknown identity, here termed NUT-variant, offering further opportunity to discover additional NUT-fusion partners that may play a key role in BRD4-NUT complex function.

METHODS. In this study, we identified 7 NCs with NUT-variant fusions. Archer® FusionPlex®, a next generation sequencing-based derivation of rapid amplification of cDNA ends (RACE), was used to identify the NUT-fusions using RNA extracted from formalin-fixed paraffin-embedded (FFPE) archival tumor tissue. Confirmatory studies included RT-PCR, Sanger sequencing, and fluorescence in situ hybridization (FISH).

RESULTS. Of the 7 NCs with unknown NUT- fusion partners, Archer® FusionPlex® identified two tumors with classic BRD4-NUT fusions, two BRD3-NUT, one non-fusion, one technical failure, and one novel ZNF592-NUT fusion. The tumor harboring the ZNF592-NUT fusion was biopsied from a pelvic tumor in an 18 year old female patient. She was initially diagnosed with NC based on positive NUT immunohistochemical staining of her pelvic tumor, and died 13 months later. RT-PCR and Sanger sequencing confirmed the ZNF592-NUT fusion detected in Archer® FusionPlex® analysis. The ZNF592-NUT fusion was further confirmed by DNA-FISH performed on the same tissue. ZNF592 is known to form a “Z4” coregulator complex with ZNF532 (a NUT-fusion partner), ZNF687 and ZMYND8. Interestingly, we have previously reported that all Z4 members co-purify with the BRD4-NUT chromatin complex. Here we show that wild type ZNF532 and wild type ZNF592 co-localize with BRD4-NUT by immunofluorescence in the BRD4-NUT NC cell line TC-797. Furthermore, knockdown of ZNF532 and ZNF592 reduced cell viability/growth of TC-797 cells, indicating these are critical factors for NC growth and/or viability.

CONCLUSION. Here we report, through the identification of a novel NUT-fusion partner, ZNF592, evidence that supports a critical role of Z4 factors in BRD4-NUT complex oncogenic function, providing the opportunity for future mechanistic insights and targeted therapeutic inhibition. Secondly, the findings suggest that Archer® FusionPlex® may in certain circumstances have a greater sensitivity for the detection of NUT-fusions than conventional FISH, and should be considered as an alternative approach for fusion-partner detection and discovery in NC.

Citation Format: Hitoshi Shiota, Janine E. Elya, Kristina Danga, Kelly Becht, Pauline Chou, Artyom A. Alekseyenko, Valentina Nardi, Christopher A. French. Novel ZNF-fusions in NUT carcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 2487.

Abstract 2655: Oncogenic chromatin factors drive cell type-specific transcription within megadomains in NUT midline carcinoma

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 find 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. Megadomains drive transcription of underlying DNA in NMC patient cells and in naïve cells induced to express BRD4-NUT. Here we characterize the constituents of BRD4-NUT chromatin complexes using a crosslinking approach, BioTAP-XL. We find many transcriptional activating proteins known to associate with BRD4, along with novel interactors including p300/CBP and a previously uncharacterized BRD4-NUT Megadomain Associated Protein (BMAP1). BMAP1 is expressed in primary NMC tissue and a subset of more common head and neck squamous cell carcinomas (HNSQC). Concurrently, we discovered a patient-derived NMC harboring a novel BMAP1-NUT fusion. BMAP-NUT blocks differentiation, and like BRD4-NUT recruits p300 to form hyperacetylated megadomains, including at the MYC locus. Thus, our proteomic and genetic approaches have converged on a novel mechanism that involves reprogramming very large regulatory regions to drive oncogenic transcription.

Citation Format: Mitzi I. Kuroda, Artyom A. Alekseyenko, Erica M. Walsh, Xin Wang, Adlai Grayson, Peter T. Hsi, Peter V. Kharchenko, Christopher A. French. Oncogenic chromatin factors drive cell type-specific transcription within megadomains in NUT midline carcinoma. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 2655.

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

In this study, Kang et al. investigate how PcG complexes form repressive chromatin domains. The findings show that Scm, a transcriptional repressor, is an important regulator of PRC1, PRC2, and transcriptional silencing and suggest that Scm coordinates PcG complexes and polymerizes, resulting in PcG silencing.

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.

BioTAP-XL: Cross-linking/Tandem Affinity Purification to Study DNA Targets, RNA, and Protein Components of Chromatin-Associated Complexes

In order to understand how chromatin complexes function in the nucleus, it is important to obtain a comprehensive picture of their protein, DNA, and RNA components, as well as their mutual interactions. This unit presents a chromatin cross-linking approach (BioTAP-XL) that utilizes a special BioTAP-tagged transgenic protein bait along with mass spectrometry to identify protein complex components, and high-throughput sequencing to identify RNA components and DNA binding sites. Full protocols are provided for Drosophila cells and for human cells in culture, along with an additional protocol for Drosophila embryos as the source material. A key element of the approach in all cases is the generation of control data from input chromatin samples.

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

Heterochromatin protein 1 (HP1a) plays conserved roles in gene silencing and heterochromatin and is also implicated in transcription, DNA replication, and repair. Using BioTAP-XL mass spectrometry and sequencing across multiple life stages of Drosophila, Alekseyenko et al. identify HP1a chromatin-associated protein and RNA interactions. They discover 13 novel candidates among the top interactions. Furthermore, HP1a selectively associates with a broad set of RNAs transcribed from repetitive regions. The validation of several novel HP1a protein interactors reveals new HP1a links to chromatin organization and function.

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.

"Jump start and gain" model for dosage compensation in Drosophila based on direct sequencing of nascent transcripts

Dosage compensation in Drosophila is mediated by the MSL complex, which increases male X-linked gene expression approximately 2-fold. The MSL complex preferentially binds the bodies of active genes on the male X, depositing H4K16ac with a 3' bias. Two models have been proposed for the influence of the MSL complex on transcription: one based on promoter recruitment of RNA polymerase II (Pol II), and a second featuring enhanced transcriptional elongation. Here, we utilize nascent RNA sequencing to document dosage compensation during transcriptional elongation. We also compare X and autosomes from published data on paused and elongating polymerase in order to assess the role of Pol II recruitment. Our results support a model for differentially regulated elongation, starting with release from 5' pausing and increasing through X-linked gene bodies. Our results highlight facilitated transcriptional elongation as a key mechanism for the coordinated regulation of a diverse set of genes.

The epigenome of evolving Drosophila neo-sex chromosomes: dosage compensation and heterochromatin formation

This study shows how young sex chromosomes have altered their chromatin structure in Drosophila, and what genomic changes have led to silencing of the Y, and hyper-transcription of the X.

Sex chromosomes originated from autosomes but have evolved a highly specialized chromatin structure. Drosophila Y chromosomes are composed entirely of silent heterochromatin, while male X chromosomes have highly accessible chromatin and are hypertranscribed as a result of dosage compensation. Here, we dissect the molecular mechanisms and functional pressures driving heterochromatin formation and dosage compensation of the recently formed neo-sex chromosomes of Drosophila miranda. We show that the onset of heterochromatin formation on the neo-Y is triggered by an accumulation of repetitive DNA. The neo-X has evolved partial dosage compensation and we find that diverse mutational paths have been utilized to establish several dozen novel binding consensus motifs for the dosage compensation complex on the neo-X, including simple point mutations at pre-binding sites, insertion and deletion mutations, microsatellite expansions, or tandem amplification of weak binding sites. Spreading of these silencing or activating chromatin modifications to adjacent regions results in massive mis-expression of neo-sex linked genes, and little correspondence between functionality of genes and their silencing on the neo-Y or dosage compensation on the neo-X. Intriguingly, the genomic regions being targeted by the dosage compensation complex on the neo-X and those becoming heterochromatic on the neo-Y show little overlap, possibly reflecting different propensities along the ancestral chromosome that formed the sex chromosome to adopt active or repressive chromatin configurations. Our findings have broad implications for current models of sex chromosome evolution, and demonstrate how mechanistic constraints can limit evolutionary adaptations. Our study also highlights how evolution can follow predictable genetic trajectories, by repeatedly acquiring the same 21-bp consensus motif for recruitment of the dosage compensation complex, yet utilizing a diverse array of random mutational changes to attain the same phenotypic outcome.

Sex chromosomes differ from non-sex chromosomes (“autosomes”) at the genomic, transcriptomic, and epigenomic level, yet the X and Y share a common evolutionary origin. The Drosophila Y chromosome is gene-poor and associated with a compact and transcriptionally inactive form of genetic material called heterochromatin. The X, in contrast, is enriched for activating chromatin marks and is consequently hyper-transcribed, a process thought to be an adaptation to decay and silencing of genes on the Y, resulting in “dosage compensation.” How sex chromosomes have altered their chromatin structure, and what genomic changes led to this dramatically different epigenetic makeup, however, has remained a mystery. By studying the genome, epigenome, and transcriptome of a species with a very recently evolved pair of sex chromosomes (the neo-X and neo-Y of a fruit fly, Drosophila miranda), we here recapitulate how both dosage compensation and heterochromatin formation evolve in Drosophila and establish several novel and important principles governing the evolution of chromatin structure. We dissect the evolutionary history of over 60 novel binding sites for the dosage compensation complex that evolved by natural selection on the neo-X within the last one million years. We show that the 21-bp consensus motifs for recruiting the dosage compensation complex were acquired by diverse molecular mechanisms along the neo-X, while the onset of heterochromatin formation is triggered by the accumulation of transposable elements, leading to silencing of adjacent neo-Y genes. We find that spreading of these chromatin modifications results in massive mis-expression of neo-sex linked genes, and that little correspondence exists between functional activity of genes on the neo-Y and whether they are dosage-compensated on the neo-X. Intriguingly, the genomic regions being targeted by the dosage compensation complex on the neo-X and those that are heterochromatic on the neo-Y show little overlap, possibly reflecting different propensities of the ancestral chromosome that formed the sex chromosome to evolve active versus repressive chromatin configurations. These findings have broad implications for current models of sex chromosome evolution.

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.

Comment on "Drosophila dosage compensation involves enhanced Pol II recruitment to male X-linked promoters"

Conrad et al. (Reports, 10 August 2012, p. 742) reported a doubling of RNA polymerase II (Pol II) occupancy at X-linked promoters to support 5' recruitment as the key mechanism for dosage compensation in Drosophila. However, they employed an erroneous data-processing step, overestimating Pol II differences. Reanalysis of the data fails to support the authors' model for dosage compensation.

Chromatin proteins captured by ChIP-mass spectrometry are linked to dosage compensation in Drosophila

X-chromosome dosage compensation by the MSL (male-specific lethal) complex is required in Drosophila melanogaster to increase gene expression from the single male X to equal that of both female X chromosomes. Instead of focusing solely on protein complexes released from DNA, here we used chromatin-interacting protein MS (ChIP-MS) to identify MSL interactions on cross-linked chromatin. We identified MSL-enriched histone modifications, including histone H4 Lys16 acetylation and histone H3 Lys36 methylation, and CG4747, a putative Lys36-trimethylated histone H3 (H3K36me3)-binding protein. CG4747 is associated with the bodies of active genes, coincident with H3K36me3, and is mislocalized in the Set2 mutant lacking H3K36me3. CG4747 loss of function in vivo results in partial mislocalization of the MSL complex to autosomes, and RNA interference experiments confirm that CG4747 and Set2 function together to facilitate targeting of the MSL complex to active genes, validating the ChIP-MS approach.

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.

Sequence-specific targeting of dosage compensation in Drosophila favors an active chromatin context

The Drosophila MSL complex mediates dosage compensation by increasing transcription of the single X chromosome in males approximately two-fold. This is accomplished through recognition of the X chromosome and subsequent acetylation of histone H4K16 on X-linked genes. Initial binding to the X is thought to occur at “entry sites” that contain a consensus sequence motif (“MSL recognition element” or MRE). However, this motif is only ∼2 fold enriched on X, and only a fraction of the motifs on X are initially targeted. Here we ask whether chromatin context could distinguish between utilized and non-utilized copies of the motif, by comparing their relative enrichment for histone modifications and chromosomal proteins mapped in the modENCODE project. Through a comparative analysis of the chromatin features in male S2 cells (which contain MSL complex) and female Kc cells (which lack the complex), we find that the presence of active chromatin modifications, together with an elevated local GC content in the surrounding sequences, has strong predictive value for functional MSL entry sites, independent of MSL binding. We tested these sites for function in Kc cells by RNAi knockdown of Sxl, resulting in induction of MSL complex. We show that ectopic MSL expression in Kc cells leads to H4K16 acetylation around these sites and a relative increase in X chromosome transcription. Collectively, our results support a model in which a pre-existing active chromatin environment, coincident with H3K36me3, contributes to MSL entry site selection. The consequences of MSL targeting of the male X chromosome include increase in nucleosome lability, enrichment for H4K16 acetylation and JIL-1 kinase, and depletion of linker histone H1 on active X-linked genes. Our analysis can serve as a model for identifying chromatin and local sequence features that may contribute to selection of functional protein binding sites in the genome.

The genomes of complex organisms encompass hundreds of millions of base pairs of DNA, and regulatory molecules must distinguish specific targets within this vast landscape. In general, regulatory factors find target genes through sequence-specific interactions with the underlying DNA. However, sequence-specific factors typically bind only a fraction of the candidate genomic regions containing their specific target sequence motif. Here we identify potential roles for chromatin environment and flanking sequence composition in helping regulatory factors find their appropriate binding sites, using targeting of the Drosophila dosage compensation complex as a model. The initial stage of dosage compensation involves binding of the Male Specific Lethal (MSL) complex to a sequence motif called the MSL recognition element [1]. Using data from a large chromatin mapping effort (the modENCODE project), we successfully identify an active chromatin environment as predictive of selective MRE binding by the MSL complex. Our study provides a framework for using genome-wide datasets to analyze and predict functional protein–DNA binding site selection.

Plasticity in patterns of histone modifications and chromosomal proteins in Drosophila heterochromatin

Eukaryotic genomes are packaged in two basic forms, euchromatin and heterochromatin. We have examined the composition and organization of Drosophila melanogaster heterochromatin in different cell types using ChIP-array analysis of histone modifications and chromosomal proteins. As anticipated, the pericentric heterochromatin and chromosome 4 are on average enriched for the "silencing" marks H3K9me2, H3K9me3, HP1a, and SU(VAR)3-9, and are generally depleted for marks associated with active transcription. The locations of the euchromatin-heterochromatin borders identified by these marks are similar in animal tissues and most cell lines, although the amount of heterochromatin is variable in some cell lines. Combinatorial analysis of chromatin patterns reveals distinct profiles for euchromatin, pericentric heterochromatin, and the 4th chromosome. Both silent and active protein-coding genes in heterochromatin display complex patterns of chromosomal proteins and histone modifications; a majority of the active genes exhibit both "activation" marks (e.g., H3K4me3 and H3K36me3) and "silencing" marks (e.g., H3K9me2 and HP1a). The hallmark of active genes in heterochromatic domains appears to be a loss of H3K9 methylation at the transcription start site. We also observe complex epigenomic profiles of intergenic regions, repeated transposable element (TE) sequences, and genes in the heterochromatic extensions. An unexpectedly large fraction of sequences in the euchromatic chromosome arms exhibits a heterochromatic chromatin signature, which differs in size, position, and impact on gene expression among cell types. We conclude that patterns of heterochromatin/euchromatin packaging show greater complexity and plasticity than anticipated. This comprehensive analysis provides a foundation for future studies of gene activity and chromosomal functions that are influenced by or dependent upon heterochromatin.

Identification of functional elements and regulatory circuits by Drosophila modENCODE

To gain insight into how genomic information is translated into cellular and developmental programs, the Drosophila model organism Encyclopedia of DNA Elements (modENCODE) project is comprehensively mapping transcripts, histone modifications, chromosomal proteins, transcription factors, replication proteins and intermediates, and nucleosome properties across a developmental time course and in multiple cell lines. We have generated more than 700 data sets and discovered protein-coding, noncoding, RNA regulatory, replication, and chromatin elements, more than tripling the annotated portion of the Drosophila genome. Correlated activity patterns of these elements reveal a functional regulatory network, which predicts putative new functions for genes, reveals stage- and tissue-specific regulators, and enables gene-expression prediction. Our results provide a foundation for directed experimental and computational studies in Drosophila and related species and also a model for systematic data integration toward comprehensive genomic and functional annotation.

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.

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.

MSL complex is attracted to genes marked by H3K36 trimethylation using a sequence-independent mechanism

In Drosophila, X chromosome dosage compensation requires the male-specific lethal (MSL) complex, which associates with actively transcribed genes on the single male X chromosome to upregulate transcription approximately 2-fold. We found that on the male X chromosome, or when MSL complex is ectopically localized to an autosome, histone H3K36 trimethylation (H3K36me3) is a strong predictor of MSL binding. We isolated mutants lacking Set2, the H3K36me3 methyltransferase, and found that Set2 is an essential gene in both sexes of Drosophila. In set2 mutant males, MSL complex maintains X specificity but exhibits reduced binding to target genes. Furthermore, recombinant MSL3 protein preferentially binds nucleosomes marked by H3K36me3 in vitro. Our results support a model in which MSL complex uses high-affinity sites to initially recognize the X chromosome and then associates with many of its targets through sequence-independent features of transcribed genes.

DNA underreplication in intercalary heterochromatin regions in polytene chromosomes of Drosophila melanogaster correlates with the formation of partial chromosomal aberrations and ectopic pairing

We studied the influence of the Suppressor of Underreplication (SuUR) gene expression on the intercalary heterochromatin (IH) regions of Drosophila melanogaster polytene chromosomes. We observed a strong positive correlation between increased SuUR expression, underreplication extent, amount of DNA truncation, and formation of ectopic contacts in IH regions. SuUR overexpression from heat shock-driven transgene results in the formation of partial chromosomal aberrations whose breakpoints map exclusively to the regions of intercalary and pericentric heterochromatin. It is important to note that all these effects are seen only if SuUR overexpression is induced during early stages of chromosome polytenization. Therefore, we developed the idea that ectopic pairing results from the joining of free DNA ends, which are formed as a consequence of underreplication.

High-resolution ChIP-chip analysis reveals that the Drosophila MSL complex selectively identifies active genes on the male X chromosome

X-chromosome dosage compensation in Drosophila requires the male-specific lethal (MSL) complex, which up-regulates gene expression from the single male X chromosome. Here, we define X-chromosome-specific MSL binding at high resolution in two male cell lines and in late-stage embryos. We find that the MSL complex is highly enriched over most expressed genes, with binding biased toward the 3' end of transcription units. The binding patterns are largely similar in the distinct cell types, with approximately 600 genes clearly bound in all three cases. Genes identified as clearly bound in one cell type and not in another indicate that attraction of MSL complex correlates with expression state. Thus, sequence alone is not sufficient to explain MSL targeting. We propose that the MSL complex recognizes most X-linked genes, but only in the context of chromatin factors or modifications indicative of active transcription. Distinguishing expressed genes from the bulk of the genome is likely to be an important function common to many chromatin organizing and modifying activities.

MSL complex associates with clusters of actively transcribed genes along the Drosophila male X chromosome

Dosage compensation in Drosophila serves as a model system for understanding the targeting of chromatin-modifying complexes to their sites of action. The MSL (male-specific lethal) complex up-regulates transcription of the single male X chromosome, thereby equalizing levels of transcription of X-linked genes between the sexes. Recruitment of the MSL complex to its binding sites on the male X chromosome requires each of the MSL proteins and at least one of the two large noncoding roX RNAs. To better understand how the MSL complex specifically targets the X chromosome, we have defined the binding using high-resolution genomic tiling arrays. Our results indicate that the MSL complex largely associates with transcribed genes that are present in clusters along the X chromosome. We hypothesize that after initial recruitment of the MSL complex to the X chromosome by unknown mechanisms, nascent transcripts or chromatin marks associated with active transcription attract the MSL complex to its final targets. Defining MSL-complex-binding sites will provide a tool for understanding functions of large noncoding RNAs that have remained elusive.

Genomic analysis of Drosophila chromosome underreplication reveals a link between replication control and transcriptional territories

In Drosophila polytene chromosomes, most late-replicating regions remain underreplicated. A loss-of-function mutant of the suppressor of underreplication [Su(UR)] gene suppresses underreplication (UR), whereas extra copies of this gene enhance the level and number of regions showing UR. By combining DNA microarray analysis with manipulation of the number of Su(UR) gene copies, we achieved genomic-scale molecular identification of 1,036 genes that are arranged in clusters located in 52 UR chromosomal regions. These regions overlap extensively (96%) but are not completely identical with late-replicating regions of mitotically dividing Kc cells in culture. Reanalysis of published gene expression profiles revealed that genomic regions defined by replication properties include clusters of coordinately expressed genes. Genomic regions that are UR in polytene chromosomes and late replicated in Kc cell chromosomes show a particularly common association with transcriptional territories that are expressed in testis/males but not ovary/females or embryos. An attractive hypothesis for future testing is that factors involved in replication control, such as SU(UR), may interact physically with those involved in epigenetic silencing of transcription territories.

Sequence-specific targeting of MSL complex regulates transcription of the roX RNA genes

In Drosophila, dosage compensation is controlled by the male-specific lethal (MSL) complex consisting of at least five proteins and two noncoding RNAs, roX1 and roX2. The roX RNAs function in targeting MSL complex to the X chromosome, and roX transgenes can nucleate spreading of the MSL complex into flanking chromatin when inserted on an autosome. An MSL-binding site (DHS, DNaseI hypersensitive site) has been identified in each roX gene. Here, we investigate the functions of the DHS using transgenic deletion analyses and reporter assays. We find that MSL interaction with the DHS counteracts constitutive repression at roX1, resulting in male-specific expression of roX1 RNA. Surprisingly, the DHS is not required for initiation of cis spreading of MSL complex, instead local transcription of roX RNAs correlates with extensive spreading.

Effect of the Suppressor of Underreplication (SuUR) Gene on Position-Effect Variegation Silencing in Drosophila melanogaster

It has been previously shown that the SuUR gene encodes a protein located in intercalary and pericentromeric heterochromatin in Drosophila melanogaster polytene chromosomes. The SuUR mutation suppresses the formation of ectopic contacts and DNA underreplication in polytene chromosomes; SuUR+ in extra doses enhances the expression of these characters. This study demonstrates that heterochromatin-dependent PEV silencing is also influenced by SuUR. The SuUR protein localizes to chromosome regions compacted as a result of PEV; the SuUR mutation suppresses DNA underreplication arising in regions of polytene chromosomes undergoing PEV. The SuUR mutation also suppresses variegation of both adult morphological characters and chromatin compaction observed in rearranged chromosomes. In contrast, SuUR+ in extra doses and its overexpression enhance variegation. Thus, SuUR affects PEV silencing in a dose-dependent manner. However, its effect is expressed weaker than that of the strong modifier Su(var)2-5.

The MSL complex levels are critical for its correct targeting to the chromosomes in Drosophila melanogaster

In Drosophila, dosage compensation requires assembly of the Male Specific Lethal (MSL) protein complex for doubling transcription of most X-linked genes in males. The recognition of the X chromosome by the MSL complex has been suggested to include initial assembly at approximately 35 chromatin entry sites and subsequent spreading of mature complexes in cis to numerous additional sites along the chromosome. To understand this process further we examined MSL patterns in a range of wild-type and mutant backgrounds producing different amounts of MSL components. Our data support a model in which MSL complex binding to the X is directed by a hierarchy of target sites that display different affinities for the MSL proteins. Chromatin entry sites differ in their ability to provide local intensive binding of complexes to adjacent regions, and need high MSL complex titers to achieve this. We also mapped a set of definite autosomal regions (approximately 70) competent to associate with the functional MSL complex in wild-type males. Overexpression of both MSL1 and MSL2 stabilizes this binding and results in inappropriate MSL binding to the chromocenter and the 4th chromosome. Thus, wild-type MSL complex titers are critical for correct targeting to the X chromosome.

Influence of the SuUR gene on intercalary heterochromatin in Drosophila melanogaster polytene chromosomes

Salivary gland polytene chromosomes of Drosophila melanogaster have a reproducible set of intercalary heterochromatin (IH) sites, characterized by late DNA replication, underreplicated DNA, breaks and frequent ectopic contacts. The SuUR mutation has been shown to suppress underreplication, and wild-type SuUR protein is found at late-replicating IH sites and in pericentric heterochromatin. Here we show that the SuUR gene influences all four IH features. The SuUR mutation leads to earlier completion of DNA replication. Using transgenic strains with two, four or six additional SuUR(+) doses (4-8xSuUR(+)) we show that wild-type SuUR is an enhancer of DNA underreplication, causing many late-replicating sites to become underreplicated. We map the underreplication sites and show that their number increases from 58 in normal strains (2xSuUR(+)) to 161 in 4-8xSuUR(+) strains. In one of these new sites (1AB) DNA polytenization decreases from 100% in the wild type to 51%-85% in the 4xSuUR (+) strain. In the 4xSuUR(+) strain, 60% of the weak points coincide with the localization of Polycomb group (PcG) proteins. At the IH region 89E1-4 (the Bithorax complex), a typical underreplication site, the degree of underreplication increases with four doses of SuUR(+) but the extent of the underreplicated region is the same as in wild type and corresponds to the region containing PcG binding sites. We conclude that the polytene chromosome regions known as IH are binding sites for SuUR protein and in many cases PcG silencing proteins. We propose that these stable silenced regions are late replicated and, in the presence of SuUR protein, become underreplicated.

Intercalary heterochromatin in Drosophila melanogaster polytene chromosomes and the problem of genetic silencing

The morphological characteristics of intercalary heterochromatin (IH) are compared with those of other types of silenced chromatin in the Drosophila melanogaster genome: pericentric heterochromatin (PH) and regions subject to position effect variegation (PEV). We conclude that IH regions in polytene chromosomes are binding sites of silencing complexes such as PcG complexes and of SuUR protein. Binding of these proteins results in the appearance of condensed chromatin and late replication of DNA, which in turn may result in DNA underreplication. IH and PH as well as regions subject to PEV have in common the condensed chromatin appearance, the localization of specific proteins, late replication, underreplication in polytene chromosomes, and ectopic pairing.

Dosage compensation and intercalary heterochromatin in X chromosomes of Drosophila melanogaster

Regions of intercalary heterochromatin (IH) are dispersed in the euchromatic arms of polytene chromosomes and share the main properties of heterochromatin, namely chromosome constrictions resulting from DNA underreplication. These constrictions are frequent on the paired X chromosomes of females, but are practically absent from the single X chromosome of males. These sex-specific differences have been proposed to reflect the different levels of transcription and chromosome compaction due to dosage compensation, which in turn may affect the degree of underreplication in IH regions. To test this hypothesis, we induced dosage compensation in females by ectopic expression of MSL-2 protein. We then measured the extent of underreplication in IH regions by determining frequencies of constrictions, or by Southern blot analysis using a fragment of the ten (a) gene which is located in IH region 11A6-9. Females transheterozygous for Sxl (fhv1)/ Sxl (f1) or carrying a constitutive msl-2 transgene are known to hypertranscribe their X chromosomes. In such females, both the frequency of constrictions and DNA underreplication were reduced. Suppression of underreplication occurs only when a complete functional MSL complex assembles on the X chromosomes. We also used three strains that carried constitutive transgenes of msl-2 with mutations in the 5' untranslated regions. These strains produced normal levels of SXL protein, but variable levels of MSL-2 protein. The SXL protein did not prevent the formation of an MSL complex in these transgenic females. We found that the extent of underreplication of ten (a) DNA in IH region 11A6-9 negatively correlates with the amount of MSL complex.

The Bithorax Complex of Drosophila melanogaster: Underreplication and morphology in polytene chromosomes

The level of polyteny of the Drosophila salivary gland chromosomes was determined throughout the chromosome region 89E1-4, the locus of the Bithorax Complex. A zone of underreplication spans the 300 kb of DNA from the Ubx to Abd-B loci. From the centromere proximal end of the complex, a 70-kb-long gradual decrease of polytenization starts with the Ubx transcription unit and, after a floor corresponding to the abd-A locus, raises gradually back to the maximum over 70 kb in the region of the Abd-B transcription unit. In flies carrying the mutation Suppressor of DNA Underreplication [Su(UR)ES], the underreplication of the Bithorax Complex is fully suppressed. In the wild type, the Bithorax Complex forms a weak point featuring thinner bands separated by clefts or constrictions. In Su(UR)ES strain in contrast, the 89E1-4 band looks like a single solid band consisting of homogenous dense material. We speculate that the wild-type Su(UR)ES protein hampers DNA replication of silenced domains and leads to their underreplication in salivary gland polytene chromosomes.

Dynamic organization of the beta-heterochromatin in the Drosophila melanogaster polytene X chromosome

Region 20 of the polytene X chromosome of Drosophila melanogaster was studied in salivary glands (SG) and pseudonurse cells (PNC) of otu mutants. In SG chromosomes the morphology of the region strongly depends on two modifiers of position effect variegation: temperature and amount of heterochromatin. It is banded in XYY males at 25 degrees C and beta-heterochromatic in X0 males at 14 degrees C, i.e. it shows dynamic transitions. In PNC chromosomes region 20 is not heterochromatic, but demonstrates a clear banding pattern. Some molecular markers of mitotic heterochromatin were localized by means of in situ hybridization on PNC chromosomes: DNA of the gene su(f) in section 20C, the nucleolar organizer and 359-bp satellite in 20F. The 359-bp satellite, which has been considered to be specific for heterochromatin of the mitotic X chromosome, was found at two additional sites on chromosome 3L, proximally to 80C. The right arm of the X chromosome in SG chromosomes was localized in the inversion In(ILR)pn2b: the telomeric HeT-A DNA and AAGAG satellite from the right arm are polytenized, having been relocated from heterochromatin to euchromatin.

Su(UR)ES: a gene suppressing DNA underreplication in intercalary and pericentric heterochromatin of Drosophila melanogaster polytene chromosomes

A genetic locus suppressing DNA underreplication in intercalary heterochromatin (IH) and pericentric heterochromatin (PH) of the polytene chromosomes of Drosophila melanogaster salivary glands, has been described. Found in the In(1)scV2 strain, the mutation, designated as Su(UR)ES, was located on chromosome 3L at position 34. 8 and cytologically mapped to region 68A3-B4. A cytological phenotype was observed in the salivary gland chromosomes of larvae homozygous and hemizygous for Su(UR)ES: (i) in the IH regions, that normally are incompletely polytenized and so they often break to form "weak points," underreplication is suppressed, breaks and ectopic contacts disappear; (ii) the degree of polytenization in PH grows higher. That is why the regions in chromosome arm basements, normally beta-heterochromatic, acquire a distinct banding pattern, i. e., become euchromatic by morphological criteria; (iii) an additional bulk of polytenized material arises between the arms of chromosome 3 to form a fragment with a typical banding pattern. Chromosome 2 PH reveals additional alpha-heterochromatin. Su(UR)ES does not affect the viability, fertility, or morphological characters of the imago, and has semidominant expression in the heterozygote and distinct maternal effect. The results obtained provide evidence that the processes leading to DNA underreplication in IH and PH are affected by the same genetic mechanism.

Alpha and beta heterochromatin in polytene chromosome 2 of Drosophila melanogaster

The formation of alpha and beta heterochromatin in chromosomes of Drosophila melanogaster was studied in salivary glands (SGs) and pseudonurse cells (PNCs). In SGs of X0, XY, XYY, XX and XXY individuals the amounts of alpha heterochromatin were similar, suggesting that the Y chromosome does not substantially contribute to alpha heterochromatin formation. Pericentric heterochromatin developed a linear sequence of blocks in PNCs, showing morphology of both alpha and beta heterochromatin. In situ hybridization with Rsp sequences (Ho clone) revealed that the most proximal heterochromatic segment of the mitotic map (region h39) formed a polytenized block in PNCs. Dot analysis showed that the clone had a hybridization rate with PNC-DNA very close to that with DNA from mainly diploid head cells, whereas the homologous SG-DNA was dramatically underrepresented. A similar increase of DNA representation in PNC was found for AAGAC satellite DNA. The mitotic region h44 was found not to polytenize in the SG chromosome, whereas in PNC chromosome 2 this region was partly polytenized and presented as an array of several blocks of alpha and beta heterochromatin. The mapping of deficiencies with proximal breakpoints in the most distal heterochromatin segments h35 in arm 2L and h46 in 2R showed that the mitotic eu-heterochromatin transitions were located in SG chromosomes distally to the polytene 40E and 41C regions, respectively. Thus, the transition zones between mitotic hetero- and euchromatin are located in banded polytene euchromatin. A scheme for dynamic organization of pericentric heterochromatin in nuclei with polytene chromosomes is proposed.