Active and repressive chromatin are interspersed without spreading in an imprinted gene cluster in the mammalian genome

The Igf2r imprinted cluster is an epigenetic silencing model in which expression of a ncRNA silences multiple genes in cis. Here, we map a 250 kb region in mouse embryonic fibroblast cells to show that histone modifications associated with expressed and silent genes are mutually exclusive and localized to discrete regions. Expressed genes were modified at promoter regions by H3K4me3 + H3K4me2 + H3K9Ac and on putative regulatory elements flanking active promoters by H3K4me2 + H3K9Ac. Silent genes showed two types of nonoverlapping profile. One type spread over large domains of tissue-specific silent genes and contained H3K27me3 alone. A second type formed localized foci on silent imprinted gene promoters and a nonexpressed pseudogene and contained H3K9me3 + H4K20me3 +/- HP1. Thus, mammalian chromosome arms contain active chromatin interspersed with repressive chromatin resembling the type of heterochromatin previously considered a feature of centromeres, telomeres, and the inactive X chromosome.     

Epigenetic regulation of nuclear lamina-associated heterochromatin by HAT1 and the acetylation of newly synthesized histones

A central component of the epigenome is the pattern of histone post-translational modifications that play a critical role in the formation of specific chromatin states. Following DNA replication, nascent chromatin is a 1:1 mixture of parental and newly synthesized histones and the transfer of modification patterns from parental histones to new histones is a fundamental step in epigenetic inheritance. Here we report that loss of HAT1, which acetylates lysines 5 and 12 of newly synthesized histone H4 during replication-coupled chromatin assembly, results in the loss of accessibility of large domains of heterochromatin, termed HAT1-dependent Accessibility Domains (HADs). HADs are mega base-scale domains that comprise ∼10% of the mouse genome. HAT1 globally represses H3 K9 me3 levels and HADs correspond to the regions of the genome that display HAT1-dependent increases in H3 K9me3 peak density. HADs display a high degree of overlap with a subset of Lamin-Associated Domains (LADs). HAT1 is required to maintain nuclear structure and integrity. These results indicate that HAT1 and the acetylation of newly synthesized histones may be critical regulators of the epigenetic inheritance of heterochromatin and suggest a new mechanism for the epigenetic regulation of nuclear lamina-heterochromatin interactions.     

β-Globin cis-elements determine differential nuclear targeting through epigenetic modifications

Increasing evidence points to nuclear compartmentalization as a contributing mechanism for gene regulation, yet mechanisms for compartmentalization remain unclear. In this paper, we use autonomous targeting of bacterial artificial chromosome (BAC) transgenes to reveal cis requirements for peripheral targeting. Three peripheral targeting regions (PTRs) within an HBB BAC bias a competition between pericentric versus peripheral heterochromatin targeting toward the nuclear periphery, which correlates with increased H3K9me3 across the β-globin gene cluster and locus control region. Targeting to both heterochromatin compartments is dependent on Suv39H-mediated H3K9me3 methylation. In different chromosomal contexts, PTRs confer no targeting, targeting to pericentric heterochromatin, or targeting to the periphery. A combination of fluorescent in situ hybridization, BAC transgenesis, and knockdown experiments reveals that peripheral tethering of the endogenous HBB locus depends both on Suv39H-mediated H3K9me3 methylation over hundreds of kilobases surrounding HBB and on G9a-mediated H3K9me2 methylation over flanking sequences in an adjacent lamin-associated domain. Our results demonstrate that multiple cis-elements regulate the overall balance of specific epigenetic marks and peripheral gene targeting.     

The visualization of large organized chromatin domains enriched in the H3K9me2 mark within a single chromosome in a single cell

Despite considerable efforts, our understanding of the organization of higher order chromatin conformations in single cells and how these relate to chromatin marks remains poor. We have earlier invented the Chromatin In Situ Proximity (ChrISP) technique to determine proximities between chromatin fibers within a single chromosome. Here we used ChrISP to identify chromosome 11-specific hubs that are enriched in the H3K9me2 mark and that project toward the nuclear membrane in finger-like structures. Conversely, chromosome 11-specfic chromatin hubs, visualized by the presence of either H3K9me1 or H3K9me3 marks, are chromosome-wide and largely absent at the nuclear periphery. As the nuclear periphery-specific chromatin hubs were lost in the induced reduction of H3K9me2 levels, they likely represent Large Organization Chromatin in Lysine Methylation (LOCK) domains, previously identified by ChIP-seq analysis. Strikingly, the downregulation of the H3K9me2/3 marks also led to the chromosome-wide compaction of chromosome 11, suggesting a pleiotropic function of these features not recognized before. The ChrISP-mediated visualization of dynamic chromatin states in single cells thus provides an analysis of chromatin structures with a resolution far exceeding that of any other light microscopic technique.     

The visualization of large organized chromatin domains enriched in the H3K9me2 mark within a single chromosome in a single cell

