Quantitative Analysis of Chip-Seq Data Uncovers Dynamic and Sustained H3k4me3 and H3k27me3 Modulation in Cancer Cells Under Hypo

Total Page:16

File Type:pdf, Size:1020Kb

Quantitative Analysis of Chip-Seq Data Uncovers Dynamic and Sustained H3k4me3 and H3k27me3 Modulation in Cancer Cells Under Hypo Adriaens et al. Epigenetics & Chromatin (2016) 9:48 DOI 10.1186/s13072-016-0090-4 Epigenetics & Chromatin RESEARCH Open Access Quantitative analysis of ChIP‑seq data uncovers dynamic and sustained H3K4me3 and H3K27me3 modulation in cancer cells under hypoxia Michiel E. Adriaens1,2*† , Peggy Prickaerts3†, Michelle Chan‑Seng‑Yue4,8, Twan van den Beucken5,6,7, Vivian E. H. Dahlmans3, Lars M. Eijssen2, Timothy Beck4,9, Bradly G. Wouters5,6,7, Jan Willem Voncken3 and Chris T. A. Evelo2 Abstract Background: A comprehensive assessment of the epigenetic dynamics in cancer cells is the key to understand‑ ing the molecular mechanisms underlying cancer and to improving cancer diagnostics, prognostics and treat‑ ment. By combining genome-wide ChIP-seq epigenomics and microarray transcriptomics, we studied the effects of oxygen deprivation and subsequent reoxygenation on histone 3 trimethylation of lysine 4 (H3K4me3) and lysine 27 (H3K27me3) in a breast cancer cell line, serving as a model for abnormal oxygenation in solid tumors. A priori, epigenetic markings and gene expression levels not only are expected to vary greatly between hypoxic and normoxic conditions, but also display a large degree of heterogeneity across the cell population. Where traditionally ChIP-seq data are often treated as dichotomous data, the model and experiment here necessitate a quantitative, data-driven analysis of both datasets. Results: We first identified genomic regions with sustained epigenetic markings, which provided a sample-specific reference enabling quantitative ChIP-seq data analysis. Sustained H3K27me3 marking was located around cen‑ tromeres and intergenic regions, while sustained H3K4me3 marking is associated with genes involved in RNA binding, translation and protein transport and localization. Dynamic marking with both H3K4me3 and H3K27me3 (hypoxia- induced bivalency) was found in CpG-rich regions at loci encoding factors that control developmental processes, congruent with observations in embryonic stem cells. Conclusions: In silico-identified epigenetically sustained and dynamic genomic regions were confirmed through ChIP-PCR in vitro, and obtained results are corroborated by published data and current insights regarding epigenetic regulation. Keywords: Epigenomics, Transcriptomics, Data normalization, ChIP-sequencing, H3K4me3, H3K27me3, Hypoxia, MCF7 *Correspondence: [email protected] †Michiel E Adriaens and Peggy Prickaerts have contributed equally to this work 1 Maastricht Centre for Systems Biology – MaCSBio, Maastricht University, Maastricht, The Netherlands Full list of author information is available at the end of the article © 2016 The Author(s). This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/ publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated. Adriaens et al. Epigenetics & Chromatin (2016) 9:48 Page 2 of 11 Background unsupervised approach is required to identify and define Epigenomics explores genome-wide epigenetic modi- epigenetically invariant genomic regions. fications, such as DNA (hydroxy)methylation or post- Second, although biological interpretation of ChIP-seq translational modification of N-terminal histone tails data by itself is legitimate [19–21], integration of tran- [1]; such marking is often dynamic and serves to support scriptomics data is pivotal for a comprehensive biological phenotypic plasticity. Well-known examples of histone interpretation [22]. We extrapolate from this observation modifications are trimethylation of histone 3 lysine 4 to propose that transcriptomics data analysis is a require- (H3K4me3) and lysine 27 (H3K27me3), which are often ment for biology-driven preprocessing of ChIP-seq data. associated with transcriptional activity and repression, Specifically, the invariant regions used for ChIP-seq data respectively. Combined chromatin immunoprecipita- normalization should be defined in conjunction with tion (ChIP) and high-throughput sequencing technology transcriptomics data. (ChIP-seq [2, 3]) has proven to be a powerful approach to Considering the challenges presented by a highly vari- create genome-wide maps of such histone modifications ant biological system, we pursued the invariant-set para- in multiple human cancer cell lines [4–7], and the key digm outlined above in a data-driven analysis approach to understanding the molecular mechanisms underlying combining genes with sustained expression and genomic cancer and to improving cancer diagnostics, prognostics regions with sustained H3K4me3 and H3K27me3 mark- and treatment [1, 8, 9]. ings across all states. The approach enables quantitative In this study we focused specifically on ChIP-seq- data comparison for H3K4me3 and H3K27me3 ChIP-seq based analysis of dynamically changing H3K4me3 and data and a comprehensive integrative analysis of epig- H3K27me3 enrichment in response to fluctuating oxygen enomics and transcriptomics data. The biological insights levels. MCF7 cells were exposed to 8 and 24 h of hypoxia, derived from the analysis are further corroborated a and subsequently reoxygenated, serving as a model for posteriori by publicly available data and current insights such dynamic epigenetic effects occurring in response to regarding epigenetic regulation. reduced oxygen availability in solid tumors [10]. Analysis of a highly dynamic biological system, where the major- Results ity of epigenetic markings and gene expression levels will Data generation be altered between physiologically distinct states, as well MCF7 cells were exposed to 8 or 24 h of hypoxia (0.02 % as display a large degree of heterogeneity across the cell oxygen) and subsequently reoxygenated for 8 h (21 % population due to epigenetic selection and stochasticity, oxygen) (see Methods section for details). Cells were har- poses a number of challenges on data analysis. vested and processed for ChIP-seq and microarray-based First, the expected heterogeneity, i.e., the observation expression analysis at 0, 8 and 24 h of hypoxia (referred that the same region can harbor a specific epigenetic to as t = 0, t = 8 and t = 24, respectively), and at 8 h of mark in some part of the cell population, while not in reoxygenation (referred to henceforth as t =+8). the other, results in (1) a quantitative range of ChIP-seq signals within a state and (2) quantitative differences Quantitative analysis of sustained and dynamic H3K27me3 between states. Although ChIP-seq data are tradition- enrichment under fluctuating oxygen levels ally treated as dichotomous (i.e., present/not present), To enable quantitative analysis, H3K27me3 datasets were here quantitative data comparison is therefore a must. normalized based on identification of regions with sus- This requires normalization of the ChIP-seq data [11]. tained H3K27me3 enrichment between all samples ana- Most commonly, a scaling approach relative to the total lyzed. The cumulative area under the curve (AUC) for all number of aligned reads is applied, based on the predic- peaks in all these regions was determined for each con- tion that the number of differences between conditions is dition, resulting in sample-specific scaling factors (see limited [12–14]. This assumption, however, does not hold Methods section for details). here, as quantitative data comparison between highly For H3K27me3, after normalization, a single cutoff variant states requires invariant measurement sets across value for all samples was set at the peak height above all experimental conditions as a reference, such as an which the H3K27me3 data correlated with the ChIP- artificial spike-in [15–17], or biologically steadily marked seq data of CBX8, a known H3K27me3-binding protein H3K4me3/H3K27me3 regions (by analogy: steadily under normoxic conditions (t = 0) [23], corresponding expressed genes in transcriptomic studies). However, to a normalized height of 5.0 on a log2 scale. For biologi- a priori definition of regions with sustained epigenetic cal interpretation analyses, we set the level of biologi- marking is challenging, as this depends on the biologi- cal significance at two times this background level (6.0 cal model as well as the studied histone modification on a log2 scale) and considered peaks below this value [13]. To accommodate study-specific approaches [18], an background noise (i.e., biologically irrelevant). To add Adriaens et al. Epigenetics & Chromatin (2016) 9:48 Page 3 of 11 further biological validity to this threshold, we defined enrichment of H3K4me3 surrounding in TSS (region a set of 1137 genes whose expression levels were above defined as− 1000 bp to +100 bp of the TSS [4]) was the 95th percentile for all samples analyzed, i.e., consist- summed, yielding sample-specific scaling factors to ently highly expressed genes. Any repressive H3K27me3 which each peak in the dataset was scaled. A number of enrichment present in these genes were considered back- genes with dynamically altered or sustained expression ground noise, and indeed, the median enrichment of were validated by conventional ChIP-PCR. Oxygen dep- highly expressed genes H3K27me3 peaks was found
Recommended publications
  • Epigenetics Analysis and Integrated Analysis of Multiomics Data, Including Epigenetic Data, Using Artificial Intelligence in the Era of Precision Medicine
    biomolecules Review Epigenetics Analysis and Integrated Analysis of Multiomics Data, Including Epigenetic Data, Using Artificial Intelligence in the Era of Precision Medicine Ryuji Hamamoto 1,2,*, Masaaki Komatsu 1,2, Ken Takasawa 1,2 , Ken Asada 1,2 and Syuzo Kaneko 1 1 Division of Molecular Modification and Cancer Biology, National Cancer Center Research Institute, 5-1-1 Tsukiji, Chuo-ku, Tokyo 104-0045, Japan; [email protected] (M.K.); [email protected] (K.T.); [email protected] (K.A.); [email protected] (S.K.) 2 Cancer Translational Research Team, RIKEN Center for Advanced Intelligence Project, 1-4-1 Nihonbashi, Chuo-ku, Tokyo 103-0027, Japan * Correspondence: [email protected]; Tel.: +81-3-3547-5271 Received: 1 December 2019; Accepted: 27 December 2019; Published: 30 December 2019 Abstract: To clarify the mechanisms of diseases, such as cancer, studies analyzing genetic mutations have been actively conducted for a long time, and a large number of achievements have already been reported. Indeed, genomic medicine is considered the core discipline of precision medicine, and currently, the clinical application of cutting-edge genomic medicine aimed at improving the prevention, diagnosis and treatment of a wide range of diseases is promoted. However, although the Human Genome Project was completed in 2003 and large-scale genetic analyses have since been accomplished worldwide with the development of next-generation sequencing (NGS), explaining the mechanism of disease onset only using genetic variation has been recognized as difficult. Meanwhile, the importance of epigenetics, which describes inheritance by mechanisms other than the genomic DNA sequence, has recently attracted attention, and, in particular, many studies have reported the involvement of epigenetic deregulation in human cancer.
    [Show full text]
  • CUT&Runtools: a Flexible Pipeline for CUT&RUN Processing and Footprint Analysis
    Zhu et al. Genome Biology (2019) 20:192 https://doi.org/10.1186/s13059-019-1802-4 SOFTWARE Open Access CUT&RUNTools: a flexible pipeline for CUT&RUN processing and footprint analysis Qian Zhu1, Nan Liu2, Stuart H. Orkin2,3* and Guo-Cheng Yuan1* Abstract We introduce CUT&RUNTools as a flexible, general pipeline for facilitating the identification of chromatin-associated protein binding and genomic footprinting analysis from antibody-targeted CUT&RUN primary cleavage data. CUT&RUNTools extracts endonuclease cut site information from sequences of short-read fragments and produces single-locus binding estimates, aggregate motif footprints, and informative visualizations to support the high-resolution mapping capability of CUT&RUN. CUT&RUNTools is available at https://bitbucket.org/qzhudfci/cutruntools/. Introduction Results Mapping the occupancy of DNA-associated proteins, in- Overview cluding transcription factors (TFs) and histones, is central CUT&RUNTools takes paired-end sequencing read to determining cellular regulatory circuits. Conventional FASTQ files as the input and performs a set of analytical ChIP sequencing (ChIP-seq) relies on the cross-linking of steps: trimming of adapter sequences, alignment to the target proteins to DNA and physical fragmentation of reference genome, peak calling, estimation of cut matrix chromatin [1]. In practice, epitope masking and insolubil- at single-nucleotide resolution, de novo motif searching, ity of protein complexes may interfere with the successful motif footprinting analysis, direct binding site identifica- use of conventional ChIP-seq for some chromatin- tion, and data visualization (Fig. 1b). The outputs of the associated proteins [2–4]. CUT&RUN is a recently de- pipeline (Fig. 1c) are (1) an aggregate footprint capturing scribed native endonuclease-based method based on the the characteristics of chromatin-associated protein bind- binding of an antibody to a chromatin-associated protein ing (Fig.
    [Show full text]
  • Differential Chromatin Binding of the Lung Lineage Transcription Factor NKX2-1 Resolves Opposing Murine Alveolar Cell Fates in Vivo ✉ Danielle R
    ARTICLE https://doi.org/10.1038/s41467-021-22817-6 OPEN Differential chromatin binding of the lung lineage transcription factor NKX2-1 resolves opposing murine alveolar cell fates in vivo ✉ Danielle R. Little1,2, Anne M. Lynch1,3, Yun Yan 2, Haruhiko Akiyama4, Shioko Kimura 5 & Jichao Chen 1 Differential transcription of identical DNA sequences leads to distinct tissue lineages and then multiple cell types within a lineage, an epigenetic process central to progenitor and stem 1234567890():,; cell biology. The associated genome-wide changes, especially in native tissues, remain insufficiently understood, and are hereby addressed in the mouse lung, where the same lineage transcription factor NKX2-1 promotes the diametrically opposed alveolar type 1 (AT1) and AT2 cell fates. Here, we report that the cell-type-specific function of NKX2-1 is attributed to its differential chromatin binding that is acquired or retained during development in coordination with partner transcriptional factors. Loss of YAP/TAZ redirects NKX2-1 from its AT1-specific to AT2-specific binding sites, leading to transcriptionally exaggerated AT2 cells when deleted in progenitors or AT1-to-AT2 conversion when deleted after fate commitment. Nkx2-1 mutant AT1 and AT2 cells gain distinct chromatin accessible sites, including those specific to the opposite fate while adopting a gastrointestinal fate, suggesting an epigenetic plasticity unexpected from transcriptional changes. Our genomic analysis of single or purified cells, coupled with precision genetics, provides an epigenetic basis for alveolar cell fate and potential, and introduces an experimental benchmark for deciphering the in vivo function of lineage transcription factors. 1 Department of Pulmonary Medicine, The University of Texas MD Anderson Cancer Center, Houston, TX, USA.
    [Show full text]
  • Chip Sequencing
    DNBseqTM SERVICE OVERVIEW ChIP Sequencing Service Description Project Workflow ChIP Sequencing is widely used to analyze protein interactions with DNA. We care for your samples from the start through It combines chromatin immunoprecipitation (ChIP) with massively parallel to the result reporting. Highly experienced DNA sequencing to identify binding sites of DNA-associated proteins, laboratory professionals follow strict quality and can be used to precisely map global binding sites for any protein of procedures to ensure high quality results. abundant information in comparison with array-based ChIP-chip. SAMPLE PREPARATION and customized bioinformatics services to suit your specific research needs. Sample QC Confidence: Great correlation with HiSeq data. Low input: As low as 5ng ChIP-ed DNA/sample for human sample. Comprehensive analysis: Correlation analysis between ChIP-Seq and RNA-seq. LIBRARY CONSTRUCTION Sequencing Service Specification The BGI ChIP-Seq Service will be performed with BGI’s DNBseq Library QC sequencing technology, featuring combinatorial probe-anchor synthesis (cPAS) and DNA Nanoballs (DNB) technology[1] for superior data quality. SEQUENCING • 50bp Single-end sequencing reads • Standard output 20 Million reads per sample Sequencing QC • Clean data and bioinformatics analysis are available in standard file formats • Available data storage and bioinformatics RAW DATA OUTPUT applications • Cloud-based data storage and delivery system Data QC Turn Around Time BIOINFORMATICS • Typical 25 working days from sample QC acceptance ANALYSIS to final data delivery • Expedited services are available, contact your local BGI specialist for details Delivery QC DNBseqTM ChIP-Seq Data Analysis In addition to clean data output, BGI oers a range of standard, advanced and customized bioinformatics pipelines for your ChIP-Seq analysis project, including the correlation analysis of dierential expression genes and peak-related genes.
    [Show full text]
  • REVIEW Genome-Wide Identification of DNA–Protein Interactions Using Chromatin Immunoprecipitation Coupled with Flow Cell Seque
    1 REVIEW Genome-wide identification of DNA–protein interactions using chromatin immunoprecipitation coupled with flow cell sequencing Brad G Hoffman and Steven J M Jones1 Department of Cancer Endocrinology, BC Cancer Research Center, 675 West 10th Avenue, Vancouver, BC, Canada V5Z 1L3 1Micheal Smith Genome Sciences Centre, BC Cancer Agency, Suite 100-570 West 7th Avenue, Vancouver, BC, Canada V5Z 4S6 (Correspondence should be addressed to S J M Jones; Email: [email protected]) Abstract The transcriptional networks underlying mammalian cell performed, the information generated by ChIP-seq is development and function are largely unknown. The recently expected to allow the development of a framework for described use of flow cell sequencing devices in combination networks describing the transcriptional regulation of cellular with chromatin immunoprecipitation (ChIP-seq) stands to development and function. However, to date, this technology revolutionize the identification of DNA–protein interactions. has been applied only to a small number of cell types, and even As such, ChIP-seq is rapidly becoming the method of choice fewer tissues, suggesting a huge potential for novel discovery for the genome-wide localization of histone modifications in this field. and transcription factor binding sites. As further studies are Journal of Endocrinology (2009) 201, 1–13 Introduction sheared chromatin. An antibody specific to the protein of interest is then added to the sonicated material and DNA The transcriptional networks driving mammalian cell fragments bound to the protein of interest isolated via development and function are only beginning to be immunoprecipitation. DNA fragments are then released by elucidated. In many tissues transcription factors critical to reversing the cross-links and the fragments purified.
    [Show full text]
  • GETTING STARTED with Chip-SEQ
    GETTING STARTED WITH ChIP-SEQ INTRODUCTION TO ChIP SEQ IN THIS GUIDE WE WILL INTRODUCE Chromatin-immunoprecipitation (ChIP) followed by ChIP SEQUENCING AND OUTLINE KEY sequencing of the immuno-precipitated DNA is a powerful STEPS OF THE EXPERIMENTAL PROCESS tool for the investigation of Protein:DNA interactions. To INCLUDING: perform ChIP-seq, chromatin is isolated from cells or tissues and fragmented. Antibodies against chromatin associated • Experimental design proteins are used to enrich for specific chromatin fragments. • Controls for ChIP-seq experiments The DNA is recovered, sequenced and aligned to a reference • Reference genome alignment of ChIP-seq reads genome to determine specific protein binding loci. ChIP (mapping) studies have increased our knowledge of transcription factor • Background estimation biology, DNA methylation and histone modifications. • Peak finding • Quality control of ChIP-seq experiments ChIP-seq was first described in 2007 (1). ChIP sequencing • Differential binding analysis (and also microRNA sequencing) was one of the first • Motif analysis methods to make use of the power of massively parallel or next-generation sequencing (NGS) to significantly ChIP-seq may have evolved from microarray analysis but it advance real-time PCR and array-based methods. ChIP- required a completely new set of analysis tools to make the seq is a counting assay that uses only short reads to align most of the platform. ChIP-seq analysis begins with mapping to the genome, but requires millions of them to provide of trimmed sequence reads to a reference genome. Next, meaningful data. Fortunately the Solexa 1G NGS gave up peaks are found using peak-calling algorithms. To further to 30M 21-35bp reads per run.
    [Show full text]
  • Histone H3k27ac Antibody, SNAP-Chip® Certified
    Histone H3K27ac Antibody, SNAP-ChIP® Certified Catalog No. 13-0045 Lot No. 21096003-01 Pack Size 50 µg Type Monoclonal Target Size 15 kDa Reactivity H, M, WR Host Mouse Format Aff. Pur. IgG Appl. ChIP, ChIP-Seq, L, IF Product Description: This antibody meets EpiCypher’s “SNAP-ChIP® Certified” criteria for specificity and efficient target enrichment in a ChIP experiment (<20% cross-reactivity across the panel, >5% recovery of target input). Histone H3 is one of the four proteins that are present in the nucleosome, the basic repeating subunit of chromatin, consisting of 147 base pairs of DNA wrapped around an octamer of core histone proteins (H2A, H2B, H3 and H4). This antibody reacts to H3K27ac and no cross reactivity Luminex Data: Histone H3K27ac antibody was assessed using a TM with other lysine acylations in the EpiCypher SNAP-ChIP K- Luminex® based approach employing dCypher Nucleosome K- AcylStat Panel (EpiCypher Catalog No. 16-9003). The panel AcylStat panel is detected. Antibody binding to H3K27ac in the comprises biotinylated designer nucleosomes (x-axis) individually context of phosphorylation at S28 (H3K27acS28ph) is inhibited coupled to uniquely identifiable Luminex MagPlex® beads. to varying degrees in Luminex and ChIP (see figures, right). Antibody binding to nucleosomes was tested in multiplex (23-plex) at a 1:4000 dilution, and detected with second layer anti-IgG*PE. Immunogen: Data was generated using a Luminex FlexMAP3D®. Data A synthetic peptide corresponding to histone H3 acetylated at normalized to relevant on-target (H3K27ac; set to 100) is shown. lysine 27. Formulation: Protein A affinity-purified antibody (1 mg/mL) in PBS pH 7.4, with 0.05% sodium azide.
    [Show full text]
  • Cpg Islands Influence Chromatin Structure Via the Cpg-Binding Protein Cfp1
    Vol 464 | 15 April 2010 | doi:10.1038/nature08924 LETTERS CpG islands influence chromatin structure via the CpG-binding protein Cfp1 John P. Thomson1*, Peter J. Skene1*, Jim Selfridge1, Thomas Clouaire1, Jacky Guy1, Shaun Webb1, Alastair R. W. Kerr1, Aime´e Deaton1, Rob Andrews2, Keith D. James2, Daniel J. Turner2, Robert Illingworth1 & Adrian Bird1 CpG islands (CGIs) are prominent in the mammalian genome to non-methylated CGIs. We first tested CXXC finger protein 1 owing to their GC-rich base composition and high density of (Cfp1), which binds to non-methylated CpG dinucleotides in vitro CpG dinucleotides1,2. Most human gene promoters are embedded by a CXXC zinc finger domain6,12. The data showed that Cfp1 is within CGIs that lack DNA methylation and coincide with sites of enriched within the CGI fraction of the genome (Fig. 1a). histone H3 lysine 4 trimethylation (H3K4me3), irrespective of Similarly, Kdm2a, an H3K36 demethylase that also contains a transcriptional activity3,4. In spite of these intriguing correlations, CXXC domain13, was enriched in the CGI fraction. the functional significance of non-methylated CGI sequences with Focusing on Cfp1, we tested its in vivo binding specificity by chro- respect to chromatin structure and transcription is unknown. By matin immunoprecipitation (ChIP) at an endogenous CGI that is performing a search for proteins that are common to all CGIs, here present in both methylated and non-methylated states. The Xist CGI we show high enrichment for Cfp1, which selectively binds to non- is mono-allelically methylated in female cells, but fully methylated in methylated CpGs in vitro5,6.
    [Show full text]
  • Histone H3k27me3 An5body, SNAP-Chip® Cer5fied, CUTANA™ CUT&RUN Compa5ble
    Histone H3K27me3 An5body, SNAP-ChIP® Cer5fied, CUTANA™ CUT&RUN Compa5ble Catalog No. 13-0030 Lot No. 18303001 Pack Size 100 µg Type Monoclonal Target Size 15 kDa Reac5vity H, M, WR Host Rabbit Format Aff. Pur. IgG Appl. ChIP, ChIP-Seq, CUT&RUN, Luminex, WB, IHC, IF Product Descrip5on: This an,body meets EpiCypher’s “SNAP-ChIP® Cer,fied” criteria for specificity and target enrichment in ChIP (<20% cross- reac,vity to related histone post-transla,onal modifica,ons and >5% recovery of target input determined using SNAP-ChIP K-MetStat Panel spike-in controls; EpiCypher Catalog No. 19-1001). Although its specificity in CUT&RUN has yet to be empirically determined in situ using spike-in controls, CUT&RUN data produced by this an,body shows a genome wide enrichment paZern characteris,c of H3K27me3 and is Figure 1: Luminex mul5plexed specificity profiling. H3K27me3 an,body highly correlated with ChIP-seq (Figures 3-5). was assessed using a Luminex® based approach employing dCypher® Nucleosome K-MetStat Panel (EpiCypher Catalog No. 16-9002). The panel Immunogen: comprises bio,nylated designer nucleosomes (x-axis) individually coupled to color coded Luminex Magplex® beads. An,body binding to the panel of A synthe,c pep,de corresponding to histone H3 16 nucleosomes was tested in mul,plex at a 1:4000 dilu,on, and detected trimethylated at lysine 27. with second layer an,-IgG*PE. Data was generated using a Luminex FlexMAP3D®. Data is normalized to target (H3K27me3; set to 100). Formula5on: Protein A affinity-purified an,body (1 mg/mL) in PBS, with 0.09% sodium azide, 1% BSA, and 50% glycerol.
    [Show full text]
  • Siq-Chip: a Reverse-Engineered Quantitative Framework for Chip-Sequencing
    bioRxiv preprint doi: https://doi.org/10.1101/672220; this version posted June 15, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. Dickson et al. RESEARCH siQ-ChIP: A reverse-engineered quantitative framework for ChIP-sequencing Bradley M Dickson*, Rochelle L Tiedemann, Alison A Chomiak, Robert M Vaughan, Evan M Cornett and Scott B Rothbart ingful “y-axis” to ChIP-seq datasets[2, 3, 4, 5, 6, 7], calls to arms have been issued[6, 8, 9, 10], and addi- tional ChIP methods presenting newer solutions are in- troduced regularly[11]. Herein, we apply physics-based Abstract mathematical modeling to derive a quantitative frame- Chromatin immunoprecipitation followed by work for ChIP. In this attempt to shake the blues, we next-generation sequencing (ChIP-seq) is a key demonstrate that ChIP-seq is (and always has been) technique for mapping the distribution and relative quantitative without the need to modify standard- abundance of histone posttranslational ized sample preparation pipelines[12] with “spike-in” modifications (PTMs) and chromatin-associated reagents. factors across genomes. There is a perceived Fundamentally, the IP step of ChIP is a competitive challenge regarding the ability to quantitatively binding reaction. As such, this step can be described plot ChIP-seq data, and as such, approaches through biophysical models used to define binding con- making use of exogenous additives, or “spike-ins” stants and explain competition reactions.[13] Upon de- have recently been developed.
    [Show full text]
  • Bisulfite Sequencing of Chromatin Immunoprecipitated DNA (Bischip-Seq) Directly Informs Methylation Status of Histone-Modified DNA
    Downloaded from genome.cshlp.org on September 28, 2021 - Published by Cold Spring Harbor Laboratory Press Method Bisulfite sequencing of chromatin immunoprecipitated DNA (BisChIP-seq) directly informs methylation status of histone-modified DNA Aaron L. Statham,1,6 Mark D. Robinson,1,2,3,6 Jenny Z. Song,1 Marcel W. Coolen,1,4 Clare Stirzaker,1,5,7 and Susan J. Clark1,5,7,8 1Cancer Program, Garvan Institute of Medical Research, Sydney 2010, New South Wales, Australia; 2Bioinformatics Division, Walter and Eliza Hall Institute, Parkville 3052, Victoria, Australia; 3Department of Medical Biology, University of Melbourne, Parkville 3050, Victoria, Australia; 4Department of Human Genetics, Nijmegen Centre for Molecular Life Sciences (NCMLS), Radboud University Nijmegen Medical Centre, 6500 HB, Nijmegen, The Netherlands; 5St. Vincent’s Clinical School, University of NSW, Sydney 2010, New South Wales, Australia The complex relationship between DNA methylation, chromatin modification, and underlying DNA sequence is often difficult to unravel with existing technologies. Here, we describe a novel technique based on high-throughput sequencing of bisulfite-treated chromatin immunoprecipitated DNA (BisChIP-seq), which can directly interrogate genetic and epi- genetic processes that occur in normal and diseased cells. Unlike most previous reports based on correlative techniques, we found using direct bisulfite sequencing of Polycomb H3K27me3-enriched DNA from normal and prostate cancer cells that DNA methylation and H3K27me3-marked histones are not always mutually exclusive, but can co-occur in a genomic region-dependent manner. Notably, in cancer, the co-dependency of marks is largely redistributed with an increase of the dual repressive marks at CpG islands and transcription start sites of silent genes.
    [Show full text]
  • 4D Nucleomes in Single Cells: What Can Computational Modeling Reveal About Spatial Chromatin Conformation? Monika Sekelja, Jonas Paulsen and Philippe Collas*
    Sekelja et al. Genome Biology (2016) 17:54 DOI 10.1186/s13059-016-0923-2 REVIEW Open Access 4D nucleomes in single cells: what can computational modeling reveal about spatial chromatin conformation? Monika Sekelja, Jonas Paulsen and Philippe Collas* frequencies are drawn to each other in the nuclear space. Abstract To improve the accuracy of 3D models of chromatin, Genome-wide sequencing technologies enable other constraints can potentially be incorporated into investigations of the structural properties of the structural models based on association of chromatin with genome in various spatial dimensions. Here, we known anchors in the nucleus, such as the nuclear enve- review computational techniques developed to lope [4, 12], nuclear pore complexes [13, 14], or nucleoli model the three-dimensional genome in single cells [15, 16]. versus ensembles of cells and assess their underlying Most 3D genome reconstructions are performed on cell assumptions. We further address approaches to study population-averaged Hi-C contact matrices [6, 8, 17–23]. the spatio-temporal aspects of genome organization The results consistently provide a hierarchical view of from single-cell data. folding of the genome, with chromatin divided into supra- megabase compartments of transcriptionally active or in- active chromatin (the so-called A and B compartments) Background [6, 8] and, within these compartments, megabase-scale Increasing evidence indicates that the spatial, three- topologically associated domains (TADs) [7, 24, 25]. TADs dimensional (3D) organization of chromatin influences show distinct boundaries, within which loci interact more gene expression and cell fate [1–8]. Chromosome frequently with one another than with loci of adjacent conformation capture (3C) techniques coupled with high- TADs.
    [Show full text]