An HMM approach to genome-wide identification of differential histone modification sites from ChIP-seq data

MOTIVATION: Epigenetic modifications are one of the critical factors to regulate gene expression and genome function. Among different epigenetic modifications, the differential histone modification sites (DHMSs) are of great interest to study the dynamic nature of epigenetic and gene expression regulations among various cell types, stages or environmental responses. To capture the histone modifications at whole genome scale, ChIP-seq technology is becoming a robust and comprehensive approach. Thus the DHMSs are potentially identifiable by comparing two ChIP-seq libraries. However, little has been addressed on this issue in literature. RESULTS: Aiming at identifying DHMSs, we propose an approach called ChIPDiff for the genome-wide comparison of histone modification sites identified by ChIP-seq. Based on the observations of ChIP fragment counts, the proposed approach employs a hidden Markov model (HMM) to infer the states of histone modification changes at each genomic location. We evaluated the performance of ChIPDiff by comparing the H3K27me3 modification sites between mouse embryonic stem cell (ESC) and neural progenitor cell (NPC). We demonstrated that the H3K27me3 DHMSs identified by our approach are of high sensitivity, specificity and technical reproducibility. ChIPDiff was further applied to uncover the differential H3K4me3 and H3K36me3 sites between different cell states. Interesting biological discoveries were achieved from such comparison in our study.     

Stress‐associated H3K4 methylation accumulates during postnatal development and aging of rhesus macaque brain

Epigenetic modifications are critical determinants of cellular and developmental states. Epigenetic changes, such as decreased H3K27me3 histone methylation on insulin/IGF1 genes, have been previously shown to modulate lifespan through gene expression regulation. However, global epigenetic changes during aging and their biological functions, if any, remain elusive. Here, we examined the histone modification H3K4 dimethylation (H3K4me2) in the prefrontal cortex of individual rhesus macaques at different ages by chromatin immunoprecipitation, followed by deep sequencing (ChIP‐seq) at the whole genome level. Through integrative analysis of the ChIP‐seq profiles with gene expression data, we found that H3K4me2 increased at promoters and enhancers globally during postnatal development and aging, and those that correspond to gene expression changes in cis are enriched for stress responses, such as the DNA damage response. This suggests that metabolic and environmental stresses experienced by an organism are associated with the progressive opening of chromatin. In support of this, we also observed increased expression of two H3K4 methyltransferases, SETD7 and DPY30, in aged macaque brain.     

Stella Regulates the Development of Female Germline Stem Cells by Modulating Chromatin Structure and DNA Methylation

Female germline stem cells (FGSCs) have the ability to self-renew and differentiate into oocytes. Stella, encoded by a maternal effect gene, plays an important role in oogenesis and early embryonic development. However, its function in FGSCs remains unclear. In this study, we showed that CRISPR/Cas9-mediated knockout of Stella promoted FGSC proliferation and reduced the level of genome-wide DNA methylation of FGSCs. Conversely, Stella overexpression led to the opposite results, and enhanced FGSC differentiation. We also performed an integrative analysis of chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq), high-throughput genome-wide chromosome conformation capture (Hi-C), and use of our published epigenetic data. Results indicated that the binding sites of STELLA and active histones H3K4me3 and H3K27ac were enriched near the TAD boundaries. Hi-C analysis showed that Stella overexpression attenuated the interaction within TADs, and interestingly enhanced the TAD boundary strength in STELLA-associated regions. Taking these findings together, our study not only reveals the role of Stella in regulating DNA methylation and chromatin structure, but also provides a better understanding of FGSC development.     

Heterochromatin Controls γH2A Localization in Neurospora crassa

