Large DNA Methylation Nadirs Anchor Chromatin Loops Maintaining Hematopoietic Stem Cell Identity

PRC2-AgeIndex as a universal biomarker of aging and rejuvenation

DNA methylation (DNAm) is one of the most reliable biomarkers of aging across mammalian tissues. While the age-dependent global loss of DNAm has been well characterized, DNAm gain is less characterized. Studies have demonstrated that CpGs which gain methylation with age are enriched in Polycomb Repressive Complex 2 (PRC2) targets. However, whole-genome examination of all PRC2 targets as well as determination of the pan-tissue or tissue-specific nature of these associations is lacking. Here, we show that low-methylated regions (LMRs) which are highly bound by PRC2 in embryonic stem cells (PRC2 LMRs) gain methylation with age in all examined somatic mitotic cells. We estimated that this epigenetic change represents around 90% of the age-dependent DNAm gain genome-wide. Therefore, we propose the "PRC2-AgeIndex," defined as the average DNAm in PRC2 LMRs, as a universal biomarker of cellular aging in somatic cells which can distinguish the effect of different anti-aging interventions.

Strong association between genomic 3D structure and CRISPR cleavage efficiency

CRISPR is a gene editing technology which enables precise in-vivo genome editing; but its potential is hampered by its relatively low specificity and sensitivity. Improving CRISPR's on-target and off-target effects requires a better understanding of its mechanism and determinants. Here we demonstrate, for the first time, the chromosomal 3D spatial structure's association with CRISPR's cleavage efficiency, and its predictive capabilities. We used high-resolution Hi-C data to estimate the 3D distance between different regions in the human genome and utilized these spatial properties to generate 3D-based features, characterizing each region's density. We evaluated these features based on empirical, in-vivo CRISPR efficiency data and compared them to 425 features used in state-of-the-art models. The 3D features ranked in the top 13% of the features, and significantly improved the predictive power of LASSO and xgboost models trained with these features. The features indicated that sites with lower spatial density demonstrated higher efficiency. Understanding how CRISPR is affected by the 3D DNA structure provides insight into CRISPR's mechanism in general and improves our ability to correctly predict CRISPR's cleavage as well as design sgRNAs for therapeutic and scientific use.

CUX1 regulates human hematopoietic stem cell chromatin accessibility via the BAF complex

CUX1 is a homeodomain-containing transcription factor that is essential for the development and differentiation of multiple tissues. CUX1 is recurrently mutated or deleted in cancer, particularly in myeloid malignancies. However, the mechanism by which CUX1 regulates gene expression and differentiation remains poorly understood, creating a barrier to understanding the tumor-suppressive functions of CUX1. Here, we demonstrate that CUX1 directs the BAF chromatin remodeling complex to DNA to increase chromatin accessibility in hematopoietic cells. CUX1 preferentially regulates lineage-specific enhancers, and CUX1 target genes are predictive of cell fate in vivo. These data indicate that CUX1 regulates hematopoietic lineage commitment and homeostasis via pioneer factor activity, and CUX1 deficiency disrupts these processes in stem and progenitor cells, facilitating transformation.

From compartments to loops: understanding the unique chromatin organization in neuronal cells

The three-dimensional organization of the genome plays a central role in the regulation of cellular functions, particularly in the human brain. This review explores the intricacies of chromatin organization, highlighting the distinct structural patterns observed between neuronal and non-neuronal brain cells. We integrate findings from recent studies to elucidate the characteristics of various levels of chromatin organization, from differential compartmentalization and topologically associating domains (TADs) to chromatin loop formation. By defining the unique chromatin landscapes of neuronal and non-neuronal brain cells, these distinct structures contribute to the regulation of gene expression specific to each cell type. In particular, we discuss potential functional implications of unique neuronal chromatin organization characteristics, such as weaker compartmentalization, neuron-specific TAD boundaries enriched with active histone marks, and an increased number of chromatin loops. Additionally, we explore the role of Polycomb group (PcG) proteins in shaping cell-type-specific chromatin patterns. This review further emphasizes the impact of variations in chromatin architecture between neuronal and non-neuronal cells on brain development and the onset of neurological disorders. It highlights the need for further research to elucidate the details of chromatin organization in the human brain in order to unravel the complexities of brain function and the genetic mechanisms underlying neurological disorders. This research will help bridge a significant gap in our comprehension of the interplay between chromatin structure and cell functions.

Tet-mediated DNA methylation dynamics affect chromosome organization

DNA Methylation is a significant epigenetic modification that can modulate chromosome states, but its role in orchestrating chromosome organization has not been well elucidated. Here we systematically assessed the effects of DNA Methylation on chromosome organization with a multi-omics strategy to capture DNA Methylation and high-order chromosome interaction simultaneously on mouse embryonic stem cells with DNA methylation dioxygenase Tet triple knock-out (Tet-TKO). Globally, upon Tet-TKO, we observed weakened compartmentalization, corresponding to decreased methylation differences between CpG island (CGI) rich and poor domains. Tet-TKO could also induce hypermethylation for the CTCF binding peaks in TAD boundaries and chromatin loop anchors. Accordingly, CTCF peak generally weakened upon Tet-TKO, which results in weakened TAD structure and depletion of long-range chromatin loops. Genes that lost enhancer-promoter looping upon Tet-TKO showed DNA hypermethylation in their gene bodies, which may compensate for the disruption of gene expression. We also observed distinct effects of Tet1 and Tet2 on chromatin organization and increased DNA methylation correlation on spatially interacted fragments upon Tet inactivation. Our work showed the broad effects of Tet inactivation and DNA methylation dynamics on chromosome organization.

Hematologic DNMT3A reduction and high-fat diet synergize to promote weight gain and tissue inflammation

During aging, blood cell production becomes dominated by a limited number of variant hematopoietic stem cell (HSC) clones. Differentiated progeny of variant HSCs are thought to mediate the detrimental effects of such clonal hematopoiesis on organismal health, but the mechanisms are poorly understood. While somatic mutations in DNA methyltransferase 3A (DNMT3A) frequently drive clonal dominance, the aging milieu also likely contributes. Here, we examined in mice the interaction between high-fat diet (HFD) and reduced DNMT3A in hematopoietic cells; strikingly, this combination led to weight gain. HFD amplified pro-inflammatory pathways and upregulated inflammation-associated genes in mutant cells along a pro-myeloid trajectory. Aberrant DNA methylation during myeloid differentiation and in response to HFD led to pro-inflammatory activation and maintenance of stemness genes. These findings suggest that reduced DNMT3A in hematopoietic cells contributes to weight gain, inflammation, and metabolic dysfunction, highlighting a role for DNMT3A loss in the development of metabolic disorders.

The impact of DNA methylation on CTCF-mediated 3D genome organization

Cytosine DNA methylation is a highly conserved epigenetic mark in eukaryotes. Although the role of DNA methylation at gene promoters and repetitive elements has been extensively studied, the function of DNA methylation in other genomic contexts remains less clear. In the nucleus of mammalian cells, the genome is spatially organized at different levels, and strongly influences myriad genomic processes. There are a number of factors that regulate the three-dimensional (3D) organization of the genome, with the CTCF insulator protein being among the most well-characterized. Pertinently, CTCF binding has been reported as being DNA methylation-sensitive in certain contexts, perhaps most notably in the process of genomic imprinting. Therefore, it stands to reason that DNA methylation may play a broader role in the regulation of chromatin architecture. Here we summarize the current understanding that is relevant to both the mammalian DNA methylation and chromatin architecture fields and attempt to assess the extent to which DNA methylation impacts the folding of the genome. The focus is in early embryonic development and cellular transitions when the epigenome is in flux, but we also describe insights from pathological contexts, such as cancer, in which the epigenome and 3D genome organization are misregulated.

