Quantifying ChIP-seq data: a spiking method providing an internal reference for sample-to-sample normalization

SpikeFlow: automated and flexible analysis of ChIP-Seq data with spike-in control

ChIP with reference exogenous genome (ChIP-Rx) is widely used to study histone modification changes across different biological conditions. A key step in the bioinformatics analysis of this data is calculating the normalization factors, which vary from the standard ChIP-seq pipelines. Choosing and applying the appropriate normalization method is crucial for interpreting the biological results. However, a comprehensive pipeline for complete ChIP-Rx data analysis is lacking. To address these challenges, we introduce SpikeFlow, an integrated Snakemake workflow that combines features from various existing tools to streamline ChIP-Rx data processing and enhance usability. SpikeFlow automates spike-in data scaling and provides multiple normalization options. It also performs peak calling and differential analysis with distinct modalities, enabling the detection of enrichment regions for histone modifications and transcription factor binding. Our workflow runs in-depth quality control at all the processing steps and generates an analysis report with tables and graphs to facilitate results interpretation. We validated the pipeline by performing a comparative analysis with DiffBind and SpikChIP, demonstrating robust performances in various biological models. By combining diverse functionalities into a single platform, SpikeFlow aims to simplify ChIP-Rx data analysis for the research community.

L3MBTL2 maintains MYCN-amplified neuroblastoma cell proliferation through silencing NRIP3 and BRME1 genes

Epigenetic alterations critically affect tumor development. Polycomb-group complexes constitute an evolutionarily conserved epigenetic machinery that regulates stem cell fate and development. They are implicated in tumorigenesis, primarily via histone modification. Polycomb repressive complex 1 (PRC1) complexes 1-6 (PRC1.1-6) mediate the ubiquitination of histone H2A on lysine 119 (H2AK119ub). Here, we studied the functional roles of a PRC1.6 molecule, L3MBTL2, in neuroblastoma (NB) cells. L3MBTL2-knockout and knockdown revealed that L3MBTL2 depletion suppressed NB cell proliferation via cell-cycle arrest and gamma-H2A.X upregulation. L3MBTL2-knockout profoundly suppressed xenograft tumor formation. Transcriptome analysis revealed suppressed cell-cycle-related and activated differentiation-related pathways. Break repair meiotic recombinase recruitment factor 1 (BRME1) and nuclear receptor interacting protein 3 (NRIP3) were notably de-repressed by L3MBTL2-knockout. The deletion of L3MBTL2 reduced enrichment of H2AK119ub and PCGF6 at transcriptional start site proximal regions of the targets. Add-back studies unveiled the importance of L3MBTL2-BRME1 and -NRIP3 axes for NB cell proliferation. We further manifested the association of MYCN with de-repression of NRIP3 in an L3MBTL2-deficient context. Therefore, this study first revealed the significance of L3MBTL2-mediated gene silencing in MYCN-amplified NB cells.

Highly quantitative measurement of differential protein-genome binding with PerCell chromatin sequencing

We report a universal strategy for 2D chromatin sequencing, to increase uniform data analyses and sharing across labs, and to facilitate highly quantitative comparisons across experimental conditions. Within our system, we provide wetlab and drylab tools for researchers to establish and analyze protein-genome binding data with PerCell ChIP-seq. Our methodology is virtually no cost and flexible, enabling rapid, quantitative, internally normalized chromatin sequencing to catalyze project development in a variety of systems, including in vivo zebrafish epigenomics and cancer cell epigenomics. While we highlight utility in these key areas, our methodology is flexible enough such that rapid comparisons of cellular spike-in versus non spike-in are possible, and generalizability to nuclease-based 2D chromatin sequencing would also be possible within the framework of our pipeline. Through the use of well-defined cellular ratios containing orthologous species' chromatin, we enable cross-species comparative epigenomics and highly quantitative low-cost chromatin sequencing with utility across a range of disciplines.

A Constitutive Heterochromatic Region Shapes Genome Organization and Impacts Gene Expression in Neurospora crassa

Genome organization is essential for proper function, including gene expression. In metazoan genome organization, chromatin loops and Topologically Associated Domains (TADs) facilitate local gene clustering, while chromosomes form distinct nuclear territories characterized by compartmentalization of silent heterochromatin at the nuclear periphery and active euchromatin in the nucleus center. A similar hierarchical organization occurs in the fungus Neurospora crassa where its seven chromosomes form a Rabl conformation, where heterochromatic centromeres and telomeres independently cluster at the nuclear membrane, while interspersed heterochromatic loci in Neurospora aggregate across megabases of linear genomic distance for forming TAD-like structures. However, the role of individual heterochromatic loci in normal genome organization and function is unknown. Here, we examined the genome organization of a Neurospora strain harboring a ∼47.4 kilobase facultative (temporarily silent) heterochromatic region deletion, as well as the genome organization of a strain deleted of a 110.6 kilobase permanently silent constitutive heterochromatic region. While the facultative heterochromatin deletion had little effect on local chromatin structure, the constitutive heterochromatin deletion altered local TAD-like structures, gene expression, and the predicted 3D genome structure by qualitatively repositioning genes into the nucleus center. Our work elucidates the role of individual heterochromatic regions for genome organization and function.

Genome-wide ATAC-see screening identifies TFDP1 as a modulator of global chromatin accessibility

Chromatin accessibility is a hallmark of active regulatory regions and is functionally linked to transcriptional networks and cell identity. However, the molecular mechanisms and networks that govern chromatin accessibility have not been thoroughly studied. Here we conducted a genome-wide CRISPR screening combined with an optimized ATAC-see protocol to identify genes that modulate global chromatin accessibility. In addition to known chromatin regulators like CREBBP and EP400, we discovered a number of previously unrecognized proteins that modulate chromatin accessibility, including TFDP1, HNRNPU, EIF3D and THAP11 belonging to diverse biological pathways. ATAC-seq analysis upon their knockouts revealed their distinct and specific effects on chromatin accessibility. Remarkably, we found that TFDP1, a transcription factor, modulates global chromatin accessibility through transcriptional regulation of canonical histones. In addition, our findings highlight the manipulation of chromatin accessibility as an approach to enhance various cell engineering applications, including genome editing and induced pluripotent stem cell reprogramming.

Quantitative transcriptomic and epigenomic data analysis: a primer

The advent of microarray and second generation sequencing technology has revolutionized the field of molecular biology, allowing researchers to quantitatively assess transcriptomic and epigenomic features in a comprehensive and cost-efficient manner. Moreover, technical advancements have pushed the resolution of these sequencing techniques to the single cell level. As a result, the bottleneck of molecular biology research has shifted from the bench to the subsequent omics data analysis. Even though most methodologies share the same general strategy, state-of-the-art literature typically focuses on data type specific approaches and already assumes expert knowledge. Here, however, we aim at providing conceptual insight in the principles of genome-wide quantitative transcriptomic and epigenomic (including open chromatin assay) data analysis by describing a generic workflow. By starting from a general framework and its assumptions, the need for alternative or additional data-analytical solutions when working with specific data types becomes clear, and are hence introduced. Thus, we aim to enable readers with basic omics expertise to deepen their conceptual and statistical understanding of general strategies and pitfalls in omics data analysis and to facilitate subsequent progression to more specialized literature.

