1
|
Liu E, Lyu H, Liu Y, Fu L, Cheng X, Yin X. Identifying TAD-like domains on single-cell Hi-C data by graph embedding and changepoint detection. Bioinformatics 2024; 40:btae138. [PMID: 38449288 PMCID: PMC10960928 DOI: 10.1093/bioinformatics/btae138] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2023] [Revised: 01/10/2024] [Accepted: 03/05/2024] [Indexed: 03/08/2024] Open
Abstract
MOTIVATION Topologically associating domains (TADs) are fundamental building blocks of 3D genome. TAD-like domains in single cells are regarded as the underlying genesis of TADs discovered in bulk cells. Understanding the organization of TAD-like domains helps to get deeper insights into their regulatory functions. Unfortunately, it remains a challenge to identify TAD-like domains on single-cell Hi-C data due to its ultra-sparsity. RESULTS We propose scKTLD, an in silico tool for the identification of TAD-like domains on single-cell Hi-C data. It takes Hi-C contact matrix as the adjacency matrix for a graph, embeds the graph structures into a low-dimensional space with the help of sparse matrix factorization followed by spectral propagation, and the TAD-like domains can be identified using a kernel-based changepoint detection in the embedding space. The results tell that our scKTLD is superior to the other methods on the sparse contact matrices, including downsampled bulk Hi-C data as well as simulated and experimental single-cell Hi-C data. Besides, we demonstrated the conservation of TAD-like domain boundaries at single-cell level apart from heterogeneity within and across cell types, and found that the boundaries with higher frequency across single cells are more enriched for architectural proteins and chromatin marks, and they preferentially occur at TAD boundaries in bulk cells, especially at those with higher hierarchical levels. AVAILABILITY AND IMPLEMENTATION scKTLD is freely available at https://github.com/lhqxinghun/scKTLD.
Collapse
Affiliation(s)
- Erhu Liu
- School of Information and Control Engineering, Xi'an University of Architecture and Technology, Xi'an 710055, China
| | - Hongqiang Lyu
- School of Automation Science and Engineering, Faculty of Electronic and Information Engineering, Xi’an Jiaotong University, Xi'an 710049, China
| | - Yuan Liu
- School of Automation Science and Engineering, Faculty of Electronic and Information Engineering, Xi’an Jiaotong University, Xi'an 710049, China
| | - Laiyi Fu
- School of Automation Science and Engineering, Faculty of Electronic and Information Engineering, Xi’an Jiaotong University, Xi'an 710049, China
| | - Xiaoliang Cheng
- Department of Pharmacy, The First Affiliated Hospital of Xi’an Jiaotong University, Xi'an 710061, China
| | - Xiaoran Yin
- Department of Oncology, The Second Affiliated Hospital of Xi’an Jiaotong University, Xi'an 710004, China
| |
Collapse
|
2
|
Korsak S, Plewczynski D. LoopSage: An energy-based Monte Carlo approach for the loop extrusion modeling of chromatin. Methods 2024; 223:106-117. [PMID: 38295892 DOI: 10.1016/j.ymeth.2024.01.015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2023] [Revised: 12/29/2023] [Accepted: 01/10/2024] [Indexed: 02/05/2024] Open
Abstract
The connection between the patterns observed in 3C-type experiments and the modeling of polymers remains unresolved. This paper presents a simulation pipeline that generates thermodynamic ensembles of 3D structures for topologically associated domain (TAD) regions by loop extrusion model (LEM). The simulations consist of two main components: a stochastic simulation phase, employing a Monte Carlo approach to simulate the binding positions of cohesins, and a dynamical simulation phase, utilizing these cohesins' positions to create 3D structures. In this approach, the system's total energy is the combined result of the Monte Carlo energy and the molecular simulation energy, which are iteratively updated. The structural maintenance of chromosomes (SMC) protein complexes are represented as loop extruders, while the CCCTC-binding factor (CTCF) locations on DNA sequence are modeled as energy minima on the Monte Carlo energy landscape. Finally, the spatial distances between DNA segments from ChIA-PET experiments are compared with the computer simulations, and we observe significant Pearson correlations between predictions and the real data. LoopSage model offers a fresh perspective on chromatin loop dynamics, allowing us to observe phase transition between sparse and condensed states in chromatin.
Collapse
Affiliation(s)
- Sevastianos Korsak
- Faculty of Mathematics and Information Science, Warsaw University of Technology, Warsaw, Poland; Center of New Technologies, University of Warsaw, Warsaw, Poland
| | - Dariusz Plewczynski
- Faculty of Mathematics and Information Science, Warsaw University of Technology, Warsaw, Poland; Center of New Technologies, University of Warsaw, Warsaw, Poland.
| |
Collapse
|
3
|
Han MH, Issagulova D, Park M. Interplay between epigenome and 3D chromatin structure. BMB Rep 2023; 56:633-644. [PMID: 38052424 PMCID: PMC10761748] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2023] [Revised: 11/28/2023] [Accepted: 12/05/2023] [Indexed: 12/07/2023] Open
Abstract
Epigenetic mechanisms, primarily mediated through histone and DNA modifications, play a pivotal role in orchestrating the functional identity of a cell and its response to environmental cues. Similarly, the spatial arrangement of chromatin within the threedimensional (3D) nucleus has been recognized as a significant factor influencing genomic function. Investigating the relationship between epigenetic regulation and 3D chromatin structure has revealed correlation and causality between these processes, from the global alignment of average chromatin structure with chromatin marks to the nuanced correlations at smaller scales. This review aims to dissect the biological significance and the interplay between the epigenome and 3D chromatin structure, while also exploring the underlying molecular mechanisms. By synthesizing insights from both experimental and modeling perspectives, we seek to provide a comprehensive understanding of cellular functions. [BMB Reports 2023; 56(12): 633-644].
Collapse
Affiliation(s)
- Man-Hyuk Han
- Department of Biological Sciences, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea
| | - Dariya Issagulova
- Department of Biological Sciences, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea
| | - Minhee Park
- Department of Biological Sciences, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea; Graduate School of Engineering Biology, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141; KAIST Institute for the BioCentury, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141; KAIST Stem Cell Center, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea
| |
Collapse
|
4
|
Li A, Zeng G, Wang H, Li X, Zhang Z. DeDoc2 Identifies and Characterizes the Hierarchy and Dynamics of Chromatin TAD-Like Domains in the Single Cells. ADVANCED SCIENCE (WEINHEIM, BADEN-WURTTEMBERG, GERMANY) 2023; 10:e2300366. [PMID: 37162225 PMCID: PMC10369259 DOI: 10.1002/advs.202300366] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/16/2023] [Revised: 04/18/2023] [Indexed: 05/11/2023]
Abstract
Topologically associating domains (TADs) are functional chromatin units with hierarchical structure. However, the existence, prevalence, and dynamics of such hierarchy in single cells remain unexplored. Here, a new generation TAD-like domain (TLD) detection algorithm, named deDoc2, to decode the hierarchy of TLDs in single cells, is reported. With dynamic programming, deDoc2 seeks genome partitions with global minimal structure entropy for both whole and local contact matrix. Notably, deDoc2 outperforms state-of-the-art tools and is one of only two tools able to identify the hierarchy of TLDs in single cells. By applying deDoc2, it is showed that the hierarchy of TLDs in single cells is highly dynamic during cell cycle, as well as among human brain cortex cells, and that it is associated with cellular identity and functions. Thus, the results demonstrate the abundance of information potentially encoded by TLD hierarchy for functional regulation. The deDoc2 can be freely accessed at https://github.com/zengguangjie/deDoc2.
Collapse
Affiliation(s)
- Angsheng Li
- State Key Laboratory of Software Development EnvironmentSchool of Computer ScienceBeihang UniversityBeijing100191P. R. China
- Zhongguancun LaboratoryBeijing100094P. R. China
| | - Guangjie Zeng
- State Key Laboratory of Software Development EnvironmentSchool of Computer ScienceBeihang UniversityBeijing100191P. R. China
| | - Haoyu Wang
- CAS Key Laboratory of Genome Sciences and InformationBeijing Institute of GenomicsChinese Academy of Sciences and China National Center for BioinformationBeijing100101China
- School of Life ScienceUniversity of Chinese Academy of SciencesBeijing101408P. R. China
| | - Xiao Li
- CAS Key Laboratory of Genome Sciences and InformationBeijing Institute of GenomicsChinese Academy of Sciences and China National Center for BioinformationBeijing100101China
- School of Life ScienceUniversity of Chinese Academy of SciencesBeijing101408P. R. China
| | - Zhihua Zhang
- CAS Key Laboratory of Genome Sciences and InformationBeijing Institute of GenomicsChinese Academy of Sciences and China National Center for BioinformationBeijing100101China
- School of Life ScienceUniversity of Chinese Academy of SciencesBeijing101408P. R. China
| |
Collapse
|
5
|
Esposito A, Abraham A, Conte M, Vercellone F, Prisco A, Bianco S, Chiariello AM. The Physics of DNA Folding: Polymer Models and Phase-Separation. Polymers (Basel) 2022; 14:1918. [PMID: 35567087 PMCID: PMC9104579 DOI: 10.3390/polym14091918] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2022] [Revised: 04/23/2022] [Accepted: 04/27/2022] [Indexed: 02/04/2023] Open
Abstract
Within cell nuclei, several biophysical processes occur in order to allow the correct activities of the genome such as transcription and gene regulation. To quantitatively investigate such processes, polymer physics models have been developed to unveil the molecular mechanisms underlying genome functions. Among these, phase-separation plays a key role since it controls gene activity and shapes chromatin spatial structure. In this paper, we review some recent experimental and theoretical progress in the field and show that polymer physics in synergy with numerical simulations can be helpful for several purposes, including the study of molecular condensates, gene-enhancer dynamics, and the three-dimensional reconstruction of real genomic regions.
Collapse
Affiliation(s)
- Andrea Esposito
- Dipartimento di Fisica, Università di Napoli Federico II, INFN Napoli, Complesso Universitario di Monte Sant’Angelo, 80126 Naples, Italy; (A.E.); (A.A.); (M.C.); (F.V.)
| | - Alex Abraham
- Dipartimento di Fisica, Università di Napoli Federico II, INFN Napoli, Complesso Universitario di Monte Sant’Angelo, 80126 Naples, Italy; (A.E.); (A.A.); (M.C.); (F.V.)
| | - Mattia Conte
- Dipartimento di Fisica, Università di Napoli Federico II, INFN Napoli, Complesso Universitario di Monte Sant’Angelo, 80126 Naples, Italy; (A.E.); (A.A.); (M.C.); (F.V.)
| | - Francesca Vercellone
- Dipartimento di Fisica, Università di Napoli Federico II, INFN Napoli, Complesso Universitario di Monte Sant’Angelo, 80126 Naples, Italy; (A.E.); (A.A.); (M.C.); (F.V.)
| | | | - Simona Bianco
- Dipartimento di Fisica, Università di Napoli Federico II, INFN Napoli, Complesso Universitario di Monte Sant’Angelo, 80126 Naples, Italy; (A.E.); (A.A.); (M.C.); (F.V.)
- Berlin Institute for Medical Systems Biology, Max-Delbrück Centre (MDC) for Molecular Medicine, 10115 Berlin, Germany
| | - Andrea M. Chiariello
- Dipartimento di Fisica, Università di Napoli Federico II, INFN Napoli, Complesso Universitario di Monte Sant’Angelo, 80126 Naples, Italy; (A.E.); (A.A.); (M.C.); (F.V.)
| |
Collapse
|
6
|
Serna-Pujol N, Salinas-Pena M, Mugianesi F, Le Dily F, Marti-Renom MA, Jordan A. Coordinated changes in gene expression, H1 variant distribution and genome 3D conformation in response to H1 depletion. Nucleic Acids Res 2022; 50:3892-3910. [PMID: 35380694 PMCID: PMC9023279 DOI: 10.1093/nar/gkac226] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2021] [Revised: 03/21/2022] [Accepted: 03/28/2022] [Indexed: 11/12/2022] Open
Abstract
Up to seven members of the histone H1 family may contribute to chromatin compaction and its regulation in human somatic cells. In breast cancer cells, knock-down of multiple H1 variants deregulates many genes, promotes the appearance of genome-wide accessibility sites and triggers an interferon response via activation of heterochromatic repeats. However, how these changes in the expression profile relate to the re-distribution of H1 variants as well as to genome conformational changes have not been yet studied. Here, we combined ChIP-seq of five endogenous H1 variants with Chromosome Conformation Capture analysis in wild-type and H1.2/H1.4 knock-down T47D cells. The results indicate that H1 variants coexist in the genome in two large groups depending on the local GC content and that their distribution is robust with respect to H1 depletion. Despite the small changes in H1 variants distribution, knock-down of H1 translated into more isolated but de-compacted chromatin structures at the scale of topologically associating domains (TADs). Such changes in TAD structure correlated with a coordinated gene expression response of their resident genes. This is the first report describing simultaneous profiling of five endogenous H1 variants and giving functional evidence of genome topology alterations upon H1 depletion in human cancer cells.
Collapse
Affiliation(s)
- Núria Serna-Pujol
- Molecular Biology Institute of Barcelona (IBMB-CSIC), Barcelona, 08028 Spain
| | - Mónica Salinas-Pena
- Molecular Biology Institute of Barcelona (IBMB-CSIC), Barcelona, 08028 Spain
| | - Francesca Mugianesi
- CNAG-CRG, Centre for Genomic Regulation, The Barcelona Institute of Science and Technology, Baldiri Reixac 4, Barcelona 08028, Spain
| | - François Le Dily
- Centre for Genomic Regulation, The Barcelona Institute for Science and Technology, Carrer del Doctor Aiguader 88, Barcelona 08003, Spain
| | - Marc A Marti-Renom
- CNAG-CRG, Centre for Genomic Regulation, The Barcelona Institute of Science and Technology, Baldiri Reixac 4, Barcelona 08028, Spain.,Centre for Genomic Regulation, The Barcelona Institute for Science and Technology, Carrer del Doctor Aiguader 88, Barcelona 08003, Spain.,Pompeu Fabra University, Doctor Aiguader 88, Barcelona 08003, Spain.,ICREA, Pg. Lluís Companys 23, 08010 Barcelona, Spain
| | - Albert Jordan
- Molecular Biology Institute of Barcelona (IBMB-CSIC), Barcelona, 08028 Spain
| |
Collapse
|
7
|
Di Stefano M, Nützmann HW. Modeling the 3D genome of plants. Nucleus 2021; 12:65-81. [PMID: 34057011 PMCID: PMC8168717 DOI: 10.1080/19491034.2021.1927503] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2020] [Revised: 04/29/2021] [Accepted: 05/05/2021] [Indexed: 12/11/2022] Open
Abstract
Chromosomes are the carriers of inheritable traits and define cell function and development. This is not only based on the linear DNA sequence of chromosomes but also on the additional molecular information they are associated with, including the transcription machinery, histone modifications, and their three-dimensional folding. The synergistic application of experimental approaches and computer simulations has helped to unveil how these organizational layers of the genome interplay in various organisms. However, such multidisciplinary approaches are still rarely explored in the plant kingdom. Here, we provide an overview of our current knowledge on plant 3D genome organization and review recent efforts to integrate cutting-edge experiments from microscopy and next-generation sequencing approaches with theoretical models. Building on these recent approaches, we propose possible avenues to extend the application of theoretical modeling in the characterization of the 3D genome organization in plants.
Collapse
Affiliation(s)
- Marco Di Stefano
- Institute of Human Genetics, Centre National de la Recherche Scientifique, University of Montpellier, Montpellier, France
| | - Hans-Wilhelm Nützmann
- The Milner Centre for Evolution, Department of Biology and Biochemistry, University of Bath, Bath, UK
| |
Collapse
|
8
|
Liang J, Perez-Rathke A. Minimalistic 3D chromatin models: Sparse interactions in single cells drive the chromatin fold and form many-body units. Curr Opin Struct Biol 2021; 71:200-214. [PMID: 34399301 DOI: 10.1016/j.sbi.2021.06.017] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2021] [Revised: 06/27/2021] [Accepted: 06/29/2021] [Indexed: 11/26/2022]
Abstract
Computational three-dimensional chromatin modeling has helped uncover principles of genome organization. Here, we discuss methods for modeling three-dimensional chromatin structures, with focus on a minimalistic polymer model which inverts population Hi-C into single-cell conformations. Utilizing only basic physical properties, this model reveals that a few specific Hi-C interactions can fold chromatin into conformations consistent with single-cell imaging, Dip-C, and FISH measurements. Aggregated single-cell chromatin conformations also reproduce Hi-C frequencies. This approach allows quantification of structural heterogeneity and discovery of many-body interaction units and has revealed additional insights, including (1) topologically associating domains as a byproduct of folding driven by specific interactions, (2) cell subpopulations with different structural scaffolds are developmental stage dependent, and (3) the functional landscape of many-body units within enhancer-rich regions. We also discuss these findings in relation to the genome structure-function relationship.
Collapse
Affiliation(s)
- Jie Liang
- Center for Bioinformatics and Quantitative Biology & Richard and Loan Hill Department of Bioengineering, University of Illinois at Chicago, Chicago, IL, 60612, USA.
| | - Alan Perez-Rathke
- Center for Bioinformatics and Quantitative Biology & Richard and Loan Hill Department of Bioengineering, University of Illinois at Chicago, Chicago, IL, 60612, USA
| |
Collapse
|
9
|
Li X, Zeng G, Li A, Zhang Z. DeTOKI identifies and characterizes the dynamics of chromatin TAD-like domains in a single cell. Genome Biol 2021; 22:217. [PMID: 34311744 PMCID: PMC8314462 DOI: 10.1186/s13059-021-02435-7] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2021] [Accepted: 07/14/2021] [Indexed: 12/19/2022] Open
Abstract
Topologically associating domains (TAD) are a key structure of the 3D mammalian genomes. However, the prevalence and dynamics of TAD-like domains in single cells remain elusive. Here we develop a new algorithm, named deTOKI, to decode TAD-like domains with single-cell Hi-C data. By non-negative matrix factorization, deTOKI seeks regions that insulate the genome into blocks with minimal chance of clustering. deTOKI outperforms competing tools and reliably identifies TAD-like domains in single cells. Finally, we find that TAD-like domains are not only prevalent, but also subject to tight regulation in single cells.
Collapse
Affiliation(s)
- Xiao Li
- CAS Key Laboratory of Genome Sciences and Information, Beijing Institute of Genomics, Chinese Academy of Sciences, and China National Center for Bioinformation, Beijing, 100101, China
- School of Life Science, University of Chinese Academy of Sciences, Beijing, People's Republic of China
- School of Artificial Intelligence, University of Chinese Academy of Sciences, Beijing, People's Republic of China
| | - Guangjie Zeng
- State Key Laboratory of Software Development Environment, School of Computer Science, Beihang University, 100083, Beijing, People's Republic of China
| | - Angsheng Li
- State Key Laboratory of Software Development Environment, School of Computer Science, Beihang University, 100083, Beijing, People's Republic of China
| | - Zhihua Zhang
- CAS Key Laboratory of Genome Sciences and Information, Beijing Institute of Genomics, Chinese Academy of Sciences, and China National Center for Bioinformation, Beijing, 100101, China.
- School of Life Science, University of Chinese Academy of Sciences, Beijing, People's Republic of China.
- School of Artificial Intelligence, University of Chinese Academy of Sciences, Beijing, People's Republic of China.
| |
Collapse
|
10
|
Polymer models are a versatile tool to study chromatin 3D organization. Biochem Soc Trans 2021; 49:1675-1684. [PMID: 34282837 DOI: 10.1042/bst20201004] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2021] [Revised: 06/21/2021] [Accepted: 06/25/2021] [Indexed: 11/17/2022]
Abstract
The development of new experimental technologies is opening the way to a deeper investigation of the three-dimensional organization of chromosomes inside the cell nucleus. Genome architecture is linked to vital functional purposes, yet a full comprehension of the mechanisms behind DNA folding is still far from being accomplished. Theoretical approaches based on polymer physics have been employed to understand the complexity of chromatin architecture data and to unveil the basic mechanisms shaping its structure. Here, we review some recent advances in the field to discuss how Polymer Physics, combined with numerical Molecular Dynamics simulation and Machine Learning based inference, can capture important aspects of genome organization, including the description of tissue-specific structural rearrangements, the detection of novel, regulatory-linked architectural elements and the structural variability of chromatin at the single-cell level.
Collapse
|
11
|
Abstract
The spatial organization of the genome in the cell nucleus is pivotal to cell function. However, how the 3D genome organization and its dynamics influence cellular phenotypes remains poorly understood. The very recent development of single-cell technologies for probing the 3D genome, especially single-cell Hi-C (scHi-C), has ushered in a new era of unveiling cell-to-cell variability of 3D genome features at an unprecedented resolution. Here, we review recent developments in computational approaches to the analysis of scHi-C, including data processing, dimensionality reduction, imputation for enhancing data quality, and the revealing of 3D genome features at single-cell resolution. While much progress has been made in computational method development to analyze single-cell 3D genomes, substantial future work is needed to improve data interpretation and multimodal data integration, which are critical to reveal fundamental connections between genome structure and function among heterogeneous cell populations in various biological contexts.
Collapse
Affiliation(s)
- Tianming Zhou
- Computational Biology Department, School of Computer Science, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA;
| | - Ruochi Zhang
- Computational Biology Department, School of Computer Science, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA;
| | - Jian Ma
- Computational Biology Department, School of Computer Science, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA;
| |
Collapse
|
12
|
Lin X, Qi Y, Latham AP, Zhang B. Multiscale modeling of genome organization with maximum entropy optimization. J Chem Phys 2021; 155:010901. [PMID: 34241389 PMCID: PMC8253599 DOI: 10.1063/5.0044150] [Citation(s) in RCA: 34] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2021] [Accepted: 04/28/2021] [Indexed: 12/15/2022] Open
Abstract
Three-dimensional (3D) organization of the human genome plays an essential role in all DNA-templated processes, including gene transcription, gene regulation, and DNA replication. Computational modeling can be an effective way of building high-resolution genome structures and improving our understanding of these molecular processes. However, it faces significant challenges as the human genome consists of over 6 × 109 base pairs, a system size that exceeds the capacity of traditional modeling approaches. In this perspective, we review the progress that has been made in modeling the human genome. Coarse-grained models parameterized to reproduce experimental data via the maximum entropy optimization algorithm serve as effective means to study genome organization at various length scales. They have provided insight into the principles of whole-genome organization and enabled de novo predictions of chromosome structures from epigenetic modifications. Applications of these models at a near-atomistic resolution further revealed physicochemical interactions that drive the phase separation of disordered proteins and dictate chromatin stability in situ. We conclude with an outlook on the opportunities and challenges in studying chromosome dynamics.
Collapse
Affiliation(s)
- Xingcheng Lin
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Yifeng Qi
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Andrew P. Latham
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Bin Zhang
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| |
Collapse
|
13
|
MacKay K, Kusalik A. Computational methods for predicting 3D genomic organization from high-resolution chromosome conformation capture data. Brief Funct Genomics 2021; 19:292-308. [PMID: 32353112 PMCID: PMC7388788 DOI: 10.1093/bfgp/elaa004] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2019] [Revised: 01/30/2020] [Accepted: 02/07/2020] [Indexed: 12/19/2022] Open
Abstract
The advent of high-resolution chromosome conformation capture assays (such as 5C, Hi-C and Pore-C) has allowed for unprecedented sequence-level investigations into the structure-function relationship of the genome. In order to comprehensively understand this relationship, computational tools are required that utilize data generated from these assays to predict 3D genome organization (the 3D genome reconstruction problem). Many computational tools have been developed that answer this need, but a comprehensive comparison of their underlying algorithmic approaches has not been conducted. This manuscript provides a comprehensive review of the existing computational tools (from November 2006 to September 2019, inclusive) that can be used to predict 3D genome organizations from high-resolution chromosome conformation capture data. Overall, existing tools were found to use a relatively small set of algorithms from one or more of the following categories: dimensionality reduction, graph/network theory, maximum likelihood estimation (MLE) and statistical modeling. Solutions in each category are far from maturity, and the breadth and depth of various algorithmic categories have not been fully explored. While the tools for predicting 3D structure for a genomic region or single chromosome are diverse, there is a general lack of algorithmic diversity among computational tools for predicting the complete 3D genome organization from high-resolution chromosome conformation capture data.
Collapse
|
14
|
Zha M, Wang N, Zhang C, Wang Z. Inferring Single-Cell 3D Chromosomal Structures Based on the Lennard-Jones Potential. Int J Mol Sci 2021; 22:ijms22115914. [PMID: 34072879 PMCID: PMC8199262 DOI: 10.3390/ijms22115914] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2021] [Revised: 05/23/2021] [Accepted: 05/28/2021] [Indexed: 11/16/2022] Open
Abstract
Reconstructing three-dimensional (3D) chromosomal structures based on single-cell Hi-C data is a challenging scientific problem due to the extreme sparseness of the single-cell Hi-C data. In this research, we used the Lennard-Jones potential to reconstruct both 500 kb and high-resolution 50 kb chromosomal structures based on single-cell Hi-C data. A chromosome was represented by a string of 500 kb or 50 kb DNA beads and put into a 3D cubic lattice for simulations. A 2D Gaussian function was used to impute the sparse single-cell Hi-C contact matrices. We designed a novel loss function based on the Lennard-Jones potential, in which the ε value, i.e., the well depth, was used to indicate how stable the binding of every pair of beads is. For the bead pairs that have single-cell Hi-C contacts and their neighboring bead pairs, the loss function assigns them stronger binding stability. The Metropolis-Hastings algorithm was used to try different locations for the DNA beads, and simulated annealing was used to optimize the loss function. We proved the correctness and validness of the reconstructed 3D structures by evaluating the models according to multiple criteria and comparing the models with 3D-FISH data.
Collapse
Affiliation(s)
- Mengsheng Zha
- School of Computing Sciences and Computer Engineering, University of Southern Mississippi, 118 College Dr, Hattiesburg, MS 39406, USA; (M.Z.); (C.Z.)
| | - Nan Wang
- Department of Computer Science, New Jersey City University, 2039 Kennedy Blvd, Jersey City, NJ 07305, USA;
| | - Chaoyang Zhang
- School of Computing Sciences and Computer Engineering, University of Southern Mississippi, 118 College Dr, Hattiesburg, MS 39406, USA; (M.Z.); (C.Z.)
| | - Zheng Wang
- Department of Computer Science, University of Miami, 1364 Memorial Drive, Coral Gables, FL 33124, USA
- Correspondence:
| |
Collapse
|
15
|
Dynamics of genome architecture and chromatin function during human B cell differentiation and neoplastic transformation. Nat Commun 2021; 12:651. [PMID: 33510161 PMCID: PMC7844026 DOI: 10.1038/s41467-020-20849-y] [Citation(s) in RCA: 66] [Impact Index Per Article: 22.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2020] [Accepted: 12/21/2020] [Indexed: 02/06/2023] Open
Abstract
To investigate the three-dimensional (3D) genome architecture across normal B cell differentiation and in neoplastic cells from different subtypes of chronic lymphocytic leukemia and mantle cell lymphoma patients, here we integrate in situ Hi-C and nine additional omics layers. Beyond conventional active (A) and inactive (B) compartments, we uncover a highly-dynamic intermediate compartment enriched in poised and polycomb-repressed chromatin. During B cell development, 28% of the compartments change, mostly involving a widespread chromatin activation from naive to germinal center B cells and a reversal to the naive state upon further maturation into memory B cells. B cell neoplasms are characterized by both entity and subtype-specific alterations in 3D genome organization, including large chromatin blocks spanning key disease-specific genes. This study indicates that 3D genome interactions are extensively modulated during normal B cell differentiation and that the genome of B cell neoplasias acquires a tumor-specific 3D genome architecture.
Collapse
|
16
|
Nussinov R, Jang H, Nir G, Tsai CJ, Cheng F. A new precision medicine initiative at the dawn of exascale computing. Signal Transduct Target Ther 2021; 6:3. [PMID: 33402669 PMCID: PMC7785737 DOI: 10.1038/s41392-020-00420-3] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2020] [Revised: 10/27/2020] [Accepted: 10/30/2020] [Indexed: 12/14/2022] Open
Abstract
Which signaling pathway and protein to select to mitigate the patient's expected drug resistance? The number of possibilities facing the physician is massive, and the drug combination should fit the patient status. Here, we briefly review current approaches and data and map an innovative patient-specific strategy to forecast drug resistance targets that centers on parallel (or redundant) proliferation pathways in specialized cells. It considers the availability of each protein in each pathway in the specific cell, its activating mutations, and the chromatin accessibility of its encoding gene. The construction of the resulting Proliferation Pathway Network Atlas will harness the emerging exascale computing and advanced artificial intelligence (AI) methods for therapeutic development. Merging the resulting set of targets, pathways, and proteins, with current strategies will augment the choice for the attending physicians to thwart resistance.
Collapse
Affiliation(s)
- Ruth Nussinov
- Computational Structural Biology Section, Frederick National Laboratory for Cancer Research in the Laboratory of Cancer Immunometabolism, National Cancer Institute, Frederick, MD, 21702, USA.
- Department of Human Molecular Genetics and Biochemistry, Sackler School of Medicine, Tel Aviv University, Tel Aviv, 69978, Israel.
| | - Hyunbum Jang
- Computational Structural Biology Section, Frederick National Laboratory for Cancer Research in the Laboratory of Cancer Immunometabolism, National Cancer Institute, Frederick, MD, 21702, USA
| | - Guy Nir
- Department of Genetics, Harvard Medical School, Boston, MA, 02115, USA
- Wyss Institute for Biologically Inspired Engineering, Harvard University, Boston, MA, 02115, USA
- Department of Biochemistry & Molecular Biology, Department of Neuroscience, Cell Biology and Anatomy, Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, TX, 77555, USA
| | - Chung-Jung Tsai
- Computational Structural Biology Section, Frederick National Laboratory for Cancer Research in the Laboratory of Cancer Immunometabolism, National Cancer Institute, Frederick, MD, 21702, USA
| | - Feixiong Cheng
- Genomic Medicine Institute, Lerner Research Institute, Cleveland Clinic, Cleveland, OH, 44106, USA
- Department of Molecular Medicine, Cleveland Clinic Lerner College of Medicine, Case Western Reserve University, Cleveland, OH, 44195, USA
| |
Collapse
|
17
|
Qi Y, Reyes A, Johnstone SE, Aryee MJ, Bernstein BE, Zhang B. Data-Driven Polymer Model for Mechanistic Exploration of Diploid Genome Organization. Biophys J 2020; 119:1905-1916. [PMID: 33086041 DOI: 10.1016/j.bpj.2020.09.009] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2020] [Revised: 08/24/2020] [Accepted: 09/08/2020] [Indexed: 12/21/2022] Open
Abstract
Chromosomes are positioned nonrandomly inside the nucleus to coordinate with their transcriptional activity. The molecular mechanisms that dictate the global genome organization and the nuclear localization of individual chromosomes are not fully understood. We introduce a polymer model to study the organization of the diploid human genome. It is data-driven because all parameters can be derived from Hi-C data; it is also a mechanistic model because the energy function is explicitly written out based on a few biologically motivated hypotheses. These two features distinguish the model from existing approaches and make it useful both for reconstructing genome structures and for exploring the principles of genome organization. We carried out extensive validations to show that simulated genome structures reproduce a wide variety of experimental measurements, including chromosome radial positions and spatial distances between homologous pairs. Detailed mechanistic investigations support the importance of both specific interchromosomal interactions and centromere clustering for chromosome positioning. We anticipate the polymer model, when combined with Hi-C experiments, to be a powerful tool for investigating large-scale rearrangements in genome structure upon cell differentiation and tumor progression.
Collapse
Affiliation(s)
- Yifeng Qi
- Departments of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts
| | - Alejandro Reyes
- Broad Institute of MIT and Harvard, Cambridge, Massachusetts; Department of Data Sciences, Dana Farber Cancer Institute, Boston, Massachusetts; Department of Biostatistics, Harvard School of Public Health, Boston, Massachusetts
| | - Sarah E Johnstone
- Broad Institute of MIT and Harvard, Cambridge, Massachusetts; Department of Pathology, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts; Center for Cancer Research, Massachusetts General Hospital, Boston, Massachusetts
| | - Martin J Aryee
- Broad Institute of MIT and Harvard, Cambridge, Massachusetts; Department of Pathology, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts; Center for Cancer Research, Massachusetts General Hospital, Boston, Massachusetts
| | - Bradley E Bernstein
- Broad Institute of MIT and Harvard, Cambridge, Massachusetts; Department of Pathology, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts; Center for Cancer Research, Massachusetts General Hospital, Boston, Massachusetts
| | - Bin Zhang
- Departments of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts.
| |
Collapse
|
18
|
Abstract
Chromatin Conformation Capture techniques have unveiled several layers of chromosome organization such as the segregation in compartments, the folding in topologically associating domains (TADs), and site-specific looping interactions. The discovery of this genome hierarchical organization emerged from the computational analysis of chromatin capture data. With the increasing availability of such data, automatic pipelines for the robust comparison, grouping, and classification of multiple experiments are needed. Here we present a pipeline based on the TADbit framework that emphasizes reproducibility, automation, quality check, and statistical robustness. This comprehensive modular pipeline covers all the steps from the sequencing products to the visualization of reconstructed 3D models of the chromatin.
Collapse
|
19
|
Li FZ, Liu ZE, Li XY, Bu LM, Bu HX, Liu H, Zhang CM. Chromatin 3D structure reconstruction with consideration of adjacency relationship among genomic loci. BMC Bioinformatics 2020; 21:272. [PMID: 32611376 PMCID: PMC7329537 DOI: 10.1186/s12859-020-03612-4] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2019] [Accepted: 06/18/2020] [Indexed: 01/19/2023] Open
Abstract
BACKGROUND Chromatin 3D conformation plays important roles in regulating gene or protein functions. High-throughout chromosome conformation capture (3C)-based technologies, such as Hi-C, have been exploited to acquire the contact frequencies among genomic loci at genome-scale. Various computational tools have been proposed to recover the underlying chromatin 3D structures from in situ Hi-C contact map data. As connected residuals in a polymer, neighboring genomic loci have intrinsic mutual dependencies in building a 3D conformation. However, current methods seldom take this feature into account. RESULTS We present a method called ShNeigh, which combines the classical MDS technique with local dependence of neighboring loci modeled by a Gaussian formula, to infer the best 3D structure from noisy and incomplete contact frequency matrices. We validated ShNeigh by comparing it to two typical distance-based algorithms, ShRec3D and ChromSDE. The comparison results on simulated Hi-C dataset showed that, while keeping the high-speed nature of classical MDS, ShNeigh can recover the true structure better than ShRec3D and ChromSDE. Meanwhile, ShNeigh is more robust to data noise. On the publicly available human GM06990 Hi-C data, we demonstrated that the structures reconstructed by ShNeigh are more reproducible between different restriction enzymes than by ShRec3D and ChromSDE, especially at high resolutions manifested by sparse contact maps, which means ShNeigh is more robust to signal coverage. CONCLUSIONS Our method can recover stable structures in high noise and sparse signal settings. It can also reconstruct similar structures from Hi-C data obtained using different restriction enzymes. Therefore, our method provides a new direction for enhancing the reconstruction quality of chromatin 3D structures.
Collapse
Affiliation(s)
- Fang-Zhen Li
- School of Computer Science and Technology, Shandong University of Finance and Economics, Jinan, China. .,Key Laboratory of Machine Learning and Financial Data Mining in Universities of Shandong, Jinan, China.
| | - Zhi-E Liu
- College of Physics and Electronic Engineering, Qilu Normal University, Jinan, China
| | - Xiu-Yuan Li
- School of Computer Science and Technology, Shandong University of Finance and Economics, Jinan, China.,Key Laboratory of Machine Learning and Financial Data Mining in Universities of Shandong, Jinan, China
| | - Li-Mei Bu
- Department of Gastroenterology, Shanghai Pudong Hospital, Fudan University Pudong Medical Center, Shanghai, China
| | - Hong-Xia Bu
- Key Laboratory of Machine Learning and Financial Data Mining in Universities of Shandong, Jinan, China
| | - Hui Liu
- School of Computer Science and Technology, Shandong University of Finance and Economics, Jinan, China.,Digital Media Technology Key Lab of Shandong Province, Jinan, China
| | - Cai-Ming Zhang
- School of Computer Science and Technology, Shandong University of Finance and Economics, Jinan, China.,Digital Media Technology Key Lab of Shandong Province, Jinan, China
| |
Collapse
|
20
|
Stik G, Vidal E, Barrero M, Cuartero S, Vila-Casadesús M, Mendieta-Esteban J, Tian TV, Choi J, Berenguer C, Abad A, Borsari B, le Dily F, Cramer P, Marti-Renom MA, Stadhouders R, Graf T. CTCF is dispensable for immune cell transdifferentiation but facilitates an acute inflammatory response. Nat Genet 2020; 52:655-661. [DOI: 10.1038/s41588-020-0643-0] [Citation(s) in RCA: 61] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2019] [Accepted: 05/08/2020] [Indexed: 11/09/2022]
|
21
|
Di Stefano M, Stadhouders R, Farabella I, Castillo D, Serra F, Graf T, Marti-Renom MA. Transcriptional activation during cell reprogramming correlates with the formation of 3D open chromatin hubs. Nat Commun 2020; 11:2564. [PMID: 32444798 PMCID: PMC7244774 DOI: 10.1038/s41467-020-16396-1] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/26/2019] [Accepted: 05/01/2020] [Indexed: 12/22/2022] Open
Abstract
Chromosome structure is a crucial regulatory factor for a wide range of nuclear processes. Chromosome conformation capture (3C)-based experiments combined with computational modelling are pivotal for unveiling 3D chromosome structure. Here, we introduce TADdyn, a tool that integrates time-course 3C data, restraint-based modelling, and molecular dynamics to simulate the structural rearrangements of genomic loci in a completely data-driven way. We apply TADdyn on in situ Hi-C time-course experiments studying the reprogramming of murine B cells to pluripotent cells, and characterize the structural rearrangements that take place upon changes in the transcriptional state of 21 genomic loci of diverse expression dynamics. By measuring various structural and dynamical properties, we find that during gene activation, the transcription starting site contacts with open and active regions in 3D chromatin domains. We propose that these 3D hubs of open and active chromatin may constitute a general feature to trigger and maintain gene transcription.
Collapse
Affiliation(s)
- Marco Di Stefano
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Baldiri i Reixac 4, 08028, Barcelona, Spain. .,Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Dr. Aiguader 88, 08003, Barcelona, Spain.
| | - Ralph Stadhouders
- Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Dr. Aiguader 88, 08003, Barcelona, Spain.,Department of Pulmonary Medicine and Department of Cell Biology, Erasmus MC, Rotterdam, the Netherlands
| | - Irene Farabella
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Baldiri i Reixac 4, 08028, Barcelona, Spain.,Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Dr. Aiguader 88, 08003, Barcelona, Spain
| | - David Castillo
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Baldiri i Reixac 4, 08028, Barcelona, Spain.,Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Dr. Aiguader 88, 08003, Barcelona, Spain
| | - François Serra
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Baldiri i Reixac 4, 08028, Barcelona, Spain.,Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Dr. Aiguader 88, 08003, Barcelona, Spain.,Computational Biology Group-Barcelona Supercomputing Center (BSC), 08034, Barcelona, Spain
| | - Thomas Graf
- Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Dr. Aiguader 88, 08003, Barcelona, Spain.
| | - Marc A Marti-Renom
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Baldiri i Reixac 4, 08028, Barcelona, Spain. .,Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Dr. Aiguader 88, 08003, Barcelona, Spain. .,Universitat Pompeu Fabra (UPF), 08002, Barcelona, Spain. .,ICREA, Pg. Lluís Companys 23, 08010, Barcelona, Spain.
| |
Collapse
|
22
|
4D Genome Rewiring during Oncogene-Induced and Replicative Senescence. Mol Cell 2020; 78:522-538.e9. [PMID: 32220303 PMCID: PMC7208559 DOI: 10.1016/j.molcel.2020.03.007] [Citation(s) in RCA: 92] [Impact Index Per Article: 23.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2019] [Revised: 12/20/2019] [Accepted: 03/04/2020] [Indexed: 12/26/2022]
Abstract
To understand the role of the extensive senescence-associated 3D genome reorganization, we generated genome-wide chromatin interaction maps, epigenome, replication-timing, whole-genome bisulfite sequencing, and gene expression profiles from cells entering replicative senescence (RS) or upon oncogene-induced senescence (OIS). We identify senescence-associated heterochromatin domains (SAHDs). Differential intra- versus inter-SAHD interactions lead to the formation of senescence-associated heterochromatin foci (SAHFs) in OIS but not in RS. This OIS-specific configuration brings active genes located in genomic regions adjacent to SAHDs in close spatial proximity and favors their expression. We also identify DNMT1 as a factor that induces SAHFs by promoting HMGA2 expression. Upon DNMT1 depletion, OIS cells transition to a 3D genome conformation akin to that of cells in replicative senescence. These data show how multi-omics and imaging can identify critical features of RS and OIS and discover determinants of acute senescence and SAHF formation. Deep multi-omics characterization of replicative and oncogene-induced senescence Senescence-associated heterochromatin domains (SAHDs) form SAHFs via 3D changes DNMT1 is required for SAHF formation via regulation of HMGA2 expression SAHF formation leads to expression of SAHF-adjacent genes via 3D chromatin contacts
Collapse
|
23
|
Fiorillo L, Bianco S, Esposito A, Conte M, Sciarretta R, Musella F, Chiariello AM. A modern challenge of polymer physics: Novel ways to study, interpret, and reconstruct chromatin structure. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1454] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Luca Fiorillo
- Dipartimento di Fisica Università di Napoli Federico II, and INFN Napoli Complesso Universitario di Monte Sant'Angelo Naples Italy
| | - Simona Bianco
- Dipartimento di Fisica Università di Napoli Federico II, and INFN Napoli Complesso Universitario di Monte Sant'Angelo Naples Italy
| | - Andrea Esposito
- Dipartimento di Fisica Università di Napoli Federico II, and INFN Napoli Complesso Universitario di Monte Sant'Angelo Naples Italy
| | - Mattia Conte
- Dipartimento di Fisica Università di Napoli Federico II, and INFN Napoli Complesso Universitario di Monte Sant'Angelo Naples Italy
| | - Renato Sciarretta
- Dipartimento di Fisica Università di Napoli Federico II, and INFN Napoli Complesso Universitario di Monte Sant'Angelo Naples Italy
| | - Francesco Musella
- Dipartimento di Fisica Università di Napoli Federico II, and INFN Napoli Complesso Universitario di Monte Sant'Angelo Naples Italy
| | - Andrea M. Chiariello
- Dipartimento di Fisica Università di Napoli Federico II, and INFN Napoli Complesso Universitario di Monte Sant'Angelo Naples Italy
| |
Collapse
|
24
|
Briand N, Collas P. Laminopathy-causing lamin A mutations reconfigure lamina-associated domains and local spatial chromatin conformation. Nucleus 2019. [PMID: 29517398 PMCID: PMC5973257 DOI: 10.1080/19491034.2018.1449498] [Citation(s) in RCA: 42] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
The nuclear lamina contributes to the regulation of gene expression and to chromatin organization. Mutations in A-type nuclear lamins cause laminopathies, some of which are associated with a loss of heterochromatin at the nuclear periphery. Until recently however, little if any information has been provided on where and how lamin A interacts with the genome and on how disease-causing lamin A mutations may rearrange genome conformation. Here, we review aspects of nuclear lamin association with the genome. We highlight recent evidence of reorganization of lamin A-chromatin interactions in cellular models of laminopathies, and implications on the 3-dimensional rearrangement of chromatin in these models, including patient cells. We discuss how a hot-spot lipodystrophic lamin A mutation alters chromatin conformation and epigenetic patterns at an anti-adipogenic locus, and conclude with remarks on links between lamin A, Polycomb and the pathophysiology of laminopathies. The recent findings presented here collectively argue towards a deregulation of large-scale and local spatial genome organization by a subset of lamin A mutations causing laminopathies.
Collapse
Affiliation(s)
- Nolwenn Briand
- a Department of Molecular Medicine , Institute of Basic Medical Sciences, Faculty of Medicine, University of Oslo , Oslo , Norway
| | - Philippe Collas
- a Department of Molecular Medicine , Institute of Basic Medical Sciences, Faculty of Medicine, University of Oslo , Oslo , Norway.,b Norwegian Center for Stem Cell Research, Department of Immunology and Transfusion Medicine , Oslo University Hospital , Oslo , Norway
| |
Collapse
|
25
|
Human pancreatic islet three-dimensional chromatin architecture provides insights into the genetics of type 2 diabetes. Nat Genet 2019; 51:1137-1148. [PMID: 31253982 DOI: 10.1038/s41588-019-0457-0] [Citation(s) in RCA: 165] [Impact Index Per Article: 33.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2018] [Accepted: 05/29/2019] [Indexed: 01/07/2023]
Abstract
Genetic studies promise to provide insight into the molecular mechanisms underlying type 2 diabetes (T2D). Variants associated with T2D are often located in tissue-specific enhancer clusters or super-enhancers. So far, such domains have been defined through clustering of enhancers in linear genome maps rather than in three-dimensional (3D) space. Furthermore, their target genes are often unknown. We have created promoter capture Hi-C maps in human pancreatic islets. This linked diabetes-associated enhancers to their target genes, often located hundreds of kilobases away. It also revealed >1,300 groups of islet enhancers, super-enhancers and active promoters that form 3D hubs, some of which show coordinated glucose-dependent activity. We demonstrate that genetic variation in hubs impacts insulin secretion heritability, and show that hub annotations can be used for polygenic scores that predict T2D risk driven by islet regulatory variants. Human islet 3D chromatin architecture, therefore, provides a framework for interpretation of T2D genome-wide association study (GWAS) signals.
Collapse
|
26
|
Collas P, Liyakat Ali TM, Brunet A, Germier T. Finding Friends in the Crowd: Three-Dimensional Cliques of Topological Genomic Domains. Front Genet 2019; 10:602. [PMID: 31275364 PMCID: PMC6593077 DOI: 10.3389/fgene.2019.00602] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2019] [Accepted: 06/05/2019] [Indexed: 12/31/2022] Open
Abstract
The mammalian genome is intricately folded in a three-dimensional topology believed to be important for the orchestration of gene expression regulating development, differentiation and tissue homeostasis. Important features of spatial genome conformation in the nucleus are promoter-enhancer contacts regulating gene expression within topologically-associated domains (TADs), short- and long-range interactions between TADs and associations of chromatin with nucleoli and nuclear speckles. In addition, anchoring of chromosomes to the nuclear lamina via lamina-associated domains (LADs) at the nuclear periphery is a key regulator of the radial distribution of chromatin. To what extent TADs and LADs act in concert as genomic organizers to shape the three-dimensional topology of chromatin has long remained unknown. A new study addressing this key question provides evidence of (i) preferred long-range associations between TADs forming TAD “cliques” which organize large heterochromatin domains, and (ii) stabilization of TAD cliques by LADs at the nuclear periphery after induction of terminal differentiation. Here, we review these findings, address the issue of whether TAD cliques exist in single cells and discuss the extent of cell-to-cell heterogeneity in higher-order chromatin conformation. The recent observations provide a first appreciation of changes in 4-dimensional higher-order genome topologies during differentiation.
Collapse
Affiliation(s)
- Philippe Collas
- Department of Molecular Medicine, Faculty of Medicine, Institute of Basic Medical Sciences, University of Oslo, Oslo, Norway.,Department of Immunology and Transfusion Medicine, Oslo University Hospital, Oslo, Norway
| | - Tharvesh M Liyakat Ali
- Department of Molecular Medicine, Faculty of Medicine, Institute of Basic Medical Sciences, University of Oslo, Oslo, Norway
| | - Annaël Brunet
- Department of Molecular Medicine, Faculty of Medicine, Institute of Basic Medical Sciences, University of Oslo, Oslo, Norway
| | - Thomas Germier
- Department of Molecular Medicine, Faculty of Medicine, Institute of Basic Medical Sciences, University of Oslo, Oslo, Norway
| |
Collapse
|
27
|
Rosenthal M, Bryner D, Huffer F, Evans S, Srivastava A, Neretti N. Bayesian Estimation of Three-Dimensional Chromosomal Structure from Single-Cell Hi-C Data. J Comput Biol 2019; 26:1191-1202. [PMID: 31211598 PMCID: PMC6856950 DOI: 10.1089/cmb.2019.0100] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
The problem of three-dimensional (3D) chromosome structure inference from Hi-C data sets is important and challenging. While bulk Hi-C data sets contain contact information derived from millions of cells and can capture major structural features shared by the majority of cells in the sample, they do not provide information about local variability between cells. Single-cell Hi-C can overcome this problem, but contact matrices are generally very sparse, making structural inference more problematic. We have developed a Bayesian multiscale approach, named Structural Inference via Multiscale Bayesian Approach, to infer 3D structures of chromosomes from single-cell Hi-C while including the bulk Hi-C data and some regularization terms as a prior. We study the landscape of solutions for each single-cell Hi-C data set as a function of prior strength and demonstrate clustering of solutions using data from the same cell.
Collapse
Affiliation(s)
- Michael Rosenthal
- Science and Technology Department, Naval Surface Warfare Center, Panama City Division, Panama City, Florida
| | - Darshan Bryner
- Science and Technology Department, Naval Surface Warfare Center, Panama City Division, Panama City, Florida
| | - Fred Huffer
- Department of Statistics, Florida State University, Tallahassee, Florida
| | - Shane Evans
- Center for Computational Molecular Biology, Brown University, Providence, Rhode Island
| | - Anuj Srivastava
- Department of Statistics, Florida State University, Tallahassee, Florida
| | - Nicola Neretti
- Center for Computational Molecular Biology, Brown University, Providence, Rhode Island.,Department of Molecular Biology, Cell Biology, and Biochemistry, Brown University, Providence, Rhode Island
| |
Collapse
|
28
|
Fu S, Zhang L, Lv J, Zhu B, Wang W, Wang X. Two main stream methods analysis and visual 3D genome architecture. Semin Cell Dev Biol 2019; 90:43-53. [DOI: 10.1016/j.semcdb.2018.07.011] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2018] [Accepted: 07/10/2018] [Indexed: 01/07/2023]
|
29
|
Oluwadare O, Highsmith M, Cheng J. An Overview of Methods for Reconstructing 3-D Chromosome and Genome Structures from Hi-C Data. Biol Proced Online 2019; 21:7. [PMID: 31049033 PMCID: PMC6482566 DOI: 10.1186/s12575-019-0094-0] [Citation(s) in RCA: 74] [Impact Index Per Article: 14.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2019] [Accepted: 04/01/2019] [Indexed: 01/08/2023] Open
Abstract
Over the past decade, methods for predicting three-dimensional (3-D) chromosome and genome structures have proliferated. This has been primarily due to the development of high-throughput, next-generation chromosome conformation capture (3C) technologies, which have provided next-generation sequencing data about chromosome conformations in order to map the 3-D genome structure. The introduction of the Hi-C technique-a variant of the 3C method-has allowed researchers to extract the interaction frequency (IF) for all loci of a genome at high-throughput and at a genome-wide scale. In this review we describe, categorize, and compare the various methods developed to map chromosome and genome structures from 3C data-particularly Hi-C data. We summarize the improvements introduced by these methods, describe the approach used for method evaluation, and discuss how these advancements shape the future of genome structure construction.
Collapse
Affiliation(s)
- Oluwatosin Oluwadare
- Department of Electrical Engineering and Computer Science, University of Missouri, Columbia, MO 65211 USA
| | - Max Highsmith
- Department of Electrical Engineering and Computer Science, University of Missouri, Columbia, MO 65211 USA
| | - Jianlin Cheng
- Department of Electrical Engineering and Computer Science, University of Missouri, Columbia, MO 65211 USA
- Informatics Institute, University of Missouri, Columbia, MO 65211 USA
| |
Collapse
|
30
|
Hierarchical Reconstruction of High-Resolution 3D Models of Large Chromosomes. Sci Rep 2019; 9:4971. [PMID: 30899036 PMCID: PMC6428844 DOI: 10.1038/s41598-019-41369-w] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2018] [Accepted: 03/07/2019] [Indexed: 11/08/2022] Open
Abstract
Eukaryotic chromosomes are often composed of components organized into multiple scales, such as nucleosomes, chromatin fibers, topologically associated domains (TAD), chromosome compartments, and chromosome territories. Therefore, reconstructing detailed 3D models of chromosomes in high resolution is useful for advancing genome research. However, the task of constructing quality high-resolution 3D models is still challenging with existing methods. Hence, we designed a hierarchical algorithm, called Hierarchical3DGenome, to reconstruct 3D chromosome models at high resolution (<=5 Kilobase (KB)). The algorithm first reconstructs high-resolution 3D models at TAD level. The TAD models are then assembled to form complete high-resolution chromosomal models. The assembly of TAD models is guided by a complete low-resolution chromosome model. The algorithm is successfully used to reconstruct 3D chromosome models at 5 KB resolution for the human B-cell (GM12878). These high-resolution models satisfy Hi-C chromosomal contacts well and are consistent with models built at lower (i.e. 1 MB) resolution, and with the data of fluorescent in situ hybridization experiments. The Java source code of Hierarchical3DGenome and its user manual are available here https://github.com/BDM-Lab/Hierarchical3DGenome .
Collapse
|
31
|
Caudai C, Salerno E, Zoppe M, Tonazzini A. Estimation of the Spatial Chromatin Structure Based on a Multiresolution Bead-Chain Model. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2019; 16:550-559. [PMID: 29994172 DOI: 10.1109/tcbb.2018.2791439] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
We present a method to infer 3D chromatin configurations from Chromosome Conformation Capture data. Quite a few methods have been proposed to estimate the structure of the nuclear dna in homogeneous populations of cells from this kind of data. Many of them transform contact frequencies into euclidean distances between pairs of chromatin fragments, and then reconstruct the structure by solving a distance-to-geometry problem. To avoid inconsistencies, our method is based on a score function that does not require any frequency-to-distance translation. We propose a multiscale chromatin model where the chromatin fiber is suitably partitioned at each scale. The partial structures are estimated independently, and connected to rebuild the whole fiber. Our score function consists of a data-fit part and a penalty part, balanced automatically at each scale and each subchain. The penalty part enforces soft geometric constraints. As many different structures can fit the data, our sampling strategy produces a set of solutions with similar scores. The procedure contains a few parameters, independent of both the scale and the genomic segment treated. The partition of the fiber, along with intrinsically parallel parts, make this method computationally efficient. Results from human genome data support the biological plausibility of our solutions.
Collapse
|
32
|
Nir G, Farabella I, Pérez Estrada C, Ebeling CG, Beliveau BJ, Sasaki HM, Lee SD, Nguyen SC, McCole RB, Chattoraj S, Erceg J, AlHaj Abed J, Martins NMC, Nguyen HQ, Hannan MA, Russell S, Durand NC, Rao SSP, Kishi JY, Soler-Vila P, Di Pierro M, Onuchic JN, Callahan SP, Schreiner JM, Stuckey JA, Yin P, Aiden EL, Marti-Renom MA, Wu CT. Walking along chromosomes with super-resolution imaging, contact maps, and integrative modeling. PLoS Genet 2018; 14:e1007872. [PMID: 30586358 PMCID: PMC6324821 DOI: 10.1371/journal.pgen.1007872] [Citation(s) in RCA: 164] [Impact Index Per Article: 27.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2018] [Revised: 01/08/2019] [Accepted: 12/04/2018] [Indexed: 12/13/2022] Open
Abstract
Chromosome organization is crucial for genome function. Here, we present a method for visualizing chromosomal DNA at super-resolution and then integrating Hi-C data to produce three-dimensional models of chromosome organization. Using the super-resolution microscopy methods of OligoSTORM and OligoDNA-PAINT, we trace 8 megabases of human chromosome 19, visualizing structures ranging in size from a few kilobases to over a megabase. Focusing on chromosomal regions that contribute to compartments, we discover distinct structures that, in spite of considerable variability, can predict whether such regions correspond to active (A-type) or inactive (B-type) compartments. Imaging through the depths of entire nuclei, we capture pairs of homologous regions in diploid cells, obtaining evidence that maternal and paternal homologous regions can be differentially organized. Finally, using restraint-based modeling to integrate imaging and Hi-C data, we implement a method-integrative modeling of genomic regions (IMGR)-to increase the genomic resolution of our traces to 10 kb.
Collapse
MESH Headings
- Cells, Cultured
- Chromosome Painting/methods
- Chromosome Structures/chemistry
- Chromosome Structures/genetics
- Chromosome Structures/ultrastructure
- Chromosome Walking/methods
- Chromosomes, Human, Pair 19/chemistry
- Chromosomes, Human, Pair 19/genetics
- Chromosomes, Human, Pair 19/ultrastructure
- Female
- Fluorescent Dyes
- Humans
- Imaging, Three-Dimensional
- In Situ Hybridization, Fluorescence/methods
- Male
- Models, Genetic
- Oligonucleotide Probes
- Pedigree
Collapse
Affiliation(s)
- Guy Nir
- Department of Genetics, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Irene Farabella
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Barcelona, Spain
| | - Cynthia Pérez Estrada
- Department of Genetics, Harvard Medical School, Boston, Massachusetts, United States of America
- Center for Genome Architecture, Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States of America
- Center for Theoretical Biological Physics, Rice University, Houston, Texas, United States of America
| | - Carl G. Ebeling
- Bruker Nano Inc., Salt Lake City, Utah, United States of America
| | - Brian J. Beliveau
- Department of Genetics, Harvard Medical School, Boston, Massachusetts, United States of America
- Wyss Institute for Biologically Inspired Engineering, Harvard University, Boston, Massachusetts, United States of America
- Department of Systems Biology, Harvard Medical School, Boston, Massachusetts, United States of America
- Department of Genome Sciences, University of Washington, Seattle, Washington, United States of America
| | - Hiroshi M. Sasaki
- Wyss Institute for Biologically Inspired Engineering, Harvard University, Boston, Massachusetts, United States of America
- Department of Systems Biology, Harvard Medical School, Boston, Massachusetts, United States of America
| | - S. Dean Lee
- Department of Genetics, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Son C. Nguyen
- Department of Genetics, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Ruth B. McCole
- Department of Genetics, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Shyamtanu Chattoraj
- Department of Genetics, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Jelena Erceg
- Department of Genetics, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Jumana AlHaj Abed
- Department of Genetics, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Nuno M. C. Martins
- Department of Genetics, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Huy Q. Nguyen
- Department of Genetics, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Mohammed A. Hannan
- Department of Genetics, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Sheikh Russell
- Center for Genome Architecture, Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States of America
| | - Neva C. Durand
- Center for Genome Architecture, Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States of America
- Broad Institute of the Massachusetts Institute of Technology and Harvard University, Cambridge, Massachusetts, United States of America
| | - Suhas S. P. Rao
- Center for Genome Architecture, Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States of America
- Center for Theoretical Biological Physics, Rice University, Houston, Texas, United States of America
- Department of Structural Biology, Stanford University School of Medicine, Stanford, California, United States of America
| | - Jocelyn Y. Kishi
- Wyss Institute for Biologically Inspired Engineering, Harvard University, Boston, Massachusetts, United States of America
- Department of Systems Biology, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Paula Soler-Vila
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Barcelona, Spain
| | - Michele Di Pierro
- Center for Theoretical Biological Physics, Rice University, Houston, Texas, United States of America
| | - José N. Onuchic
- Center for Theoretical Biological Physics, Rice University, Houston, Texas, United States of America
| | | | | | - Jeff A. Stuckey
- Bruker Nano Inc., Middleton, Wisconsin, United States of America
| | - Peng Yin
- Wyss Institute for Biologically Inspired Engineering, Harvard University, Boston, Massachusetts, United States of America
- Department of Systems Biology, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Erez Lieberman Aiden
- Center for Genome Architecture, Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States of America
- Center for Theoretical Biological Physics, Rice University, Houston, Texas, United States of America
- Broad Institute of the Massachusetts Institute of Technology and Harvard University, Cambridge, Massachusetts, United States of America
- Departments of Computer Science and Computational and Applied Mathematics, Rice University, Houston, Texas, United States of America
| | - Marc A. Marti-Renom
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Barcelona, Spain
- Gene Regulation, Stem Cells and Cancer Program, Centre for Genomic Regulation (CRG), Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
- ICREA, Barcelona, Spain
| | - C.-ting Wu
- Department of Genetics, Harvard Medical School, Boston, Massachusetts, United States of America
- Wyss Institute for Biologically Inspired Engineering, Harvard University, Boston, Massachusetts, United States of America
| |
Collapse
|
33
|
Le Treut G, Képès F, Orland H. A Polymer Model for the Quantitative Reconstruction of Chromosome Architecture from HiC and GAM Data. Biophys J 2018; 115:2286-2294. [PMID: 30527448 DOI: 10.1016/j.bpj.2018.10.032] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2018] [Revised: 10/03/2018] [Accepted: 10/26/2018] [Indexed: 01/03/2023] Open
Abstract
It is widely believed that the folding of the chromosome in the nucleus has a major effect on genetic expression. For example, coregulated genes in several species have been shown to colocalize in space despite being far away on the DNA sequence. In this manuscript, we present a new, to our knowledge, method to model the three-dimensional structure of the chromosome in live cells based on DNA-DNA interactions measured in high-throughput chromosome conformation capture experiments and genome architecture mapping. Our approach incorporates a polymer model and directly uses the contact probabilities measured in high-throughput chromosome conformation capture experiments and genome architecture mapping experiments rather than estimates of average distances between genomic loci. Specifically, we model the chromosome as a Gaussian polymer with harmonic interactions and extract the coupling coefficients best reproducing the experimental contact probabilities. In contrast to existing methods, we give an exact expression of the contact probabilities at thermodynamic equilibrium. The Gaussian effective model reconstructed with our method reproduces experimental contacts with high accuracy. We also show how Brownian dynamics simulations of our reconstructed Gaussian effective model can be used to study chromatin organization and possibly give some clue about its dynamics.
Collapse
Affiliation(s)
- Guillaume Le Treut
- Department of Physics, University of California San Diego, La Jolla, California.
| | - François Képès
- institute of Systems and Synthetic Biology, Genopole, CNRS, UEVE, Université Paris-Saclay, Évry, France
| | - Henri Orland
- Institut de Physique Théorique, CEA, CNRS-URA 2306, Gif-sur-Yvette, France; Beijing Computational Science Research Center, Beijing, China
| |
Collapse
|
34
|
Diament A, Tuller T. Modeling three-dimensional genomic organization in evolution and pathogenesis. Semin Cell Dev Biol 2018; 90:78-93. [PMID: 30030143 DOI: 10.1016/j.semcdb.2018.07.008] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2018] [Accepted: 07/08/2018] [Indexed: 12/17/2022]
Abstract
The regulation of gene expression is mediated via the complex three-dimensional (3D) conformation of the genetic material and its interactions with various intracellular factors. Various experimental and computational approaches have been developed in recent years for understating the relation between the 3D conformation of the genome and the phenotypes of cells in normal condition and diseases. In this review, we will discuss novel approaches for analyzing and modeling the 3D genomic conformation, focusing on deciphering disease-causing mutations that affect gene expression. We conclude that as this is a very challenging mission, an important direction should involve the comparative analysis of various 3D models from various organisms or cells.
Collapse
Affiliation(s)
- Alon Diament
- Dept. of Biomedical Engineering, Tel Aviv University, Tel Aviv 6997801, Israel
| | - Tamir Tuller
- Dept. of Biomedical Engineering, Tel Aviv University, Tel Aviv 6997801, Israel; The Sagol School of Neuroscience, Tel-Aviv University, Tel Aviv 6997801, Israel.
| |
Collapse
|
35
|
Waldispühl J, Zhang E, Butyaev A, Nazarova E, Cyr Y. Storage, visualization, and navigation of 3D genomics data. Methods 2018; 142:74-80. [PMID: 29792917 DOI: 10.1016/j.ymeth.2018.05.008] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2017] [Revised: 05/07/2018] [Accepted: 05/09/2018] [Indexed: 01/27/2023] Open
Abstract
The field of 3D genomics grew at increasing rates in the last decade. The volume and complexity of 2D and 3D data produced is progressively outpacing the capacities of the technology previously used for distributing genome sequences. The emergence of new technologies provides also novel opportunities for the development of innovative approaches. In this paper, we review the state-of-the-art computing technology, as well as the solutions adopted by the platforms currently available.
Collapse
Affiliation(s)
| | - Eric Zhang
- School of Computer Science, McGill University, Montréal, Canada
| | | | - Elena Nazarova
- School of Computer Science, McGill University, Montréal, Canada
| | - Yan Cyr
- Beam Me Up Labs, Montréal, Canada
| |
Collapse
|
36
|
Caudai C, Salerno E, Zoppe M, Merelli I, Tonazzini A. ChromStruct 4: A Python Code to Estimate the Chromatin Structure from Hi-C Data. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2018:1-1. [PMID: 29993555 DOI: 10.1109/tcbb.2018.2838669] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
A method and a stand-alone Python(TM) code to estimate the 3D chromatin structure from chromosome conformation capture data are presented. The method is based on a multiresolution, modified-bead-chain chromatin model, evolved through quaternion operators in a Monte Carlo sampling. The solution space to be sampled is generated by a score function with a data-fit part and a constraint part where the available prior knowledge is implicitly coded. The final solution is a set of 3D configurations that are compatible with both the data and the prior knowledge. The iterative code, provided here as additional material, is equipped with a graphical user interface and stores its results in standard-format files for 3D visualization. We describe the mathematical-computational aspects of the method and explain the details of the code. Some experimental results are reported, with a demonstration of their fit to the data.
Collapse
|
37
|
Abstract
Motivation Recent experiments have provided Hi-C data at resolution as high as 1 kbp. However, 3D structural inference from high-resolution Hi-C datasets is often computationally unfeasible using existing methods. Results We have developed miniMDS, an approximation of multidimensional scaling (MDS) that partitions a Hi-C dataset, performs high-resolution MDS separately on each partition, and then reassembles the partitions using low-resolution MDS. miniMDS is faster, more accurate, and uses less memory than existing methods for inferring the human genome at high resolution (10 kbp). Availability and implementation A Python implementation of miniMDS is available on GitHub: https://github.com/seqcode/miniMDS. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Lila Rieber
- Department of Biochemistry and Molecular Biology and Center for Eukaryotic Gene Regulation, The Pennsylvania State University, University Park, PA, USA
| | - Shaun Mahony
- Department of Biochemistry and Molecular Biology and Center for Eukaryotic Gene Regulation, The Pennsylvania State University, University Park, PA, USA
| |
Collapse
|
38
|
Abstract
Chrom3D is a computational platform for 3D genome modeling that simulates the spatial positioning of chromosome domains relative to each other and relative to the nuclear periphery. In Chrom3D, chromosomes are modeled as chains of contiguous beads, in which each bead represents a genomic domain. In this protocol, a bead represents a topologically associated domain (TAD) mapped from ensemble Hi-C data. Chrom3D takes as input data significant pairwise TAD-TAD interactions determined from a Hi-C contact matrix, and TAD interactions with the nuclear periphery, determined by ChIP-sequencing of nuclear lamins to define lamina-associated domains (LADs). Chrom3D is based on Monte Carlo simulations initiated from a starting random bead configuration. During the optimization process, TAD-TAD interactions constrain bead positions relative to each other, whereas LAD information constrains the corresponding bead toward the nuclear periphery. Optimization can be repeated many times to generate an ensemble of 3D genome models. Analyses of the models enable estimations of the radial positioning of genomic sites in the nucleus across cells in a population. Chrom3D provides opportunities to reveal spatial relationships between TADs and LADs. More generally, predictions from Chrom3D models can be experimentally tested in the laboratory. We describe the entire Chrom3D protocol for modeling a 3D diploid human genome, from the creation of input files to the final rendering of 3D genome structures. The procedure takes ∼18 h. Chrom3D is freely available on GitHub.
Collapse
|
39
|
Abstract
Chromosome conformation capture technologies such as Hi-C are widely used to investigate the spatial organization of genomes. Because genome structures can vary considerably between individual cells of a population, interpreting ensemble-averaged Hi-C data can be challenging, in particular for long-range and interchromosomal interactions. We pioneered a probabilistic approach for the generation of a population of distinct diploid 3D genome structures consistent with all the chromatin-chromatin interaction probabilities from Hi-C experiments. Each structure in the population is a physical model of the genome in 3D. Analysis of these models yields new insights into the causes and the functional properties of the genome's organization in space and time. We provide a user-friendly software package, called PGS, which runs on local machines (for practice runs) and high-performance computing platforms. PGS takes a genome-wide Hi-C contact frequency matrix, along with information about genome segmentation, and produces an ensemble of 3D genome structures entirely consistent with the input. The software automatically generates an analysis report, and provides tools to extract and analyze the 3D coordinates of specific domains. Basic Linux command-line knowledge is sufficient for using this software. A typical running time of the pipeline is ∼3 d with 300 cores on a computer cluster to generate a population of 1,000 diploid genome structures at topological-associated domain (TAD)-level resolution.
Collapse
|
40
|
Oluwadare O, Zhang Y, Cheng J. A maximum likelihood algorithm for reconstructing 3D structures of human chromosomes from chromosomal contact data. BMC Genomics 2018; 19:161. [PMID: 29471801 PMCID: PMC5824572 DOI: 10.1186/s12864-018-4546-8] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2017] [Accepted: 02/13/2018] [Indexed: 01/07/2023] Open
Abstract
Background The development of chromosomal conformation capture techniques, particularly, the Hi-C technique, has made the analysis and study of the spatial conformation of a genome an important topic in bioinformatics and computational biology. Aided by high-throughput next generation sequencing techniques, the Hi-C technique can generate genome-wide, large-scale intra- and inter-chromosomal interaction data capable of describing in details the spatial interactions within a genome. These data can be used to reconstruct 3D structures of chromosomes that can be used to study DNA replication, gene regulation, genome interaction, genome folding, and genome function. Results Here, we introduce a maximum likelihood algorithm called 3DMax to construct the 3D structure of a chromosome from Hi-C data. 3DMax employs a maximum likelihood approach to infer the 3D structures of a chromosome, while automatically re-estimating the conversion factor (α) for converting Interaction Frequency (IF) to distance. Our results show that the models generated by 3DMax from a simulated Hi-C dataset match the true models better than most of the existing methods. 3DMax is more robust to structural variability and noise. Compared on a real Hi-C dataset, 3DMax constructs chromosomal models that fit the data better than most methods, and it is faster than all other methods. The models reconstructed by 3DMax were consistent with fluorescent in situ hybridization (FISH) experiments and existing knowledge about the organization of human chromosomes, such as chromosome compartmentalization. Conclusions 3DMax is an effective approach to reconstructing 3D chromosomal models. The results, and the models generated for the simulated and real Hi-C datasets are available here: http://sysbio.rnet.missouri.edu/bdm_download/3DMax/. The source code is available here: https://github.com/BDM-Lab/3DMax. A short video demonstrating how to use 3DMax can be found here: https://youtu.be/ehQUFWoHwfo.
Collapse
Affiliation(s)
- Oluwatosin Oluwadare
- Electrical Engineering & Computer Science Department, University of Missouri, Columbia, MO, 65211, USA
| | - Yuxiang Zhang
- Electrical Engineering & Computer Science Department, University of Missouri, Columbia, MO, 65211, USA
| | - Jianlin Cheng
- Electrical Engineering & Computer Science Department, University of Missouri, Columbia, MO, 65211, USA. .,Informatics Institute, University of Missouri, Columbia, MO, 65211, USA.
| |
Collapse
|
41
|
Trieu T, Cheng J. 3D genome structure modeling by Lorentzian objective function. Nucleic Acids Res 2017; 45:1049-1058. [PMID: 28180292 PMCID: PMC5430849 DOI: 10.1093/nar/gkw1155] [Citation(s) in RCA: 38] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2016] [Revised: 11/01/2016] [Accepted: 11/04/2016] [Indexed: 12/19/2022] Open
Abstract
The 3D structure of the genome plays a vital role in biological processes such as gene interaction, gene regulation, DNA replication and genome methylation. Advanced chromosomal conformation capture techniques, such as Hi-C and tethered conformation capture, can generate chromosomal contact data that can be used to computationally reconstruct 3D structures of the genome. We developed a novel restraint-based method that is capable of reconstructing 3D genome structures utilizing both intra-and inter-chromosomal contact data. Our method was robust to noise and performed well in comparison with a panel of existing methods on a controlled simulated data set. On a real Hi-C data set of the human genome, our method produced chromosome and genome structures that are consistent with 3D FISH data and known knowledge about the human chromosome and genome, such as, chromosome territories and the cluster of small chromosomes in the nucleus center with the exception of the chromosome 18. The tool and experimental data are available at https://missouri.box.com/v/LorDG.
Collapse
Affiliation(s)
- Tuan Trieu
- Computer Science Department, University of Missouri-Columbia, MO, USA
| | - Jianlin Cheng
- Computer Science Department, University of Missouri-Columbia, MO, USA.,Informatics Institute, University of Missouri-Columbia, MO, USA
| |
Collapse
|
42
|
Integrative modelling of cellular assemblies. Curr Opin Struct Biol 2017; 46:102-109. [PMID: 28735107 DOI: 10.1016/j.sbi.2017.07.001] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2017] [Revised: 07/01/2017] [Accepted: 07/04/2017] [Indexed: 02/06/2023]
Abstract
A wide variety of experimental techniques can be used for understanding the precise molecular mechanisms underlying the activities of cellular assemblies. The inherent limitations of a single experimental technique often requires integration of data from complementary approaches to gain sufficient insights into the assembly structure and function. Here, we review popular computational approaches for integrative modelling of cellular assemblies, including protein complexes and genomic assemblies. We provide recent examples of integrative models generated for such assemblies by different experimental techniques, especially including data from 3D electron microscopy (3D-EM) and chromosome conformation capture experiments, respectively. We highlight general concepts in integrative modelling and discuss the need for careful formulation and merging of different types of information.
Collapse
|
43
|
Automatic analysis and 3D-modelling of Hi-C data using TADbit reveals structural features of the fly chromatin colors. PLoS Comput Biol 2017; 13:e1005665. [PMID: 28723903 PMCID: PMC5540598 DOI: 10.1371/journal.pcbi.1005665] [Citation(s) in RCA: 174] [Impact Index Per Article: 24.9] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2017] [Revised: 08/02/2017] [Accepted: 07/03/2017] [Indexed: 11/19/2022] Open
Abstract
The sequence of a genome is insufficient to understand all genomic processes carried out in the cell nucleus. To achieve this, the knowledge of its three-dimensional architecture is necessary. Advances in genomic technologies and the development of new analytical methods, such as Chromosome Conformation Capture (3C) and its derivatives, provide unprecedented insights in the spatial organization of genomes. Here we present TADbit, a computational framework to analyze and model the chromatin fiber in three dimensions. Our package takes as input the sequencing reads of 3C-based experiments and performs the following main tasks: (i) pre-process the reads, (ii) map the reads to a reference genome, (iii) filter and normalize the interaction data, (iv) analyze the resulting interaction matrices, (v) build 3D models of selected genomic domains, and (vi) analyze the resulting models to characterize their structural properties. To illustrate the use of TADbit, we automatically modeled 50 genomic domains from the fly genome revealing differential structural features of the previously defined chromatin colors, establishing a link between the conformation of the genome and the local chromatin composition. TADbit provides three-dimensional models built from 3C-based experiments, which are ready for visualization and for characterizing their relation to gene expression and epigenetic states. TADbit is an open-source Python library available for download from https://github.com/3DGenomes/tadbit.
Collapse
|
44
|
Le Dily F, Serra F, Marti-Renom MA. 3D modeling of chromatin structure: is there a way to integrate and reconcile single cell and population experimental data? WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1308] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- François Le Dily
- Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology; Dr. Aiguader 88; Barcelona Spain
- Universitat Pompeu Fabra (UPF); Barcelona Spain
| | - François Serra
- Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology; Dr. Aiguader 88; Barcelona Spain
- Universitat Pompeu Fabra (UPF); Barcelona Spain
- Structural Genomic Group, CNAG-CRG, Centre for Genomic Regulation (CRG); The Barcelona Institute of Science and Technology, Baldiri Reixac 4; Barcelona Spain
| | - Marc A. Marti-Renom
- Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology; Dr. Aiguader 88; Barcelona Spain
- Universitat Pompeu Fabra (UPF); Barcelona Spain
- Structural Genomic Group, CNAG-CRG, Centre for Genomic Regulation (CRG); The Barcelona Institute of Science and Technology, Baldiri Reixac 4; Barcelona Spain
- Institució Catalana de Recerca i Estudis Avançats (ICREA), Pg. Lluis Companys 23; Barcelona Spain
| |
Collapse
|
45
|
Paulsen J, Sekelja M, Oldenburg AR, Barateau A, Briand N, Delbarre E, Shah A, Sørensen AL, Vigouroux C, Buendia B, Collas P. Chrom3D: three-dimensional genome modeling from Hi-C and nuclear lamin-genome contacts. Genome Biol 2017; 18:21. [PMID: 28137286 PMCID: PMC5278575 DOI: 10.1186/s13059-016-1146-2] [Citation(s) in RCA: 134] [Impact Index Per Article: 19.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2016] [Accepted: 12/23/2016] [Indexed: 01/09/2023] Open
Abstract
Current three-dimensional (3D) genome modeling platforms are limited by their inability to account for radial placement of loci in the nucleus. We present Chrom3D, a user-friendly whole-genome 3D computational modeling framework that simulates positions of topologically-associated domains (TADs) relative to each other and to the nuclear periphery. Chrom3D integrates chromosome conformation capture (Hi-C) and lamin-associated domain (LAD) datasets to generate structure ensembles that recapitulate radial distributions of TADs detected in single cells. Chrom3D reveals unexpected spatial features of LAD regulation in cells from patients with a laminopathy-causing lamin mutation. Chrom3D is freely available on github.
Collapse
Affiliation(s)
- Jonas Paulsen
- Department of Molecular Medicine, Institute of Basic Medical Sciences, Faculty of Medicine, University of Oslo, Oslo, Norway
| | - Monika Sekelja
- Department of Molecular Medicine, Institute of Basic Medical Sciences, Faculty of Medicine, University of Oslo, Oslo, Norway
| | - Anja R Oldenburg
- Department of Molecular Medicine, Institute of Basic Medical Sciences, Faculty of Medicine, University of Oslo, Oslo, Norway
| | | | - Nolwenn Briand
- Department of Molecular Medicine, Institute of Basic Medical Sciences, Faculty of Medicine, University of Oslo, Oslo, Norway
| | - Erwan Delbarre
- Department of Molecular Medicine, Institute of Basic Medical Sciences, Faculty of Medicine, University of Oslo, Oslo, Norway
| | - Akshay Shah
- Department of Molecular Medicine, Institute of Basic Medical Sciences, Faculty of Medicine, University of Oslo, Oslo, Norway
| | - Anita L Sørensen
- Department of Molecular Medicine, Institute of Basic Medical Sciences, Faculty of Medicine, University of Oslo, Oslo, Norway
| | - Corinne Vigouroux
- INSERM, UMR S938, Centre de Recherches Saint-Antoine, Paris, France.,UPMC Université Paris 6 UMR S938, Paris, France.,ICAN, Paris, France.,AP-HP Hôpital Tenon, Paris, France
| | | | - Philippe Collas
- Department of Molecular Medicine, Institute of Basic Medical Sciences, Faculty of Medicine, University of Oslo, Oslo, Norway. .,Norwegian Center for Stem Cell Research, Oslo University Hospital, Oslo, Norway.
| |
Collapse
|
46
|
Computational inference of physical spatial organization of eukaryotic genomes. QUANTITATIVE BIOLOGY 2016. [DOI: 10.1007/s40484-016-0082-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
47
|
Szałaj P, Tang Z, Michalski P, Pietal MJ, Luo OJ, Sadowski M, Li X, Radew K, Ruan Y, Plewczynski D. An integrated 3-Dimensional Genome Modeling Engine for data-driven simulation of spatial genome organization. Genome Res 2016; 26:1697-1709. [PMID: 27789526 PMCID: PMC5131821 DOI: 10.1101/gr.205062.116] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2016] [Accepted: 10/20/2016] [Indexed: 02/03/2023]
Abstract
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.
Collapse
Affiliation(s)
- Przemysław Szałaj
- Centre of New Technologies, Warsaw University, 02-097 Warsaw, Poland.,Centre for Innovative Research, Medical University of Bialystok, 15-089 Białystok, Poland.,I-BioStat, Hasselt University, BE3590 Hasselt, Belgium
| | - Zhonghui Tang
- The Jackson Laboratory for Genomic Medicine, Farmington, Connecticut 06032, USA
| | - Paul Michalski
- The Jackson Laboratory for Genomic Medicine, Farmington, Connecticut 06032, USA
| | - Michal J Pietal
- Centre of New Technologies, Warsaw University, 02-097 Warsaw, Poland
| | - Oscar J Luo
- The Jackson Laboratory for Genomic Medicine, Farmington, Connecticut 06032, USA
| | - Michał Sadowski
- Centre of New Technologies, Warsaw University, 02-097 Warsaw, Poland
| | - Xingwang Li
- The Jackson Laboratory for Genomic Medicine, Farmington, Connecticut 06032, USA
| | - Kamen Radew
- Centre of New Technologies, Warsaw University, 02-097 Warsaw, Poland
| | - Yijun Ruan
- The Jackson Laboratory for Genomic Medicine, Farmington, Connecticut 06032, USA.,Department of Genetics and Genome Sciences, UConn Health, Farmington, Connecticut 06032, USA
| | - Dariusz Plewczynski
- Centre of New Technologies, Warsaw University, 02-097 Warsaw, Poland.,Centre for Innovative Research, Medical University of Bialystok, 15-089 Białystok, Poland.,Faculty of Pharmacy, Medical University of Warsaw, 02-097 Warsaw, Poland
| |
Collapse
|
48
|
Three-dimensional reconstruction of single-cell chromosome structure using recurrence plots. Sci Rep 2016; 6:34982. [PMID: 27725694 PMCID: PMC5057099 DOI: 10.1038/srep34982] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2016] [Accepted: 09/21/2016] [Indexed: 11/21/2022] Open
Abstract
Single-cell analysis of the three-dimensional (3D) chromosome structure can reveal cell-to-cell variability in genome activities. Here, we propose to apply recurrence plots, a mathematical method of nonlinear time series analysis, to reconstruct the 3D chromosome structure of a single cell based on information of chromosomal contacts from genome-wide chromosome conformation capture (Hi-C) data. This recurrence plot-based reconstruction (RPR) method enables rapid reconstruction of a unique structure in single cells, even from incomplete Hi-C information.
Collapse
|
49
|
Zhang B, Wolynes PG. Shape Transitions and Chiral Symmetry Breaking in the Energy Landscape of the Mitotic Chromosome. PHYSICAL REVIEW LETTERS 2016; 116:248101. [PMID: 27367409 DOI: 10.1103/physrevlett.116.248101] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/10/2015] [Indexed: 05/18/2023]
Abstract
We derive an unbiased information theoretic energy landscape for chromosomes at metaphase using a maximum entropy approach that accurately reproduces the details of the experimentally measured pairwise contact probabilities between genomic loci. Dynamical simulations using this landscape lead to cylindrical, helically twisted structures reflecting liquid crystalline order. These structures are similar to those arising from a generic ideal homogenized chromosome energy landscape. The helical twist can be either right or left handed so chiral symmetry is broken spontaneously. The ideal chromosome landscape when augmented by interactions like those leading to topologically associating domain formation in the interphase chromosome reproduces these behaviors. The phase diagram of this landscape shows that the helical fiber order and the cylindrical shape persist at temperatures above the onset of chiral symmetry breaking, which is limited by the topologically associating domain interaction strength.
Collapse
Affiliation(s)
- Bin Zhang
- Department of Chemistry, Rice University, Houston, Texas 77005, USA
- Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, USA
| | - Peter G Wolynes
- Department of Chemistry, Rice University, Houston, Texas 77005, USA
- Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, USA
- Department of Physics and Astronomy, Rice University, Houston, Texas 77005, USA
| |
Collapse
|
50
|
Szalaj P, Michalski PJ, Wróblewski P, Tang Z, Kadlof M, Mazzocco G, Ruan Y, Plewczynski D. 3D-GNOME: an integrated web service for structural modeling of the 3D genome. Nucleic Acids Res 2016; 44:W288-93. [PMID: 27185892 PMCID: PMC4987952 DOI: 10.1093/nar/gkw437] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2016] [Accepted: 05/07/2016] [Indexed: 11/13/2022] Open
Abstract
Recent advances in high-throughput chromosome conformation capture (3C) technology, such as Hi-C and ChIA-PET, have demonstrated the importance of 3D genome organization in development, cell differentiation and transcriptional regulation. There is now a widespread need for computational tools to generate and analyze 3D structural models from 3C data. Here we introduce our 3D GeNOme Modeling Engine (3D-GNOME), a web service which generates 3D structures from 3C data and provides tools to visually inspect and annotate the resulting structures, in addition to a variety of statistical plots and heatmaps which characterize the selected genomic region. Users submit a bedpe (paired-end BED format) file containing the locations and strengths of long range contact points, and 3D-GNOME simulates the structure and provides a convenient user interface for further analysis. Alternatively, a user may generate structures using published ChIA-PET data for the GM12878 cell line by simply specifying a genomic region of interest. 3D-GNOME is freely available at http://3dgnome.cent.uw.edu.pl/.
Collapse
Affiliation(s)
- Przemyslaw Szalaj
- Centre of New Technologies, University of Warsaw, 02-097 Warsaw, Poland Center for Bioinformatics and Data Analysis, Medical University of Bialystok, 15-089 Bialystok, Poland I-BioStat, Hasselt University, 3500 Hasselt, Belgium
| | - Paul J Michalski
- The Jackson Laboratory for Genomic Medicine, Farmington, CT 06032, USA
| | | | - Zhonghui Tang
- The Jackson Laboratory for Genomic Medicine, Farmington, CT 06032, USA
| | - Michal Kadlof
- Centre of New Technologies, University of Warsaw, 02-097 Warsaw, Poland
| | - Giovanni Mazzocco
- Centre of New Technologies, University of Warsaw, 02-097 Warsaw, Poland
| | - Yijun Ruan
- The Jackson Laboratory for Genomic Medicine, Farmington, CT 06032, USA Department of Genetics and Genome Sciences, UConn Health, Farmington, CT 06030-6403, USA
| | - Dariusz Plewczynski
- Centre of New Technologies, University of Warsaw, 02-097 Warsaw, Poland Center for Bioinformatics and Data Analysis, Medical University of Bialystok, 15-089 Bialystok, Poland Faculty of Pharmacy, Medical University of Warsaw, 02-097 Warsaw, Poland
| |
Collapse
|