Polycomb domain formation depends on short and long distance regulatory cues

Integrated bioinformatics analysis reveals that EZH2-rich domains promote transcriptional repression in cervical cancer

Cervical cancer is the third female cancer most common worldwide. The carcinogenic process involves an alteration of the mechanisms associated with transcription. Several studies have reported an oncogenic role of the polycomb complex subunit, EZH2. However, the role of EZH2 in cervical cancer is unknown. Hence, the objective of this study was to determine the role of EZH2 in transcriptional regulation in cervical cancer. The EZH2 expression and the methylation status of its promoter were analyzed in The Cancer Genome Atlas. The EZH2 enrichment profile was analyzed using chromatin immunoprecipitation with massively parallel DNA sequencing data provided by ENCODE project. The chromatin compartments were identified in the 4D Nucleome Data Portal. The functional annotation was examined in Enrichr. We report that EZH2 expression is increased in cervical cancer which is associated with hypomethylation of its promoter. EZH2 is enriched at promoter and distal intergenic regions. We identified that EZH2 defines chromatin domains enriched with H3K27me3 within repressive compartments in the HeLa-S3 cell line. Additionally, high EZH2 expression is associated with the repression of the senescent phenotype in cervical cancer patients. Our results suggest the participation of EZH2 in the generation of domains with a silencer function in cervical cancer, which regulate the expression of genes associated with cellular senescence.

Long-Distance Repression by Human Silencers: Chromatin Interactions and Phase Separation in Silencers

Three-dimensional genome organization represents an additional layer in the epigenetic regulation of gene expression. Active transcription controlled by enhancers or super-enhancers has been extensively studied. Enhancers or super-enhancers can recruit activators or co-activators to activate target gene expression through long-range chromatin interactions. Chromatin interactions and phase separation play important roles in terms of enhancer or super-enhancer functioning. Silencers are another major type of cis-regulatory element that can mediate gene regulation by turning off or reducing gene expression. However, compared to active transcription, silencer studies are still in their infancy. This review covers the current knowledge of human silencers, especially the roles of chromatin interactions and phase separation in silencers. This review also proposes future directions for human silencer studies.

H3K27me3-rich genomic regions can function as silencers to repress gene expression via chromatin interactions

The mechanisms underlying gene repression and silencers are poorly understood. Here we investigate the hypothesis that H3K27me3-rich regions of the genome, defined from clusters of H3K27me3 peaks, may be used to identify silencers that can regulate gene expression via proximity or looping. We find that H3K27me3-rich regions are associated with chromatin interactions and interact preferentially with each other. H3K27me3-rich regions component removal at interaction anchors by CRISPR leads to upregulation of interacting target genes, altered H3K27me3 and H3K27ac levels at interacting regions, and altered chromatin interactions. Chromatin interactions did not change at regions with high H3K27me3, but regions with low H3K27me3 and high H3K27ac levels showed changes in chromatin interactions. Cells with H3K27me3-rich regions knockout also show changes in phenotype associated with cell identity, and altered xenograft tumor growth. Finally, we observe that H3K27me3-rich regions-associated genes and long-range chromatin interactions are susceptible to H3K27me3 depletion. Our results characterize H3K27me3-rich regions and their mechanisms of functioning via looping.

Defining the Boundaries of Polycomb Domains in Drosophila

Polycomb group (PcG) proteins are an important group of transcriptional repressors that act by modifying chromatin. PcG target genes are covered by the repressive chromatin mark H3K27me3. Polycomb repressive complex 2 (PRC2) is a multiprotein complex that is responsible for generating H3K27me3. In Drosophila, PRC2 is recruited by Polycomb Response Elements (PREs) and then trimethylates flanking nucleosomes, spreading the H3K27me3 mark over large regions of the genome, the "Polycomb domains." What defines the boundary of a Polycomb domain? There is experimental evidence that insulators, PolII, and active transcription can all form the boundaries of Polycomb domains. Here we divide the boundaries of larval Polycomb domains into six different categories. In one category, genes are transcribed toward the Polycomb domain, where active transcription is thought to stop the spreading of H3K27me3. In agreement with this, we show that introducing a transcriptional terminator into such a transcription unit causes an extension of the Polycomb domain. Additional data suggest that active transcription of a boundary gene may restrict the range of enhancer activity of a Polycomb-regulated gene.

