Capture Hi-C Research Articles

Algorithmic considerations when analysing capture Hi-C data.

Chromosome conformation capture methodologies have provided insight into the effect of 3D genomic architecture on gene regulation. Capture Hi-C (CHi-C) is a recent extension of Hi-C that improves the effective resolution of chromatin interactions by enriching for defined regions of biological relevance. The varying targeting efficiency between capture regions, however, introduces bias not present in conventional Hi-C, making analysis more complicated. Here we consider salient features of an algorithm that should be considered in evaluating the performance of a program used to analyse CHi-C data in order to infer meaningful interactions. We use the program CHICAGO to analyse promotor capture Hi-C data generated on 28 different cell lines as a case study.

3DIV update for 2021: a comprehensive resource of 3D genome and 3D cancer genome.

Three-dimensional (3D) genome organization is tightly coupled with gene regulation in various biological processes and diseases. In cancer, various types of large-scale genomic rearrangements can disrupt the 3D genome, leading to oncogenic gene expression. However, unraveling the pathogenicity of the 3D cancer genome remains a challenge since closer examinations have been greatly limited due to the lack of appropriate tools specialized for disorganized higher-order chromatin structure. Here, we updated a 3D-genome Interaction Viewer and database named 3DIV by uniformly processing ∼230 billion raw Hi-C reads to expand our contents to the 3D cancer genome. The updates of 3DIV are listed as follows: (i) the collection of 401 samples including 220 cancer cell line/tumor Hi-C data, 153 normal cell line/tissue Hi-C data, and 28 promoter capture Hi-C data, (ii) the live interactive manipulation of the 3D cancer genome to simulate the impact of structural variations and (iii) the reconstruction of Hi-C contact maps by user-defined chromosome order to investigate the 3D genome of the complex genomic rearrangement. In summary, the updated 3DIV will be the most comprehensive resource to explore the gene regulatory effects of both the normal and cancer 3D genome. ‘3DIV’ is freely available at http://3div.kr.

A computational framework for identifying the transcription factors involved in enhancer-promoter loop formation

The pairwise interaction between transcription factors (TFs) plays an important role in enhancer-promoter loop formation. Although thousands of TFs in the human genome have been found, only a few TF pairs have been demonstrated to be related to loop formation. It is still a challenge to determine which TF pairs could be involved in the enhancer-promoter regulation network. This work describes a computational framework to identify TF pairs in enhancer-promoter regulation. By integrating different levels of data derived from Promoter Capture Hi-C, chromatin immunoprecipitation sequencing (ChIP-seq) of histone marks, RNA-seq, protein-protein interaction (PPI), and TF motif, we identified 361 significant TF pairs and constructed a TF interaction network. From the network, we found several hub-TFs, which may have important roles in the regulation of long-range interactions. Our studies extended TF pairs identified in other experimental and computational approaches. These findings will help the further study of long-range interactions between enhancers and promoters.

Altered transcriptional and chromatin responses to rhinovirus in bronchial epithelial cells from adults with asthma

There is a life-long relationship between rhinovirus (RV) infection and the development and clinical manifestations of asthma. In this study we demonstrate that cultured primary bronchial epithelial cells from adults with asthma (n = 9) show different transcriptional and chromatin responses to RV infection compared to those without asthma (n = 9). Both the number and magnitude of transcriptional and chromatin responses to RV were muted in cells from asthma cases compared to controls. Pathway analysis of the transcriptionally responsive genes revealed enrichments of apoptotic pathways in controls but inflammatory pathways in asthma cases. Using promoter capture Hi-C we tethered regions of RV-responsive chromatin to RV-responsive genes and showed enrichment of these regions and genes at asthma GWAS loci. Taken together, our studies indicate a delayed or prolonged inflammatory state in cells from asthma cases and highlight genes that may contribute to genetic risk for asthma.

Algorithmic considerations when analysing capture Hi-C data

Algorithmic considerations when analysing capture Hi-C data

Algorithmic considerations when analysing capture Hi-C data

Algorithmic considerations when analysing capture Hi-C data

Chromatin interactions in differentiating keratinocytes reveal novel atopic dermatitis– and psoriasis-associated genes

