An integrated 3-Dimensional Genome Modeling Engine for data-driven simulation of spatial genome organization

• Published in: Genome research 2016-12, Vol.26 (12), p.1697-1709
• Main Authors: Szałaj, Przemysław, Tang, Zhonghui, Michalski, Paul, Pietal, Michal J, Luo, Oscar J, Sadowski, Michał, Li, Xingwang, Radew, Kamen, Ruan, Yijun, Plewczynski, Dariusz
• Format: Article
• Language: English
• Subjects: B-Lymphocytes - cytology; Cells, Cultured; Chromatin - genetics; Chromosomes; Chromosomes, Human - genetics; Computational Biology - methods; Computer Simulation; Gene Expression Regulation; Humans; Imaging, Three-Dimensional - methods; Information Storage and Retrieval; Method; Models, Genetic
• ISSN: 1088-9051; 1549-5469
• DOI: 10.1101/gr.205062.116

ChIA-PET is a high-throughput mapping technology that reveals long-range chromatin interactions and provides insights into the basic principles of spatial genome organization and gene regulation mediated by specific protein factors. Recently, we showed that a single ChIA-PET experiment provides information at all genomic scales of interest, from the high-resolution locations of binding sites and enriched chromatin interactions mediated by specific protein factors, to the low resolution of nonenriched interactions that reflect topological neighborhoods of higher-order chromosome folding. This multilevel nature of ChIA-PET data offers an opportunity to use multiscale 3D models to study structural-functional relationships at multiple length scales, but doing so requires a structural modeling platform. Here, we report the development of 3D-GNOME (3-Dimensional Genome Modeling Engine), a complete computational pipeline for 3D simulation using ChIA-PET data. 3D-GNOME consists of three integrated components: a graph-distance-based heat map normalization tool, a 3D modeling platform, and an interactive 3D visualization tool. Using ChIA-PET and Hi-C data derived from human B-lymphocytes, we demonstrate the effectiveness of 3D-GNOME in building 3D genome models at multiple levels, including the entire genome, individual chromosomes, and specific segments at megabase (Mb) and kilobase (kb) resolutions of single average and ensemble structures. Further incorporation of CTCF-motif orientation and high-resolution looping patterns in 3D simulation provided additional reliability of potential biologically plausible topological structures.