The multiscale effects of polycomb mechanisms on 3D chromatin folding

Polycomb group (PcG) proteins silence master regulatory genes required to properly confer cell identity during the development of both Drosophila and mammals. They may act through chromatin compaction and higher-order folding of chromatin inside the cell nucleus. During the last decade, analysis on interphase chromosome architecture discovered self-interacting regions named topologically associated domains (TADs). TADs result from the 3D chromatin folding of a succession of transcribed and repressed epigenomic domains and from loop extrusion mediated by cohesin/CTCF in mammals. Polycomb silenced chromatin constitutes one type of repressed epigenomic domains which form compacted nano-compartments inside cell nuclei. Recruitment of canonical PcG proteins on chromatin relies on initial binding to discrete elements and further spreading into large chromatin domains covered with H3K27me3. Some of these discrete elements have a bivalent nature both in mammals and Drosophila and are dynamically regulated during development. Loops can occur between them, suggesting that their interaction plays both functional and structural roles. Formation of large chromatin domains covered by H3K27me3 seems crucial for PcG silencing and PcG proteins might exert their function through compaction of these domains in both mammals and flies, rather than by directly controlling the nucleosomal accessibility of discrete regulatory elements. In addition, PcG chromatin domains interact over long genomic distances, shaping a higher-order chromatin network. Therefore, PcG silencing might rely on multiscale chromatin folding to maintain cell identity during differentiation.

Separate Polycomb Response Elements control chromatin state and activation of the vestigial gene

Patterned expression of many developmental genes is specified by transcription factor gene expression, but is thought to be refined by chromatin-mediated repression. Regulatory DNA sequences called Polycomb Response Elements (PREs) are required to repress some developmental target genes, and are widespread in genomes, suggesting that they broadly affect developmental programs. While PREs in transgenes can nucleate trimethylation on lysine 27 of the histone H3 tail (H3K27me3), none have been demonstrated to be necessary at endogenous chromatin domains. This failure is thought to be due to the fact that most endogenous H3K27me3 domains contain many PREs, and individual PREs may be redundant. In contrast to these ideas, we show here that PREs near the wing selector gene vestigial have distinctive roles at their endogenous locus, even though both PREs are repressors in transgenes. First, a PRE near the promoter is required for vestigial activation and not for repression. Second, only the distal PRE contributes to H3K27me3, but even removal of both PREs does not eliminate H3K27me3 across the vestigial domain. Thus, endogenous chromatin domains appear to be intrinsically marked by H3K27me3, and PREs appear required to enhance this chromatin modification to high levels at inactive genes.

The Ubx Polycomb response element bypasses an unpaired Fab-8 insulator via cis transvection in Drosophila

Chromatin insulators or boundary elements protect genes from regulatory activities from neighboring genes or chromatin domains. In the Drosophila Abdominal-B (Abd-B) locus, the deletion of such elements, such as Frontabdominal-7 (Fab-7) or Fab-8 led to dominant gain of function phenotypes, presumably due to the loss of chromatin barriers. Homologous chromosomes are paired in Drosophila, creating a number of pairing dependent phenomena including transvection, and whether transvection may affect the function of Polycomb response elements (PREs) and thus contribute to the phenotypes are not known. Here, we studied the chromatin barrier activity of Fab-8 and how it is affected by the zygosity of the transgene, and found that Fab-8 is able to block the silencing effect of the Ubx PRE on the DsRed reporter gene in a CTCF binding sites dependent manner. However, the blocking also depends on the zygosity of the transgene in that the barrier activity is present when the transgene is homozygous, but absent when the transgene is heterozygous. To analyze this effect, we performed chromatin immunoprecipitation and quantitative PCR (ChIP-qPCR) experiments on homozygous transgenic embryos, and found that H3K27me3 and H3K9me3 marks are restricted by Fab-8, but they spread beyond Fab-8 into the DsRed gene when the two CTCF binding sites within Fab-8 were mutated. Consistent with this, the mutation reduced H3K4me3 and RNA Pol II binding to the DsRed gene, and consequently, DsRed expression. Importantly, in heterozygous embryos, Fab-8 is unable to prevent the spread of H3K27me3 and H3K9me3 marks from crossing Fab-8 into DsRed, suggesting an insulator bypass. These results suggest that in the Abd-B locus, deletion of the insulator in one copy of the chromosome could lead to the loss of insulator activity on the homologous chromosome, and in other loci where chromosomal deletion created hemizygous regions of the genome, the chromatin barrier could be compromised. This study highlights a role of homologous chromosome pairing in the regulation of gene expression in the Drosophila genome.