Hundreds of variants associated with atopic dermatitis (AD) and psoriasis, 2 common inflammatory skin disorders, have previously been discovered through genome-wide association studies (GWASs). The majority of these variants are in noncoding regions, and their target genes remain largely unclear. We sought to understand the effects of these noncoding variants on the development of AD and psoriasis by linking them to the genes that they regulate. We constructed genomic 3-dimensional maps of human keratinocytes during differentiation by using targeted chromosome conformation capture (Capture Hi-C) targeting more than 20,000 promoters and 214 GWAS variants and combined these data with transcriptome and epigenomic data sets. We validated our results with reporter assays, clustered regularly interspaced short palindromic repeats activation, and examination of patient gene expression from previous studies. We identified 118 target genes of 82 AD and psoriasis GWAS variants. Differential expression of 58 of the 118 target genes (49%) occurred in either AD or psoriatic lesions, many of which were not previously linked to any skin disease. We highlighted the genes AFG1L, CLINT1, ADO, LINC00302, and RP1-140J1.1 and provided further evidence for their potential roles in AD and psoriasis. Our work focused on skin barrier pathology through investigation of the interaction profile of GWAS variants during keratinocyte differentiation. We have provided a catalogue of candidate genes that could modulate the risk of AD and psoriasis. Given that only 35% of the target genes are the gene nearest to the known GWAS variants, we expect that our work will contribute to the discovery of novel pathways involved in AD and psoriasis.

Single-Cell Transcriptome Analysis Reveals Dynamic Cell Populations and Differential Gene Expression Patterns in Control and Aneurysmal Human Aortic Tissue.

Ascending thoracic aortic aneurysm (ATAA) is caused by the progressive weakening and dilatation of the aortic wall and can lead to aortic dissection, rupture, and other life-threatening complications. To improve our understanding of ATAA pathogenesis, we aimed to comprehensively characterize the cellular composition of the ascending aortic wall and to identify molecular alterations in each cell population of human ATAA tissues. We performed single-cell RNA sequencing analysis of ascending aortic tissues from 11 study participants, including 8 patients with ATAA (4 women and 4 men) and 3 control subjects (2 women and 1 man). Cells extracted from aortic tissue were analyzed and categorized with single-cell RNA sequencing data to perform cluster identification. ATAA-related changes were then examined by comparing the proportions of each cell type and the gene expression profiles between ATAA and control tissues. We also examined which genes may be critical for ATAA by performing the integrative analysis of our single-cell RNA sequencing data with publicly available data from genome-wide association studies. We identified 11 major cell types in human ascending aortic tissue; the high-resolution reclustering of these cells further divided them into 40 subtypes. Multiple subtypes were observed for smooth muscle cells, macrophages, and T lymphocytes, suggesting that these cells have multiple functional populations in the aortic wall. In general, ATAA tissues had fewer nonimmune cells and more immune cells, especially T lymphocytes, than control tissues did. Differential gene expression data suggested the presence of extensive mitochondrial dysfunction in ATAA tissues. In addition, integrative analysis of our single-cell RNA sequencing data with public genome-wide association study data and promoter capture Hi-C data suggested that the erythroblast transformation-specific related gene(ERG) exerts an important role in maintaining normal aortic wall function. Our study provides a comprehensive evaluation of the cellular composition of the ascending aortic wall and reveals how the gene expression landscape is altered in human ATAA tissue. The information from this study makes important contributions to our understanding of ATAA formation and progression.

Epigenomics and transcriptomics of systemic sclerosis CD4+ T cells reveal long-range dysregulation of key inflammatory pathways mediated by disease-associated susceptibility loci

BackgroundSystemic sclerosis (SSc) is a genetically complex autoimmune disease mediated by the interplay between genetic and epigenetic factors in a multitude of immune cells, with CD4+ T lymphocytes as one of the principle drivers of pathogenesis.MethodsDNA samples exacted from CD4+ T cells of 48 SSc patients and 16 healthy controls were hybridized on MethylationEPIC BeadChip array. In parallel, gene expression was interrogated by hybridizing total RNA on Clariom™ S array. Downstream bioinformatics analyses were performed to identify correlating differentially methylated CpG positions (DMPs) and differentially expressed genes (DEGs), which were then confirmed utilizing previously published promoter capture Hi-C (PCHi-C) data.ResultsWe identified 9112 and 3929 DMPs and DEGs, respectively. These DMPs and DEGs are enriched in functional categories related to inflammation and T cell biology. Furthermore, correlation analysis identified 17,500 possible DMP-DEG interaction pairs within a window of 5 Mb, and utilizing PCHi-C data, we observed that 212 CD4+ T cell-specific pairs of DMP-DEG also formed part of three-dimensional promoter-enhancer networks, potentially involving CTCF. Finally, combining PCHi-C data with SSc GWAS data, we identified four important SSc-associated susceptibility loci, TNIP1 (rs3792783), GSDMB (rs9303277), IL12RB1 (rs2305743), and CSK (rs1378942), that could potentially interact with DMP-DEG pairs cg17239269-ANXA6, cg19458020-CCR7, cg10808810-JUND, and cg11062629-ULK3, respectively.ConclusionOur study unveils a potential link between genetic, epigenetic, and transcriptional deregulation in CD4+ T cells of SSc patients, providing a novel integrated view of molecular components driving SSc pathogenesis.

