ChIPseqSpikeInFree: a ChIP-seq normalization approach to reveal global changes in histone modifications without spike-in

Histone lysine demethylase 4 family proteins maintain the transcriptional program and adrenergic cellular state of MYCN-amplified neuroblastoma

Neuroblastoma with MYCN amplification (MNA) is a high-risk disease that has a poor survival rate. Neuroblastoma displays cellular heterogeneity, including more differentiated (adrenergic) and more primitive (mesenchymal) cellular states. Here, we demonstrate that MYCN oncoprotein promotes a cellular state switch in mesenchymal cells to an adrenergic state, accompanied by induction of histone lysine demethylase 4 family members (KDM4A-C) that act in concert to control the expression of MYCN and adrenergic core regulatory circulatory (CRC) transcription factors. Pharmacologic inhibition of KDM4 blocks expression of MYCN and the adrenergic CRC transcriptome with genome-wide induction of transcriptionally repressive H3K9me3, resulting in potent anticancer activity against neuroblastomas with MNA by inducing neuroblastic differentiation and apoptosis. Furthermore, a short-term KDM4 inhibition in combination with conventional, cytotoxic chemotherapy results in complete tumor responses of xenografts with MNA. Thus, KDM4 blockade may serve as a transformative strategy to target the adrenergic CRC dependencies in MNA neuroblastomas.

Single-cell chromatin profiling reveals genetic programs activating proregenerative states in nonmyocyte cells

Cardiac regeneration requires coordinated participation of multiple cell types whereby their communications result in transient activation of proregenerative cell states. Although the molecular characteristics and lineage origins of these activated cell states and their contribution to cardiac regeneration have been studied, the extracellular signaling and the intrinsic genetic program underlying the activation of the transient functional cell states remain largely unexplored. In this study, we delineated the chromatin landscapes of the noncardiomyocytes (nonCMs) of the regenerating heart at the single-cell level and inferred the cis-regulatory architectures and trans-acting factors that control cell type-specific gene expression programs. Moreover, further motif analysis and cell-specific genetic manipulations suggest that the macrophage-derived inflammatory signal tumor necrosis factor-α, acting via its downstream transcription factor complex activator protein-1, functions cooperatively with discrete transcription regulators to activate respective nonCM cell types critical for cardiac regeneration. Thus, our study defines the regulatory architectures and intercellular communication principles in zebrafish heart regeneration.

Type I MADS-box transcription factor TaMADS-GS regulates grain size by stabilizing cytokinin signalling during endosperm cellularization in wheat

Grain size is one of the important traits in wheat breeding programs aimed at improving yield, and cytokinins, mainly involved in cell division, have a positive impact on grain size. Here, we identified a novel wheat gene TaMADS-GS encoding type I MADS-box transcription factor, which regulates the cytokinins signalling pathway during early stages of grain development to modulate grain size and weight in wheat. TaMADS-GS is exclusively expressed in grains at early stage of seed development and its knockout leads to delayed endosperm cellularization, smaller grain size and lower grain weight. TaMADS-GS protein interacts with the Polycomb Repressive Complex 2 (PRC2) and leads to repression of genes encoding cytokinin oxidase/dehydrogenases (CKXs) stimulating cytokinins inactivation by mediating accumulation of the histone H3 trimethylation at lysine 27 (H3K27me3). Through the screening of a large wheat germplasm collection, an elite allele of the TaMADS-GS exhibits higher ability to repress expression of genes inactivating cytokinins and a positive correlation with grain size and weight, thus representing a novel marker for breeding programs in wheat. Overall, these findings support the relevance of TaMADS-GS as a key regulator of wheat grain size and weight.

Reduced female fertility due to sequestration of RNA Pol II by pervasive transcription in exosome RNase-depleted oocytes

Perturbing the transcriptome of mammalian oocytes results in meiotic failure. We previously reported that RNA-exosome-associated RNase, EXOSC10, degrades unwanted protein-coding RNA and processes ribosomal RNA to ensure proper oocyte maturation. Here, we establish oocyte-specific knockout mice of another RNA-exosome-associated RNase, DIS3. Mutant females (Dis3cKO) exhibit significantly reduced fertility because oocytes arrest after the growth phase. Single-oocyte RNA sequencing (RNA-seq) and CUT&Tag analyses show that DIS3 degrades intergenic RNA and mediates transcription silencing that is essential for chromatin condensation and resumption of meiosis. Dis3cKO oocytes exhibit elevated H3K27me3 in a pre-defined manner due to insufficient demethylation. During oocyte growth, EXOSC10 functions with DIS3 to degrade intergenic RNA. Double-knockout oocytes have earlier growth defects and more accumulated transcripts. We conclude that RNA exosomes synergistically degrade unwanted RNA and mediate transcription termination to ensure transcriptome integrity during oocyte development.