How Communication Between Nucleosomes Enables Spreading and Epigenetic Memory of Histone Modifications

Nucleosomes "talk" to each other about their modification state to form extended domains of modified histones independently of the underlying DNA sequence. At the same time, DNA elements promote modification of nucleosomes in their vicinity. How do these site-specific and histone-based activities act together to regulate spreading of histone modifications along the genome? How do they enable epigenetic memory to preserve cell identity? Many models for the dynamics of repressive histone modifications emphasize the role of strong positive feedback loops, which reinforce histone modifications by recruiting histone modifiers to preexisting modifications. Recent experiments question that repressive histone modifications are self-sustained independently of their genomic context, thereby indicating that histone-based feedback is relatively weak. In the present review, current models for the dynamics of histone modifications are compared and it is suggested that limitation of histone-based feedback is key to intrinsic confinement of spreading and coexistence of short- and long-term memory at different genomic loci. See also the video abstract here: https://youtu.be/3bxr_xDEZfQ.

Formation of a Polycomb-Domain in the Absence of Strong Polycomb Response Elements

Polycomb group response elements (PREs) in Drosophila are DNA-elements that recruit Polycomb proteins (PcG) to chromatin and regulate gene expression. PREs are easily recognizable in the Drosophila genome as strong peaks of PcG-protein binding over discrete DNA fragments; many small but statistically significant PcG peaks are also observed in PcG domains. Surprisingly, in vivo deletion of the four characterized strong PREs from the PcG regulated invected-engrailed (inv-en) gene complex did not disrupt the formation of the H3K27me3 domain and did not affect inv-en expression in embryos or larvae suggesting the presence of redundant PcG recruitment mechanism. Further, the 3D-structure of the inv-en domain was only minimally altered by the deletion of the strong PREs. A reporter construct containing a 7.5kb en fragment that contains three weak peaks but no large PcG peaks forms an H3K27me3 domain and is PcG-regulated. Our data suggests a model for the recruitment of PcG-complexes to Drosophila genes via interactions with multiple, weak PREs spread throughout an H3K27me3 domain.

Generalized nucleation and looping model for epigenetic memory of histone modifications

Histone modifications can redistribute along the genome in a sequence-independent manner, giving rise to chromatin position effects and epigenetic memory. The underlying mechanisms shape the endogenous chromatin landscape and determine its response to ectopically targeted histone modifiers. Here, we simulate linear and looping-driven spreading of histone modifications and compare both models to recent experiments on histone methylation in fission yeast. We find that a generalized nucleation-and-looping mechanism describes key observations on engineered and endogenous methylation domains including intrinsic spatial confinement, independent regulation of domain size and memory, variegation in the absence of antagonists, and coexistence of short- and long-term memory at loci with weak and strong constitutive nucleation. These findings support a straightforward relationship between the biochemical properties of chromatin modifiers and the spatiotemporal modification pattern. The proposed mechanism gives rise to a phase diagram for cellular memory that may be generally applicable to explain epigenetic phenomena across different species.

Regulation of Genome Architecture and Function by Polycomb Proteins

Polycomb group (PcG) proteins dynamically define cellular identities through the epigenetic repression of key developmental regulatory genes. PcG proteins are recruited to specific regulatory elements to modify the chromatin surrounding them. In addition, they regulate the organization of their target genes in the 3D space of the nucleus, and this regulatory function of the 3D genome architecture is involved in cell differentiation and the maintenance of cellular memory. In this review we discuss recent advances in our understanding of how PcG proteins are recruited to chromatin to induce local and global changes in chromosome conformation and regulate their target genes.

Spectral identification of topological domains

MOTIVATION