Analysis of chromatin organization and gene expression in T cells identifies functional genes for rheumatoid arthritis

Genome-wide association studies have identified genetic variation contributing to complex disease risk. However, assigning causal genes and mechanisms has been more challenging because disease-associated variants are often found in distal regulatory regions with cell-type specific behaviours. Here, we collect ATAC-seq, Hi-C, Capture Hi-C and nuclear RNA-seq data in stimulated CD4+ T cells over 24 h, to identify functional enhancers regulating gene expression. We characterise changes in DNA interaction and activity dynamics that correlate with changes in gene expression, and find that the strongest correlations are observed within 200 kb of promoters. Using rheumatoid arthritis as an example of T cell mediated disease, we demonstrate interactions of expression quantitative trait loci with target genes, and confirm assigned genes or show complex interactions for 20% of disease associated loci, including FOXO1, which we confirm using CRISPR/Cas9.

A supervised learning framework for chromatin loop detection in genome-wide contact maps

Accurately predicting chromatin loops from genome-wide interaction matrices such as Hi-C data is critical to deepening our understanding of proper gene regulation. Current approaches are mainly focused on searching for statistically enriched dots on a genome-wide map. However, given the availability of orthogonal data types such as ChIA-PET, HiChIP, Capture Hi-C, and high-throughput imaging, a supervised learning approach could facilitate the discovery of a comprehensive set of chromatin interactions. Here, we present Peakachu, a Random Forest classification framework that predicts chromatin loops from genome-wide contact maps. We compare Peakachu with current enrichment-based approaches, and find that Peakachu identifies a unique set of short-range interactions. We show that our models perform well in different platforms, across different sequencing depths, and across different species. We apply this framework to predict chromatin loops in 56 Hi-C datasets, and release the results at the 3D Genome Browser.

Cohesin-Dependent and -Independent Mechanisms Mediate Chromosomal Contacts between Promoters and Enhancers

SummaryIt is currently assumed that 3D chromosomal organization plays a central role in transcriptional control. However, depletion of cohesin and CTCF affects the steady-state levels of only a minority of transcripts. Here, we use high-resolution Capture Hi-C to interrogate the dynamics of chromosomal contacts of all annotated human gene promoters upon degradation of cohesin and CTCF. We show that a majority of promoter-anchored contacts are lost in these conditions, but many contacts with distinct properties are maintained, and some new ones are gained. The rewiring of contacts between promoters and active enhancers upon cohesin degradation associates with rapid changes in target gene transcription as detected by SLAM sequencing (SLAM-seq). These results provide a mechanistic explanation for the limited, but consistent, effects of cohesin and CTCF depletion on steady-state transcription and suggest the existence of both cohesin-dependent and -independent mechanisms of enhancer-promoter pairing.

Promoter anchored interaction landscape of THP-1 macrophages captures early immune response processes

Macrophages are highly plastic immune cells with temporally distinct transcriptome changes upon lipopolysaccride (LPS) activation. However, to what extent transcriptome reprogramming is mediated via spatial chromatin looping is not well studied. We generated high resolution chromatin interaction maps for LPS-stimulated THP-1 macrophages (0 and 2 h) using capture Hi-C. Success of LPS stimulation was validated with transcriptome sequencing. Circa 2900 genes changed their interaction profile upon LPS stimulation and those gaining interactions were enriched for LPS response relevant processes, suggesting a substantial role for distal regulation. Immune and cardiovascular risk variants were enriched within the interacting regions, thereby providing insights into macrophage biology.

