scholarly journals Chromatin 3D structure reconstruction with consideration of adjacency relationship among genomic loci

2019 ◽  
Author(s):  
Fang-Zhen Li ◽  
Zhi-E Liu ◽  
Xiu-Yuan Li ◽  
Li-Mei Bu ◽  
Hong-Xia Bu ◽  
...  

AbstractChromatin 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. We present a method called ShNeigh, which combines the classical MDS technique with local dependence of neighboring loci modelled by a Gaussian formula, to infer the best 3D structure from noisy and incomplete contact frequency matrices. The results obtained on simulations and real Hi-C data showed, while keeping the high-speed nature of classical MDS, ShNeigh is more accurate and robust than existing methods, especially for sparse contact maps. A Matlab implementation of the proposed method is available at https://github.com/fangzhen-li/ShNeigh.Author summaryWe propose a new method to infer a consensus 3D genome structure from a Hi-C contact map. The novelty of our method is that it takes into accounts the adjacency of genomic loci along chromosomes. Specifically, the proposed method penalizes the optimization problem of the classical multidimensional scaling method with a smoothness constraint weighted by a function of the genomic distance between the pairs of genomic loci. We demonstrate this optimization problem can still be solved efficiently by a classical multidimensional scaling method. We then show that the method can recover stable structures in high noise settings. We also show that it can reconstruct similar structures from data obtained using different restriction enzymes.


2019 ◽  
Author(s):  
Oluwatosin Oluwadare ◽  
Max Highsmith ◽  
Jianlin Cheng

ABSTRACTAdvances in the study of chromosome conformation capture (3C) technologies, such as Hi-C technique - capable of capturing chromosomal interactions in a genome-wide scale - have led to the development of three-dimensional (3D) chromosome and genome structure reconstruction methods from Hi-C data. The 3D genome structure is important because it plays a role in a variety of important biological activities such as DNA replication, gene regulation, genome interaction, and gene expression. In recent years, numerous Hi-C datasets have been generated, and likewise, a number of genome structure construction algorithms have been developed. However, until now, there has been no freely available repository for 3D chromosome structures. In this work, we outline the construction of a novel Genome Structure Database (GSDB) to create a comprehensive repository that contains 3D structures for Hi-C datasets constructed by a variety of 3D structure reconstruction tools. GSDB contains over 50,000 structures constructed by 12 state-of-the-art chromosome and genome structure prediction methods for publicly used Hi-C datasets with varying resolution. The database is useful for the community to study the function of genome from a 3D perspective. GSDB is accessible at http://sysbio.rnet.missouri.edu/3dgenome/GSDB



2020 ◽  
Vol 21 (1) ◽  
Author(s):  
Priyojit Das ◽  
Tongye Shen ◽  
Rachel Patton McCord

Abstract Background The nonrandom radial organization of eukaryotic chromosome territories (CTs) inside the nucleus plays an important role in nuclear functional compartmentalization. Increasingly, chromosome conformation capture (Hi-C) based approaches are being used to characterize the genome structure of many cell types and conditions. Computational methods to extract 3D arrangements of CTs from this type of pairwise contact data will thus increase our ability to analyze CT organization in a wider variety of biological situations. Results A number of full-scale polymer models have successfully reconstructed the 3D structure of chromosome territories from Hi-C. To supplement such methods, we explore alternative, direct, and less computationally intensive approaches to capture radial CT organization from Hi-C data. We show that we can infer relative chromosome ordering using PCA on a thresholded inter-chromosomal contact matrix. We simulate an ensemble of possible CT arrangements using a force-directed network layout algorithm and propose an approach to integrate additional chromosome properties into our predictions. Our CT radial organization predictions have a high correlation with microscopy imaging data for various cell nucleus geometries (lymphoblastoid, skin fibroblast, and breast epithelial cells), and we can capture previously documented changes in senescent and progeria cells. Conclusions Our analysis approaches provide rapid and modular approaches to screen for alterations in CT organization across widely available Hi-C data. We demonstrate which stages of the approach can extract meaningful information, and also describe limitations of pairwise contacts alone to predict absolute 3D positions.



Genes ◽  
2021 ◽  
Vol 12 (11) ◽  
pp. 1671
Author(s):  
Ryo Ishibashi ◽  
Y-h. Taguchi

The positions of enhancers and promoters on genomic DNA remain poorly understood. Chromosomes cannot be observed during the cell division cycle because the genome forms a chromatin structure and spreads within the nucleus. However, high-throughput chromosome conformation capture (Hi-C) measures the physical interactions of genomes. In previous studies, DNA extrusion loops were directly derived from Hi-C heat maps. Multidimensional Scaling (MDS) is used in this assessment to more precisely locate enhancers and promoters. MDS is a multivariate analysis method that reproduces the original coordinates from the distance matrix between elements. We used Hi-C data of cultured osteosarcoma cells and applied MDS as the distance matrix of the genome. In addition, we selected columns 2 and 3 of the orthogonal matrix U as the desired structure. Overall, the DNA loops from the reconstructed genome structure contained bioprocesses involved in transcription, such as the pre-transcriptional initiation complex and RNA polymerase II initiation complex, and transcription factors involved in cancer, such as Foxm1 and CREB3. Therefore, our results are consistent with the biological findings. Our method is suitable for identifying enhancers and promoters in the genome.