PR-DUB safeguards Polycomb repression through H2AK119ub1 restriction

Polycomb group (PcG) proteins are critical chromatin regulators for cell fate control. The mono-ubiquitylation on histone H2AK119 (H2AK119ub1) is one of the well-recognized mechanisms for Polycomb repressive complex 1 (PRC1)-mediated transcription repression. Unexpectedly, the specific H2AK119 deubiquitylation complex composed by additional sex comb-like proteins and BAP1 has also been genetically characterized as Polycomb repressive deubiquitnase (PR-DUB) for unclear reasons. However, it remains a mystery whether and how PR-DUB deficiency affects chromatin states and cell fates through impaired PcG silencing. Here through a careful epigenomic analysis, we demonstrate that a bulk of H2AK119ub1 is diffusely distributed away from promoter regions and their enrichment is positively correlated with PRC1 occupancy. Upon deletion of Asxl2 in mouse embryonic stem cells (ESCs), a pervasive gain of H2AK119ub1 is coincident with increased PRC1 sampling at chromatin. Accordingly, PRC1 is significantly lost from a subset of highly occupied promoters, leading to impaired silencing of associated genes before and after lineage differentiation of Asxl2-null ESCs. Therefore, our study highlights the importance of genome-wide H2AK119ub1 restriction by PR-DUB in safeguarding robust PRC1 deposition and its roles in developmental regulation.

Regulatory architecture of housekeeping genes is driven by promoter assemblies

Genes that are key to cell identity are generally regulated by cell-type-specific enhancer elements bound by transcription factors, some of which facilitate looping to distant gene promoters. In contrast, genes that encode housekeeping functions, whose regulation is essential for normal cell metabolism and growth, generally lack interactions with distal enhancers. We find that Ronin (Thap11) assembles multiple promoters of housekeeping and metabolic genes to regulate gene expression. This behavior is analogous to how enhancers are brought together with promoters to regulate cell identity genes. Thus, Ronin-dependent promoter assemblies provide a mechanism to explain why housekeeping genes can forgo distal enhancer elements and why Ronin is important for cellular metabolism and growth control. We propose that clustering of regulatory elements is a mechanism common to cell identity and housekeeping genes but is accomplished by different factors binding distinct control elements to establish enhancer-promoter or promoter-promoter interactions, respectively.

Extracellular acidosis restricts one-carbon metabolism and preserves T cell stemness

The accumulation of acidic metabolic waste products within the tumor microenvironment inhibits effector functions of tumor-infiltrating lymphocytes (TILs). However, it remains unclear how an acidic environment affects T cell metabolism and differentiation. Here we show that prolonged exposure to acid reprograms T cell intracellular metabolism and mitochondrial fitness and preserves T cell stemness. Mechanistically, elevated extracellular acidosis impairs methionine uptake and metabolism via downregulation of SLC7A5, therefore altering H3K27me3 deposition at the promoters of key T cell stemness genes. These changes promote the maintenance of a 'stem-like memory' state and improve long-term in vivo persistence and anti-tumor efficacy in mice. Our findings not only reveal an unexpected capacity of extracellular acidosis to maintain the stem-like properties of T cells, but also advance our understanding of how methionine metabolism affects T cell stemness.

Dynamic chromatin regulatory programs during embryogenesis of hexaploid wheat

BACKGROUND

Plant and animal embryogenesis have conserved and distinct features. Cell fate transitions occur during embryogenesis in both plants and animals. The epigenomic processes regulating plant embryogenesis remain largely elusive.

RESULTS