The CUT&RUN greenlist: genomic regions of consistent noise are effective normalizing factors for quantitative epigenome mapping

Cleavage Under Targets and Release Using Nuclease (CUT&RUN) is a recent development for epigenome mapping, but its unique methodology can hamper proper quantitative analyses. As traditional normalization approaches have been shown to be inaccurate, we sought to determine endogenous normalization factors based on the human genome regions of constant nonspecific signal. This constancy was determined by applying Shannon's information entropy, and the set of normalizer regions, which we named the 'Greenlist', was extensively validated using publicly available datasets. We demonstrate here that the greenlist normalization outperforms the current top standards, and remains consistent across different experimental setups, cell lines and antibodies; the approach can even be applied to different species or to CUT&Tag. Requiring no additional experimental steps and no added cost, this approach can be universally applied to CUT&RUN experiments to greatly minimize the interference of technical variation over the biological epigenome changes of interest.

Selective gene expression maintains human tRNA anticodon pools during differentiation

Transfer RNAs are essential for translating genetic information into proteins. The human genome contains hundreds of predicted tRNA genes, many in multiple copies. How their expression is regulated to control tRNA repertoires is unknown. Here we combined quantitative tRNA profiling and chromatin immunoprecipitation with sequencing to measure tRNA expression following the differentiation of human induced pluripotent stem cells into neuronal and cardiac cells. We find that tRNA transcript levels vary substantially, whereas tRNA anticodon pools, which govern decoding rates, are more stable among cell types. Mechanistically, RNA polymerase III transcribes a wide range of tRNA genes in human induced pluripotent stem cells but on differentiation becomes constrained to a subset we define as housekeeping tRNAs. This shift is mediated by decreased mTORC1 signalling, which activates the RNA polymerase III repressor MAF1. Our data explain how tRNA anticodon pools are buffered to maintain decoding speed across cell types and reveal that mTORC1 drives selective tRNA expression during differentiation.

Loss of cohesin regulator PDS5A reveals repressive role of Polycomb loops

Polycomb Repressive Complexes 1 and 2 (PRC1, PRC2) are conserved epigenetic regulators that promote transcriptional gene silencing. PRC1 and PRC2 converge on shared targets, catalyzing repressive histone modifications. Additionally, a subset of PRC1/PRC2 targets engage in long-range interactions whose functions in gene silencing are poorly understood. Using a CRISPR screen in mouse embryonic stem cells, we found that the cohesin regulator PDS5A links transcriptional silencing by Polycomb and 3D genome organization. PDS5A deletion impairs cohesin unloading and results in derepression of a subset of endogenous PRC1/PRC2 target genes. Importantly, derepression is not linked to loss of Polycomb chromatin domains. Instead, PDS5A removal causes aberrant cohesin activity leading to ectopic insulation sites, which disrupt the formation of ultra-long Polycomb loops. We show that these loops are important for robust silencing at a subset of PRC1/PRC2 target genes and that maintenance of cohesin-dependent genome architecture is critical for Polycomb regulation.

Contrasting effects of whole-body and hepatocyte-specific deletion of the RNA polymerase III repressor Maf1 in the mouse

MAF1 is a nutrient-sensitive, TORC1-regulated repressor of RNA polymerase III (Pol III). MAF1 downregulation leads to increased lipogenesis in Drosophila melanogaster, Caenorhabditis elegans, and mice. However, Maf1 -/- mice are lean as increased lipogenesis is counterbalanced by futile pre-tRNA synthesis and degradation, resulting in increased energy expenditure. We compared Chow-fed Maf1 -/- mice with Chow- or High Fat (HF)-fed Maf1 hep-/- mice that lack MAF1 specifically in hepatocytes. Unlike Maf1 -/- mice, Maf1 hep-/- mice become heavier and fattier than control mice with old age and much earlier under a HF diet. Liver ChIPseq, RNAseq and proteomics analyses indicate increased Pol III occupancy at Pol III genes, very few differences in mRNA accumulation, and protein accumulation changes consistent with increased lipogenesis. Futile pre-tRNA synthesis and degradation in the liver, as likely occurs in Maf1 hep-/- mice, thus seems insufficient to counteract increased lipogenesis. Indeed, RNAseq and metabolite profiling indicate that liver phenotypes of Maf1 -/- mice are strongly influenced by systemic inter-organ communication. Among common changes in the three phenotypically distinct cohorts, Angiogenin downregulation is likely linked to increased Pol III occupancy of tRNA genes in the Angiogenin promoter.

Histone deacetylation and cytosine methylation compartmentalize heterochromatic regions in the genome organization of Neurospora crassa

Chromosomes must correctly fold in eukaryotic nuclei for proper genome function. Eukaryotic organisms hierarchically organize their genomes, including in the fungus Neurospora crassa, where chromatin fiber loops compact into Topologically Associated Domain-like structures formed by heterochromatic region aggregation. However, insufficient data exist on how histone posttranslational modifications (PTMs), including acetylation, affect genome organization. In Neurospora, the HCHC complex [composed of the proteins HDA-1, CDP-2 (Chromodomain Protein-2), Heterochromatin Protein-1, and CHAP (CDP-2 and HDA-1 Associated Protein)] deacetylates heterochromatic nucleosomes, as loss of individual HCHC members increases centromeric acetylation, and alters the methylation of cytosines in DNA. Here, we assess whether the HCHC complex affects genome organization by performing Hi-C in strains deleted of the cdp-2 or chap genes. CDP-2 loss increases intra- and interchromosomal heterochromatic region interactions, while loss of CHAP decreases heterochromatic region compaction. Individual HCHC mutants exhibit different patterns of histone PTMs genome-wide, as CDP-2 deletion increases heterochromatic H4K16 acetylation, yet smaller heterochromatic regions lose H3K9 trimethylation and gain interheterochromatic region interactions; CHAP loss produces minimal acetylation changes but increases heterochromatic H3K9me3 enrichment. Loss of both CDP-2 and the DIM-2 DNA methyltransferase causes extensive genome disorder as heterochromatic-euchromatic contacts increase despite additional H3K9me3 enrichment. Our results highlight how the increased cytosine methylation in HCHC mutants ensures genome compartmentalization when heterochromatic regions become hyperacetylated without HDAC activity.

Functional dissection of PRC1 subunits RYBP and YAF2 during neural differentiation of embryonic stem cells