A distant global control region is essential for normal expression of anterior HOXA genes during mouse and human craniofacial development

Craniofacial abnormalities account for approximately one third of birth defects. The regulatory programs that build the face require precisely controlled spatiotemporal gene expression, achieved through tissue-specific enhancers. Clusters of coactivated enhancers and their target genes, known as superenhancers, are important in determining cell identity but have been largely unexplored in development. In this study we identified superenhancer regions unique to human embryonic craniofacial tissue. To demonstrate the importance of such regions in craniofacial development and disease, we focused on an ~600 kb noncoding region located between NPVF and NFE2L3. We identified long range interactions with this region in both human and mouse embryonic craniofacial tissue with the anterior portion of the HOXA gene cluster. Mice lacking this superenhancer exhibit perinatal lethality, and present with highly penetrant skull defects and orofacial clefts phenocopying Hoxa2-/- mice. Moreover, we identified two cases of de novo copy number changes of the superenhancer in humans both with severe craniofacial abnormalities. This evidence suggests we have identified a critical noncoding locus control region that specifically regulates anterior HOXA genes and copy number changes are pathogenic in human patients.

H3K27me3 spreading organizes canonical PRC1 chromatin architecture to regulate developmental programs

Polycomb Repressive Complex 2 (PRC2)-mediated histone H3K27 tri-methylation (H3K27me3) recruits canonical PRC1 (cPRC1) to maintain heterochromatin. In early development, polycomb-regulated genes are connected through long-range 3D interactions which resolve upon differentiation. Here, we report that polycomb looping is controlled by H3K27me3 spreading and regulates target gene silencing and cell fate specification. Using glioma-derived H3 Lys-27-Met (H3K27M) mutations as tools to restrict H3K27me3 deposition, we show that H3K27me3 confinement concentrates the chromatin pool of cPRC1, resulting in heightened 3D interactions mirroring chromatin architecture of pluripotency, and stringent gene repression that maintains cells in progenitor states to facilitate tumor development. Conversely, H3K27me3 spread in pluripotent stem cells, following neural differentiation or loss of the H3K36 methyltransferase NSD1, dilutes cPRC1 concentration and dissolves polycomb loops. These results identify the regulatory principles and disease implications of polycomb looping and nominate histone modification-guided distribution of reader complexes as an important mechanism for nuclear compartment organization.

Highlights

The confinement of H3K27me3 at PRC2 nucleation sites without its spreading correlates with increased 3D chromatin interactions.The H3K27M oncohistone concentrates canonical PRC1 that anchors chromatin loop interactions in gliomas, silencing developmental programs.Stem and progenitor cells require factors promoting H3K27me3 confinement, including H3K36me2, to maintain cPRC1 loop architecture.The cPRC1-H3K27me3 interaction is a targetable driver of aberrant self-renewal in tumor cells.

Lineage specific 3D genome structure in the adult human brain and neurodevelopmental changes in the chromatin interactome

The human brain is a complex organ comprised of distinct cell types, and the contribution of the 3D genome to lineage specific gene expression remains poorly understood. To decipher cell type specific genome architecture, and characterize fine scale changes in the chromatin interactome across neural development, we compared the 3D genome of the human fetal cortical plate to that of neurons and glia isolated from the adult prefrontal cortex. We found that neurons have weaker genome compartmentalization compared to glia, but stronger TADs, which emerge during fetal development. Furthermore, relative to glia, the neuronal genome shifts more strongly towards repressive compartments. Neurons have differential TAD boundaries that are proximal to active promoters involved in neurodevelopmental processes. CRISPRi on CNTNAP2 in hIPSC-derived neurons reveals that transcriptional inactivation correlates with loss of insulation at the differential boundary. Finally, re-wiring of chromatin loops during neural development is associated with transcriptional and functional changes. Importantly, differential loops in the fetal cortex are associated with autism GWAS loci, suggesting a neuropsychiatric disease mechanism affecting the chromatin interactome. Furthermore, neural development involves gaining enhancer-promoter loops that upregulate genes that control synaptic activity. Altogether, our study provides multi-scale insights on the 3D genome in the human brain.

Domains Rearranged Methylase 2 maintains DNA methylation at large DNA hypomethylated shores and long-range chromatin interactions in rice

DNA methylation plays an important role in gene regulation and genomic stability. However, large DNA hypomethylated regions known as DNA methylation valleys (DMVs) or canyons have also been suggested to serve unique regulatory functions, largely unknown in rice (Oryza sativa). Here, we describe the DMVs in rice seedlings, which were highly enriched with developmental and transcription regulatory genes. Further detailed analysis indicated that grand DMVs (gDMVs) might be derived from nuclear integrants of organelle DNA (NORGs). Furthermore, Domains Rearranged Methylase 2 (OsDRM2) maintained DNA methylation at short DMV (sDMV) shores. Epigenetic maps indicated that sDMVs were marked with H3K4me3 and/or H3K27me3, although the loss of DNA methylation had a negligible effect on histone modification within these regions. In addition, we constructed H3K27me3-associated interaction maps for homozygous T-DNA insertion mutant of the gene (osdrm2) and wild type (WT). From a global perspective, most (90%) compartments were stable between osdrm2 and WT plants. At a high resolution, we observed a dramatic loss of long-range chromatin loops in osdrm2, which suffered an extensive loss of non-CG (CHG and CHH, H = A, T, or C) methylation. From another viewpoint, the loss of non-CG methylation at sDMV shores in osdrm2 could disrupt H3K27me3-mediated chromatin interaction networks. Overall, our results demonstrated that DMVs are a key genomic feature in rice and are precisely regulated by epigenetic modifications, including DNA methylation and histone modifications. OsDRM2 maintained DNA methylation at sDMV shores, while OsDRM2 deficiency strongly affected three-dimensional (3D) genome architectures.

Integrome signatures of lentiviral gene therapy for SCID-X1 patients

Lentiviral vector (LV)-based gene therapy holds promise for a broad range of diseases. Analyzing more than 280,000 vector integration sites (VISs) in 273 samples from 10 patients with X-linked severe combined immunodeficiency (SCID-X1), we discovered shared LV integrome signatures in 9 of 10 patients in relation to the genomics, epigenomics, and 3D structure of the human genome. VISs were enriched in the nuclear subcompartment A1 and integrated into super-enhancers close to nuclear pore complexes. These signatures were validated in T cells transduced with an LV encoding a CD19-specific chimeric antigen receptor. Intriguingly, the one patient whose VISs deviated from the identified integrome signatures had a distinct clinical course. Comparison of LV and gamma retrovirus integromes regarding their 3D genome signatures identified differences that might explain the lower risk of insertional mutagenesis in LV-based gene therapy. Our findings suggest that LV integrome signatures, shaped by common features such as genome organization, may affect the efficacy of LV-based cellular therapies.

TLR8 escapes X chromosome inactivation in human monocytes and CD4+ T cells

BACKGROUND

Human endosomal Toll-like receptors TLR7 and TLR8 recognize self and non-self RNA ligands, and are important mediators of innate immunity and autoimmune pathogenesis. TLR7 and TLR8 are, respectively, encoded by adjacent X-linked genes. We previously established that TLR7 evades X chromosome inactivation (XCI) in female immune cells. Whether TLR8 also evades XCI, however, has not yet been explored.

METHOD

