Spatial clustering for identification of ChIP-enriched regions (SICER) to map regions of histone methylation patterns in embryonic stem cells

CTCF/RAD21 organize the ground state of chromatin-nuclear speckle association

Recent findings indicate that nuclear speckles, a distinct type of nuclear body, interact with certain chromatin regions in a ground state. Here, we report that the chromatin structural factors CTCF and cohesin are required for full ground-state association between DNA and nuclear speckles. We identified a putative speckle-targeting motif (STM) within cohesin subunit RAD21 and demonstrated that the STM is required for chromatin-nuclear speckle association, disruption of which also impaired induction of speckle-associated genes. Depletion of the cohesin-releasing factor WAPL, which stabilizes cohesin on chromatin, resulted in reinforcement of DNA-speckle contacts and enhanced inducibility of speckle-associated genes. Additionally, we observed disruption of chromatin-nuclear speckle association in patient-derived cells with Cornelia de Lange syndrome, a congenital neurodevelopmental disorder involving defective cohesin pathways. In summary, our findings reveal a mechanism for establishing the ground state of chromatin-speckle association and promoting gene inducibility, with relevance to human disease.

Copy number normalization distinguishes differential signals driven by copy number differences in ATAC-seq and ChIP-seq

A common objective across ATAC-seq and ChIP-seq analyses is to identify differential signals across contrasted conditions. However, in differential analyses, the impact of copy number variation is often overlooked. Here, we demonstrated copy number differences among samples could drive, if not dominate, differential signals. To address this, we propose a pipeline featuring copy number normalization. By comparing the averaged signal per gene copy, it effectively segregates differential signals driven by copy number from other factors. Further applying it to Down syndrome unveiled distinct dosage-dependent and -independent changes on chromosome 21. Thus, we recommend copy number normalization as a general approach.

VRN2-PRC2 facilitates light-triggered repression of PIF signaling to coordinate growth in Arabidopsis

VERNALIZATION2 (VRN2) is a flowering plant-specific subunit of the polycomb-repressive complex 2 (PRC2), a conserved eukaryotic holoenzyme that represses gene expression by depositing the histone H3 lysine 27 trimethylation (H3K27me3) mark in chromatin. Previous work established VRN2 as an oxygen-regulated target of the N-degron pathway that may function as a sensor subunit connecting PRC2 activity to the perception of endogenous and environmental cues. Here, we show that VRN2 is enriched in the hypoxic shoot apex and emerging leaves of Arabidopsis, where it negatively regulates growth by establishing a stable and conditionally repressed chromatin state in key PHYTOCHROME INTERACTING FACTOR (PIF)-regulated genes that promote cell expansion. This function is required to keep these genes poised for repression via a light-responsive signaling cascade later in leaf development. Thus, we identify VRN2-PRC2 as a core component of a developmentally and spatially encoded epigenetic mechanism that coordinates plant growth through facilitating the signal-dependent suppression of PIF signaling.

Hao-Fountain syndrome protein USP7 controls neuronal differentiation via BCOR-ncPRC1.1

Pathogenic variants in the ubiquitin-specific protease 7 (USP7) gene cause a neurodevelopmental disorder called Hao-Fountain syndrome. However, it remains unclear which of USP7's pleiotropic functions are relevant for neurodevelopment. Here, we present a combination of quantitative proteomics, transcriptomics, and epigenomics to define the USP7 regulatory circuitry during neuronal differentiation. USP7 activity is required for the transcriptional programs that direct both the differentiation of embryonic stem cells into neural stem cells and the neuronal differentiation of SH-SY5Y neuroblastoma cells. USP7 controls the dosage of the Polycomb monubiquitylated histone H2A lysine 119 (H2AK119ub1) ubiquitin ligase complexes ncPRC1.1 and ncPRC1.6. Loss-of-function experiments revealed that BCOR-ncPRC1.1, but not ncPRC1.6, is a key effector of USP7 during neuronal differentiation. Indeed, BCOR-ncPRC1.1 mediates a major portion of USP7-dependent gene regulation during this process. Besides providing a detailed map of the USP7 regulome during neurodifferentiation, our results suggest that USP7- and ncPRC1.1-associated neurodevelopmental disorders involve dysregulation of a shared epigenetic network.

Annotation of cis-regulatory-associated histone modifications in the genomes of two Thoroughbred stallions

The Functional Annotation of Animal Genomes (FAANG) consortium aims to annotate animal genomes across species, and work in the horse has substantially contributed to that goal. As part of this initiative, chromatin immunoprecipitation with sequencing (ChIP-seq) was performed to identify histone modifications corresponding to enhancers (H3K4me1), promoters (H3K4me3), activators (H3K27ac), and repressors (H3K27me3) in eight tissues from two Thoroughbred stallions: adipose, parietal cortex, heart, lamina, liver, lung, skeletal muscle, and testis. The average genome coverage of peaks identified by MACS2 for H3K4me1, H3K4me3, and H3K27ac was 6.2%, 2.2%, and 4.1%, respectively. Peaks were called for H3K27me3, a broad mark, using both MACS2 and SICERpy, with MACS2 identifying a greater average number of peaks (158K; 10.4% genome coverage) than SICERpy (32K; 24.3% genome coverage). Tissue-unique peaks were identified with BEDTools, and 1%-47% of peaks were unique to a tissue for a given histone modification. However, correlations among usable reads, total peak number, and unique peak number ranged from 0.01 to 0.92, indicating additional data collection is necessary to parse technical from true biological differences. These publicly available data expand a growing resource available for identifying regulatory regions within the equine genome, and they serve as a reference for genome regulation across healthy tissues of the adult Thoroughbred stallion.

Evaluating the Performance of Peak Calling Algorithms Available for Intracellular G-Quadruplex Sequencing

DNA G-quadruplexes (G4) are non-canonical DNA structures that play key roles in various biological processes. Antibody-dependent sequencing is an important tool for identifying intracellularly formed DNA G4s, and peak calling is a crucial step in processing the sequencing data. As the applicability of existing peak calling algorithms to intracellular G4 data has not been previously assessed, we systematically compared and evaluated these algorithms to determine those best suited for G4 detection. We selected seven representative candidates from 43 published peak calling algorithms for detailed evaluation. The performance of each candidate on six published intracellular G4 sequencing datasets (GSE107690, GSE145090, GSE133379, GSE178668ChIP-seq, GSE178668CUT&Tag, GSE221437) were assessed by precision and recall against customized benchmarks integrating results from multiple algorithms, as well as consistency with known G4 information (pG4 predicted by pqsfinder, oG4 from GSE63874, and multi-cell-line conserved G4s) and epigenetic signals. We identified MACS2, PeakRanger, and GoPeaks as the most effective algorithms for analyzing intracellular G4 sequencing data, and attributed their superior performance partially to the distribution model of sequencing reads/fragments used in the hypothesis testing step of the peak calling procedures. These findings provide guidance and rationale for selecting peak callers appropriate for intracellular G4 data.

Interrogation of the interplay between DNA N6-methyladenosine (6mA) and hypoxia-induced chromatin accessibility by a randomized empirical model (EnrichShuf)

N 6-Methyladenosine (6mA) is an epigenetic mark in eukaryotes regulating development, stress response and tumor progression. METTL4 has been reported as a 6mA methyltransferase induced by hypoxia. The detection and annotation of 6mA signals in mammalian cells have been hampered by the techniques and analytical methods developed so far. Here we developed a 6mA-ChIP-exo-5.1-seq to improve the sensitivity of detecting 6mAs in human cell lines. Furthermore, an EnrichShuf analysis tool for comprehensively comparing 6mA-ChIP-exo-5.1-seq, ATAC-seq, ChIP-seq and RNA-seq has been developed to annotate the functional relevance of 6mA in relation to chromatin accessibility and histone marks. Using a hypoxia-induced 6mA induction system as a model, we showed that hypoxic 6mA signals positively correlated with accessible chromatin regions. These 6mA signals correlate with their regulation by METTL4 under hypoxia, consistent with previous results. 6mAs also co-exist with H3K4me1, a histone mark for enhancers. Further analysis of enhancers using an ABC (active-by-contact) model shows that hypoxia-inducible factor-1α-induced H3K4me3 surrounds the 6mA/H3K4me1 site to augment active enhancers. These results suggest that correlation between 6mA and accessible chromatin regions plays a significant role in enhancer-promoter interactions during hypoxia-induced gene expression.

