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

HBO1 functions as an epigenetic barrier to hepatocyte plasticity and reprogramming during liver injury

Hepatocytes can reprogram into biliary epithelial cells (BECs) during liver injury, but the underlying epigenetic mechanisms remain poorly understood. Here, we define the chromatin dynamics of this process using single-cell ATAC-seq and identify YAP/TEAD activation as a key driver of chromatin remodeling. An in vivo CRISPR screen highlights the histone acetyltransferase HBO1 as a critical barrier to reprogramming. HBO1 is recruited by YAP to target loci, where it promotes histone H3 lysine 14 acetylation (H3K14ac) and engages the chromatin reader zinc-finger MYND-type containing 8 (ZMYND8) to suppress YAP/TEAD-driven transcription. Loss of HBO1 accelerates chromatin remodeling, enhances YAP binding, and enables a more complete hepatocyte-to-BEC transition. Our findings position HBO1 as an epigenetic brake that restrains YAP-mediated reprogramming, suggesting that targeting HBO1 may enhance hepatocyte plasticity for liver regeneration.

Dynamic regulation of murine RNA polymerase III transcription during heat shock stress

Cells respond to many different types of stresses by overhauling gene expression patterns, both at the transcriptional and translational levels. Under heat stress, global transcription and translation are inhibited, while the expression of chaperone proteins is preferentially favored. As the direct link between mRNA transcription and protein translation, transfer RNA (tRNA) expression is intricately regulated during the stress response. Despite extensive research into the heat shock response (HSR), the regulation of tRNA expression by RNA polymerase III (Pol III) transcription has yet to be fully elucidated in mammalian cells. Here, we examine the regulation of Pol III transcription during different stages of heat shock stress in mouse embryonic stem cells. We observe that Pol III transcription is downregulated after 30 min of heat shock, followed by an overall increase in transcription after 60 min of heat shock. This effect is more evident in tRNAs, although other Pol III gene targets are also similarly affected. Notably, we show that the downregulation at 30 min of heat shock is independent of HSF1, the master transcription factor of the HSR, but that the subsequent increase in expression at 60 min requires HSF1. Taken together, these results demonstrate an adaptive RNA Pol III response to heat stress and an intricate relationship between the canonical HSR and tRNA expression.

NSD2-epigenomic reprogramming and maintenance of plasma cell phenotype in t(4;14) myeloma

Overexpression of the H3K36 histone methyltransferase NSD2 in t(4;14) multiple myeloma (MM) is an early, oncogenic event, and understanding its impact on genomic organisation and expression is relevant to understanding MM biology. We performed epigenetic, transcriptional and phenotypic profiling of the t(4;14) KMS11 myeloma cell line and its isogenic translocation knock out (TKO) to characterise the sequelae of NSD2 overexpression. We found a marked global impact of NSD2 on gene expression and DNA organisation implicating cell identity genes; notably the early lymphocyte regulator, LAIR1 and MM cell surface markers, including CD38, a classical marker of plasma cells which was reduced in TKO cells. Plasma cell transcription factors such as PRDM1, IRF4 and XBP1 were unaffected, suggesting a downstream direct gene effect of NSD2 on cell identity. Changes in cell surface markers suggest an altered surface immunophenotype. Our findings suggest a role for NSD2 in maintaining MM cell identity, with potential implications for future therapeutic strategies based on targeting of NSD2.

Transcriptional stress induces the overexpression of novel lncRNAs that regulate the BRCA1 locus

Long non-coding RNAs (lncRNAs) have been shown to play a role during transcriptional regulation in response to stress. However, their function under stress caused by transcriptional inhibition has not yet been addressed. Using genome-wide assays to elucidate the transcriptional response in human cells caused by RNA polymerase II transcription inhibition, we found three novel regulatory lncRNAs, TILR-1 , TILR-2 , and LINC00910 , that are upregulated as a response to this transcriptional stress. Knockdown experiments showed that the expression of these RNAs is interdependent, and together, they regulate transcription of the nearby BRCA1 locus. The lack of these novel regulatory transcripts also resulted in an increase in cellular proliferation and survival. Public transcriptomic data from different cell lines treated with a variety of transcriptional inhibitors or with heat shock and arsenic stress showed that TILR-1 , TILR-2 , and LINC00910 are commonly upregulated in a broad array of stress conditions. Evolutionary analysis showed that TILR-1 , TILR-2 , and LINC00910 are highly conserved among primates, and their emergence correlates with the duplication of the bidirectional promoter of BRCA1 and NBR1 . We conclude that that coordinate transcription of TILR-1 , TILR-2 , and LINC00910 is stimulated generally by stress and the resulting lncRNAs are novel, functionally-conserved regulators of the BRCA1 locus.