In the current study, we used RNA fluorescence in situ hybridization (RNA FISH) to directly visualize, on a single-cell basis, primary transcripts of TLR7 and TLR8 relative to X chromosome territories in CD14+ monocytes and CD4+ T lymphocytes from women, Klinefelter syndrome (KS) men, and euploid men. To assign X chromosome territories in cells lacking robust expression of a XIST compartment, we designed probes specific for X-linked genes that do not escape XCI and therefore robustly label the active X chromosome. We also assessed whether XCI escape of TLR8 was associated with sexual dimorphism in TLR8 protein expression by western blot and flow cytometry.

RESULTS

Using RNA FISH, we show that TLR8, like TLR7, evades XCI in immune cells, and that cells harboring simultaneously TLR7 and TLR8 transcript foci are more frequent in women and KS men than in euploid men, resulting in a sevenfold difference in frequency. This transcriptional bias was again observable when comparing the single X of XY males with the active X of cells from females or KS males. Interestingly, TLR8 protein expression was significantly higher in female mononuclear blood cells, including all monocyte subsets, than in male cells.

CONCLUSIONS

TLR8, mirroring TLR7, escapes XCI in human monocytes and CD4+ T cells. Co-dependent transcription from the active X chromosome and escape from XCI could both contribute to higher TLR8 protein abundance in female cells, which may have implications for the response to viruses and bacteria, and the risk of developing inflammatory and autoimmune diseases.

Widespread contribution of transposable elements to the rewiring of mammalian 3D genomes

Transposable elements (TEs) are major contributors of genetic material in mammalian genomes. These often include binding sites for architectural proteins, including the multifarious master protein, CTCF, which shapes the 3D genome by creating loops, domains, compartment borders, and RNA-DNA interactions. These play a role in the compact packaging of DNA and have the potential to facilitate regulatory function. In this study, we explore the widespread contribution of TEs to mammalian 3D genomes by quantifying the extent to which they give rise to loops and domain border differences across various cell types and species using several 3D genome mapping technologies. We show that specific families and subfamilies of TEs have contributed to lineage-specific 3D chromatin structures across mammalian species. In many cases, these loops may facilitate sustained interaction between distant cis-regulatory elements and target genes, and domains may segregate chromatin state to impact gene expression in a lineage-specific manner. An experimental validation of our analytical findings using CRISPR-Cas9 to delete a candidate TE resulted in disruption of species-specific 3D chromatin structure. Taken together, we comprehensively quantify and selectively validate our finding that TEs contribute to shaping 3D genome organization and may, in some cases, impact gene regulation during the course of mammalian evolution.

The aberrant epigenome of DNMT3B-mutated ICF1 patient iPSCs is amenable to correction, with the exception of a subset of regions with H3K4me3- and/or CTCF-based epigenetic memory

Bi-allelic hypomorphic mutations in DNMT3B disrupt DNA methyltransferase activity and lead to immunodeficiency, centromeric instability, facial anomalies syndrome, type 1 (ICF1). Although several ICF1 phenotypes have been linked to abnormally hypomethylated repetitive regions, the unique genomic regions responsible for the remaining disease phenotypes remain largely uncharacterized. Here we explored two ICF1 patient-derived induced pluripotent stem cells (iPSCs) and their CRISPR-Cas9-corrected clones to determine whether DNMT3B correction can globally overcome DNA methylation defects and related changes in the epigenome. Hypomethylated regions throughout the genome are highly comparable between ICF1 iPSCs carrying different DNMT3B variants, and significantly overlap with those in ICF1 patient peripheral blood and lymphoblastoid cell lines. These regions include large CpG island domains, as well as promoters and enhancers of several lineage-specific genes, in particular immune-related, suggesting that they are premarked during early development. CRISPR-corrected ICF1 iPSCs reveal that the majority of phenotype-related hypomethylated regions reacquire normal DNA methylation levels following editing. However, at the most severely hypomethylated regions in ICF1 iPSCs, which also display the highest increases in H3K4me3 levels and/or abnormal CTCF binding, the epigenetic memory persists, and hypomethylation remains uncorrected. Overall, we demonstrate that restoring the catalytic activity of DNMT3B can reverse the majority of the aberrant ICF1 epigenome. However, a small fraction of the genome is resilient to this rescue, highlighting the challenge of reverting disease states that are due to genome-wide epigenetic perturbations. Uncovering the basis for the persistent epigenetic memory will promote the development of strategies to overcome this obstacle.

A genetic disorder reveals a hematopoietic stem cell regulatory network co-opted in leukemia

The molecular regulation of human hematopoietic stem cell (HSC) maintenance is therapeutically important, but limitations in experimental systems and interspecies variation have constrained our knowledge of this process. Here, we have studied a rare genetic disorder due to MECOM haploinsufficiency, characterized by an early-onset absence of HSCs in vivo. By generating a faithful model of this disorder in primary human HSCs and coupling functional studies with integrative single-cell genomic analyses, we uncover a key transcriptional network involving hundreds of genes that is required for HSC maintenance. Through our analyses, we nominate cooperating transcriptional regulators and identify how MECOM prevents the CTCF-dependent genome reorganization that occurs as HSCs differentiate. We show that this transcriptional network is co-opted in high-risk leukemias, thereby enabling these cancers to acquire stem cell properties. Collectively, we illuminate a regulatory network necessary for HSC self-renewal through the study of a rare experiment of nature.

Isocitrate dehydrogenase 1 mutation drives leukemogenesis by PDGFRA activation due to insulator disruption in acute myeloid leukemia (AML)

Acute myeloid leukemia (AML) is characterized by complex molecular alterations and driver mutations. Elderly patients show increased frequencies of IDH mutations with high chemoresistance and relapse rates despite recent therapeutic advances. Besides being associated with global promoter hypermethylation, IDH1 mutation facilitated changes in 3D DNA-conformation by CTCF-anchor methylation and upregulated oncogene expression in glioma, correlating with poor prognosis. Here, we investigated the role of IDH1 p.R132H mutation in altering 3D DNA-architecture and subsequent oncogene activation in AML. Using public RNA-Seq data, we identified upregulation of tyrosine kinase PDGFRA in IDH1-mutant patients, correlating with poor prognosis. DNA methylation analysis identified CpG hypermethylation within a CTCF-anchor upstream of PDGFRA in IDH1-mutant patients. Increased PDGFRA expression, PDGFRA-CTCF methylation and decreased CTCF binding were confirmed in AML CRISPR cells with heterozygous IDH1 p.R132H mutation and upon exogenous 2-HG treatment. IDH1-mutant cells showed higher sensitivity to tyrosine kinase inhibitor dasatinib, which was supported by reduced blast count in a patient with refractory IDH1-mutant AML after dasatinib treatment. Our data illustrate that IDH1 p.R132H mutation leads to CTCF hypermethylation, disrupting DNA-looping and insulation of PDGFRA, resulting in PDGFRA upregulation in IDH1-mutant AML. Treatment with dasatinib may offer a novel treatment strategy for IDH1-mutant AML.

Diverse silent chromatin states modulate genome compartmentalization and loop extrusion barriers

The relationships between chromosomal compartmentalization, chromatin state and function are poorly understood. Here by profiling long-range contact frequencies in HCT116 colon cancer cells, we distinguish three silent chromatin states, comprising two types of heterochromatin and a state enriched for H3K9me2 and H2A.Z that exhibits neutral three-dimensional interaction preferences and which, to our knowledge, has not previously been characterized. We find that heterochromatin marked by H3K9me3, HP1α and HP1β correlates with strong compartmentalization. We demonstrate that disruption of DNA methyltransferase activity greatly remodels genome compartmentalization whereby domains lose H3K9me3-HP1α/β binding and acquire the neutrally interacting state while retaining late replication timing. Furthermore, we show that H3K9me3-HP1α/β heterochromatin is permissive to loop extrusion by cohesin but refractory to CTCF binding. Together, our work reveals a dynamic structural and organizational diversity of the silent portion of the genome and establishes connections between the regulation of chromatin state and chromosome organization, including an interplay between DNA methylation, compartmentalization and loop extrusion.