Ezh2 Delays Activation of Differentiation Genes During Normal Cerebellar Granule Neuron Development and in Medulloblastoma

Medulloblastoma (MB) is the most common malignant brain tumour in children. The Sonic Hedgehog (SHH)-medulloblastoma subtype arises from the cerebellar granule neuron lineage. Terminally differentiated neurons are incapable of undergoing further cell division, so an effective treatment for this tumour could be to force neuronal differentiation. Differentiation therapy provides a potential alternative for patients with medulloblastoma who harbor mutations that impair cell death pathways (TP53), which is associated a with high mortality. To this end, our goal was to explore epigenetic regulation of cerebellar granule neuron differentiation in medulloblastoma cells. Key regulators were discovered using chromatin immunoprecipitation with high-throughput sequencing. DNA-bound protein and chromatin protein modifications were investigated across all genes. We discovered that Ezh2-mediated tri-methylation of the H3 histone (H3K27me3), occurred on more than half of the 787 genes whose transcription normally increases as granule neurons terminally differentiate. Conditional knockout of Ezh2 led to early initiation of differentiation in granule neuron precursors (GNPs), but only after cell cycle exit had occurred. Similarly, in MB cells, neuronal differentiation could be induced by preventing H3K27me3 modifications using an Ezh2 inhibitor (UNC1999), but only when UNC1999 was combined with forced cell cycle exit driven by a CDK4/6 inhibitor (Palbociclib). Ezh2 emerges as a powerful restraint upon post-mitotic differentiation during normal GNP development and combination of Ezh2 inhibition with cell cycle exit leads to MB cell differentiation.

PML is a constitutive component of chromatin domains enriched in repetitive elements and duplicated gene clusters in cancer cells

Heterochromatin is a pivotal element in the functional organization of genomes. In our study, we delve into the heterochromatin pattern of association by the PML (promyelocytic leukemia) protein. By using PML chromatin immunoprecipitation and sequencing data and comparing computational methodologies to depict PML chromatin association, we describe PML-associated domains or PADs as large heterochromatic regions that exhibit similar genomic features across cancer cell lines. We show that PADs are specifically enriched in non-coding genes, duplicated gene clusters, and repetitive DNA elements. Moreover, we find enriched binding motifs of KZFPs, which are involved in orchestrating epigenetic repression at repetitive DNA elements. Hence, our findings suggest that PML conservatively associates to heterochromatic domains enriched in repetitive DNA elements and duplicated gene clusters in cancer. These findings contribute to a broader understanding of the complex regulatory framework of genome organization by heterochromatin in cancer.

Phosphorylation of the DNA damage repair factor 53BP1 by ATM kinase controls neurodevelopmental programs in cortical brain organoids

53BP1 is a well-established DNA damage repair factor that has recently emerged to critically regulate gene expression for tumor suppression and neural development. However, its precise function and regulatory mechanisms remain unclear. Here, we showed that phosphorylation of 53BP1 at serine 25 by ATM is required for neural progenitor cell proliferation and neuronal differentiation in cortical brain organoids. Dynamic phosphorylation of 53BP1-serine 25 controls 53BP1 target genes governing neuronal differentiation and function, cellular response to stress, and apoptosis. Mechanistically, ATM and RNF168 govern 53BP1's binding to gene loci to directly affect gene regulation, especially at genes for neuronal differentiation and maturation. 53BP1 serine 25 phosphorylation effectively impedes its binding to bivalent or H3K27me3-occupied promoters, especially at genes regulating H3K4 methylation, neuronal functions, and cell proliferation. Beyond 53BP1, ATM-dependent phosphorylation displays wide-ranging effects, regulating factors in neuronal differentiation, cytoskeleton, p53 regulation, as well as key signaling pathways such as ATM, BDNF, and WNT during cortical organoid differentiation. Together, our data suggest that the interplay between 53BP1 and ATM orchestrates essential genetic programs for cell morphogenesis, tissue organization, and developmental pathways crucial for human cortical development.

Foxg1 regulates translation of neocortical neuronal genes, including the main NMDA receptor subunit gene, Grin1

BACKGROUND

Mainly known as a transcription factor patterning the rostral brain and governing its histogenesis, FOXG1 has been also detected outside the nucleus; however, biological meaning of that has been only partially clarified.

RESULTS

Prompted by FOXG1 expression in cytoplasm of pallial neurons, we investigated its implication in translational control. We documented the impact of FOXG1 on ribosomal recruitment of Grin1-mRNA, encoding for the main subunit of NMDA receptor. Next, we showed that FOXG1 increases GRIN1 protein level by enhancing the translation of its mRNA, while not increasing its stability. Molecular mechanisms underlying this activity included FOXG1 interaction with EIF4E and, possibly, Grin1-mRNA. Besides, we found that, within murine neocortical cultures, de novo synthesis of GRIN1 undergoes a prominent and reversible, homeostatic regulation and FOXG1 is instrumental to that. Finally, by integrated analysis of multiple omic data, we inferred that FOXG1 is implicated in translational control of hundreds of neuronal genes, modulating ribosome engagement and progression. In a few selected cases, we experimentally verified such inference.

CONCLUSIONS

These findings point to FOXG1 as a key effector, potentially crucial to multi-scale temporal tuning of neocortical pyramid activity, an issue with profound physiological and neuropathological implications.

Association with TFIIIC limits MYCN localisation in hubs of active promoters and chromatin accumulation of non-phosphorylated RNA polymerase II

MYC family oncoproteins regulate the expression of a large number of genes and broadly stimulate elongation by RNA polymerase II (RNAPII). While the factors that control the chromatin association of MYC proteins are well understood, much less is known about how interacting proteins mediate MYC's effects on transcription. Here, we show that TFIIIC, an architectural protein complex that controls the three-dimensional chromatin organisation at its target sites, binds directly to the amino-terminal transcriptional regulatory domain of MYCN. Surprisingly, TFIIIC has no discernible role in MYCN-dependent gene expression and transcription elongation. Instead, MYCN and TFIIIC preferentially bind to promoters with paused RNAPII and globally limit the accumulation of non-phosphorylated RNAPII at promoters. Consistent with its ubiquitous role in transcription, MYCN broadly participates in hubs of active promoters. Depletion of TFIIIC further increases MYCN localisation to these hubs. This increase correlates with a failure of the nuclear exosome and BRCA1, both of which are involved in nascent RNA degradation, to localise to active promoters. Our data suggest that MYCN and TFIIIC exert an censoring function in early transcription that limits promoter accumulation of inactive RNAPII and facilitates promoter-proximal degradation of nascent RNA.

Ageing-Related Changes to H3K4me3, H3K27ac, and H3K27me3 in Purified Mouse Neurons

Neurons are central to lifelong learning and memory, but ageing disrupts their morphology and function, leading to cognitive decline. Although epigenetic mechanisms are known to play crucial roles in learning and memory, neuron-specific genome-wide epigenetic maps into old age remain scarce, often being limited to whole-brain homogenates and confounded by glial cells. Here, we mapped H3K4me3, H3K27ac, and H3K27me3 in mouse neurons across their lifespan. This revealed stable H3K4me3 and global losses of H3K27ac and H3K27me3 into old age. We observed patterns of synaptic function gene deactivation, regulated through the loss of the active mark H3K27ac, but not H3K4me3. Alongside this, embryonic development loci lost repressive H3K27me3 in old age. This suggests a loss of a highly refined neuronal cellular identity linked to global chromatin reconfiguration. Collectively, these findings indicate a key role for epigenetic regulation in neurons that is inextricably linked with ageing.

