Single-cell 3D genome structure reveals distinct human pluripotent states.

BACKGROUND: Pluripotent states of embryonic stem cells (ESCs) with distinct transcriptional profiles affect ESC differentiative capacity and therapeutic potential. Although single-cell RNA sequencing has revealed additional subpopulations and specific features of naive and primed human pluripotent stem cells (hPSCs), the underlying mechanisms that regulate their specific transcription and that control their pluripotent states remain elusive. RESULTS: By single-cell analysis of high-resolution, three-dimensional (3D) genomic structure, we herein demonstrate that remodeling of genomic structure is highly associated with the pluripotent states of human ESCs (hESCs). The naive pluripotent state is featured with specialized 3D genomic structures and clear chromatin compartmentalization that is distinct from the primed state. The naive pluripotent state is achieved by remodeling the active euchromatin compartment and reducing chromatin interactions at the nuclear center. This unique genomic organization is linked to enhanced chromatin accessibility on enhancers and elevated expression levels of naive pluripotent genes localized to this region. In contradistinction, the primed state exhibits intermingled genomic organization. Moreover, active euchromatin and primed pluripotent genes are distributed at the nuclear periphery, while repressive heterochromatin is densely concentrated at the nuclear center, reducing chromatin accessibility and the transcription of naive genes. CONCLUSIONS: Our data provide insights into the chromatin structure of ESCs in their naive and primed states, and we identify specific patterns of modifications in transcription and chromatin structure that might explain the genes that are differentially expressed between naive and primed hESCs. Thus, the inversion or relocation of heterochromatin to euchromatin via compartmentalization is related to the regulation of chromatin accessibility, thereby defining pluripotent states and cellular identity.

ISWI chromatin remodeling complexes recruit NSD2 and H3K36me2 in pericentromeric heterochromatin.

Histone H3 lysine36 dimethylation (H3K36me2) is generally distributed in the gene body and euchromatic intergenic regions. However, we found that H3K36me2 is enriched in pericentromeric heterochromatin in some mouse cell lines. We here revealed the mechanism of heterochromatin targeting of H3K36me2. Among several H3K36 methyltransferases, NSD2 was responsible for inducing heterochromatic H3K36me2. Depletion and overexpression analyses of NSD2-associating proteins revealed that NSD2 recruitment to heterochromatin was mediated through the imitation switch (ISWI) chromatin remodeling complexes, such as BAZ1B-SMARCA5 (WICH), which directly binds to AT-rich DNA via a BAZ1B domain-containing AT-hook-like motifs. The abundance and stoichiometry of NSD2, SMARCA5, and BAZ1B could determine the localization of H3K36me2 in different cell types. In mouse embryos, H3K36me2 heterochromatin localization was observed at the two- to four-cell stages, suggesting its physiological relevance.

Active transcription and epigenetic reactions synergistically regulate meso-scale genomic organization.

In interphase nuclei, chromatin forms dense domains of characteristic sizes, but the influence of transcription and histone modifications on domain size is not understood. We present a theoretical model exploring this relationship, considering chromatin-chromatin interactions, histone modifications, and chromatin extrusion. We predict that the size of heterochromatic domains is governed by a balance among the diffusive flux of methylated histones sustaining them and the acetylation reactions in the domains and the process of loop extrusion via supercoiling by RNAPII at their periphery, which contributes to size reduction. Super-resolution and nano-imaging of five distinct cell lines confirm the predictions indicating that the absence of transcription leads to larger heterochromatin domains. Furthermore, the model accurately reproduces the findings regarding how transcription-mediated supercoiling loss can mitigate the impacts of excessive cohesin loading. Our findings shed light on the role of transcription in genome organization, offering insights into chromatin dynamics and potential therapeutic targets.

H1 restricts euchromatin-associated methylation pathways from heterochromatic encroachment.