Multiplexed chromatin immunoprecipitation sequencing for quantitative study of histone modifications and chromatin factors

ChIP-seq is a widely used technique for studying histone post-translational modifications and DNA-binding proteins. DNA fragments associated with a specific protein or histone modification epitope are captured by using antibodies, sequenced and mapped to a reference genome. Albeit versatile and popular, performing many parallel ChIP-seq experiments to compare different conditions, replicates and epitopes is laborious, is prone to experimental variation and does not allow quantitative comparisons unless adequate spike-in chromatin is included. We present a detailed protocol for performing and analyzing a multiplexed quantitative chromatin immunoprecipitation-sequencing experiment (MINUTE-ChIP), in which multiple samples are profiled against multiple epitopes in a single workflow. Multiplexing not only dramatically increases the throughput of ChIP-seq experiments (e.g., profiling 12 samples against multiple histone modifications or DNA-binding proteins in a single experiment), but also enables accurate quantitative comparisons. The protocol consists of four parts: sample preparation (i.e., lysis, chromatin fragmentation and barcoding of native or formaldehyde-fixed material), pooling and splitting of the barcoded chromatin into parallel immunoprecipitation reactions, preparation of next-generation sequencing libraries from input and immunoprecipitated DNA and data analysis using our dedicated analysis pipeline. This pipeline autonomously generates quantitatively scaled ChIP-seq tracks for downstream analysis and visualization, alongside necessary quality control indicators. The entire workflow requires basic knowledge in molecular biology and bioinformatics and can be completed in 1 week. MINUTE-ChIP empowers biologists to perform every ChIP-seq experiment with an appropriate number of replicates and control conditions, delivering more statistically robust, exquisitely quantitative and biologically meaningful results.

ERH Enables Early Embryonic Differentiation and Overlays H3K9me3 Heterochromatin on a Cryptic Pluripotency H3K9me3 Landscape in Somatic Cells

Enhancer of Rudimentary Homolog (ERH) is an evolutionarily conserved protein originally characterized in fission yeast 1 and recently shown to maintain H3K9me3 in human fibroblasts 2 . Here, we find that ERH depletion in fibroblasts reverts the H3K9me3 landscape to an embryonic stem cell (ESC) state and enables activation of naïve and pluripotency genes and transposable elements during induced pluripotent stem cell (iPSC) reprogramming. We find that ERH similarly represses totipotent and alternative lineage programs during mouse preimplantation development and is required for proper segregation of the inner cell mass and trophectoderm cell lineages. During human ESC differentiation into germ layer lineages, ERH silences naïve and pluripotency genes, transposable elements, and alternative lineage somatic genes. As in fission yeast, we find that mammalian ERH interacts with RNA-binding proteins to engage and repress its chromatin targets. Our findings reveal a fundamental role for ERH in cell fate specification via the initiation and maintenance of early developmental gene repression.

Epigenetic characterization of adult rhesus monkey spermatogonial stem cells identifies key regulators of stem cell homeostasis

Spermatogonial stem cells (SSCs) play crucial roles in the preservation of male fertility. However, successful ex vivo expansion of authentic human SSCs remains elusive due to the inadequate understanding of SSC homeostasis regulation. Using rhesus monkeys (Macaca mulatta) as a representative model, we characterized SSCs and progenitor subsets through single-cell RNA sequencing using a cell-specific network approach. We also profiled chromatin status and major histone modifications (H3K4me1, H3K4me3, H3K27ac, H3K27me3 and H3K9me3), and subsequently mapped promoters and active enhancers in TSPAN33+ putative SSCs. Comparing the epigenetic changes between fresh TSPAN33+ cells and cultured TSPAN33+ cells (resembling progenitors), we identified the regulatory elements with higher activity in SSCs, and the potential transcription factors and signaling pathways implicated in SSC regulation. Specifically, TGF-β signaling is activated in monkey putative SSCs. We provided evidence supporting its role in promoting self-renewal of monkey SSCs in culture. Overall, this study outlines the epigenetic landscapes of monkey SSCs and provides clues for optimization of the culture condition for primate SSCs expansion.

Genome-wide mapping of native co-localized G4s and R-loops in living cells

The interplay between G4s and R-loops are emerging in regulating DNA repair, replication, and transcription. A comprehensive picture of native co-localized G4s and R-loops in living cells is currently lacking. Here, we describe the development of HepG4-seq and an optimized HBD-seq methods, which robustly capture native G4s and R-loops, respectively, in living cells. We successfully employed these methods to establish comprehensive maps of native co-localized G4s and R-loops in human HEK293 cells and mouse embryonic stem cells (mESCs). We discovered that co-localized G4s and R-loops are dynamically altered in a cell type-dependent manner and are largely localized at active promoters and enhancers of transcriptional active genes. We further demonstrated the helicase Dhx9 as a direct and major regulator that modulates the formation and resolution of co-localized G4s and R-loops. Depletion of Dhx9 impaired the self-renewal and differentiation capacities of mESCs by altering the transcription of co-localized G4s and R-loops -associated genes. Taken together, our work established that the endogenous co-localized G4s and R-loops are prevalently persisted in the regulatory regions of active genes and are involved in the transcriptional regulation of their linked genes, opening the door for exploring broader roles of co-localized G4s and R-loops in development and disease.