Single-housing-induced islet epigenomic changes are related to polymorphisms in diabetic KK mice

A lack of social relationships is increasingly recognized as a type 2 diabetes (T2D) risk. To investigate the underlying mechanism, we used male KK mice, an inbred strain with spontaneous diabetes. Given the association between living alone and T2D risk in humans, we divided the non-diabetic mice into singly housed (KK-SH) and group-housed control mice. Around the onset of diabetes in KK-SH mice, we compared H3K27ac ChIP-Seq with RNA-Seq using pancreatic islets derived from each experimental group, revealing a positive correlation between single-housing-induced changes in H3K27ac and gene expression levels. In particular, single-housing-induced H3K27ac decreases revealed a significant association with islet cell functions and GWAS loci for T2D and related diseases, with significant enrichment of binding motifs for transcription factors representative of human diabetes. Although these H3K27ac regions were preferentially localized to a polymorphic genomic background, SNVs and indels did not cause sequence disruption of enriched transcription factor motifs in most of these elements. These results suggest alternative roles of genetic variants in environment-dependent epigenomic changes and provide insights into the complex mode of disease inheritance.

Dysregulated expression of cholesterol biosynthetic genes in Alzheimer's disease alters epigenomic signatures of hippocampal neurons

Aging is the main risk factor of cognitive neurodegenerative diseases such as Alzheimer's disease, with epigenome alterations as a contributing factor. Here, we compared transcriptomic/epigenomic changes in the hippocampus, modified by aging and by tauopathy, an AD-related feature. We show that the cholesterol biosynthesis pathway is severely impaired in hippocampal neurons of tauopathic but not of aged mice pointing to vulnerability of these neurons in the disease. At the epigenomic level, histone hyperacetylation was observed at neuronal enhancers associated with glutamatergic regulations only in the tauopathy. Lastly, a treatment of tau mice with the CSP-TTK21 epi-drug that restored expression of key cholesterol biosynthesis genes counteracted hyperacetylation at neuronal enhancers and restored object memory. As acetyl-CoA is the primary substrate of both pathways, these data suggest that the rate of the cholesterol biosynthesis in hippocampal neurons may trigger epigenetic-driven changes, that may compromise the functions of hippocampal neurons in pathological conditions.

Merkel cell polyomavirus small tumor antigen contributes to immune evasion by interfering with type I interferon signaling

Merkel cell polyomavirus (MCPyV) is the causative agent of the majority of Merkel cell carcinomas (MCC). The virus has limited coding capacity, with its early viral proteins, large T (LT) and small T (sT), being multifunctional and contributing to infection and transformation. A fundamental difference in early viral gene expression between infection and MCPyV-driven tumorigenesis is the expression of a truncated LT (LTtr) in the tumor. In contrast, sT is expressed in both conditions and contributes significantly to oncogenesis. Here, we identified novel functions of early viral proteins by performing genome-wide transcriptome and chromatin studies in primary human fibroblasts. Due to current limitations in infection and tumorigenesis models, we mimic these conditions by ectopically expressing sT, LT or LTtr, individually or in combination, at different time points. In addition to its known function in cell cycle and inflammation modulation, we reveal a fundamentally new function of sT. We show that sT regulates the type I interferon (IFN) response downstream of the type I interferon receptor (IFNAR) by interfering with the interferon-stimulated gene factor 3 (ISGF3)-induced interferon-stimulated gene (ISG) response. Expression of sT leads to a reduction in the expression of interferon regulatory factor 9 (IRF9) which is a central component of the ISGF3 complex. We further show that this function of sT is conserved in BKPyV. We provide a first mechanistic understanding of which early viral proteins trigger and control the type I IFN response, which may influence MCPyV infection, persistence and, during MCC progression, regulation of the tumor microenvironment.

Endonucleosis mediates internalization of cytoplasm into the nucleus

Setd8 regulates transcription elongation, mitotic DNA condensation, DNA damage response and replication licensing. Here we show that, in mitogen-stimulated liver-specific Setd8-KO mice, most of the hepatocytes are eliminated by necrosis but a significant number of them survive via entering a stage exhibiting several senescence-related features. Setd8-deficient hepatocytes had enlarged nuclei, chromosomal hyperploidy and nuclear engulfments progressing to the formation of intranuclear vesicles surrounded by nuclear lamina. These vesicles contain glycogen, cytoplasmic proteins and even entire organelles. We term this process "endonucleosis". Intranuclear vesicles are absent in hepatocytes of Setd8/Atg5 knockout mice, suggesting that the process requires the function of the canonical autophagy machinery. Endonucleosis and hyperploidization are temporary, early events in the surviving Setd8-deficient cells. Larger vesicles break down into microvesicles over time and are eventually eliminated. The results reveal sequential events in cells with extensive DNA damage, which function as part of survival mechanisms to prevent necrotic death.

Reinforcement of repressive marks in the chicken primordial germ cell epigenetic signature: divergence from basal state resetting in mammals

BACKGROUND

In mammals, primordial germ cells (PGCs), the embryonic precursors of the germline, arise from embryonic or extra-embryonic cells upon induction by the surrounding tissues during gastrulation, according to mechanisms which are elucidated in mice but remain controversial in primates. They undergo genome-wide epigenetic reprogramming, consisting of extensive DNA demethylation and histone post-translational modification (PTM) changes, toward a basal, euchromatinized state. In contrast, chicken PGCs are specified by preformation before gastrulation based on maternally-inherited factors. They can be isolated from the bloodstream during their migration to the genital ridges. Our prior research highlighted differences in the global epigenetic profile of cultured chicken PGCs compared with chicken somatic cells and mammalian PGCs. This study investigates the acquisition and evolution of this profile during development.

RESULTS

Quantitative analysis of global DNA methylation and histone PTMs, including their distribution, during key stages of chicken early development revealed divergent PGC epigenetic changes compared with mammals. Unlike mammalian PGCs, chicken PGCs do not undergo genome-wide DNA demethylation or exhibit a decrease in histone H3 lysine 9 dimethylation. However, chicken PGCs show 5‑hydroxymethylcytosine loss, macroH2A redistribution, and chromatin decompaction, mirroring mammalian processes. Chicken PGCs initiate their epigenetic signature during migration, progressively accumulating high global levels of H3K9me3, with preferential enrichment in inactive genome regions. Despite apparent global chromatin decompaction, abundant heterochromatin marks, including repressive histone PTMs, HP1 variants, and DNA methylation, persists in chicken PGCs, contrasting with mammalian PGCs.

CONCLUSIONS

Chicken PGCs' epigenetic signature does not align with the basal chromatin state observed in mammals, suggesting a departure from extensive epigenetic reprogramming. Despite disparities in early PGC development, the persistence of several epigenetic features shared with mammals implies their involvement in chromatin-regulated germ cell properties, with the distinctive elevation of chicken-specific H3K9me3 potentially participating in these processes.

Genome-wide profiling of histone (H3) lysine 4 (K4) tri-methylation (me3) under drought, heat, and combined stresses in switchgrass

BACKGROUND

Switchgrass (Panicum virgatum L.) is a warm-season perennial (C4) grass identified as an important biofuel crop in the United States. It is well adapted to the marginal environment where heat and moisture stresses predominantly affect crop growth. However, the underlying molecular mechanisms associated with heat and drought stress tolerance still need to be fully understood in switchgrass. The methylation of H3K4 is often associated with transcriptional activation of genes, including stress-responsive. Therefore, this study aimed to analyze genome-wide histone H3K4-tri-methylation in switchgrass under heat, drought, and combined stress.

