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

A Multiplexed System for Quantitative Comparisons of Chromatin Landscapes

Genome-wide profiling of histone modifications can provide systematic insight into the regulatory elements and programs engaged in a given cell type. However, conventional chromatin immunoprecipitation and sequencing (ChIP-seq) does not capture quantitative information on histone modification levels, requires large amounts of starting material, and involves tedious processing of each individual sample. Here, we address these limitations with a technology that leverages DNA barcoding to profile chromatin quantitatively and in multiplexed format. We concurrently map relative levels of multiple histone modifications across multiple samples, each comprising as few as a thousand cells. We demonstrate the technology by monitoring dynamic changes following inhibition of p300, EZH2, or KDM5, by linking altered epigenetic landscapes to chromatin regulator mutations, and by mapping active and repressive marks in purified human hematopoietic stem cells. Hence, this technology enables quantitative studies of chromatin state dynamics across rare cell types, genotypes, environmental conditions, and drug treatments.

Developmental Dynamics of X-Chromosome Dosage Compensation by the DCC and H4K20me1 in C. elegans

In Caenorhabditis elegans, the dosage compensation complex (DCC) specifically binds to and represses transcription from both X chromosomes in hermaphrodites. The DCC is composed of an X-specific condensin complex that interacts with several proteins. During embryogenesis, DCC starts localizing to the X chromosomes around the 40-cell stage, and is followed by X-enrichment of H4K20me1 between 100-cell to comma stage. Here, we analyzed dosage compensation of the X chromosome between sexes, and the roles of dpy-27 (condensin subunit), dpy-21 (non-condensin DCC member), set-1 (H4K20 monomethylase) and set-4 (H4K20 di-/tri-methylase) in X chromosome repression using mRNA-seq and ChIP-seq analyses across several developmental time points. We found that the DCC starts repressing the X chromosomes by the 40-cell stage, but X-linked transcript levels remain significantly higher in hermaphrodites compared to males through the comma stage of embryogenesis. Dpy-27 and dpy-21 are required for X chromosome repression throughout development, but particularly in early embryos dpy-27 and dpy-21 mutations produced distinct expression changes, suggesting a DCC independent role for dpy-21. We previously hypothesized that the DCC increases H4K20me1 by reducing set-4 activity on the X chromosomes. Accordingly, in the set-4 mutant, H4K20me1 increased more from the autosomes compared to the X, equalizing H4K20me1 level between X and autosomes. H4K20me1 increase on the autosomes led to a slight repression, resulting in a relative effect of X derepression. H4K20me1 depletion in the set-1 mutant showed greater X derepression compared to equalization of H4K20me1 levels between X and autosomes in the set-4 mutant, indicating that H4K20me1 level is important, but X to autosomal balance of H4K20me1 contributes slightly to X-repression. Thus H4K20me1 is not only a downstream effector of the DCC [corrected].In summary, X chromosome dosage compensation starts in early embryos as the DCC localizes to the X, and is strengthened in later embryogenesis by H4K20me1.

Biological chromodynamics: a general method for measuring protein occupancy across the genome by calibrating ChIP-seq

Sequencing DNA fragments associated with proteins following in vivo cross-linking with formaldehyde (known as ChIP-seq) has been used extensively to describe the distribution of proteins across genomes. It is not widely appreciated that this method merely estimates a protein's distribution and cannot reveal changes in occupancy between samples. To do this, we tagged with the same epitope orthologous proteins in Saccharomyces cerevisiae and Candida glabrata, whose sequences have diverged to a degree that most DNA fragments longer than 50 bp are unique to just one species. By mixing defined numbers of C. glabrata cells (the calibration genome) with S. cerevisiae samples (the experimental genomes) prior to chromatin fragmentation and immunoprecipitation, it is possible to derive a quantitative measure of occupancy (the occupancy ratio - OR) that enables a comparison of occupancies not only within but also between genomes. We demonstrate for the first time that this 'internal standard' calibration method satisfies the sine qua non for quantifying ChIP-seq profiles, namely linearity over a wide range. Crucially, by employing functional tagged proteins, our calibration process describes a method that distinguishes genuine association within ChIP-seq profiles from background noise. Our method is applicable to any protein, not merely highly conserved ones, and obviates the need for the time consuming, expensive, and technically demanding quantification of ChIP using qPCR, which can only be performed on individual loci. As we demonstrate for the first time in this paper, calibrated ChIP-seq represents a major step towards documenting the quantitative distributions of proteins along chromosomes in different cell states, which we term biological chromodynamics.

Loss of the RNA polymerase III repressor MAF1 confers obesity resistance

MAF1 is a global repressor of RNA polymerase III transcription that regulates the expression of highly abundant noncoding RNAs in response to nutrient availability and cellular stress. Thus, MAF1 function is thought to be important for metabolic economy. Here we show that a whole-body knockout of Maf1 in mice confers resistance to diet-induced obesity and nonalcoholic fatty liver disease by reducing food intake and increasing metabolic inefficiency. Energy expenditure in Maf1(-/-) mice is increased by several mechanisms. Precursor tRNA synthesis was increased in multiple tissues without significant effects on mature tRNA levels, implying increased turnover in a futile tRNA cycle. Elevated futile cycling of hepatic lipids was also observed. Metabolite profiling of the liver and skeletal muscle revealed elevated levels of many amino acids and spermidine, which links the induction of autophagy in Maf1(-/-) mice with their extended life span. The increase in spermidine was accompanied by reduced levels of nicotinamide N-methyltransferase, which promotes polyamine synthesis, enables nicotinamide salvage to regenerate NAD(+), and is associated with obesity resistance. Consistent with this, NAD(+) levels were increased in muscle. The importance of MAF1 for metabolic economy reveals the potential for MAF1 modulators to protect against obesity and its harmful consequences.