2020 ◽  
Author(s):  
Liwei Liu ◽  
Huili Yao

AbstractIn recent years, with the development of high-throughput chromosome conformation capture (Hi-C) technology and the reduction of high-throughput sequencing cost, the data volume of whole-genome interaction has increased rapidly, and the resolution of interaction map keeps improving. Great progress has been made in the research of 3D structure modeling of chromosomes and genomes. Several methods have been proposed to construct the chromosome structure from chromosome conformation capture data. Based on the Hi-C data, this paper analyses the relevant literature of chromosome 3D structure reconstruction and it summarizes the principle of 3DMAX, which is a classical algorithm to construct the 3D structure of a chromosome. In this paper, we introduce a new gradient ascent optimization algorithm called XNadam that is a variant of Nadam optimization method. When XNadam is applied to 3DMax algorithm, the performance of 3DMax algorithm can be improved, which can be used to predict the three-dimensional structure of a chromosome.Author summaryThe exploration of the three-dimensional structure of chromosomes has gradually become a necessary means to understand the relationship between genome function and gene regulation. An important problem in the construction of three-dimensional model is how to use the interaction map. Usually, the interaction frequency can be transformed into the spatial distance according to the deterministic or non-deterministic function relationship, and the interaction frequency can be weighted as weight in the objective function of the optimization problem. When the frequency of interaction is weighted as weight in the objective function of the optimization problem, what kind of optimization method is used to optimize the objective function is the problem we consider. In order to solve this problem, we provide an improved stochastic gradient ascent optimization algorithm(XNadam). The XNadam optimization algorithm combined with maximum likelihood algorithm is applied to high resolution Hi-C data set to infer 3D chromosome structure.



2019 ◽  
Author(s):  
Priyojit Das ◽  
Tongye Shen ◽  
Rachel Patton McCord

AbstractBackgroundThe nonrandom radial organization of eukaryotic chromosome territories (CTs) inside the nucleus plays an important role in nuclear functional compartmentalization. Increasingly, chromosome conformation capture (Hi-C) based approaches are being used to characterize the genome structure of many cell types and conditions. Computational methods to extract 3D arrangements of CTs from this type of pairwise contact data will thus increase our ability to analyze CT organization in a wider variety of biological situations.ResultsA number of full-scale polymer models have successfully reconstructed the 3D structure of chromosome territories from Hi-C. To supplement such methods, we explore alternative, direct, and less computationally intensive approaches to capture radial CT organization from Hi-C data. We show that we can infer relative chromo-some ordering using PCA on a thresholded inter-chromosomal contact matrix. We simulate an ensemble of possible CT arrangements using a force-directed network layout algorithm and propose an approach to integrate additional chromosome properties into our predictions. Our CT radial organization predictions have a high correlation with microscopy imaging data for various cell nucleus geometries (lymphoblastoid, skin fibroblast, and breast epithelial cells), and we can capture previously documented changes in senescent and progeria cells.ConclusionsOur analysis approaches provide rapid and modular approaches to screen for alterations in CT organization across widely available Hi-C data. We demon-strate which stages of the approach can extract meaningful information, and also de-scribe limitations of pairwise contacts alone to predict absolute 3D positions.



2017 ◽  
Vol 20 (4) ◽  
pp. 1205-1214
Author(s):  
Jincheol Park ◽  
Shili Lin

Abstract How chromosomes fold and how distal genomic elements interact with one another at a genomic scale have been actively pursued in the past decade following the seminal work describing the Chromosome Conformation Capture (3C) assay. Essentially, 3C-based technologies produce two-dimensional (2D) contact maps that capture interactions between genomic fragments. Accordingly, a plethora of analytical methods have been proposed to take a 2D contact map as input to recapitulate the underlying whole genome three-dimensional (3D) structure of the chromatin. However, their performance in terms of several factors, including data resolution and ability to handle contact map features, have not been sufficiently evaluated. This task is taken up in this article, in which we consider several recent and/or well-regarded methods, both optimization-based and model-based, for their aptness of producing 3D structures using contact maps generated based on a population of cells. These methods are evaluated and compared using both simulated and real data. Several criteria have been used. For simulated data sets, the focus is on accurate recapitulation of the entire structure given the existence of the gold standard. For real data sets, comparison with distances measured by Florescence in situ Hybridization and consistency with several genomic features of known biological functions are examined.



2021 ◽  
Author(s):  
Van Hovenga ◽  
Oluwatosin Oluwadare ◽  
Jugal Kalita