Topological domains have been proposed as the backbone of interphase chromosome structure. They are regions of high local contact frequency separated by sharp boundaries. Genes within a domain often have correlated transcription. In this paper, we present a computational efficient spectral algorithm to identify topological domains from chromosome conformation data (Hi-C data). We consider the genome as a weighted graph with vertices defined by loci on a chromosome and the edge weights given by interaction frequency between two loci. Laplacian-based graph segmentation is then applied iteratively to obtain the domains at the given compactness level. Comparison with algorithms in the literature shows the advantage of the proposed strategy.

RESULTS

An efficient algorithm is presented to identify topological domains from the Hi-C matrix.

AVAILABILITY AND IMPLEMENTATION

The Matlab source code and illustrative examples are available at http://bionetworks.ccmb.med.umich.edu/

CONTACT

: indikar@med.umich.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

One, Two, Three: Polycomb Proteins Hit All Dimensions of Gene Regulation

Polycomb group (PcG) proteins contribute to the formation and maintenance of a specific repressive chromatin state that prevents the expression of genes in a particular space and time. Polycomb repressive complexes (PRCs) consist of several PcG proteins with specific regulatory or catalytic properties. PRCs are recruited to thousands of target genes, and various recruitment factors, including DNA-binding proteins and non-coding RNAs, are involved in the targeting. PcG proteins contribute to a multitude of biological processes by altering chromatin features at different scales. PcG proteins mediate both biochemical modifications of histone tails and biophysical modifications (e.g., chromatin fiber compaction and three-dimensional (3D) chromatin conformation). Here, we review the role of PcG proteins in nuclear architecture, describing their impact on the structure of the chromatin fiber, on chromatin interactions, and on the spatial organization of the genome in nuclei. Although little is known about the role of plant PcG proteins in nuclear organization, much is known in the animal field, and we highlight similarities and differences in the roles of PcG proteins in 3D gene regulation in plants and animals.

The Arabidopsis GAGA-Binding Factor BASIC PENTACYSTEINE6 Recruits the POLYCOMB-REPRESSIVE COMPLEX1 Component LIKE HETEROCHROMATIN PROTEIN1 to GAGA DNA Motifs

Polycomb-repressive complexes (PRCs) play key roles in development by repressing a large number of genes involved in various functions. Much, however, remains to be discovered about PRC-silencing mechanisms as well as their targeting to specific genomic regions. Besides other mechanisms, GAGA-binding factors in animals can guide PRC members in a sequence-specific manner to Polycomb-responsive DNA elements. Here, we show that the Arabidopsis (Arabidopsis thaliana) GAGA-motif binding factor protein basic pentacysteine6 (BPC6) interacts with like heterochromatin protein1 (LHP1), a PRC1 component, and associates with vernalization2 (VRN2), a PRC2 component, in vivo. By using a modified DNA-protein interaction enzyme-linked immunosorbant assay, we could show that BPC6 was required and sufficient to recruit LHP1 to GAGA motif-containing DNA probes in vitro. We also found that LHP1 interacts with VRN2 and, therefore, can function as a possible scaffold between BPC6 and VRN2. The lhp1-4 bpc4 bpc6 triple mutant displayed a pleiotropic phenotype, extreme dwarfism and early flowering, which disclosed synergistic functions of LHP1 and group II plant BPC members. Transcriptome analyses supported this synergy and suggested a possible function in the concerted repression of homeotic genes, probably through histone H3 lysine-27 trimethylation. Hence, our findings suggest striking similarities between animal and plant GAGA-binding factors in the recruitment of PRC1 and PRC2 components to Polycomb-responsive DNA element-like GAGA motifs, which must have evolved through convergent evolution.

A balancing act: heterochromatin protein 1a and the Polycomb group coordinate their levels to silence chromatin in Drosophila

BACKGROUND

The small non-histone protein Heterochromatin protein 1a (HP1a) plays a vital role in packaging chromatin, most notably in forming constitutive heterochromatin at the centromeres and telomeres. A second major chromatin regulating system is that of the Polycomb/trithorax groups of genes which, respectively, maintain the repressed/activated state of euchromatin. Recent analyses suggest they affect the expression of a multitude of genes, beyond the homeotics whose alteration in expression lead to their initial discovery.

RESULTS