Active and repressed biosynthetic gene clusters have spatially distinct chromosome states

While colocalization within a bacterial operon enables coexpression of the constituent genes, the mechanistic logic of clustering of nonhomologous monocistronic genes in eukaryotes is not immediately obvious. Biosynthetic gene clusters that encode pathways for specialized metabolites are an exception to the classical eukaryote rule of random gene location and provide paradigmatic exemplars with which to understand eukaryotic cluster dynamics and regulation. Here, using 3C, Hi-C, and Capture Hi-C (CHi-C) organ-specific chromosome conformation capture techniques along with high-resolution microscopy, we investigate how chromosome topology relates to transcriptional activity of clustered biosynthetic pathway genes in Arabidopsis thaliana Our analyses reveal that biosynthetic gene clusters are embedded in local hot spots of 3D contacts that segregate cluster regions from the surrounding chromosome environment. The spatial conformation of these cluster-associated domains differs between transcriptionally active and silenced clusters. We further show that silenced clusters associate with heterochromatic chromosomal domains toward the periphery of the nucleus, while transcriptionally active clusters relocate away from the nuclear periphery. Examination of chromosome structure at unrelated clusters in maize, rice, and tomato indicates that integration of clustered pathway genes into distinct topological domains is a common feature in plant genomes. Our results shed light on the potential mechanisms that constrain coexpression within clusters of nonhomologous eukaryotic genes and suggest that gene clustering in the one-dimensional chromosome is accompanied by compartmentalization of the 3D chromosome.

242-OR: Rare Genetic Variants in the Islet Regulome Contribute to the Heritability of Type 2 Diabetes Risk

Genome-wide association studies of type 2 diabetes (T2D) show common variant heritability is enriched within gene regulatory regions of pancreatic islets. Whole genome sequence (WGS) data in cohorts allows us to determine the contribution of the complete spectrum of genetic variation to the heritability of T2D risk, including rare alleles that have recently arisen in the human population. A variant-based heritability analysis can partition disease risk into environmental and genetic components, and the genetic components can be further decomposed by variant frequency (a proxy of variant age) and pancreatic islets regulatory annotation (a proxy of functional relevance to T2D risk). We hypothesized that rare, non-protein coding variants would have a significant contribution to the heritability of T2D risk and that further restricting to gene regulatory regions active in pancreatic islets would also produce significant heritability estimates. We conducted a WGS heritability analysis (GREML method in GCTA software) of T2D risk in 15109 unrelated individuals (2215 with T2D; 12894 controls) of European ancestry from the NHLBI’s TOPMed study. We partitioned variants with allele count > 5 by minor allele frequency (MAF): rare variants [MAF<0.1%], and annotation: all non-coding variants, and non-coding variants in pancreatic islets regulatory regions. We further partitioned variants by linkage disequilibrium properties to reduce estimate bias. One of the rare, "all non-coding" variants sets showed a relatively large liability-scale heritability estimate (h2=0.6, P< 0.001). Further significant estimates in rare variants were from active transcription start sites (h2=0.1, P<0.01); active promoters (three variant sets with h2=0.07, 0.08 and 0.1, P<0.05, after body mass index adjustment) and active enhancers from promotor capture Hi-C models (h2=0.07, P<0.05). These findings show a role for rare, non-coding variants in the genetic architecture of type 2 diabetes. Disclosure A. Manning: None. Funding National Institutes of Health (K01DK107836)

Mapping DNA interaction landscapes in psoriasis susceptibility loci highlights KLF4 as a target gene in 9q31

BackgroundGenome-wide association studies (GWAS) have uncovered many genetic risk loci for psoriasis, yet many remain uncharacterised in terms of the causal gene and their biological mechanism in disease. This is largely a result of the findings that over 90% of GWAS variants map outside of protein-coding DNA and instead are enriched in cell type- and stimulation-specific gene regulatory regions.ResultsHere, we use a disease-focused Capture Hi-C (CHi-C) experiment to link psoriasis-associated variants with their target genes in psoriasis-relevant cell lines (HaCaT keratinocytes and My-La CD8+ T cells). We confirm previously assigned genes, suggest novel candidates and provide evidence for complexity at psoriasis GWAS loci. For one locus, uniquely, we combine further epigenomic evidence to demonstrate how a psoriasis-associated region forms a functional interaction with the distant (> 500 kb) KLF4 gene. This interaction occurs between the gene and active enhancers in HaCaT cells, but not in My-La cells. We go on to investigate this long-distance interaction further with Cas9 fusion protein-mediated chromatin modification (CRISPR activation) coupled with RNA-seq, demonstrating how activation of the psoriasis-associated enhancer upregulates KLF4 and its downstream targets, relevant to skin cells and apoptosis.ConclusionsThis approach utilises multiple functional genomic techniques to follow up GWAS-associated variants implicating relevant cell types and causal genes in each locus; these are vital next steps for the translation of genetic findings into clinical benefit.