Despite considerable efforts, our understanding of the organization of higher order chromatin conformations in single cells and how these relate to chromatin marks remains poor. We have earlier invented the Chromatin In Situ Proximity (ChrISP) technique to determine proximities between chromatin fibers within a single chromosome. Here we used ChrISP to identify chromosome 11-specific hubs that are enriched in the H3K9me2 mark and that project toward the nuclear membrane in finger-like structures. Conversely, chromosome 11-specfic chromatin hubs, visualized by the presence of either H3K9me1 or H3K9me3 marks, are chromosome-wide and largely absent at the nuclear periphery. As the nuclear periphery-specific chromatin hubs were lost in the induced reduction of H3K9me2 levels, they likely represent Large Organization Chromatin in Lysine Methylation (LOCK) domains, previously identified by ChIP-seq analysis. Strikingly, the downregulation of the H3K9me2/3 marks also led to the chromosome-wide compaction of chromosome 11, suggesting a pleiotropic function of these features not recognized before. The ChrISP-mediated visualization of dynamic chromatin states in single cells thus provides an analysis of chromatin structures with a resolution far exceeding that of any other light microscopic technique.     

Long ncRNA expression associates with tissue-specific enhancers

Long non-coding RNAs (ncRNA) have recently been demonstrated to be expressed from a subset of enhancers and to be required for the distant regulation of gene expression. Several approaches to predict enhancers have been developed based on various chromatin marks and occupancy of enhancer-binding proteins. Despite the rapid advances in the field, no consensus how to define tissue specific enhancers yet exists. Here, we identify 2,695 long ncRNAs annotated by ENCODE (corresponding to 28% of all ENCODE annotated long ncRNAs) that overlap tissue-specific enhancers. We use a recently developed algorithm to predict tissue-specific enhancers, PreSTIGE, that is based on the H3K4me1 mark and tissue specific expression of mRNAs. The expression of the long ncRNAs overlapping enhancers is significantly higher when the enhancer is predicted as active in a specific cell line, suggesting a general interdependency of active enhancers and expression of long ncRNAs. This dependency is not identified using previous enhancer prediction algorithms that do not account for expression of their downstream targets. The predicted enhancers that overlap annotated long ncRNAs generally have a lower ratio of H3K4me1 to H3K4me3, suggesting that enhancers expressing long ncRNAs might be associated with specific epigenetic marks. In conclusion, we demonstrate the tissue-specific predictive power of PreSTIGE and provide evidence for thousands of long ncRNAs that are expressed from active tissue-specific enhancers, suggesting a particularly important functional relationship between long ncRNAs and enhancer activity in determining tissue-specific gene expression.     

Long ncRNA expression associates with tissue-specific enhancers

Long non-coding RNAs (ncRNA) have recently been demonstrated to be expressed from a subset of enhancers and to be required for the distant regulation of gene expression. Several approaches to predict enhancers have been developed based on various chromatin marks and occupancy of enhancer-binding proteins. Despite the rapid advances in the field, no consensus how to define tissue specific enhancers yet exists. Here, we identify 2,695 long ncRNAs annotated by ENCODE (corresponding to 28% of all ENCODE annotated long ncRNAs) that overlap tissue-specific enhancers. We use a recently developed algorithm to predict tissue-specific enhancers, PreSTIGE, that is based on the H3K4me1 mark and tissue specific expression of mRNAs. The expression of the long ncRNAs overlapping enhancers is significantly higher when the enhancer is predicted as active in a specific cell line, suggesting a general interdependency of active enhancers and expression of long ncRNAs. This dependency is not identified using previous enhancer prediction algorithms that do not account for expression of their downstream targets. The predicted enhancers that overlap annotated long ncRNAs generally have a lower ratio of H3K4me1 to H3K4me3, suggesting that enhancers expressing long ncRNAs might be associated with specific epigenetic marks. In conclusion, we demonstrate the tissue-specific predictive power of PreSTIGE and provide evidence for thousands of long ncRNAs that are expressed from active tissue-specific enhancers, suggesting a particularly important functional relationship between long ncRNAs and enhancer activity in determining tissue-specific gene expression.     

Epigenomic Landscape of Human Fetal Brain,Heart, and Liver

The epigenetic regulation of spatiotemporal gene expression is crucial for human development. Here, we present whole-genome chromatin immunoprecipitation followed by high throughput DNA sequencing (ChIP-seq) analyses of a wide variety of histone markers in the brain, heart, and liver of early human embryos shortly after their formation. We identified 40,181 active enhancers, with a large portion showing tissue-specific and developmental stage-specific patterns, pointing to their roles in controlling the ordered spatiotemporal expression of the developmental genes in early human embryos. Moreover, using sequential ChIP-seq, we showed that all three organs have hundreds to thousands of bivalent domains that are marked by both H3K4me3 and H3K27me3, probably to keep the progenitor cells in these organs ready for immediate differentiation into diverse cell types during subsequent developmental processes. Our work illustrates the potentially critical roles of tissue-specific and developmental stage-specific epigenomes in regulating the spatiotemporal expression of developmental genes during early human embryonic development.     

