Genome-wide evaluation of histone methylation changes associated with leaf senescence in Arabidopsis

Leaf senescence is the orderly dismantling of older tissue that allows recycling of nutrients to developing portions of the plant and is accompanied by major changes in gene expression. Histone modifications correlate to levels of gene expression, and this study utilizes ChIP-seq to classify activating H3K4me3 and silencing H3K27me3 marks on a genome-wide scale for soil-grown mature and naturally senescent Arabidopsis leaves. ChIPnorm was used to normalize data sets and identify genomic regions with significant differences in the two histone methylation patterns, and the differences were correlated to changes in gene expression. Genes that showed an increase in the H3K4me3 mark in older leaves were senescence up-regulated, while genes that showed a decrease in the H3K4me3 mark in the older leaves were senescence down-regulated. For the H3K27me3 modification, genes that lost the H3K27me3 mark in older tissue were senescence up-regulated. Only a small number of genes gained the H3K27me3 mark, and these were senescence down-regulated. Approximately 50% of senescence up-regulated genes lacked the H3K4me3 mark in both mature and senescent leaf tissue. Two of these genes, SAG12 and At1g73220, display strong senescence up-regulation without the activating H3K4me3 histone modification. This study provides an initial epigenetic framework for the developmental transition into senescence.     

Investigating repetitively matching short sequencing reads: The enigmatic nature of H3K9me3

Most histone modifications can easily be characterized as either activating or repressive. For example, histone3, lysine 4 trimethylation (H3K4me3) is generally considered a distinct sign of actively transcribed promoters while H3K27me3 is generally found at repressed genes. This is not the case for H3K9me3, the subject of this communication, which is a modification that has traditionally been considered a mark of constitutive heterochromatin, but has also been found in significant levels in expressed genes. We therefore sought to use new high-throughput genome-wide maps of H3K9me3 localization to investigate the conflicting hypotheses concerning the nature of this modification. Before we could accurately analyze the locations of H3K9me3 along the genome, and especially in repetitive locations, we developed a method for accurately utilizing short sequencing reads that do not map uniquely to a location in the genome. Investigating the locations of H3K9me3 along the genome allowed us to determine that, while there are high levels of H3K9me3 outside of genes, this modification is not absent from genes. Therefore, we suggest that H3K9me3 may have a role in chromatin organization rather than being directly related to gene expression. In addition, we have found that there is a need to include repetitively matching reads in any high-throughput sequencing experiment.     

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.     

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.     

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.