Defining bacterial regulons using ChIP-seq

Chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) is a powerful method that identifies protein-DNA binding sites in vivo. Recent studies have illustrated the value of ChIP-seq in studying transcription factor binding in various bacterial species under a variety of growth conditions. These results show that in addition to identifying binding sites, correlation of ChIP-seq data with expression data can reveal important information about bacterial regulons and regulatory networks. In this chapter, we provide an overview of the current state of knowledge about ChIP-seq methodology in bacteria, from sample preparation to raw data analysis. We also describe visualization and various bioinformatic analyses of processed ChIP-seq data.

Calibrating ChIP-Seq with Nucleosomal Internal Standards to Measure Histone Modification Density Genome Wide

Chromatin immunoprecipitation (ChIP) serves as a central experimental technique in epigenetics research, yet there are serious drawbacks: it is a relative measurement, which untethered to any external scale obscures fair comparison among experiments; it employs antibody reagents that have differing affinities and specificities for target epitopes that vary in abundance; and it is frequently not reproducible. To address these problems, we developed Internal Standard Calibrated ChIP (ICeChIP), wherein a native chromatin sample is spiked with nucleosomes reconstituted from recombinant and semisynthetic histones on barcoded DNA prior to immunoprecipitation. ICeChIP measures local histone modification densities on a biologically meaningful scale, enabling unbiased trans-experimental comparisons, and reveals unique insight into the nature of bivalent domains. This technology provides in situ assessment of the immunoprecipitation step, accommodating for many experimental pitfalls as well as providing a critical examination of untested assumptions inherent to conventional ChIP.

Genome-wide assays that identify and quantify modified cytosines in human disease studies

The number of different assays that has been published to study DNA methylation is extensive, complemented by recently described assays that test modifications of cytosine other than the most abundant 5-methylcytosine (5mC) variant. In this review, we describe the considerations involved in choosing how to study 5mC throughout the genome, with an emphasis on the common application of testing for epigenetic dysregulation in human disease. While microarray studies of 5mC continue to be commonly used, these lack the additional qualitative information from sequencing-based approaches that is increasingly recognized to be valuable. When we test the representation of functional elements in the human genome by several current assay types, we find that no survey approach interrogates anything more than a small minority of the nonpromoter cis-regulatory sites where DNA methylation variability is now appreciated to influence gene expression and to be associated with human disease. However, whole-genome bisulphite sequencing (WGBS) adds a substantial representation of loci at which DNA methylation changes are unlikely to be occurring with transcriptional consequences. Our assessment is that the most effective approach to DNA methylation studies in human diseases is to use targeted bisulphite sequencing of the cis-regulatory loci in a cell type of interest, using a capture-based or comparable system, and that no single design of a survey approach will be suitable for all cell types.

Aggregate and Heatmap Representations of Genome-Wide Localization Data Using VAP, a Versatile Aggregate Profiler

In the analysis of experimental data corresponding to the signal enrichment of chromatin features such as histone modifications throughout the genome, it is often useful to represent the signal over known regions of interest, such as genes, using aggregate or individual profiles. In the present chapter, we describe and explain the best practices on how to generate such profiles as well as other usages of the versatile aggregate profiler (VAP) tool (Coulombe et al., Nucleic Acids Res 42:W485-W493, 2014), with a particular focus on the new functionalities introduced in version 1.1.0 of VAP.

Quantitative ChIP-Seq normalization reveals global modulation of the epigenome

Epigenomic profiling by chromatin immunoprecipitation coupled with massively parallel DNA sequencing (ChIP-seq) is a prevailing methodology used to investigate chromatin-based regulation in biological systems such as human disease, but the lack of an empirical methodology to enable normalization among experiments has limited the precision and usefulness of this technique. Here, we describe a method called ChIP with reference exogenous genome (ChIP-Rx) that allows one to perform genome-wide quantitative comparisons of histone modification status across cell populations using defined quantities of a reference epigenome. ChIP-Rx enables the discovery and quantification of dynamic epigenomic profiles across mammalian cells that would otherwise remain hidden using traditional normalization methods. We demonstrate the utility of this method for measuring epigenomic changes following chemical perturbations and show how reference normalization of ChIP-seq experiments enables the discovery of disease-relevant changes in histone modification occupancy.

Identifying and mitigating bias in next-generation sequencing methods for chromatin biology

Next-generation sequencing (NGS) technologies have been used in diverse ways to investigate various aspects of chromatin biology by identifying genomic loci that are bound by transcription factors, occupied by nucleosomes or accessible to nuclease cleavage, or loci that physically interact with remote genomic loci. However, reaching sound biological conclusions from such NGS enrichment profiles requires many potential biases to be taken into account. In this Review, we discuss common ways in which biases may be introduced into NGS chromatin profiling data, approaches to diagnose these biases and analytical techniques to mitigate their effect.