Nuclear envelope, chromatin organizers, histones, and DNA: The many achilles heels exploited across cancers

In eukaryotic cells, the genome is organized in the form of chromatin composed of DNA and histones that organize and regulate gene expression. The dysregulation of chromatin remodeling, including the aberrant incorporation of histone variants and their consequent post-translational modifications, is prevalent across cancers. Additionally, nuclear envelope proteins are often deregulated in cancers, which impacts the 3D organization of the genome. Altered nuclear morphology, genome organization, and gene expression are defining features of cancers. With advances in single-cell sequencing, imaging technologies, and high-end data mining approaches, we are now at the forefront of designing appropriate small molecules to selectively inhibit the growth and proliferation of cancer cells in a genome- and epigenome-specific manner. Here, we review recent advances and the emerging significance of aberrations in nuclear envelope proteins, histone variants, and oncohistones in deregulating chromatin organization and gene expression in oncogenesis.

Subtype-specific 3D genome alteration in acute myeloid leukaemia

Acute myeloid leukaemia (AML) represents a set of heterogeneous myeloid malignancies, and hallmarks include mutations in epigenetic modifiers, transcription factors and kinases1-5. The extent to which mutations in AML drive alterations in chromatin 3D structure and contribute to myeloid transformation is unclear. Here we use Hi-C and whole-genome sequencing to analyse 25 samples from patients with AML and 7 samples from healthy donors. Recurrent and subtype-specific alterations in A/B compartments, topologically associating domains and chromatin loops were identified. RNA sequencing, ATAC with sequencing and CUT&Tag for CTCF, H3K27ac and H3K27me3 in the same AML samples also revealed extensive and recurrent AML-specific promoter-enhancer and promoter-silencer loops. We validated the role of repressive loops on their target genes by CRISPR deletion and interference. Structural variation-induced enhancer-hijacking and silencer-hijacking events were further identified in AML samples. Hijacked enhancers play a part in AML cell growth, as demonstrated by CRISPR screening, whereas hijacked silencers have a downregulating role, as evidenced by CRISPR-interference-mediated de-repression. Finally, whole-genome bisulfite sequencing of 20 AML and normal samples revealed the delicate relationship between DNA methylation, CTCF binding and 3D genome structure. Treatment of AML cells with a DNA hypomethylating agent and triple knockdown of DNMT1, DNMT3A and DNMT3B enabled the manipulation of DNA methylation to revert 3D genome organization and gene expression. Overall, this study provides a resource for leukaemia studies and highlights the role of repressive loops and hijacked cis elements in human diseases.

Single-cell multi-omics of human clonal hematopoiesis reveals that DNMT3A R882 mutations perturb early progenitor states through selective hypomethylation

Somatic mutations in cancer genes have been detected in clonal expansions across healthy human tissue, including in clonal hematopoiesis. However, because mutated and wild-type cells are admixed, we have limited ability to link genotypes with phenotypes. To overcome this limitation, we leveraged multi-modality single-cell sequencing, capturing genotype, transcriptomes and methylomes in progenitors from individuals with DNMT3A R882 mutated clonal hematopoiesis. DNMT3A mutations result in myeloid over lymphoid bias, and an expansion of immature myeloid progenitors primed toward megakaryocytic-erythroid fate, with dysregulated expression of lineage and leukemia stem cell markers. Mutated DNMT3A leads to preferential hypomethylation of polycomb repressive complex 2 targets and a specific CpG flanking motif. Notably, the hypomethylation motif is enriched in binding motifs of key hematopoietic transcription factors, serving as a potential mechanistic link between DNMT3A mutations and aberrant transcriptional phenotypes. Thus, single-cell multi-omics paves the road to defining the downstream consequences of mutations that drive clonal mosaicism.

Three-dimensional genome organization in immune cell fate and function

Immune cell development and activation demand the precise and coordinated control of transcriptional programmes. Three-dimensional (3D) organization of the genome has emerged as an important regulator of chromatin state, transcriptional activity and cell identity by facilitating or impeding long-range genomic interactions among regulatory elements and genes. Chromatin folding thus enables cell type-specific and stimulus-specific transcriptional responses to extracellular signals, which are essential for the control of immune cell fate, for inflammatory responses and for generating a diverse repertoire of antigen receptor specificities. Here, we review recent findings connecting 3D genome organization to the control of immune cell differentiation and function, and discuss how alterations in genome folding may lead to immune dysfunction and malignancy.

DYSF promotes monocyte activation in atherosclerotic cardiovascular disease as a DNA methylation-driven gene

Dysferlin (DYSF) has drawn much attention due to its involvement in dysferlinopathy and was reported to affect monocyte functions in recent studies. However, the role of DYSF in the pathogenesis of atherosclerotic cardiovascular diseases (ASCVD) and the regulation mechanism of DYSF expression have not been fully studied. In this study, Gene Expression Omnibus (GEO) database and epigenome-wide association study (EWAS) literatures were searched to find the DNA methylation-driven genes (including DYSF) of ASCVD. The hub genes related to DYSF were also identified through weighted correlation network analysis (WGCNA). Regulation of DYSF expression through its promoter methylation status was verified using peripheral blood leucocytes (PBLs) from ASCVD patients and normal controls, and experiments on THP1 cells and Apoe^-/- mice. Similarly, the expressions of DYSF related hub genes, mainly contained SELL, STAT3 and TMX1, were also validated. DYSF functions were then evaluated by phagocytosis, transwell and adhesion assays in DYSF knock-down and overexpressed THP1 cells. The results showed that DYSF promoter hypermethylation up-regulated its expression in clinical samples, THP1 cells and Apoe^-/- mice, confirming DYSF as a DNA methylation-driven gene. The combination of DYSF expression and methylation status in PBLs had a considerable prediction value for ASCVD. Besides, DYSF could enhance the phagocytosis, migration and adhesion ability of THP1 cells. Among DYSF related hub genes, SELL was proven to be the downstream target of DYSF by wet experiments. In conclusion, DYSF promoter hypermethylation upregulated its expression and promoted monocytes activation, which further participated in the pathogenesis of ASCVD.

DNA methylation in transposable elements buffers the connection between three-dimensional chromatin organization and gene transcription upon rice genome duplication

INTRODUCTION

Polyploidy is a major force in plant evolution and the domestication of cultivated crops.

OBJECTIVES

The study aimed to explore the relationship and underlying mechanism between three-dimensional (3D) chromatin organization and gene transcription upon rice genome duplication.

METHODS

The 3D chromatin structures between diploid (2C) and autotetraploid (4C) rice were compared using high-throughput chromosome conformation capture (Hi-C) analysis. The study combined genetics, transcriptomics, whole-genome bisulfite sequencing (WGBS-seq) and 3D genomics approaches to uncover the mechanism for DNA methylation in modulating gene transcription through 3D chromatin architectures upon rice genome duplication.

RESULTS

We found that 4C rice presents weakened intra-chromosomal interactions compared to its 2C progenitor in some chromosomes. In addition, we found that changes of 3D chromatin organizations including chromatin compartments, topologically associating domains (TADs), and loops, are uncorrelated with gene transcription. Moreover, DNA methylations in the regulatory sequences of genes in compartment A/B switched regions and TAD boundaries are unrelated to their expression. Importantly, although there was no significant difference in the methylation levels in transposable elements (TEs) in differentially expressed gene (DEG) and non-DEG promoters between 2C and 4C rice, we found that the hypermethylated TEs across genes in compartment A/B switched regions and TAD boundaries may suppress the expression of these genes.