Silencing pathways prevent transposable element (TE) proliferation and help to maintain genome integrity through cell division. Silenced genomic regions can be classified as either euchromatic or heterochromatic, and are targeted by genetically separable epigenetic pathways. In plants, the RNA-directed DNA methylation (RdDM) pathway targets mostly euchromatic regions, while CMT DNA methyltransferases are mainly associated with heterochromatin. However, many epigenetic features - including DNA methylation patterning - are largely indistinguishable between these regions, so how the functional separation is maintained is unclear. The linker histone H1 is preferentially localized to heterochromatin and has been proposed to restrict RdDM from encroachment. To test this hypothesis, we followed RdDM genomic localization in an h1 mutant by performing ChIP-seq on the largest subunit, NRPE1, of the central RdDM polymerase, Pol V. Loss of H1 resulted in NRPE1 enrichment predominantly in heterochromatic TEs. Increased NRPE1 binding was associated with increased chromatin accessibility in h1, suggesting that H1 restricts NRPE1 occupancy by compacting chromatin. However, RdDM occupancy did not impact H1 localization, demonstrating that H1 hierarchically restricts RdDM positioning. H1 mutants experience major symmetric (CG and CHG) DNA methylation gains, and by generating an h1/nrpe1 double mutant, we demonstrate these gains are largely independent of RdDM. However, loss of NRPE1 occupancy from a subset of euchromatic regions in h1 corresponded to the loss of methylation in all sequence contexts, while at ectopically bound heterochromatic loci, NRPE1 deposition correlated with increased methylation specifically in the CHH context. Additionally, we found that H1 similarly restricts the occupancy of the methylation reader, SUVH1, and polycomb-mediated H3K27me3. Together, the results support a model whereby H1 helps maintain the exclusivity of heterochromatin by preventing encroachment from other competing pathways.

Downregulation of Histone H3.3 Induces p53-Dependent Cellular Senescence in Human Diploid Fibroblasts.

Cellular senescence is an irreversible growth arrest that acts as a barrier to cancer initiation and progression. Histone alteration is one of the major events during replicative senescence. However, little is known about the function of H3.3 in cellular senescence. Here we found that the downregulation of H3.3 induced growth suppression with senescence-like phenotypes such as senescence-associated heterochromatin foci (SAHF) and ß-galactosidase (SA-ß-gal) activity. Furthermore, H3.3 depletion induced senescence-like phenotypes with the p53/p21-depedent pathway. In addition, we identified miR-22-3p, tumor suppressive miRNA, as an upstream regulator of the H3F3B (H3 histone, family 3B) gene which is the histone variant H3.3 and replaces conventional H3 in active genes. Therefore, our results reveal for the first time the molecular mechanisms for cellular senescence which are regulated by H3.3 abundance. Taken together, our studies suggest that H3.3 exerts functional roles in regulating cellular senescence and is a promising target for cancer therapy.

The SUN-family protein Sad1 mediates heterochromatin spatial organization through interaction with histone H2A-H2B.

Heterochromatin is generally associated with the nuclear periphery, but how the spatial organization of heterochromatin is regulated to ensure epigenetic silencing remains unclear. Here we found that Sad1, an inner nuclear membrane SUN-family protein in fission yeast, interacts with histone H2A-H2B but not H3-H4. We solved the crystal structure of the histone binding motif (HBM) of Sad1 in complex with H2A-H2B, revealing the intimate contacts between Sad1HBM and H2A-H2B. Structure-based mutagenesis studies revealed that the H2A-H2B-binding activity of Sad1 is required for the dynamic distribution of Sad1 throughout the nuclear envelope (NE). The Sad1-H2A-H2B complex mediates tethering telomeres and the mating-type locus to the NE. This complex is also important for heterochromatin silencing. Mechanistically, H2A-H2B enhances the interaction between Sad1 and HDACs, including Clr3 and Sir2, to maintain epigenetic identity of heterochromatin. Interestingly, our results suggest that Sad1 exhibits the histone-enhanced liquid-liquid phase separation property, which helps recruit heterochromatin factors to the NE. Our results uncover an unexpected role of SUN-family proteins in heterochromatin regulation and suggest a nucleosome-independent role of H2A-H2B in regulating Sad1's functionality.