Here, we elucidate chromatin and transcriptomic dynamics during embryogenesis of the most cultivated crop, hexaploid wheat. Time-series analysis reveals stage-specific and proximal-distal distinct chromatin accessibility and dynamics concordant with transcriptome changes. Following fertilization, the remodeling kinetics of H3K4me3, H3K27ac, and H3K27me3 differ from that in mammals, highlighting considerable species-specific epigenomic dynamics during zygotic genome activation. Polycomb repressive complex 2 (PRC2)-mediated H3K27me3 deposition is important for embryo establishment. Later H3K27ac, H3K27me3, and chromatin accessibility undergo dramatic remodeling to establish a permissive chromatin environment facilitating the access of transcription factors to cis-elements for fate patterning. Embryonic maturation is characterized by increasing H3K27me3 and decreasing chromatin accessibility, which likely participates in restricting totipotency while preventing extensive organogenesis. Finally, epigenomic signatures are correlated with biased expression among homeolog triads and divergent expression after polyploidization, revealing an epigenomic contributor to subgenome diversification in an allohexaploid genome.

CONCLUSIONS

Collectively, we present an invaluable resource for comparative and mechanistic analysis of the epigenomic regulation of crop embryogenesis.

Statistical Analysis in ChIP-seq-Related Applications

Chromatin immunoprecipitation sequencing (ChIP-seq) has been widely performed to identify protein binding information along the genome. The sequencing protocol is quite flexible and mature to measure different types of protein binding as long as sequencing parameters are properly tailored to accommodate protein features. Two distinct types of protein binding are point-source-like binding by transcription factors and diffused-distribution binding by histone modifications. Consequently, statistical approaches have been proposed to address ChIP-seq-related questions according to different protein features. In this chapter, we briefly summarize statistical principles, approaches, and tools that are widely implemented in modeling ChIP-seq data, from raw data quality control to final result reporting. We discuss the key solutions in addressing eight routine questions in ChIP-seq applications. We also include discussion on approaches fitting unique data features in different ChIP-seq types. We hope this chapter will serve as a brief guide, especially for ChIP-seq beginners, to provide them with a high-level overview to understand and design processing plans for their ChIP-seq experiments.

Targeting KDM4 for treating PAX3-FOXO1-driven alveolar rhabdomyosarcoma

Chimeric transcription factors drive lineage-specific oncogenesis but are notoriously difficult to target. Alveolar rhabdomyosarcoma (RMS) is an aggressive childhood soft tissue sarcoma transformed by the pathognomonic Paired Box 3-Forkhead Box O1 (PAX3-FOXO1) fusion protein, which governs a core regulatory circuitry transcription factor network. Here, we show that the histone lysine demethylase 4B (KDM4B) is a therapeutic vulnerability for PAX3-FOXO1⁺ RMS. Genetic and pharmacologic inhibition of KDM4B substantially delayed tumor growth. Suppression of KDM4 proteins inhibited the expression of core oncogenic transcription factors and caused epigenetic alterations of PAX3-FOXO1-governed superenhancers. Combining KDM4 inhibition with cytotoxic chemotherapy led to tumor regression in preclinical PAX3-FOXO1⁺ RMS subcutaneous xenograft models. In summary, we identified a targetable mechanism required for maintenance of the PAX3-FOXO1-related transcription factor network, which may translate to a therapeutic approach for fusion-positive RMS.

Reorganization of postmitotic neuronal chromatin accessibility for maturation of serotonergic identity

Assembly of transcriptomes encoding unique neuronal identities requires selective accessibility of transcription factors to cis-regulatory sequences in nucleosome-embedded postmitotic chromatin. Yet, the mechanisms controlling postmitotic neuronal chromatin accessibility are poorly understood. Here, we show that unique distal enhancers define the Pet1 neuron lineage that generates serotonin (5-HT) neurons in mice. Heterogeneous single-cell chromatin landscapes are established early in postmitotic Pet1 neurons and reveal the putative regulatory programs driving Pet1 neuron subtype identities. Distal enhancer accessibility is highly dynamic as Pet1 neurons mature, suggesting the existence of regulatory factors that reorganize postmitotic neuronal chromatin. We find that Pet1 and Lmx1b control chromatin accessibility to select Pet1-lineage-specific enhancers for 5-HT neurotransmission. Additionally, these factors are required to maintain chromatin accessibility during early maturation suggesting that postmitotic neuronal open chromatin is unstable and requires continuous regulatory input. Together, our findings reveal postmitotic transcription factors that reorganize accessible chromatin for neuron specialization.