CONCLUSION

The study proposed that the rice genome doubling might modulate TE methylation to buffer the effects of chromatin architecture on gene transcription in compartment A/B switched regions and TAD boundaries, resulting in the disconnection between 3D chromatin structure alteration and gene transcription upon rice genome duplication.

Polycomb-lamina antagonism partitions heterochromatin at the nuclear periphery

The genome can be divided into two spatially segregated compartments, A and B, which partition active and inactive chromatin states. While constitutive heterochromatin is predominantly located within the B compartment near the nuclear lamina, facultative heterochromatin marked by H3K27me3 spans both compartments. How epigenetic modifications, compartmentalization, and lamina association collectively maintain heterochromatin architecture remains unclear. Here we develop Lamina-Inducible Methylation and Hi-C (LIMe-Hi-C) to jointly measure chromosome conformation, DNA methylation, and lamina positioning. Through LIMe-Hi-C, we identify topologically distinct sub-compartments with high levels of H3K27me3 and differing degrees of lamina association. Inhibition of Polycomb repressive complex 2 (PRC2) reveals that H3K27me3 is essential for sub-compartment segregation. Unexpectedly, PRC2 inhibition promotes lamina association and constitutive heterochromatin spreading into H3K27me3-marked B sub-compartment regions. Consistent with this repositioning, genes originally marked with H3K27me3 in the B compartment, but not the A compartment, remain largely repressed, suggesting that constitutive heterochromatin spreading can compensate for H3K27me3 loss at a transcriptional level. These findings demonstrate that Polycomb sub-compartments and their antagonism with lamina association are fundamental features of genome structure. More broadly, by jointly measuring nuclear position and Hi-C contacts, our study demonstrates how compartmentalization and lamina association represent distinct but interdependent modes of heterochromatin regulation.

Epigenetic modifier SMCHD1 maintains a normal pool of long-term hematopoietic stem cells

SMCHD1 (structural maintenance of chromosomes hinge domain containing 1) is a noncanonical SMC protein that mediates long-range repressive chromatin structures. SMCHD1 is required for X chromosome inactivation in female cells and repression of imprinted and clustered autosomal genes, with SMCHD1 mutations linked to human diseases facioscapulohumeral muscular dystrophy (FSHD) and bosma arhinia and micropthalmia syndrome (BAMS). We used a conditional mouse model to investigate SMCHD1 in hematopoiesis. Smchd1-deleted mice maintained steady-state hematopoiesis despite showing an impaired reconstitution capacity in competitive bone marrow transplantations and age-related hematopoietic stem cell (HSC) loss. This phenotype was more pronounced in Smchd1-deleted females, which showed a loss of quiescent HSCs and fewer B cells. Gene expression profiling of Smchd1-deficient HSCs and B cells revealed known and cell-type-specific SMCHD1-sensitive genes and significant disruption to X-linked gene expression in female cells. These data show SMCHD1 is a regulator of HSCs whose effects are more profound in females.

•

SMCHD1 is not required to maintain steady-state hematopoiesis

•

Smchd1-deletion leads to loss of adult hematopoietic stem cells

•

Smchd1-deleted female mice are more severely affected than males

•

SMCHD1 maintains cellular quiescence in female hematopoietic stem cells

Hematopoietic differentiation persists in human iPSCs defective in de novo DNA methylation

BACKGROUND

DNA methylation is involved in the epigenetic regulation of gene expression during developmental processes and is primarily established by the DNA methyltransferase 3A (DNMT3A) and 3B (DNMT3B). DNMT3A is one of the most frequently mutated genes in clonal hematopoiesis and leukemia, indicating that it plays a crucial role for hematopoietic differentiation. However, the functional relevance of Dnmt3a for hematopoietic differentiation and hematological malignancies has mostly been analyzed in mice, with the specific role for human hematopoiesis remaining elusive. In this study, we therefore investigated if DNMT3A is essential for hematopoietic differentiation of human induced pluripotent stem cells (iPSCs).

RESULTS

We generated iPSC lines with knockout of either exon 2, 19, or 23 and analyzed the impact of different DNMT3A exon knockouts on directed differentiation toward mesenchymal and hematopoietic lineages. Exon 19-/- and 23-/- lines displayed an almost entire absence of de novo DNA methylation during mesenchymal and hematopoietic differentiation. Yet, differentiation efficiency was only slightly reduced in exon 19-/- and rather increased in exon 23-/- lines, while there was no significant impact on gene expression in hematopoietic progenitors (iHPCs). Notably, DNMT3A-/- iHPCs recapitulate some DNA methylation patterns of acute myeloid leukemia (AML) with DNMT3A mutations. Furthermore, multicolor genetic barcoding revealed growth advantage of exon 23-/- iHPCs in a syngeneic competitive differentiation assay.

CONCLUSIONS

Our results demonstrate that iPSCs with homozygous knockout of different exons of DNMT3A remain capable of mesenchymal and hematopoietic differentiation-and exon 23-/- iHPCs even gained growth advantage-despite loss of almost the entire de novo DNA methylation. Partial recapitulation of DNA methylation patterns of AML with DNMT3A mutations by our DNMT3A knockout iHPCs indicates that our model system can help to elucidate mechanisms of clonal hematopoiesis.

Super-Enhancers, Phase-Separated Condensates, and 3D Genome Organization in Cancer

3D chromatin organization plays an important role in transcription regulation and gene expression. The 3D genome is highly maintained by several architectural proteins, such as CTCF, Yin Yang 1, and cohesin complex. This structural organization brings regulatory DNA elements in close proximity to their target promoters. In this review, we discuss the 3D chromatin organization of super-enhancers and their relationship to phase-separated condensates. Super-enhancers are large clusters of DNA elements. They can physically contact with their target promoters by chromatin looping during transcription. Multiple transcription factors can bind to enhancer and promoter sequences and recruit a complex array of transcriptional co-activators and RNA polymerase II to effect transcriptional activation. Phase-separated condensates of transcription factors and transcriptional co-activators have been implicated in assembling the transcription machinery at particular enhancers. Cancer cells can hijack super-enhancers to drive oncogenic transcription to promote cell survival and proliferation. These dysregulated transcriptional programs can cause cancer cells to become highly dependent on transcriptional regulators, such as Mediator and BRD4. Moreover, the expression of oncogenes that are driven by super-enhancers is sensitive to transcriptional perturbation and often occurs in phase-separated condensates, supporting therapeutic rationales of targeting SE components, 3D genome organization, or dysregulated condensates in cancer.

Polycomb-mediated genome architecture enables long-range spreading of H3K27 methylation