Chromatin phase separation and nuclear shape fluctuations are correlated in a polymer model of the nucleus.

Abnormal cell nuclear shapes are hallmarks of diseases, including progeria, muscular dystrophy, and many cancers. Experiments have shown that disruption of heterochromatin and increases in euchromatin lead to nuclear deformations, such as blebs and ruptures. However, the physical mechanisms through which chromatin governs nuclear shape are poorly understood. To investigate how heterochromatin and euchromatin might govern nuclear morphology, we studied chromatin microphase separation in a composite coarse-grained polymer and elastic shell simulation model. By varying chromatin density, heterochromatin composition, and heterochromatin-lamina interactions, we show how the chromatin phase organization may perturb nuclear shape. Increasing chromatin density stabilizes the lamina against large fluctuations. However, increasing heterochromatin levels or heterochromatin-lamina interactions enhances nuclear shape fluctuations by a "wetting"-like interaction. In contrast, fluctuations are insensitive to heterochromatin's internal structure. Our simulations suggest that peripheral heterochromatin accumulation could perturb nuclear morphology, while nuclear shape stabilization likely occurs through mechanisms other than chromatin microphase organization.

Phase separation and inheritance of repressive chromatin domains.

Polycomb-associated chromatin and pericentromeric heterochromatin form genomic domains important for the epigenetic regulation of gene expression. Both Polycomb complexes and heterochromatin factors rely on 'read and write' mechanisms, which, on their own, are not sufficient to explain the formation and the maintenance of these epigenetic domains. Microscopy has revealed that they form specific nuclear compartments separated from the rest of the genome. Recently, some subunits of these molecular machineries have been shown to undergo phase separation, both in vitro and in vivo, suggesting that phase separation might play important roles in the formation and the function of these two kinds of repressive chromatin. In this review, we will present the recent advances in the field of facultative and constitutive heterochromatin formation and maintenance through phase separation.

MeCP2 binds to methylated DNA independently of phase separation and heterochromatin organisation.

Correlative evidence has suggested that the methyl-CpG-binding protein MeCP2 contributes to the formation of heterochromatin condensates via liquid-liquid phase separation. This interpretation has been reinforced by the observation that heterochromatin, DNA methylation and MeCP2 co-localise within prominent foci in mouse cells. The findings presented here revise this view. MeCP2 localisation is independent of heterochromatin as MeCP2 foci persist even when heterochromatin organisation is disrupted. Additionally, MeCP2 foci fail to show hallmarks of phase separation in live cells. Importantly, we find that mouse cellular models are highly atypical as MeCP2 distribution is diffuse in most mammalian species, including humans. Notably, MeCP2 foci are absent in Mus spretus which is a mouse subspecies lacking methylated satellite DNA repeats. We conclude that MeCP2 has no intrinsic tendency to form condensates and its localisation is independent of heterochromatin. Instead, the distribution of MeCP2 in the nucleus is primarily determined by global DNA methylation patterns.

Suv39h-catalyzed H3K9me3 is critical for euchromatic genome organization and the maintenance of gene transcription.

H3K9me3-dependent heterochromatin is critical for the silencing of repeat-rich pericentromeric regions and also has key roles in repressing lineage-inappropriate protein-coding genes in differentiation and development. Here, we investigate the molecular consequences of heterochromatin loss in cells deficient in both SUV39H1 and SUV39H2 (Suv39DKO), the major mammalian histone methyltransferase enzymes that catalyze heterochromatic H3K9me3 deposition. We reveal a paradoxical repression of protein-coding genes in Suv39DKO cells, with these differentially expressed genes principally in euchromatic (Tn5-accessible, H3K4me3- and H3K27ac-marked) rather than heterochromatic (H3K9me3-marked) or polycomb (H3K27me3-marked) regions. Examination of the three-dimensional (3D) nucleome reveals that transcriptomic dysregulation occurs in euchromatic regions close to the nuclear periphery in 3D space. Moreover, this transcriptomic dysregulation is highly correlated with altered 3D genome organization in Suv39DKO cells. Together, our results suggest that the nuclear lamina-tethering of Suv39-dependent H3K9me3 domains provides an essential scaffold to support euchromatic genome organization and the maintenance of gene transcription for healthy cellular function.