An adult-stage transcriptional program for survival of serotonergic connectivity

Neurons must function for decades of life, but how these non-dividing cells are preserved is poorly understood. Using mouse serotonin (5-HT) neurons as a model, we report an adult-stage transcriptional program specialized to ensure the preservation of neuronal connectivity. We uncover a switch in Lmx1b and Pet1 transcription factor function from controlling embryonic axonal growth to sustaining a transcriptomic signature of 5-HT connectivity comprising functionally diverse synaptic and axonal genes. Adult-stage deficiency of Lmx1b and Pet1 causes slowly progressing degeneration of 5-HT synapses and axons, increased susceptibility of 5-HT axons to neurotoxic injury, and abnormal stress responses. Axon degeneration occurs in a die back pattern and is accompanied by accumulation of α-synuclein and amyloid precursor protein in spheroids and mitochondrial fragmentation without cell body loss. Our findings suggest that neuronal connectivity is transcriptionally protected by maintenance of connectivity transcriptomes; progressive decay of such transcriptomes may contribute to age-related diseases of brain circuitry.

PROSER1 mediates TET2 O-GlcNAcylation to regulate DNA demethylation on UTX-dependent enhancers and CpG islands

PROSER1 promotes the interaction between TET2 and the glycosyltransferase OGT to regulate TET2 O-GlcNAcylation and stability on genomic elements that depend on the activity of the MLL3/4 complexes.

DNA methylation at enhancers and CpG islands usually leads to gene repression, which is counteracted by DNA demethylation through the TET protein family. However, how TET enzymes are recruited and regulated at these genomic loci is not fully understood. Here, we identify TET2, the glycosyltransferase OGT and a previously undescribed proline and serine rich protein, PROSER1 as interactors of UTX, a component of the enhancer-associated MLL3/4 complexes. We find that PROSER1 mediates the interaction between OGT and TET2, thus promoting TET2 O-GlcNAcylation and protein stability. In addition, PROSER1, UTX, TET1/2, and OGT colocalize on many genomic elements genome-wide. Loss of PROSER1 results in lower enrichment of UTX, TET1/2, and OGT at enhancers and CpG islands, with a concomitant increase in DNA methylation and transcriptional down-regulation of associated target genes and increased DNA hypermethylation encroachment at H3K4me1-predisposed CpG islands. Furthermore, we provide evidence that PROSER1 acts as a more general regulator of OGT activity by controlling O-GlcNAcylation of multiple other chromatin signaling pathways. Taken together, this study describes for the first time a regulator of TET2 O-GlcNAcylation and its implications in mediating DNA demethylation at UTX-dependent enhancers and CpG islands and supports an important role for PROSER1 in regulating the function of various chromatin-associated proteins via OGT-mediated O-GlcNAcylation.

CHIPIN: ChIP-seq inter-sample normalization based on signal invariance across transcriptionally constant genes

Background

Multiple studies rely on ChIP-seq experiments to assess the effect of gene modulation and drug treatments on protein binding and chromatin structure. However, most methods commonly used for the normalization of ChIP-seq binding intensity signals across conditions, e.g., the normalization to the same number of reads, either assume a constant signal-to-noise ratio across conditions or base the estimates of correction factors on genomic regions with intrinsically different signals between conditions. Inaccurate normalization of ChIP-seq signal may, in turn, lead to erroneous biological conclusions.

Results

We developed a new R package, CHIPIN, that allows normalizing ChIP-seq signals across different conditions/samples when spike-in information is not available, but gene expression data are at hand. Our normalization technique is based on the assumption that, on average, no differences in ChIP-seq signals should be observed in the regulatory regions of genes whose expression levels are constant across samples/conditions. In addition to normalizing ChIP-seq signals, CHIPIN provides as output a number of graphs and calculates statistics allowing the user to assess the efficiency of the normalization and qualify the specificity of the antibody used. In addition to ChIP-seq, CHIPIN can be used without restriction on open chromatin ATAC-seq or DNase hypersensitivity data. We validated the CHIPIN method on several ChIP-seq data sets and documented its superior performance in comparison to several commonly used normalization techniques.

Conclusions

The CHIPIN method provides a new way for ChIP-seq signal normalization across conditions when spike-in experiments are not available. The method is implemented in a user-friendly R package available on GitHub: https://github.com/BoevaLab/CHIPIN

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04320-3.