Polycomb repressive complex 1 (PRC1) comprises two different complexes: CBX-containing canonical PRC1 (cPRC1) and RYBP/YAF2-containing variant PRC1 (vPRC1). RYBP-vPRC1 or YAF2-vPRC1 catalyzes H2AK119ub through a positive-feedback model; however, whether RYBP and YAF2 have different regulatory functions is still unclear. Here, we show that the expression of RYBP and YAF2 decreases and increases, respectively, during neural differentiation of embryonic stem cells (ESCs). Rybp knockout impairs neural differentiation by activating Wnt signaling and derepressing nonneuroectoderm-associated genes. However, Yaf2 knockout promotes neural differentiation and leads to redistribution of RYBP binding, increases enrichment of RYBP and H2AK119ub on the RYBP-YAF2 cotargeted genes, and prevents ectopic derepression of nonneuroectoderm-associated genes in neural-differentiated cells. Taken together, this study reveals that RYBP and YAF2 function differentially in regulating mESC neural differentiation.

MYC function and regulation in physiological perspective

MYC, a key member of the Myc-proto-oncogene family, is a universal transcription amplifier that regulates almost every physiological process in a cell including cell cycle, proliferation, metabolism, differentiation, and apoptosis. MYC interacts with several cofactors, chromatin modifiers, and regulators to direct gene expression. MYC levels are tightly regulated, and deregulation of MYC has been associated with numerous diseases including cancer. Understanding the comprehensive biology of MYC under physiological conditions is an utmost necessity to demark biological functions of MYC from its pathological functions. Here we review the recent advances in biological mechanisms, functions, and regulation of MYC. We also emphasize the role of MYC as a global transcription amplifier.

Histone deacetylation and cytosine methylation compartmentalize heterochromatic regions in the genome organization of Neurospora crassa

Chromosomes must correctly fold in eukaryotic nuclei for proper genome function. Eukaryotic organisms hierarchically organize their genomes, including in the fungus Neurospora crassa, where chromatin fiber loops compact into Topologically Associated Domain (TAD)-like structures formed by heterochromatic region aggregation. However, insufficient data exists on how histone post-translational modifications, including acetylation, affect genome organization. In Neurospora, the HCHC complex (comprised of the proteins HDA-1, CDP-2, HP1, and CHAP) deacetylates heterochromatic nucleosomes, as loss of individual HCHC members increases centromeric acetylation and alters the methylation of cytosines in DNA. Here, we assess if the HCHC complex affects genome organization by performing Hi-C in strains deleted of the cdp-2 or chap genes. CDP-2 loss increases intra- and inter-chromosomal heterochromatic region interactions, while loss of CHAP decreases heterochromatic region compaction. Individual HCHC mutants exhibit different patterns of histone post-translational modifications genome-wide: without CDP-2, heterochromatic H4K16 acetylation is increased, yet smaller heterochromatic regions lose H3K9 trimethylation and gain inter-heterochromatic region interactions; CHAP loss produces minimal acetylation changes but increases heterochromatic H3K9me3 enrichment. Loss of both CDP-2 and the DIM-2 DNA methyltransferase causes extensive genome disorder, as heterochromatic-euchromatic contacts increase despite additional H3K9me3 enrichment. Our results highlight how the increased cytosine methylation in HCHC mutants ensures genome compartmentalization when heterochromatic regions become hyperacetylated without HDAC activity.

Globally elevated levels of histone H3 lysine 9 trimethylation in early infancy are associated with poor growth trajectory in Bangladeshi children

BACKGROUND

Stunting is a global health problem affecting hundreds of millions of children worldwide and contributing to 45% of deaths in children under the age of five. Current therapeutic interventions have limited efficacy. Understanding the epigenetic changes underlying stunting will elucidate molecular mechanisms and likely lead to new therapies.

RESULTS

We profiled the repressive mark histone H3 lysine 9 trimethylation (H3K9me3) genome-wide in peripheral blood mononuclear cells (PBMCs) from 18-week-old infants (n = 15) and mothers (n = 14) enrolled in the PROVIDE study established in an urban slum in Bangladesh. We associated H3K9me3 levels within individual loci as well as genome-wide with anthropometric measurements and other biomarkers of stunting and performed functional annotation of differentially affected regions. Despite the relatively small number of samples from this vulnerable population, we observed globally elevated H3K9me3 levels were associated with poor linear growth between birth and one year of age. A large proportion of the differentially methylated genes code for proteins targeting viral mRNA and highly significant regions were enriched in transposon elements with potential regulatory roles in immune system activation and cytokine production. Maternal data show a similar trend with child's anthropometry; however, these trends lack statistical significance to infer an intergenerational relationship.

CONCLUSIONS

We speculate that high H3K9me3 levels may result in poor linear growth by repressing genes involved in immune system activation. Importantly, changes to H3K9me3 were detectable before the overt manifestation of stunting and therefore may be valuable as new biomarkers of stunting.

Streamlined quantitative analysis of histone modification abundance at nucleosome-scale resolution with siQ-ChIP version 2.0

We recently introduced an absolute and physical quantitative scale for chromatin immunoprecipitation followed by sequencing (ChIP-seq). The scale itself was determined directly from measurements routinely made on sequencing samples without additional reagents or spike-ins. We called this approach sans spike-in quantitative ChIP, or siQ-ChIP. Herein, we extend those results in several ways. First, we simplified the calculations defining the quantitative scale, reducing practitioner burden. Second, we reveal a normalization constraint implied by the quantitative scale and introduce a new scheme for generating 'tracks'. The constraint requires that tracks are probability distributions so that quantified ChIP-seq is analogous to a mass distribution. Third, we introduce some whole-genome analyses that allow us, for example, to project the IP mass (immunoprecipitated mass) onto the genome to evaluate how much of any genomic interval was captured in the IP. We applied siQ-ChIP to p300/CBP inhibition and compare our results to those of others. We detail how the same data-level observations are misinterpreted in the literature when tracks are not understood as probability densities and are compared without correct quantitative scaling, and we offer new interpretations of p300/CBP inhibition outcomes.

The ENCODE Imputation Challenge: a critical assessment of methods for cross-cell type imputation of epigenomic profiles

A promising alternative to comprehensively performing genomics experiments is to, instead, perform a subset of experiments and use computational methods to impute the remainder. However, identifying the best imputation methods and what measures meaningfully evaluate performance are open questions. We address these questions by comprehensively analyzing 23 methods from the ENCODE Imputation Challenge. We find that imputation evaluations are challenging and confounded by distributional shifts from differences in data collection and processing over time, the amount of available data, and redundancy among performance measures. Our analyses suggest simple steps for overcoming these issues and promising directions for more robust research.

Incomplete erasure of histone marks during epigenetic reprogramming in medaka early development