Polycomb-group proteins play critical roles in gene silencing through the deposition of histone H3 lysine 27 trimethylation (H3K27me3) and chromatin compaction. This process is essential for embryonic stem cell (ESC) pluripotency, differentiation, and development. Polycomb repressive complex 2 (PRC2) can both read and write H3K27me3, enabling progressive spreading of H3K27me3 on the linear genome. Long-range Polycomb-associated DNA contacts have also been described, but their regulation and role in gene silencing remain unclear. Here, we apply H3K27me3 HiChIP, a protein-directed chromosome conformation method, and optical reconstruction of chromatin architecture to profile long-range Polycomb-associated DNA loops that span tens to hundreds of megabases across multiple topological associated domains in mouse ESCs and human induced pluripotent stem cells. We find that H3K27me3 loop anchors are enriched for Polycomb nucleation points and coincide with key developmental genes. Genetic deletion of H3K27me3 loop anchors results in disruption of spatial contact between distant loci and altered H3K27me3 in cis, both locally and megabases away on the same chromosome. In mouse embryos, loop anchor deletion leads to ectopic activation of the partner gene, suggesting that Polycomb-associated loops control gene silencing during development. Further, we find that alterations in PRC2 occupancy resulting from an RNA binding–deficient EZH2 mutant are accompanied by loss of Polycomb-associated DNA looping. Together, these results suggest PRC2 uses RNA binding to enhance long-range chromosome folding and H3K27me3 spreading. Developmental gene loci have unique roles in Polycomb spreading, emerging as important architectural elements of the epigenome.

Expression regulation of genes is linked to their CpG density distributions around transcription start sites

The CpG dinucleotide and its methylation behaviors play vital roles in gene regulation. Previous studies have divided genes into several categories based on the CpG intensity around transcription starting sites and found that housekeeping genes tend to possess high CpG density, whereas tissue-specific genes are generally characterized by low CpG density. In this study, we investigated how the CpG density distribution of a gene affects its transcription and regulation pattern. Based on the CpG density distribution around transcription starting site, by means of a semi-supervised neural network we designed, which took data augmentation into account, we divided the human genes into three categories, and genes within each cluster shared similar CpG density distribution. Not only sequence properties, these different clusters exhibited distinctly different structural features, regulatory mechanisms, correlation patterns between the expression level and CpG/TpG density, and expression and epigenetic mark variations during tumorigenesis. For instance, the activation of cluster 3 genes relies more on 3D genome reorganization, compared with cluster 1 and 2 genes, whereas cluster 2 genes showed the strongest correlation between gene expression and H3K27me3. Genes exhibiting uncoupled correlation between gene regulation and histone modifications are mainly in cluster 3. These results emphasized that the usage of epigenetic marks in gene regulation is partially rooted in the sequence property of genes such as their CpG density distribution and explained to some extent why the relation between epigenetic marks and gene expression is controversial.

Large-Scale Chromatin Rearrangements in Cancer

Simple Summary

Cancers have many genetic mutations such as nucleotide changes, deletions, amplifications, and chromosome gains or losses. Some of these genetic alterations directly contribute to the initiation and progression of tumors. In parallel to these genetic changes, cancer cells acquire modifications to their chromatin landscape, i.e., to the marks that are carried by DNA and the histone proteins it is associated with. These “epimutations” have consequences for gene expression and genome stability, and also contribute to tumoral initiation and progression. Some of these chromatin changes are very local, affecting just one or a few genes. In contrast, some chromatin alterations observed in cancer are more widespread and affect a large part of the genome. In this review, we present different types of large-scale chromatin rearrangements in cancer, explain how they may occur, and why they are relevant for cancer diagnosis and treatment.

Abstract

Epigenetic abnormalities are extremely widespread in cancer. Some of them are mere consequences of transformation, but some actively contribute to cancer initiation and progression; they provide powerful new biological markers, as well as new targets for therapies. In this review, we examine the recent literature and focus on one particular aspect of epigenome deregulation: large-scale chromatin changes, causing global changes of DNA methylation or histone modifications. After a brief overview of the one-dimension (1D) and three-dimension (3D) epigenome in healthy cells and of its homeostasis mechanisms, we use selected examples to describe how many different events (mutations, changes in metabolism, and infections) can cause profound changes to the epigenome and fuel cancer. We then present the consequences for therapies and briefly discuss the role of single-cell approaches for the future progress of the field.

The disordered N-terminal domain of DNMT3A recognizes H2AK119ub and is required for postnatal development

DNA methyltransferase 3a (DNMT3A) plays a crucial role during mammalian development. Two isoforms of DNMT3A are differentially expressed from stem cells to somatic tissues, but their individual functions remain largely uncharacterized. Here we report that the long isoform DNMT3A1, but not the short DNMT3A2, is essential for mouse postnatal development. DNMT3A1 binds to and regulates bivalent neurodevelopmental genes in the brain. Strikingly, Dnmt3a1 knockout perinatal lethality could be partially rescued by DNMT3A1 restoration in the nervous system. We further show that the intrinsically disordered N terminus of DNMT3A1 is required for normal development and DNA methylation at DNMT3A1-enriched regions. Mechanistically, a ubiquitin-interacting motif embedded in a putative α-helix within the N terminus binds to mono-ubiquitinated histone H2AK119, probably mediating recruitment of DNMT3A1 to Polycomb-regulated regions. These data demonstrate an isoform-specific role for DNMT3A1 in mouse postnatal development and reveal the N terminus as a necessary regulatory domain for DNMT3A1 chromatin occupancy and functions in the nervous system.

Anopheles mosquitoes reveal new principles of 3D genome organization in insects

Chromosomes are hierarchically folded within cell nuclei into territories, domains and subdomains, but the functional importance and evolutionary dynamics of these hierarchies are poorly defined. Here, we comprehensively profile genome organizations of five Anopheles mosquito species and show how different levels of chromatin architecture influence each other. Patterns observed on Hi-C maps are associated with known cytological structures, epigenetic profiles, and gene expression levels. Evolutionary analysis reveals conservation of chromatin architecture within synteny blocks for tens of millions of years and enrichment of synteny breakpoints in regions with increased genomic insulation. However, in-depth analysis shows a confounding effect of gene density on both insulation and distribution of synteny breakpoints, suggesting limited causal relationship between breakpoints and regions with increased genomic insulation. At the level of individual loci, we identify specific, extremely long-ranged looping interactions, conserved for ~100 million years. We demonstrate that the mechanisms underlying these looping contacts differ from previously described Polycomb-dependent interactions and clustering of active chromatin.

Lamin B1 deletion in myeloid neoplasms causes nuclear anomaly and altered hematopoietic stem cell function

Abnormal nuclear morphology is a hallmark of malignant cells widely used in cancer diagnosis. Pelger-Huët anomaly (PHA) is a common abnormality of neutrophil nuclear morphology of unknown molecular etiology in myeloid neoplasms (MNs). We show that loss of nuclear lamin B1 (LMNB1) encoded on chromosome 5q, which is frequently deleted in MNs, induces defects in nuclear morphology and human hematopoietic stem cell (HSC) function associated with malignancy. LMNB1 deficiency alters genome organization inducing in vitro and in vivo expansion of HSCs, myeloid-biased differentiation with impaired lymphoid commitment, and genome instability due to defective DNA damage repair. Nuclear dysmorphology of neutrophils in patients with MNs is associated with 5q deletions spanning the LMNB1 locus, and lamin B1 loss is both necessary and sufficient to cause PHA in normal and 5q-deleted neutrophils. LMNB1 loss thus causes acquired PHA and links abnormal nuclear morphology with HSCs and progenitor cell fate determination via genome organization.

Epigenome rewiring in human pluripotent stem cells

The epigenome plays a crucial role in modulating the activity of regulatory elements, thereby orchestrating diverse transcriptional programs during embryonic development. Human (h)PSC stepwise differentiation provides an excellent platform for capturing dynamic epigenomic events during lineage transition in human development. Here we discuss how recent technological advances, from epigenomic mapping to targeted perturbation, are providing a more comprehensive appreciation of remodeling of the chromatin landscape during human development with implications for aberrant rewiring in disease. We predict that the continuous innovation of hPSC differentiation methods, epigenome mapping, and CRISPR (clustered regularly interspaced short palindromic repeats) perturbation technologies will allow researchers to build toward not only a comprehensive understanding of the epigenomic mechanisms governing development, but also a highly flexible way to model diseases with opportunities for translation.