Heterochromatic 3D genome organization is directed by HP1a- and H3K9-dependent and independent mechanisms.

Whether and how histone post-translational modifications and the proteins that bind them drive 3D genome organization remains unanswered. Here, we evaluate the contribution of H3K9-methylated constitutive heterochromatin to 3D genome organization in Drosophila tissues. We find that the predominant organizational feature of wild-type tissues is the segregation of euchromatic chromosome arms from heterochromatic pericentromeres. Reciprocal perturbation of HP1a⋅H3K9me binding, using a point mutation in the HP1a chromodomain or replacement of the replication-dependent histone H3 with H3K9R mutant histones, revealed that HP1a binding to methylated H3K9 in constitutive heterochromatin is required to limit contact frequency between pericentromeres and chromosome arms and regulate the distance between arm and pericentromeric regions. Surprisingly, the self-association of pericentromeric regions is largely preserved despite the loss of H3K9 methylation and HP1a occupancy. Thus, the HP1a⋅H3K9 interaction contributes to but does not solely drive the segregation of euchromatin and heterochromatin inside the nucleus.

Mathematical model for the role of multiple pericentromeric repeats on heterochromatin assembly.

Although the length and constituting sequences for pericentromeric repeats are highly variable across eukaryotes, the presence of multiple pericentromeric repeats is one of the conserved features of the eukaryotic chromosomes. Pericentromeric heterochromatin is often misregulated in human diseases, with the expansion of pericentromeric repeats in human solid cancers. In this article, we have developed a mathematical model of the RNAi-dependent methylation of H3K9 in the pericentromeric region of fission yeast. Our model, which takes copy number as an explicit parameter, predicts that the pericentromere is silenced only if there are many copies of repeats. It becomes bistable or desilenced if the copy number of repeats is reduced. This suggests that the copy number of pericentromeric repeats alone can determine the fate of heterochromatin silencing in fission yeast. Through sensitivity analysis, we identified parameters that favor bistability and desilencing. Stochastic simulation shows that faster cell division and noise favor the desilenced state. These results show the unexpected role of pericentromeric repeat copy number in gene silencing and provide a quantitative basis for how the copy number allows or protects repetitive and unique parts of the genome from heterochromatin silencing, respectively.

Differentiation shifts from a reversible to an irreversible heterochromatin state at the DM1 locus.

Epigenetic defects caused by hereditary or de novo mutations are implicated in various human diseases. It remains uncertain whether correcting the underlying mutation can reverse these defects in patient cells. Here we show by the analysis of myotonic dystrophy type 1 (DM1)-related locus that in mutant human embryonic stem cells (hESCs), DNA methylation and H3K9me3 enrichments are completely abolished by repeat excision (CTG2000 expansion), whereas in patient myoblasts (CTG2600 expansion), repeat deletion fails to do so. This distinction between undifferentiated and differentiated cells arises during cell differentiation, and can be reversed by reprogramming of gene-edited myoblasts. We demonstrate that abnormal methylation in DM1 is distinctively maintained in the undifferentiated state by the activity of the de novo DNMTs (DNMT3b in tandem with DNMT3a). Overall, the findings highlight a crucial difference in heterochromatin maintenance between undifferentiated (sequence-dependent) and differentiated (sequence-independent) cells, thus underscoring the role of differentiation as a locking mechanism for repressive epigenetic modifications at the DM1 locus.

Interplay between charge distribution and DNA in shaping HP1 paralog phase separation and localization.