Chromosome conformation capture (3C) is a method of measuring chromosome topology in terms of loci interaction. The Hi-C method is a derivative of 3C that allows for genome wide quantification of chromosome interaction. From such interaction data, it is possible to infer the three-dimensional (3D) structure of the underlying chromosome. In this paper, we use a node embedding algorithm and a graph neural network to predict the 3D coordinates of each genomic loci from the corresponding Hi-C contact data. Unlike other chromosome structure prediction methods, our method can generalize a single model across Hi-C resolutions, multiple restriction enzymes, and multiple cell populations while maintaining reconstruction accuracy. We derive these results using three separate Hi-C data sets from the GM12878, GM06990, and K562 cell lines. We also compare the reconstruction accuracy of our method to four other existing methods and show that our method yields superior performance. Our algorithm outperforms the state-of-the-art methods in the accuracy of prediction and introduces a novel method for 3D structure prediction from Hi-C data.



Author(s):  
Zijian Guo ◽  
Tanghong Liu ◽  
Wenhui Li ◽  
Yutao Xia

The present work focuses on the aerodynamic problems resulting from a high-speed train (HST) passing through a tunnel. Numerical simulations were employed to obtain the numerical results, and they were verified by a moving-model test. Two responses, [Formula: see text] (coefficient of the peak-to-peak pressure of a single fluctuation) and[Formula: see text] (pressure value of micro-pressure wave), were studied with regard to the three building parameters of the portal-hat buffer structure of the tunnel entrance and exit. The MOPSO (multi-objective particle swarm optimization) method was employed to solve the optimization problem in order to find the minimum [Formula: see text] and[Formula: see text]. Results showed that the effects of the three design parameters on [Formula: see text] were not monotonous, and the influences of[Formula: see text] (the oblique angle of the portal) and [Formula: see text] (the height of the hat structure) were more significant than that of[Formula: see text] (the angle between the vertical line of the portal and the hat). Monotonically decreasing responses were found in [Formula: see text] for [Formula: see text] and[Formula: see text]. The Pareto front of [Formula: see text] and[Formula: see text]was obtained. The ideal single-objective optimums for each response located at the ends of the Pareto front had values of 1.0560 for [Formula: see text] and 101.8 Pa for[Formula: see text].



2018 ◽  
Vol 2018 ◽  
pp. 1-8
Author(s):  
Lingyun Lu ◽  
Tian Wang ◽  
Wei Ni ◽  
Kai Li ◽  
Bo Gao

This paper presents a suboptimal approach for resource allocation of massive MIMO-OFDMA systems for high-speed train (HST) applications. An optimization problem is formulated to alleviate the severe Doppler effect and maximize the energy efficiency (EE) of the system. We propose to decouple the problem between the allocations of antennas, subcarriers, and transmit powers and solve the problem by carrying out the allocations separately and iteratively in an alternating manner. Fast convergence can be achieved for the proposed approach within only several iterations. Simulation results show that the proposed algorithm is superior to existing techniques in terms of system EE and throughput in different system configurations of HST applications.



2017 ◽  
Author(s):  
◽  
Tuan Anh Trieu

[ACCESS RESTRICTED TO THE UNIVERSITY OF MISSOURI AT AUTHOR'S REQUEST.] Different cell types of an organism have the same DNA sequence, but they can function differently because their difference in 3D organization allows them to express different genes and has different cellular functions. Understanding the 3D organization of the genome is the key to understand functions of the cell. Chromosome conformation capture techniques like Hi-C and TCC that can capture interactions between proximal chromosome fragments have allowed the study of 3D genome organization in high resolution and high through-put. My work focuses on developing computational methods to reconstruct 3D genome structures from Hi-C data. I presented three methods to reconstruct 3D genome and chromosome structures. The first method can build 3D genome models from soft constraints of contacts and non-contacts. This method utilizes the concept of contact and non-contact to reconstruct 3D models without translating interaction frequencies into physical distances. The translation is commonly used by other methods even though it makes a strong assumption about the relationship between interaction frequencies and physical distances. In synthetic dataset, when the relationship was known, my method performed comparably with other methods assuming the relationship. This shows the potential of my method for real Hi-C datasets where the relationship is unknown. The limitation of the method is that it has parameters requiring manual adjustment. I developed the second method to reconstruct 3D genome models. This method utilizes a commonly used function to translate interaction frequencies to physical distances to build 3D models. I proposed a novel way to derive soft constraints to handle inconsistency in the data and to make the method robust. Building 3D models at high resolution is a more challenging problem as the number of constraints is small and the feasible space is larger. I introduced a third method to build 3D chromosome models at high resolution. The method reconstructs models at low resolution and then uses them to guide the reconstruction of models at high resolution. The last part of my work is the development of a comprehensive tool with intuitive graphic user interface to analyze Hi-C data, reconstruct and analyze 3D models.



Sign in / Sign up

Export Citation Format

Share Document