Focal disruption of DNA methylation dynamics at enhancers in IDH-mutant AML cells

Recurrent mutations in IDH1 or IDH2 in acute myeloid leukemia (AML) are associated with increased DNA methylation, but the genome-wide patterns of this hypermethylation phenotype have not been comprehensively studied in AML samples. We analyzed whole-genome bisulfite sequencing data from 15 primary AML samples with IDH1 or IDH2 mutations, which identified ~4000 focal regions that were uniquely hypermethylated in IDH^mut samples vs. normal CD34+ cells and other AMLs. These regions had modest hypermethylation in AMLs with biallelic TET2 mutations, and levels of 5-hydroxymethylation that were diminished in IDH and TET-mutant samples, indicating that this hypermethylation results from inhibition of TET-mediated demethylation. Focal hypermethylation in IDH^mut AMLs occurred at regions with low methylation in CD34+ cells, implying that DNA methylation and demethylation are active at these loci. AML samples containing IDH and DNMT3A^R882 mutations were significantly less hypermethylated, suggesting that IDH^mut-associated hypermethylation is mediated by DNMT3A. IDH^mut-specific hypermethylation was highly enriched for enhancers that form direct interactions with genes involved in normal hematopoiesis and AML, including MYC and ETV6. These results suggest that focal hypermethylation in IDH-mutant AML occurs by altering the balance between DNA methylation and demethylation, and that disruption of these pathways at enhancers may contribute to AML pathogenesis.

Redefining the biological basis of lineage-ambiguous leukemia through genomics: BCL11B deregulation in acute leukemias of ambiguous lineage

Acute leukemias of ambiguous lineage (ALAL), including mixed phenotype acute leukemia (MPAL) and related entities such as early T-cell precursor acute leukemia (ETP-ALL), remain diagnostic and clinical challenges due to limited understanding of pathogenesis, reliance of immunophenotyping to classify disease, and the lack of a rational approach to guide selection of appropriate therapy. Recent studies utilizing genomic sequencing and complementary approaches have provided key insights that are changing the way in which such leukemias are classified, and potentially, treated. Several recurrent genomic alterations define leukemias that straddle immunophenotypic entities, such as ZNF384-rearranged childhood B-ALL and B/myeloid MPAL, and BCL11B-rearranged T/myeloid MPAL, ETP-ALL and AML. In contrast, some cases of MPAL represent canonical ALL/AML entities exhibiting lineage aberrancy. For many cases of ALAL, experimental approaches indicate lineage aberrancy arises from acquisition of a founding genetic alteration into a hematopoietic stem or progenitor cell. Determination of optimal therapeutic approach requires genomic characterization of uniformly treated ALAL patients in prospective studies, but several approaches, including kinase inhibitors and BH3 mimetics may be efficacious in subsets of ALAL.

Enhancer Hijacking Drives Oncogenic BCL11B Expression in Lineage-Ambiguous Stem Cell Leukemia

Lineage-ambiguous leukemias are high-risk malignancies of poorly understood genetic basis. Here, we describe a distinct subgroup of acute leukemia with expression of myeloid, T lymphoid, and stem cell markers driven by aberrant allele-specific deregulation of BCL11B, a master transcription factor responsible for thymic T-lineage commitment and specification. Mechanistically, this deregulation was driven by chromosomal rearrangements that juxtapose BCL11B to superenhancers active in hematopoietic progenitors, or focal amplifications that generate a superenhancer from a noncoding element distal to BCL11B. Chromatin conformation analyses demonstrated long-range interactions of rearranged enhancers with the expressed BCL11B allele and association of BCL11B with activated hematopoietic progenitor cell cis-regulatory elements, suggesting BCL11B is aberrantly co-opted into a gene regulatory network that drives transformation by maintaining a progenitor state. These data support a role for ectopic BCL11B expression in primitive hematopoietic cells mediated by enhancer hijacking as an oncogenic driver of human lineage-ambiguous leukemia. SIGNIFICANCE: Lineage-ambiguous leukemias pose significant diagnostic and therapeutic challenges due to a poorly understood molecular and cellular basis. We identify oncogenic deregulation of BCL11B driven by diverse structural alterations, including de novo superenhancer generation, as the driving feature of a subset of lineage-ambiguous leukemias that transcend current diagnostic boundaries.This article is highlighted in the In This Issue feature, p. 2659.

Making Lemonade: Putting the Wisdom of the Genome to Work in Atopic Dermatitis

Getting from a GWAS hit to an actionable gene remains a challenge in complex disease genetics. In a new article of the Journal of Investigative Dermatology, Sobczyk et al. (2021) use a wide variety of genomic data to generate a prioritization algorithm to tackle this problem in atopic dermatitis, calling on the wisdom of the genome to generate promising results.

T-BET and EOMES Accelerate and Enhance Functional Differentiation of Human Natural Killer Cells

T-bet and Eomes are transcription factors that are known to be important in maturation and function of murine natural killer (NK) cells. Reduced T-BET and EOMES expression results in dysfunctional NK cells and failure to control tumor growth. In contrast to mice, the current knowledge on the role of T-BET and EOMES in human NK cells is rudimentary. Here, we ectopically expressed either T-BET or EOMES in human hematopoietic progenitor cells. Combined transcriptome, chromatin accessibility and protein expression analyses revealed that T-BET or EOMES epigenetically represses hematopoietic stem cell quiescence and non-NK lineage differentiation genes, while activating an NK cell-specific transcriptome and thereby drastically accelerating NK cell differentiation. In this model, the effects of T-BET and EOMES are largely overlapping, yet EOMES shows a superior role in early NK cell maturation and induces faster NK receptor and enhanced CD16 expression. T-BET particularly controls transcription of terminal maturation markers and epigenetically controls strong induction of KIR expression. Finally, NK cells generated upon T-BET or EOMES overexpression display improved functionality, including increased IFN-γ production and killing, and especially EOMES overexpression NK cells have enhanced antibody-dependent cellular cytotoxicity. Our findings reveal novel insights on the regulatory role of T-BET and EOMES in human NK cell maturation and function, which is essential to further understand human NK cell biology and to optimize adoptive NK cell therapies.

The concurrence of DNA methylation and demethylation is associated with transcription regulation

The mammalian DNA methylome is formed by two antagonizing processes, methylation by DNA methyltransferases (DNMT) and demethylation by ten-eleven translocation (TET) dioxygenases. Although the dynamics of either methylation or demethylation have been intensively studied in the past decade, the direct effects of their interaction on gene expression remain elusive. Here, we quantify the concurrence of DNA methylation and demethylation by the percentage of unmethylated CpGs within a partially methylated read from bisulfite sequencing. After verifying ‘methylation concurrence’ by its strong association with the co-localization of DNMT and TET enzymes, we observe that methylation concurrence is strongly correlated with gene expression. Notably, elevated methylation concurrence in tumors is associated with the repression of 40~60% of tumor suppressor genes, which cannot be explained by promoter hypermethylation alone. Furthermore, methylation concurrence can be used to stratify large undermethylated regions with negligible differences in average methylation into two subgroups with distinct chromatin accessibility and gene regulation patterns. Together, methylation concurrence represents a unique methylation metric important for transcription regulation and is distinct from conventional metrics, such as average methylation and methylation variation.

The global pattern of the mammalian methylome is formed by changes in methylation and demethylation. Here the authors describe a metric methylation concurrence that measures the ratio of unmethylated CpGs inside the partially methylated reads and show that methylation concurrence is associated with epigenetically regulated tumour suppressor genes.