RESULTS

In total, ~ 1.3 million H3K4me3 peaks were identified in this study using SICER. Among them, 7,342; 6,510; and 8,536 peaks responded under drought (DT), drought and heat (DTHT), and heat (HT) stresses, respectively. Most DT and DTHT peaks spanned 0 to + 2000 bases from the transcription start site [TSS]. By comparing differentially marked peaks with RNA-Seq data, we identified peaks associated with genes: 155 DT-responsive peaks with 118 DT-responsive genes, 121 DTHT-responsive peaks with 110 DTHT-responsive genes, and 175 HT-responsive peaks with 136 HT-responsive genes. We have identified various transcription factors involved in DT, DTHT, and HT stresses. Gene Ontology analysis using the AgriGO revealed that most genes belonged to biological processes. Most annotated peaks belonged to metabolite interconversion, RNA metabolism, transporter, protein modifying, defense/immunity, membrane traffic protein, transmembrane signal receptor, and transcriptional regulator protein families. Further, we identified significant peaks associated with TFs, hormones, signaling, fatty acid and carbohydrate metabolism, and secondary metabolites. qRT-PCR analysis revealed the relative expressions of six abiotic stress-responsive genes (transketolase, chromatin remodeling factor-CDH3, fatty-acid desaturase A, transmembrane protein 14C, beta-amylase 1, and integrase-type DNA binding protein genes) that were significantly (P < 0.05) marked during drought, heat, and combined stresses by comparing stress-induced against un-stressed and input controls.

CONCLUSION

Our study provides a comprehensive and reproducible epigenomic analysis of drought, heat, and combined stress responses in switchgrass. Significant enrichment of H3K4me3 peaks downstream of the TSS of protein-coding genes was observed. In addition, the cost-effective experimental design, modified ChIP-Seq approach, and analyses presented here can serve as a prototype for other non-model plant species for conducting stress studies.

Genomic profiling of six human somatic histone H1 variants denotes that H1X accumulates at recently incorporated transposable elements

Histone H1, a vital component in chromatin structure, binds to linker DNA and regulates nuclear processes. We have investigated the distribution of histone H1 variants in a breast cancer cell line using ChIP-Seq. Two major groups of variants are identified: H1.2, H1.3, H1.5 and H1.0 are abundant in low GC regions (B compartment), while H1.4 and H1X preferentially localize in high GC regions (A compartment). Examining their abundance within transposable elements (TEs) reveals that H1X and H1.4 are enriched in recently-incorporated TEs (SVA and SINE-Alu), while H1.0/H1.2/H1.3/H1.5 are more abundant in older elements. Notably, H1X is particularly enriched in SVA families, while H1.4 shows the highest abundance in young AluY elements. Although low GC variants are generally enriched in LINE, LTR and DNA repeats, H1X and H1.4 are also abundant in a subset of recent LINE-L1 and LTR repeats. H1X enrichment at SVA and Alu is consistent across multiple cell lines. Further, H1X depletion leads to TE derepression, suggesting its role in maintaining TE repression. Overall, this study provides novel insights into the differential distribution of histone H1 variants among repetitive elements, highlighting the potential involvement of H1X in repressing TEs recently incorporated within the human genome.

Acute myeloid leukemias with UBTF tandem duplications are sensitive to menin inhibitors

ABSTRACT

UBTF tandem duplications (UBTF-TDs) have recently emerged as a recurrent alteration in pediatric and adult acute myeloid leukemia (AML). UBTF-TD leukemias are characterized by a poor response to conventional chemotherapy and a transcriptional signature that mirrors NUP98-rearranged and NPM1-mutant AMLs, including HOX-gene dysregulation. However, the mechanism by which UBTF-TD drives leukemogenesis remains unknown. In this study, we investigated the genomic occupancy of UBTF-TD in transformed cord blood CD34+ cells and patient-derived xenograft models. We found that UBTF-TD protein maintained genomic occupancy at ribosomal DNA loci while also occupying genomic targets commonly dysregulated in UBTF-TD myeloid malignancies, such as the HOXA/HOXB gene clusters and MEIS1. These data suggest that UBTF-TD is a gain-of-function alteration that results in mislocalization to genomic loci dysregulated in UBTF-TD leukemias. UBTF-TD also co-occupies key genomic loci with KMT2A and menin, which are known to be key partners involved in HOX-dysregulated leukemias. Using a protein degradation system, we showed that stemness, proliferation, and transcriptional signatures are dependent on sustained UBTF-TD localization to chromatin. Finally, we demonstrate that primary cells from UBTF-TD leukemias are sensitive to the menin inhibitor SNDX-5613, resulting in markedly reduced in vitro and in vivo tumor growth, myeloid differentiation, and abrogation of the UBTF-TD leukemic expression signature. These findings provide a viable therapeutic strategy for patients with this high-risk AML subtype.

CDCA7-associated global aberrant DNA hypomethylation translates to localized, tissue-specific transcriptional responses

Disruption of cell division cycle associated 7 (CDCA7) has been linked to aberrant DNA hypomethylation, but the impact of DNA methylation loss on transcription has not been investigated. Here, we show that CDCA7 is critical for maintaining global DNA methylation levels across multiple tissues in vivo. A pathogenic Cdca7 missense variant leads to the formation of large, aberrantly hypomethylated domains overlapping with the B genomic compartment but without affecting the deposition of H3K9 trimethylation (H3K9me3). CDCA7-associated aberrant DNA hypomethylation translated to localized, tissue-specific transcriptional dysregulation that affected large gene clusters. In the brain, we identify CDCA7 as a transcriptional repressor and epigenetic regulator of clustered protocadherin isoform choice. Increased protocadherin isoform expression frequency is accompanied by DNA methylation loss, gain of H3K4 trimethylation (H3K4me3), and increased binding of the transcriptional regulator CCCTC-binding factor (CTCF). Overall, our in vivo work identifies a key role for CDCA7 in safeguarding tissue-specific expression of gene clusters via the DNA methylation pathway.

Loss of NAT10 disrupts enhancer organization via p300 mislocalization and suppresses transcription of genes necessary for metastasis progression

Acetylation of protein and RNA represent a critical event for development and cancer progression. NAT10 is the only known RNA acetylase that catalyzes the N4-actylcytidine (ac4C) modification of RNAs. Here, we show that the loss of NAT10 significantly decreases lung metastasis in allograft and genetically engineered mouse models of breast cancer. NAT10 interacts with a mechanosensitive, metastasis susceptibility protein complex at the nuclear pore. In addition to its canonical role in RNA acetylation, we find that NAT10 interacts with p300 at gene enhancers. NAT10 loss is associated with p300 mislocalization into heterochromatin regions. NAT10 depletion disrupts enhancer organization, leading to alteration of gene transcription necessary for metastatic progression, including reduced myeloid cell-recruiting chemokines that results in a less metastasis-prone tumor microenvironment. Our study uncovers a distinct role of NAT10 in enhancer organization of metastatic tumor cells and suggests its involvement in the tumor-immune crosstalk dictating metastatic outcomes.

Disruption of H3K36 methylation provokes cellular plasticity to drive aberrant glandular formation and squamous carcinogenesis

Chromatin organization is essential for maintaining cell-fate trajectories and developmental programs. Here, we find that disruption of H3K36 methylation dramatically impairs normal epithelial differentiation and development, which promotes increased cellular plasticity and enrichment of alternative cell fates. Specifically, we observe a striking increase in the aberrant generation of excessive epithelial glandular tissues, including hypertrophic salivary, sebaceous, and meibomian glands, as well as enhanced squamous tumorigenesis. These phenotypic and gene expression manifestations are associated with loss of H3K36me2 and rewiring of repressive H3K27me3, changes we also observe in human patients with glandular hyperplasia. Collectively, these results have identified a critical role for H3K36 methylation in both in vivo epithelial cell-fate decisions and the prevention of squamous carcinogenesis and suggest that H3K36 methylation modulation may offer new avenues for the treatment of numerous common disorders driven by altered glandular function, which collectively affect large segments of the human population.