Epigenetic modifications undergo drastic erasure and reestablishment after fertilization. This reprogramming is required for proper embryonic development and cell differentiation. In mammals, some histone modifications are not completely reprogrammed and play critical roles in later development. In contrast, in nonmammalian vertebrates, most histone modifications are thought to be more intensively erased and reestablished by the stage of zygotic genome activation (ZGA). However, histone modifications that escape reprogramming in nonmammalian vertebrates and their potential functional roles remain unknown. Here, we quantitatively and comprehensively analyzed histone modification dynamics during epigenetic reprogramming in Japanese killifish, medaka (Oryzias latipes) embryos. Our data revealed that H3K27ac, H3K27me3, and H3K9me3 escape complete reprogramming, whereas H3K4 methylation is completely erased during cleavage stage. Furthermore, we experimentally showed the functional roles of such retained modifications at early stages: (i) H3K27ac premarks promoters during the cleavage stage, and inhibition of histone acetyltransferases disrupts proper patterning of H3K4 and H3K27 methylation at CpG-dense promoters, but does not affect chromatin accessibility after ZGA; (ii) H3K9me3 is globally erased but specifically retained at telomeric regions, which is required for maintenance of genomic stability during the cleavage stage. These results expand the understanding of diversity and conservation of reprogramming in vertebrates, and unveil previously uncharacterized functions of histone modifications retained during epigenetic reprogramming.

Quantitative Comparison of Multiple Chromatin Immunoprecipitation-Sequencing (ChIP-seq) Experiments with spikChIP

The chromatin immunoprecipitation coupled with the next-generation sequencing (ChIP-seq) is a powerful technique that enables to characterize the genomic distribution of chromatin-associated proteins, histone posttranslational modifications, and histone variants. However, in the absence of a reference control for monitoring experimental and biological variations, the standard ChIP-seq scheme is unable to accurately assess changes in the abundance of chromatin targets across different experimental samples. To overcome this limitation, the combination of external spike-in material with the experimental chromatin is offered as an effective solution for quantitative comparison of ChIP-seq data across different conditions. Here, we detail (i) the experimental protocol for preparing quality control spike-in chromatin from Drosophila melanogaster cells and (ii) the computational protocol to compare ChIP-seq samples with spike-in based on the use of the spikChIP software.

CBP Is Required for Establishing Adaptive Gene Programs in the Adult Mouse Brain

Environmental factors and life experiences impinge on brain circuits triggering adaptive changes. Epigenetic regulators contribute to this neuroadaptation by enhancing or suppressing specific gene programs. The paralogous transcriptional coactivators and lysine acetyltransferases CREB binding protein (CBP) and p300 are involved in brain plasticity and stimulus-dependent transcription, but their specific roles in neuroadaptation are not fully understood. Here we investigated the impact of eliminating either CBP or p300 in excitatory neurons of the adult forebrain of mice from both sexes using inducible and cell type-restricted knock-out strains. The elimination of CBP, but not p300, reduced the expression and chromatin acetylation of plasticity genes, dampened activity-driven transcription, and caused memory deficits. The defects became more prominent in elderly mice and in paradigms that involved enduring changes in transcription, such as kindling and environmental enrichment, in which CBP loss interfered with the establishment of activity-induced transcriptional and epigenetic changes in response to stimulus or experience. These findings further strengthen the link between CBP deficiency in excitatory neurons and etiopathology in the nervous system.SIGNIFICANCE STATEMENT How environmental conditions and life experiences impinge on mature brain circuits to elicit adaptive responses that favor the survival of the organism remains an outstanding question in neurosciences. Epigenetic regulators are thought to contribute to neuroadaptation by initiating or enhancing adaptive gene programs. In this article, we examined the role of CREB binding protein (CBP) and p300, two paralogous transcriptional coactivators and histone acetyltransferases involved in cognitive processes and intellectual disability, in neuroadaptation in adult hippocampal circuits. Our experiments demonstrate that CBP, but not its paralog p300, plays a highly specific role in the epigenetic regulation of neuronal plasticity gene programs in response to stimulus and provide unprecedented insight into the molecular mechanisms underlying neuroadaptation.

DNA damage-induced transcription stress triggers the genome-wide degradation of promoter-bound Pol II

The precise regulation of RNA Polymerase II (Pol II) transcription after genotoxic stress is crucial for proper execution of the DNA damage-induced stress response. While stalling of Pol II on transcription-blocking lesions (TBLs) blocks transcript elongation and initiates DNA repair in cis, TBLs additionally elicit a response in trans that regulates transcription genome-wide. Here we uncover that, after an initial elongation block in cis, TBLs trigger the genome-wide VCP-mediated proteasomal degradation of promoter-bound, P-Ser5-modified Pol II in trans. This degradation is mechanistically distinct from processing of TBL-stalled Pol II, is signaled via GSK3, and contributes to the TBL-induced transcription block, even in transcription-coupled repair-deficient cells. Thus, our data reveal the targeted degradation of promoter-bound Pol II as a critical pathway that allows cells to cope with DNA damage-induced transcription stress and enables the genome-wide adaptation of transcription to genotoxic stress.

Cohesin-dependent chromosome loop extrusion is limited by transcription and stalled replication forks

Genome function depends on regulated chromosome folding, and loop extrusion by the protein complex cohesin is essential for this multilayered organization. The chromosomal positioning of cohesin is controlled by transcription, and the complex also localizes to stalled replication forks. However, the role of transcription and replication in chromosome looping remains unclear. Here, we show that reduction of chromosome-bound RNA polymerase weakens normal cohesin loop extrusion boundaries, allowing cohesin to form new long-range chromosome cis interactions. Stress response genes induced by transcription inhibition are also shown to act as new loop extrusion boundaries. Furthermore, cohesin loop extrusion during early S phase is jointly controlled by transcription and replication units. Together, the results reveal that replication and transcription machineries are chromosome-folding regulators that block the progression of loop-extruding cohesin, opening for new perspectives on cohesin's roles in genome function and stability.

CTCF: a misguided jack-of-all-trades in cancer cells

The emergence and progression of cancers is accompanied by a dysregulation of transcriptional programs. The three-dimensional (3D) organization of the human genome has emerged as an important multi-level mediator of gene transcription and regulation. In cancer cells, this organization can be restructured, providing a framework for the deregulation of gene activity. The CTCF protein, initially identified as the product from a tumor suppressor gene, is a jack-of-all-trades for the formation of 3D genome organization in normal cells. Here, we summarize how CTCF is involved in the multi-level organization of the human genome and we discuss emerging insights into how perturbed CTCF function and DNA binding causes the activation of oncogenes in cancer cells, mostly through a process of enhancer hijacking. Moreover, we highlight non-canonical functions of CTCF that can be relevant for the emergence of cancers as well. Finally, we provide guidelines for the computational identification of perturbed CTCF binding and reorganized 3D genome structure in cancer cells.

Prediction of histone post-translational modification patterns based on nascent transcription data

The role of histone modifications in transcription remains incompletely understood. Here, we examine the relationship between histone modifications and transcription using experimental perturbations combined with sensitive machine-learning tools. Transcription predicted the variation in active histone marks and complex chromatin states, like bivalent promoters, down to single-nucleosome resolution and at an accuracy that rivaled the correspondence between independent ChIP-seq experiments. Blocking transcription rapidly removed two punctate marks, H3K4me3 and H3K27ac, from chromatin indicating that transcription is required for active histone modifications. Transcription was also required for maintenance of H3K27me3, consistent with a role for RNA in recruiting PRC2. A subset of DNase-I-hypersensitive sites were refractory to prediction, precluding models where transcription initiates pervasively at any open chromatin. Our results, in combination with past literature, support a model in which active histone modifications serve a supportive, rather than an essential regulatory, role in transcription.