Loss of grand histone H3 lysine 27 trimethylation domains mediated transcriptional activation in esophageal squamous cell carcinoma

Trimethylation of histone H3 lysine 27 trimethylation (H3K27me3) may be recruited by repressive Polycomb complexes to mediate gene silencing, which is critical for maintaining embryonic stem cell pluripotency and differentiation. However, the roles of aberrant H3K27me3 patterns in tumorigenesis are not fully understood. Here, we discovered that grand silencer domains (breadth > 50 kb) for H3K27me3 were significantly associated with epithelial cell differentiation and exhibited high gene essentiality and conservation in human esophageal epithelial cells. These grand H3K27me3 domains exhibited high modification signals involved in gene silencing, and preferentially occupied the entirety of topologically associating domains and interact with each other. We found that widespread loss of the grand H3K27me3 domains in of esophageal squamous cell carcinomas (ESCCs) were enriched in genes involved in epithelium and endothelium differentiation, which were significantly associated with overexpression with increase of active modifications of H3K4me3, H3K4me1, and H3K27ac marks, as well as DNA hypermethylation in the gene bodies. A total of 208 activated genes with loss of grand H3K27me3 domains in ESCC were identified, where the higher expression and mutation of T-box transcription factor 20 (TBX20) were associated with worse patients’ outcomes. Our results showed that knockdown of TBX20 may have led to a striking defect in esophageal cancer cell growth and carcinogenesis-related pathway, including cell cycle and homologous recombination. Together, our results reveal that loss of grand H3K27me3 domains represent a catalog of remarkable activating regulators involved in carcinogenesis.

STAG2 loss rewires oncogenic and developmental programs to promote metastasis in Ewing sarcoma

The core cohesin subunit STAG2 is recurrently mutated in Ewing sarcoma but its biological role is less clear. Here, we demonstrate that cohesin complexes containing STAG2 occupy enhancer and polycomb repressive complex (PRC2)-marked regulatory regions. Genetic suppression of STAG2 leads to a compensatory increase in cohesin-STAG1 complexes, but not in enhancer-rich regions, and results in reprogramming of cis-chromatin interactions. Strikingly, in STAG2 knockout cells the oncogenic genetic program driven by the fusion transcription factor EWS/FLI1 was highly perturbed, in part due to altered enhancer-promoter contacts. Moreover, loss of STAG2 also disrupted PRC2-mediated regulation of gene expression. Combined, these transcriptional changes converged to modulate EWS/FLI1, migratory, and neurodevelopmental programs. Finally, consistent with clinical observations, functional studies revealed that loss of STAG2 enhances the metastatic potential of Ewing sarcoma xenografts. Our findings demonstrate that STAG2 mutations can alter chromatin architecture and transcriptional programs to promote an aggressive cancer phenotype.

The interplay between DNA and histone methylation: molecular mechanisms and disease implications

Methylation of cytosine in CpG dinucleotides and histone lysine and arginine residues is a chromatin modification that critically contributes to the regulation of genome integrity, replication, and accessibility. A strong correlation exists between the genome-wide distribution of DNA and histone methylation, suggesting an intimate relationship between these epigenetic marks. Indeed, accumulating literature reveals complex mechanisms underlying the molecular crosstalk between DNA and histone methylation. These in vitro and in vivo discoveries are further supported by the finding that genes encoding DNA- and histone-modifying enzymes are often mutated in overlapping human diseases. Here, we summarize recent advances in understanding how DNA and histone methylation cooperate to maintain the cellular epigenomic landscape. We will also discuss the potential implication of these insights for understanding the etiology of, and developing biomarkers and therapies for, human congenital disorders and cancers that are driven by chromatin abnormalities.

The push and pull of DNA methylation

The control of DNA methylation guides cell-fate decisions in development and disease

Protein phase separation and its role in chromatin organization and diseases

In the physical sciences, solid, liquid, and gas are the most familiar phase states, whose essence is their existence reflecting the different spatial distribution of molecular components. The biological molecules in the living cell also have differences in spatial distribution. The molecules organized in the form of membrane-bound organelles are well recognized. However, the biomolecules organized in membraneless compartments called biomolecular condensates remain elusive. The liquid-liquid phase separation (LLPS), as a new emerging scientific breakthrough, describes the biomolecules assembled in special distribution and appeared as membraneless condensates in the form of a new "phase" compared with the surrounding liquid milieu. LLPS provides an important theoretical basis for explaining the composition of biological molecules and related biological reactions. Mounting evidence has emerged recently that phase-separated condensates participate in various biological activities. This article reviews the occurrence of LLPS and underlying regulatory mechanisms for understanding how multivalent molecules drive phase transitions to form the biomolecular condensates. And, it also summarizes recent major progress in elucidating the roles of LLPS in chromatin organization and provides clues for the development of new innovative therapeutic strategies for related diseases.

Clonal hematopoiesis: mechanisms driving dominance of stem cell clones

The discovery of clonal hematopoiesis (CH) in older individuals has changed the way hematologists and stem cell biologists view aging. Somatic mutations accumulate in stem cells over time. While most mutations have no impact, some result in subtle functional differences that ultimately manifest in distinct stem cell behaviors. With a large pool of stem cells and many decades to compete, some of these differences confer advantages under specific contexts. Approximately 20 genes are recurrently found as mutated in CH, indicating they confer some advantage. The impact of these mutations has begun to be analyzed at a molecular level by modeling in cell lines and in mice. Mutations in epigenetic regulators such as DNMT3A and TET2 confer an advantage by enhancing self-renewal of stem and progenitor cells and inhibiting their differentiation. Mutations in other genes involved in the DNA damage response may simply enhance cell survival. Here, we review proposed mechanisms that lead to CH, specifically in the context of stem cell biology, based on our current understanding of the function of some of the CH-associated genes.

BART3D: inferring transcriptional regulators associated with differential chromatin interactions from Hi-C data

SUMMARY

Identification of functional transcriptional regulators (TRs) associated with chromatin interactions is an important problem in studies of 3-dimensional genome organization and gene regulation. Direct inference of TR binding has been limited by the resolution of Hi-C data. Here, we present BART3D, a computational method for inferring TRs associated with genome-wide differential chromatin interactions by comparing Hi-C maps from two states, leveraging public ChIP-seq data for human and mouse. We demonstrate that BART3D can detect relevant TRs from dynamic Hi-C profiles with TR perturbation or cell differentiation. BART3D can be a useful tool in 3D genome data analysis and functional genomics research.

AVAILABILITY AND IMPLEMENTATION

BART3D is implemented in Python and the source code is available at https://github.com/zanglab/bart3d.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The acquisition of molecular drivers in pediatric therapy-related myeloid neoplasms

Pediatric therapy-related myeloid neoplasms (tMN) occur in children after exposure to cytotoxic therapy and have a dismal prognosis. The somatic and germline genomic alterations that drive these myeloid neoplasms in children and how they arise have yet to be comprehensively described. We use whole exome, whole genome, and/or RNA sequencing to characterize the genomic profile of 84 pediatric tMN cases (tMDS: n = 28, tAML: n = 56). Our data show that Ras/MAPK pathway mutations, alterations in RUNX1 or TP53, and KMT2A rearrangements are frequent somatic drivers, and we identify cases with aberrant MECOM expression secondary to enhancer hijacking. Unlike adults with tMN, we find no evidence of pre-existing minor tMN clones (including those with TP53 mutations), but rather the majority of cases are unrelated clones arising as a consequence of cytotoxic therapy. These studies also uncover rare cases of lineage switch disease rather than true secondary neoplasms.