3DeFDR: Identifying cell type-specific looping interactions with empirical false discovery rate guided thresholding

The mammalian genome is connected into tens of thousands of long-range looping interactions critically linked to spatiotemporal gene expression regulation. An important unanswered question is to what extent looping interactions change across developmental models, genetic perturbations, drug treatments, and disease states. Although methods exist for calling loops in single biological conditions, there is a severe shortage of computational tools for rigorous assessment of cell type-specific looping interactions across multiple biological conditions. Here we present 3DeFDR, a simple and effective statistical tool for classifying dynamic looping interactions across biological conditions from Chromosome-Conformation-Capture-Carbon-Copy (5C) data. 3DeFDR parses chromatin contacts into invariant and cell type-specific classes by thresholding on differences in modeled interaction strength signal across two or three cellular states. Thresholds are iteratively adjusted based on a target empirical false discovery rate computed between real and simulated 5C maps. 3DeFDR enables the sensitive detection of high-confidence looping interactions and markedly reduces false positives when benchmarked against a classic analysis of variance (ANOVA) test, our newly formulated parametric likelihood ratio test (3DLRT), and the leading Hi-C differential interaction caller diffHic. 3DeFDR also sensitively and specifically calls loops in Mb-scale genomic regions parsed from Hi-C data. Our work provides a statistical framework and an open-source coding library for identifying dynamic long-range looping interactions in high-resolution 5C data from multiple cellular conditions.

[1]  Michael Q. Zhang,et al.  FIND: difFerential chromatin INteractions Detection using a spatial Poisson process , 2018, Genome research.

[2]  J. Dekker,et al.  Capturing Chromosome Conformation , 2002, Science.

[3]  S. Mundlos,et al.  Formation of new chromatin domains determines pathogenicity of genomic duplications , 2016, Nature.

[4]  Neva C. Durand,et al.  Juicer Provides a One-Click System for Analyzing Loop-Resolution Hi-C Experiments. , 2016, Cell systems.

[5]  Michael P Snyder,et al.  Static and dynamic DNA loops form AP-1 bound activation hubs during macrophage development , 2017, bioRxiv.

[6]  L. Mirny,et al.  Biological techniques: Chromosomes captured one by one , 2013, Nature.

[7]  Jing Liang,et al.  Chromatin architecture reorganization during stem cell differentiation , 2015, Nature.

[8]  Job Dekker,et al.  Invariant TAD Boundaries Constrain Cell-Type-Specific Looping Interactions between Promoters and Distal Elements around the CFTR Locus. , 2016, American journal of human genetics.

[9]  S. Bicciato,et al.  Comparison of computational methods for Hi-C data analysis , 2017, Nature Methods.

[10]  Nuno A. Fonseca,et al.  Two independent modes of chromatin organization revealed by cohesin removal , 2017, Nature.

[11]  M. Gobbi,et al.  Analysis of hundreds of cis-regulatory landscapes at high resolution in a single, high-throughput experiment , 2014, Nature Genetics.

[12]  J. Dekker,et al.  The long-range interaction landscape of gene promoters , 2012, Nature.

[13]  Borbala Mifsud,et al.  GOTHiC, a probabilistic model to resolve complex biases and to identify real interactions in Hi-C data , 2017, PloS one.

[14]  Thomas G. Gilgenast,et al.  Systematic Evaluation of Statistical Methods for Identifying Looping Interactions in 5C Data. , 2019, Cell systems.

[15]  P. Flicek,et al.  Molecular maps of the reorganization of genome-nuclear lamina interactions during differentiation. , 2010, Molecular cell.

[16]  Peter J. Park,et al.  The 4D Nucleome Project , 2017 .

[17]  Michael P Snyder,et al.  ChIA-PET2: a versatile and flexible pipeline for ChIA-PET data analysis , 2016, Nucleic acids research.

[18]  J. Telenius,et al.  Multiplexed analysis of chromosome conformation at vastly improved sensitivity , 2015, Nature Methods.

[19]  Qi Zheng,et al.  HIPPIE: a high-throughput identification pipeline for promoter interacting enhancer elements , 2015, Bioinform..

[20]  Shawn M. Gillespie,et al.  Insulator dysfunction and oncogene activation in IDH mutant gliomas , 2015, Nature.

[21]  Michael Q. Zhang,et al.  MICC: an R package for identifying chromatin interactions from ChIA-PET data , 2015, Bioinform..

[22]  A. Tanay,et al.  Probabilistic modeling of Hi-C contact maps eliminates systematic biases to characterize global chromosomal architecture , 2011, Nature Genetics.

[23]  B. Lenhard,et al.  GenomicInteractions: An R/Bioconductor package for manipulating and investigating chromatin interaction data , 2015, BMC Genomics.