Differential Analysis of Protein-DNA Binding Using ChIP-Seq Data

Chromatin immunoprecipitation in combination with next-generation sequencing (ChIP-Seq) allows probing of protein-DNA binding in a rapid and genome-wide fashion. Herein we describe the required steps to preprocess ChIP-Seq data and to analyze the differential binding of proteins to DNA for perturbation experiments. In these experiments, different conditions are compared to find the underlying biological mechanisms caused by the stimulus or treatment. In addition, we provide a sample analysis using the steps outlined in the chapter.

Insights on the nuclear shuttling of H2A-H2B histone chaperones

All cellular processes that involve the unwinding of DNA also lead to the systematic shuttling of histones. Histone shuttling across the nuclear membrane is facilitated by a class of proteins known as - histone chaperones. Histone chaperones are classified based on their binding to H3/H4 histones or H2A/H2B histones. During the shuttling process, two types of signals - NLS and NES are recognized by the nuclear transport proteins. However, this is the nuclear transport protein and the mechanism of signal recognition by the protein is still unknown. Thus, in this piece of work, the NLS and NES signals are predicted on important H2A/H2B binding histone chaperones. In addition, cellular localization and potential DNA binding regions of histone chaperones are predicted. Mapping of predicted regions on the histone chaperone's structure suggested that the critical binding regions mainly lie on the disordered region of the histone chaperones. NLS and NES are present in the N- and C-terminal of the histone chaperones. Most histone chaperones contain bipartiate NLS signals. This article sheds light on the crucial aspect that in addition of being directly engaged in nucleosome synthesis and disassembly in vivo, histone chaperone also performs various specific roles via histone binding activity.

Human histone H1 variants impact splicing outcome by controlling RNA polymerase II elongation

Histones shape chromatin structure and the epigenetic landscape. H1, the most diverse histone in the human genome, has 11 variants. Due to the high structural similarity between the H1s, their unique functions in transferring information from the chromatin to mRNA-processing machineries have remained elusive. Here, we generated human cell lines lacking up to five H1 subtypes, allowing us to characterize the genomic binding profiles of six H1 variants. Most H1s bind to specific sites, and binding depends on multiple factors, including GC content. The highly expressed H1.2 has a high affinity for exons, whereas H1.3 binds intronic sequences. H1s are major splicing regulators, especially of exon skipping and intron retention events, through their effects on the elongation of RNA polymerase II (RNAPII). Thus, H1 variants determine splicing fate by modulating RNAPII elongation.

Asymmetric distribution of parental H3K9me3 in S phase silences L1 elements

In eukaryotes, repetitive DNA sequences are transcriptionally silenced through histone H3 lysine 9 trimethylation (H3K9me3). Loss of silencing of the repeat elements leads to genome instability and human diseases, including cancer and ageing1-3. Although the role of H3K9me3 in the establishment and maintenance of heterochromatin silencing has been extensively studied4-6, the pattern and mechanism that underlie the partitioning of parental H3K9me3 at replicating DNA strands are unknown. Here we report that H3K9me3 is preferentially transferred onto the leading strands of replication forks, which occurs predominantly at long interspersed nuclear element (LINE) retrotransposons (also known as LINE-1s or L1s) that are theoretically transcribed in the head-on direction with replication fork movement. Mechanistically, the human silencing hub (HUSH) complex interacts with the leading-strand DNA polymerase Pol ε and contributes to the asymmetric segregation of H3K9me3. Cells deficient in Pol ε subunits (POLE3 and POLE4) or the HUSH complex (MPP8 and TASOR) show compromised H3K9me3 asymmetry and increased LINE expression. Similar results were obtained in cells expressing a MPP8 mutant defective in H3K9me3 binding and in TASOR mutants with reduced interactions with Pol ε. These results reveal an unexpected mechanism whereby the HUSH complex functions with Pol ε to promote asymmetric H3K9me3 distribution at head-on LINEs to suppress their expression in S phase.

CTCF/cohesin organize the ground state of chromatin-nuclear speckle association

The interchromatin space in the cell nucleus contains various membrane-less nuclear bodies. Recent findings indicate that nuclear speckles, comprising a distinct nuclear body, exhibit interactions with certain chromatin regions in a ground state. Key questions are how this ground state of chromatin-nuclear speckle association is established and what are the gene regulatory roles of this layer of nuclear organization. We report here that chromatin structural factors CTCF and cohesin are required for full ground state association between DNA and nuclear speckles. Disruption of ground state DNA-speckle contacts via either CTCF depletion or cohesin depletion had minor effects on basal level expression of speckle-associated genes, however we show strong negative effects on stimulus-dependent induction of speckle-associated genes. We identified a putative speckle targeting motif (STM) within cohesin subunit RAD21 and demonstrated that the STM is required for chromatin-nuclear speckle association. In contrast to reduction of CTCF or RAD21, depletion of the cohesin releasing factor WAPL stabilized cohesin on chromatin and DNA-speckle contacts, resulting in enhanced inducibility of speckle-associated genes. In addition, we observed disruption of chromatin-nuclear speckle association in patient derived cells with Cornelia de Lange syndrome (CdLS), a congenital neurodevelopmental diagnosis involving defective cohesin pathways, thus revealing nuclear speckles as an avenue for therapeutic inquiry. In summary, our findings reveal a mechanism to establish the ground organizational state of chromatin-speckle association, to promote gene inducibility, and with relevance to human disease.

Discovery of a highly potent, selective, orally bioavailable inhibitor of KAT6A/B histone acetyltransferases with efficacy against KAT6A-high ER+ breast cancer

KAT6A, and its paralog KAT6B, are histone lysine acetyltransferases (HAT) that acetylate histone H3K23 and exert an oncogenic role in several tumor types including breast cancer where KAT6A is frequently amplified/overexpressed. However, pharmacologic targeting of KAT6A to achieve therapeutic benefit has been a challenge. Here we describe identification of a highly potent, selective, and orally bioavailable KAT6A/KAT6B inhibitor CTx-648 (PF-9363), derived from a benzisoxazole series, which demonstrates anti-tumor activity in correlation with H3K23Ac inhibition in KAT6A over-expressing breast cancer. Transcriptional and epigenetic profiling studies show reduced RNA Pol II binding and downregulation of genes involved in estrogen signaling, cell cycle, Myc and stem cell pathways associated with CTx-648 anti-tumor activity in ER-positive (ER+) breast cancer. CTx-648 treatment leads to potent tumor growth inhibition in ER+ breast cancer in vivo models, including models refractory to endocrine therapy, highlighting the potential for targeting KAT6A in ER+ breast cancer.

Lineage Plasticity and Stemness Phenotypes in Prostate Cancer: Harnessing the Power of Integrated "Omics" Approaches to Explore Measurable Metrics

Prostate cancer (PCa), the most frequent and second most lethal cancer type in men in developed countries, is a highly heterogeneous disease. PCa heterogeneity, therapy resistance, stemness, and lethal progression have been attributed to lineage plasticity, which refers to the ability of neoplastic cells to undergo phenotypic changes under microenvironmental pressures by switching between developmental cell states. What remains to be elucidated is how to identify measurements of lineage plasticity, how to implement them to inform preclinical and clinical research, and, further, how to classify patients and inform therapeutic strategies in the clinic. Recent research has highlighted the crucial role of next-generation sequencing technologies in identifying potential biomarkers associated with lineage plasticity. Here, we review the genomic, transcriptomic, and epigenetic events that have been described in PCa and highlight those with significance for lineage plasticity. We further focus on their relevance in PCa research and their benefits in PCa patient classification. Finally, we explore ways in which bioinformatic analyses can be used to determine lineage plasticity based on large omics analyses and algorithms that can shed light on upstream and downstream events. Most importantly, an integrated multiomics approach may soon allow for the identification of a lineage plasticity signature, which would revolutionize the molecular classification of PCa patients.