The heterochromatin protein 1 (HP1) family is a crucial component of heterochromatin with diverse functions in gene regulation, cell cycle control, and cell differentiation. In humans, there are three paralogs, HP1α, HP1ß, and HP1γ, which exhibit remarkable similarities in their domain architecture and sequence properties. Nevertheless, these paralogs display distinct behaviors in liquid-liquid phase separation (LLPS), a process linked to heterochromatin formation. Here, we employ a coarse-grained simulation framework to uncover the sequence features responsible for the observed differences in LLPS. We highlight the significance of the net charge and charge patterning along the sequence in governing paralog LLPS propensities. We also show that both highly conserved folded and less-conserved disordered domains contribute to the observed differences. Furthermore, we explore the potential co-localization of different HP1 paralogs in multicomponent assemblies and the impact of DNA on this process. Importantly, our study reveals that DNA can significantly reshape the stability of a minimal condensate formed by HP1 paralogs due to competitive interactions of HP1α with HP1ß and HP1γ versus DNA. In conclusion, our work highlights the physicochemical nature of interactions that govern the distinct phase-separation behaviors of HP1 paralogs and provides a molecular framework for understanding their role in chromatin organization.

Two-way feedback between chromatin compaction and histone modification state explains Saccharomyces cerevisiae heterochromatin bistability.

Compact chromatin is closely linked with gene silencing in part by sterically masking access to promoters, inhibiting transcription factor binding and preventing polymerase from efficiently transcribing a gene. However, a broader hypothesis suggests that chromatin compaction can be both a cause and a consequence of the locus histone modification state, with a tight bidirectional interaction underpinning bistable transcriptional states. To rigorously test this hypothesis, we developed a mathematical model for the dynamics of the HMR locus in Saccharomyces cerevisiae, that incorporates activating histone modifications, silencing proteins, and a dynamic, acetylation-dependent, three-dimensional locus size. Chromatin compaction enhances silencer protein binding, which in turn feeds back to remove activating histone modifications, leading to further compaction. The bistable output of the model was in good agreement with prior quantitative data, including switching rates from expressed to silent states (and vice versa), and protein binding/histone modification levels within the locus. We then tested the model by predicting changes in switching rates as the genetic length of the locus was increased, which were then experimentally verified. Such bidirectional feedback between chromatin compaction and the histone modification state may be a widespread and important regulatory mechanism given the hallmarks of many heterochromatic regions: physical chromatin compaction and dimerizing (or multivalent) silencing proteins.

Lifelong persistence of nuclear RNAs in the mouse brain.

Genomic DNA that resides in the nuclei of mammalian neurons can be as old as the organism itself. The life span of nuclear RNAs, which are critical for proper chromatin architecture and transcription regulation, has not been determined in adult tissues. In this work, we identified and characterized nuclear RNAs that do not turn over for at least 2 years in a subset of postnatally born cells in the mouse brain. These long-lived RNAs were stably retained in nuclei in a neural cell type-specific manner and were required for the maintenance of heterochromatin. Thus, the life span of neural cells may depend on both the molecular longevity of DNA for the storage of genetic information and also the extreme stability of RNA for the functional organization of chromatin.

Heterochromatin Is Not the Only Place for satDNAs: The High Diversity of satDNAs in the Euchromatin of the Beetle Chrysolina americana (Coleoptera, Chrysomelidae).

The satellitome of the beetle Chrysolina americana Linneo, 1758 has been characterized through chromosomal analysis, genomic sequencing, and bioinformatics tools. C-banding reveals the presence of constitutive heterochromatin blocks enriched in A+T content, primarily located in pericentromeric regions. Furthermore, a comprehensive satellitome analysis unveils the extensive diversity of satellite DNA families within the genome of C. americana. Using fluorescence in situ hybridization techniques and the innovative CHRISMAPP approach, we precisely map the localization of satDNA families on assembled chromosomes, providing insights into their organization and distribution patterns. Among the 165 identified satDNA families, only three of them exhibit a remarkable amplification and accumulation, forming large blocks predominantly in pericentromeric regions. In contrast, the remaining, less abundant satDNA families are dispersed throughout euchromatic regions, challenging the traditional association of satDNA with heterochromatin. Overall, our findings underscore the complexity of repetitive DNA elements in the genome of C. americana and emphasize the need for further exploration to elucidate their functional significance and evolutionary implications.