The “toolbox” for chromatin conformation capture

<p indent=0mm>Studies on the 3D genomics could facilitate the understanding of relationship between genome structure and functions. Based on different research purposes, a variety of chromatin three-dimensional conformation capture methods have been developed, greatly promoting the progress of three-dimensional genomics research. This review aims to systematically introduce and summarize various 3D genomics research methods. The methods are generally divided into three categories. The first category is based on proximity ligation and sequencing, including 3C, Hi-C, <italic>in situ</italic> Hi-C, capture Hi-C, single-cell Hi-C, DLO Hi-C, and ChIA-PET. The other category is not based on proximity ligation, such as ChIA-Drop, SPRITE, and GAM. The last category is based on microscopy, including seqFISH, CRISPRainbow, and SABER amplifies FISH. Finally, the perspective of conformation capture technologies is discussed.

Promoter-anchored chromatin interactions predicted from genetic analysis of epigenomic data

Promoter-anchored chromatin interactions (PAIs) play a pivotal role in transcriptional regulation. Current high-throughput technologies for detecting PAIs, such as promoter capture Hi-C, are not scalable to large cohorts. Here, we present an analytical approach that uses summary-level data from cohort-based DNA methylation (DNAm) quantitative trait locus (mQTL) studies to predict PAIs. Using mQTL data from human peripheral blood (n = 1980), we predict 34,797 PAIs which show strong overlap with the chromatin contacts identified by previous experimental assays. The promoter-interacting DNAm sites are enriched in enhancers or near expression QTLs. Genes whose promoters are involved in PAIs are more actively expressed, and gene pairs with promoter-promoter interactions are enriched for co-expression. Integration of the predicted PAIs with GWAS data highlight interactions among 601 DNAm sites associated with 15 complex traits. This study demonstrates the use of mQTL data to predict PAIs and provides insights into the role of PAIs in complex trait variation.

Identification of cis Elements for Spatio-temporal Control of DNA Replication

The temporal order of DNA replication (replication timing, RT) is highly coupled with genome architecture, but cis-elements regulating spatio-temporal control of replicatio have remained elusive. We performed an extensive series of CRISPR mediated deletions and inversions and high-resolution capture Hi-C of a pluripotency associated domain (DppA2/4) in mouse embryonic stem cells. Whereas CTCF mediated loops and chromatin domain boundaries were dispensable, deletion of three intra-domain prominent CTCF-independent 3D contact sites caused a domain-wide delay in RT, shift in sub-nuclear chromatin compartment and loss of transcriptional activity, These “early replication control elements” (ERCEs) display prominent chromatin features resembling enhancers/promoters and individual and pair-wise deletions of the ERCEs confirmed their partial redundancy and interdependency in controlling domain-wide RT and transcription. Our results demonstrate that discrete cis-regulatory elements mediate domain-wide RT, chromatin compartmentalization, and transcription, representing a major advance in dissecting the relationship between genome structure and function.

Galaxy HiCExplorer 3: a web server for reproducible Hi-C, capture Hi-C and single-cell Hi-C data analysis, quality control and visualization.

The Galaxy HiCExplorer provides a web service at https://hicexplorer.usegalaxy.eu. It enables the integrative analysis of chromosome conformation by providing tools and computational resources to pre-process, analyse and visualize Hi-C, Capture Hi-C (cHi-C) and single-cell Hi-C (scHi-C) data. Since the last publication, Galaxy HiCExplorer has been expanded considerably with new tools to facilitate the analysis of cHi-C and to provide an in-depth analysis of Hi-C data. Moreover, it supports the analysis of scHi-C data by offering a broad range of tools. With the help of the standard graphical user interface of Galaxy, presented workflows, extensive documentation and tutorials, novices as well as Hi-C experts are supported in their Hi-C data analysis with Galaxy HiCExplorer.