M33 condenses chromatin through nuclear body formation and methylation of both histone H3 lysine 9 and lysine 27

Heterochromatin, a type of condensed DNA in eukaryotic cells, has two main categories: Constitutive heterochromatin, which contains H3K9 methylation, and facultative heterochromatin, which contains H3K27 methylation. Methylated H3K9 and H3K27 serve as docking sites for chromodomain-containing proteins that compact chromatin. M33 (also known as CBX2) is a chromodomain-containing protein that binds H3K27me3 and compacts chromatin in vitro. However, whether M33 mediates chromatin compaction in cellulo remains unknown. Here we show that M33 compacts chromatin into DAPI-intense heterochromatin domains in cells. The formation of these heterochromatin domains requires H3K27me3, which recruits M33 to form nuclear bodies. G9a and SUV39H1 are sequentially recruited into M33 nuclear bodies to create H3K9 methylated chromatin in a process that is independent of HP1α. Finally, M33 decreases progerin-induced nuclear envelope disruption caused by loss of heterochromatin. Our findings demonstrate that M33 mediates the formation of condensed chromatin by forming nuclear bodies containing both H3K27me3 and H3K9me3. Our model of M33-dependent chromatin condensation suggests H3K27 methylation corroborates with H3K9 methylation during the formation of facultative heterochromatin and provides the theoretical basis for developing novel therapies to treat heterochromatin-related diseases.     

Dynamics and memory of heterochromatin in living cells

Posttranslational histone modifications are important for gene regulation, yet the mode of propagation and the contribution to heritable gene expression states remains controversial. To address these questions, we developed a chromatin in vivo assay (CiA) system employing chemically induced proximity to initiate and terminate chromatin modifications in living cells. We selectively recruited HP1α to induce H3K9me3-dependent gene silencing and describe the kinetics and extent of chromatin modifications at the Oct4 locus in fibroblasts and pluripotent cells. H3K9me3 propagated symmetrically and continuously at average rates of ~0.18 nucleosomes/hr to produce domains of up to 10 kb. After removal of the HP1α stimulus, heterochromatic domains were heritably transmitted, undiminished through multiple cell generations. Our data enabled quantitative modeling of reaction kinetics, which revealed that dynamic competition between histone marking and turnover, determines the boundaries and stability of H3K9me3 domains. This framework predicts the steady-state dynamics and spatial features of the majority of euchromatic H3K9me3 domains over the genome.     

Hdac3, Setdb1, and Kap1 mark H3K9me3/H3K14ac bivalent regions in young and aged liver

Post‐translational modifications of histone tails play a crucial role in gene regulation. Here, we performed chromatin profiling by quantitative targeted mass spectrometry to assess all possible modifications of the core histones. We identified a bivalent combination, a dually marked H3K9me3/H3K14ac modification in the liver, that is significantly decreased in old hepatocytes. Subsequent sequential ChIP‐Seq identified dually marked single nucleosome regions, with reduced number of sites and decreased signal in old livers, confirming mass spectrometry results. We detected H3K9me3 and H3K14ac bulk ChIP‐Seq signal in reChIP nucleosome regions, suggesting a correlation between H3K9me3/H3K14ac bulk bivalent genomic regions and dually marked single nucleosomes. Histone H3K9 deacetylase Hdac3, as well as H3K9 methyltransferase Setdb1, found in complex Kap1, occupied both bulk and single nucleosome bivalent regions in both young and old livers, correlating to presence of H3K9me3. Expression of genes associated with bivalent regions in young liver, including those regulating cholesterol secretion and triglyceride synthesis, is upregulated in old liver once the bivalency is lost. Hence, H3K9me3/H3K14ac dually marked regions define a poised inactive state that is resolved with loss of one or both of the chromatin marks, which subsequently leads to change in gene expression.     

Chromosomal loop anchorage of the kappa immunoglobulin gene occurs next to the enhancer in a region containing topoisomerase II sites

Introduction of torsional stress into active chromatin domains requires that linear DNA molecules be anchored in vivo to impede free rotation. While searching for these anchorage elements, we have localized a nuclear matrix association region (MAR) within the mouse immunoglobulin kappa gene that contains two topoisomerase II sites and is adjacent to the tissue-specific enhancer. The same matrix contact occurs when the kappa locus is in germ-line (inactive) or rear-ranged (transcribed) configurations. This constitutive anchorage site partitions the gene into V-J and C region chromatin domains. We demonstrate that at least 10,000 similar and evolutionarily conserved MAR binding sites exist in the nucleus. We propose that these sites, in association with topoisomerase II and possibly in conjunction with enhancers, play fundamental roles in the functional organization of chromatin loop domains.