reChIP-seq reveals widespread bivalency of H3K4me3 and H3K27me3 in CD4+ memory T cells

The combinatorial action of co-localizing chromatin modifications and regulators determines chromatin structure and function. However, identifying co-localizing chromatin features in a high-throughput manner remains a technical challenge. Here we describe a novel reChIP-seq approach and tailored bioinformatic analysis tool, normR that allows for the sequential enrichment and detection of co-localizing DNA-associated proteins in an unbiased and genome-wide manner. We illustrate the utility of the reChIP-seq method and normR by identifying H3K4me3 or H3K27me3 bivalently modified nucleosomes in primary human CD4+ memory T cells. We unravel widespread bivalency at hypomethylated CpG-islands coinciding with inactive promoters of developmental regulators. reChIP-seq additionally uncovered heterogeneous bivalency in the population, which was undetectable by intersecting H3K4me3 and H3K27me3 ChIP-seq tracks. Finally, we provide evidence that bivalency is established and stabilized by an interplay between the genome and epigenome. Our reChIP-seq approach augments conventional ChIP-seq and is broadly applicable to unravel combinatorial modes of action.

[1]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[2]  Ivan G. Costa,et al.  Detecting differential peaks in ChIP-seq signals with ODIN , 2015, Bioinform..

[3]  Dustin E. Schones,et al.  High-Resolution Profiling of Histone Methylations in the Human Genome , 2007, Cell.

[4]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[5]  I. Ial,et al.  Nature Communications , 2010, Nature Cell Biology.

[6]  C. Allis,et al.  The language of covalent histone modifications , 2000, Nature.

[7]  T. Mikkelsen,et al.  Genome-wide maps of chromatin state in pluripotent and lineage-committed cells , 2007, Nature.

[8]  Benjamin A. Garcia,et al.  Asymmetrically Modified Nucleosomes , 2012, Cell.

[9]  Guo-Liang Xu,et al.  The N-terminus of histone H3 is required for de novo DNA methylation in chromatin , 2009, Proceedings of the National Academy of Sciences.

[10]  P. Soloway,et al.  Single molecule and single cell epigenomics. , 2015, Methods.

[11]  K. Struhl,et al.  Analysis of protein co-occupancy by quantitative sequential chromatin immunoprecipitation. , 2004, Current protocols in molecular biology.

[12]  Bronwen L. Aken,et al.  GENCODE: The reference human genome annotation for The ENCODE Project , 2012, Genome research.

[13]  Yuka Kanno,et al.  Global mapping of H3K4me3 and H3K27me3 reveals specificity and plasticity in lineage fate determination of differentiating CD4+ T cells. , 2009, Immunity.

[14]  Cole Trapnell,et al.  TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions , 2013, Genome Biology.

[15]  Suresh Cuddapah,et al.  The genomic landscape of histone modifications in human T cells , 2006, Proceedings of the National Academy of Sciences.

[16]  B. Fierz,et al.  Histone H3K27 Trimethylation Inhibits H3 Binding and Function of SET1-Like H3K4 Methyltransferase Complexes , 2013, Molecular and Cellular Biology.

[17]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[18]  Isaac Dialsingh,et al.  Estimating the proportion of true null hypotheses when the statistics are discrete , 2015, Bioinform..

[19]  Tao Liu,et al.  Using MACS to Identify Peaks from ChIP‐Seq Data , 2011, Current protocols in bioinformatics.

[20]  Kevin Struhl,et al.  Quantitative sequential chromatin immunoprecipitation, a method for analyzing co-occupancy of proteins at genomic regions in vivo. , 2004, Nucleic acids research.

[21]  P. Laird,et al.  Bis-SNP: Combined DNA methylation and SNP calling for Bisulfite-seq data , 2012, Genome Biology.

[22]  James A. Cuff,et al.  A Bivalent Chromatin Structure Marks Key Developmental Genes in Embryonic Stem Cells , 2006, Cell.

[23]  Michael B. Stadler,et al.  Lineage-specific polycomb targets and de novo DNA methylation define restriction and potential of neuronal progenitors. , 2008, Molecular cell.

[24]  John D. Storey.  A direct approach to false discovery rates , 2002 .

[25]  Feng Lin,et al.  An HMM approach to genome-wide identification of differential histone modification sites from ChIP-seq data , 2008, Bioinform..

[26]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[27]  R. Young,et al.  A Chromatin Landmark and Transcription Initiation at Most Promoters in Human Cells , 2007, Cell.

[28]  M. Frommer,et al.  CpG islands in vertebrate genomes. , 1987, Journal of molecular biology.

[29]  C. Allis,et al.  Translating the Histone Code , 2001, Science.

[30]  T. Richmond,et al.  Crystal structure of the nucleosome core particle at 2.8 Å resolution , 1997, Nature.

[31]  R. Kornberg.  Chromatin structure: a repeating unit of histones and DNA. , 1974, Science.

[32]  M. Furlan-Magaril,et al.  Sequential chromatin immunoprecipitation protocol: ChIP-reChIP. , 2009, Methods in molecular biology.

[33]  A. Bauer,et al.  Histone methylation by PRC2 is inhibited by active chromatin marks. , 2011, Molecular cell.

[34]  Christodoulos A Floudas,et al.  High Throughput Characterization of Combinatorial Histone Codes* , 2009, Molecular & Cellular Proteomics.

[35]  Scott B. Dewell,et al.  Engineering of a Histone-Recognition Domain in Dnmt3a Alters the Epigenetic Landscape and Phenotypic Features of Mouse ESCs. , 2015, Molecular cell.

[36]  C. Allis,et al.  DNMT3L connects unmethylated lysine 4 of histone H3 to de novo methylation of DNA , 2007, Nature.

[37]  G. Mizuguchi,et al.  Stepwise Histone Replacement by SWR1 Requires Dual Activation with Histone H2A.Z and Canonical Nucleosome , 2010, Cell.

[38]  H. Ng,et al.  Uniform, optimal signal processing of mapped deep-sequencing data , 2013, Nature Biotechnology.

[39]  A. Bird,et al.  Synthetic CpG islands reveal DNA sequence determinants of chromatin structure , 2014, eLife.

[40]  Kairong Cui,et al.  H3.3/H2A.Z double variant-containing nucleosomes mark ‘nucleosome-free regions’ of active promoters and other regulatory regions in the human genome , 2009, Nature Genetics.

[41]  Bradley E. Bernstein,et al.  GC-Rich Sequence Elements Recruit PRC2 in Mammalian ES Cells , 2010, PLoS genetics.

[42]  A. Truax,et al.  ChIP and Re-ChIP assays: investigating interactions between regulatory proteins, histone modifications, and the DNA sequences to which they bind. , 2012, Methods in molecular biology.

[43]  Sarah L Vowler,et al.  Cooperative interaction between retinoic acid receptor-alpha and estrogen receptor in breast cancer. , 2010, Genes & development.

[44]  Clifford A. Meyer,et al.  Model-based Analysis of ChIP-Seq (MACS) , 2008, Genome Biology.

[45]  Xianhua Dai,et al.  BGDB: a database of bivalent genes , 2013, Database J. Biol. Databases Curation.