SNIP1 and PRC2 coordinate cell fates of neural progenitors during brain development

Stem cell survival versus death is a developmentally programmed process essential for morphogenesis, sizing, and quality control of genome integrity and cell fates. Cell death is pervasive during development, but its programming is little known. Here, we report that Smad nuclear interacting protein 1 (SNIP1) promotes neural progenitor cell survival and neurogenesis and is, therefore, integral to brain development. The SNIP1-depleted brain exhibits dysplasia with robust induction of caspase 9-dependent apoptosis. Mechanistically, SNIP1 regulates target genes that promote cell survival and neurogenesis, and its activities are influenced by TGFβ and NFκB signaling pathways. Further, SNIP1 facilitates the genomic occupancy of Polycomb complex PRC2 and instructs H3K27me3 turnover at target genes. Depletion of PRC2 is sufficient to reduce apoptosis and brain dysplasia and to partially restore genetic programs in the SNIP1-depleted brain in vivo. These findings suggest a loci-specific regulation of PRC2 and H3K27 marks to toggle cell survival and death in the developing brain.

Mutant FUS induces chromatin reorganization in the hippocampus and alters memory processes

Cytoplasmic mislocalization of the nuclear Fused in Sarcoma (FUS) protein is associated to amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Cytoplasmic FUS accumulation is recapitulated in the frontal cortex and spinal cord of heterozygous Fus∆NLS/+ mice. Yet, the mechanisms linking FUS mislocalization to hippocampal function and memory formation are still not characterized. Herein, we show that in these mice, the hippocampus paradoxically displays nuclear FUS accumulation. Multi-omic analyses showed that FUS binds to a set of genes characterized by the presence of an ETS/ELK-binding motifs, and involved in RNA metabolism, transcription, ribosome/mitochondria and chromatin organization. Importantly, hippocampal nuclei showed a decompaction of the neuronal chromatin at highly expressed genes and an inappropriate transcriptomic response was observed after spatial training of Fus∆NLS/+ mice. Furthermore, these mice lacked precision in a hippocampal-dependent spatial memory task and displayed decreased dendritic spine density. These studies shows that mutated FUS affects epigenetic regulation of the chromatin landscape in hippocampal neurons, which could participate in FTD/ALS pathogenic events. These data call for further investigation in the neurological phenotype of FUS-related diseases and open therapeutic strategies towards epigenetic drugs.

Next-Generation Sequencing Technology: Current Trends and Advancements

The advent of next-generation sequencing (NGS) has brought about a paradigm shift in genomics research, offering unparalleled capabilities for analyzing DNA and RNA molecules in a high-throughput and cost-effective manner. This transformative technology has swiftly propelled genomics advancements across diverse domains. NGS allows for the rapid sequencing of millions of DNA fragments simultaneously, providing comprehensive insights into genome structure, genetic variations, gene expression profiles, and epigenetic modifications. The versatility of NGS platforms has expanded the scope of genomics research, facilitating studies on rare genetic diseases, cancer genomics, microbiome analysis, infectious diseases, and population genetics. Moreover, NGS has enabled the development of targeted therapies, precision medicine approaches, and improved diagnostic methods. This review provides an insightful overview of the current trends and recent advancements in NGS technology, highlighting its potential impact on diverse areas of genomic research. Moreover, the review delves into the challenges encountered and future directions of NGS technology, including endeavors to enhance the accuracy and sensitivity of sequencing data, the development of novel algorithms for data analysis, and the pursuit of more efficient, scalable, and cost-effective solutions that lie ahead.

Impaired histone inheritance promotes tumor progression

Faithful inheritance of parental histones is essential to maintain epigenetic information and cellular identity during cell division. Parental histones are evenly deposited onto the replicating DNA of sister chromatids in a process dependent on the MCM2 subunit of DNA helicase. However, the impact of aberrant parental histone partition on human disease such as cancer is largely unknown. In this study, we construct a model of impaired histone inheritance by introducing MCM2-2A mutation (defective in parental histone binding) in MCF-7 breast cancer cells. The resulting impaired histone inheritance reprograms the histone modification landscapes of progeny cells, especially the repressive histone mark H3K27me3. Lower H3K27me3 levels derepress the expression of genes associated with development, cell proliferation, and epithelial to mesenchymal transition. These epigenetic changes confer fitness advantages to some newly emerged subclones and consequently promote tumor growth and metastasis after orthotopic implantation. In summary, our results indicate that impaired inheritance of parental histones can drive tumor progression.

ChroKit: a Shiny-based framework for interactive analysis, visualization and integration of genomic data

We developed ChroKit (the Chromatin toolKit), an interactive web-based framework written in R that enables intuitive exploration, multidimensional analyses, and visualization of genomic data from ChIP-Seq, DNAse-Seq or any other NGS experiment that reports the enrichment of aligned reads over genomic regions. This program takes preprocessed NGS data and performs operations on genomic regions of interest, including resetting their boundaries, their annotation based on proximity to genomic features, the association to gene ontologies, and signal enrichment calculations. Genomic regions can be further refined or subsetted by user-defined logical operations and unsupervised classification algorithms. ChroKit generates a full range of plots that are easily manipulated by point and click operations, thus allowing 'on the fly' re-analysis and fast exploration of the data. Working sessions can be exported for reproducibility, accountability, and easy sharing within the bioinformatics community. ChroKit is multiplatform and can be deployed on a server to enhance computational speed and provide simultaneous access by multiple users. ChroKit is a fast and intuitive genomic analysis tool suited for a wide range of users due to its architecture and its user-friendly graphical interface. ChroKit source code is available at https://github.com/ocroci/ChroKit and the Docker image at https://hub.docker.com/r/ocroci/chrokit.

Phosphorylation of 53BP1 by ATM enforce neurodevelopmental programs in cortical organoids

53BP1 is a well-established DNA damage repair factor recently shown to regulate gene expression and critically influence tumor suppression and neural development. For gene regulation, how 53BP1 is regulated remains unclear. Here, we showed that 53BP1-serine 25 phosphorylation by ATM is required for neural progenitor cell proliferation and neuronal differentiation in cortical organoids. 53BP1-serine 25 phosphorylation dynamics controls 53BP1 target genes for neuronal differentiation and function, cellular response to stress, and apoptosis. Beyond 53BP1, ATM is required for phosphorylation of factors in neuronal differentiation, cytoskeleton, p53 regulation, and ATM, BNDF, and WNT signaling pathways for cortical organoid differentiation. Overall, our data suggest that 53BP1 and ATM control key genetic programs required for human cortical development.

High-throughput methods for the analysis of transcription factors and chromatin modifications: Low input, single cell and spatial genomic technologies

Genome-wide analysis of transcription factors and epigenomic features is instrumental to shed light on DNA-templated regulatory processes such as transcription, cellular differentiation or to monitor cellular responses to environmental cues. Two decades of technological developments have led to a rich set of approaches progressively pushing the limits of epigenetic profiling towards single cells. More recently, disruptive technologies using innovative biochemistry came into play. Assays such as CUT&RUN, CUT&Tag and variations thereof show considerable potential to survey multiple TFs or histone modifications in parallel from a single experiment and in native conditions. These are in the path to become the dominant assays for genome-wide analysis of TFs and chromatin modifications in bulk, single-cell, and spatial genomic applications. The principles together with pros and cons are discussed.