[24]  Jennifer E. Phillips-Cremins,et al.  Architectural Protein Subclasses Shape 3D Organization of Genomes during Lineage Commitment , 2013, Cell.

[25]  Daniel S. Day,et al.  YY1 Is a Structural Regulator of Enhancer-Promoter Loops , 2017, Cell.

[26]  B. Steensel,et al.  Nuclear organization of active and inactive chromatin domains uncovered by chromosome conformation capture–on-chip (4C) , 2006, Nature Genetics.

[27]  A. Tanay,et al.  Multiscale 3D Genome Rewiring during Mouse Neural Development , 2017, Cell.

[28]  L. Mirny,et al.  Iterative Correction of Hi-C Data Reveals Hallmarks of Chromosome Organization , 2012, Nature Methods.

[29]  Neva C. Durand,et al.  A 3D Map of the Human Genome at Kilobase Resolution Reveals Principles of Chromatin Looping , 2014, Cell.

[30]  Daniel S. Day,et al.  Activation of proto-oncogenes by disruption of chromosome neighborhoods , 2015, Science.

[31]  Lusy Handoko,et al.  Dynamic Reorganization of Extremely Long-Range Promoter-Promoter Interactions between Two States of Pluripotency. , 2015, Cell stem cell.

[32]  A. Tanay,et al.  Single cell Hi-C reveals cell-to-cell variability in chromosome structure , 2013, Nature.

[33]  Michael P. Snyder,et al.  Mango: a bias-correcting ChIA-PET analysis pipeline , 2015, Bioinform..

[34]  L. Wessels,et al.  Domain organization of human chromosomes revealed by mapping of nuclear lamina interactions , 2008, Nature.

[35]  J. Dekker,et al.  Genomics tools for the unraveling of chromosome architecture , 2010, Nature Biotechnology.

[36]  Geir Kjetil Sandve,et al.  HiBrowse: multi-purpose statistical analysis of genome-wide chromatin 3D organization , 2014, Bioinform..

[37]  Jesse R. Dixon,et al.  Topological Domains in Mammalian Genomes Identified by Analysis of Chromatin Interactions , 2012, Nature.

[38]  Jill M Dowen,et al.  Control of Cell Identity Genes Occurs in Insulated Neighborhoods in Mammalian Chromosomes , 2014, Cell.

[39]  Wouter de Laat,et al.  3C-based technologies to study the shape of the genome. , 2012, Methods.

[40]  K. Sandhu,et al.  Circular chromosome conformation capture (4C) uncovers extensive networks of epigenetically regulated intra- and interchromosomal interactions , 2006, Nature Genetics.

[41]  Martin J. Aryee,et al.  hichipper: a preprocessing pipeline for calling DNA loops from HiChIP data , 2018, Nature Methods.

[42]  Pedro P. Rocha,et al.  CTCF establishes discrete functional chromatin domains at the Hox clusters during differentiation , 2015, Science.

[43]  Gaelen T. Hess,et al.  Static and dynamic DNA loops form AP-1 bound activation hubs during macrophage development , 2017, bioRxiv.

[44]  William Stafford Noble,et al.  Statistical confidence estimation for Hi-C data reveals regulatory chromatin contacts , 2014, Genome research.

[45]  L. Mirny,et al.  Targeted Degradation of CTCF Decouples Local Insulation of Chromosome Domains from Genomic Compartmentalization , 2017, Cell.

[46]  Jonathan M. Cairns,et al.  Lineage-Specific Genome Architecture Links Enhancers and Non-coding Disease Variants to Target Gene Promoters , 2016, Cell.

[47]  C. Nusbaum,et al.  Chromosome Conformation Capture Carbon Copy (5C): a massively parallel solution for mapping interactions between genomic elements. , 2006, Genome research.

[48]  I. Amit,et al.  Comprehensive mapping of long range interactions reveals folding principles of the human genome , 2011 .

[49]  Erez Lieberman Aiden,et al.  Cohesin Loss Eliminates All Loop Domains , 2017, Cell.

[50]  Philip A. Ewels,et al.  Mapping long-range promoter contacts in human cells with high-resolution capture Hi-C , 2015, Nature Genetics.

[51]  Jennifer E. Phillips-Cremins,et al.  YY1 and CTCF orchestrate a 3D chromatin looping switch during early neural lineage commitment , 2017, Genome research.

[52]  Aaron T. L. Lun,et al.  diffHic: a Bioconductor package to detect differential genomic interactions in Hi-C data , 2015, BMC Bioinformatics.

[53]  A. Visel,et al.  Disruptions of Topological Chromatin Domains Cause Pathogenic Rewiring of Gene-Enhancer Interactions , 2015, Cell.

[54]  Martin J. Aryee,et al.  diffloop: a computational framework for identifying and analyzing differential DNA loops from sequencing data , 2016, bioRxiv.