Metabolically controlled histone H4K5 acylation/acetylation ratio drives BRD4 genomic distribution

In addition to acetylation, histones are modified by a series of competing longer-chain acylations. Most of these acylation marks are enriched and co-exist with acetylation on active gene regulatory elements. Their seemingly redundant functions hinder our understanding of histone acylations' specific roles. Here, by using an acute lymphoblastic leukemia (ALL) cell model and blasts from individuals with B-precusor ALL (B-ALL), we demonstrate a role of mitochondrial activity in controlling the histone acylation/acetylation ratio, especially at histone H4 lysine 5 (H4K5). An increase in the ratio of non-acetyl acylations (crotonylation or butyrylation) over acetylation on H4K5 weakens bromodomain containing protein 4 (BRD4) bromodomain-dependent chromatin interaction and enhances BRD4 nuclear mobility and availability for binding transcription start site regions of active genes. Our data suggest that the metabolism-driven control of the histone acetylation/longer-chain acylation(s) ratio could be a common mechanism regulating the bromodomain factors' functional genomic distribution.

A histone H3.3K36M mutation in mice causes an imbalance of histone modifications and defects in chondrocyte differentiation

Histone lysine-to-methionine (K-to-M) mutations have been identified as driver mutations in human cancers. Interestingly, these 'oncohistone' mutations inhibit the activity of histone methyltransferases. Therefore, they can potentially be used as versatile tools to investigate the roles of histone modifications. In this study, we generated a genetically engineered mouse line in which an H3.3K36M mutation could be induced in the endogenous H3f3b gene. Since H3.3K36M has been identified as a causative mutation of human chondroblastoma, we induced this mutation in the chondrocyte lineage in mouse embryonic limbs. We found that H3.3K36M causes a global reduction in H3K36me2 and defects in chondrocyte differentiation. Importantly, the reduction of H3K36me2 was accompanied by a collapse of normal H3K27me3 distribution. Furthermore, the changes in H3K27me3, especially the loss of H3K27me3 at gene regulatory elements, were associated with the mis-regulated expression of a set of genes important for limb development, including HoxA cluster genes. Thus, through the in vivo induction of the H3.3K36M mutation, we reveal the importance of maintaining the balance between H3K36me2 and H3K27me3 during chondrocyte differentiation and limb development.

Dynamic chromatin accessibility profiling reveals changes in host genome organization in response to baculovirus infection

DNA viruses can hijack and manipulate the host chromatin state to facilitate their infection. Multiple lines of evidences reveal that DNA virus infection results in the host chromatin relocation, yet there is little known about the effects of viral infection on the architecture of host chromatin. Here, a combination of epigenomic, transcriptomic and biochemical assays were conducted to investigate the temporal dynamics of chromatin accessibility in response to Bombyx mori nucleopolyhedrovirus (BmNPV) infection. The high-quality ATAC-seq data indicated that progressive chromatin remodeling took place following BmNPV infection. Viral infection resulted in a more open chromatin architecture, along with the marginalization of host genome and nucleosome disassembly. Moreover, our results revealed that chromatin accessibility in uninfected cells was regulated by euchromatic modifications, whereas the viral-induced highly accessible chromatin regions were originally associated with facultative heterochromatic modification. Overall, our findings illustrate for the first time the organization and accessibility of host chromatin in BmNPV-infected cells, which lay the foundation for future studies on epigenomic regulation mediated by DNA viruses.

As a well-studied arthropod-specific double-stranded DNA virus, Bombyx mori nucleopolyhedrovirus (BmNPV) is a representative member of baculoviruses. BmNPV infection results in significant host chromatin marginalization, which has also been found in most DNA viruses. However, the effects of baculovirus infection on the organization and accessibility of host chromatin are poorly understood. Here, by using ATAC-seq, we show that DNA virus BmNPV infection gradually remodels the accessibility of host chromatin. ATAC-seq data reveal that the marginalized host chromatin is a more accessible architecture along with the depletion of multi-nucleosome depositions. Moreover, our findings suggest the increased accessibility regions are regulated by the facultative heterochromatic modification. Overall, we provide a novel understanding of molecular mechanisms by which baculovirus and DNA viruses alter the organization of host chromatin in epigenomic regulation.