Loss of Ezh2 function remodels the DNA replication initiation landscape

In metazoan cells, DNA replication initiates from thousands of genomic loci scattered throughout the genome called DNA replication origins. Origins are strongly associated with euchromatin, particularly open genomic regions such as promoters and enhancers. However, over a third of transcriptionally silent genes are associated with DNA replication initiation. Most of these genes are bound and repressed by the Polycomb repressive complex-2 (PRC2) through the repressive H3K27me3 mark. This is the strongest overlap observed for a chromatin regulator with replication origin activity. Here, we asked whether Polycomb-mediated gene repression is functionally involved in recruiting DNA replication origins to transcriptionally silent genes. We show that the absence of EZH2, the catalytic subunit of PRC2, results in increased DNA replication initiation, specifically in the vicinity of EZH2 binding sites. The increase in DNA replication initiation does not correlate with transcriptional de-repression or the acquisition of activating histone marks but does correlate with loss of H3K27me3 from bivalent promoters.

NSD2 E1099K drives relapse in pediatric acute lymphoblastic leukemia by disrupting 3D chromatin organization

BACKGROUND

The NSD2 p.E1099K (EK) mutation is shown to be enriched in patients with relapsed acute lymphoblastic leukemia (ALL), indicating a role in clonal evolution and drug resistance.

RESULTS

To uncover 3D chromatin architecture-related mechanisms underlying drug resistance, we perform Hi-C on three B-ALL cell lines heterozygous for NSD2 EK. The NSD2 mutation leads to widespread remodeling of the 3D genome, most dramatically in terms of compartment changes with a strong bias towards A compartment shifts. Systematic integration of the Hi-C data with previously published ATAC-seq, RNA-seq, and ChIP-seq data show an expansion in H3K36me2 and a shrinkage in H3K27me3 within A compartments as well as increased gene expression and chromatin accessibility. These results suggest that NSD2 EK plays a prominent role in chromatin decompaction through enrichment of H3K36me2. In contrast, we identify few changes in intra-topologically associating domain activity. While compartment changes vary across cell lines, a common core of decompacting loci are shared, driving the expression of genes/pathways previously implicated in drug resistance. We further perform RNA sequencing on a cohort of matched diagnosis/relapse ALL patients harboring the relapse-specific NSD2 EK mutation. Changes in patient gene expression upon relapse significantly correlate with core compartment changes, further implicating the role of NSD2 EK in genome decompaction.

CONCLUSIONS

In spite of cell-context-dependent changes mediated by EK, there appears to be a shared transcriptional program dependent on compartment shifts which could explain phenotypic differences across EK cell lines. This core program is an attractive target for therapeutic intervention.

Epigenetic signals that direct cell type-specific interferon beta response in mouse cells

The antiviral response induced by type I interferon (IFN) via the JAK-STAT signaling cascade activates hundreds of IFN-stimulated genes (ISGs) across human and mouse tissues but varies between cell types. However, the links between the underlying epigenetic features and the ISG profile are not well understood. We mapped ISGs, binding sites of the STAT1 and STAT2 transcription factors, chromatin accessibility, and histone H3 lysine modification by acetylation (ac) and mono-/tri-methylation (me1, me3) in mouse embryonic stem cells and fibroblasts before and after IFNβ treatment. A large fraction of ISGs and STAT-binding sites was cell type specific with promoter binding of a STAT1/2 complex being a key driver of ISGs. Furthermore, STAT1/2 binding to putative enhancers induced ISGs as inferred from a chromatin co-accessibility analysis. STAT1/2 binding was dependent on the chromatin context and positively correlated with preexisting H3K4me1 and H3K27ac marks in an open chromatin state, whereas the presence of H3K27me3 had an inhibitory effect. Thus, chromatin features present before stimulation represent an additional regulatory layer for the cell type-specific antiviral response.

Adolescent binge ethanol impacts H3K36me3 regulation of synaptic genes

Adolescence is marked in part by the ongoing development of the prefrontal cortex (PFC). Binge ethanol use during this critical stage in neurodevelopment induces significant structural changes to the PFC, as well as cognitive and behavioral deficits that can last into adulthood. Previous studies showed that adolescent binge ethanol causes lasting deficits in working memory, decreases in the expression of chromatin remodeling genes responsible for the methylation of histone 3 lysine 36 (H3K36), and global decreases in H3K36 in the PFC. H3K36me3 is present within the coding region of actively-transcribed genes, and safeguards against aberrant, cryptic transcription by RNA Polymerase II. We hypothesize that altered methylation of H3K36 could play a role in adolescent binge ethanol-induced memory deficits. To investigate this at the molecular level, ethanol (4 g/kg, i.g.) or water was administered intermittently to adolescent mice. RNA-and ChIP-sequencing were then performed within the same tissue to determine gene expression changes and identify genes and loci where H3K36me3 was disrupted by ethanol. We further assessed ethanol-induced changes at the transcription level with differential exon-use and cryptic transcription analysis - a hallmark of decreased H3K36me3. Here, we found ethanol-induced changes to the gene expression and H3K36me3-regulation of synaptic-related genes in all our analyses. Notably, H3K36me3 was differentially trimethylated between ethanol and control conditions at synaptic-related genes, and Snap25 and Cplx1 showed evidence of cryptic transcription in males and females treated with ethanol during adolescence. Our results provide preliminary evidence that ethanol-induced changes to H3K36me3 during adolescent neurodevelopment may be linked to synaptic dysregulation at the transcriptional level, which may explain the reported ethanol-induced changes to PFC synaptic function.

Lysine methyltransferase Kmt2d regulates naive CD8+ T cell activation-induced survival

Lysine specific methyltransferase 2D (Kmt2d) catalyzes the mono-methylation of histone 3 lysine 4 (H3K4me1) and plays a critical role in regulatory T cell generation via modulating Foxp3 gene expression. Here we report a role of Kmt2d in naïve CD8⁺ T cell generation and survival. In the absence of Kmt2d, the number of CD8⁺ T cells, particularly naïve CD8⁺ T cells (CD62L^hi/CD44^lo), in spleen was greatly decreased and in vitro activation-related death significantly increased from Kmt2d ^fl/flCD4cre⁺ (KO) compared to Kmt2d ^fl/flCD4cre^- (WT) mice. Furthermore, analyses by ChIPseq, RNAseq, and scRNAseq showed reduced H3K4me1 levels in enhancers and reduced expression of apoptosis-related genes in activated naïve CD8⁺ T cells in the absence of Kmt2d. Finally, we confirmed the activation-induced death of antigen-specific naïve CD8⁺ T cells in vivo in Kmt2d KO mice upon challenge with Listeria monocytogenes infection. These findings reveal that Kmt2d regulates activation-induced naïve CD8⁺ T cell survival via modulating H3K4me1 levels in enhancer regions of apoptosis and immune function-related genes.

USP7 regulates the ncPRC1 Polycomb axis to stimulate genomic H2AK119ub1 deposition uncoupled from H3K27me3

Ubiquitin-specific protease 7 (USP7) has been implicated in cancer progression and neurodevelopment. However, its molecular targets remain poorly characterized. We combined quantitative proteomics, transcriptomics, and epigenomics to define the core USP7 network. Our multi-omics analysis reveals USP7 as a control hub that links genome regulation, tumor suppression, and histone H2A ubiquitylation (H2AK119ub1) by noncanonical Polycomb-repressive complexes (ncPRC1s). USP7 strongly stabilizes ncPRC1.6 and, to a lesser extent, ncPRC1.1. Moreover, USP7 represses expression of AUTS2, which suppresses H2A ubiquitylation by ncPRC1.3/5. Collectively, these USP7 activities promote the genomic deposition of H2AK119ub1 by ncPRC1, especially at transcriptionally repressed loci. Notably, USP7-dependent changes in H2AK119ub1 levels are uncoupled from H3K27me3. Even complete loss of the PRC1 catalytic core and H2AK119ub1 has only a limited effect on H3K27me3. Besides defining the USP7 regulome, our results reveal that H2AK119ub1 dosage is largely disconnected from H3K27me3.