Chromodomain mutation in S. pombe Kat5/Mst1 affects centromere dynamics and DNA repair.

KAT5 (S. pombe Mst1, human TIP60) is a MYST family histone acetyltransferase conserved from yeast to humans that is involved in multiple cellular activities. This family is characterized in part by containing a chromodomain, a motif associated with binding methylated histones. We show that a chromodomain mutation in the S. pombe Kat5, mst1-W66R, has defects in pericentromere silencing. mst1-W66R is sensitive to camptothecin (CPT) but only at an increased temperature of 36°C, although it is proficient for growth at this temperature. We also describe a de-silencing effect at the pericentromere by CPT that is independent of RNAi and methylation machinery. We also show that mst1-W66R disrupts recruitment of proteins to repair foci in response to camptothecin-induced DNA damage. Our data suggest a function of Mst1 chromodomain in centromere heterochromatin formation and a separate role in genome-wide damage repair in CPT.

Linker histone H1 regulates homeostasis of heterochromatin-associated cRNAs.

Chromatin-associated RNAs (cRNAs) are a poorly characterized fraction of cellular RNAs that co-purify with chromatin. Their full complexity and the mechanisms regulating their packaging and chromatin association remain poorly understood. Here, we address these questions in Drosophila. We find that cRNAs constitute a heterogeneous group of RNA species that is abundant in heterochromatic transcripts. We show that heterochromatic cRNAs interact with the heterogeneous nuclear ribonucleoproteins (hnRNP) hrp36/hrp48 and that depletion of linker histone dH1 impairs this interaction. dH1 depletion induces the accumulation of RNA::DNA hybrids (R-loops) in heterochromatin and, as a consequence, increases retention of heterochromatic cRNAs. These effects correlate with increased RNA polymerase II (RNAPII) occupancy at heterochromatin. Notably, impairing cRNA assembly by depletion of hrp36/hrp48 mimics heterochromatic R-loop accumulation induced by dH1 depletion. We also show that dH1 depletion alters nucleosome organization, increasing accessibility of heterochromatin. Altogether, these perturbations facilitate annealing of cRNAs to the DNA template, enhancing R-loop formation and cRNA retention at heterochromatin.

Accelerated plasma-cell differentiation in Bach2-deficient mouse B cells is caused by altered IRF4 functions.

Transcription factors BACH2 and IRF4 are both essential for antibody class-switch recombination (CSR) in activated B lymphocytes, while they oppositely regulate the differentiation of plasma cells (PCs). Here, we investigated how BACH2 and IRF4 interact during CSR and plasma-cell differentiation. We found that BACH2 organizes heterochromatin formation of target gene loci in mouse splenic B cells, including targets of IRF4 activation such as Aicda, an inducer of CSR, and Prdm1, a master plasma-cell regulator. Release of these gene loci from heterochromatin in response to B-cell receptor stimulation was coupled to AKT-mTOR pathway activation. In Bach2-deficient B cells, PC genes' activation depended on IRF4 protein accumulation, without an increase in Irf4 mRNA. Mechanistically, a PU.1-IRF4 heterodimer in activated B cells promoted BACH2 function by inducing gene expression of Bach2 and Pten, a negative regulator of AKT signaling. Elevated AKT activity in Bach2-deficient B cells resulted in IRF4 protein accumulation. Thus, BACH2 and IRF4 mutually modulate the activity of each other, and BACH2 inhibits PC differentiation by both the repression of PC genes and the restriction of IRF4 protein accumulation.