Our data suggest that early in Drosophila development, HP1a collaborates with the Polycomb/trithorax groups of proteins to regulate gene expression and that the two chromatin systems do not act separately as convention describes. HP1a affects the levels of both the Polycomb complexes and RNA polymerase II at promoters, as assayed by chromatin immunoprecipitation analysis. Deposition of both the repressive (H3K27me3) and activating (H3K4me3) marks promoted by the Polycomb/trithorax group genes at gene promoters is affected. Additionally, depending on which parent contributes the null mutation of the HP1a gene, the levels of the H3K27me3 and H3K9me3 silencing marks at both promoters and heterochromatin are different. Changes in levels of the H3K27me3 and H3K9me3 repressive marks show a mostly reciprocal nature. The time around the mid-blastula transition, when the zygotic genome begins to be actively transcribed, appears to be a transition/decision point for setting the levels.

CONCLUSIONS

We find that HP1a, which is normally critical for the formation of constitutive heterochromatin, also affects the generation of the epigenetic marks of the Polycomb/trithorax groups of proteins, chromatin modifiers which are key to maintaining gene expression in euchromatin. At gene promoters, deposition of both the repressive H3K27me3 and activating H3K4me3 marks of histone modifications shows a dependence on HP1a. Around the mid-blastula transition, when the zygotic genome begins to be actively transcribed, a pivotal decision for the level of silencing appears to take place. This is also when the embryo organizes its genome into heterochromatin and euchromatin. A balance between the HP1a and Polycomb group silencing systems appears to be set for the chromatin types that each system will primarily regulate.

Expanding the roles of chromatin insulators in nuclear architecture, chromatin organization and genome function

Of the numerous classes of elements involved in modulating eukaryotic chromosome structure and function, chromatin insulators arguably remain the most poorly understood in their contribution to these processes in vivo. Indeed, our view of chromatin insulators has evolved dramatically since their chromatin boundary and enhancer blocking properties were elucidated roughly a quarter of a century ago as a result of recent genome-wide, high-throughput methods better suited to probing the role of these elements in their native genomic contexts. The overall theme that has emerged from these studies is that chromatin insulators function as general facilitators of higher-order chromatin loop structures that exert both physical and functional constraints on the genome. In this review, we summarize the result of recent work that supports this idea as well as a number of other studies linking these elements to a diverse array of nuclear processes, suggesting that chromatin insulators exert master control over genome organization and behavior.

Insulator function and topological domain border strength scale with architectural protein occupancy

BACKGROUND

Chromosome conformation capture studies suggest that eukaryotic genomes are organized into structures called topologically associating domains. The borders of these domains are highly enriched for architectural proteins with characterized roles in insulator function. However, a majority of architectural protein binding sites localize within topological domains, suggesting sites associated with domain borders represent a functionally different subclass of these regulatory elements. How topologically associating domains are established and what differentiates border-associated from non-border architectural protein binding sites remain unanswered questions.

RESULTS

By mapping the genome-wide target sites for several Drosophila architectural proteins, including previously uncharacterized profiles for TFIIIC and SMC-containing condensin complexes, we uncover an extensive pattern of colocalization in which architectural proteins establish dense clusters at the borders of topological domains. Reporter-based enhancer-blocking insulator activity as well as endogenous domain border strength scale with the occupancy level of architectural protein binding sites, suggesting co-binding by architectural proteins underlies the functional potential of these loci. Analyses in mouse and human stem cells suggest that clustering of architectural proteins is a general feature of genome organization, and conserved architectural protein binding sites may underlie the tissue-invariant nature of topologically associating domains observed in mammals.

CONCLUSIONS

We identify a spectrum of architectural protein occupancy that scales with the topological structure of chromosomes and the regulatory potential of these elements. Whereas high occupancy architectural protein binding sites associate with robust partitioning of topologically associating domains and robust insulator function, low occupancy sites appear reserved for gene-specific regulation within topological domains.

Are we there yet? Initial targeting of the Male-Specific Lethal and Polycomb group chromatin complexes in Drosophila

