1
|
龚 海, 张 司, 张 晓. [An identification method of chromatin topological associated domains based on spatial density clustering]. SHENG WU YI XUE GONG CHENG XUE ZA ZHI = JOURNAL OF BIOMEDICAL ENGINEERING = SHENGWU YIXUE GONGCHENGXUE ZAZHI 2024; 41:552-559. [PMID: 38932542 PMCID: PMC11208664 DOI: 10.7507/1001-5515.202311059] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Subscribe] [Scholar Register] [Received: 11/28/2023] [Revised: 03/05/2024] [Indexed: 06/28/2024]
Abstract
The rapid development of high-throughput chromatin conformation capture (Hi-C) technology provides rich genomic interaction data between chromosomal loci for chromatin structure analysis. However, existing methods for identifying topologically associated domains (TADs) based on Hi-C data suffer from low accuracy and sensitivity to parameters. In this context, a TAD identification method based on spatial density clustering was designed and implemented in this paper. The method preprocessed the raw Hi-C data to obtain normalized Hi-C contact matrix data. Then, it computed the distance matrix between loci, generated a reachability graph based on the core distance and reachability distance of loci, and extracted clustering clusters. Finally, it extracted TAD boundaries based on clustering results. This method could identify TAD structures with higher coherence, and TAD boundaries were enriched with more ChIP-seq factors. Experimental results demonstrate that our method has advantages such as higher accuracy and practical significance in TAD identification.
Collapse
Affiliation(s)
- 海燕 龚
- 北京科技大学 新材料技术研究院(北京 100083)Institute for Advanced Materials and Technology, University of Science and Technology Beijing, Beijing 100083, P. R. China
- 北京科技大学 计算机与通信工程学院(北京 100083)School of Computer and Communication Engineering, University of Science and Technology Beijing, Beijing 100083, P. R. China
| | - 司臣 张
- 北京科技大学 新材料技术研究院(北京 100083)Institute for Advanced Materials and Technology, University of Science and Technology Beijing, Beijing 100083, P. R. China
| | - 晓彤 张
- 北京科技大学 新材料技术研究院(北京 100083)Institute for Advanced Materials and Technology, University of Science and Technology Beijing, Beijing 100083, P. R. China
- 北京科技大学 计算机与通信工程学院(北京 100083)School of Computer and Communication Engineering, University of Science and Technology Beijing, Beijing 100083, P. R. China
| |
Collapse
|
2
|
Chen B, Ren C, Ouyang Z, Xu J, Xu K, Li Y, Guo H, Bai X, Tian M, Xu X, Wang Y, Li H, Bo X, Chen H. Stratifying TAD boundaries pinpoints focal genomic regions of regulation, damage, and repair. Brief Bioinform 2024; 25:bbae306. [PMID: 38935071 PMCID: PMC11210073 DOI: 10.1093/bib/bbae306] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2024] [Revised: 06/01/2024] [Accepted: 06/13/2024] [Indexed: 06/28/2024] Open
Abstract
Advances in chromatin mapping have exposed the complex chromatin hierarchical organization in mammals, including topologically associating domains (TADs) and their substructures, yet the functional implications of this hierarchy in gene regulation and disease progression are not fully elucidated. Our study delves into the phenomenon of shared TAD boundaries, which are pivotal in maintaining the hierarchical chromatin structure and regulating gene activity. By integrating high-resolution Hi-C data, chromatin accessibility, and DNA double-strand breaks (DSBs) data from various cell lines, we systematically explore the complex regulatory landscape at high-level TAD boundaries. Our findings indicate that these boundaries are not only key architectural elements but also vibrant hubs, enriched with functionally crucial genes and complex transcription factor binding site-clustered regions. Moreover, they exhibit a pronounced enrichment of DSBs, suggesting a nuanced interplay between transcriptional regulation and genomic stability. Our research provides novel insights into the intricate relationship between the 3D genome structure, gene regulation, and DNA repair mechanisms, highlighting the role of shared TAD boundaries in maintaining genomic integrity and resilience against perturbations. The implications of our findings extend to understanding the complexities of genomic diseases and open new avenues for therapeutic interventions targeting the structural and functional integrity of TAD boundaries.
Collapse
Affiliation(s)
- Bijia Chen
- Academy of Military Medical Sciences, Beijing 100850, China
| | - Chao Ren
- Academy of Military Medical Sciences, Beijing 100850, China
| | - Zhangyi Ouyang
- Academy of Military Medical Sciences, Beijing 100850, China
| | - Jingxuan Xu
- Key Laboratory of Carcinogenesis and Translational Research (Ministry of Education/Beijing), Department of Gastrointestinal Surgery, Peking University Cancer Hospital & Institute, Beijing 100142, China
| | - Kang Xu
- School of Software, Shandong University, Jinan 250101, China
| | - Yaru Li
- Academy of Military Medical Sciences, Beijing 100850, China
| | - Hejiang Guo
- Academy of Military Medical Sciences, Beijing 100850, China
| | - Xuemei Bai
- Academy of Military Medical Sciences, Beijing 100850, China
| | - Mengge Tian
- The First Affiliated Hospital of Harbin Medical University, Harbin 150001, China
| | - Xiang Xu
- Academy of Military Medical Sciences, Beijing 100850, China
| | - Yuyang Wang
- College of Computer and Data Science, Fuzhou University, Fuzhou 350108, China
| | - Hao Li
- Academy of Military Medical Sciences, Beijing 100850, China
| | - Xiaochen Bo
- Academy of Military Medical Sciences, Beijing 100850, China
| | - Hebing Chen
- Academy of Military Medical Sciences, Beijing 100850, China
| |
Collapse
|
3
|
Xu J, Xu X, Huang D, Luo Y, Lin L, Bai X, Zheng Y, Yang Q, Cheng Y, Huang A, Shi J, Bo X, Gu J, Chen H. A comprehensive benchmarking with interpretation and operational guidance for the hierarchy of topologically associating domains. Nat Commun 2024; 15:4376. [PMID: 38782890 PMCID: PMC11116433 DOI: 10.1038/s41467-024-48593-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2023] [Accepted: 05/03/2024] [Indexed: 05/25/2024] Open
Abstract
Topologically associating domains (TADs), megabase-scale features of chromatin spatial architecture, are organized in a domain-within-domain TAD hierarchy. Within TADs, the inner and smaller subTADs not only manifest cell-to-cell variability, but also precisely regulate transcription and differentiation. Although over 20 TAD callers are able to detect TAD, their usability in biomedicine is confined by a disagreement of outputs and a limit in understanding TAD hierarchy. We compare 13 computational tools across various conditions and develop a metric to evaluate the similarity of TAD hierarchy. Although outputs of TAD hierarchy at each level vary among callers, data resolutions, sequencing depths, and matrices normalization, they are more consistent when they have a higher similarity of larger TADs. We present comprehensive benchmarking of TAD hierarchy callers and operational guidance to researchers of life science researchers. Moreover, by simulating the mixing of different types of cells, we confirm that TAD hierarchy is generated not simply from stacking Hi-C heatmaps of heterogeneous cells. Finally, we propose an air conditioner model to decipher the role of TAD hierarchy in transcription.
Collapse
Affiliation(s)
- Jingxuan Xu
- Key Laboratory of Carcinogenesis and Translational Research (Ministry of Education/Beijing), Department of Gastrointestinal Surgery, Peking University Cancer Hospital & Institute, Beijing, 100142, China
| | - Xiang Xu
- Academy of Military Medical Science, Beijing, 100850, China
| | - Dandan Huang
- Department of Oncology, Peking University Shougang Hospital, Beijing, China
- Center for Precision Diagnosis and Treatment of Colorectal Cancer and Inflammatory Diseases, Peking University Health Science Center, Beijing, China
| | - Yawen Luo
- Academy of Military Medical Science, Beijing, 100850, China
| | - Lin Lin
- Academy of Military Medical Science, Beijing, 100850, China
- School of Computer Science and Information Technology& KLAS, Northeast Normal University, Changchun, China
| | - Xuemei Bai
- Academy of Military Medical Science, Beijing, 100850, China
| | - Yang Zheng
- Academy of Military Medical Science, Beijing, 100850, China
| | - Qian Yang
- Academy of Military Medical Science, Beijing, 100850, China
| | - Yu Cheng
- Key Laboratory of Carcinogenesis and Translational Research (Ministry of Education/Beijing), Department of Gastrointestinal Surgery, Peking University Cancer Hospital & Institute, Beijing, 100142, China
| | - An Huang
- Key Laboratory of Carcinogenesis and Translational Research (Ministry of Education/Beijing), Department of Gastrointestinal Surgery, Peking University Cancer Hospital & Institute, Beijing, 100142, China
| | - Jingyi Shi
- Key Laboratory of Carcinogenesis and Translational Research (Ministry of Education/Beijing), Department of Gastrointestinal Surgery, Peking University Cancer Hospital & Institute, Beijing, 100142, China
| | - Xiaochen Bo
- Academy of Military Medical Science, Beijing, 100850, China.
| | - Jin Gu
- Key Laboratory of Carcinogenesis and Translational Research (Ministry of Education/Beijing), Department of Gastrointestinal Surgery, Peking University Cancer Hospital & Institute, Beijing, 100142, China.
- Department of Oncology, Peking University Shougang Hospital, Beijing, China.
- Center for Precision Diagnosis and Treatment of Colorectal Cancer and Inflammatory Diseases, Peking University Health Science Center, Beijing, China.
- Peking-Tsinghua Center for Life Sciences, Peking University, Beijing, China.
- Peking University International Cancer Institute, Beijing, China.
| | - Hebing Chen
- Academy of Military Medical Science, Beijing, 100850, China.
| |
Collapse
|
4
|
Zhang B, Long Y, Pei L, Huang X, Li B, Han B, Zhang M, Lindsey K, Zhang X, Wang M, Yang X. Drought response revealed by chromatin organization variation and transcriptional regulation in cotton. BMC Biol 2024; 22:114. [PMID: 38764013 PMCID: PMC11103878 DOI: 10.1186/s12915-024-01906-0] [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: 11/14/2023] [Accepted: 04/29/2024] [Indexed: 05/21/2024] Open
Abstract
BACKGROUND Cotton is a major world cash crop and an important source of natural fiber, oil, and protein. Drought stress is becoming a restrictive factor affecting cotton production. To facilitate the development of drought-tolerant cotton varieties, it is necessary to study the molecular mechanism of drought stress response by exploring key drought-resistant genes and related regulatory factors. RESULTS In this study, two cotton varieties, ZY007 (drought-sensitive) and ZY168 (drought-tolerant), showing obvious phenotypic differences under drought stress, were selected. A total of 25,898 drought-induced genes were identified, exhibiting significant enrichment in pathways related to plant stress responses. Under drought induction, At subgenome expression bias was observed at the whole-genome level, which may be due to stronger inhibition of Dt subgenome expression. A gene co-expression module that was significantly associated with drought resistance was identified. About 90% of topologically associating domain (TAD) boundaries were stable, and 6613 TAD variation events were identified between the two varieties under drought. We identified 92 genes in ZY007 and 98 in ZY168 related to chromatin 3D structural variation and induced by drought stress. These genes are closely linked to the cotton response to drought stress through canonical hormone-responsive pathways, modulation of kinase and phosphatase activities, facilitation of calcium ion transport, and other related molecular mechanisms. CONCLUSIONS These results lay a foundation for elucidating the molecular mechanism of the cotton drought response and provide important regulatory locus and gene resources for the future molecular breeding of drought-resistant cotton varieties.
Collapse
Affiliation(s)
- Boyang Zhang
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Yuexuan Long
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Liuling Pei
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Xianhui Huang
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Baoqi Li
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Bei Han
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Mengmeng Zhang
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Keith Lindsey
- Department of Biosciences, Durham University, South Road, Durham, DH1 3LE, UK
| | - Xianlong Zhang
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Maojun Wang
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China.
| | - Xiyan Yang
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China.
| |
Collapse
|
5
|
Yan Z, Yang J, Wei WT, Zhou ML, Mo DX, Wan X, Ma R, Wu MM, Huang JH, Liu YJ, Lv FH, Li MH. A time-resolved multi-omics atlas of transcriptional regulation in response to high-altitude hypoxia across whole-body tissues. Nat Commun 2024; 15:3970. [PMID: 38730227 PMCID: PMC11087590 DOI: 10.1038/s41467-024-48261-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2023] [Accepted: 04/23/2024] [Indexed: 05/12/2024] Open
Abstract
High-altitude hypoxia acclimatization requires whole-body physiological regulation in highland immigrants, but the underlying genetic mechanism has not been clarified. Here we use sheep as an animal model for low-to-high altitude translocation. We generate multi-omics data including whole-genome sequences, time-resolved bulk RNA-Seq, ATAC-Seq and single-cell RNA-Seq from multiple tissues as well as phenotypic data from 20 bio-indicators. We characterize transcriptional changes of all genes in each tissue, and examine multi-tissue temporal dynamics and transcriptional interactions among genes. Particularly, we identify critical functional genes regulating the short response to hypoxia in each tissue (e.g., PARG in the cerebellum and HMOX1 in the colon). We further identify TAD-constrained cis-regulatory elements, which suppress the transcriptional activity of most genes under hypoxia. Phenotypic and transcriptional evidence indicate that antenatal hypoxia could improve hypoxia tolerance in offspring. Furthermore, we provide time-series expression data of candidate genes associated with human mountain sickness (e.g., BMPR2) and high-altitude adaptation (e.g., HIF1A). Our study provides valuable resources and insights for future hypoxia-related studies in mammals.
Collapse
Affiliation(s)
- Ze Yan
- State Key Laboratory of Animal Biotech Breeding, China Agricultural University, Beijing, 100193, China
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Ji Yang
- State Key Laboratory of Animal Biotech Breeding, China Agricultural University, Beijing, 100193, China
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Wen-Tian Wei
- State Key Laboratory of Animal Biotech Breeding, China Agricultural University, Beijing, 100193, China
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Ming-Liang Zhou
- Sichuan Academy of Grassland Science, Chengdu, 611743, China
| | - Dong-Xin Mo
- State Key Laboratory of Animal Biotech Breeding, China Agricultural University, Beijing, 100193, China
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Xing Wan
- State Key Laboratory of Animal Biotech Breeding, China Agricultural University, Beijing, 100193, China
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Rui Ma
- State Key Laboratory of Animal Biotech Breeding, China Agricultural University, Beijing, 100193, China
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Mei-Ming Wu
- State Key Laboratory of Animal Biotech Breeding, China Agricultural University, Beijing, 100193, China
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Jia-Hui Huang
- State Key Laboratory of Animal Biotech Breeding, China Agricultural University, Beijing, 100193, China
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Ya-Jing Liu
- State Key Laboratory of Animal Biotech Breeding, China Agricultural University, Beijing, 100193, China
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Feng-Hua Lv
- State Key Laboratory of Animal Biotech Breeding, China Agricultural University, Beijing, 100193, China
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Meng-Hua Li
- State Key Laboratory of Animal Biotech Breeding, China Agricultural University, Beijing, 100193, China.
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China.
| |
Collapse
|
6
|
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
|
7
|
Zhang L, Zhao R, Liang J, Cai X, Zhang L, Guo H, Zhang Z, Wu J, Wang X. BL-Hi-C reveals the 3D genome structure of Brassica crops with high sensitivity. HORTICULTURE RESEARCH 2024; 11:uhae017. [PMID: 38464474 PMCID: PMC10923644 DOI: 10.1093/hr/uhae017] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/01/2023] [Accepted: 01/03/2024] [Indexed: 03/12/2024]
Abstract
High-throughput Chromatin Conformation Capture (Hi-C) technologies can be used to investigate the three-dimensional genomic structure of plants. However, the practical utility of these technologies is impeded by significant background noise, hindering their capability in detecting fine 3D genomic structures. In this study, we optimized the Bridge Linker Hi-C technology (BL-Hi-C) to comprehensively investigate the 3D chromatin landscape of Brassica rapa and Brassica oleracea. The Bouquet configuration of both B. rapa and B. oleracea was elucidated through the construction of a 3D genome simulation. The optimized BL-Hi-C exhibited lower background noise compared to conventional Hi-C methods. Taking this advantage, we used BL-Hi-C to identify FLC gene loops in Arabidopsis, B. rapa, and B. oleracea. We observed that gene loops of FLC2 exhibited conservation across Arabidopsis, B. rapa, and B. oleracea. While gene loops of syntenic FLCs exhibited conservation across B. rapa and B. oleracea, variations in gene loops were evident among multiple paralogs FLCs within the same species. Collectively, our findings highlight the high sensitivity of optimized BL-Hi-C as a powerful tool for investigating the fine 3D genomic organization.
Collapse
Affiliation(s)
- Lupeng Zhang
- State Key Laboratory of Vegetable Biobreeding, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing 100081, China
| | - Ranze Zhao
- State Key Laboratory of Vegetable Biobreeding, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing 100081, China
| | - Jianli Liang
- State Key Laboratory of Vegetable Biobreeding, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing 100081, China
| | - Xu Cai
- State Key Laboratory of Vegetable Biobreeding, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing 100081, China
| | - Lei Zhang
- State Key Laboratory of Vegetable Biobreeding, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing 100081, China
| | - Huiling Guo
- State Key Laboratory of Vegetable Biobreeding, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing 100081, China
| | - Zhicheng Zhang
- State Key Laboratory of Vegetable Biobreeding, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing 100081, China
| | - Jian Wu
- State Key Laboratory of Vegetable Biobreeding, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing 100081, China
| | - Xiaowu Wang
- State Key Laboratory of Vegetable Biobreeding, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing 100081, China
| |
Collapse
|
8
|
Xu MRX, Liao ZY, Brock JR, Du K, Li GY, Chen ZQ, Wang YH, Gao ZN, Agarwal G, Wei KHC, Shao F, Pang S, Platts AE, van de Velde J, Lin HM, Teresi SJ, Bird K, Niederhuth CE, Xu JG, Yu GH, Yang JY, Dai SF, Nelson A, Braasch I, Zhang XG, Schartl M, Edger PP, Han MJ, Zhang HH. Maternal dominance contributes to subgenome differentiation in allopolyploid fishes. Nat Commun 2023; 14:8357. [PMID: 38102128 PMCID: PMC10724154 DOI: 10.1038/s41467-023-43740-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2023] [Accepted: 11/17/2023] [Indexed: 12/17/2023] Open
Abstract
Teleost fishes, which are the largest and most diverse group of living vertebrates, have a rich history of ancient and recent polyploidy. Previous studies of allotetraploid common carp and goldfish (cyprinids) reported a dominant subgenome, which is more expressed and exhibits biased gene retention. However, the underlying mechanisms contributing to observed 'subgenome dominance' remains poorly understood. Here we report high-quality genomes of twenty-one cyprinids to investigate the origin and subsequent subgenome evolution patterns following three independent allopolyploidy events. We identify the closest extant relatives of the diploid progenitor species, investigate genetic and epigenetic differences among subgenomes, and conclude that observed subgenome dominance patterns are likely due to a combination of maternal dominance and transposable element densities in each polyploid. These findings provide an important foundation to understanding subgenome dominance patterns observed in teleost fishes, and ultimately the role of polyploidy in contributing to evolutionary innovations.
Collapse
Affiliation(s)
- Min-Rui-Xuan Xu
- College of Pharmacy and Life Science, Jiujiang University, Jiujiang, China
| | - Zhen-Yang Liao
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
| | - Jordan R Brock
- Department of Horticulture, Michigan State University, East Lansing, MI, USA
| | - Kang Du
- The Xiphophorus Genetic Stock Center, Texas State University, San Marcos, TX, USA
| | - Guo-Yin Li
- College of Life Science and Agronomy, Zhoukou Normal University, Zhoukou, Henan, China
| | | | - Ying-Hao Wang
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
| | - Zhong-Nan Gao
- College of Pharmacy and Life Science, Jiujiang University, Jiujiang, China
| | - Gaurav Agarwal
- Department of Plant Biology, Michigan State University, East Lansing, MI, USA
| | - Kevin H-C Wei
- Department of Integrative Biology, University of California Berkeley, Berkeley, CA, USA
- Department of Zoology, University of British Columbia, Vancouver, British Columbia, Canada
| | - Feng Shao
- Key Laboratory of Freshwater Fish Reproduction and Development (Ministry of Education), Southwest University, School of Life Sciences, Chongqing, China
| | | | - Adrian E Platts
- Department of Horticulture, Michigan State University, East Lansing, MI, USA
| | - Jozefien van de Velde
- Department of Chromosome Biology, Max Planck Institute for Plant Breeding Research, Cologne, Germany
| | - Hong-Min Lin
- College of Pharmacy and Life Science, Jiujiang University, Jiujiang, China
| | - Scott J Teresi
- Department of Horticulture, Michigan State University, East Lansing, MI, USA
| | - Kevin Bird
- Department of Horticulture, Michigan State University, East Lansing, MI, USA
| | - Chad E Niederhuth
- Department of Plant Biology, Michigan State University, East Lansing, MI, USA
| | - Jin-Gen Xu
- Jiujiang Academy of Agricultural Sciences, Jiujiang, China
| | - Guo-Hua Yu
- College of Pharmacy and Life Science, Jiujiang University, Jiujiang, China
| | - Jian-Yuan Yang
- College of Pharmacy and Life Science, Jiujiang University, Jiujiang, China
| | - Si-Fa Dai
- College of Pharmacy and Life Science, Jiujiang University, Jiujiang, China
| | | | - Ingo Braasch
- Department of Integrative Biology, Michigan State University, East Lansing, MI, USA
| | - Xiao-Gu Zhang
- College of Pharmacy and Life Science, Jiujiang University, Jiujiang, China.
| | - Manfred Schartl
- The Xiphophorus Genetic Stock Center, Texas State University, San Marcos, TX, USA.
- Developmental Biochemistry, Biocenter, University of Würzburg, Würzburg, Bayern, Germany.
| | - Patrick P Edger
- Department of Horticulture, Michigan State University, East Lansing, MI, USA.
| | - Min-Jin Han
- State Key Laboratory of Resource Insects, Key Laboratory for Sericulture Functional Genomics and Biotechnology of Agricultural Ministry, Southwest University, Chongqing, China.
| | - Hua-Hao Zhang
- College of Pharmacy and Life Science, Jiujiang University, Jiujiang, China.
| |
Collapse
|
9
|
Shen Y, Yu L, Qiu Y, Zhang T, Kingsford C. Improving Hi-C contact matrices using genome graphs. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.11.08.566275. [PMID: 37986943 PMCID: PMC10659349 DOI: 10.1101/2023.11.08.566275] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/22/2023]
Abstract
Three-dimensional chromosome structure plays an important role in fundamental genomic functions. Hi-C, a high-throughput, sequencing-based technique, has drastically expanded our comprehension of 3D chromosome structures. The first step of Hi-C analysis pipeline involves mapping sequencing reads from Hi-C to linear reference genomes. However, the linear reference genome does not incorporate genetic variation information, which can lead to incorrect read alignments, especially when analyzing samples with substantial genomic differences from the reference such as cancer samples. Using genome graphs as the reference facilitates more accurate mapping of reads, however, new algorithms are required for inferring linear genomes from Hi-C reads mapped on genome graphs and constructing corresponding Hi-C contact matrices, which is a prerequisite for the subsequent steps of the Hi-C analysis such as identifying topologically associated domains and calling chromatin loops. We introduce the problem of genome sequence inference from Hi-C data mediated by genome graphs. We formalize this problem, show the hardness of solving this problem, and introduce a novel heuristic algorithm specifically tailored to this problem. We provide a theoretical analysis to evaluate the efficacy of our algorithm. Finally, our empirical experiments indicate that the linear genomes inferred from our method lead to the creation of improved Hi-C contact matrices. These enhanced matrices show a reduction in erroneous patterns caused by structural variations and are more effective in accurately capturing the structures of topologically associated domains.
Collapse
Affiliation(s)
- Yihang Shen
- Computational Biology Department, School of Computer Science, Carnegie Mellon University, 5000 Forbes Avenue, Pittsburgh, PA
| | - Lingge Yu
- Computational Biology Department, School of Computer Science, Carnegie Mellon University, 5000 Forbes Avenue, Pittsburgh, PA
| | - Yutong Qiu
- Computational Biology Department, School of Computer Science, Carnegie Mellon University, 5000 Forbes Avenue, Pittsburgh, PA
| | - Tianyu Zhang
- Department of Statistics and Data Science, Carnegie Mellon University, 5000 Forbes Avenue, Pittsburgh, PA
| | - Carl Kingsford
- Computational Biology Department, School of Computer Science, Carnegie Mellon University, 5000 Forbes Avenue, Pittsburgh, PA
| |
Collapse
|
10
|
Lee H, Seo P. Accessible gene borders establish a core structural unit for chromatin architecture in Arabidopsis. Nucleic Acids Res 2023; 51:10261-10277. [PMID: 37884483 PMCID: PMC10602878 DOI: 10.1093/nar/gkad710] [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: 03/22/2023] [Revised: 08/08/2023] [Accepted: 08/16/2023] [Indexed: 10/28/2023] Open
Abstract
Three-dimensional (3D) chromatin structure is linked to transcriptional regulation in multicellular eukaryotes including plants. Taking advantage of high-resolution Hi-C (high-throughput chromatin conformation capture), we detected a small structural unit with 3D chromatin architecture in the Arabidopsis genome, which lacks topologically associating domains, and also in the genomes of tomato, maize, and Marchantia polymorpha. The 3D folding domain unit was usually established around an individual gene and was dependent on chromatin accessibility at the transcription start site (TSS) and transcription end site (TES). We also observed larger contact domains containing two or more neighboring genes, which were dependent on accessible border regions. Binding of transcription factors to accessible TSS/TES regions formed these gene domains. We successfully simulated these Hi-C contact maps via computational modeling using chromatin accessibility as input. Our results demonstrate that gene domains establish basic 3D chromatin architecture units that likely contribute to higher-order 3D genome folding in plants.
Collapse
Affiliation(s)
- Hongwoo Lee
- Department of Chemistry, Seoul National University, Seoul 08826, Korea
| | - Pil Joon Seo
- Department of Chemistry, Seoul National University, Seoul 08826, Korea
- Plant Genomics and Breeding Institute, Seoul National University, Seoul 08826, Korea
| |
Collapse
|
11
|
Gong H, Zhang D, Zhang X. TOAST: A novel method for identifying topologically associated domains based on graph auto-encoders and clustering. Comput Struct Biotechnol J 2023; 21:4759-4768. [PMID: 37822562 PMCID: PMC10562672 DOI: 10.1016/j.csbj.2023.09.019] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2023] [Revised: 09/16/2023] [Accepted: 09/16/2023] [Indexed: 10/13/2023] Open
Abstract
Topologically associated domains (TADs) play a pivotal role in disease detection. This study introduces a novel TADs recognition approach named TOAST, leveraging graph auto-encoders and clustering techniques. TOAST conceptualizes each genomic bin as a node of a graph and employs the Hi-C contact matrix as the graph's adjacency matrix. By employing graph auto-encoders, TOAST generates informative embeddings as features. Subsequently, the unsupervised clustering algorithm HDBSCAN is utilized to assign labels to each genomic bin, facilitating the identification of contiguous regions with the same label as TADs. Our experimental analysis of several simulated Hi-C data sets shows that TOAST can quickly and accurately identify TADs from different types of simulated Hi-C contact matrices, outperforming existing algorithms. We also determined the anchoring ratio of TAD boundaries by analyzing different TAD recognition algorithms, and obtained an average ratio of anchoring CTCF, SMC3, RAD21, POLR2A, H3K36me3, H3K9me3, H3K4me3, H3K4me1, Enhancer, and Promoters of 0.66, 0.47, 0.54, 0.27, 0.24, 0.12, 0.32, 0.41, 0.26, and 0.13, respectively. In conclusion, TOAST is a method that can quickly identify TAD boundary parameters that are easy to understand and have important biological significance. The TOAST web server can be accessed via http://223.223.185.189:4005/. The code of TOAST is available online at https://github.com/ghaiyan/TOAST.
Collapse
Affiliation(s)
- Haiyan Gong
- Institute for Advanced Materials and Technology, University of Science and Technology Beijing, Beijing, 100083, China
- School of Computer and Communication Engineering, Beijing Advanced Innovation Center for Materials Genome Engineering, University of Science and Technology Beijing, Beijing, 100083, China
- Shunde innovation School, University of Science and Technology Beijing, Foshan, 528399, Guangdong, China
| | - Dawei Zhang
- Institute for Advanced Materials and Technology, University of Science and Technology Beijing, Beijing, 100083, China
- Shunde innovation School, University of Science and Technology Beijing, Foshan, 528399, Guangdong, China
| | - Xiaotong Zhang
- School of Computer and Communication Engineering, Beijing Advanced Innovation Center for Materials Genome Engineering, University of Science and Technology Beijing, Beijing, 100083, China
- Shunde innovation School, University of Science and Technology Beijing, Foshan, 528399, Guangdong, China
| |
Collapse
|
12
|
Raffo A, Paulsen J. The shape of chromatin: insights from computational recognition of geometric patterns in Hi-C data. Brief Bioinform 2023; 24:bbad302. [PMID: 37646128 PMCID: PMC10516369 DOI: 10.1093/bib/bbad302] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2023] [Revised: 07/05/2023] [Accepted: 08/03/2023] [Indexed: 09/01/2023] Open
Abstract
The three-dimensional organization of chromatin plays a crucial role in gene regulation and cellular processes like deoxyribonucleic acid (DNA) transcription, replication and repair. Hi-C and related techniques provide detailed views of spatial proximities within the nucleus. However, data analysis is challenging partially due to a lack of well-defined, underpinning mathematical frameworks. Recently, recognizing and analyzing geometric patterns in Hi-C data has emerged as a powerful approach. This review provides a summary of algorithms for automatic recognition and analysis of geometric patterns in Hi-C data and their correspondence with chromatin structure. We classify existing algorithms on the basis of the data representation and pattern recognition paradigm they make use of. Finally, we outline some of the challenges ahead and promising future directions.
Collapse
Affiliation(s)
- Andrea Raffo
- Department of Biosciences, University of Oslo, 0316 Oslo, Norway
| | - Jonas Paulsen
- Department of Biosciences, University of Oslo, 0316 Oslo, Norway
- Centre for Bioinformatics, Department of Informatics, University of Oslo, 0316 Oslo, Norway
| |
Collapse
|
13
|
Song T, Yao M, Yang Y, Liu Z, Zhang L, Li W. Integrative Identification by Hi-C Revealed Distinct Advanced Structural Variations in Lung Adenocarcinoma Tissue. PHENOMICS (CHAM, SWITZERLAND) 2023; 3:390-407. [PMID: 37589026 PMCID: PMC10425312 DOI: 10.1007/s43657-023-00103-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/16/2021] [Revised: 03/01/2023] [Accepted: 03/03/2023] [Indexed: 08/18/2023]
Abstract
Advanced three-dimensional structure variations of chromatin in large genome fragments, such as conversion of A/B compartment, topologically associated domains (TADs) and chromatin loops are related closely to occurrence of malignant tumors. However, the structural characteristics of lung cancer still remain uncovered. In this study, we used high-throughput chromosome (Hi-C) conformation capture technology to detect the advanced structural variations in chromatin of two non-smoking lung adenocarcinoma (LUAD) tumor and paired normal tissues. The results indicate that significant chromatin variations are detected in tumor tissues compared with normal tissues. At compartment scale, the main conversion type of compartment is A → B in tumor tissues, which are concentrated mainly on chromosome 3 (Chr3) (33.6%). A total of 216 tumor-specific TADs are identified in tumor tissues, which are distributed mainly in Chr1 (19), Chr2 (15) and Chr3 (17). Forty-one distinct enhancer-promoter loops are observed in tumor tissue, which are associated closely to tumor-related pathways including mitogen-activated protein kinase (MAPK), Phosphatidylinositol-3-kinase-Protein kinase B (PI3K-AKT), Ras, Wnt and Ras1. The most important observation in this study is that we identify five important genes (SYT16, NCEH1, NXPE3, MB21D2, and DZIP1L), which are detected in both A → B compartment, TADs and chromatin loops in tumor samples, and four of these genes (NCEH1, NXPE3, MB21D2, and DZIP1L) locate on q arm of Chr3. Further gene expression and invasion experiment analysis show that NCEH1, MB21D2 and SYT16 are involved in the tumor development. Thus, we provide a comprehensive overview of advanced structures in LUAD for the first time and provide a basis for further research on the genetic variation of this tumor.
Collapse
Affiliation(s)
- Tingting Song
- Department of Respiratory and Critical Care Medicine, Center of Precision Medicine, Precision Medicine Key Laboratory of Sichuan Province, Frontiers Science Center for Disease-Related Molecular Network, West China Hospital, Sichuan University, Chengdu, 610041 Sichuan Province China
| | - Menglin Yao
- Department of Respiratory and Critical Care Medicine, Center of Precision Medicine, Precision Medicine Key Laboratory of Sichuan Province, Frontiers Science Center for Disease-Related Molecular Network, West China Hospital, Sichuan University, Chengdu, 610041 Sichuan Province China
| | - Ying Yang
- Department of Respiratory and Critical Care Medicine, Center of Precision Medicine, Precision Medicine Key Laboratory of Sichuan Province, Frontiers Science Center for Disease-Related Molecular Network, West China Hospital, Sichuan University, Chengdu, 610041 Sichuan Province China
| | - Zhiqiang Liu
- Department of Respiratory and Critical Care Medicine, Center of Precision Medicine, Precision Medicine Key Laboratory of Sichuan Province, Frontiers Science Center for Disease-Related Molecular Network, West China Hospital, Sichuan University, Chengdu, 610041 Sichuan Province China
| | - Li Zhang
- Department of Respiratory and Critical Care Medicine, Center of Precision Medicine, Precision Medicine Key Laboratory of Sichuan Province, Frontiers Science Center for Disease-Related Molecular Network, West China Hospital, Sichuan University, Chengdu, 610041 Sichuan Province China
| | - Weimin Li
- Department of Respiratory and Critical Care Medicine, Center of Precision Medicine, Precision Medicine Key Laboratory of Sichuan Province, Frontiers Science Center for Disease-Related Molecular Network, West China Hospital, Sichuan University, Chengdu, 610041 Sichuan Province China
| |
Collapse
|
14
|
Dang D, Zhang SW, Duan R, Zhang S. Defining the separation landscape of topological domains for decoding consensus domain organization of the 3D genome. Genome Res 2023; 33:386-400. [PMID: 36894325 PMCID: PMC10078287 DOI: 10.1101/gr.277187.122] [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: 08/08/2022] [Accepted: 02/23/2023] [Indexed: 03/11/2023]
Abstract
Topologically associating domains (TADs) have emerged as basic structural and functional units of genome organization and have been determined by many computational methods from Hi-C contact maps. However, the TADs obtained by different methods vary greatly, which makes the accurate determination of TADs a challenging issue and hinders subsequent biological analyses about their organization and functions. Obvious inconsistencies among the TADs identified by different methods indeed make the statistical and biological properties of TADs overly depend on the chosen method rather than on the data. To this end, we use the consensus structural information captured by these methods to define the TAD separation landscape for decoding the consensus domain organization of the 3D genome. We show that the TAD separation landscape could be used to compare domain boundaries across multiple cell types for discovering conserved and divergent topological structures, decipher three types of boundary regions with diverse biological features, and identify consensus TADs (ConsTADs). We illustrate that these analyses could deepen our understanding of the relationships between the topological domains and chromatin states, gene expression, and DNA replication timing.
Collapse
Affiliation(s)
- Dachang Dang
- Key Laboratory of Information Fusion Technology of Ministry of Education, School of Automation, Northwestern Polytechnical University, Xi'an 710072, China
| | - Shao-Wu Zhang
- Key Laboratory of Information Fusion Technology of Ministry of Education, School of Automation, Northwestern Polytechnical University, Xi'an 710072, China;
| | - Ran Duan
- Department of Software Engineering, Yunnan University, Kunming 650500, China
| | - Shihua Zhang
- NCMIS, CEMS, RCSDS, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing 100190, China;
- School of Mathematical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
- Center for Excellence in Animal Evolution and Genetics, Chinese Academy of Sciences, Kunming 650223, China
- Key Laboratory of Systems Health Science of Zhejiang Province, School of Life Science, Hangzhou Institute for Advanced Study, University of Chinese Academy of Sciences, Hangzhou 310024, China
| |
Collapse
|
15
|
Rosen J, Lee L, Abnousi A, Chen J, Wen J, Hu M, Li Y. HPTAD: A computational method to identify topologically associating domains from HiChIP and PLAC-seq datasets. Comput Struct Biotechnol J 2023; 21:931-939. [PMID: 38213897 PMCID: PMC10782010 DOI: 10.1016/j.csbj.2023.01.003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2022] [Revised: 01/06/2023] [Accepted: 01/07/2023] [Indexed: 01/11/2023] Open
Abstract
High-throughput chromatin conformation capture technologies, such as Hi-C and Micro-C, have enabled genome-wide view of chromatin spatial organization. Most recently, Hi-C-derived enrichment-based technologies, including HiChIP and PLAC-seq, offer attractive alternatives due to their high signal-to-noise ratio and low cost. While a series of computational tools have been developed for Hi-C data, methods tailored for HiChIP and PLAC-seq data are still under development. Here we present HPTAD, a computational method to identify topologically associating domains (TADs) from HiChIP and PLAC-seq data. We performed comprehensive benchmark analysis to demonstrate its superior performance over existing TAD callers designed for Hi-C data. HPTAD is freely available at https://github.com/yunliUNC/HPTAD.
Collapse
Affiliation(s)
- Jonathan Rosen
- Department of Genetics, University of North Carolina, Chapel Hill, NC, USA
| | - Lindsay Lee
- Department of Quantitative Health Sciences, Lerner Research Institute, Cleveland Clinic Foundation, Cleveland, OH, USA
| | - Armen Abnousi
- Department of Quantitative Health Sciences, Lerner Research Institute, Cleveland Clinic Foundation, Cleveland, OH, USA
| | - Jiawen Chen
- Department of Biostatistics, University of North Carolina, Chapel Hill, NC, USA
| | - Jia Wen
- Department of Genetics, University of North Carolina, Chapel Hill, NC, USA
| | - Ming Hu
- Department of Quantitative Health Sciences, Lerner Research Institute, Cleveland Clinic Foundation, Cleveland, OH, USA
| | - Yun Li
- Department of Genetics, University of North Carolina, Chapel Hill, NC, USA
- Department of Biostatistics, University of North Carolina, Chapel Hill, NC, USA
- Department of Computer Science, University of North Carolina, Chapel Hill, NC, USA
| |
Collapse
|
16
|
Liu K, Li HD, Li Y, Wang J, Wang J. A Comparison of Topologically Associating Domain Callers Based on Hi-C Data. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2023; 20:15-29. [PMID: 35104223 DOI: 10.1109/tcbb.2022.3147805] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
Topologically associating domains (TADs) are local chromatin interaction domains, which have been shown to play an important role in gene expression regulation. TADs were originally discovered in the investigation of 3D genome organization based on High-throughput Chromosome Conformation Capture (Hi-C) data. Continuous considerable efforts have been dedicated to developing methods for detecting TADs from Hi-C data. Different computational methods for TADs identification vary in their assumptions and criteria in calling TADs. As a consequence, the TADs called by these methods differ in their similarities and biological features they are enriched in. In this work, we performed a systematic comparison of twenty-six TAD callers. We first compared the TADs and gaps between adjacent TADs across different methods, resolutions, and sequencing depths. We then assessed the quality of TADs and TAD boundaries according to three criteria: the decay of contact frequencies over the genomic distance, enrichment and depletion of regulatory elements around TAD boundaries, and reproducibility of TADs and TAD boundaries in replicate samples. Last, due to the lack of a gold standard of TADs, we also evaluated the performance of the methods on synthetic datasets. We discussed the key principles of TAD callers, and pinpointed current situation in the detection of TADs. We provide a concise, comprehensive, and systematic framework for evaluating the performance of TAD callers, and expect our work will provide useful guidance in choosing suitable approaches for the detection and evaluation of TADs.
Collapse
|
17
|
Reference panel guided topological structure annotation of Hi-C data. Nat Commun 2022; 13:7426. [PMID: 36460680 PMCID: PMC9718747 DOI: 10.1038/s41467-022-35231-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2022] [Accepted: 11/22/2022] [Indexed: 12/03/2022] Open
Abstract
Accurately annotating topological structures (e.g., loops and topologically associating domains) from Hi-C data is critical for understanding the role of 3D genome organization in gene regulation. This is a challenging task, especially at high resolution, in part due to the limited sequencing coverage of Hi-C data. Current approaches focus on the analysis of individual Hi-C data sets of interest, without taking advantage of the facts that (i) several hundred Hi-C contact maps are publicly available, and (ii) the vast majority of topological structures are conserved across multiple cell types. Here, we present RefHiC, an attention-based deep learning framework that uses a reference panel of Hi-C datasets to facilitate topological structure annotation from a given study sample. We compare RefHiC against tools that do not use reference samples and find that RefHiC outperforms other programs at both topological associating domain and loop annotation across different cell types, species, and sequencing depths.
Collapse
|
18
|
Wang M, Li J, Qi Z, Long Y, Pei L, Huang X, Grover CE, Du X, Xia C, Wang P, Liu Z, You J, Tian X, Ma Y, Wang R, Chen X, He X, Fang DD, Sun Y, Tu L, Jin S, Zhu L, Wendel JF, Zhang X. Genomic innovation and regulatory rewiring during evolution of the cotton genus Gossypium. Nat Genet 2022; 54:1959-1971. [PMID: 36474047 DOI: 10.1038/s41588-022-01237-2] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2022] [Accepted: 10/20/2022] [Indexed: 12/13/2022]
Abstract
Phenotypic diversity and evolutionary innovation ultimately trace to variation in genomic sequence and rewiring of regulatory networks. Here, we constructed a pan-genome of the Gossypium genus using ten representative diploid genomes. We document the genomic evolutionary history and the impact of lineage-specific transposon amplification on differential genome composition. The pan-3D genome reveals evolutionary connections between transposon-driven genome size variation and both higher-order chromatin structure reorganization and the rewiring of chromatin interactome. We linked changes in chromatin structures to phenotypic differences in cotton fiber and identified regulatory variations that decode the genetic basis of fiber length, the latter enabled by sequencing 1,005 transcriptomes during fiber development. We showcase how pan-genomic, pan-3D genomic and genetic regulatory data serve as a resource for delineating the evolutionary basis of spinnable cotton fiber. Our work provides insights into the evolution of genome organization and regulation and will inform cotton improvement by enabling regulome-based approaches.
Collapse
Affiliation(s)
- Maojun Wang
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Jianying Li
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Zhengyang Qi
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Yuexuan Long
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Liuling Pei
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Xianhui Huang
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Corrinne E Grover
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University, Ames, IA, USA
| | - Xiongming Du
- State Key Laboratory of Cotton Biology, Institute of Cotton Research, Chinese Academy of Agricultural Sciences, Anyang, China
| | - Chunjiao Xia
- College of Life Science and Technology, Huazhong Agricultural University, Wuhan, China
| | - Pengcheng Wang
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Zhenping Liu
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Jiaqi You
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Xuehan Tian
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Yizan Ma
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Ruipeng Wang
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Xinyuan Chen
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Xin He
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - David D Fang
- Cotton Fiber Bioscience Research Unit, USDA-ARS, Southern Regional Research Center, New Orleans, LA, USA
| | - Yuqiang Sun
- Zhejiang Sci-Tech University College of Life Sciences, Zhejiang, Hangzhou, China
| | - Lili Tu
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Shuangxia Jin
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Longfu Zhu
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China
| | - Jonathan F Wendel
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University, Ames, IA, USA.
| | - Xianlong Zhang
- National Key Laboratory of Crop Genetic Improvement, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, China.
| |
Collapse
|
19
|
Yu D, Chen G, Wang Y, Wang Y, Lin R, Liu N, Zhu P, Liu H, Hu T, Feng R, Feng H, Lan F, Cai J, Chen H. Regulation of cohesin-mediated chromosome folding by PDS5 in mammals. EMBO Rep 2022; 23:e54853. [PMID: 36129789 PMCID: PMC9638874 DOI: 10.15252/embr.202254853] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2022] [Revised: 09/04/2022] [Accepted: 09/08/2022] [Indexed: 09/23/2023] Open
Abstract
Cohesin regulates sister chromatid cohesion but also contributes to chromosome folding by promoting the formation of chromatin loops, a process mediated by loop extrusion. Although PDS5 regulates cohesin dynamics on chromatin, the exact function of PDS5 in cohesin-mediated chromatin looping remains unclear. Two paralogs of PDS5 exist in vertebrates, PDS5A and PDS5B. Here we show that PDS5A and PDS5B co-localize with RAD21 and CTCF at loop anchors. Rapid PDS5A or PDS5B degradation in liver cancer cells using an inducible degron system reduces chromatin loops and increases loop size. RAD21 enrichment at loop anchors is decreased upon depletion of PDS5A or PDS5B. PDS5B loss also reduces CTCF signals at loop anchors and has a stronger effect on loop enlargement compared with PDS5A. Co-depletion of PDS5A and PDS5B reduces RAD21 levels at loop anchors although the amount of cohesin on chromatin is increased. Our study provides insight into how PDS5 proteins regulate cohesin-mediated chromatin looping.
Collapse
Affiliation(s)
- Dingdang Yu
- Key Laboratory of Medical Epigenetics and Metabolism, Institutes of Biomedical SciencesFudan UniversityShanghaiChina
- Department of Human Cell Biology and Genetics, School of MedicineSouthern University of Science and TechnologyShenzhenChina
| | - Guoyu Chen
- State Key Laboratory of Oncogenes and Related Genes, Shanghai Cancer InstituteRenji HospitalShanghaiChina
- Renji‐Med X Clinical Stem Cell Research Center, Ren Ji HospitalShanghai Jiao Tong University School of MedicineShanghaiChina
| | - Yuci Wang
- Department of Human Cell Biology and Genetics, School of MedicineSouthern University of Science and TechnologyShenzhenChina
| | - Yining Wang
- Key Laboratory of Carcinogenesis and Cancer Invasion, Ministry of Education, Liver Cancer Institute, Zhongshan HospitalFudan UniversityShanghaiChina
| | - Risheng Lin
- Department of Human Cell Biology and Genetics, School of MedicineSouthern University of Science and TechnologyShenzhenChina
| | - Nanbo Liu
- Department of Cardiac Surgery, and Department of Medical Sciences, Guangdong Cardiovascular Institute, Guangdong Provincial People's HospitalGuangdong Academy of Medical SciencesGuangzhouChina
| | - Ping Zhu
- Department of Cardiac Surgery, and Department of Medical Sciences, Guangdong Cardiovascular Institute, Guangdong Provincial People's HospitalGuangdong Academy of Medical SciencesGuangzhouChina
| | - Hang Liu
- Key Laboratory of Medical Epigenetics and Metabolism, Institutes of Biomedical SciencesFudan UniversityShanghaiChina
| | - Tao Hu
- Key Laboratory of Medical Epigenetics and Metabolism, Institutes of Biomedical SciencesFudan UniversityShanghaiChina
| | - Rui Feng
- Department of Human Cell Biology and Genetics, School of MedicineSouthern University of Science and TechnologyShenzhenChina
| | - Haizhong Feng
- State Key Laboratory of Oncogenes and Related Genes, Shanghai Cancer InstituteRenji HospitalShanghaiChina
- Renji‐Med X Clinical Stem Cell Research Center, Ren Ji HospitalShanghai Jiao Tong University School of MedicineShanghaiChina
| | - Fei Lan
- Key Laboratory of Medical Epigenetics and Metabolism, Institutes of Biomedical SciencesFudan UniversityShanghaiChina
- Key Laboratory of Carcinogenesis and Cancer Invasion, Ministry of Education, Liver Cancer Institute, Zhongshan HospitalFudan UniversityShanghaiChina
| | - Jiabin Cai
- Key Laboratory of Medical Epigenetics and Metabolism, Institutes of Biomedical SciencesFudan UniversityShanghaiChina
- Key Laboratory of Carcinogenesis and Cancer Invasion, Ministry of Education, Liver Cancer Institute, Zhongshan HospitalFudan UniversityShanghaiChina
| | - Hao Chen
- Department of Human Cell Biology and Genetics, School of MedicineSouthern University of Science and TechnologyShenzhenChina
| |
Collapse
|
20
|
3D chromatin remodeling potentiates transcriptional programs driving cell invasion. Proc Natl Acad Sci U S A 2022; 119:e2203452119. [PMID: 36037342 PMCID: PMC9457068 DOI: 10.1073/pnas.2203452119] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The contribution of deregulated chromatin architecture, including topologically associated domains (TADs), to cancer progression remains ambiguous. CCCTC-binding factor (CTCF) is a central regulator of higher-order chromatin structure that undergoes copy number loss in over half of all breast cancers, but the impact of this defect on epigenetic programming and chromatin architecture remains unclear. We find that under physiological conditions, CTCF organizes subTADs to limit the expression of oncogenic pathways, including phosphatidylinositol 3-kinase (PI3K) and cell adhesion networks. Loss of a single CTCF allele potentiates cell invasion through compromised chromatin insulation and a reorganization of chromatin architecture and histone programming that facilitates de novo promoter-enhancer contacts. However, this change in the higher-order chromatin landscape leads to a vulnerability to inhibitors of mTOR. These data support a model whereby subTAD reorganization drives both modification of histones at de novo enhancer-promoter contacts and transcriptional up-regulation of oncogenic transcriptional networks.
Collapse
|
21
|
Su C, Gao L, May CL, Pippin JA, Boehm K, Lee M, Liu C, Pahl MC, Golson ML, Naji A, Grant SFA, Wells AD, Kaestner KH. 3D chromatin maps of the human pancreas reveal lineage-specific regulatory architecture of T2D risk. Cell Metab 2022; 34:1394-1409.e4. [PMID: 36070683 PMCID: PMC9664375 DOI: 10.1016/j.cmet.2022.08.014] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/23/2022] [Revised: 06/03/2022] [Accepted: 08/17/2022] [Indexed: 12/20/2022]
Abstract
Three-dimensional (3D) chromatin organization maps help dissect cell-type-specific gene regulatory programs. Furthermore, 3D chromatin maps contribute to elucidating the pathogenesis of complex genetic diseases by connecting distal regulatory regions and genetic risk variants to their respective target genes. To understand the cell-type-specific regulatory architecture of diabetes risk, we generated transcriptomic and 3D epigenomic profiles of human pancreatic acinar, alpha, and beta cells using single-cell RNA-seq, single-cell ATAC-seq, and high-resolution Hi-C of sorted cells. Comparisons of these profiles revealed differential A/B (open/closed) chromatin compartmentalization, chromatin looping, and transcriptional factor-mediated control of cell-type-specific gene regulatory programs. We identified a total of 4,750 putative causal-variant-to-target-gene pairs at 194 type 2 diabetes GWAS signals using pancreatic 3D chromatin maps. We found that the connections between candidate causal variants and their putative target effector genes are cell-type stratified and emphasize previously underappreciated roles for alpha and acinar cells in diabetes pathogenesis.
Collapse
Affiliation(s)
- Chun Su
- Division of Human Genetics and Endocrinology & Diabetes, The Children's Hospital of Philadelphia, Philadelphia, PA, USA; Center for Spatial and Functional Genomics, The Children's Hospital of Philadelphia, Philadelphia, PA, USA
| | - Long Gao
- Institute for Diabetes, Obesity, and Metabolism, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, USA; Department of Genetics, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, USA
| | - Catherine L May
- Institute for Diabetes, Obesity, and Metabolism, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, USA; Department of Genetics, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, USA
| | - James A Pippin
- Division of Human Genetics and Endocrinology & Diabetes, The Children's Hospital of Philadelphia, Philadelphia, PA, USA; Center for Spatial and Functional Genomics, The Children's Hospital of Philadelphia, Philadelphia, PA, USA
| | - Keith Boehm
- Division of Human Genetics and Endocrinology & Diabetes, The Children's Hospital of Philadelphia, Philadelphia, PA, USA; Center for Spatial and Functional Genomics, The Children's Hospital of Philadelphia, Philadelphia, PA, USA
| | - Michelle Lee
- Institute for Diabetes, Obesity, and Metabolism, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, USA; Department of Genetics, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, USA
| | - Chengyang Liu
- Institute for Diabetes, Obesity, and Metabolism, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, USA; Department of Surgery, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, USA
| | - Matthew C Pahl
- Division of Human Genetics and Endocrinology & Diabetes, The Children's Hospital of Philadelphia, Philadelphia, PA, USA; Center for Spatial and Functional Genomics, The Children's Hospital of Philadelphia, Philadelphia, PA, USA
| | - Maria L Golson
- Institute for Diabetes, Obesity, and Metabolism, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, USA; Department of Genetics, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, USA
| | - Ali Naji
- Institute for Diabetes, Obesity, and Metabolism, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, USA; Department of Surgery, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, USA
| | - Struan F A Grant
- Division of Human Genetics and Endocrinology & Diabetes, The Children's Hospital of Philadelphia, Philadelphia, PA, USA; Center for Spatial and Functional Genomics, The Children's Hospital of Philadelphia, Philadelphia, PA, USA; Institute for Diabetes, Obesity, and Metabolism, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, USA; Department of Genetics, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, USA; Department of Pediatrics, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, USA.
| | - Andrew D Wells
- Center for Spatial and Functional Genomics, The Children's Hospital of Philadelphia, Philadelphia, PA, USA; Department of Pathology and Laboratory Medicine, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, USA; Institute for Immunology, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, USA.
| | - Klaus H Kaestner
- Institute for Diabetes, Obesity, and Metabolism, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, USA; Department of Genetics, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, USA.
| |
Collapse
|
22
|
Yu W, Zhong Q, Wen Z, Zhang W, Huang Y. Genome architecture plasticity underlies DNA replication timing dynamics in cell differentiation. Front Genet 2022; 13:961612. [PMID: 36118849 PMCID: PMC9478753 DOI: 10.3389/fgene.2022.961612] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2022] [Accepted: 07/15/2022] [Indexed: 12/04/2022] Open
Abstract
During the S-phase of eukaryotic cell cycle, DNA is replicated in a dedicatedly regulated temporal order, with regions containing active and inactive genes replicated early and late, respectively. Recent advances in sequencing technology allow us to explore the connection between replication timing (RT), histone modifications, and three-dimensional (3D) chromatin structure in diverse cell types. To characterize the dynamics during cell differentiation, corresponding sequencing data for human embryonic stem cells and four differentiated cell types were collected. By comparing RT and its extent of conservation before and after germ layer specification, the human genome was partitioned into distinct categories. Each category is then subject to comparisons on genomic, epigenetic, and chromatin 3D structural features. As expected, while constitutive early and late replication regions showed active and inactive features, respectively, dynamic regions with switched RT showed intermediate features. Surprisingly, although early-to-late replication and late-to-early replication regions showed similar histone modification patterns in hESCs, their structural preferences were opposite. Specifically, in hESCs, early-to-late replication regions tended to appear in the B compartment and large topologically associated domains, while late-to-early replication regions showed the opposite. Our results uncover the coordinated regulation of RT and 3D genome structure that underlies the loss of pluripotency and lineage commitment and indicate the importance and potential roles of genome architecture in biological processes.
Collapse
Affiliation(s)
- Wenjun Yu
- Center for Genetics and Developmental Systems Biology, Department of Obstetrics and Gynecology, Nanfang Hospital, Southern Medical University, Guangzhou, China
- *Correspondence: Wenjun Yu,
| | - Quan Zhong
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, China
| | - Zi Wen
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, China
| | - Weihan Zhang
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, China
| | - Yanrong Huang
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, China
| |
Collapse
|
23
|
Lebeau B, Zhao K, Jangal M, Zhao T, Guerra M, Greenwood CMT, Witcher M. Single base-pair resolution analysis of DNA binding motif with MoMotif reveals an oncogenic function of CTCF zinc-finger 1 mutation. Nucleic Acids Res 2022; 50:8441-8458. [PMID: 35947648 PMCID: PMC9410893 DOI: 10.1093/nar/gkac658] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2022] [Accepted: 07/21/2022] [Indexed: 12/24/2022] Open
Abstract
Defining the impact of missense mutations on the recognition of DNA motifs is highly dependent on bioinformatic tools that define DNA binding elements. However, classical motif analysis tools remain limited in their capacity to identify subtle changes in complex binding motifs between distinct conditions. To overcome this limitation, we developed a new tool, MoMotif, that facilitates a sensitive identification, at the single base-pair resolution, of complex, or subtle, alterations to core binding motifs, discerned from ChIP-seq data. We employed MoMotif to define the previously uncharacterized recognition motif of CTCF zinc-finger 1 (ZF1), and to further define the impact of CTCF ZF1 mutation on its association with chromatin. Mutations of CTCF ZF1 are exclusive to breast cancer and are associated with metastasis and therapeutic resistance, but the underlying mechanisms are unclear. Using MoMotif, we identified an extension of the CTCF core binding motif, necessitating a functional ZF1 to bind appropriately. Using a combination of ChIP-Seq and RNA-Seq, we discover that the inability to bind this extended motif drives an altered transcriptional program associated with the oncogenic phenotypes observed clinically. Our study demonstrates that MoMotif is a powerful new tool for comparative ChIP-seq analysis and characterising DNA-protein contacts.
Collapse
Affiliation(s)
| | | | - Maika Jangal
- Lady Davis Institute, Jewish General Hospital, Montréal, Québec H3T 1E2, Canada
| | - Tiejun Zhao
- Lady Davis Institute, Jewish General Hospital, Montréal, Québec H3T 1E2, Canada
| | - Maria Guerra
- Lady Davis Institute, Jewish General Hospital, Montréal, Québec H3T 1E2, Canada
| | - Celia M T Greenwood
- Correspondence may also be addressed to Celia Greenwood. Tel: +1 514 340 8222 (Ext 28397);
| | - Michael Witcher
- To whom correspondence should be addressed. Tel: +1 514 340 8222 (Ext 23363);
| |
Collapse
|
24
|
Liu E, Lyu H, Peng Q, Liu Y, Wang T, Han J. TADfit is a multivariate linear regression model for profiling hierarchical chromatin domains on replicate Hi-C data. Commun Biol 2022; 5:608. [PMID: 35725901 PMCID: PMC9209495 DOI: 10.1038/s42003-022-03546-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2021] [Accepted: 05/31/2022] [Indexed: 11/26/2022] Open
Abstract
Topologically associating domains (TADs) are fundamental building blocks of three dimensional genome, and organized into complex hierarchies. Identifying hierarchical TADs on Hi-C data helps to understand the relationship between genome architectures and gene regulation. Herein we propose TADfit, a multivariate linear regression model for profiling hierarchical chromatin domains, which tries to fit the interaction frequencies in Hi-C contact matrix with and without replicates using all-possible hierarchical TADs, and the significant ones can be determined by the regression coefficients obtained with the help of an online learning solver called Follow-The-Regularized-Leader (FTRL). Beyond the existing methods, TADfit has an ability to handle multiple contact matrix replicates and find partially overlapping TADs on them, which helps to find the comprehensive underlying TADs across replicates from different experiments. The comparative results tell that TADfit has better accuracy and reproducibility, and the hierarchical TADs called by it exhibit a reasonable biological relevance. TADfit is a computational method that can identify hierarchical or partially overlapping TADs from Hi-C data, in part by using information from multiple replicates to improve detection power.
Collapse
Affiliation(s)
- Erhu Liu
- School of Automation Science and Engineering, Faculty of Electronic and Information Engineering, Xi'an Jiaotong University, Shaanxi, 710049, China
| | - Hongqiang Lyu
- School of Automation Science and Engineering, Faculty of Electronic and Information Engineering, Xi'an Jiaotong University, Shaanxi, 710049, China. .,Guangdong Artificial Intelligence and Digital Economy Laboratory, Guangdong, 510335, China.
| | - Qinke Peng
- School of Automation Science and Engineering, Faculty of Electronic and Information Engineering, Xi'an Jiaotong University, Shaanxi, 710049, China
| | - Yuan Liu
- School of Automation Science and Engineering, Faculty of Electronic and Information Engineering, Xi'an Jiaotong University, Shaanxi, 710049, China
| | - Tian Wang
- Institute of Artificial Intelligence, Beihang University, Beijing, 100191, China
| | - Jiuqiang Han
- School of Automation Science and Engineering, Faculty of Electronic and Information Engineering, Xi'an Jiaotong University, Shaanxi, 710049, China.,Guangdong Artificial Intelligence and Digital Economy Laboratory, Guangdong, 510335, China
| |
Collapse
|
25
|
Wang X, Luan Y, Yue F. EagleC: A deep-learning framework for detecting a full range of structural variations from bulk and single-cell contact maps. SCIENCE ADVANCES 2022; 8:eabn9215. [PMID: 35704579 PMCID: PMC9200291 DOI: 10.1126/sciadv.abn9215] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/31/2021] [Accepted: 04/28/2022] [Indexed: 05/11/2023]
Abstract
The Hi-C technique has been shown to be a promising method to detect structural variations (SVs) in human genomes. However, algorithms that can use Hi-C data for a full-range SV detection have been severely lacking. Current methods can only identify interchromosomal translocations and long-range intrachromosomal SVs (>1 Mb) at less-than-optimal resolution. Therefore, we develop EagleC, a framework that combines deep-learning and ensemble-learning strategies to predict a full range of SVs at high resolution. We show that EagleC can uniquely capture a set of fusion genes that are missed by whole-genome sequencing or nanopore. Furthermore, EagleC also effectively captures SVs in other chromatin interaction platforms, such as HiChIP, Chromatin interaction analysis with paired-end tag sequencing (ChIA-PET), and capture Hi-C. We apply EagleC in more than 100 cancer cell lines and primary tumors and identify a valuable set of high-quality SVs. Last, we demonstrate that EagleC can be applied to single-cell Hi-C and used to study the SV heterogeneity in primary tumors.
Collapse
Affiliation(s)
- Xiaotao Wang
- Department of Biochemistry and Molecular Genetics, Feinberg School of Medicine Northwestern University, Chicago, IL, USA
| | - Yu Luan
- Department of Biochemistry and Molecular Genetics, Feinberg School of Medicine Northwestern University, Chicago, IL, USA
| | - Feng Yue
- Department of Biochemistry and Molecular Genetics, Feinberg School of Medicine Northwestern University, Chicago, IL, USA
- Robert H. Lurie Comprehensive Cancer Center of Northwestern University, Chicago, IL, USA
| |
Collapse
|
26
|
Sefer E. A comparison of topologically associating domain callers over mammals at high resolution. BMC Bioinformatics 2022; 23:127. [PMID: 35413815 PMCID: PMC9006547 DOI: 10.1186/s12859-022-04674-2] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2021] [Accepted: 04/07/2022] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Topologically associating domains (TADs) are locally highly-interacting genome regions, which also play a critical role in regulating gene expression in the cell. TADs have been first identified while investigating the 3D genome structure over High-throughput Chromosome Conformation Capture (Hi-C) interaction dataset. Substantial degree of efforts have been devoted to develop techniques for inferring TADs from Hi-C interaction dataset. Many TAD-calling methods have been developed which differ in their criteria and assumptions in TAD inference. Correspondingly, TADs inferred via these callers vary in terms of both similarities and biological features they are enriched in. RESULT We have carried out a systematic comparison of 27 TAD-calling methods over mammals. We use Micro-C, a recent high-resolution variant of Hi-C, to compare TADs at a very high resolution, and classify the methods into 3 categories: feature-based methods, Clustering methods, Graph-partitioning methods. We have evaluated TAD boundaries, gaps between adjacent TADs, and quality of TADs across various criteria. We also found particularly CTCF and Cohesin proteins to be effective in formation of TADs with corner dots. We have also assessed the callers performance on simulated datasets since a gold standard for TADs is missing. TAD sizes and numbers change remarkably between TAD callers and dataset resolutions, indicating that TADs are hierarchically-organized domains, instead of disjoint regions. A core subset of feature-based TAD callers regularly perform the best while inferring reproducible domains, which are also enriched for TAD related biological properties. CONCLUSION We have analyzed the fundamental principles of TAD-calling methods, and identified the existing situation in TAD inference across high resolution Micro-C interaction datasets over mammals. We come up with a systematic, comprehensive, and concise framework to evaluate the TAD-calling methods performance across Micro-C datasets. Our research will be useful in selecting appropriate methods for TAD inference and evaluation based on available data, experimental design, and biological question of interest. We also introduce our analysis as a benchmarking tool with publicly available source code.
Collapse
Affiliation(s)
- Emre Sefer
- Department of Computer Science, Ozyegin University, Istanbul, Turkey.
| |
Collapse
|
27
|
Mapping nucleosome and chromatin architectures: A survey of computational methods. Comput Struct Biotechnol J 2022; 20:3955-3962. [PMID: 35950186 PMCID: PMC9340519 DOI: 10.1016/j.csbj.2022.07.037] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2022] [Revised: 07/22/2022] [Accepted: 07/22/2022] [Indexed: 11/21/2022] Open
Abstract
With ever-growing genomic sequencing data, the data variabilities and the underlying biases of the sequencing technologies pose significant computational challenges ranging from the need for accurately detecting the nucleosome positioning or chromatin interaction to the need for developing normalization methods to eliminate systematic biases. This review mainly surveys the computational methods for mapping the higher-resolution nucleosome and higher-order chromatin architectures. While a detailed discussion of the underlying algorithms is beyond the scope of our survey, we have discussed the methods and tools that can detect the nucleosomes in the genome, then demonstrated the computational methods for identifying 3D chromatin domains and interactions. We further illustrated computational approaches for integrating multi-omics data with Hi-C data and the advance of single-cell (sc)Hi-C data analysis. Our survey provides a comprehensive and valuable resource for biomedical scientists interested in studying nucleosome organization and chromatin structures as well as for computational scientists who are interested in improving upon them.
Collapse
|
28
|
Xu L, Yin L, Qi Y, Tan X, Gao M, Peng J. 3D disorganization and rearrangement of genome provide insights into pathogenesis of NAFLD by integrated Hi-C, Nanopore, and RNA sequencing. Acta Pharm Sin B 2021; 11:3150-3164. [PMID: 34729306 PMCID: PMC8546856 DOI: 10.1016/j.apsb.2021.03.022] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2020] [Revised: 01/30/2021] [Accepted: 02/07/2021] [Indexed: 12/12/2022] Open
Abstract
The three-dimensional (3D) conformation of chromatin is integral to the precise regulation of gene expression. The 3D genome and genomic variations in non-alcoholic fatty liver disease (NAFLD) are largely unknown, despite their key roles in cellular function and physiological processes. High-throughput chromosome conformation capture (Hi-C), Nanopore sequencing, and RNA-sequencing (RNA-seq) assays were performed on the liver of normal and NAFLD mice. A high-resolution 3D chromatin interaction map was generated to examine different 3D genome hierarchies including A/B compartments, topologically associated domains (TADs), and chromatin loops by Hi-C, and whole genome sequencing identifying structural variations (SVs) and copy number variations (CNVs) by Nanopore sequencing. We identified variations in thousands of regions across the genome with respect to 3D chromatin organization and genomic rearrangements, between normal and NAFLD mice, and revealed gene dysregulation frequently accompanied by these variations. Candidate target genes were identified in NAFLD, impacted by genetic rearrangements and spatial organization disruption. Our data provide a high-resolution 3D genome interaction resource for NAFLD investigations, revealed the relationship among genetic rearrangements, spatial organization disruption, and gene regulation, and identified candidate genes associated with these variations implicated in the pathogenesis of NAFLD. The newly findings offer insights into novel mechanisms of NAFLD pathogenesis and can provide a new conceptual framework for NAFLD therapy.
Collapse
Key Words
- 3C, chromosome conformation capture
- 3D genome
- 3D, three-dimensional
- ALT, alanine aminotransferase
- AST, aspartate aminotransferase
- Abcg5, ATP-binding cassette sub-family G member 5
- BWA, Burrows-Wheeler Aligner
- CNV, copy number variation
- Camk1d, calcium/calmodulin-dependent protein kinase type 1D
- Chr, chromosome
- Chromatin looping
- DEG, differentially expressed gene
- DEL, deletion
- DI, directionality index
- DUP, duplication
- Elovl6, elongation of very long chain fatty acids protein 6
- FDR, false discovery rate
- FFA, free fatty acid
- Fgfr2, fibroblast growth factor receptor 2
- GCKR, glucokinase regulator
- GO, gene ontology
- GSH, glutathione
- Gadd45g, growth arrest and DNA damage-inducible protein GADD45 gamma
- Grm8, metabotropic glutamate receptor 8
- Gsta1, glutathione S-transferase A1
- H&E, hematoxylin-eosin
- HFD, high-fat diet
- HSD17B13, hydroxysteroid 17-beta dehydrogenase 13
- Hi-C, high-throughput chromosome conformation capture
- IDE, interaction decay exponent
- INS, insertion
- INV, inversion
- IR, inclusion ratio
- IRGM, immunity related GTPase M
- IRS4, insulin receptor substrate 4
- KEGG, Kyoto Encyclopedia of Genes and Genomes
- Kcnma1, calcium-activated potassium channel subunit alpha-1
- LPIN1, lipin 1
- MBOAT7, membrane bound O-acyltransferase domain containing 7
- MDA, malondialdehyde
- NAFLD, non-alcoholic fatty liver disease
- NF1, neurofibromin 1
- NGS, next-generation sequencing
- NOTCH1, notch receptor 1
- ONT, Oxford Nanopore Technologies
- PCA, principal component analysis
- PNPLA3, patatin like phospholipase domain containing 3
- PPP1R3B, protein phosphatase 1 regulatory subunit 3B
- PTEN, phosphatase and tensin homolog
- Pde4b, phosphodiesterase 4B
- Plce1, 1-phosphat-idylinositol 4,5-bisphosphate phosphodiesterase epsilon-1
- Plxnb1, Plexin-B1
- RB1, RB transcriptional corepressor 1
- RNA-seq, RNA-sequencing
- SD, standard deviation
- SOD, superoxide dismutase
- SV, structural variation
- Scd1, acyl-CoA desaturase 1
- Sugct, succinate-hydroxymethylglutarate CoA-transferase
- TAD, topologically associated domain
- TC, total cholesterol
- TG, triglyceride
- TM6SF2, transmembrane 6 superfamily member 2
- TP53, tumor protein p53
- TRA, translocation
- Topologically associated domain
- Transcriptome
- WGS, whole-genome sequencing
- Whole-genome sequencing
Collapse
|
29
|
Wang M, Li J, Wang P, Liu F, Liu Z, Zhao G, Xu Z, Pei L, Grover CE, Wendel JF, Wang K, Zhang X. Comparative Genome Analyses Highlight Transposon-Mediated Genome Expansion and the Evolutionary Architecture of 3D Genomic Folding in Cotton. Mol Biol Evol 2021. [PMID: 33973633 DOI: 10.21203/rs.3.rs-93594/v1] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/08/2023] Open
Abstract
Transposable element (TE) amplification has been recognized as a driving force mediating genome size expansion and evolution, but the consequences for shaping 3D genomic architecture remains largely unknown in plants. Here, we report reference-grade genome assemblies for three species of cotton ranging 3-fold in genome size, namely Gossypium rotundifolium (K2), G. arboreum (A2), and G. raimondii (D5), using Oxford Nanopore Technologies. Comparative genome analyses document the details of lineage-specific TE amplification contributing to the large genome size differences (K2, 2.44 Gb; A2, 1.62 Gb; D5, 750.19 Mb) and indicate relatively conserved gene content and synteny relationships among genomes. We found that approximately 17% of syntenic genes exhibit chromatin status change between active ("A") and inactive ("B") compartments, and TE amplification was associated with the increase of the proportion of A compartment in gene regions (∼7,000 genes) in K2 and A2 relative to D5. Only 42% of topologically associating domain (TAD) boundaries were conserved among the three genomes. Our data implicate recent amplification of TEs following the formation of lineage-specific TAD boundaries. This study sheds light on the role of transposon-mediated genome expansion in the evolution of higher-order chromatin structure in plants.
Collapse
Affiliation(s)
- Maojun Wang
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan, Hubei, China
| | - Jianying Li
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan, Hubei, China
| | - Pengcheng Wang
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan, Hubei, China
| | - Fang Liu
- State Key Laboratory of Cotton Biology, Institute of Cotton Research, Chinese Academy of Agricultural Sciences, Anyang, Henan, China
| | - Zhenping Liu
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan, Hubei, China
| | - Guannan Zhao
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan, Hubei, China
| | - Zhongping Xu
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan, Hubei, China
| | - Liuling Pei
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan, Hubei, China
| | - Corrinne E Grover
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University, Ames, USA
| | - Jonathan F Wendel
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University, Ames, USA
| | - Kunbo Wang
- State Key Laboratory of Cotton Biology, Institute of Cotton Research, Chinese Academy of Agricultural Sciences, Anyang, Henan, China
| | - Xianlong Zhang
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan, Hubei, China
| |
Collapse
|
30
|
Long Y, Liu Z, Wang P, Yang H, Wang Y, Zhang S, Zhang X, Wang M. Disruption of topologically associating domains by structural variations in tetraploid cottons. Genomics 2021; 113:3405-3414. [PMID: 34311045 DOI: 10.1016/j.ygeno.2021.07.023] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2021] [Revised: 07/12/2021] [Accepted: 07/21/2021] [Indexed: 01/18/2023]
Abstract
Structural variations (SVs) are recognized to have an important role in transcriptional regulation, especially in the light of resolved 3D genome structure using high-throughput chromosome conformation capture (Hi-C) technology in mammals. However, the effect of SVs on 3D genome organization in plants remains rarely understood. In this study, we identified 295,496 SVs and 5251 topologically associating domains (TADs) in two diploid and two tetraploid cottons. We observed that approximately 16% of SVs occurred in TAD boundary regions that were called boundary affecting-structural variations (BA-SVs), and had a large effect on disrupting TAD organization. Nevertheless, SVs preferred occurring in TAD interior instead of TAD boundary, probably associated with the relaxed evolutionary selection pressure. We noticed the biased evolution of the At and Dt subgenomes of tetraploid cottons, in terms of SV-mediated disruption of 3D genome structure relative to diploids. In addition, we provide evidence showing that both SVs and TAD disruption could lead to expression difference of orthologous genes. This study advances our understanding of the effect of SVs on 3D genome organization and gene expression regulation in plants.
Collapse
Affiliation(s)
- Yuexuan Long
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan 430070, Hubei, China
| | - Zhenping Liu
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan 430070, Hubei, China
| | - Pengcheng Wang
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan 430070, Hubei, China
| | - Hang Yang
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan 430070, Hubei, China
| | - Yuejin Wang
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan 430070, Hubei, China
| | - Sainan Zhang
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan 430070, Hubei, China
| | - Xianlong Zhang
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan 430070, Hubei, China
| | - Maojun Wang
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan 430070, Hubei, China.
| |
Collapse
|
31
|
Meng T, Huang R, Jin J, Gao J, Liu F, Wei Z, Xu X, Chang Z, Lin J, Ta N, Huang Z, Yin H, Zhou W, Song D. The comparative integrated multi-omics analysis identifies CA2 as a novel target for chordoma. Neuro Oncol 2021; 23:1709-1722. [PMID: 34214167 DOI: 10.1093/neuonc/noab156] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
BACKGROUND Chordoma is a rare mesenchymal malignancy, with a high recurrence rate and unclear tumorigenic mechanism. Genetic alterations, epigenetic regulators, and chromatin spatial organization play crucial roles in the initiation and progression of chordoma. In the current study, we aim to uncover the novel therapeutical targets for chordoma via using integrated multi-omics analysis. METHODS The RNA-sequencing (RNA-seq), assay for transposable accessible chromatin by high throughput sequencing (ATAC-seq) and Hi-C were performed between chordoma and human nucleus pulposus (HNP), along with imageological examination and clinical information. The expressions of identified targets were validated by clinical samples and their function were further evaluated by cell and animal experiments via gene knockdown and inhibitors. RESULTS The integrated multi-omics analysis revealed the important roles of bone microenvironment in chordoma tumorigenesis. By comparing the hierarchical structures, CA2 and THNSL2 were identified in the switched compartments, cell-specific boundaries and loops. Additionally, CA2 was highly expressed in chordoma, but barely found in HNP. The cell growth and migration of chordoma cells were dramatically suppressed via inhibition of CA2 either with genetic deletion or pharmaceutical treatment with Dorzolamide HCl. Furthermore, Dorzolamide HCl also regulated the bone microenvironment by blocking the osteoclast differentiation of bone marrow monocytes. CONCLUSION This study uncovers the roles of bone microenvironment in the chordoma tumorigenesis and identifies CA2 as a novel therapeutic target for chordoma. Besides, our findings suggest Dorzolamide HCl as a promising therapeutic option for chordoma.
Collapse
Affiliation(s)
- Tong Meng
- Department of Orthopedics, Shanghai General Hospital, School of Medicine, Shanghai Jiaotong University, Shanghai, China.,Tongji University Cancer Center, Shanghai Tenth People's Hospital, School of Medicine, Tongji University, Shanghai, China
| | - Runzhi Huang
- Department of Orthopedics, Tongji Hospital, Tongji University School of Medicine, Tongji University, Shanghai, China.,Department of Orthopedics, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
| | - Jiali Jin
- Tongji University Cancer Center, Shanghai Tenth People's Hospital, School of Medicine, Tongji University, Shanghai, China
| | - Jianxuan Gao
- Department of Orthopedics, Shanghai General Hospital, School of Medicine, Shanghai Jiaotong University, Shanghai, China
| | - Fuyan Liu
- Biomarker Technologies Corporation, Beijing, China
| | - Ziheng Wei
- Department of Orthopedics, Shanghai General Hospital, School of Medicine, Shanghai Jiaotong University, Shanghai, China
| | - Xiaowen Xu
- Department of Medical Imaging, Tongji Hospital, Tongji University School of Medicine, Tongji University, Shanghai, China
| | - Zhengyan Chang
- Department of Pathology, Shanghai Tenth People's Hospital, Tongji University School of Medicine, Shanghai, China
| | - Jun Lin
- Department of Pathology, Shanghai General Hospital, School of Medicine, Shanghai Jiaotong University, Shanghai, China
| | - Na Ta
- Department of Pathology, Shanghai Changhai Hospital, Navy Medical University, Shanghai, China
| | - Zongqiang Huang
- Department of Orthopedics, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
| | - Huabin Yin
- Department of Orthopedics, Shanghai General Hospital, School of Medicine, Shanghai Jiaotong University, Shanghai, China
| | - Wang Zhou
- Departments of Neurovascular Center, Shanghai Changhai Hospital, Navy Medical University, Shanghai, China.,The Musculoskeletal laboratory, Institute of Biotechnology, University of Shanghai for Science and Technology, Shanghai, China
| | - Dianwen Song
- Department of Orthopedics, Shanghai General Hospital, School of Medicine, Shanghai Jiaotong University, Shanghai, China
| |
Collapse
|
32
|
Wang X, Xu J, Zhang B, Hou Y, Song F, Lyu H, Yue F. Genome-wide detection of enhancer-hijacking events from chromatin interaction data in rearranged genomes. Nat Methods 2021; 18:661-668. [PMID: 34092790 PMCID: PMC8191102 DOI: 10.1038/s41592-021-01164-w] [Citation(s) in RCA: 68] [Impact Index Per Article: 22.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2020] [Accepted: 04/22/2021] [Indexed: 02/06/2023]
Abstract
Recent efforts have shown that structural variations (SVs) can disrupt three-dimensional genome organization and induce enhancer hijacking, yet no computational tools exist to identify such events from chromatin interaction data. Here, we develop NeoLoopFinder, a computational framework to identify the chromatin interactions induced by SVs, including interchromosomal translocations, large deletions and inversions. Our framework can automatically resolve complex SVs, reconstruct local Hi-C maps surrounding the breakpoints, normalize copy number variation and allele effects and predict chromatin loops induced by SVs. We applied NeoLoopFinder in Hi-C data from 50 cancer cell lines and primary tumors and identified tens of recurrent genes associated with enhancer hijacking. To experimentally validate NeoLoopFinder, we deleted the hijacked enhancers in prostate adenocarcinoma cells using CRISPR-Cas9, which significantly reduced expression of the target oncogene. In summary, NeoLoopFinder enables identification of critical oncogenic regulatory elements that can potentially reveal therapeutic targets.
Collapse
Affiliation(s)
- Xiaotao Wang
- Department of Biochemistry and Molecular Genetics, Feinberg School of Medicine Northwestern University, Chicago, Illinois, USA
| | - Jie Xu
- Department of Biochemistry and Molecular Genetics, Feinberg School of Medicine Northwestern University, Chicago, Illinois, USA
| | - Baozhen Zhang
- Department of Biochemistry and Molecular Genetics, Feinberg School of Medicine Northwestern University, Chicago, Illinois, USA
| | - Ye Hou
- Department of Biochemistry and Molecular Genetics, Feinberg School of Medicine Northwestern University, Chicago, Illinois, USA
| | - Fan Song
- Department of Biochemistry and Molecular Genetics, Feinberg School of Medicine Northwestern University, Chicago, Illinois, USA
| | - Huijue Lyu
- Department of Biochemistry and Molecular Genetics, Feinberg School of Medicine Northwestern University, Chicago, Illinois, USA
| | - Feng Yue
- Department of Biochemistry and Molecular Genetics, Feinberg School of Medicine Northwestern University, Chicago, Illinois, USA.,Robert H. Lurie Comprehensive Cancer Center of Northwestern University, Chicago, Illinois, USA.,Correspondence:
| |
Collapse
|
33
|
Wang M, Li J, Wang P, Liu F, Liu Z, Zhao G, Xu Z, Pei L, Grover CE, Wendel JF, Wang K, Zhang X. Comparative genome analyses highlight transposon-mediated genome expansion and the evolutionary architecture of 3D genomic folding in cotton. Mol Biol Evol 2021; 38:3621-3636. [PMID: 33973633 PMCID: PMC8382922 DOI: 10.1093/molbev/msab128] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2020] [Revised: 04/19/2021] [Accepted: 04/28/2021] [Indexed: 12/12/2022] Open
Abstract
Transposable element (TE) amplification has been recognized as a driving force mediating genome size expansion and evolution, but the consequences for shaping 3D genomic architecture remains largely unknown in plants. Here, we report reference-grade genome assemblies for three species of cotton ranging 3-fold in genome size, namely Gossypium rotundifolium (K2), G. arboreum (A2), and G. raimondii (D5), using Oxford Nanopore Technologies. Comparative genome analyses document the details of lineage-specific TE amplification contributing to the large genome size differences (K2, 2.44 Gb; A2, 1.62 Gb; D5, 750.19 Mb) and indicate relatively conserved gene content and synteny relationships among genomes. We found that approximately 17% of syntenic genes exhibit chromatin status change between active (“A”) and inactive (“B”) compartments, and TE amplification was associated with the increase of the proportion of A compartment in gene regions (∼7,000 genes) in K2 and A2 relative to D5. Only 42% of topologically associating domain (TAD) boundaries were conserved among the three genomes. Our data implicate recent amplification of TEs following the formation of lineage-specific TAD boundaries. This study sheds light on the role of transposon-mediated genome expansion in the evolution of higher-order chromatin structure in plants.
Collapse
Affiliation(s)
- Maojun Wang
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan 430070, Hubei, China
| | - Jianying Li
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan 430070, Hubei, China
| | - Pengcheng Wang
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan 430070, Hubei, China
| | - Fang Liu
- State Key Laboratory of Cotton Biology, Institute of Cotton Research, Chinese Academy of Agricultural Sciences, Anyang 455000, Henan, China
| | - Zhenping Liu
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan 430070, Hubei, China
| | - Guannan Zhao
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan 430070, Hubei, China
| | - Zhongping Xu
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan 430070, Hubei, China
| | - Liuling Pei
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan 430070, Hubei, China
| | - Corrinne E Grover
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University, Ames, IA 50011, USA
| | - Jonathan F Wendel
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University, Ames, IA 50011, USA
| | - Kunbo Wang
- State Key Laboratory of Cotton Biology, Institute of Cotton Research, Chinese Academy of Agricultural Sciences, Anyang 455000, Henan, China
| | - Xianlong Zhang
- National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan 430070, Hubei, China
| |
Collapse
|
34
|
Wang Y, Liu Y, Xu Q, Xu Y, Cao K, Deng N, Wang R, Zhang X, Zheng R, Li G, Fang Y. TAD boundary and strength prediction by integrating sequence and epigenetic profile information. Brief Bioinform 2021; 22:6236056. [PMID: 33866359 DOI: 10.1093/bib/bbab139] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/25/2020] [Revised: 03/05/2021] [Accepted: 03/23/2021] [Indexed: 11/14/2022] Open
Abstract
Topologically associated domains (TADs) are one of the important higher order chromatin structures with various sizes in the eukaryotic genomes. TAD boundaries, as the flanking regions between adjacent domains, can restrict the interactions of regulatory elements, including enhancers and promoters, and are generally dynamic and variable in different cells. However, the influence of sequence and epigenetic profile-based features in the identification of TAD boundaries is largely unknown. In this work, we proposed a method called pTADS (prediction of TAD boundary and strength), to predict TAD boundaries and boundary strength across multiple cell lines with DNA sequence and epigenetic profile information. The performance was assessed in seven cell lines and three TAD calling methods. The results demonstrate that the TAD boundary can be well predicted by the selected shared features across multiple cell lines. Especially, the model can be transferable to predict the TAD boundary from one cell line to other cell lines. The boundary strength can be characterized by boundary score with good performance. The predicted TAD boundary and TAD boundary strength are further confirmed by three Hi-C contact matrix-based methods across multiple cell lines. The codes and datasets are available at https://github.com/chrom3DEpi/pTADS.
Collapse
Affiliation(s)
- Yunlong Wang
- College of Informatics, Huazhong Agricultural University, China
| | - Yaqi Liu
- College of Informatics, Huazhong Agricultural University, China
| | - Qian Xu
- College of Informatics, Huazhong Agricultural University, China
| | - Yao Xu
- College of Informatics, Huazhong Agricultural University, China
| | - Kai Cao
- College of Informatics, Huazhong Agricultural University, China
| | - Nan Deng
- College of Informatics, Huazhong Agricultural University, China
| | - Ruimin Wang
- College of Informatics, Huazhong Agricultural University, China
| | - Xueying Zhang
- College of Informatics, Huazhong Agricultural University, China
| | - Ruiqin Zheng
- College of Informatics, Huazhong Agricultural University, China
| | - Guoliang Li
- College of Informatics, Huazhong Agricultural University, Wuhan, 430070, P. R. China
| | - Yaping Fang
- College of Informatics, Huazhong Agricultural University, Wuhan, 430070, P. R. China
| |
Collapse
|
35
|
Gong H, Yang Y, Zhang S, Li M, Zhang X. Application of Hi-C and other omics data analysis in human cancer and cell differentiation research. Comput Struct Biotechnol J 2021; 19:2070-2083. [PMID: 33995903 PMCID: PMC8086027 DOI: 10.1016/j.csbj.2021.04.016] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2020] [Revised: 04/04/2021] [Accepted: 04/04/2021] [Indexed: 02/07/2023] Open
Abstract
With the development of 3C (chromosome conformation capture) and its derivative technology Hi-C (High-throughput chromosome conformation capture) research, the study of the spatial structure of the genomic sequence in the nucleus helps researchers understand the functions of biological processes such as gene transcription, replication, repair, and regulation. In this paper, we first introduce the research background and purpose of Hi-C data visualization analysis. After that, we discuss the Hi-C data analysis methods from genome 3D structure, A/B compartment, TADs (topologically associated domain), and loop detection. We also discuss how to apply genome visualization technologies to the identification of chromosome feature structures. We continue with a review of correlation analysis differences among multi-omics data, and how to apply Hi-C and other omics data analysis into cancer and cell differentiation research. Finally, we summarize the various problems in joint analyses based on Hi-C and other multi-omics data. We believe this review can help researchers better understand the progress and applications of 3D genome technology.
Collapse
Affiliation(s)
- Haiyan Gong
- Department of Computer Science and Technology, University of Science and Technology Beijing, Beijing 100083, China
- Beijing Advanced Innovation Center for Materials Genome Engineering, University of Science and Technology Beijing, Beijing 100083, China
- Beijing Key Laboratory of Knowledge Engineering for Materials Science, Beijing 100083, China
- Shunde Graduate School of University of Science and Technology Beijing, Foshan 528000, China
| | - Yi Yang
- Department of Computer Science and Technology, University of Science and Technology Beijing, Beijing 100083, China
| | - Sichen Zhang
- Department of Computer Science and Technology, University of Science and Technology Beijing, Beijing 100083, China
| | - Minghong Li
- Department of Computer Science and Technology, University of Science and Technology Beijing, Beijing 100083, China
| | - Xiaotong Zhang
- Department of Computer Science and Technology, University of Science and Technology Beijing, Beijing 100083, China
- Beijing Advanced Innovation Center for Materials Genome Engineering, University of Science and Technology Beijing, Beijing 100083, China
- Beijing Key Laboratory of Knowledge Engineering for Materials Science, Beijing 100083, China
- Shunde Graduate School of University of Science and Technology Beijing, Foshan 528000, China
| |
Collapse
|
36
|
Deschamps S, Crow JA, Chaidir N, Peterson-Burch B, Kumar S, Lin H, Zastrow-Hayes G, May GD. Chromatin loop anchors contain core structural components of the gene expression machinery in maize. BMC Genomics 2021; 22:23. [PMID: 33407087 PMCID: PMC7789236 DOI: 10.1186/s12864-020-07324-0] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2020] [Accepted: 12/14/2020] [Indexed: 12/28/2022] Open
Abstract
BACKGROUND Three-dimensional chromatin loop structures connect regulatory elements to their target genes in regions known as anchors. In complex plant genomes, such as maize, it has been proposed that loops span heterochromatic regions marked by higher repeat content, but little is known on their spatial organization and genome-wide occurrence in relation to transcriptional activity. RESULTS Here, ultra-deep Hi-C sequencing of maize B73 leaf tissue was combined with gene expression and open chromatin sequencing for chromatin loop discovery and correlation with hierarchical topologically-associating domains (TADs) and transcriptional activity. A majority of all anchors are shared between multiple loops from previous public maize high-resolution interactome datasets, suggesting a highly dynamic environment, with a conserved set of anchors involved in multiple interaction networks. Chromatin loop interiors are marked by higher repeat contents than the anchors flanking them. A small fraction of high-resolution interaction anchors, fully embedded in larger chromatin loops, co-locate with active genes and putative protein-binding sites. Combinatorial analyses indicate that all anchors studied here co-locate with at least 81.5% of expressed genes and 74% of open chromatin regions. Approximately 38% of all Hi-C chromatin loops are fully embedded within hierarchical TAD-like domains, while the remaining ones share anchors with domain boundaries or with distinct domains. Those various loop types exhibit specific patterns of overlap for open chromatin regions and expressed genes, but no apparent pattern of gene expression. In addition, up to 63% of all unique variants derived from a prior public maize eQTL dataset overlap with Hi-C loop anchors. Anchor annotation suggests that < 7% of all loops detected here are potentially devoid of any genes or regulatory elements. The overall organization of chromatin loop anchors in the maize genome suggest a loop modeling system hypothesized to resemble phase separation of repeat-rich regions. CONCLUSIONS Sets of conserved chromatin loop anchors mapping to hierarchical domains contains core structural components of the gene expression machinery in maize. The data presented here will be a useful reference to further investigate their function in regard to the formation of transcriptional complexes and the regulation of transcriptional activity in the maize genome.
Collapse
Affiliation(s)
| | - John A. Crow
- Corteva Agriscience, 8325 NW, 62nd Avenue, Johnston, Iowa, 50131 USA
| | - Nadia Chaidir
- Corteva Agriscience, 8325 NW, 62nd Avenue, Johnston, Iowa, 50131 USA
| | | | - Sunil Kumar
- Corteva Agriscience, The V-Acendas, Atria Block, 12th Floor, Plot No.17, Madhapur, Hyderabad, Telangana 500081 India
| | - Haining Lin
- Corteva Agriscience, 8325 NW, 62nd Avenue, Johnston, Iowa, 50131 USA
| | | | - Gregory D. May
- Corteva Agriscience, 8325 NW, 62nd Avenue, Johnston, Iowa, 50131 USA
| |
Collapse
|
37
|
Order and stochasticity in the folding of individual Drosophila genomes. Nat Commun 2021; 12:41. [PMID: 33397980 PMCID: PMC7782554 DOI: 10.1038/s41467-020-20292-z] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2020] [Accepted: 11/23/2020] [Indexed: 02/06/2023] Open
Abstract
Mammalian and Drosophila genomes are partitioned into topologically associating domains (TADs). Although this partitioning has been reported to be functionally relevant, it is unclear whether TADs represent true physical units located at the same genomic positions in each cell nucleus or emerge as an average of numerous alternative chromatin folding patterns in a cell population. Here, we use a single-nucleus Hi-C technique to construct high-resolution Hi-C maps in individual Drosophila genomes. These maps demonstrate chromatin compartmentalization at the megabase scale and partitioning of the genome into non-hierarchical TADs at the scale of 100 kb, which closely resembles the TAD profile in the bulk in situ Hi-C data. Over 40% of TAD boundaries are conserved between individual nuclei and possess a high level of active epigenetic marks. Polymer simulations demonstrate that chromatin folding is best described by the random walk model within TADs and is most suitably approximated by a crumpled globule build of Gaussian blobs at longer distances. We observe prominent cell-to-cell variability in the long-range contacts between either active genome loci or between Polycomb-bound regions, suggesting an important contribution of stochastic processes to the formation of the Drosophila 3D genome. Genomes are partitioned into topologically associating domains (TADs). Here the authors present single-nucleus Hi-C maps in Drosophila at 10 kb resolution, demonstrating the presence of chromatin compartments in individual nuclei, and partitioning of the genome into non-hierarchical TADs at the scale of 100 kb, which resembles population TAD profiles.
Collapse
|
38
|
A Brief Review of Current 3D Genomics Research. SYSTEMS MEDICINE 2021. [DOI: 10.1016/b978-0-12-801238-3.11594-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022] Open
|
39
|
Luo H, Li X, Fu H, Peng C. HiCHap: a package to correct and analyze the diploid Hi-C data. BMC Genomics 2020; 21:746. [PMID: 33109075 PMCID: PMC7590616 DOI: 10.1186/s12864-020-07165-x] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2020] [Accepted: 10/20/2020] [Indexed: 12/18/2022] Open
Abstract
Background In diploid cells, it is important to construct maternal and paternal Hi-C contact maps respectively since the two homologous chromosomes can differ in chromatin three-dimensional (3D) organization. Though previous softwares could construct diploid (maternal and paternal) Hi-C contact maps by using phased genetic variants, they all neglected the systematic biases in diploid Hi-C contact maps caused by variable genetic variant density in the genome. In addition, few of softwares provided quantitative analyses on allele-specific chromatin 3D organization, including compartment, topological domain and chromatin loop. Results In this work, we revealed the feature of allele-assignment bias caused by the variable genetic variant density, and then proposed a novel strategy to correct the systematic biases in diploid Hi-C contact maps. Based on the bias correction, we developed an integrated tool, called HiCHap, to perform read mapping, contact map construction, whole-genome identification of compartments, topological domains and chromatin loops, and allele-specific testing for diploid Hi-C data. Our results show that the correction on allele-assignment bias in HiCHap does significantly improve the quality of diploid Hi-C contact maps, which subsequently facilitates the whole-genome identification of diploid chromatin 3D organization, including compartments, topological domains and chromatin loops. Finally, HiCHap also supports the data analysis for haploid Hi-C maps without distinguishing two homologous chromosomes. Conclusions We provided an integrated package HiCHap to perform the data processing, bias correction and structural analysis for diploid Hi-C data. The source code and tutorial of software HiCHap are freely available at https://pypi.org/project/HiCHap/. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-020-07165-x.
Collapse
Affiliation(s)
- Han Luo
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, 430070, China
| | - Xinxin Li
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, 430070, China
| | - Haitao Fu
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, 430070, China
| | - Cheng Peng
- Center for Life Sciences, School of Life Sciences, Yunnan University, Kunming, 650500, China.
| |
Collapse
|
40
|
Highly interconnected enhancer communities control lineage-determining genes in human mesenchymal stem cells. Nat Genet 2020; 52:1227-1238. [PMID: 33020665 DOI: 10.1038/s41588-020-0709-z] [Citation(s) in RCA: 45] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2019] [Accepted: 09/01/2020] [Indexed: 12/20/2022]
Abstract
Adipocyte differentiation is driven by waves of transcriptional regulators that reprogram the enhancer landscape and change the wiring of the promoter interactome. Here, we use high-throughput chromosome conformation enhancer capture to interrogate the role of enhancer-to-enhancer interactions during differentiation of human mesenchymal stem cells. We find that enhancers form an elaborate network that is dynamic during differentiation and coupled with changes in enhancer activity. Transcription factors (TFs) at baited enhancers amplify TF binding at target enhancers, a phenomenon we term cross-interaction stabilization of TFs. Moreover, highly interconnected enhancers (HICE) act as integration hubs orchestrating differentiation by the formation of three-dimensional enhancer communities, inside which, HICE, and other enhancers, converge on phenotypically important gene promoters. Collectively, these results indicate that enhancer interactions play a key role in the regulation of enhancer function, and that HICE are important for both signal integration and compartmentalization of the genome.
Collapse
|
41
|
Lyu H, Li L, Wu Z, Wang T, Zheng J, Wang H. TADBD: a sensitive and fast method for detection of typologically associated domain boundaries. Biotechniques 2020; 69:376-383. [DOI: 10.2144/btn-2019-0165] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022] Open
Abstract
A topologically associated domain (TAD) is a self-interacting genomic block. Detection of TAD boundaries on Hi-C contact matrix is one of the most important issues in the analysis of 3D genome architecture at TAD level. Here, we present TAD boundary detection (TADBD), a sensitive and fast computational method for detection of TAD boundaries on Hi-C contact matrix. This method implements a Haar-based algorithm by considering Haar diagonal template, acceleration via a compact integrogram, multi-scale aggregation at template size and statistical filtering. In most cases, comparison results from simulated and experimental data show that TADBD outperforms the other five methods. In addition, a new R package for TADBD is freely available online.
Collapse
Affiliation(s)
- Hongqiang Lyu
- School of Electronic & Information Engineering, Xi’an Jiaotong University, Xi’an 710049, China
| | - Lin Li
- School of Electronic & Information Engineering, Xi’an Jiaotong University, Xi’an 710049, China
| | - Zhifang Wu
- School of Electronic & Information Engineering, Xi’an Jiaotong University, Xi’an 710049, China
| | - Tian Wang
- School of Electronic & Information Engineering, Xi’an Jiaotong University, Xi’an 710049, China
- School of Automation Science & Electrical Engineering, Beihang University, Beijing 100191, China
| | - Jiguang Zheng
- School of Electronic & Information Engineering, Xi’an Jiaotong University, Xi’an 710049, China
| | - Hongda Wang
- School of Electronic & Information Engineering, Xi’an Jiaotong University, Xi’an 710049, China
| |
Collapse
|
42
|
Advances in technologies for 3D genomics research. SCIENCE CHINA-LIFE SCIENCES 2020; 63:811-824. [PMID: 32394244 DOI: 10.1007/s11427-019-1704-2] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/03/2020] [Accepted: 04/21/2020] [Indexed: 01/08/2023]
Abstract
The spatial structure of the orderly organized chromatin in the nucleus has important roles in maintaining normal cell function and in regulation of gene expression, and the high-throughput Hi-C and ChIA-PET methods have been widely used in various biological studies for determining potential spatial genome structures and their functions. However, there are still great difficulties and challenges in three-dimensional (3D) genomics research. More efficient, economical, and unbiased approaches to studying 3D genomics need to be developed for more widespread and easier applications. Here, we review the most recent studies on new 3D genomics research technologies, such as improvements of the traditional Hi-C and ChIA-PET methods, new approaches based on non-proximal-ligation strategies, and imaging-based methods improved in recent years. Especially, we review the CRISPR-based methods for functional validations in 3D genomics, which could be the forthcoming directions. We hope this review can show some insights into the potential improvements for future 3D genomics.
Collapse
|
43
|
Wang G, Meng Q, Xia B, Zhang S, Lv J, Zhao D, Li Y, Wang X, Zhang L, Cooke JP, Cao Q, Chen K. TADsplimer reveals splits and mergers of topologically associating domains for epigenetic regulation of transcription. Genome Biol 2020; 21:84. [PMID: 32241291 PMCID: PMC7114812 DOI: 10.1186/s13059-020-01992-7] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2019] [Accepted: 03/09/2020] [Indexed: 12/19/2022] Open
Abstract
We present TADsplimer, the first computational tool to systematically detect topologically associating domain (TAD) splits and mergers across the genome between Hi-C samples. TADsplimer recaptures splits and mergers of TADs with high accuracy in simulation analyses and defines hundreds of TAD splits and mergers between pairs of different cell types, such as endothelial cells and fibroblasts. Our work reveals a key role for TAD remodeling in epigenetic regulation of transcription and delivers the first tool for the community to perform dynamic analysis of TAD splits and mergers in numerous biological and disease models.
Collapse
Affiliation(s)
- Guangyu Wang
- Center for Bioinformatics and Computational Biology, Houston Methodist Research Institute, Houston, TX, 77030, USA
- Center for Cardiovascular Regeneration, Houston Methodist Research Institute, Houston, TX, 77030, USA
- Department of Cardiothoracic Surgery, Weill Cornell Medicine, Cornell University, New York, NY, 10065, USA
| | - Qingshu Meng
- Department of Urology, Northwestern University Feinberg School of Medicine, Chicago, IL, 60611, USA
- Robert H. Lurie Comprehensive Cancer Center, Northwestern University Feinberg School of Medicine, Chicago, IL, 60611, USA
| | - Bo Xia
- Center for Bioinformatics and Computational Biology, Houston Methodist Research Institute, Houston, TX, 77030, USA
- Center for Cardiovascular Regeneration, Houston Methodist Research Institute, Houston, TX, 77030, USA
- Department of Cardiothoracic Surgery, Weill Cornell Medicine, Cornell University, New York, NY, 10065, USA
| | - Shuo Zhang
- Center for Bioinformatics and Computational Biology, Houston Methodist Research Institute, Houston, TX, 77030, USA
- Center for Cardiovascular Regeneration, Houston Methodist Research Institute, Houston, TX, 77030, USA
- Department of Cardiothoracic Surgery, Weill Cornell Medicine, Cornell University, New York, NY, 10065, USA
| | - Jie Lv
- Center for Bioinformatics and Computational Biology, Houston Methodist Research Institute, Houston, TX, 77030, USA
- Center for Cardiovascular Regeneration, Houston Methodist Research Institute, Houston, TX, 77030, USA
- Department of Cardiothoracic Surgery, Weill Cornell Medicine, Cornell University, New York, NY, 10065, USA
| | - Dongyu Zhao
- Center for Bioinformatics and Computational Biology, Houston Methodist Research Institute, Houston, TX, 77030, USA
- Center for Cardiovascular Regeneration, Houston Methodist Research Institute, Houston, TX, 77030, USA
- Department of Cardiothoracic Surgery, Weill Cornell Medicine, Cornell University, New York, NY, 10065, USA
| | - Yanqiang Li
- Center for Bioinformatics and Computational Biology, Houston Methodist Research Institute, Houston, TX, 77030, USA
- Center for Cardiovascular Regeneration, Houston Methodist Research Institute, Houston, TX, 77030, USA
- Department of Cardiothoracic Surgery, Weill Cornell Medicine, Cornell University, New York, NY, 10065, USA
| | - Xin Wang
- Center for Bioinformatics and Computational Biology, Houston Methodist Research Institute, Houston, TX, 77030, USA
- Center for Cardiovascular Regeneration, Houston Methodist Research Institute, Houston, TX, 77030, USA
- Department of Cardiothoracic Surgery, Weill Cornell Medicine, Cornell University, New York, NY, 10065, USA
| | - Lili Zhang
- Center for Cardiovascular Regeneration, Houston Methodist Research Institute, Houston, TX, 77030, USA
- Department of Cardiothoracic Surgery, Weill Cornell Medicine, Cornell University, New York, NY, 10065, USA
| | - John P Cooke
- Center for Cardiovascular Regeneration, Houston Methodist Research Institute, Houston, TX, 77030, USA
- Department of Cardiothoracic Surgery, Weill Cornell Medicine, Cornell University, New York, NY, 10065, USA
| | - Qi Cao
- Department of Urology, Northwestern University Feinberg School of Medicine, Chicago, IL, 60611, USA.
- Robert H. Lurie Comprehensive Cancer Center, Northwestern University Feinberg School of Medicine, Chicago, IL, 60611, USA.
| | - Kaifu Chen
- Center for Bioinformatics and Computational Biology, Houston Methodist Research Institute, Houston, TX, 77030, USA.
- Center for Cardiovascular Regeneration, Houston Methodist Research Institute, Houston, TX, 77030, USA.
- Department of Cardiothoracic Surgery, Weill Cornell Medicine, Cornell University, New York, NY, 10065, USA.
| |
Collapse
|
44
|
Dong P, Tu X, Li H, Zhang J, Grierson D, Li P, Zhong S. Tissue-specific Hi-C analyses of rice, foxtail millet and maize suggest non-canonical function of plant chromatin domains. JOURNAL OF INTEGRATIVE PLANT BIOLOGY 2020; 62:201-217. [PMID: 30920762 DOI: 10.1111/jipb.12809] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/19/2019] [Accepted: 03/21/2019] [Indexed: 06/09/2023]
Abstract
Chromatins are not randomly packaged in the nucleus and their organization plays important roles in transcription regulation, which is best studied in the mammalian models. Using in situ Hi-C, we have compared the 3D chromatin architectures of rice mesophyll and endosperm, foxtail millet bundle sheath and mesophyll, and maize bundle sheath, mesophyll and endosperm tissues. We found that their global A/B compartment partitions are stable across tissues, while local A/B compartment has tissue-specific dynamic associated with differential gene expression. Plant domains are largely stable across tissues, while new domain border formations are often associated with transcriptional activation in the region. Genes inside plant domains are not conserved across species, and lack significant co-expression behavior unlike those in mammalian TADs. Although we only observed chromatin loops between gene islands in the large genomes, the maize loop gene pairs' syntenic orthologs have shorter physical distances in small genome monocots, suggesting that loops instead of domains might have conserved biological function. Our study showed that plants' chromatin features might not have conserved biological functions as the mammalian ones.
Collapse
Affiliation(s)
- Pengfei Dong
- State Key Laboratory of Agrobiotechnology, School of Life Sciences, The Chinese University of Hong Kong, Hong Kong, China
| | - Xiaoyu Tu
- State Key Laboratory of Agrobiotechnology, School of Life Sciences, The Chinese University of Hong Kong, Hong Kong, China
| | - Haoxuan Li
- State Key Laboratory of Agrobiotechnology, School of Life Sciences, The Chinese University of Hong Kong, Hong Kong, China
| | - Jianhua Zhang
- State Key Laboratory of Agrobiotechnology, School of Life Sciences, The Chinese University of Hong Kong, Hong Kong, China
| | - Donald Grierson
- Zhejiang Provincial Key Laboratory of Horticultural Plant Integrative Biology, Zhejiang University, Hangzhou, Zhejiang, China
| | - Pinghua Li
- State Key Laboratory of Crop Biology, College of Agronomy, Shandong Agricultural University, Tai'an, Shandong, China
| | - Silin Zhong
- State Key Laboratory of Agrobiotechnology, School of Life Sciences, The Chinese University of Hong Kong, Hong Kong, China
| |
Collapse
|
45
|
Razin SV, Ulianov SV, Gavrilov AA. 3D Genomics. Mol Biol 2019; 53:802-812. [DOI: 10.1134/s0026893319060153] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2019] [Revised: 06/01/2019] [Accepted: 06/03/2019] [Indexed: 08/30/2023]
|
46
|
Muñiz Moreno MDM, Brault V, Birling MC, Pavlovic G, Herault Y. Modeling Down syndrome in animals from the early stage to the 4.0 models and next. PROGRESS IN BRAIN RESEARCH 2019; 251:91-143. [PMID: 32057313 DOI: 10.1016/bs.pbr.2019.08.001] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
The genotype-phenotype relationship and the physiopathology of Down Syndrome (DS) have been explored in the last 20 years with more and more relevant mouse models. From the early age of transgenesis to the new CRISPR/CAS9-derived chromosomal engineering and the transchromosomic technologies, mouse models have been key to identify homologous genes or entire regions homologous to the human chromosome 21 that are necessary or sufficient to induce DS features, to investigate the complexity of the genetic interactions that are involved in DS and to explore therapeutic strategies. In this review we report the new developments made, how genomic data and new genetic tools have deeply changed our way of making models, extended our panel of animal models, and increased our understanding of the neurobiology of the disease. But even if we have made an incredible progress which promises to make DS a curable condition, we are facing new research challenges to nurture our knowledge of DS pathophysiology as a neurodevelopmental disorder with many comorbidities during ageing.
Collapse
Affiliation(s)
- Maria Del Mar Muñiz Moreno
- Université de Strasbourg, CNRS, INSERM, Institut de Génétique et de Biologie Moléculaire et Cellulaire, Illkirch, France
| | - Véronique Brault
- Université de Strasbourg, CNRS, INSERM, Institut de Génétique et de Biologie Moléculaire et Cellulaire, Illkirch, France
| | - Marie-Christine Birling
- Université de Strasbourg, CNRS, INSERM, PHENOMIN Institut Clinique de la Souris, Illkirch, France
| | - Guillaume Pavlovic
- Université de Strasbourg, CNRS, INSERM, PHENOMIN Institut Clinique de la Souris, Illkirch, France
| | - Yann Herault
- Université de Strasbourg, CNRS, INSERM, Institut de Génétique et de Biologie Moléculaire et Cellulaire, Illkirch, France; Université de Strasbourg, CNRS, INSERM, PHENOMIN Institut Clinique de la Souris, Illkirch, France.
| |
Collapse
|
47
|
Malik L, Patro R. Rich Chromatin Structure Prediction from Hi-C Data. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2019; 16:1448-1458. [PMID: 29994683 DOI: 10.1109/tcbb.2018.2851200] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Recent studies involving the 3-dimensional conformation of chromatin have revealed the important role it has to play in different processes within the cell. These studies have also led to the discovery of densely interacting segments of the chromosome, called topologically associating domains. The accurate identification of these domains from Hi-C interaction data is an interesting and important computational problem for which numerous methods have been proposed. Unfortunately, most existing algorithms designed to identify these domains assume that they are non-overlapping whereas there is substantial evidence to believe a nested structure exists. We present a methodology to predict hierarchical chromatin domains using chromatin conformation capture data. Our method predicts domains at different resolutions, calculated using intrinsic properties of the chromatin data, and effectively clusters these to construct the hierarchy. At each individual level, the domains are non-overlapping in such a way that the intra-domain interaction frequencies are maximized. We show that our predicted structure is highly enriched for actively transcribing housekeeping genes and various chromatin markers, including CTCF, around the domain boundaries. We also show that large-scale domains, at multiple resolutions within our hierarchy, are conserved across cell types and species. We also provide comparisons against existing tools for extracting hierarchical domains. Our software, Matryoshka, is written in C++11 and licensed under GPL v3; it is available at https://github.com/COMBINE-lab/matryoshka.
Collapse
|
48
|
Xie T, Zhang FG, Zhang HY, Wang XT, Hu JH, Wu XM. Biased gene retention during diploidization in Brassica linked to three-dimensional genome organization. NATURE PLANTS 2019; 5:822-832. [PMID: 31383969 DOI: 10.1038/s41477-019-0479-8] [Citation(s) in RCA: 36] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/06/2019] [Accepted: 06/19/2019] [Indexed: 05/07/2023]
Abstract
The non-random three-dimensional (3D) organization of the genome in the nucleus is critical to gene regulation and genome function. Using high-throughput chromatin conformation capture, we generated chromatin interaction maps for Brassica rapa and Brassica oleracea at a high resolution and characterized the conservation and divergence of chromatin organization in these two species. Large-scale chromatin structures, including A/B compartments and topologically associating domains, are notably conserved between B. rapa and B. oleracea, yet their KNOT structures are highly divergent. We found that genes retained in less fractionated subgenomes exhibited stronger interaction strengths, and diploidization-resistant duplicates retained in pairs or triplets are more likely to be colocalized in both B. rapa and B. oleracea. These observations suggest that spatial constraint in duplicated genes is correlated to their biased retention in the diploidization process. In addition, we found strong similarities in the epigenetic modification and Gene Ontology terms of colocalized paralogues, which were largely conserved across B. rapa and B. oleracea, indicating functional constraints on their 3D positioning in the nucleus. This study presents an investigation of the spatial organization of genomes in Brassica and provides insights on the role of 3D organization in the genome evolution of this genus.
Collapse
Affiliation(s)
- Ting Xie
- Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agriculture and Rural Affairs, Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences, Wuhan, China.
| | - Fu-Gui Zhang
- Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agriculture and Rural Affairs, Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences, Wuhan, China
| | - Hong-Yu Zhang
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, China
| | - Xiao-Tao Wang
- Department of Biochemistry and Molecular Biology, College of Medicine, Pennsylvania State University, Hershey, PA, USA
| | - Ji-Hong Hu
- Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agriculture and Rural Affairs, Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences, Wuhan, China
| | - Xiao-Ming Wu
- Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agriculture and Rural Affairs, Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences, Wuhan, China
| |
Collapse
|
49
|
Chen F, Li G, Zhang MQ, Chen Y. HiCDB: a sensitive and robust method for detecting contact domain boundaries. Nucleic Acids Res 2019; 46:11239-11250. [PMID: 30184171 PMCID: PMC6265446 DOI: 10.1093/nar/gky789] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2018] [Accepted: 08/27/2018] [Indexed: 12/20/2022] Open
Abstract
Contact domains are closely linked to gene regulation and lineage commitment, while current understanding of contact domains and their boundaries is still limited. Here, we present a novel method HiCDB, which is constructively based on local relative insulation metric and multi-scale aggregation approach to detect contact domain boundaries (CDBs) on Hi-C maps. Compared with other ‘state-of-art’ methods, HiCDB shows improved sensitivity and specificity in determining CDBs at various Hi-C resolutions. The superiority of HiCDB enabled us to study the epigenetic features of detected CDBs and showed enrichment of architectural proteins and cell-type-specific transcription factor binding sites at CDBs. The further comparison of GM12878 and IMR90 Hi-C datasets suggested that cell-type-specific CDBs are marked by active regulatory signals and correlate with activation of nearby cell identity genes.
Collapse
Affiliation(s)
- Fengling Chen
- MOE Key Laboratory of Bioinformatics, Center for Synthetic and Systems Biology, Bioinformatics Division, BNRist, Department of Automation, Tsinghua University, Beijing 100084, China
| | - Guipeng Li
- MOE Key Laboratory of Bioinformatics, Center for Synthetic and Systems Biology, Bioinformatics Division, BNRist, Department of Automation, Tsinghua University, Beijing 100084, China.,Department of Biology, Medi-X Institute, SUSTech Academy for Advanced Interdisciplinary Studies, Southern University of Science and Technology, Shenzhen 518055, China
| | - Michael Q Zhang
- MOE Key Laboratory of Bioinformatics, Center for Synthetic and Systems Biology, Bioinformatics Division, BNRist, Department of Automation, Tsinghua University, Beijing 100084, China.,MOE Key Laboratory of Bioinformatics, Center for Synthetic and Systems Biology, School of Medicine, Tsinghua University, Beijing 100084, China.,Department of Biological Sciences, Center for Systems Biology, The University of Texas, Dallas 800 West Campbell Road, RL11, Richardson, TX 75080-3021, USA
| | - Yang Chen
- MOE Key Laboratory of Bioinformatics, Center for Synthetic and Systems Biology, Bioinformatics Division, BNRist, Department of Automation, Tsinghua University, Beijing 100084, China.,MOE Key Laboratory of Bioinformatics, Center for Synthetic and Systems Biology, School of Medicine, Tsinghua University, Beijing 100084, China
| |
Collapse
|
50
|
Abstract
In the epigenetics field, large-scale functional genomics datasets of ever-increasing size and complexity have been produced using experimental techniques based on high-throughput sequencing. In particular, the study of the 3D organization of chromatin has raised increasing interest, thanks to the development of advanced experimental techniques. In this context, Hi-C has been widely adopted as a high-throughput method to measure pairwise contacts between virtually any pair of genomic loci, thus yielding unprecedented challenges for analyzing and handling the resulting complex datasets. In this review, we focus on the increasing complexity of available Hi-C datasets, which parallels the adoption of novel protocol variants. We also review the complexity of the multiple data analysis steps required to preprocess Hi-C sequencing reads and extract biologically meaningful information. Finally, we discuss solutions for handling and visualizing such large genomics datasets.
Collapse
|