Bioinformatics Methods for ChIP-seq Histone Analysis

The field of genomics and genome-wide analysis has exploded since around 2008 with the development of high-throughput omics approaches, largely driven by the emergence of the next-generation sequencing technologies. Among the different biological applications supported by recent sequencing technologies, ChIP-seq (Chromatin ImmunoPrecipitation followed by Sequencing) is one of the most powerful techniques which has dramatically changed our view of the epigenetics landscape of cells.In this chapter, I will present and discuss the main steps of bioinformatic and biostatistical analysis of ChIP-seq data (Fig. 1). While this technique has been widely used to study transcription factor binding sites, I will focus here on the analysis of histone modifications.

Productive visualization of high-throughput sequencing data using the SeqCode open portable platform

Large-scale sequencing techniques to chart genomes are entirely consolidated. Stable computational methods to perform primary tasks such as quality control, read mapping, peak calling, and counting are likewise available. However, there is a lack of uniform standards for graphical data mining, which is also of central importance. To fill this gap, we developed SeqCode, an open suite of applications that analyzes sequencing data in an elegant but efficient manner. Our software is a portable resource written in ANSI C that can be expected to work for almost all genomes in any computational configuration. Furthermore, we offer a user-friendly front-end web server that integrates SeqCode functions with other graphical analysis tools. Our analysis and visualization toolkit represents a significant improvement in terms of performance and usability as compare to other existing programs. Thus, SeqCode has the potential to become a key multipurpose instrument for high-throughput professional analysis; further, it provides an extremely useful open educational platform for the world-wide scientific community. SeqCode website is hosted at http://ldicrocelab.crg.eu, and the source code is freely distributed at https://github.com/eblancoga/seqcode.

PRC1 drives Polycomb-mediated gene repression by controlling transcription initiation and burst frequency

The Polycomb repressive system plays a fundamental role in controlling gene expression during mammalian development. To achieve this, Polycomb repressive complexes 1 and 2 (PRC1 and PRC2) bind target genes and use histone modification-dependent feedback mechanisms to form Polycomb chromatin domains and repress transcription. The inter-relatedness of PRC1 and PRC2 activity at these sites has made it difficult to discover the specific components of Polycomb chromatin domains that drive gene repression and to understand mechanistically how this is achieved. Here, by exploiting rapid degron-based approaches and time-resolved genomics, we kinetically dissect Polycomb-mediated repression and discover that PRC1 functions independently of PRC2 to counteract RNA polymerase II binding and transcription initiation. Using single-cell gene expression analysis, we reveal that PRC1 acts uniformly within the cell population and that repression is achieved by controlling transcriptional burst frequency. These important new discoveries provide a mechanistic and conceptual framework for Polycomb-dependent transcriptional control.

Histone H3 lysine 27 acetylation profile undergoes two global shifts in undernourished children and suggests altered one-carbon metabolism

BACKGROUND

Stunting is a condition in which a child does not reach their full growth potential due to chronic undernutrition. It arises during the first 2 years of a child's life and is associated with developmental deficiencies and life-long health problems. Current interventions provide some benefit, but new approaches to prevention and treatment grounded in a molecular understanding of stunting are needed. Epigenetic analyses are critical as they can provide insight into how signals from a poor environment lead to changes in cell function.

RESULTS

Here we profiled histone H3 acetylation on lysine 27 (H3K27ac) in peripheral blood mononuclear cells (PBMCs) of 18-week-old (n = 14) and 1-year-old children (n = 22) living in an urban slum in Dhaka, Bangladesh. We show that 18-week-old children destined to become stunted have elevated levels of H3K27ac overall, functional analysis of which indicates activation of the immune system and stress response pathways as a primary response to a poor environment with high pathogen load. Conversely, overt stunting at 1-year-of age is associated with globally reduced H3K27ac that is indicative of metabolic rewiring and downregulation of the immune system and DNA repair pathways that are likely secondary responses to chronic exposure to a poor environment with limited nutrients. Among processes altered in 1-year-old children, we identified one-carbon metabolism, the significance of which is supported by integrative analysis with results from histone H3 trimethylation on lysine 4 (H3K4me3). Together, these results suggest altered one-carbon metabolism in this population of stunted children.

CONCLUSIONS

The epigenomes of stunted children undergo two global changes in H3K27ac within their first year of life, which are associated with probable initial hyperactive immune responses followed by reduced metabolic capacity. Limitation of one-carbon metabolites may play a key role in the development of stunting. Trial registration ClinicalTrials.gov NCT01375647. Registered 17 June 2011, retrospectively registered, https://clinicaltrials.gov/ct2/show/NCT01375647 .

SpikChIP: a novel computational methodology to compare multiple ChIP-seq using spike-in chromatin