Chromatin-binding proteins must navigate the complex nuclear milieu to find their sites of action, and a constellation of protein factors and other properties are likely to influence targeting specificity. Despite considerable progress, the precise rules by which binding specificity is achieved have remained elusive. Here, we consider early targeting events for two groups of chromatin-binding complexes in Drosophila: the Male-Specific Lethal (MSL) and the Polycomb group (PcG) complexes. These two serve as models for understanding targeting, because they have been extensively studied and play vital roles in Drosophila, and their targets have been documented at high resolution. Furthermore, the proteins and biochemical properties of both complexes are largely conserved in multicellular organisms, including humans. While the MSL complex increases gene expression and PcG members repress genes, the two groups share many similarities such as the ability to modify their chromatin environment to create active or repressive domains, respectively. With legacies of in-depth genetic, biochemical and now genomic approaches, the MSL and PcG complexes will continue to provide tractable systems for understanding the recruitment of multiprotein chromatin complexes to their target loci.

YY1 DNA binding and interaction with YAF2 is essential for Polycomb recruitment

Polycomb Group (PcG) proteins are crucial for epigenetic inheritance of cell identity and are functionally conserved from Drosophila to humans. PcG proteins regulate expression of homeotic genes and are essential for axial body patterning during development. Earlier we showed that transcription factor YY1 functions as a PcG protein. YY1 also physically interacts with YAF2, a homolog of RYBP. Here we characterize the mechanism and physiologic relevance of this interaction. We found phenotypic and biochemical correction of dRYBP mutant flies by mouse YAF2 demonstrating functional conservation across species. Further biochemical analysis revealed that YAF2 bridges interaction between YY1 and the PRC1 complex. ChIP assays in HeLa cells showed that YAF2 is responsible for PcG recruitment to DNA, which is mediated by YY1 DNA binding. Knock-down of YY1 abrogated PcG recruitment, which was not compensated by exogenous YAF2 demonstrating that YY1 DNA binding is a priori necessary for Polycomb assembly on chromatin. Finally, we found that although YAF2 and RYBP regulate a similar number of Polycomb target genes, there are very few genes that are regulated by both implying functional distinction between the two proteins. We present a model of YAF2-dependent and independent PcG DNA recruitment by YY1.

Surviving an identity crisis: a revised view of chromatin insulators in the genomics era

The control of complex, developmentally regulated loci and partitioning of the genome into active and silent domains is in part accomplished through the activity of DNA-protein complexes termed chromatin insulators. Together, the multiple, well-studied classes of insulators in Drosophila melanogaster appear to be generally functionally conserved. In this review, we discuss recent genomic-scale experiments and attempt to reconcile these newer findings in the context of previously defined insulator characteristics based on classical genetic analyses and transgenic approaches. Finally, we discuss the emerging understanding of mechanisms of chromatin insulator regulation. This article is part of a Special Issue entitled: Chromatin and epigenetic regulation of animal development.

The trithorax group proteins Kismet and ASH1 promote H3K36 dimethylation to counteract Polycomb group repression in Drosophila

Members of the Polycomb group of repressors and trithorax group of activators maintain heritable states of transcription by modifying nucleosomal histones or remodeling chromatin. Although tremendous progress has been made toward defining the biochemical activities of Polycomb and trithorax group proteins, much remains to be learned about how they interact with each other and the general transcription machinery to maintain on or off states of gene expression. The trithorax group protein Kismet (KIS) is related to the SWI/SNF and CHD families of chromatin remodeling factors. KIS promotes transcription elongation, facilitates the binding of the trithorax group histone methyltransferases ASH1 and TRX to active genes, and counteracts repressive methylation of histone H3 on lysine 27 (H3K27) by Polycomb group proteins. Here, we sought to clarify the mechanism of action of KIS and how it interacts with ASH1 to antagonize H3K27 methylation in Drosophila. We present evidence that KIS promotes transcription elongation and counteracts Polycomb group repression via distinct mechanisms. A chemical inhibitor of transcription elongation, DRB, had no effect on ASH1 recruitment or H3K27 methylation. Conversely, loss of ASH1 function had no effect on transcription elongation. Mutations in kis cause a global reduction in the di- and tri-methylation of histone H3 on lysine 36 (H3K36) - modifications that antagonize H3K27 methylation in vitro. Furthermore, loss of ASH1 significantly decreases H3K36 dimethylation, providing further evidence that ASH1 is an H3K36 dimethylase in vivo. These and other findings suggest that KIS antagonizes Polycomb group repression by facilitating ASH1-dependent H3K36 dimethylation.