In response to genotoxic stress, ATR and ATM kinases phosphorylate H2A in fungi and H2AX in animals on a C-terminal serine. The resulting modified histone, called γH2A, recruits chromatin-binding proteins that stabilize stalled replication forks or promote DNA double-strand-break repair. To identify genomic loci that might be prone to replication fork stalling or DNA breakage in Neurospora crassa, we performed chromatin immunoprecipitation (ChIP) of γH2A followed by next-generation sequencing (ChIP-seq). γH2A-containing nucleosomes are enriched in Neurospora heterochromatin domains. These domains are comprised of A·T-rich repetitive DNA sequences associated with histone H3 methylated at lysine-9 (H3K9me), the H3K9me-binding protein heterochromatin protein 1 (HP1), and DNA cytosine methylation. H3K9 methylation, catalyzed by DIM-5, is required for normal γH2A localization. In contrast, γH2A is not required for H3K9 methylation or DNA methylation. Normal γH2A localization also depends on HP1 and a histone deacetylase, HDA-1, but is independent of the DNA methyltransferase DIM-2. γH2A is globally induced in dim-5 mutants under normal growth conditions, suggesting that the DNA damage response is activated in these mutants in the absence of exogenous DNA damage. Together, these data suggest that heterochromatin formation is essential for normal DNA replication or repair.     

An Assessment of Fixed and Native Chromatin Preparation Methods to Study Histone Post-Translational Modifications at a Whole Genome Scale in Skeletal Muscle Tissue

Background

Genomic loci associated with histone marks are typically analyzed by immunoprecipitation of the chromatin followed by quantitative-PCR (ChIP-qPCR) or high throughput sequencing (ChIP-seq). Chromatin can be either cross-linked (X-ChIP) or used in the native state (N-ChIP). Cross-linking of DNA and proteins helps stabilizing their interactions before analysis. Despite X-ChIP is the most commonly used method, muscle tissue fixation is known to be relatively inefficient. Moreover, no protocol described a simple and reliable preparation of skeletal muscle chromatin of sufficient quality for subsequent high-throughput sequencing. Here we aimed to set-up and compare both chromatin preparation methods for a genome-wide analysis of H3K27me3, a broad-peak histone mark, using chicken P. major muscle tissue.

Results

Fixed and unfixed chromatin were prepared from chicken muscle tissues (Pectoralis major). Chromatin fixation, shearing by sonication or digestion and immunoprecipitation performed equivalently. High-quality Illumina reads were obtained (q30 > 93%). The bioinformatic analysis of the data was performed using epic, a tool based on SICER, and MACS2. Forty millions of reads were analyzed for both X-ChIP-seq and N-ChIP-seq experiments. Surprisingly, H3K27me3 X-ChIP-seq analysis led to the identification of only 2000 enriched regions compared to about 15,000 regions identified in the case of N-ChIP-seq. N-ChIP-seq peaks were more consistent between replicates compared to X-ChIP-seq. Higher N-ChIP-seq enrichments were confirmed by ChIP-qPCR at the PAX5 and SOX2 loci known to be enriched for H3K27me3 in myotubes and at the loci of common regions of enrichment identified in this study.

Conclusions

Our findings suggest that the preparation of muscle chromatin for ChIP-seq in cross-linked conditions can compromise the systematic analysis of broad histone marks. Therefore, native chromatin preparation should be preferred to cross-linking when a ChIP experiment has to be performed on skeletal muscle tissue, particularly when a broad source signal is considered.

     

The Arabidopsis SWI2/SNF2 Chromatin Remodeler BRAHMA Regulates Polycomb Function during Vegetative Development and Directly Activates the Flowering Repressor Gene SVP