In order to evaluate cell- and disease-specific changes in the interacting strength of chromatin targets, ChIP-seq signal across multiple conditions must undergo robust normalization. However, this is not possible using the standard ChIP-seq scheme, which lacks a reference for the control of biological and experimental variabilities. While several studies have recently proposed different solutions to circumvent this problem, substantial analytical differences among methodologies could hamper the experimental reproducibility and quantitative accuracy. Here, we propose a computational method to accurately compare ChIP-seq experiments, with exogenous spike-in chromatin, across samples in a genome-wide manner by using a local regression strategy (spikChIP). In contrast to the previous methodologies, spikChIP reduces the influence of sequencing noise of spike-in material during ChIP-seq normalization, while minimizes the overcorrection of non-occupied genomic regions in the experimental ChIP-seq. We demonstrate the utility of spikChIP with both histone and non-histone chromatin protein, allowing us to monitor for experimental reproducibility and the accurate ChIP-seq comparison of distinct experimental schemes. spikChIP software is available on GitHub (https://github.com/eblancoga/spikChIP).

CTCF looping is established during gastrulation in medaka embryos

Chromatin looping plays an important role in genome regulation. However, because ChIP-seq and loop-resolution Hi-C (DNA-DNA proximity ligation) are extremely challenging in mammalian early embryos, the developmental stage at which cohesin-mediated loops form remains unknown. Here, we study early development in medaka (the Japanese killifish, Oryzias latipes) at 12 time points before, during, and after gastrulation (the onset of cell differentiation) and characterize transcription, protein binding, and genome architecture. We find that gastrulation is associated with drastic changes in genome architecture, including the formation of the first loops between sites bound by the insulator protein CTCF and a large increase in the size of contact domains. In contrast, the binding of the CTCF is fixed throughout embryogenesis. Loops form long after genome-wide transcriptional activation, and long after domain formation seen in mouse embryos. These results suggest that, although loops may play a role in differentiation, they are not required for zygotic transcription. When we repeated our experiments in zebrafish, loops did not emerge until gastrulation, that is, well after zygotic genome activation. We observe that loop positions are highly conserved in synteny blocks of medaka and zebrafish, indicating that the 3D genome architecture has been maintained for >110–200 million years of evolution.

Genome-Wide Chromatin Analysis of FFPE Tissues Using a Dual-Arm Robot with Clinical Potential

Although chromatin immunoprecipitation and next-generation sequencing (ChIP-seq) using formalin-fixed paraffin-embedded tissue (FFPE) has been reported, it remained elusive whether they retained accurate transcription factor binding. Here, we developed a method to identify the binding sites of the insulator transcription factor CTCF and the genome-wide distribution of histone modifications involved in transcriptional activation. Importantly, we provide evidence that the ChIP-seq datasets obtained from FFPE samples are similar to or even better than the data for corresponding fresh-frozen samples, indicating that FFPE samples are compatible with ChIP-seq analysis. H3K27ac ChIP-seq analyses of 69 FFPE samples using a dual-arm robot revealed that driver mutations in EGFR were distinguishable from pan-negative cases and were relatively homogeneous as a group in lung adenocarcinomas. Thus, our results demonstrate that FFPE samples are an important source for epigenomic research, enabling the study of histone modifications, nuclear chromatin structure, and clinical data.

The Drosophila Y Chromosome Affects Heterochromatin Integrity Genome-Wide

The Drosophila Y chromosome is gene poor and mainly consists of silenced, repetitive DNA. Nonetheless, the Y influences expression of hundreds of genes genome-wide, possibly by sequestering key components of the heterochromatin machinery away from other positions in the genome. To test the influence of the Y chromosome on the genome-wide chromatin landscape, we assayed the genomic distribution of histone modifications associated with gene activation (H3K4me3) or heterochromatin (H3K9me2 and H3K9me3) in fruit flies with varying sex chromosome complements (X0, XY, and XYY males; XX and XXY females). Consistent with the general deficiency of active chromatin modifications on the Y, we find that Y gene dose has little influence on the genomic distribution of H3K4me3. In contrast, both the presence and the number of Y chromosomes strongly influence genome-wide enrichment patterns of repressive chromatin modifications. Highly repetitive regions such as the pericentromeres, the dot, and the Y chromosome (if present) are enriched for heterochromatic modifications in wildtype males and females, and even more strongly in X0 flies. In contrast, the additional Y chromosome in XYY males and XXY females diminishes the heterochromatic signal in these normally silenced, repeat-rich regions, which is accompanied by an increase in expression of Y-linked repeats. We find hundreds of genes that are expressed differentially between individuals with aberrant sex chromosome karyotypes, many of which also show sex-biased expression in wildtype Drosophila. Thus, Y chromosomes influence heterochromatin integrity genome-wide, and differences in the chromatin landscape of males and females may also contribute to sex-biased gene expression and sexual dimorphisms.

Toxic Y chromosome: Increased repeat expression and age-associated heterochromatin loss in male Drosophila with a young Y chromosome

Sex-specific differences in lifespan are prevalent across the tree of life and influenced by heteromorphic sex chromosomes. In species with XY sex chromosomes, females often outlive males. Males and females can differ in their overall repeat content due to the repetitive Y chromosome, and repeats on the Y might lower survival of the heterogametic sex (toxic Y effect). Here, we take advantage of the well-assembled young Y chromosome of Drosophila miranda to study the sex-specific dynamics of chromatin structure and repeat expression during aging in male and female flies. Male D. miranda have about twice as much repetitive DNA compared to females, and live shorter than females. Heterochromatin is crucial for silencing of repetitive elements, yet old D. miranda flies lose H3K9me3 modifications in their pericentromere, with heterochromatin loss being more severe during aging in males than females. Satellite DNA becomes de-repressed more rapidly in old vs. young male flies relative to females. In contrast to what is observed in D. melanogaster, we find that transposable elements (TEs) are expressed at higher levels in male D. miranda throughout their life. We show that epigenetic silencing via heterochromatin formation is ineffective on the TE-rich neo-Y chromosome, presumably due to active transcription of a large number of neo-Y linked genes, resulting in up-regulation of Y-linked TEs already in young males. This is consistent with an interaction between the evolutionary age of the Y chromosome and the genomic effects of aging. Our data support growing evidence that "toxic Y chromosomes" can diminish male fitness and a reduction in heterochromatin can contribute to sex-specific aging.

A new approach for quantifying epigenetic landscapes

ChIP-Seq is a widespread experimental method for determining the global enrichment of chromatin modifications and genome-associated factors. Whereas it is straightforward to compare the relative genomic distribution of these epigenetic features, researchers have also made efforts to compare their signal strength using external references for normalization. New work now suggests that these "spike-ins" could lead to inaccurate conclusions due to intrinsic issues of the methodology and instead calls for new criteria of experimental reporting that may permit internal standardization when certain parameters are fulfilled.

Dynamics of transcription-mediated conversion from euchromatin to facultative heterochromatin at the Xist promoter by Tsix

The fine-scale dynamics from euchromatin (EC) to facultative heterochromatin (fHC) has remained largely unclear. Here, we focus on Xist and its silencing initiator Tsix as a paradigm of transcription-mediated conversion from EC to fHC. In mouse epiblast stem cells, induction of Tsix recapitulates the conversion at the Xist promoter. Investigating the dynamics reveals that the conversion proceeds in a stepwise manner. Initially, a transient opened chromatin structure is observed. In the second step, gene silencing is initiated and dependent on Tsix, which is reversible and accompanied by simultaneous changes in multiple histone modifications. At the last step, maintenance of silencing becomes independent of Tsix and irreversible, which correlates with occupation of the -1 position of the transcription start site by a nucleosome and initiation of DNA methylation introduction. This study highlights the hierarchy of multiple chromatin events upon stepwise gene silencing establishment.

Methods for ChIP-seq analysis: A practical workflow and advanced applications

Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is a central method in epigenomic research. Genome-wide analysis of histone modifications, such as enhancer analysis and genome-wide chromatin state annotation, enables systematic analysis of how the epigenomic landscape contributes to cell identity, development, lineage specification, and disease. In this review, we first present a typical ChIP-seq analysis workflow, from quality assessment to chromatin-state annotation. We focus on practical, rather than theoretical, approaches for biological studies. Next, we outline various advanced ChIP-seq applications and introduce several state-of-the-art methods, including prediction of gene expression level and chromatin loops from epigenome data and data imputation. Finally, we discuss recently developed single-cell ChIP-seq analysis methodologies that elucidate the cellular diversity within complex tissues and cancers.

Induction of spontaneous human neocentromere formation and long-term maturation

Human centromeres form primarily on α-satellite DNA but sporadically arise de novo at naive ectopic loci, creating neocentromeres. Centromere inheritance is driven primarily by chromatin containing the histone H3 variant CENP-A. Here, we report a chromosome engineering system for neocentromere formation in human cells and characterize the first experimentally induced human neocentromere at a naive locus. The spontaneously formed neocentromere spans a gene-poor 100-kb domain enriched in histone H3 lysine 9 trimethylated (H3K9me3). Long-read sequencing revealed this neocentromere was formed by purely epigenetic means and assembly of a functional kinetochore correlated with CENP-A seeding, eviction of H3K9me3 and local accumulation of mitotic cohesin and RNA polymerase II. At formation, the young neocentromere showed markedly reduced chromosomal passenger complex (CPC) occupancy and poor sister chromatin cohesion. However, long-term tracking revealed increased CPC assembly and low-level transcription providing evidence for centromere maturation over time.

Analysis of Myc Chromatin Binding by Calibrated ChIP-Seq Approach

Here, we present a strategy to map and quantify the interactions between Myc and chromatin using a calibrated Myc ChIP-seq approach. We recommend the use of an internal spike-in control for post-sequencing normalization to enable detection of broad changes in Myc binding as can occur under conditions with varied Myc abundance. We also highlight a range of bioinformatic analyses that can dissect the downstream effects of Myc binding. These methods include peak calling, mapping Myc onto an integrated metagenome, juxtaposing ChIP-seq data with matching RNA-seq data, and identifying gene ontologies enriched for genes with high Myc binding. Our aim is to provide a guided strategy, from cell harvest through to bioinformatic analysis, to elucidate the global effects of Myc on transcription.

Histone Recycling by FACT and Spt6 during Transcription Prevents the Scrambling of Histone Modifications

Genomic DNA is framed by additional layers of information, referred to as the epigenome. Epigenomic marks such as DNA methylation, histone modifications, and histone variants are concentrated on specific genomic sites, where they can both instruct and reflect gene expression. How this information is maintained, notably in the face of transcription, is not completely understood. Specifically, the extent to which modified histones themselves are retained through RNA polymerase II passage is unclear. Here, we show that several histone modifications are mislocalized when the transcription-coupled histone chaperones FACT or Spt6 are disrupted in Saccharomyces cerevisiae. In the absence of functional FACT or Spt6, transcription generates nucleosome loss, which is partially compensated for by the increased activity of non-transcription-coupled histone chaperones. The random incorporation of transcription-evicted modified histones scrambles epigenomic information. Our work highlights the importance of local recycling of modified histones by FACT and Spt6 during transcription in the maintenance of the epigenomic landscape.

Targeted inhibition of KDM6 histone demethylases eradicates tumor-initiating cells via enhancer reprogramming in colorectal cancer

Tumor-initiating cells (TICs) maintain heterogeneity within tumors and seed metastases at distant sites, contributing to therapeutic resistance and disease recurrence. In colorectal cancer (CRC), strategy that effectively eradicates TICs and is of potential value for clinical use still remains in need. Methods: The anti-tumorigenic activity of a small-molecule inhibitor of KDM6 histone demethylases named GSK-J4 in CRC was evaluated by in vitro assays and in vivo imaging of xenografted tumors. Sphere formation, flow cytometry analysis of cell surface markers and intestinal organoid formation were performed to examine the impact of GSK-J4 on TIC properties. Transcriptome analysis and global profiling of H3K27ac, H3K27me3, and KDM6A levels by ChIP-seq were conducted to elucidate how KDM6 inhibition reshapes epigenetic landscape and thereby eliminating TICs. Results: GSK-J4 alleviated the malignant phenotypes of CRC cells in vitro and in vivo, sensitized them to chemotherapeutic treatment, and strongly repressed TIC properties and stemness-associated gene signatures in these cells. Mechanistically, KDM6 inhibition induced global enhancer reprogramming with a preferential impact on super-enhancer-associated genes, including some key genes that control stemness in CRC such as ID1. Besides, expression of both Kdm6a and Kdm6b was more abundant in mouse intestinal crypt when compared with upper villus and inhibition of their activities blocked intestinal organoid formation. Finally, we unveiled the power of KDM6B in predicting both the overall survival outcome and recurrence of CRC patients. Conclusions: Our study provides a novel rational strategy to eradicate TICs through reshaping epigenetic landscape in CRC, which might also be beneficial for optimizing current therapeutics.

Chromatin Hyperacetylation Impacts Chromosome Folding by Forming a Nuclear Subcompartment

Delineating how chromosomes fold at length scales beyond one megabase remains obscure relative to smaller-scale folding into TADs, loops, and nucleosomes. We find that rather than simply unfolding chromatin, histone hyperacetylation results in interactions between distant genomic loci separated by tens to hundreds of megabases, even in the absence of transcription. These hyperacetylated "megadomains" are formed by the BRD4-NUT fusion oncoprotein, interact both within and between chromosomes, and form a specific nuclear subcompartment that has elevated gene activity with respect to other subcompartments. Pharmacological degradation of BRD4-NUT results in collapse of megadomains and attenuation of the interactions between them. In contrast, these interactions persist and contacts between newly acetylated regions are formed after inhibiting RNA polymerase II initiation. Our structure-function approach thus reveals that broad chromatin domains of identical biochemical composition, independent of transcription, form nuclear subcompartments, and also indicates the potential of altering chromosome structure for treating human disease.

MAF1 is a chronic repressor of RNA polymerase III transcription in the mouse

Maf1^−/− mice are lean, obesity-resistant and metabolically inefficient. Their increased energy expenditure is thought to be driven by a futile RNA cycle that reprograms metabolism to meet an increased demand for nucleotides stemming from the deregulation of RNA polymerase (pol) III transcription. Metabolic changes consistent with this model have been reported in both fasted and refed mice, however the impact of the fasting-refeeding-cycle on pol III function has not been examined. Here we show that changes in pol III occupancy in the liver of fasted versus refed wild-type mice are largely confined to low and intermediate occupancy genes; high occupancy genes are unchanged. However, in Maf1^−/− mice, pol III occupancy of the vast majority of active loci in liver and the levels of specific precursor tRNAs in this tissue and other organs are higher than wild-type in both fasted and refed conditions. Thus, MAF1 functions as a chronic repressor of active pol III loci and can modulate transcription under different conditions. Our findings support the futile RNA cycle hypothesis, elaborate the mechanism of pol III repression by MAF1 and demonstrate a modest effect of MAF1 on global translation via reduced mRNA levels and translation efficiencies for several ribosomal proteins.

The Y chromosome may contribute to sex-specific ageing in Drosophila

Heterochromatin suppresses repetitive DNA, and a loss of heterochromatin has been observed in aged cells of several species, including humans and Drosophila. Males often contain substantially more heterochromatic DNA than females, due to the presence of a large, repeat-rich Y chromosome, and male flies generally have shorter average life spans than females. Here we show that repetitive DNA becomes de-repressed more rapidly in old male flies relative to females, and repeats on the Y chromosome are disproportionally mis-expressed during aging. This is associated with a loss of heterochromatin at repetitive elements during aging in male flies, and a general loss of repressive chromatin in aged males away from pericentromeric regions and the Y. By generating flies with different sex chromosome karyotypes (XXY females; X0 and XYY males), we show that repeat de-repression and average lifespan is correlated with the number of Y chromosomes. This suggests that sex-specific chromatin differences may contribute to sex-specific aging in flies.

PRC1 Catalytic Activity Is Central to Polycomb System Function

The Polycomb repressive system is an essential chromatin-based regulator of gene expression. Despite being extensively studied, how the Polycomb system selects its target genes is poorly understood, and whether its histone-modifying activities are required for transcriptional repression remains controversial. Here, we directly test the requirement for PRC1 catalytic activity in Polycomb system function. To achieve this, we develop a conditional mutation system in embryonic stem cells that completely removes PRC1 catalytic activity. Using this system, we demonstrate that catalysis by PRC1 drives Polycomb chromatin domain formation and long-range chromatin interactions. Furthermore, we show that variant PRC1 complexes with DNA-binding activities occupy target sites independently of PRC1 catalytic activity, providing a putative mechanism for Polycomb target site selection. Finally, we discover that Polycomb-mediated gene repression requires PRC1 catalytic activity. Together these discoveries provide compelling evidence that PRC1 catalysis is central to Polycomb system function and gene regulation.

User-Friendly and Interactive Analysis of ChIP-Seq Data Using EaSeq

ChIP-seq is a central method to gain understanding of the regulatory networks in the genome of stem cells and during differentiation. Exploration and analysis of such genome-wide data often leads to unexpected discoveries and new hypotheses. It therefore accelerates and improves the discovery phase, when scientists with biological understanding are enabled to analyze and visualize data. EaSeq ( http://easeq.net ) offers integrated exploration of genome-wide data in a visual, versatile, user-friendly, and interactive manner that connects abstract interpretations to the signal distribution at the underlying loci. Here we introduce the interface, data types, and acquisition, and guide the reader through two example workflows. These workflows will enable the reader to perform genome-wide analysis and visualization of transcription factor binding sites and histone marks. This includes making basic plots; finding, annotating, sorting, and filtering of peaks; using EaSeq as a genome browser; measuring ChIP-seq signal and calculating ratios; as well as data import and export.

KDM2 proteins constrain transcription from CpG island gene promoters independently of their histone demethylase activity

CpG islands (CGIs) are associated with the majority of mammalian gene promoters and function to recruit chromatin modifying enzymes. It has therefore been proposed that CGIs regulate gene expression through chromatin-based mechanisms, however in most cases this has not been directly tested. Here, we reveal that the histone H3 lysine 36 (H3K36) demethylase activity of the CGI-binding KDM2 proteins contributes only modestly to the H3K36me2-depleted state at CGI-associated gene promoters and is dispensable for normal gene expression. Instead, we discover that KDM2 proteins play a widespread and demethylase-independent role in constraining gene expression from CGI-associated gene promoters. We further show that KDM2 proteins shape RNA Polymerase II occupancy but not chromatin accessibility at CGI-associated promoters. Together this reveals a demethylase-independent role for KDM2 proteins in transcriptional repression and uncovers a new function for CGIs in constraining gene expression.

Native internally calibrated chromatin immunoprecipitation for quantitative studies of histone post-translational modifications

Chromatin immunoprecipitation coupled to next-generation sequencing (ChIP-seq) has served as the central method for the study of histone modifications for the past decade. In ChIP-seq analyses, antibodies selectively capture nucleosomes bearing a modification of interest and the associated DNA is then mapped to the genome to determine the distribution of the mark. This approach has several important drawbacks: (i) ChIP interpretation necessitates the assumption of perfect antibody specificity, despite growing evidence that this is often not the case. (ii) Common methods for evaluating antibody specificity in other formats have little or no bearing on specificity within a ChIP experiment. (iii) Uncalibrated ChIP is reported as relative enrichment, which is biologically meaningless outside the experimental reference frame defined by a discrete immunoprecipitation (IP), thus preventing facile comparison across experimental conditions or modifications. (iv) Differential library amplification and loading onto next-generation sequencers, as well as computational normalization, can further compromise quantitative relationships that may exist between samples. Consequently, the researcher is presented with a series of potential pitfalls and is blind to nearly all of them. Here we provide a detailed protocol for internally calibrated ChIP (ICeChIP), a method we recently developed to resolve these problems by spike-in of defined nucleosomal standards within a ChIP procedure. This protocol is optimized for specificity and quantitative power, allowing for measurement of antibody specificity and absolute measurement of histone modification density (HMD) at genomic loci on a biologically meaningful scale enabling unambiguous comparisons. We provide guidance on optimal conditions for next-generation sequencing (NGS) and instructions for data analysis. This protocol takes between 17 and 18 h, excluding time for sequencing or bioinformatic analysis. The ICeChIP procedure enables accurate measurement of histone post-translational modifications (PTMs) genome-wide in mammalian cells as well as Drosophila melanogaster and Caenorhabditis elegans, indicating suitability for use in eukaryotic cells more broadly.

Promoter-specific dynamics of TATA-binding protein association with the human genome

Transcription factor binding to target sites in vivo is a dynamic process that involves cycles of association and dissociation, with individual proteins differing in their binding dynamics. The dynamics at individual sites on a genomic scale have been investigated in yeast cells, but comparable experiments have not been done in multicellular eukaryotes. Here, we describe a tamoxifen-inducible, time-course ChIP-seq approach to measure transcription factor binding dynamics at target sites throughout the human genome. As observed in yeast cells, the TATA-binding protein (TBP) typically displays rapid turnover at RNA polymerase (Pol) II-transcribed promoters, slow turnover at Pol III promoters, and very slow turnover at the Pol I promoter. Turnover rates vary widely among Pol II promoters in a manner that does not correlate with the level of TBP occupancy. Human Pol II promoters with slow TBP dissociation preferentially contain a TATA consensus motif, support high transcriptional activity of downstream genes, and are linked with specific activators and chromatin remodelers. These properties of human promoters with slow TBP turnover differ from those of yeast promoters with slow turnover. These observations suggest that TBP binding dynamics differentially affect promoter function and gene expression, possibly at the level of transcriptional reinitiation/bursting.

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

MOTIVATION

The traditional reads per million normalization method is inappropriate for the evaluation of ChIP-seq data when treatments or mutations have global effects. Changes in global levels of histone modifications can be detected with exogenous reference spike-in controls. However, most ChIP-seq studies overlook the normalization that must be corrected with spike-in. A method that retrospectively renormalizes datasets without spike-in is lacking.

RESULTS

ChIPseqSpikeInFree is a novel ChIP-seq normalization method to effectively determine scaling factors for samples across various conditions and treatments, which does not rely on exogenous spike-in chromatin or peak detection to reveal global changes in histone modification occupancy. Application of ChIPseqSpikeInFree on five datasets demonstrates that this in silico approach reveals a similar magnitude of global changes as the spike-in method does.

AVAILABILITY AND IMPLEMENTATION

St. Jude Cloud (https://pecan.stjude.cloud/permalink/spikefree) and St. Jude Github ( https://github.com/stjude/ChIPseqSpikeInFree).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.