1
|
Beagrie RA, Thieme CJ, Annunziatella C, Baugher C, Zhang Y, Schueler M, Kukalev A, Kempfer R, Chiariello AM, Bianco S, Li Y, Davis T, Scialdone A, Welch LR, Nicodemi M, Pombo A. Multiplex-GAM: genome-wide identification of chromatin contacts yields insights overlooked by Hi-C. Nat Methods 2023:10.1038/s41592-023-01903-1. [PMID: 37336949 PMCID: PMC10333126 DOI: 10.1038/s41592-023-01903-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2021] [Accepted: 05/01/2023] [Indexed: 06/21/2023]
Abstract
Technology for measuring 3D genome topology is increasingly important for studying gene regulation, for genome assembly and for mapping of genome rearrangements. Hi-C and other ligation-based methods have become routine but have specific biases. Here, we develop multiplex-GAM, a faster and more affordable version of genome architecture mapping (GAM), a ligation-free technique that maps chromatin contacts genome-wide. We perform a detailed comparison of multiplex-GAM and Hi-C using mouse embryonic stem cells. When examining the strongest contacts detected by either method, we find that only one-third of these are shared. The strongest contacts specifically found in GAM often involve 'active' regions, including many transcribed genes and super-enhancers, whereas in Hi-C they more often contain 'inactive' regions. Our work shows that active genomic regions are involved in extensive complex contacts that are currently underestimated in ligation-based approaches, and highlights the need for orthogonal advances in genome-wide contact mapping technologies.
Collapse
Affiliation(s)
- Robert A Beagrie
- Max-Delbrück-Center for Molecular Medicine in the Helmholtz Association (MDC), Berlin Institute for Medical Systems Biology (BIMSB), Epigenetic Regulation and Chromatin Architecture Group, Berlin, Germany
- Laboratory of Gene Regulation, Weatherall Institute of Molecular Medicine, Oxford, UK
- Chromatin and Disease Group, Wellcome Centre for Human Genetics, Oxford, UK
| | - Christoph J Thieme
- Max-Delbrück-Center for Molecular Medicine in the Helmholtz Association (MDC), Berlin Institute for Medical Systems Biology (BIMSB), Epigenetic Regulation and Chromatin Architecture Group, Berlin, Germany
| | - Carlo Annunziatella
- Dipartimento di Fisica, Università di Napoli Federico II, and INFN Napoli, CNR-SPIN, Complesso Universitario di Monte Sant'Angelo, Naples, Italy
| | - Catherine Baugher
- School of Electrical Engineering and Computer Science, Ohio University, Athens, OH, USA
| | - Yingnan Zhang
- School of Electrical Engineering and Computer Science, Ohio University, Athens, OH, USA
| | - Markus Schueler
- Max-Delbrück-Center for Molecular Medicine in the Helmholtz Association (MDC), Berlin Institute for Medical Systems Biology (BIMSB), Epigenetic Regulation and Chromatin Architecture Group, Berlin, Germany
| | - Alexander Kukalev
- Max-Delbrück-Center for Molecular Medicine in the Helmholtz Association (MDC), Berlin Institute for Medical Systems Biology (BIMSB), Epigenetic Regulation and Chromatin Architecture Group, Berlin, Germany
| | - Rieke Kempfer
- Max-Delbrück-Center for Molecular Medicine in the Helmholtz Association (MDC), Berlin Institute for Medical Systems Biology (BIMSB), Epigenetic Regulation and Chromatin Architecture Group, Berlin, Germany
- Humboldt-Universität zu Berlin, Berlin, Germany
| | - Andrea M Chiariello
- Dipartimento di Fisica, Università di Napoli Federico II, and INFN Napoli, CNR-SPIN, Complesso Universitario di Monte Sant'Angelo, Naples, Italy
| | - Simona Bianco
- Max-Delbrück-Center for Molecular Medicine in the Helmholtz Association (MDC), Berlin Institute for Medical Systems Biology (BIMSB), Epigenetic Regulation and Chromatin Architecture Group, Berlin, Germany
- Dipartimento di Fisica, Università di Napoli Federico II, and INFN Napoli, CNR-SPIN, Complesso Universitario di Monte Sant'Angelo, Naples, Italy
| | - Yichao Li
- School of Electrical Engineering and Computer Science, Ohio University, Athens, OH, USA
| | - Trenton Davis
- School of Electrical Engineering and Computer Science, Ohio University, Athens, OH, USA
| | - Antonio Scialdone
- Institute of Epigenetics and Stem Cells, Helmholtz Zentrum München - German Research Center for Environmental Health, Munich, Germany
- Institute of Functional Epigenetics, Helmholtz Zentrum München - German Research Center for Environmental Health, Neuherberg, Germany
- Institute of Computational Biology, Helmholtz Zentrum München - German Research Center for Environmental Health, Neuherberg, Germany
| | - Lonnie R Welch
- School of Electrical Engineering and Computer Science, Ohio University, Athens, OH, USA.
| | - Mario Nicodemi
- Max-Delbrück-Center for Molecular Medicine in the Helmholtz Association (MDC), Berlin Institute for Medical Systems Biology (BIMSB), Epigenetic Regulation and Chromatin Architecture Group, Berlin, Germany.
- Dipartimento di Fisica, Università di Napoli Federico II, and INFN Napoli, CNR-SPIN, Complesso Universitario di Monte Sant'Angelo, Naples, Italy.
- Berlin Institute of Health (BIH), MDC-Berlin, Berlin, Germany.
| | - Ana Pombo
- Max-Delbrück-Center for Molecular Medicine in the Helmholtz Association (MDC), Berlin Institute for Medical Systems Biology (BIMSB), Epigenetic Regulation and Chromatin Architecture Group, Berlin, Germany.
- Humboldt-Universität zu Berlin, Berlin, Germany.
| |
Collapse
|
2
|
Weigert R, Hetzel S, Bailly N, Haggerty C, Ilik IA, Yung PYK, Navarro C, Bolondi A, Kumar AS, Anania C, Brändl B, Meierhofer D, Lupiáñez DG, Müller FJ, Aktas T, Elsässer SJ, Kretzmer H, Smith ZD, Meissner A. Dynamic antagonism between key repressive pathways maintains the placental epigenome. Nat Cell Biol 2023; 25:579-591. [PMID: 37024684 PMCID: PMC10104784 DOI: 10.1038/s41556-023-01114-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2022] [Accepted: 02/21/2023] [Indexed: 04/08/2023]
Abstract
DNA and Histone 3 Lysine 27 methylation typically function as repressive modifications and operate within distinct genomic compartments. In mammals, the majority of the genome is kept in a DNA methylated state, whereas the Polycomb repressive complexes regulate the unmethylated CpG-rich promoters of developmental genes. In contrast to this general framework, the extra-embryonic lineages display non-canonical, globally intermediate DNA methylation levels, including disruption of local Polycomb domains. Here, to better understand this unusual landscape's molecular properties, we genetically and chemically perturbed major epigenetic pathways in mouse trophoblast stem cells. We find that the extra-embryonic epigenome reflects ongoing and dynamic de novo methyltransferase recruitment, which is continuously antagonized by Polycomb to maintain intermediate, locally disordered methylation. Despite its disorganized molecular appearance, our data point to a highly controlled equilibrium between counteracting repressors within extra-embryonic cells, one that can seemingly persist indefinitely without bistable features typically seen for embryonic forms of epigenetic regulation.
Collapse
Affiliation(s)
- Raha Weigert
- Department of Genome Regulation, Max Planck Institute for Molecular Genetics, Berlin, Germany
- Department of Medical Biotechnology, Technische Universität Berlin, Berlin, Germany
| | - Sara Hetzel
- Department of Genome Regulation, Max Planck Institute for Molecular Genetics, Berlin, Germany
| | - Nina Bailly
- Department of Genome Regulation, Max Planck Institute for Molecular Genetics, Berlin, Germany
- Department of Biology, Chemistry and Pharmacy, Freie Universität Berlin, Berlin, Germany
| | - Chuck Haggerty
- Department of Genome Regulation, Max Planck Institute for Molecular Genetics, Berlin, Germany
- Department of Biology, Chemistry and Pharmacy, Freie Universität Berlin, Berlin, Germany
| | - Ibrahim A Ilik
- Otto Warburg Laboratories, Max Planck Institute for Molecular Genetics, Berlin, Germany
| | - Philip Yuk Kwong Yung
- Science for Life Laboratory, Department of Medical Biochemistry and Biophysics, Karolinska Institutet, Stockholm, Sweden
- Ming Wai Lau Centre for Reparative Medicine, Stockholm node, Karolinska Institutet, Stockholm, Sweden
| | - Carmen Navarro
- Science for Life Laboratory, Department of Medical Biochemistry and Biophysics, Karolinska Institutet, Stockholm, Sweden
- Ming Wai Lau Centre for Reparative Medicine, Stockholm node, Karolinska Institutet, Stockholm, Sweden
| | - Adriano Bolondi
- Department of Genome Regulation, Max Planck Institute for Molecular Genetics, Berlin, Germany
- Department of Biology, Chemistry and Pharmacy, Freie Universität Berlin, Berlin, Germany
| | - Abhishek Sampath Kumar
- Department of Genome Regulation, Max Planck Institute for Molecular Genetics, Berlin, Germany
| | - Chiara Anania
- Epigenetics and Sex Development Group, Berlin Institute for Medical Systems Biology, Max-Delbrück Center for Molecular Medicine, Berlin-Buch, Germany
| | - Björn Brändl
- Department of Genome Regulation, Max Planck Institute for Molecular Genetics, Berlin, Germany
- Universitätsklinikum Schleswig-Holstein Campus Kiel, Zentrum für Integrative Psychiatrie gGmbH, Kiel, Germany
| | - David Meierhofer
- Mass Spectrometry Joint Facilities Scientific Service, Max Planck Institute for Molecular Genetics, Berlin, Germany
| | - Darío G Lupiáñez
- Epigenetics and Sex Development Group, Berlin Institute for Medical Systems Biology, Max-Delbrück Center for Molecular Medicine, Berlin-Buch, Germany
| | - Franz-Josef Müller
- Department of Genome Regulation, Max Planck Institute for Molecular Genetics, Berlin, Germany
- Universitätsklinikum Schleswig-Holstein Campus Kiel, Zentrum für Integrative Psychiatrie gGmbH, Kiel, Germany
| | - Tugce Aktas
- Otto Warburg Laboratories, Max Planck Institute for Molecular Genetics, Berlin, Germany
| | - Simon J Elsässer
- Science for Life Laboratory, Department of Medical Biochemistry and Biophysics, Karolinska Institutet, Stockholm, Sweden
- Ming Wai Lau Centre for Reparative Medicine, Stockholm node, Karolinska Institutet, Stockholm, Sweden
| | - Helene Kretzmer
- Department of Genome Regulation, Max Planck Institute for Molecular Genetics, Berlin, Germany
| | - Zachary D Smith
- Department of Genetics, Yale Stem Cell Center, Yale School of Medicine, New Haven, CT, USA.
| | - Alexander Meissner
- Department of Genome Regulation, Max Planck Institute for Molecular Genetics, Berlin, Germany.
- Department of Biology, Chemistry and Pharmacy, Freie Universität Berlin, Berlin, Germany.
- Broad Institute of MIT and Harvard, Cambridge, MA, USA.
- Department of Stem Cell and Regenerative Biology, Harvard University, Cambridge, MA, US.
| |
Collapse
|
3
|
Naik M, Rueda L, Vasighizaker A. Identification of Enriched Regions in ChIP-Seq Data via a Linear-Time Multi-Level Thresholding Algorithm. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2022; 19:2842-2850. [PMID: 34398762 DOI: 10.1109/tcbb.2021.3104734] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Chromatin immunoprecipitation (ChIP-Seq) has emerged as a superior alternative to microarray technology as it provides higher resolution, less noise, greater coverage and wider dynamic range. While ChIP-Seq enables probing of DNA-protein interaction over the entire genome, it requires the use of sophisticated tools to recognize hidden patterns and extract meaningful data. Over the years, various attempts have resulted in several algorithms making use of different heuristics to accurately determine individual peaks corresponding to unique DNA-protein. However, finding all the significant peaks with high accuracy in a reasonable time is still a challenge. In this work, we propose the use of Multi-level thresholding algorithm, which we call LinMLTBS, used to identify the enriched regions on ChIP-Seq data. Although various suboptimal heuristics have been proposed for multi-level thresholding, we emphasize on the use of an algorithm capable of obtaining an optimal solution, while maintaining linear-time complexity. Testing various algorithm on various ENCODE project datasets shows that our approach attains higher accuracy relative to previously proposed peak finders while retaining a reasonable processing speed.
Collapse
|
4
|
Jing F, Zhang SW, Zhang S. Prediction of the transcription factor binding sites with meta-learning. Methods 2022; 203:207-213. [DOI: 10.1016/j.ymeth.2022.04.010] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2021] [Revised: 04/01/2022] [Accepted: 04/17/2022] [Indexed: 11/26/2022] Open
|
5
|
Gonzalez-Teran B, Pittman M, Felix F, Thomas R, Richmond-Buccola D, Hüttenhain R, Choudhary K, Moroni E, Costa MW, Huang Y, Padmanabhan A, Alexanian M, Lee CY, Maven BEJ, Samse-Knapp K, Morton SU, McGregor M, Gifford CA, Seidman JG, Seidman CE, Gelb BD, Colombo G, Conklin BR, Black BL, Bruneau BG, Krogan NJ, Pollard KS, Srivastava D. Transcription factor protein interactomes reveal genetic determinants in heart disease. Cell 2022; 185:794-814.e30. [PMID: 35182466 PMCID: PMC8923057 DOI: 10.1016/j.cell.2022.01.021] [Citation(s) in RCA: 34] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2020] [Revised: 08/20/2021] [Accepted: 01/25/2022] [Indexed: 02/08/2023]
Abstract
Congenital heart disease (CHD) is present in 1% of live births, yet identification of causal mutations remains challenging. We hypothesized that genetic determinants for CHDs may lie in the protein interactomes of transcription factors whose mutations cause CHDs. Defining the interactomes of two transcription factors haplo-insufficient in CHD, GATA4 and TBX5, within human cardiac progenitors, and integrating the results with nearly 9,000 exomes from proband-parent trios revealed an enrichment of de novo missense variants associated with CHD within the interactomes. Scoring variants of interactome members based on residue, gene, and proband features identified likely CHD-causing genes, including the epigenetic reader GLYR1. GLYR1 and GATA4 widely co-occupied and co-activated cardiac developmental genes, and the identified GLYR1 missense variant disrupted interaction with GATA4, impairing in vitro and in vivo function in mice. This integrative proteomic and genetic approach provides a framework for prioritizing and interrogating genetic variants in heart disease.
Collapse
Affiliation(s)
- Barbara Gonzalez-Teran
- Gladstone Institutes, San Francisco, CA, USA; Roddenberry Center for Stem Cell Biology and Medicine at Gladstone, San Francisco, CA, USA
| | - Maureen Pittman
- Gladstone Institutes, San Francisco, CA, USA; Department of Epidemiology & Biostatistics, Institute for Computational Health Sciences, and Institute for Human Genetics, University of California San Francisco, San Francisco, CA, USA
| | - Franco Felix
- Gladstone Institutes, San Francisco, CA, USA; Roddenberry Center for Stem Cell Biology and Medicine at Gladstone, San Francisco, CA, USA
| | | | - Desmond Richmond-Buccola
- Gladstone Institutes, San Francisco, CA, USA; Roddenberry Center for Stem Cell Biology and Medicine at Gladstone, San Francisco, CA, USA
| | - Ruth Hüttenhain
- Gladstone Institutes, San Francisco, CA, USA; Department of Cellular and Molecular Pharmacology, University of California San Francisco, San Francisco, CA, USA; Quantitative Biosciences Institute (QBI), University of California San Francisco, San Francisco, CA, USA
| | | | | | - Mauro W Costa
- Gladstone Institutes, San Francisco, CA, USA; Roddenberry Center for Stem Cell Biology and Medicine at Gladstone, San Francisco, CA, USA
| | - Yu Huang
- Gladstone Institutes, San Francisco, CA, USA; Roddenberry Center for Stem Cell Biology and Medicine at Gladstone, San Francisco, CA, USA
| | - Arun Padmanabhan
- Gladstone Institutes, San Francisco, CA, USA; Roddenberry Center for Stem Cell Biology and Medicine at Gladstone, San Francisco, CA, USA; Division of Cardiology, Department of Medicine, University of California, San Francisco, CA, USA
| | - Michael Alexanian
- Gladstone Institutes, San Francisco, CA, USA; Roddenberry Center for Stem Cell Biology and Medicine at Gladstone, San Francisco, CA, USA
| | - Clara Youngna Lee
- Gladstone Institutes, San Francisco, CA, USA; Roddenberry Center for Stem Cell Biology and Medicine at Gladstone, San Francisco, CA, USA
| | - Bonnie E J Maven
- Gladstone Institutes, San Francisco, CA, USA; Roddenberry Center for Stem Cell Biology and Medicine at Gladstone, San Francisco, CA, USA; Developmental and Stem Cell Biology Graduate Program, University of California San Francisco, San Francisco, CA, USA
| | - Kaitlen Samse-Knapp
- Gladstone Institutes, San Francisco, CA, USA; Roddenberry Center for Stem Cell Biology and Medicine at Gladstone, San Francisco, CA, USA
| | - Sarah U Morton
- Division of Newborn Medicine, Department of Medicine, Boston Children's Hospital, Boston, MA, USA; Department of Pediatrics, Harvard Medical School, Boston, MA, USA
| | - Michael McGregor
- Gladstone Institutes, San Francisco, CA, USA; Department of Cellular and Molecular Pharmacology, University of California San Francisco, San Francisco, CA, USA; Quantitative Biosciences Institute (QBI), University of California San Francisco, San Francisco, CA, USA
| | - Casey A Gifford
- Gladstone Institutes, San Francisco, CA, USA; Roddenberry Center for Stem Cell Biology and Medicine at Gladstone, San Francisco, CA, USA
| | - J G Seidman
- Department of Genetics, Harvard Medical School, Boston, MA, USA
| | - Christine E Seidman
- Department of Genetics, Harvard Medical School, Boston, MA, USA; Howard Hughes Medical Institute, Harvard Medical School, Boston, MA, USA; Cardiovascular Division, Brigham and Women's Hospital, Boston, MA, USA
| | - Bruce D Gelb
- Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, USA; Mindich Child Health and Development Institute, Icahn School of Medicine at Mount Sinai, New York, NY, USA; Department of Pediatrics, Icahn School of Medicine at Mount Sinai, New York, NY, USA
| | | | - Bruce R Conklin
- Gladstone Institutes, San Francisco, CA, USA; Roddenberry Center for Stem Cell Biology and Medicine at Gladstone, San Francisco, CA, USA
| | - Brian L Black
- Cardiovascular Research Institute, University of California San Francisco, San Francisco, CA, USA
| | - Benoit G Bruneau
- Gladstone Institutes, San Francisco, CA, USA; Roddenberry Center for Stem Cell Biology and Medicine at Gladstone, San Francisco, CA, USA; Cardiovascular Research Institute, University of California San Francisco, San Francisco, CA, USA; Division of Cardiology, Department of Pediatrics, UCSF School of Medicine, San Francisco, CA, USA
| | - Nevan J Krogan
- Gladstone Institutes, San Francisco, CA, USA; Department of Cellular and Molecular Pharmacology, University of California San Francisco, San Francisco, CA, USA; Quantitative Biosciences Institute (QBI), University of California San Francisco, San Francisco, CA, USA
| | - Katherine S Pollard
- Gladstone Institutes, San Francisco, CA, USA; Chan Zuckerberg Biohub, San Francisco, CA, USA; Department of Epidemiology & Biostatistics, Institute for Computational Health Sciences, and Institute for Human Genetics, University of California San Francisco, San Francisco, CA, USA.
| | - Deepak Srivastava
- Gladstone Institutes, San Francisco, CA, USA; Roddenberry Center for Stem Cell Biology and Medicine at Gladstone, San Francisco, CA, USA; Division of Cardiology, Department of Pediatrics, UCSF School of Medicine, San Francisco, CA, USA; Department of Biochemistry and Biophysics, University of California San Francisco, San Francisco, CA, USA.
| |
Collapse
|
6
|
Brahma safeguards canalization of cardiac mesoderm differentiation. Nature 2022; 602:129-134. [PMID: 35082446 DOI: 10.1038/s41586-021-04336-y] [Citation(s) in RCA: 22] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2020] [Accepted: 12/08/2021] [Indexed: 12/26/2022]
Abstract
Differentiation proceeds along a continuum of increasingly fate-restricted intermediates, referred to as canalization1,2. Canalization is essential for stabilizing cell fate, but the mechanisms that underlie robust canalization are unclear. Here we show that the BRG1/BRM-associated factor (BAF) chromatin-remodelling complex ATPase gene Brm safeguards cell identity during directed cardiogenesis of mouse embryonic stem cells. Despite the establishment of a well-differentiated precardiac mesoderm, Brm-/- cells predominantly became neural precursors, violating germ layer assignment. Trajectory inference showed a sudden acquisition of a non-mesodermal identity in Brm-/- cells. Mechanistically, the loss of Brm prevented de novo accessibility of primed cardiac enhancers while increasing the expression of neurogenic factor POU3F1, preventing the binding of the neural suppressor REST and shifting the composition of BRG1 complexes. The identity switch caused by the Brm mutation was overcome by increasing BMP4 levels during mesoderm induction. Mathematical modelling supports these observations and demonstrates that Brm deletion affects cell fate trajectory by modifying saddle-node bifurcations2. In the mouse embryo, Brm deletion exacerbated mesoderm-deleted Brg1-mutant phenotypes, severely compromising cardiogenesis, and reveals an in vivo role for Brm. Our results show that Brm is a compensable safeguard of the fidelity of mesoderm chromatin states, and support a model in which developmental canalization is not a rigid irreversible path, but a highly plastic trajectory.
Collapse
|
7
|
Bayat F, Libbrecht M. VSS: Variance-stabilized signals for sequencing-based genomic signals. Bioinformatics 2021; 37:4383-4391. [PMID: 34165492 PMCID: PMC8652025 DOI: 10.1093/bioinformatics/btab457] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2020] [Revised: 04/28/2021] [Accepted: 06/17/2021] [Indexed: 11/12/2022] Open
Abstract
MOTIVATION A sequencing-based genomic assay such as ChIP-seq outputs a real-valued signal for each position in the genome that measures the strength of activity at that position. Most genomic signals lack the property of variance stabilization. That is, a difference between 0 and 100 reads usually has a very different statistical importance from a difference between 1,000 and 1,100 reads. A statistical model such as a negative binomial distribution can account for this pattern, but learning these models is computationally challenging. Therefore, many applications - including imputation and segmentation and genome annotation (SAGA) - instead use Gaussian models and use a transformation such as log or inverse hyperbolic sine (asinh) to stabilize variance. RESULTS We show here that existing transformations do not fully stabilize variance in genomic data sets. To solve this issue, we propose VSS, a method that produces variance-stabilized signals for sequencing-based genomic signals. VSS learns the empirical relationship between the mean and variance of a given signal data set and produces transformed signals that normalize for this dependence. We show that VSS successfully stabilizes variance and that doing so improves downstream applications such as SAGA. VSS will eliminate the need for downstream methods to implement complex mean-variance relationship models, and will enable genomic signals to be easily understood by eye. AVAILABILITY https://github.com/faezeh-bayat/VSS. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Faezeh Bayat
- Department of Computing Science, Simon Fraser University, Burnaby, BC V5A 1S6, Canada
| | - Maxwell Libbrecht
- Department of Computing Science, Simon Fraser University, Burnaby, BC V5A 1S6, Canada
| |
Collapse
|
8
|
Zheng A, Lamkin M, Qiu Y, Ren K, Goren A, Gymrek M. A flexible ChIP-sequencing simulation toolkit. BMC Bioinformatics 2021; 22:201. [PMID: 33879052 PMCID: PMC8056602 DOI: 10.1186/s12859-021-04097-5] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2021] [Accepted: 03/22/2021] [Indexed: 11/17/2022] Open
Abstract
Background A major challenge in evaluating quantitative ChIP-seq analyses, such as peak calling and differential binding, is a lack of reliable ground truth data. Accurate simulation of ChIP-seq data can mitigate this challenge, but existing frameworks are either too cumbersome to apply genome-wide or unable to model a number of important experimental conditions in ChIP-seq. Results We present ChIPs, a toolkit for rapidly simulating ChIP-seq data using statistical models of key experimental steps. We demonstrate how ChIPs can be used for a range of applications, including benchmarking analysis tools and evaluating the impact of various experimental parameters. ChIPs is implemented as a standalone command-line program written in C++ and is available from https://github.com/gymreklab/chips. Conclusions ChIPs is an efficient ChIP-seq simulation framework that generates realistic datasets over a flexible range of experimental conditions. It can serve as an important component in various ChIP-seq analyses where ground truth data are needed. Supplementary Information The online version contains supplementary material available at 10.1186/s12859-021-04097-5.
Collapse
Affiliation(s)
- An Zheng
- Department of Computer Science and Engineering, University of California San Diego, 9500 Gilman Drive, La Jolla, CA, 92093, USA
| | - Michael Lamkin
- Department of Bioengineering, University of California San Diego, 9500 Gilman Drive, La Jolla, CA, 92093, USA
| | - Yutong Qiu
- Department of Computer Science and Engineering, University of California San Diego, 9500 Gilman Drive, La Jolla, CA, 92093, USA.,School of Computer Science, Carnegie Mellon University, 5000 Forbes Avenue, Pittsburgh, PA, 15213, USA
| | - Kevin Ren
- Department of Mathematics, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA, 02139, USA
| | - Alon Goren
- Department of Medicine, University of California San Diego, 9500 Gilman Drive, La Jolla, CA, 92093, USA.
| | - Melissa Gymrek
- Department of Computer Science and Engineering, University of California San Diego, 9500 Gilman Drive, La Jolla, CA, 92093, USA. .,Department of Medicine, University of California San Diego, 9500 Gilman Drive, La Jolla, CA, 92093, USA.
| |
Collapse
|
9
|
Xing H, Wu Y, Zhang MQ, Chen Y. Deciphering hierarchical organization of topologically associated domains through change-point testing. BMC Bioinformatics 2021; 22:183. [PMID: 33838653 PMCID: PMC8037919 DOI: 10.1186/s12859-021-04113-8] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2020] [Accepted: 03/30/2021] [Indexed: 12/20/2022] Open
Abstract
Background The nucleus of eukaryotic cells spatially packages chromosomes into a hierarchical and distinct segregation that plays critical roles in maintaining transcription regulation. High-throughput methods of chromosome conformation capture, such as Hi-C, have revealed topologically associating domains (TADs) that are defined by biased chromatin interactions within them. Results We introduce a novel method, HiCKey, to decipher hierarchical TAD structures in Hi-C data and compare them across samples. We first derive a generalized likelihood-ratio (GLR) test for detecting change-points in an interaction matrix that follows a negative binomial distribution or general mixture distribution. We then employ several optimal search strategies to decipher hierarchical TADs with p values calculated by the GLR test. Large-scale validations of simulation data show that HiCKey has good precision in recalling known TADs and is robust against random collisions of chromatin interactions. By applying HiCKey to Hi-C data of seven human cell lines, we identified multiple layers of TAD organization among them, but the vast majority had no more than four layers. In particular, we found that TAD boundaries are significantly enriched in active chromosomal regions compared to repressed regions. Conclusions HiCKey is optimized for processing large matrices constructed from high-resolution Hi-C experiments. The method and theoretical result of the GLR test provide a general framework for significance testing of similar experimental chromatin interaction data that may not fully follow negative binomial distributions but rather more general mixture distributions. Supplementary Information The online version contains supplementary material available at 10.1186/s12859-021-04113-8.
Collapse
Affiliation(s)
- Haipeng Xing
- Department of Applied Mathematics and Statistics, State University of New York at Stony Brook, 100 Nicolls Rd, Stony Brook, NY, 11794, USA
| | - Yingru Wu
- Department of Applied Mathematics and Statistics, State University of New York at Stony Brook, 100 Nicolls Rd, Stony Brook, NY, 11794, USA
| | - Michael Q Zhang
- Center for System Biology, University of Texas at Dallas, 800 W Campbell Rd, Richardson, TX, 75080, USA
| | - Yong Chen
- Department of Molecular and Cellular Biosciences, Rowan University, 201 Mullica Hill Rd, Glassboro, NJ, 08028, USA.
| |
Collapse
|
10
|
Chitpin JG, Awdeh A, Perkins TJ. RECAP reveals the true statistical significance of ChIP-seq peak calls. Bioinformatics 2020; 35:3592-3598. [PMID: 30824903 PMCID: PMC6761936 DOI: 10.1093/bioinformatics/btz150] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2018] [Revised: 01/18/2019] [Accepted: 02/27/2019] [Indexed: 12/29/2022] Open
Abstract
Motivation Chromatin Immunopreciptation (ChIP)-seq is used extensively to identify sites of transcription factor binding or regions of epigenetic modifications to the genome. A key step in ChIP-seq analysis is peak calling, where genomic regions enriched for ChIP versus control reads are identified. Many programs have been designed to solve this task, but nearly all fall into the statistical trap of using the data twice—once to determine candidate enriched regions, and again to assess enrichment by classical statistical hypothesis testing. This double use of the data invalidates the statistical significance assigned to enriched regions, thus the true significance or reliability of peak calls remains unknown. Results Using simulated and real ChIP-seq data, we show that three well-known peak callers, MACS, SICER and diffReps, output biased P-values and false discovery rate estimates that can be many orders of magnitude too optimistic. We propose a wrapper algorithm, RECAP, that uses resampling of ChIP-seq and control data to estimate a monotone transform correcting for biases built into peak calling algorithms. When applied to null hypothesis data, where there is no enrichment between ChIP-seq and control, P-values recalibrated by RECAP are approximately uniformly distributed. On data where there is genuine enrichment, RECAP P-values give a better estimate of the true statistical significance of candidate peaks and better false discovery rate estimates, which correlate better with empirical reproducibility. RECAP is a powerful new tool for assessing the true statistical significance of ChIP-seq peak calls. Availability and implementation The RECAP software is available through www.perkinslab.ca or on github at https://github.com/theodorejperkins/RECAP. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Justin G Chitpin
- Translational and Molecular Medicine Program, University of Ottawa, Ottawa, ON K1H8M5, Canada.,Regenerative Medicine Program, Ottawa Hospital Research Institute, Ottawa, ON K1H8L6, Canada
| | - Aseel Awdeh
- Regenerative Medicine Program, Ottawa Hospital Research Institute, Ottawa, ON K1H8L6, Canada.,School of Electrical Engineering and Computer Science, University of Ottawa, Ottawa, ON K1N6N5, Canada
| | - Theodore J Perkins
- Regenerative Medicine Program, Ottawa Hospital Research Institute, Ottawa, ON K1H8L6, Canada.,School of Electrical Engineering and Computer Science, University of Ottawa, Ottawa, ON K1N6N5, Canada.,Department of Biochemistry, Microbiology and Immunology, University of Ottawa, Ottawa, ON K1H8M5, Canada
| |
Collapse
|
11
|
Höllbacher B, Balázs K, Heinig M, Uhlenhaut NH. Seq-ing answers: Current data integration approaches to uncover mechanisms of transcriptional regulation. Comput Struct Biotechnol J 2020; 18:1330-1341. [PMID: 32612756 PMCID: PMC7306512 DOI: 10.1016/j.csbj.2020.05.018] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2020] [Revised: 05/21/2020] [Accepted: 05/23/2020] [Indexed: 02/06/2023] Open
Abstract
Advancements in the field of next generation sequencing lead to the generation of ever-more data, with the challenge often being how to combine and reconcile results from different OMICs studies such as genome, epigenome and transcriptome. Here we provide an overview of the standard processing pipelines for ChIP-seq and RNA-seq as well as common downstream analyses. We describe popular multi-omics data integration approaches used to identify target genes and co-factors, and we discuss how machine learning techniques may predict transcriptional regulators and gene expression.
Collapse
Affiliation(s)
- Barbara Höllbacher
- Institute for Diabetes and Cancer IDC, Helmholtz Zentrum Muenchen (HMGU) and German Center for Diabetes Research (DZD), Munich 85764, Neuherberg, Germany.,Institute of Computational Biology ICB, Helmholtz Zentrum Muenchen (HMGU) and German Center for Diabetes Research (DZD), Munich 85764, Neuherberg, Germany.,Department of Informatics, TUM, Munich 85748, Garching, Germany
| | - Kinga Balázs
- Institute for Diabetes and Cancer IDC, Helmholtz Zentrum Muenchen (HMGU) and German Center for Diabetes Research (DZD), Munich 85764, Neuherberg, Germany
| | - Matthias Heinig
- Institute of Computational Biology ICB, Helmholtz Zentrum Muenchen (HMGU) and German Center for Diabetes Research (DZD), Munich 85764, Neuherberg, Germany.,Department of Informatics, TUM, Munich 85748, Garching, Germany
| | - N Henriette Uhlenhaut
- Institute for Diabetes and Cancer IDC, Helmholtz Zentrum Muenchen (HMGU) and German Center for Diabetes Research (DZD), Munich 85764, Neuherberg, Germany.,Metabolic Programming, TUM School of Life Sciences Weihenstephan, Munich 85354, Freising, Germany
| |
Collapse
|
12
|
Baldoni PL, Rashid NU, Ibrahim JG. Improved detection of epigenomic marks with mixed-effects hidden Markov models. Biometrics 2019; 75:1401-1413. [PMID: 31081192 PMCID: PMC6851437 DOI: 10.1111/biom.13083] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2018] [Accepted: 05/03/2019] [Indexed: 11/30/2022]
Abstract
Chromatin immunoprecipitation followed by next-generation sequencing (ChIP-seq) is a technique to detect genomic regions containing protein-DNA interaction, such as transcription factor binding sites or regions containing histone modifications. One goal of the analysis of ChIP-seq experiments is to identify genomic loci enriched for sequencing reads pertaining to DNA bound to the factor of interest. The accurate identification of such regions aids in the understanding of epigenomic marks and gene regulatory mechanisms. Given the reduction of massively parallel sequencing costs, methods to detect consensus regions of enrichment across multiple samples are of interest. Here, we present a statistical model to detect broad consensus regions of enrichment from ChIP-seq technical or biological replicates through a class of zero-inflated mixed-effects hidden Markov models. We show that the proposed model outperforms existing methods for consensus peak calling in common epigenomic marks by accounting for the excess zeros and sample-specific biases. We apply our method to data from the Encyclopedia of DNA Elements and Roadmap Epigenomics projects and also from an extensive simulation study.
Collapse
Affiliation(s)
- Pedro L. Baldoni
- Department of Biostatistics, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, U.S.A
| | - Naim U. Rashid
- Department of Biostatistics, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, U.S.A
| | - Joseph G. Ibrahim
- Department of Biostatistics, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, U.S.A
| |
Collapse
|
13
|
Gheorghe M, Sandve GK, Khan A, Chèneby J, Ballester B, Mathelier A. A map of direct TF-DNA interactions in the human genome. Nucleic Acids Res 2019; 47:e21. [PMID: 30517703 PMCID: PMC6393237 DOI: 10.1093/nar/gky1210] [Citation(s) in RCA: 48] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2018] [Revised: 10/31/2018] [Accepted: 11/20/2018] [Indexed: 12/11/2022] Open
Abstract
Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is the most popular assay to identify genomic regions, called ChIP-seq peaks, that are bound in vivo by transcription factors (TFs). These regions are derived from direct TF-DNA interactions, indirect binding of the TF to the DNA (through a co-binding partner), nonspecific binding to the DNA, and noise/bias/artifacts. Delineating the bona fide direct TF-DNA interactions within the ChIP-seq peaks remains challenging. We developed a dedicated software, ChIP-eat, that combines computational TF binding models and ChIP-seq peaks to automatically predict direct TF-DNA interactions. Our work culminated with predicted interactions covering >4% of the human genome, obtained by uniformly processing 1983 ChIP-seq peak data sets from the ReMap database for 232 unique TFs. The predictions were a posteriori assessed using protein binding microarray and ChIP-exo data, and were predominantly found in high quality ChIP-seq peaks. The set of predicted direct TF-DNA interactions suggested that high-occupancy target regions are likely not derived from direct binding of the TFs to the DNA. Our predictions derived co-binding TFs supported by protein-protein interaction data and defined cis-regulatory modules enriched for disease- and trait-associated SNPs. We provide this collection of direct TF-DNA interactions and cis-regulatory modules through the UniBind web-interface (http://unibind.uio.no).
Collapse
Affiliation(s)
- Marius Gheorghe
- Centre for Molecular Medicine Norway (NCMM), University of Oslo, Oslo, Norway
| | | | - Aziz Khan
- Centre for Molecular Medicine Norway (NCMM), University of Oslo, Oslo, Norway
| | - Jeanne Chèneby
- Aix Marseille Université, INSERM, TAGC, Marseille, France
| | | | - Anthony Mathelier
- Centre for Molecular Medicine Norway (NCMM), University of Oslo, Oslo, Norway.,Department of Cancer Genetics, Institute for Cancer Research, Radiumhospitalet, Oslo, Norway
| |
Collapse
|
14
|
Noguchi YT, Nakamura M, Hino N, Nogami J, Tsuji S, Sato T, Zhang L, Tsujikawa K, Tanaka T, Izawa K, Okada Y, Doi T, Kokubo H, Harada A, Uezumi A, Gessler M, Ohkawa Y, Fukada SI. Cell-autonomous and redundant roles of Hey1 and HeyL in muscle stem cells: HeyL requires Hes1 to bind diverse DNA sites. Development 2019; 146:dev.163618. [DOI: 10.1242/dev.163618] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2018] [Accepted: 01/31/2019] [Indexed: 12/20/2022]
Abstract
The undifferentiated state of muscle stem (satellite) cells (MuSCs) is maintained by the canonical Notch pathway. Although three bHLH transcriptional factors, Hey1, HeyL, and Hes1, are considered to be potential effectors of the Notch pathway exerting anti-myogenic effects, neither HeyL nor Hes1 inhibits myogenic differentiation of myogenic cell lines. Furthermore, whether these factors work redundantly or cooperatively is unknown. Here, we showed cell-autonomous functions of Hey1 and HeyL in MuSCs using conditional and genetic null mice. Analysis of cultured MuSCs revealed anti-myogenic activity of both HeyL and Hes1. We found that HeyL forms heterodimeric complexes with Hes1 in living cells. Moreover, our ChIP-Seq experiments demonstrated that, compared with HeyL alone, HeyL-Hes1 heterodimer bound with high affinity to specific sites in the chromatin including the cis-element of Hey1. Finally, the analyses of myogenin promoter activity showed that HeyL and Hes1 acted synergistically to suppress myogenic differentiation. Collectively, those results suggest that HeyL and Hey1 function redundantly in MuSCs, and that HeyL requires Hes1 for effective DNA binding and biological activity.
Collapse
Affiliation(s)
- Yu-taro Noguchi
- Laboratory of Molecular and Cellular Physiology, Graduate School of Pharmaceutical Sciences, Osaka University, 1-6 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Miki Nakamura
- Laboratory of Molecular and Cellular Physiology, Graduate School of Pharmaceutical Sciences, Osaka University, 1-6 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Nobumasa Hino
- Laboratory of Molecular Medicine, Graduate School of Pharmaceutical Sciences, Osaka University, 1-6 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Jumpei Nogami
- Division of Transcriptomics, Medical Institute of Bioregulation, Kyushu University, Fukuoka 812-8582, Japan
| | - Sayaka Tsuji
- Laboratory of Molecular Medicine, Graduate School of Pharmaceutical Sciences, Osaka University, 1-6 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Takahiko Sato
- Department of Ophthalmology, Kyoto Prefectural University of Medicine, Kyoto 602-8566, Japan
| | - Lidan Zhang
- Laboratory of Molecular and Cellular Physiology, Graduate School of Pharmaceutical Sciences, Osaka University, 1-6 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Kazutake Tsujikawa
- Laboratory of Molecular and Cellular Physiology, Graduate School of Pharmaceutical Sciences, Osaka University, 1-6 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Toru Tanaka
- Laboratory of Molecular Medicine, Graduate School of Pharmaceutical Sciences, Osaka University, 1-6 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Kohei Izawa
- Laboratory of Molecular Medicine, Graduate School of Pharmaceutical Sciences, Osaka University, 1-6 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Yoshiaki Okada
- Laboratory of Molecular Medicine, Graduate School of Pharmaceutical Sciences, Osaka University, 1-6 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Takefumi Doi
- Laboratory of Molecular Medicine, Graduate School of Pharmaceutical Sciences, Osaka University, 1-6 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Hiroki Kokubo
- Department of Cardiovascular Physiology and Medicine, Graduate School of Biomedical and Health Sciences, Hiroshima University, 1-2-3 Kasumi, Minamiku, Hiroshima 734-8551, Japan
| | - Akihito Harada
- Division of Transcriptomics, Medical Institute of Bioregulation, Kyushu University, Fukuoka 812-8582, Japan
| | - Akiyoshi Uezumi
- Department of Geriatric Medicine, Tokyo Metropolitan Institute of Gerontology, Itabashi-ku, Tokyo 173-0015, Japan
| | - Manfred Gessler
- Developmental Biochemistry, Theodor-Boveri-Institute / Biocenter, and Comprehensive Cancer Center Mainfranken, University of Wuerzburg, 97074 Wuerzburg, Germany
| | - Yasuyuki Ohkawa
- Division of Transcriptomics, Medical Institute of Bioregulation, Kyushu University, Fukuoka 812-8582, Japan
| | - So-ichiro Fukada
- Laboratory of Molecular and Cellular Physiology, Graduate School of Pharmaceutical Sciences, Osaka University, 1-6 Yamadaoka, Suita, Osaka 565-0871, Japan
| |
Collapse
|
15
|
Wu C, Jiao Y, Shen M, Pan C, Cheng G, Jia D, Zhu J, Zhang L, Zheng M, Jia J. Clustering-local-unique-enriched-signals (CLUES) promotes identification of novel regulators of ES cell self-renewal and pluripotency. PLoS One 2018; 13:e0206844. [PMID: 30399165 PMCID: PMC6219791 DOI: 10.1371/journal.pone.0206844] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2018] [Accepted: 10/19/2018] [Indexed: 01/08/2023] Open
Abstract
BACKGROUND Key regulators of developmental processes can be prioritized through integrated analysis of ChIP-Seq data of master transcriptional factors (TFs) such as Nanog and Oct4, active histone modifications (HMs) such as H3K4me3 and H3K27ac, and repressive HMs such as H3K27me3. Recent studies show that broad enrichment signals such as super-enhancers and broad H3K4me3 enrichment signals play more dominant roles than short enrichment signals of the master TFs and H3K4me3 in epigenetic regulatory mechanism. Besides the broad enrichment signals, up to ten thousands of short enrichment signals of these TFs and HMs exist in genome. Prioritization of these broad enrichment signals from ChIP-Seq data is a prerequisite for such integrated analysis. RESULTS Here, we present a method named Clustering-Local-Unique-Enriched-Signals (CLUES), which uses an adaptive-size-windows strategy to identify enriched regions (ERs) and cluster them into broad enrichment signals. Tested on 62 ENCODE ChIP-Seq datasets of Ctcf and Nrsf, CLUES performs equally well as MACS2 regarding prioritization of ERs with the TF's motif. Tested on 165 ENCODE ChIP-Seq datasets of H3K4me3, H3K27me3, and H3K36me3, CLUES performs better than existing algorithms on prioritizing broad enrichment signals implicating cell functions influenced by epigenetic regulatory mechanism in cells. Most importantly, CLUES helps to confirm several novel regulators of mouse ES cell self-renewal and pluripotency through integrated analysis of prioritized broad enrichment signals of H3K4me3, H3K27me3, Nanog and Oct4 with the support of a CRISPR/Cas9 negative selection genetic screen. CONCLUSIONS CLUES holds promise for prioritizing broad enrichment signals from ChIP-Seq data. The download site for CLUES is https://github.com/Wuchao1984/CLUESv1.
Collapse
Affiliation(s)
- Chao Wu
- Life Sciences Institute and Innovation Center for Cell Signaling Network, Zhejiang University, Hangzhou, Zhejiang, PRC
- * E-mail: (MZ); (CW); (JJ)
| | - Yang Jiao
- Life Sciences Institute and Innovation Center for Cell Signaling Network, Zhejiang University, Hangzhou, Zhejiang, PRC
| | - Manli Shen
- Life Sciences Institute and Innovation Center for Cell Signaling Network, Zhejiang University, Hangzhou, Zhejiang, PRC
| | - Chen Pan
- Life Sciences Institute and Innovation Center for Cell Signaling Network, Zhejiang University, Hangzhou, Zhejiang, PRC
| | - Guo Cheng
- Life Sciences Institute and Innovation Center for Cell Signaling Network, Zhejiang University, Hangzhou, Zhejiang, PRC
| | - Danmei Jia
- Life Sciences Institute and Innovation Center for Cell Signaling Network, Zhejiang University, Hangzhou, Zhejiang, PRC
| | - Jing Zhu
- Beijing Ming-tian Genetics Ltd., Beijing, PRC
| | - Long Zhang
- Life Sciences Institute and Innovation Center for Cell Signaling Network, Zhejiang University, Hangzhou, Zhejiang, PRC
| | - Min Zheng
- Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, Zhejiang University, Hangzhou, Zhejiang, PRC
- State Key Laboratory for Diagnosis and Treatment of Infectious Diseases, The First Affiliated Hospital, Zhejiang University, Hangzhou, Zhejiang, PRC
- * E-mail: (MZ); (CW); (JJ)
| | - Junling Jia
- Life Sciences Institute and Innovation Center for Cell Signaling Network, Zhejiang University, Hangzhou, Zhejiang, PRC
- Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, Zhejiang University, Hangzhou, Zhejiang, PRC
- State Key Laboratory for Diagnosis and Treatment of Infectious Diseases, The First Affiliated Hospital, Zhejiang University, Hangzhou, Zhejiang, PRC
- * E-mail: (MZ); (CW); (JJ)
| |
Collapse
|
16
|
Fiziev P, Ernst J. ChromTime: modeling spatio-temporal dynamics of chromatin marks. Genome Biol 2018; 19:109. [PMID: 30097049 PMCID: PMC6085762 DOI: 10.1186/s13059-018-1485-2] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2017] [Accepted: 07/17/2018] [Indexed: 11/17/2022] Open
Abstract
To model spatial changes of chromatin mark peaks over time we develop and apply ChromTime, a computational method that predicts peaks to be either expanding, contracting, or holding steady between time points. Predicted expanding and contracting peaks can mark regulatory regions associated with transcription factor binding and gene expression changes. Spatial dynamics of peaks provide information about gene expression changes beyond localized signal density changes. ChromTime detects asymmetric expansions and contractions, which for some marks associate with the direction of transcription. ChromTime facilitates the analysis of time course chromatin data in a range of biological systems.
Collapse
Affiliation(s)
- Petko Fiziev
- Bioinformatics Interdepartmental Program, University of California, Los Angeles, CA, USA.,Department of Biological Chemistry, University of California, Los Angeles, CA, USA.,Eli and Edythe Broad Center of Regenerative Medicine and Stem Cell Research at UCLA, Los Angeles, CA, USA
| | - Jason Ernst
- Bioinformatics Interdepartmental Program, University of California, Los Angeles, CA, USA. .,Department of Biological Chemistry, University of California, Los Angeles, CA, USA. .,Eli and Edythe Broad Center of Regenerative Medicine and Stem Cell Research at UCLA, Los Angeles, CA, USA. .,Computer Science Department, University of California, Los Angeles, CA, USA. .,Jonsson Comprehensive Cancer Center, University of California, Los Angeles, CA, USA. .,Molecular Biology Institute, University of California, Los Angeles, CA, USA.
| |
Collapse
|
17
|
Dai Z, Mentch SJ, Gao X, Nichenametla SN, Locasale JW. Methionine metabolism influences genomic architecture and gene expression through H3K4me3 peak width. Nat Commun 2018; 9:1955. [PMID: 29769529 PMCID: PMC5955993 DOI: 10.1038/s41467-018-04426-y] [Citation(s) in RCA: 83] [Impact Index Per Article: 13.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2017] [Accepted: 04/27/2018] [Indexed: 01/07/2023] Open
Abstract
Nutrition and metabolism are known to influence chromatin biology and epigenetics through post-translational modifications, yet how this interaction influences genomic architecture and connects to gene expression is unknown. Here we consider, as a model, the metabolically-driven dynamics of H3K4me3, a histone methylation mark that is known to encode information about active transcription, cell identity, and tumor suppression. We analyze the genome-wide changes in H3K4me3 and gene expression in response to alterations in methionine availability in both normal mouse physiology and human cancer cells. Surprisingly, we find that the location of H3K4me3 peaks is largely preserved under methionine restriction, while the response of H3K4me3 peak width encodes almost all aspects of H3K4me3 biology including changes in expression levels, and the presence of cell identity and cancer-associated genes. These findings may reveal general principles for how nutrient availability modulates specific aspects of chromatin dynamics to mediate biological function.
Collapse
Affiliation(s)
- Ziwei Dai
- Department of Pharmacology and Cancer Biology, Duke Molecular Physiology Institute, Duke Cancer Institute, Duke University School of Medicine, Durham, NC, 27710, USA
| | - Samantha J Mentch
- Department of Pharmacology and Cancer Biology, Duke Molecular Physiology Institute, Duke Cancer Institute, Duke University School of Medicine, Durham, NC, 27710, USA
| | - Xia Gao
- Department of Pharmacology and Cancer Biology, Duke Molecular Physiology Institute, Duke Cancer Institute, Duke University School of Medicine, Durham, NC, 27710, USA
| | | | - Jason W Locasale
- Department of Pharmacology and Cancer Biology, Duke Molecular Physiology Institute, Duke Cancer Institute, Duke University School of Medicine, Durham, NC, 27710, USA.
| |
Collapse
|
18
|
Wang X, Sterr M, Burtscher I, Chen S, Hieronimus A, Machicao F, Staiger H, Häring HU, Lederer G, Meitinger T, Cernilogar FM, Schotta G, Irmler M, Beckers J, Hrabě de Angelis M, Ray M, Wright CVE, Bakhti M, Lickert H. Genome-wide analysis of PDX1 target genes in human pancreatic progenitors. Mol Metab 2018; 9:57-68. [PMID: 29396371 PMCID: PMC5870105 DOI: 10.1016/j.molmet.2018.01.011] [Citation(s) in RCA: 56] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 11/16/2017] [Revised: 01/05/2018] [Accepted: 01/16/2018] [Indexed: 12/14/2022] Open
Abstract
Objective Homozygous loss-of-function mutations in the gene coding for the homeobox transcription factor (TF) PDX1 leads to pancreatic agenesis, whereas heterozygous mutations can cause Maturity-Onset Diabetes of the Young 4 (MODY4). Although the function of Pdx1 is well studied in pre-clinical models during insulin-producing β-cell development and homeostasis, it remains elusive how this TF controls human pancreas development by regulating a downstream transcriptional program. Also, comparative studies of PDX1 binding patterns in pancreatic progenitors and adult β-cells have not been conducted so far. Furthermore, many studies reported the association between single nucleotide polymorphisms (SNPs) and T2DM, and it has been shown that islet enhancers are enriched in T2DM-associated SNPs. Whether regions, harboring T2DM-associated SNPs are PDX1 bound and active at the pancreatic progenitor stage has not been reported so far. Methods In this study, we have generated a novel induced pluripotent stem cell (iPSC) line that efficiently differentiates into human pancreatic progenitors (PPs). Furthermore, PDX1 and H3K27ac chromatin immunoprecipitation sequencing (ChIP-seq) was used to identify PDX1 transcriptional targets and active enhancer and promoter regions. To address potential differences in the function of PDX1 during development and adulthood, we compared PDX1 binding profiles from PPs and adult islets. Moreover, combining ChIP-seq and GWAS meta-analysis data we identified T2DM-associated SNPs in PDX1 binding sites and active chromatin regions. Results ChIP-seq for PDX1 revealed a total of 8088 PDX1-bound regions that map to 5664 genes in iPSC-derived PPs. The PDX1 target regions include important pancreatic TFs, such as PDX1 itself, RFX6, HNF1B, and MEIS1, which were activated during the differentiation process as revealed by the active chromatin mark H3K27ac and mRNA expression profiling, suggesting that auto-regulatory feedback regulation maintains PDX1 expression and initiates a pancreatic TF program. Remarkably, we identified several PDX1 target genes that have not been reported in the literature in human so far, including RFX3, required for ciliogenesis and endocrine differentiation in mouse, and the ligand of the Notch receptor DLL1, which is important for endocrine induction and tip-trunk patterning. The comparison of PDX1 profiles from PPs and adult human islets identified sets of stage-specific target genes, associated with early pancreas development and adult β-cell function, respectively. Furthermore, we found an enrichment of T2DM-associated SNPs in active chromatin regions from iPSC-derived PPs. Two of these SNPs fall into PDX1 occupied sites that are located in the intronic regions of TCF7L2 and HNF1B. Both of these genes are key transcriptional regulators of endocrine induction and mutations in cis-regulatory regions predispose to diabetes. Conclusions Our data provide stage-specific target genes of PDX1 during in vitro differentiation of stem cells into pancreatic progenitors that could be useful to identify pathways and molecular targets that predispose for diabetes. In addition, we show that T2DM-associated SNPs are enriched in active chromatin regions at the pancreatic progenitor stage, suggesting that the susceptibility to T2DM might originate from imperfect execution of a β-cell developmental program. PDX1 ChIP-seq analysis revealed 5664 target genes in human pancreatic progenitors, including unreported target genes. Comparison of PDX1 profiles from PPs and adult human islets identified stage-specific PDX1 target gene sets. T2DM-associated SNPs are enriched in active chromatin regions from iPSC-derived PPs. Three SNPs fall into PDX1 occupied sites, located in intronic regions of the developmental regulatory TFs TCF7L2 and HNF1B.
Collapse
Affiliation(s)
- Xianming Wang
- Institute of Diabetes and Regeneration Research, Helmholtz Zentrum München, Parkring 11, 85748, Garching, Germany; Institute of Stem Cell Research, Helmholtz Zentrum München, 85764 Neuherberg, Germany; Chair of ß-Cell Biology, Technische Universität München, Ismaningerstraße 22, 81675 München, Germany
| | - Michael Sterr
- Institute of Diabetes and Regeneration Research, Helmholtz Zentrum München, Parkring 11, 85748, Garching, Germany; Institute of Stem Cell Research, Helmholtz Zentrum München, 85764 Neuherberg, Germany; Chair of ß-Cell Biology, Technische Universität München, Ismaningerstraße 22, 81675 München, Germany
| | - Ingo Burtscher
- Institute of Diabetes and Regeneration Research, Helmholtz Zentrum München, Parkring 11, 85748, Garching, Germany; Institute of Stem Cell Research, Helmholtz Zentrum München, 85764 Neuherberg, Germany
| | - Shen Chen
- iPS and Cancer Research Unit, Department of Histology and Embryology, Zhongshan School of Medicine, Sun Yat-sen University, Guangzhou 510080, China
| | - Anja Hieronimus
- Institute for Diabetes Research and Metabolic Diseases of the Helmholtz Zentrum München at the University of Tübingen, 72076 Tübingen, Germany; Department of Internal Medicine, Division of Endocrinology, Diabetology, Vascular Disease, Nephrology and Clinical Chemistry, University of Tübingen, 72076 Tübingen, Germany; German Center for Diabetes Research (DZD), 85764 Neuherberg, Germany
| | - Fausto Machicao
- Institute of Experimental Genetics, Helmholtz Zentrum München, 85764 Neuherberg, Germany; German Center for Diabetes Research (DZD), 85764 Neuherberg, Germany
| | - Harald Staiger
- Institute for Diabetes Research and Metabolic Diseases of the Helmholtz Zentrum München at the University of Tübingen, 72076 Tübingen, Germany; Institute of Pharmaceutical Sciences, Department of Pharmacy and Biochemistry, Eberhard Karls University Tübingen, 72076 Tübingen, Germany; German Center for Diabetes Research (DZD), 85764 Neuherberg, Germany
| | - Hans-Ulrich Häring
- Institute for Diabetes Research and Metabolic Diseases of the Helmholtz Zentrum München at the University of Tübingen, 72076 Tübingen, Germany; Department of Internal Medicine, Division of Endocrinology, Diabetology, Vascular Disease, Nephrology and Clinical Chemistry, University of Tübingen, 72076 Tübingen, Germany; German Center for Diabetes Research (DZD), 85764 Neuherberg, Germany
| | - Gabriele Lederer
- Institute of Human Genetics, Helmholtz Zentrum München, 85764 Neuherberg, Germany
| | - Thomas Meitinger
- Institute of Human Genetics, Helmholtz Zentrum München, 85764 Neuherberg, Germany
| | - Filippo M Cernilogar
- Biomedical Center and Center for Integrated Protein Science Munich, Ludwig-Maximilians-University, 82152 Planegg-Martinsried, Germany
| | - Gunnar Schotta
- Biomedical Center and Center for Integrated Protein Science Munich, Ludwig-Maximilians-University, 82152 Planegg-Martinsried, Germany
| | - Martin Irmler
- Institute of Experimental Genetics, Helmholtz Zentrum München, 85764 Neuherberg, Germany
| | - Johannes Beckers
- Institute of Experimental Genetics, Helmholtz Zentrum München, 85764 Neuherberg, Germany; German Center for Diabetes Research (DZD), 85764 Neuherberg, Germany; Chair of Experimental Genetics, School of Life Sciences Weihenstephan, Technische Universität München, 85354 Freising, Germany
| | - Martin Hrabě de Angelis
- Institute of Experimental Genetics, Helmholtz Zentrum München, 85764 Neuherberg, Germany; German Center for Diabetes Research (DZD), 85764 Neuherberg, Germany; Chair of Experimental Genetics, School of Life Sciences Weihenstephan, Technische Universität München, 85354 Freising, Germany
| | - Michael Ray
- Vanderbilt University Program in Developmental Biology, Department of Cell and Developmental Biology, Vanderbilt University, Nashville, TN 37232, USA; Vanderbilt Center for Stem Cell Biology, Vanderbilt University, Nashville, TN 37232, USA
| | - Christopher V E Wright
- Vanderbilt University Program in Developmental Biology, Department of Cell and Developmental Biology, Vanderbilt University, Nashville, TN 37232, USA; Vanderbilt Center for Stem Cell Biology, Vanderbilt University, Nashville, TN 37232, USA
| | - Mostafa Bakhti
- Institute of Diabetes and Regeneration Research, Helmholtz Zentrum München, Parkring 11, 85748, Garching, Germany; Institute of Stem Cell Research, Helmholtz Zentrum München, 85764 Neuherberg, Germany; German Center for Diabetes Research (DZD), 85764 Neuherberg, Germany
| | - Heiko Lickert
- Institute of Diabetes and Regeneration Research, Helmholtz Zentrum München, Parkring 11, 85748, Garching, Germany; Institute of Stem Cell Research, Helmholtz Zentrum München, 85764 Neuherberg, Germany; Chair of ß-Cell Biology, Technische Universität München, Ismaningerstraße 22, 81675 München, Germany; German Center for Diabetes Research (DZD), 85764 Neuherberg, Germany.
| |
Collapse
|
19
|
Nakato R, Shirahige K. Recent advances in ChIP-seq analysis: from quality management to whole-genome annotation. Brief Bioinform 2017; 18:279-290. [PMID: 26979602 PMCID: PMC5444249 DOI: 10.1093/bib/bbw023] [Citation(s) in RCA: 78] [Impact Index Per Article: 11.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2015] [Indexed: 02/06/2023] Open
Abstract
Chromatin immunoprecipitation followed by sequencing (ChIP-seq) analysis can detect protein/DNA-binding and histone-modification sites across an entire genome. Recent advances in sequencing technologies and analyses enable us to compare hundreds of samples simultaneously; such large-scale analysis has potential to reveal the high-dimensional interrelationship level for regulatory elements and annotate novel functional genomic regions de novo. Because many experimental considerations are relevant to the choice of a method in a ChIP-seq analysis, the overall design and quality management of the experiment are of critical importance. This review offers guiding principles of computation and sample preparation for ChIP-seq analyses, highlighting the validity and limitations of the state-of-the-art procedures at each step. We also discuss the latest challenges of single-cell analysis that will encourage a new era in this field.
Collapse
Affiliation(s)
- Ryuichiro Nakato
- Research Center for Epigenetic Disease, Institute of Molecular and Cellular Biosciences, The University of Tokyo, Tokyo, Japan
| | - Katsuhiko Shirahige
- Research Center for Epigenetic Disease, Institute of Molecular and Cellular Biosciences, The University of Tokyo, Tokyo, Japan.,Core Research for Evolutional Science and Technology (CREST), Japan Science and Technology Agency, Kawaguchi, Japan
| |
Collapse
|
20
|
Ferrai C, Torlai Triglia E, Risner-Janiczek JR, Rito T, Rackham OJ, de Santiago I, Kukalev A, Nicodemi M, Akalin A, Li M, Ungless MA, Pombo A. RNA polymerase II primes Polycomb-repressed developmental genes throughout terminal neuronal differentiation. Mol Syst Biol 2017; 13:946. [PMID: 29038337 PMCID: PMC5658700 DOI: 10.15252/msb.20177754] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
Polycomb repression in mouse embryonic stem cells (ESCs) is tightly associated with promoter co‐occupancy of RNA polymerase II (RNAPII) which is thought to prime genes for activation during early development. However, it is unknown whether RNAPII poising is a general feature of Polycomb repression, or is lost during differentiation. Here, we map the genome‐wide occupancy of RNAPII and Polycomb from pluripotent ESCs to non‐dividing functional dopaminergic neurons. We find that poised RNAPII complexes are ubiquitously present at Polycomb‐repressed genes at all stages of neuronal differentiation. We observe both loss and acquisition of RNAPII and Polycomb at specific groups of genes reflecting their silencing or activation. Strikingly, RNAPII remains poised at transcription factor genes which are silenced in neurons through Polycomb repression, and have major roles in specifying other, non‐neuronal lineages. We conclude that RNAPII poising is intrinsically associated with Polycomb repression throughout differentiation. Our work suggests that the tight interplay between RNAPII poising and Polycomb repression not only instructs promoter state transitions, but also may enable promoter plasticity in differentiated cells.
Collapse
Affiliation(s)
- Carmelo Ferrai
- Epigenetic Regulation and Chromatin Architecture, Max Delbrück Center for Molecular Medicine, Berlin, Germany .,Genome Function, MRC London Institute of Medical Sciences (previously MRC Clinical Sciences Centre), London, UK.,Institute of Clinical Sciences (ICS), Faculty of Medicine, Imperial College London, London, UK
| | - Elena Torlai Triglia
- Epigenetic Regulation and Chromatin Architecture, Max Delbrück Center for Molecular Medicine, Berlin, Germany
| | - Jessica R Risner-Janiczek
- Institute of Clinical Sciences (ICS), Faculty of Medicine, Imperial College London, London, UK.,Stem Cell Neurogenesis, MRC London Institute of Medical Sciences (previously MRC Clinical Sciences Centre), London, UK.,Neurophysiology Group, MRC London Institute of Medical Sciences (previously MRC Clinical Sciences Centre), London, UK
| | - Tiago Rito
- Epigenetic Regulation and Chromatin Architecture, Max Delbrück Center for Molecular Medicine, Berlin, Germany
| | | | - Inês de Santiago
- Genome Function, MRC London Institute of Medical Sciences (previously MRC Clinical Sciences Centre), London, UK.,Institute of Clinical Sciences (ICS), Faculty of Medicine, Imperial College London, London, UK
| | - Alexander Kukalev
- Epigenetic Regulation and Chromatin Architecture, Max Delbrück Center for Molecular Medicine, Berlin, Germany
| | - Mario Nicodemi
- Dipartimento di Fisica, Università di Napoli Federico II and INFN Napoli, Complesso Universitario di Monte Sant'Angelo, Naples, Italy
| | - Altuna Akalin
- Scientific Bioinformatics Platform, Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine, Berlin, Germany
| | - Meng Li
- Institute of Clinical Sciences (ICS), Faculty of Medicine, Imperial College London, London, UK.,Stem Cell Neurogenesis, MRC London Institute of Medical Sciences (previously MRC Clinical Sciences Centre), London, UK
| | - Mark A Ungless
- Institute of Clinical Sciences (ICS), Faculty of Medicine, Imperial College London, London, UK .,Neurophysiology Group, MRC London Institute of Medical Sciences (previously MRC Clinical Sciences Centre), London, UK
| | - Ana Pombo
- Epigenetic Regulation and Chromatin Architecture, Max Delbrück Center for Molecular Medicine, Berlin, Germany .,Genome Function, MRC London Institute of Medical Sciences (previously MRC Clinical Sciences Centre), London, UK.,Institute of Clinical Sciences (ICS), Faculty of Medicine, Imperial College London, London, UK.,Institute for Biology, Humboldt-Universität zu Berlin, Berlin Germany
| |
Collapse
|
21
|
Conrad RJ, Fozouni P, Thomas S, Sy H, Zhang Q, Zhou MM, Ott M. The Short Isoform of BRD4 Promotes HIV-1 Latency by Engaging Repressive SWI/SNF Chromatin-Remodeling Complexes. Mol Cell 2017; 67:1001-1012.e6. [PMID: 28844864 DOI: 10.1016/j.molcel.2017.07.025] [Citation(s) in RCA: 82] [Impact Index Per Article: 11.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2016] [Revised: 05/24/2017] [Accepted: 07/25/2017] [Indexed: 11/18/2022]
Abstract
BET proteins commonly activate cellular gene expression, yet inhibiting their recruitment paradoxically reactivates latent HIV-1 transcription. Here we identify the short isoform of BET family member BRD4 (BRD4S) as a corepressor of HIV-1 transcription. We found that BRD4S was enriched in chromatin fractions of latently infected T cells, and it was more rapidly displaced from chromatin upon BET inhibition than the long isoform. BET inhibition induced marked nucleosome remodeling at the latent HIV-1 promoter, which was dependent on the activity of BRG1-associated factors (BAF), an SWI/SNF chromatin-remodeling complex with known repressive functions in HIV-1 transcription. BRD4S directly bound BRG1, a catalytic subunit of BAF, via its bromodomain and extraterminal (ET) domain, and this isoform was necessary for BRG1 recruitment to latent HIV-1 chromatin. Using chromatin immunoprecipitation sequencing (ChIP-seq) combined with assay for transposase-accessible chromatin coupled to high-throughput sequencing (ATAC-seq) data, we found that the latent HIV-1 promoter phenotypically resembles endogenous long terminal repeat (LTR) sequences, pointing to a select role of BRD4S-BRG1 complexes in genomic silencing of invasive retroelements.
Collapse
Affiliation(s)
- Ryan J Conrad
- Gladstone Institutes, University of California, San Francisco, San Francisco, CA 94158, USA; Graduate Program in Pharmaceutical Sciences and Pharmacogenomics, University of California, San Francisco, San Francisco, CA 94158, USA; Department of Medicine, University of California, San Francisco, San Francisco, CA 94158, USA
| | - Parinaz Fozouni
- Gladstone Institutes, University of California, San Francisco, San Francisco, CA 94158, USA; Department of Medicine, University of California, San Francisco, San Francisco, CA 94158, USA
| | - Sean Thomas
- Gladstone Institutes, University of California, San Francisco, San Francisco, CA 94158, USA
| | - Hendrik Sy
- Gladstone Institutes, University of California, San Francisco, San Francisco, CA 94158, USA
| | - Qiang Zhang
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA
| | - Ming-Ming Zhou
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA
| | - Melanie Ott
- Gladstone Institutes, University of California, San Francisco, San Francisco, CA 94158, USA; Graduate Program in Pharmaceutical Sciences and Pharmacogenomics, University of California, San Francisco, San Francisco, CA 94158, USA; Department of Medicine, University of California, San Francisco, San Francisco, CA 94158, USA.
| |
Collapse
|
22
|
Ranciati S, Viroli C, Wit EC. Mixture model with multiple allocations for clustering spatially correlated observations in the analysis of ChIP-Seq data. Biom J 2017; 59:1301-1316. [PMID: 28664629 DOI: 10.1002/bimj.201600131] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2016] [Revised: 02/13/2017] [Accepted: 05/01/2017] [Indexed: 11/10/2022]
Abstract
Model-based clustering is a technique widely used to group a collection of units into mutually exclusive groups. There are, however, situations in which an observation could in principle belong to more than one cluster. In the context of next-generation sequencing (NGS) experiments, for example, the signal observed in the data might be produced by two (or more) different biological processes operating together and a gene could participate in both (or all) of them. We propose a novel approach to cluster NGS discrete data, coming from a ChIP-Seq experiment, with a mixture model, allowing each unit to belong potentially to more than one group: these multiple allocation clusters can be flexibly defined via a function combining the features of the original groups without introducing new parameters. The formulation naturally gives rise to a 'zero-inflation group' in which values close to zero can be allocated, acting as a correction for the abundance of zeros that manifest in this type of data. We take into account the spatial dependency between observations, which is described through a latent conditional autoregressive process that can reflect different dependency patterns. We assess the performance of our model within a simulation environment and then we apply it to ChIP-seq real data.
Collapse
Affiliation(s)
- Saverio Ranciati
- Department of Statistical Sciences "Paolo Fortunati", University of Bologna, via delle Belle Arti 41, 40126 Bologna, Italy
| | - Cinzia Viroli
- Department of Statistical Sciences "Paolo Fortunati", University of Bologna, via delle Belle Arti 41, 40126 Bologna, Italy
| | - Ernst C Wit
- Johann Bernoulli Institute for Mathematics and Computer Science, University of Groningen, Nijenborgh 9, 9747 AG Groningen, The Netherlands
| |
Collapse
|
23
|
Kar G, Kim JK, Kolodziejczyk AA, Natarajan KN, Torlai Triglia E, Mifsud B, Elderkin S, Marioni JC, Pombo A, Teichmann SA. Flipping between Polycomb repressed and active transcriptional states introduces noise in gene expression. Nat Commun 2017; 8:36. [PMID: 28652613 PMCID: PMC5484669 DOI: 10.1038/s41467-017-00052-2] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2016] [Accepted: 04/28/2017] [Indexed: 11/09/2022] Open
Abstract
Polycomb repressive complexes (PRCs) are important histone modifiers, which silence gene expression; yet, there exists a subset of PRC-bound genes actively transcribed by RNA polymerase II (RNAPII). It is likely that the role of Polycomb repressive complex is to dampen expression of these PRC-active genes. However, it is unclear how this flipping between chromatin states alters the kinetics of transcription. Here, we integrate histone modifications and RNAPII states derived from bulk ChIP-seq data with single-cell RNA-sequencing data. We find that Polycomb repressive complex-active genes have greater cell-to-cell variation in expression than active genes, and these results are validated by knockout experiments. We also show that PRC-active genes are clustered on chromosomes in both two and three dimensions, and interactions with active enhancers promote a stabilization of gene expression noise. These findings provide new insights into how chromatin regulation modulates stochastic gene expression and transcriptional bursting, with implications for regulation of pluripotency and development.Polycomb repressive complexes modify histones but it is unclear how changes in chromatin states alter kinetics of transcription. Here, the authors use single-cell RNAseq and ChIPseq to find that actively transcribed genes with Polycomb marks have greater cell-to-cell variation in expression.
Collapse
Affiliation(s)
- Gozde Kar
- European Molecular Biology Laboratory-European Bioinformatics Institute (EMBL-EBI), Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SD, UK
| | - Jong Kyoung Kim
- European Molecular Biology Laboratory-European Bioinformatics Institute (EMBL-EBI), Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SD, UK
- Department of New Biology, DGIST, Daegu, 42988, Republic of Korea
| | - Aleksandra A Kolodziejczyk
- European Molecular Biology Laboratory-European Bioinformatics Institute (EMBL-EBI), Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SD, UK
- Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA, UK
| | - Kedar Nath Natarajan
- European Molecular Biology Laboratory-European Bioinformatics Institute (EMBL-EBI), Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SD, UK
- Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA, UK
| | - Elena Torlai Triglia
- Epigenetic Regulation and Chromatin Architecture Group, Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine, Robert Roessle Strasse, Berlin-Buch, 13125, Germany
| | - Borbala Mifsud
- Cancer Research UK London Research Institute, 44 Lincoln's Inn Fields, London, WC2A 3LY, UK
- Department of Genetics, Evolution and Environment, University College London, Gower Street, London, WC1E 6BT, UK
- William Harvey Research Institute, Queen Mary University London, Charterhouse Square, London, EC1M 6BQ, UK
| | - Sarah Elderkin
- Nuclear Dynamics Programme, The Babraham Institute, Babraham Research Campus, Cambridge, CB22 3AT, UK
| | - John C Marioni
- European Molecular Biology Laboratory-European Bioinformatics Institute (EMBL-EBI), Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SD, UK
- Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA, UK
- Cancer Research UK Cambridge Institute, University of Cambridge, Li Ka Shing Centre, Robinson Way, Cambridge, CB2 0RE, UK
| | - Ana Pombo
- Epigenetic Regulation and Chromatin Architecture Group, Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine, Robert Roessle Strasse, Berlin-Buch, 13125, Germany
| | - Sarah A Teichmann
- European Molecular Biology Laboratory-European Bioinformatics Institute (EMBL-EBI), Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SD, UK.
- Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA, UK.
| |
Collapse
|
24
|
Nora EP, Goloborodko A, Valton AL, Gibcus JH, Uebersohn A, Abdennur N, Dekker J, Mirny LA, Bruneau BG. Targeted Degradation of CTCF Decouples Local Insulation of Chromosome Domains from Genomic Compartmentalization. Cell 2017; 169:930-944.e22. [PMID: 28525758 PMCID: PMC5538188 DOI: 10.1016/j.cell.2017.05.004] [Citation(s) in RCA: 1040] [Impact Index Per Article: 148.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2016] [Revised: 03/14/2017] [Accepted: 05/02/2017] [Indexed: 01/17/2023]
Abstract
The molecular mechanisms underlying folding of mammalian chromosomes remain poorly understood. The transcription factor CTCF is a candidate regulator of chromosomal structure. Using the auxin-inducible degron system in mouse embryonic stem cells, we show that CTCF is absolutely and dose-dependently required for looping between CTCF target sites and insulation of topologically associating domains (TADs). Restoring CTCF reinstates proper architecture on altered chromosomes, indicating a powerful instructive function for CTCF in chromatin folding. CTCF remains essential for TAD organization in non-dividing cells. Surprisingly, active and inactive genome compartments remain properly segregated upon CTCF depletion, revealing that compartmentalization of mammalian chromosomes emerges independently of proper insulation of TADs. Furthermore, our data support that CTCF mediates transcriptional insulator function through enhancer blocking but not as a direct barrier to heterochromatin spreading. Beyond defining the functions of CTCF in chromosome folding, these results provide new fundamental insights into the rules governing mammalian genome organization.
Collapse
Affiliation(s)
- Elphège P Nora
- Gladstone Institute of Cardiovascular Disease, San Francisco, CA 94158, USA; Roddenberry Center for Stem Cell Biology and Medicine at Gladstone, San Francisco, CA 94158, USA.
| | - Anton Goloborodko
- Institute for Medical Engineering and Science and Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Anne-Laure Valton
- Howard Hughes Medical Institute, Program in Systems Biology, Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, MA 01605-0103, USA
| | - Johan H Gibcus
- Howard Hughes Medical Institute, Program in Systems Biology, Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, MA 01605-0103, USA
| | - Alec Uebersohn
- Gladstone Institute of Cardiovascular Disease, San Francisco, CA 94158, USA; Roddenberry Center for Stem Cell Biology and Medicine at Gladstone, San Francisco, CA 94158, USA
| | - Nezar Abdennur
- Institute for Medical Engineering and Science and Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Job Dekker
- Howard Hughes Medical Institute, Program in Systems Biology, Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, MA 01605-0103, USA
| | - Leonid A Mirny
- Institute for Medical Engineering and Science and Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Benoit G Bruneau
- Gladstone Institute of Cardiovascular Disease, San Francisco, CA 94158, USA; Roddenberry Center for Stem Cell Biology and Medicine at Gladstone, San Francisco, CA 94158, USA; Department of Pediatrics, University of California, San Francisco, CA 94143, USA; Cardiovascular Research Institute, University of California, San Francisco, CA 94158, USA.
| |
Collapse
|
25
|
Thomas R, Thomas S, Holloway AK, Pollard KS. Features that define the best ChIP-seq peak calling algorithms. Brief Bioinform 2017; 18:441-450. [PMID: 27169896 PMCID: PMC5429005 DOI: 10.1093/bib/bbw035] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2016] [Revised: 03/01/2016] [Indexed: 12/20/2022] Open
Abstract
Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is an important tool for studying gene regulatory proteins, such as transcription factors and histones. Peak calling is one of the first steps in the analysis of these data. Peak calling consists of two sub-problems: identifying candidate peaks and testing candidate peaks for statistical significance. We surveyed 30 methods and identified 12 features of the two sub-problems that distinguish methods from each other. We picked six methods GEM, MACS2, MUSIC, BCP, Threshold-based method (TM) and ZINBA] that span this feature space and used a combination of 300 simulated ChIP-seq data sets, 3 real data sets and mathematical analyses to identify features of methods that allow some to perform better than the others. We prove that methods that explicitly combine the signals from ChIP and input samples are less powerful than methods that do not. Methods that use windows of different sizes are more powerful than the ones that do not. For statistical testing of candidate peaks, methods that use a Poisson test to rank their candidate peaks are more powerful than those that use a Binomial test. BCP and MACS2 have the best operating characteristics on simulated transcription factor binding data. GEM has the highest fraction of the top 500 peaks containing the binding motif of the immunoprecipitated factor, with 50% of its peaks within 10 base pairs of a motif. BCP and MUSIC perform best on histone data. These findings provide guidance and rationale for selecting the best peak caller for a given application.
Collapse
Affiliation(s)
| | - Sean Thomas
- Gladstone Institutes, San Francisco, CA, USA
- Division of Biostatistics, University of California, San Francisco, CA, USA
| | - Alisha K Holloway
- Gladstone Institutes, San Francisco, CA, USA
- Division of Biostatistics, University of California, San Francisco, CA, USA
- Phylos Biosciences, Portland, OR, USA
| | - Katherine S Pollard
- Gladstone Institutes, San Francisco, CA, USA
- Division of Biostatistics, University of California, San Francisco, CA, USA
- Institute for Human Genetics and Institute for Computational Health Sciences, University of California, San Francisco, CA, USA
| |
Collapse
|
26
|
Taguchi H, Xie Y, Horikoshi N, Maehara K, Harada A, Nogami J, Sato K, Arimura Y, Osakabe A, Kujirai T, Iwasaki T, Semba Y, Tachibana T, Kimura H, Ohkawa Y, Kurumizaka H. Crystal Structure and Characterization of Novel Human Histone H3 Variants, H3.6, H3.7, and H3.8. Biochemistry 2017; 56:2184-2196. [DOI: 10.1021/acs.biochem.6b01098] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Hiroyuki Taguchi
- Laboratory of Structural
Biology,
Graduate School of Advanced Science and Engineering, Research Institute
for Science and Engineering, and Institute for Medical-oriented Structural
Biology, Waseda University, 2-2 Wakamatsu-cho, Shinjuku-ku, Tokyo 162-8480, Japan
| | - Yan Xie
- Laboratory of Structural
Biology,
Graduate School of Advanced Science and Engineering, Research Institute
for Science and Engineering, and Institute for Medical-oriented Structural
Biology, Waseda University, 2-2 Wakamatsu-cho, Shinjuku-ku, Tokyo 162-8480, Japan
| | - Naoki Horikoshi
- Laboratory of Structural
Biology,
Graduate School of Advanced Science and Engineering, Research Institute
for Science and Engineering, and Institute for Medical-oriented Structural
Biology, Waseda University, 2-2 Wakamatsu-cho, Shinjuku-ku, Tokyo 162-8480, Japan
| | - Kazumitsu Maehara
- Division of Transcriptomics, Medical
Institute of Bioregulation, Kyushu University, 3-1-1 Maidashi,
Higashi-ku, Fukuoka 812-8582, Japan
| | - Akihito Harada
- Division of Transcriptomics, Medical
Institute of Bioregulation, Kyushu University, 3-1-1 Maidashi,
Higashi-ku, Fukuoka 812-8582, Japan
| | - Jumpei Nogami
- Division of Transcriptomics, Medical
Institute of Bioregulation, Kyushu University, 3-1-1 Maidashi,
Higashi-ku, Fukuoka 812-8582, Japan
| | - Koichi Sato
- Laboratory of Structural
Biology,
Graduate School of Advanced Science and Engineering, Research Institute
for Science and Engineering, and Institute for Medical-oriented Structural
Biology, Waseda University, 2-2 Wakamatsu-cho, Shinjuku-ku, Tokyo 162-8480, Japan
| | - Yasuhiro Arimura
- Laboratory of Structural
Biology,
Graduate School of Advanced Science and Engineering, Research Institute
for Science and Engineering, and Institute for Medical-oriented Structural
Biology, Waseda University, 2-2 Wakamatsu-cho, Shinjuku-ku, Tokyo 162-8480, Japan
| | - Akihisa Osakabe
- Laboratory of Structural
Biology,
Graduate School of Advanced Science and Engineering, Research Institute
for Science and Engineering, and Institute for Medical-oriented Structural
Biology, Waseda University, 2-2 Wakamatsu-cho, Shinjuku-ku, Tokyo 162-8480, Japan
| | - Tomoya Kujirai
- Laboratory of Structural
Biology,
Graduate School of Advanced Science and Engineering, Research Institute
for Science and Engineering, and Institute for Medical-oriented Structural
Biology, Waseda University, 2-2 Wakamatsu-cho, Shinjuku-ku, Tokyo 162-8480, Japan
| | - Takeshi Iwasaki
- Division of Transcriptomics, Medical
Institute of Bioregulation, Kyushu University, 3-1-1 Maidashi,
Higashi-ku, Fukuoka 812-8582, Japan
| | - Yuichiro Semba
- Division of Transcriptomics, Medical
Institute of Bioregulation, Kyushu University, 3-1-1 Maidashi,
Higashi-ku, Fukuoka 812-8582, Japan
| | - Taro Tachibana
- Department of Bioengineering, Graduate
School of Engineering, Osaka City University, Sugimoto, Sumiyoshi-ku, Osaka 558-8585, Japan
| | - Hiroshi Kimura
- Cell Biology Unit,
Institute of
Innovative Research, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Yokohama 226-8501, Japan
| | - Yasuyuki Ohkawa
- Division of Transcriptomics, Medical
Institute of Bioregulation, Kyushu University, 3-1-1 Maidashi,
Higashi-ku, Fukuoka 812-8582, Japan
| | - Hitoshi Kurumizaka
- Laboratory of Structural
Biology,
Graduate School of Advanced Science and Engineering, Research Institute
for Science and Engineering, and Institute for Medical-oriented Structural
Biology, Waseda University, 2-2 Wakamatsu-cho, Shinjuku-ku, Tokyo 162-8480, Japan
| |
Collapse
|
27
|
Cormier N, Kolisnik T, Bieda M. Reusable, extensible, and modifiable R scripts and Kepler workflows for comprehensive single set ChIP-seq analysis. BMC Bioinformatics 2016; 17:270. [PMID: 27377783 PMCID: PMC4932705 DOI: 10.1186/s12859-016-1125-3] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2016] [Accepted: 06/07/2016] [Indexed: 11/10/2022] Open
Abstract
Background There has been an enormous expansion of use of chromatin immunoprecipitation followed by sequencing (ChIP-seq) technologies. Analysis of large-scale ChIP-seq datasets involves a complex series of steps and production of several specialized graphical outputs. A number of systems have emphasized custom development of ChIP-seq pipelines. These systems are primarily based on custom programming of a single, complex pipeline or supply libraries of modules and do not produce the full range of outputs commonly produced for ChIP-seq datasets. It is desirable to have more comprehensive pipelines, in particular ones addressing common metadata tasks, such as pathway analysis, and pipelines producing standard complex graphical outputs. It is advantageous if these are highly modular systems, available as both turnkey pipelines and individual modules, that are easily comprehensible, modifiable and extensible to allow rapid alteration in response to new analysis developments in this growing area. Furthermore, it is advantageous if these pipelines allow data provenance tracking. Results We present a set of 20 ChIP-seq analysis software modules implemented in the Kepler workflow system; most (18/20) were also implemented as standalone, fully functional R scripts. The set consists of four full turnkey pipelines and 16 component modules. The turnkey pipelines in Kepler allow data provenance tracking. Implementation emphasized use of common R packages and widely-used external tools (e.g., MACS for peak finding), along with custom programming. This software presents comprehensive solutions and easily repurposed code blocks for ChIP-seq analysis and pipeline creation. Tasks include mapping raw reads, peakfinding via MACS, summary statistics, peak location statistics, summary plots centered on the transcription start site (TSS), gene ontology, pathway analysis, and de novo motif finding, among others. Conclusions These pipelines range from those performing a single task to those performing full analyses of ChIP-seq data. The pipelines are supplied as both Kepler workflows, which allow data provenance tracking, and, in the majority of cases, as standalone R scripts. These pipelines are designed for ease of modification and repurposing.
Collapse
Affiliation(s)
- Nathan Cormier
- Department of Biochemistry and Molecular Biology, University of Calgary Cumming School of Medicine, Rm HSC1151, 3330 Hospital Dr. NW, Calgary, AB, T2N4N1, Canada
| | - Tyler Kolisnik
- Department of Biochemistry and Molecular Biology, University of Calgary Cumming School of Medicine, Rm HSC1151, 3330 Hospital Dr. NW, Calgary, AB, T2N4N1, Canada
| | - Mark Bieda
- Department of Biochemistry and Molecular Biology, University of Calgary Cumming School of Medicine, Rm HSC1151, 3330 Hospital Dr. NW, Calgary, AB, T2N4N1, Canada.
| |
Collapse
|
28
|
Chen X, Jung JG, Shajahan-Haq AN, Clarke R, Shih IM, Wang Y, Magnani L, Wang TL, Xuan J. ChIP-BIT: Bayesian inference of target genes using a novel joint probabilistic model of ChIP-seq profiles. Nucleic Acids Res 2016; 44:e65. [PMID: 26704972 PMCID: PMC4838354 DOI: 10.1093/nar/gkv1491] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2015] [Revised: 11/16/2015] [Accepted: 12/09/2015] [Indexed: 11/16/2022] Open
Abstract
Chromatin immunoprecipitation with massively parallel DNA sequencing (ChIP-seq) has greatly improved the reliability with which transcription factor binding sites (TFBSs) can be identified from genome-wide profiling studies. Many computational tools are developed to detect binding events or peaks, however the robust detection of weak binding events remains a challenge for current peak calling tools. We have developed a novel Bayesian approach (ChIP-BIT) to reliably detect TFBSs and their target genes by jointly modeling binding signal intensities and binding locations of TFBSs. Specifically, a Gaussian mixture model is used to capture both binding and background signals in sample data. As a unique feature of ChIP-BIT, background signals are modeled by a local Gaussian distribution that is accurately estimated from the input data. Extensive simulation studies showed a significantly improved performance of ChIP-BIT in target gene prediction, particularly for detecting weak binding signals at gene promoter regions. We applied ChIP-BIT to find target genes from NOTCH3 and PBX1 ChIP-seq data acquired from MCF-7 breast cancer cells. TF knockdown experiments have initially validated about 30% of co-regulated target genes identified by ChIP-BIT as being differentially expressed in MCF-7 cells. Functional analysis on these genes further revealed the existence of crosstalk between Notch and Wnt signaling pathways.
Collapse
Affiliation(s)
- Xi Chen
- Bradley Department of Electrical and Computer Engineering, Virginia Polytechnic Institute and State University, 900 North Glebe Road, Arlington, VA 22203, USA
| | - Jin-Gyoung Jung
- Department of Pathology, Johns Hopkins Medical Institutions, 1550 Orleans Street, CRB-II, Baltimore, MD 21231, USA
| | - Ayesha N Shajahan-Haq
- Department of Oncology, Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, 3970 Reservoir Road NW, Washington, DC 20057, USA
| | - Robert Clarke
- Department of Oncology, Lombardi Comprehensive Cancer Center, Georgetown University Medical Center, 3970 Reservoir Road NW, Washington, DC 20057, USA
| | - Ie-Ming Shih
- Department of Pathology, Johns Hopkins Medical Institutions, 1550 Orleans Street, CRB-II, Baltimore, MD 21231, USA
| | - Yue Wang
- Bradley Department of Electrical and Computer Engineering, Virginia Polytechnic Institute and State University, 900 North Glebe Road, Arlington, VA 22203, USA
| | - Luca Magnani
- Department of Surgery and Cancer, Imperial College London, ICTEM building, Hammersmith Hospital, DuCane Road, London W120NN, UK
| | - Tian-Li Wang
- Department of Pathology, Johns Hopkins Medical Institutions, 1550 Orleans Street, CRB-II, Baltimore, MD 21231, USA
| | - Jianhua Xuan
- Bradley Department of Electrical and Computer Engineering, Virginia Polytechnic Institute and State University, 900 North Glebe Road, Arlington, VA 22203, USA
| |
Collapse
|
29
|
Starmer J, Magnuson T. Detecting broad domains and narrow peaks in ChIP-seq data with hiddenDomains. BMC Bioinformatics 2016; 17:144. [PMID: 27009150 PMCID: PMC4806451 DOI: 10.1186/s12859-016-0991-z] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2015] [Accepted: 03/15/2016] [Indexed: 01/22/2023] Open
Abstract
Background Correctly identifying genomic regions enriched with histone modifications and transcription factors is key to understanding their regulatory and developmental roles. Conceptually, these regions are divided into two categories, narrow peaks and broad domains, and different algorithms are used to identify each one. Datasets that span these two categories are often analyzed with a single program for peak calling combined with an ad hoc method for domains. Results We developed hiddenDomains, which identifies both peaks and domains, and compare it to the leading algorithms using H3K27me3, H3K36me3, GABP, ESR1 and FOXA ChIP-seq datasets. The output from the programs was compared to qPCR-validated enriched and depleted sites, predicted transcription factor binding sites, and highly-transcribed gene bodies. With every method, hiddenDomains, performed as well as, if not better than algorithms dedicated to a specific type of analysis. Conclusions hiddenDomains performs as well as the best domain and peak calling algorithms, making it ideal for analyzing ChIP-seq datasets, especially those that contain a mixture of peaks and domains. Electronic supplementary material The online version of this article (doi:10.1186/s12859-016-0991-z) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Joshua Starmer
- Department of Genetics, The University of North Carolina at Chapel Hill, Chapel Hill, NC, 27599, USA. .,Lineberger Comprehensive Cancer Center, The University of North Carolina at Chapel Hill, Chapel Hill, NC, 27599, USA.
| | - Terry Magnuson
- Department of Genetics, The University of North Carolina at Chapel Hill, Chapel Hill, NC, 27599, USA.,Lineberger Comprehensive Cancer Center, The University of North Carolina at Chapel Hill, Chapel Hill, NC, 27599, USA
| |
Collapse
|
30
|
Dias JD, Rito T, Torlai Triglia E, Kukalev A, Ferrai C, Chotalia M, Brookes E, Kimura H, Pombo A. Methylation of RNA polymerase II non-consensus Lysine residues marks early transcription in mammalian cells. eLife 2015; 4. [PMID: 26687004 PMCID: PMC4758952 DOI: 10.7554/elife.11215] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2015] [Accepted: 12/18/2015] [Indexed: 12/16/2022] Open
Abstract
Dynamic post-translational modification of RNA polymerase II (RNAPII) coordinates the co-transcriptional recruitment of enzymatic complexes that regulate chromatin states and processing of nascent RNA. Extensive phosphorylation of serine residues at the largest RNAPII subunit occurs at its structurally-disordered C-terminal domain (CTD), which is composed of multiple heptapeptide repeats with consensus sequence Y1-S2-P3-T4-S5-P6-S7. Serine-5 and Serine-7 phosphorylation mark transcription initiation, whereas Serine-2 phosphorylation coincides with productive elongation. In vertebrates, the CTD has eight non-canonical substitutions of Serine-7 into Lysine-7, which can be acetylated (K7ac). Here, we describe mono- and di-methylation of CTD Lysine-7 residues (K7me1 and K7me2). K7me1 and K7me2 are observed during the earliest transcription stages and precede or accompany Serine-5 and Serine-7 phosphorylation. In contrast, K7ac is associated with RNAPII elongation, Serine-2 phosphorylation and mRNA expression. We identify an unexpected balance between RNAPII K7 methylation and acetylation at gene promoters, which fine-tunes gene expression levels.
Collapse
Affiliation(s)
- João D Dias
- Epigenetic Regulation and Chromatin Architecture Group, Berlin Institute for Medical Systems Biology, Max-Delbrück Centre for Molecular Medicine, Berlin, Germany.,Genome Function Group, MRC Clinical Sciences Centre, Imperial College London, London, United Kingdom.,Graduate Program in Areas of Basic and Applied Biology, University of Porto, Porto, Portugal
| | - Tiago Rito
- Epigenetic Regulation and Chromatin Architecture Group, Berlin Institute for Medical Systems Biology, Max-Delbrück Centre for Molecular Medicine, Berlin, Germany
| | - Elena Torlai Triglia
- Epigenetic Regulation and Chromatin Architecture Group, Berlin Institute for Medical Systems Biology, Max-Delbrück Centre for Molecular Medicine, Berlin, Germany
| | - Alexander Kukalev
- Epigenetic Regulation and Chromatin Architecture Group, Berlin Institute for Medical Systems Biology, Max-Delbrück Centre for Molecular Medicine, Berlin, Germany
| | - Carmelo Ferrai
- Epigenetic Regulation and Chromatin Architecture Group, Berlin Institute for Medical Systems Biology, Max-Delbrück Centre for Molecular Medicine, Berlin, Germany.,Genome Function Group, MRC Clinical Sciences Centre, Imperial College London, London, United Kingdom
| | - Mita Chotalia
- Genome Function Group, MRC Clinical Sciences Centre, Imperial College London, London, United Kingdom
| | - Emily Brookes
- Genome Function Group, MRC Clinical Sciences Centre, Imperial College London, London, United Kingdom
| | - Hiroshi Kimura
- Graduate School of Bioscience and Biotechnology, Tokyo Institute of Technology, Yokohama, Japan
| | - Ana Pombo
- Epigenetic Regulation and Chromatin Architecture Group, Berlin Institute for Medical Systems Biology, Max-Delbrück Centre for Molecular Medicine, Berlin, Germany.,Genome Function Group, MRC Clinical Sciences Centre, Imperial College London, London, United Kingdom
| |
Collapse
|
31
|
Harmanci A, Rozowsky J, Gerstein M. MUSIC: identification of enriched regions in ChIP-Seq experiments using a mappability-corrected multiscale signal processing framework. Genome Biol 2015; 15:474. [PMID: 25292436 PMCID: PMC4234855 DOI: 10.1186/s13059-014-0474-3] [Citation(s) in RCA: 53] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2014] [Indexed: 12/20/2022] Open
Abstract
We present MUSIC, a signal processing approach for identification of enriched regions in ChIP-Seq data, available at music.gersteinlab.org. MUSIC first filters the ChIP-Seq read-depth signal for systematic noise from non-uniform mappability, which fragments enriched regions. Then it performs a multiscale decomposition, using median filtering, identifying enriched regions at multiple length scales. This is useful given the wide range of scales probed in ChIP-Seq assays. MUSIC performs favorably in terms of accuracy and reproducibility compared with other methods. In particular, analysis of RNA polymerase II data reveals a clear distinction between the stalled and elongating forms of the polymerase.
Collapse
Affiliation(s)
- Arif Harmanci
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT 06520, USA
| | | | | |
Collapse
|
32
|
Hansen P, Hecht J, Ibrahim DM, Krannich A, Truss M, Robinson PN. Saturation analysis of ChIP-seq data for reproducible identification of binding peaks. Genome Res 2015; 25:1391-400. [PMID: 26163319 PMCID: PMC4561497 DOI: 10.1101/gr.189894.115] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2015] [Accepted: 07/06/2015] [Indexed: 11/24/2022]
Abstract
Chromatin immunoprecipitation coupled with next-generation sequencing (ChIP-seq) is a powerful technology to identify the genome-wide locations of transcription factors and other DNA binding proteins. Computational ChIP-seq peak calling infers the location of protein–DNA interactions based on various measures of enrichment of sequence reads. In this work, we introduce an algorithm, Q, that uses an assessment of the quadratic enrichment of reads to center candidate peaks followed by statistical analysis of saturation of candidate peaks by 5′ ends of reads. We show that our method not only is substantially faster than several competing methods but also demonstrates statistically significant advantages with respect to reproducibility of results and in its ability to identify peaks with reproducible binding site motifs. We show that Q has superior performance in the delineation of double RNAPII and H3K4me3 peaks surrounding transcription start sites related to a better ability to resolve individual peaks. The method is implemented in C+l+ and is freely available under an open source license.
Collapse
Affiliation(s)
- Peter Hansen
- Institute for Medical and Human Genetics, Charité-Universitätsmedizin Berlin, 13353 Berlin, Germany; Berlin Brandenburg Center for Regenerative Therapies (BCRT), Charité-Universitätsmedizin Berlin, 13353 Berlin, Germany
| | - Jochen Hecht
- Institute for Medical and Human Genetics, Charité-Universitätsmedizin Berlin, 13353 Berlin, Germany; Berlin Brandenburg Center for Regenerative Therapies (BCRT), Charité-Universitätsmedizin Berlin, 13353 Berlin, Germany; Max Planck Institute for Molecular Genetics, 14195 Berlin, Germany
| | - Daniel M Ibrahim
- Institute for Medical and Human Genetics, Charité-Universitätsmedizin Berlin, 13353 Berlin, Germany; Max Planck Institute for Molecular Genetics, 14195 Berlin, Germany
| | - Alexander Krannich
- Department of Biostatistics, Clinical Research Unit, Berlin Institute of Health, Charité-Universitätsmedizin Berlin, 13353 Berlin, Germany
| | - Matthias Truss
- Labor für Pädiatrische Molekularbiologie, Charité-Universitätsmedizin Berlin, 10117, Berlin, Germany
| | - Peter N Robinson
- Institute for Medical and Human Genetics, Charité-Universitätsmedizin Berlin, 13353 Berlin, Germany; Berlin Brandenburg Center for Regenerative Therapies (BCRT), Charité-Universitätsmedizin Berlin, 13353 Berlin, Germany; Max Planck Institute for Molecular Genetics, 14195 Berlin, Germany; Institute for Bioinformatics, Department of Mathematics and Computer Science, Freie Universität Berlin, 14195 Berlin, Germany
| |
Collapse
|
33
|
Abstract
The Forkhead box O (Foxo) family of transcription factors has a critical role in controlling the development, differentiation, and function of T cells. However, the direct target genes of Foxo transcription factors in T cells have not been well characterized. In this study, we focused on mapping the genome wide Foxo1-binding sites in naïve CD4+ T cells, CD8+ T cells, and Foxp3+ regulatory T (Treg) cells. By using chromatin immunoprecipitation coupled with deep sequencing (ChIP-Seq), we identified Foxo1 binding sites that were shared among or specific to the three T cell populations. Here we describe the experiments, quality controls, as well as the deep sequencing data. Part of the data analysis has been published by Ouyang W et al. in Nature 2012[1] and Kim MV et al. in Immunity 2013[2], and the associated data set were uploaded to NCBI Gene Expression Omnibus.
Collapse
|
34
|
Ibrahim MM, Lacadie SA, Ohler U. JAMM: a peak finder for joint analysis of NGS replicates. ACTA ACUST UNITED AC 2014; 31:48-55. [PMID: 25223640 DOI: 10.1093/bioinformatics/btu568] [Citation(s) in RCA: 47] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
MOTIVATION Although peak finding in next-generation sequencing (NGS) datasets has been addressed extensively, there is no consensus on how to analyze and process biological replicates. Furthermore, most peak finders do not focus on accurate determination of enrichment site widths and are not widely applicable to different types of datasets. RESULTS We developed JAMM (Joint Analysis of NGS replicates via Mixture Model clustering): a peak finder that can integrate information from biological replicates, determine enrichment site widths accurately and resolve neighboring narrow peaks. JAMM is a universal peak finder that is applicable to different types of datasets. We show that JAMM is among the best performing peak finders in terms of site detection accuracy and in terms of accurate determination of enrichment sites widths. In addition, JAMM's replicate integration improves peak spatial resolution, sorting and peak finding accuracy. AVAILABILITY AND IMPLEMENTATION JAMM is available for free and can run on Linux machines through the command line: http://code.google.com/p/jamm-peak-finder.
Collapse
Affiliation(s)
- Mahmoud M Ibrahim
- Department of Biology, Humboldt University, Invalidenstrasse 43, D-10115 Berlin, Germany and The Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine Berlin-Buch, Robert Rössle Str. 10, Berlin 13125, Germany Department of Biology, Humboldt University, Invalidenstrasse 43, D-10115 Berlin, Germany and The Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine Berlin-Buch, Robert Rössle Str. 10, Berlin 13125, Germany
| | - Scott A Lacadie
- Department of Biology, Humboldt University, Invalidenstrasse 43, D-10115 Berlin, Germany and The Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine Berlin-Buch, Robert Rössle Str. 10, Berlin 13125, Germany
| | - Uwe Ohler
- Department of Biology, Humboldt University, Invalidenstrasse 43, D-10115 Berlin, Germany and The Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine Berlin-Buch, Robert Rössle Str. 10, Berlin 13125, Germany Department of Biology, Humboldt University, Invalidenstrasse 43, D-10115 Berlin, Germany and The Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine Berlin-Buch, Robert Rössle Str. 10, Berlin 13125, Germany
| |
Collapse
|
35
|
|
36
|
Chan HP, Lai TL. Discussion on “Change-Points: From Sequential Detection to Biology and Back” by David Siegmund. Seq Anal 2013. [DOI: 10.1080/07474946.2013.751840] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|