The chromatin remodeler BRAHMA (BRM) is a Trithorax Group (TrxG) protein that antagonizes the functions of Polycomb Group (PcG) proteins in fly and mammals. Recent studies also implicate such a role for Arabidopsis (Arabidopsis thaliana) BRM but the molecular mechanisms underlying the antagonism are unclear. To understand the interplay between BRM and PcG during plant development, we performed a genome-wide analysis of trimethylated histone H3 lysine 27 (H3K27me3) in brm mutant seedlings by chromatin immunoprecipitation followed by next generation sequencing (ChIP-seq). Increased H3K27me3 deposition at several hundred genes was observed in brm mutants and this increase was partially supressed by removal of the H3K27 methyltransferase CURLY LEAF (CLF) or SWINGER (SWN). ChIP experiments demonstrated that BRM directly binds to a subset of the genes and prevents the inappropriate association and/or activity of PcG proteins at these loci. Together, these results indicate a crucial role of BRM in restricting the inappropriate activity of PcG during plant development. The key flowering repressor gene SHORT VEGETATIVE PHASE (SVP) is such a BRM target. In brm mutants, elevated PcG occupancy at SVP accompanies a dramatic increase in H3K27me3 levels at this locus and a concomitant reduction of SVP expression. Further, our gain- and loss-of-function genetic evidence establishes that BRM controls flowering time by directly activating SVP expression. This work reveals a genome-wide functional interplay between BRM and PcG and provides new insights into the impacts of these proteins in plant growth and development.     

Quantitative assessment of chromatin immunoprecipitation grade antibodies directed against histone modifications reveals patterns of co-occurring marks on histone protein molecules

The defining step in most chromatin immunoprecipitation (ChIP) assays is the use of an antibody to enrich for a particular protein or histone modification state associated with segments of chromatin. The specificity of the antibody is critical to the interpretation of the experiment, yet this property is rarely reported. Here, we present a quantitative method using mass spectrometry to characterize the specificity of key histone H3 modification-targeting antibodies that have previously been used to characterize the "histone code." We further extend the use of these antibody reagents to the observation of long range correlations among disparate histone modifications. Using purified human histones representing the mixture of chromatin states present in living cells, we were able to quantify the degree of target enrichment and the specificity of several commonly used, commercially available ChIP grade antibodies. We found significant differences in enrichment efficiency among various reagents directed against four frequently studied chromatin marks: H3K4me2, H3K4me3, H3K9me3, and H3K27me3. For some antibodies, we also detected significant off target enrichment of alternate modifications at the same site (i.e., enrichment of H3K4me2 by an antibody directed against H3K4me3). Through cluster analysis, we were able to recognize patterns of co-enrichment of marks at different sites on the same histone protein. Surprisingly, these co-enrichments corresponded well to "canonical" chromatin states that are exemplary of activated and repressed regions of chromatin. Altogether, our findings suggest that 1) the results of ChIP experiments need to be evaluated with caution given the potential for cross-reactivity of the commonly used histone modification recognizing antibodies, 2) multiple marks with consistent biological interpretation exist on the same histone protein molecule, and 3) some components of the histone code may be transduced on single proteins in living cells.     

ChIPnorm: A Statistical Method for Normalizing and Identifying Differential Regions in Histone Modification ChIP-seq Libraries

The advent of high-throughput technologies such as ChIP-seq has made possible the study of histone modifications. A problem of particular interest is the identification of regions of the genome where different cell types from the same organism exhibit different patterns of histone enrichment. This problem turns out to be surprisingly difficult, even in simple pairwise comparisons, because of the significant level of noise in ChIP-seq data. In this paper we propose a two-stage statistical method, called ChIPnorm, to normalize ChIP-seq data, and to find differential regions in the genome, given two libraries of histone modifications of different cell types. We show that the ChIPnorm method removes most of the noise and bias in the data and outperforms other normalization methods. We correlate the histone marks with gene expression data and confirm that histone modifications H3K27me3 and H3K4me3 act as respectively a repressor and an activator of genes. Compared to what was previously reported in the literature, we find that a substantially higher fraction of bivalent marks in ES cells for H3K27me3 and H3K4me3 move into a K27-only state. We find that most of the promoter regions in protein-coding genes have differential histone-modification sites. The software for this work can be downloaded from http://lcbb.epfl.ch/software.html.