Capture of RNA-binding proteins across mouse tissues using HARD-AP

RNA-binding proteins (RBPs) modulate all aspects of RNA metabolism, but a comprehensive picture of RBP expression across tissues is lacking. Here, we describe our development of the method we call HARD-AP that robustly retrieves RBPs and tightly associated RNA regulatory complexes from cultured cells and fresh tissues. We successfully use HARD-AP to establish a comprehensive atlas of RBPs across mouse primary organs. We then systematically map RNA-binding sites of these RBPs using machine learning-based modeling. Notably, the modeling reveals that the LIM domain as an RNA-binding domain in many RBPs. We validate the LIM-domain-only protein Csrp1 as a tissue-dependent RNA binding protein. Taken together, HARD-AP is a powerful approach that can be used to identify RBPomes from any type of sample, allowing comprehensive and physiologically relevant networks of RNA-protein interactions.

Evaluation of barley genotypes for drought adaptability: based on stress indices and comprehensive evaluation as criteria

The prevalence of drought events worldwide emphasizes the importance of screening and cultivating drought-adapted crops. In this study, 206 germplasm resources were used as materials, dry weight as target trait, and two genotyping methods as criteria to evaluate drought adaptability at the seedling establishment stage. The results showed a significant decrease in average dry weight of the tested germplasm resources (from 746.90 mg to 285.40 mg) and rich variation in the responses of dry weight among each genotype to drought (CV=61.14%). In traditional evaluation method, drought resistance coefficient (DC), geometric mean productivity index (GMP), mean productivity index (MP), stress susceptibility index (SSI), stress tolerance index (STI), and tolerance index (TOL) also exhibited diversity in tested genotypes (CV>30%). However, these indices showed varying degrees of explanation for dry weight under stress and non-stress environments and failed to differentiate drought adaptability among genotypes clearly. In new evaluation method, four stress indices were developed to quantify barley seedling production and stability capacities. Compared to traditional stress indices, the stress production index (SI) explained dry weight more comprehensively under stress conditions (R2 = 0.98), while the ideal production index (II) explained dry weight better under non-stress conditions (R2 = 0.89). Furthermore, the potential index (PI) and elasticity index (EI) eliminated disparities in traditional stress indices and comprehensively clarified the contribution of elasticity and potential to production capacity under drought stress. Ultimately, through grading evaluation and cluster analysis, the tested germplasm resources were effectively categorized, and 11 genotypes were identified as suitable for cultivation in arid areas. Overall, the comprehensive evaluation method based on the newly developed stress indices surpasses the traditional method in screening drought adaptability of crops and serves as a vital tool for identifying high-stability and high-production capacities genotypes in various environments, which is expected to provide practical guidance for barley planting and breeding in arid areas.

Oestrogen receptor positive breast cancer and its embedded mechanism: breast cancer resistance to conventional drugs and related therapies, a review

Traditional medication and alternative therapies have long been used to treat breast cancer. One of the main problems with current treatments is that there is an increase in drug resistance in the cancer cells owing to genetic differences such as mutational changes, epigenetic changes and miRNA (microRNA) alterations such as miR-1246, miR-298, miR-27b and miR-33a, along with epigenetic modifications, such as Histone3 acetylation and CCCTC-Binding Factor (CTCF) hypermethylation for drug resistance in breast cancer cell lines. Certain forms of conventional drug resistance have been linked to genetic changes in genes such as ABCB1, AKT, S100A8/A9, TAGLN2 and NPM. This review aims to explore the current approaches to counter breast cancer, the action mechanism, along with novel therapeutic methods endowing potential drug resistance. The investigation of novel therapeutic approaches sheds light on the phenomenon of drug resistance including genetic variations that impact distinct forms of oestrogen receptor (ER) cancer, genetic changes, epigenetics-reported resistance and their identification in patients. Long-term effective therapy for breast cancer includes selective oestrogen receptor modulators, selective oestrogen receptor degraders and genetic variations, such as mutations in nuclear genes, epigenetic modifications and miRNA alterations in target proteins. Novel research addressing combinational therapies including maytansine, photodynamic therapy, guajadiol, talazoparib, COX2 inhibitors and miRNA 1246 inhibitors have been developed to improve patient survival rates.

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.