LncRNA CCTT-mediated RNA-DNA and RNA-protein interactions facilitate the recruitment of CENP-C to centromeric DNA during kinetochore assembly

Kinetochore assembly on centromeres is central for chromosome segregation, and defects in this process cause mitotic errors and aneuploidy. Besides the well-established protein network, emerging evidence suggests the involvement of regulatory RNA in kinetochore assembly; however, it has remained elusive about the identity of such RNA, let alone its mechanism of action in this critical process. Here, we report CCTT, a previously uncharacterized long non-coding RNA (lncRNA) transcribed from the arm of human chromosome 17, which plays a vital role in kinetochore assembly. We show that CCTT highly localizes to all centromeres via the formation of RNA-DNA triplex and specifically interacts with CENP-C to help engage this blueprint protein in centromeres, and consequently, CCTT loss triggers extensive mitotic errors and aneuploidy. These findings uncover a non-centromere-derived lncRNA that recruits CENP-C to centromeres and shed critical lights on the function of centromeric DNA sequences as anchor points for kinetochore assembly.

Redistribution of lamina-associated domains reshapes binding of pioneer factor FOXA2 in development of nonalcoholic fatty liver disease

Nonalcoholic fatty liver disease (NAFLD) is highly prevalent in type 2 diabetes mellitus and the elderly, impacting 40% of individuals over 70. Regulation of heterochromatin at the nuclear lamina has been associated with aging and age-dependent metabolic changes. We previously showed that changes at the lamina in aged hepatocytes and laminopathy models lead to redistribution of lamina-associated domains (LADs), opening of repressed chromatin, and up-regulation of genes regulating lipid synthesis and storage, culminating in fatty liver. Here, we test the hypothesis that change in the expression of lamina-associated proteins and nuclear shape leads to redistribution of LADs, followed by altered binding of pioneer factor FOXA2 and by up-regulation of lipid synthesis and storage, culminating in steatosis in younger NAFLD patients (aged 21-51). Changes in nuclear morphology alter LAD partitioning and reduced lamin B1 signal correlate with increased FOXA2 binding before severe steatosis in young mice placed on a western diet. Nuclear shape is also changed in younger NAFLD patients. LADs are redistrubted and lamin B1 signal decreases similarly in mild and severe steatosis. In contrast, FOXA2 binding is similar in normal and NAFLD patients with moderate steatosis and is repositioned only in NAFLD patients with more severe lipid accumulation. Hence, changes at the nuclear lamina reshape FOXA2 binding with progression of the disease. Our results suggest a role for nuclear lamina in etiology of NAFLD, irrespective of aging, with potential for improved stratification of patients and novel treatments aimed at restoring nuclear lamina function.

Exploitation of epigenetic variation of crop wild relatives for crop improvement and agrobiodiversity preservation

Crop wild relatives (CWRs) are recognized as the best potential source of traits for crop improvement. However, successful crop improvement using CWR relies on identifying variation in genes controlling desired traits in plant germplasms and subsequently incorporating them into cultivars. Epigenetic diversity may provide an additional layer of variation within CWR and can contribute novel epialleles for key traits for crop improvement. There is emerging evidence that epigenetic variants of functional and/or agronomic importance exist in CWR gene pools. This provides a rationale for the conservation of epigenotypes of interest, thus contributing to agrobiodiversity preservation through conservation and (epi)genetic monitoring. Concepts and techniques of classical and modern breeding should consider integrating recent progress in epigenetics, initially by identifying their association with phenotypic variations and then by assessing their heritability and stability in subsequent generations. New tools available for epigenomic analysis offer the opportunity to capture epigenetic variation and integrate it into advanced (epi)breeding programmes. Advances in -omics have provided new insights into the sources and inheritance of epigenetic variation and enabled the efficient introduction of epi-traits from CWR into crops using epigenetic molecular markers, such as epiQTLs.

A noncoding single-nucleotide polymorphism at 8q24 drives IDH1-mutant glioma formation

Establishing causal links between inherited polymorphisms and cancer risk is challenging. Here, we focus on the single-nucleotide polymorphism rs55705857, which confers a sixfold greater risk of isocitrate dehydrogenase (IDH)-mutant low-grade glioma (LGG). We reveal that rs55705857 itself is the causal variant and is associated with molecular pathways that drive LGG. Mechanistically, we show that rs55705857 resides within a brain-specific enhancer, where the risk allele disrupts OCT2/4 binding, allowing increased interaction with the Myc promoter and increased Myc expression. Mutating the orthologous mouse rs55705857 locus accelerated tumor development in an Idh1R132H-driven LGG mouse model from 472 to 172 days and increased penetrance from 30% to 75%. Our work reveals mechanisms of the heritable predisposition to lethal glioma in ~40% of LGG patients.

PBRM1 presents a potential prognostic marker and therapeutic target in duodenal papillary carcinoma

BACKGROUND

Due to its rarity, duodenal papillary carcinoma (DPC) is seldom studied as a unique disease and no specific molecular features or treatment guidelines are provided.

METHODS

Whole-exome sequencing was performed to gain new insights into the DPC mutation landscape and to identify potential signalling pathways and therapeutic targets. Mechanistically, immunohistochemistry (IHC), immunofluorescence, RNA-seq, ATAC-seq and in vitro cell function experiments were performed to confirm the underlying mechanisms.

RESULTS

We described the mutational landscape of DPC for the first time as a group of rare tumours with a high frequency of dysregulation in the chromatin remodelling pathway, particularly PBRM1-inactivating mutations that are significantly higher than duodenal adenocarcinomas and ampullary adenocarcinoma (27% vs. 0% vs. 7%, p < .01). In vitro cell experiments showed that downregulation of PBRM1 expression could significantly promote the cancer progression and epithelial-to-mesenchymal transition via the PBRM1-c-JUN-VIM axis. The IHC data indicated that PBRM1 deficiency (p = .047) and c-JUN expression (p < .001) were significantly associated with poor prognosis. Meanwhile, the downregulation of PBRM1 expression in HUTU-80 cells was sensitive to radiation, which may be due to the suppression of c-JUN by irradiation.

CONCLUSIONS

Our findings define a novel molecular subgroup of PBRM1-inactivating mutations in DPC. PBRM1 play an important role in DPC progression and may serve as a potential therapeutic target and prognostic indicator.

Altered activity-regulated H3K9 acetylation at TGF-beta signaling genes during egocentric memory in Huntington's disease

Molecular mechanisms underlying cognitive deficits in Huntington's disease (HD), a striatal neurodegenerative disorder, are unknown. Here, we generated ChIPseq, 4Cseq and RNAseq data on striatal tissue of HD and control mice during striatum-dependent egocentric memory process. Multi-omics analyses showed altered activity-dependent epigenetic gene reprogramming of neuronal and glial genes regulating striatal plasticity in HD mice, which correlated with memory deficit. First, our data reveal that spatial chromatin re-organization and transcriptional induction of BDNF-related markers, regulating neuronal plasticity, were reduced since memory acquisition in the striatum of HD mice. Second, our data show that epigenetic memory implicating H3K9 acetylation, which established during late phase of memory process (e.g. during consolidation/recall) and contributed to glia-mediated, TGFβ-dependent plasticity, was compromised in HD mouse striatum. Specifically, memory-dependent regulation of H3K9 acetylation was impaired at genes controlling extracellular matrix and myelination. Our study investigating the interplay between epigenetics and memory identifies H3K9 acetylation and TGFβ signaling as new targets of striatal plasticity, which might offer innovative leads to improve HD.