CWL-Based Analysis Pipeline for Hi-C Data: From FASTQ Files to Matrices

Overview

Journal Methods Mol Biol

Specialty Molecular Biology

Date 2024 Sep 16

PMID 39283448

Authors

Hisashi Miura

Rory T Cerbus

Izumi Noda

Ichiro Hiratani

Affiliations

Soon will be listed here.

Abstract

Over a decade has passed since the development of the Hi-C method for genome-wide analysis of 3D genome organization. Hi-C utilizes next-generation sequencing (NGS) technology to generate large-scale chromatin interaction data, which has accumulated across a diverse range of species and cell types, particularly in eukaryotes. There is thus a growing need to streamline the process of Hi-C data analysis to utilize these data sets effectively. Hi-C generates data that are much larger compared to other NGS techniques such as chromatin immunoprecipitation sequencing (ChIP-seq) or RNA-seq, making the data reanalysis process computationally expensive. In an effort to bridge this resource gap, the 4D Nucleome (4DN) Data Portal has reanalyzed approximately 600 Hi-C data sets, allowing users to access and utilize the analyzed data. In this chapter, we provide detailed instructions for the implementation of the common workflow language (CWL)-based Hi-C analysis pipeline adopted by the 4DN Data Portal ecosystem. This reproducible and portable pipeline generates standard Hi-C contact matrices in formats such as .hic or .mcool from FASTQ files. It enables users to output their own Hi-C data in the same format as those registered in the 4DN Data portal, facilitating comparative analysis using data registered in the portal. Our custom-made scripts are available on GitHub at https://github.com/kuzobuta/4dn_cwl_pipeline .

References

Ghosh R, Meyer B . Spatial Organization of Chromatin: Emergence of Chromatin Structure During Development. Annu Rev Cell Dev Biol. 2021; 37:199-232. PMC: 8664233. DOI: 10.1146/annurev-cellbio-032321-035734. View

Marchal C, Sima J, Gilbert D . Control of DNA replication timing in the 3D genome. Nat Rev Mol Cell Biol. 2019; 20(12):721-737. PMC: 11567694. DOI: 10.1038/s41580-019-0162-y. View

Stanic M, Mekhail K . Integration of DNA damage responses with dynamic spatial genome organization. Trends Genet. 2021; 38(3):290-304. DOI: 10.1016/j.tig.2021.08.016. View

Lieberman-Aiden E, van Berkum N, Williams L, Imakaev M, Ragoczy T, Telling A . Comprehensive mapping of long-range interactions reveals folding principles of the human genome. Science. 2009; 326(5950):289-93. PMC: 2858594. DOI: 10.1126/science.1181369. View

Dixon J, Selvaraj S, Yue F, Kim A, Li Y, Shen Y . Topological domains in mammalian genomes identified by analysis of chromatin interactions. Nature. 2012; 485(7398):376-80. PMC: 3356448. DOI: 10.1038/nature11082. View

Nora E, Lajoie B, Schulz E, Giorgetti L, Okamoto I, Servant N . Spatial partitioning of the regulatory landscape of the X-inactivation centre. Nature. 2012; 485(7398):381-5. PMC: 3555144. DOI: 10.1038/nature11049. View

Rao S, Huntley M, Durand N, Stamenova E, Bochkov I, Robinson J . A 3D map of the human genome at kilobase resolution reveals principles of chromatin looping. Cell. 2014; 159(7):1665-80. PMC: 5635824. DOI: 10.1016/j.cell.2014.11.021. View

Jin F, Li Y, Dixon J, Selvaraj S, Ye Z, Lee A . A high-resolution map of the three-dimensional chromatin interactome in human cells. Nature. 2013; 503(7475):290-4. PMC: 3838900. DOI: 10.1038/nature12644. View

Nichols M, Corces V . Principles of 3D compartmentalization of the human genome. Cell Rep. 2021; 35(13):109330. PMC: 8265014. DOI: 10.1016/j.celrep.2021.109330. View

10.

Abdulla A, Salari H, Tortora M, Vaillant C, Jost D . 4D epigenomics: deciphering the coupling between genome folding and epigenomic regulation with biophysical modeling. Curr Opin Genet Dev. 2023; 79:102033. DOI: 10.1016/j.gde.2023.102033. View

11.

Ryba T, Hiratani I, Lu J, Itoh M, Kulik M, Zhang J . Evolutionarily conserved replication timing profiles predict long-range chromatin interactions and distinguish closely related cell types. Genome Res. 2010; 20(6):761-70. PMC: 2877573. DOI: 10.1101/gr.099655.109. View

12.

Arnould C, Rocher V, Saur F, Bader A, Muzzopappa F, Collins S . Chromatin compartmentalization regulates the response to DNA damage. Nature. 2023; 623(7985):183-192. PMC: 10620078. DOI: 10.1038/s41586-023-06635-y. View

13.

Spracklin G, Abdennur N, Imakaev M, Chowdhury N, Pradhan S, Mirny L . Diverse silent chromatin states modulate genome compartmentalization and loop extrusion barriers. Nat Struct Mol Biol. 2022; 30(1):38-51. PMC: 9851908. DOI: 10.1038/s41594-022-00892-7. View

14.

Berlivet S, Paquette D, Dumouchel A, Langlais D, Dostie J, Kmita M . Clustering of tissue-specific sub-TADs accompanies the regulation of HoxA genes in developing limbs. PLoS Genet. 2014; 9(12):e1004018. PMC: 3873244. DOI: 10.1371/journal.pgen.1004018. View

15.

Phillips-Cremins J, Sauria M, Sanyal A, Gerasimova T, Lajoie B, Bell J . Architectural protein subclasses shape 3D organization of genomes during lineage commitment. Cell. 2013; 153(6):1281-95. PMC: 3712340. DOI: 10.1016/j.cell.2013.04.053. View

16.

Vian L, Pekowska A, Rao S, Kieffer-Kwon K, Jung S, Baranello L . The Energetics and Physiological Impact of Cohesin Extrusion. Cell. 2018; 173(5):1165-1178.e20. PMC: 6065110. DOI: 10.1016/j.cell.2018.03.072. View

17.

Kraft K, Magg A, Heinrich V, Riemenschneider C, Schopflin R, Markowski J . Serial genomic inversions induce tissue-specific architectural stripes, gene misexpression and congenital malformations. Nat Cell Biol. 2019; 21(3):305-310. DOI: 10.1038/s41556-019-0273-x. View

18.

Guo Y, Al-Jibury E, Garcia-Millan R, Ntagiantas K, King J, Nash A . Chromatin jets define the properties of cohesin-driven in vivo loop extrusion. Mol Cell. 2022; 82(20):3769-3780.e5. DOI: 10.1016/j.molcel.2022.09.003. View

19.

Tashiro S, Lanctot C . The International Nucleome Consortium. Nucleus. 2015; 6(2):89-92. PMC: 4615172. DOI: 10.1080/19491034.2015.1022703. View

20.

Cremer T, Cremer M, Hubner B, Strickfaden H, Smeets D, Popken J . The 4D nucleome: Evidence for a dynamic nuclear landscape based on co-aligned active and inactive nuclear compartments. FEBS Lett. 2015; 589(20 Pt A):2931-43. DOI: 10.1016/j.febslet.2015.05.037. View