High-throughput Chromatin Conformation Capture Research Articles

Modification of the hi-c technology for molecular genetic analysis of formalin-fixed paraffin-embedded sections of tumor tissues

Molecular genetic analysis of tumor tissues is the most important step towards understanding the mechanisms of cancer development, and it is also necessary for the choice of targeted therapy. Hi-C (high-throughput chromatin conformation capture) technology can be used to detect various types of genomic variants, including balanced chromosomal rearrangements: inversions and translocations. In this work, we propose a modification of the Hi-C method for the analysis of chromatin contacts in formalin-fixed paraffin-embedded (FFPE) sections of tumor tissues. Our results demonstrate that this protocol allow to generate high-quality Hi-C data and detect all types of chromosomal rearrangements. We have analyzed various databases to compile a comprehensive list of translocations that hold clinical importance for targeted therapy selection. The practical value of molecular genetic testing is its ability to influence patient treatment strategies and provide prognostic insights. Detecting specific chromosomal rearrangements can guide the choice of targeted therapies, which is a critical aspect of personalized medicine in oncology.

Snow alga Sanguina aurantia as revealed through de novo genome assembly and annotation.

To thrive on melting alpine and polar snow, some Chlorophytes produce an abundance of astaxanthin, causing red blooms, often dominated by genus Sanguina. The red cells have not been cultured, but we recently grew a green biciliate conspecific with Sanguina aurantia from a sample of watermelon snow. This culture provided source material for Oxford Nanopore Technology and Illumina sequencing. Our assembly pipeline exemplifies the value of a hybrid long- and short-read approach for the complexities of working with a culture grown from a field sample. Using bioinformatic tools, we separated assembled contigs into 2 genomic pools based on a difference in GC content (57.5 and 55.1%). We present the data as 2 assemblies of S. aurantia variants but explore other possibilities. High-throughput chromatin conformation capture analysis (Hi-C sequencing) was used to scaffold the assemblies into a 96-Mb genome designated as "A" and a 102-Mb genome designated as "B." Both assemblies are highly contiguous: genome A consists of 38 scaffolds with an N50 of 5.4 Mb, while genome B has 50 scaffolds with an N50 of 6.4 Mb. RNA sequencing was used to improve gene annotation.

Genome Restructuring around Innate Immune Genes in Monocytes in Alcohol-associated Hepatitis.

Many inflammatory genes in the immune system are clustered in the genome. The 3D genome architecture of these clustered genes likely plays a critical role in their regulation and alterations to this structure may contribute to diseases where inflammation is poorly controlled. Alcohol-associated hepatitis (AH) is a severe inflammatory disease that contributes significantly to morbidity in alcohol associated liver disease. Monocytes in AH are hyper-responsive to inflammatory stimuli and contribute significantly to inflammation. We performed high throughput chromatin conformation capture (Hi-C) technology on monocytes isolated from 4 AH patients and 4 healthy controls to better understand how genome structure is altered in AH. Most chromosomes from AH and healthy controls were significantly dissimilar from each other. Comparing AH to HC, many regions of the genome contained significant changes in contact frequency. While there were alterations throughout the genome, there were a number of hotspots containing a higher density of changes in structure. A few of these hotspots contained genes involved in innate immunity including the NK-gene receptor complex and the CXC-chemokines. Finally, we compare these results to scRNA-seq data from patients with AH challenged with LPS to predict how chromatin conformation impacts transcription of clustered immune genes. Together, these results reveal changes in the chromatin structure of monocytes from AH patients that perturb expression of highly clustered proinflammatory genes.

Chromosome-level scaffolding of haplotype-resolved assemblies using Hi-C data without reference genomes.

Scaffolding is crucial for constructing most chromosome-level genomes. The high-throughput chromatin conformation capture (Hi-C) technology has become the primary scaffolding strategy due to its convenience and cost-effectiveness. As sequencing technologies and assembly algorithms advance, constructing haplotype-resolved genomes is increasingly preferred because haplotypes can provide additional genetic information on allelic and non-allelic variations. ALLHiC is a widely used allele-aware scaffolding tool designed for this purpose. However, its dependence on chromosome-level reference genomes and a higher chromosome misassignment rate still impede the unravelling of haplotype-resolved genomes. Here we present HapHiC, a reference-independent allele-aware scaffolding tool with superior performance on chromosome assignment as well as contig ordering and orientation. In addition, we provide new insights into the challenges in allele-aware scaffolding by conducting comprehensive analyses on various adverse factors. Finally, with the help of HapHiC, we constructed the haplotype-resolved allotriploid genome for Miscanthus × giganteus, an important lignocellulosic bioenergy crop.

Chromosome-scale reference genome of an ancient landrace: unveiling the genetic basis of seed weight in the food legume crop pigeonpea (Cajanus cajan).

Pigeonpea (Cajanus cajan) is a nutrient-rich and versatile food legume crop of tropical and subtropical regions. In this study, we describe the de novo assembly of a high-quality genome for the ancient pigeonpea landrace 'D30', achieved through a combination of Pacific Biosciences high-fidelity (PacBio HiFi) and high-throughput chromatin conformation capture (Hi-C) sequencing technologies. The assembled 'D30' genome has a size of 813.54Mb, with a contig N50 of 10.74Mb, a scaffold N50 of 73.07Mb, and a GC content of 35.67%. Genomic evaluation revealed that the 'D30' genome contains 99.2% of Benchmarking Universal Single-Copy Orthologs (BUSCO) and achieves a 29.06 long terminal repeat (LTR) assembly index (LAI). Genome annotation indicated that 'D30' encompasses 431.37Mb of repeat elements (53.02% of the genome) and 37 977 protein-coding genes. Identification of single-nucleotide polymorphisms (SNPs), insertions/deletions (indels), and structural variations between 'D30' and the published genome of pigeonpea cultivar 'Asha' suggests that genes affected by these variations may play important roles in biotic and abiotic stress responses. Further investigation of genomic regions under selection highlights genes enriched in starch and sucrose metabolism, with 42.11% of these genes highly expressed in seeds. Finally, we conducted genome-wide association studies (GWAS) to facilitate the identification of 28 marker-trait associations for six agronomic traits of pigeonpea. Notably, we discovered a calmodulin-like protein (CcCML) that harbors a dominant haplotype associated with the 100-seed weight of pigeonpea. Our study provides a foundational resource for developing genomics-assisted breeding programs in pigeonpea.

Identifying topologically associating domains using differential kernels.

Chromatin is a polymer complex of DNA and proteins that regulates gene expression. The three-dimensional (3D) structure and organization of chromatin controls DNA transcription and replication. High-throughput chromatin conformation capture techniques generate Hi-C maps that can provide insight into the 3D structure of chromatin. Hi-C maps can be represented as a symmetric matrix [Formula: see text], where each element represents the average contact probability or number of contacts between chromatin loci i and j. Previous studies have detected topologically associating domains (TADs), or self-interacting regions in [Formula: see text] within which the contact probability is greater than that outside the region. Many algorithms have been developed to identify TADs within Hi-C maps. However, most TAD identification algorithms are unable to identify nested or overlapping TADs and for a given Hi-C map there is significant variation in the location and number of TADs identified by different methods. We develop a novel method to identify TADs, KerTAD, using a kernel-based technique from computer vision and image processing that is able to accurately identify nested and overlapping TADs. We benchmark this method against state-of-the-art TAD identification methods on both synthetic and experimental data sets. We find that the new method consistently has higher true positive rates (TPR) and lower false discovery rates (FDR) than all tested methods for both synthetic and manually annotated experimental Hi-C maps. The TPR for KerTAD is also largely insensitive to increasing noise and sparsity, in contrast to the other methods. We also find that KerTAD is consistent in the number and size of TADs identified across replicate experimental Hi-C maps for several organisms. Thus, KerTAD will improve automated TAD identification and enable researchers to better correlate changes in TADs to biological phenomena, such as enhancer-promoter interactions and disease states.

An identification method of chromatin topological associated domains based on spatial density clustering

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.

Transcriptional Hubs Within Cliques in Ensemble Hi-C Chromatin Interaction Networks.

Chromatin conformation capture technologies permit the study of chromatin spatial organization on a genome-wide scale at a variety of resolutions. Despite the increasing precision and resolution of high-throughput chromatin conformation capture (Hi-C) methods, it remains challenging to conclusively link transcriptional activity to spatial organizational phenomena. We have developed a clique-based approach for analyzing Hi-C data that helps identify chromosomal hotspots that feature considerable enrichment of chromatin annotations for transcriptional start sites and, building on previously published work, show that these chromosomal hotspots are not only significantly enriched in RNA polymerase II binding sites as identified by the ENCODE project, but also identify a noticeable increase in FANTOM5 and GTEx transcription within our identified cliques across a variety of tissue types. From the obtained data, we surmise that our cliques are a suitable method for identifying transcription factories in Hi-C data, and outline further extensions to the method that may make it useful for locating regions of increased transcriptional activity in datasets where in-depth expression or polymerase data may not be available.

First telomere-to-telomere gapless assembly of the rice blast fungus Pyricularia oryzae

Rice blast caused by Pyricularia oryzae (syn., Magnaporthe oryzae) was one of the most destructive diseases of rice throughout the world. Genome assembly was fundamental to genetic variation identification and critically impacted the understanding of its ability to overcome host resistance. Here, we report a gapless genome assembly of rice blast fungus P. oryzae strain P131 using PacBio, Illumina and high throughput chromatin conformation capture (Hi-C) sequencing data. This assembly contained seven complete chromosomes (43,237,743 bp) and a circular mitochondrial genome (34,866 bp). Approximately 14.31% of this assembly carried repeat sequences, significantly greater than its previous assembled version. This assembly had a 99.9% complement in BUSCO evaluation. A total of 14,982 genes protein-coding genes were predicted. In summary, we assembled the first telomere-to-telomere gapless genome of P. oryzae, which would be a valuable genome resource for future research on the genome evolution and host adaptation.

Modification of the Hi-C Technology for Molecular Genetic Analysis of Formalin-Fixed Paraffin-Embedded Sections of Tumor Tissues.

Molecular genetic analysis of tumor tissues is the most important step towards understanding the mechanisms of cancer development; it is also necessary for the choice of targeted therapy. The Hi-C (high-throughput chromatin conformation capture) technology can be used to detect various types of genomic variants, including balanced chromosomal rearrangements, such as inversions and translocations. We propose a modification of the Hi-C method for the analysis of chromatin contacts in formalin-fixed paraffin-embedded (FFPE) sections of tumor tissues. The developed protocol allows to generate high-quality Hi-C data and detect all types of chromosomal rearrangements. We have analyzed various databases to compile a comprehensive list of translocations that hold clinical importance for the targeted therapy selection. The practical value of molecular genetic testing is its ability to influence the treatment strategies and to provide prognostic insights. Detecting specific chromosomal rearrangements can guide the choice of the targeted therapies, which is a critical aspect of personalized medicine in oncology.

BL-Hi-C reveals the 3D genome structure of Brassica crops with high sensitivity.

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.

Chromosome-level genome assembly of the Stoliczka’s Asian trident bat (Aselliscus stoliczkanus)

Stoliczka’s Asian trident bat (Aselliscus stoliczkanus) is a small-bodied species and very sensitive to climate change. Here, we presented a chromosome-level genome assembly of A. stoliczkanus by combining Illumina sequencing, Nanopore sequencing and high-throughput chromatin conformation capture (Hi-C) sequencing technology. The genome assembly was 2.18 Gb in size with 98.26% of the genome sequences anchored onto 14 autosomes and two sex chromosomes (X and Y). The quality of the genome assembly is very high with a contig and scaffold N50 of 72.98 and 162 Mb, respectively, Benchmarking Universal Single-Copy Orthologs (BUSCO) score of 96.6%, and the consensus quality value (QV) of 47.44. A total of 20,567 genes were predicted and 98.8% of these genes were functionally annotated. Syntenic blocks between A. stoliczkanus and Homo sapiens, together with previous comparative cytogenetic studies, provide valuable foundations for further comparative genomic and cytogenetic studies in mammals. The reference-quality genome of A. stoliczkanus contributes an important resource for conservative genomics and landscape genomics in predicting adaptation and vulnerability to climate change.

Hi-C sequencing unravels dynamic three-dimensional chromatin interactions in muntjac lineage: insights from chromosome fusions in Fea's muntjac genome.

Eukaryotes have varying numbers and structures of characteristic chromosomes across lineages or species. The evolutionary trajectory of species may have been affected by spontaneous genome rearrangements. Chromosome fusion drastically alters karyotypes. However, the mechanisms and consequences of chromosome fusions, particularly in muntjac species, are poorly understood. Recent research-based advancements in three-dimensional (3D) genomics, particularly high-throughput chromatin conformation capture (Hi-C) sequencing, have allowed for the identification of chromosome fusions and provided mechanistic insights into three muntjac species: Muntiacus muntjak, M. reevesi, and M. crinifrons. This study aimed to uncover potential genome rearrangement patterns in the threatened species Fea's muntjac (Muntiacus feae), which have not been previously examined for such characteristics. Deep Hi-C sequencing (31.42 × coverage) was performed to reveal the 3D chromatin architecture of the Fea's muntjac genome. Patterns of repeated chromosome fusions that were potentially mediated by high-abundance transposable elements were identified. Comparative Hi-C maps demonstrated linkage homology between the sex chromosomes in Fea's muntjac and autosomes in M. reevesi, indicating that fusions may have played a crucial role in the evolution of the sex chromosomes of the lineage. The species-level dynamics of topologically associated domains (TADs) suggest that TAD organization could be altered by differential chromosome interactions owing to repeated chromosome fusions. However, research on the effect of TADs on muntjac genome evolution is insufficient. This study generated Hi-C data for the Fea's muntjac, providing a genomic resource for future investigations of the evolutionary patterns of chromatin conformation at the chromosomal level.

Accessible gene borders establish a core structural unit for chromatin architecture in Arabidopsis.

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.

SnapHiC-D: a computational pipeline to identify differential chromatin contacts from single-cell Hi-C data.

Single-cell high-throughput chromatin conformation capture technologies (scHi-C) has been used to map chromatin spatial organization in complex tissues. However, computational tools to detect differential chromatin contacts (DCCs) from scHi-C datasets in development and through disease pathogenesis are still lacking. Here, we present SnapHiC-D, a computational pipeline to identify DCCs between two scHi-C datasets. Compared to methods designed for bulk Hi-C data, SnapHiC-D detects DCCs with high sensitivity and accuracy. We used SnapHiC-D to identify cell-type-specific chromatin contacts at 10 Kb resolution in mouse hippocampal and human prefrontal cortical tissues, demonstrating that DCCs detected in the hippocampal and cortical cell types are generally associated with cell-type-specific gene expression patterns and epigenomic features. SnapHiC-D is freely available at https://github.com/HuMingLab/SnapHiC-D.

HiBrowser: an interactive and dynamic browser for synchronous Hi-C data visualization.

With the development of chromosome conformation capture technology, the genome-wide investigation of higher-order chromatin structure by using high-throughput chromatin conformation capture (Hi-C) technology is emerging as an important component for understanding the mechanism of gene regulation. Considering genetic and epigenetic differences are typically used to explore the pathological reasons on the chromosome and gene level, visualizing multi-omics data and performing an intuitive analysis by using an interactive browser become a powerful and welcomed way. In this paper, we develop an effective sequence and chromatin interaction data display browser called HiBrowser for visualizing and analyzing Hi-C data and their associated genetic and epigenetic annotations. The advantages of HiBrowser are flexible multi-omics navigation, novel multidimensional synchronization comparisons and dynamic interaction system. In particular, HiBrowser first provides an out of the box web service and allows flexible and dynamic reconstruction of custom annotation tracks on demand during running. In order to conveniently and intuitively analyze the similarities and differences among multiple samples, such as visual comparisons of normal and tumor tissue samples, and pan genomes of multiple (consanguineous) species, HiBrowser develops a clone mode to synchronously display the genome coordinate positions or the same regions of multiple samples on the same page of visualization. HiBrowser also supports a pluralistic and precise search on correlation data of distal cis-regulatory elements and navigation to any region on Hi-C heatmap of interest according to the searching results. HiBrowser is a no-build tool, and could be easily deployed in local server. The source code is available at https://github.com/lyotvincent/HiBrowser.

3D genome structural variations play important roles in regulating seed oil content of Brassica napus

Dissecting the complex regulatory mechanism of seed oil content (SOC) is one of the main research goals in Brassica napus. Increasing evidence suggests that genome architecture is linked to multiple biological functions. However, the effect of genome architecture on SOC regulation remains unclear. Here, we used high-throughput chromatin conformation capture to characterize differences in the three-dimensional (3D) landscape of genome architecture of seeds from two B. napus lines, N53-2 (with high SOC) and Ken-C8 (with low SOC). Bioinformatics analysis demonstrated that differentially accessible regions and differentially expressed genes between N53-2 and Ken-C8 were preferentially enriched in regions with quantitative trait loci (QTLs)/associated genomic regions (AGRs) for SOC. A multi-omics analysis demonstrated that expression of SOC-related genes was tightly correlated with genome structural variations in QTLs/AGRs of B. napus. The candidate gene BnaA09g48250D, which showed structural variation in a QTL/AGR on chrA09, was identified by fine-mapping of a KN double-haploid population derived from hybridization of N53-2 and Ken-C8. Overexpression and knockout of BnaA09g48250D led to significant increases and decreases in SOC, respectively, in the transgenic lines. Taken together, our results reveal the 3D genome architecture of B. napus seeds and the roles of genome structural variations in SOC regulation, enriching our understanding of the molecular mechanisms of SOC regulation from the perspective of spatial chromatin structure.

Spatio-temporal transcriptome dynamics coordinate rapid transition of core crop functions in 'lactating' pigeon.

Pigeons (Columba livia) are among a select few avian species that have developed a specialized reproductive mode wherein the parents produce a 'milk' in their crop to feed newborn squabs. Nonetheless, the transcriptomic dynamics and role in the rapid transition of core crop functions during 'lactation' remain largely unexplored. Here, we generated a de novo pigeon genome assembly to construct a high resolution spatio-temporal transcriptomic landscape of the crop epithelium across the entire breeding stage. This multi-omics analysis identified a set of 'lactation'-related genes involved in lipid and protein metabolism, which contribute to the rapid functional transitions in the crop. Analysis of in situ high-throughput chromatin conformation capture (Hi-C) sequencing revealed extensive reorganization of promoter-enhancer interactions linked to the dynamic expression of these 'lactation'-related genes between stages. Moreover, their expression is spatially localized in specific epithelial layers, and can be correlated with phenotypic changes in the crop. These results illustrate the preferential de novo synthesis of 'milk' lipids and proteins in the crop, and provides candidate enhancer loci for further investigation of the regulatory elements controlling pigeon 'lactation'.

Abstract A054: Accurate prediction of cohesin-mediated 3D genome organization using 2D chromatin features

Abstract The three-dimensional (3D) genome organization influences diverse nuclear processes. Chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) and High-throughput chromatin conformation capture (Hi-C) are powerful methods used to study the 3D genome organization. However, these two experiments are costly, time-consuming, require tens to hundreds of millions of cells, and challenging to optimize and analyze. Predicting ChIA-PET/Hi-C data using cheaper ChIP-Seq data and other easily obtainable features is a useful alternative. Also, it is well-established that the cohesin protein complex is a key determinant of 3D genome organization. Considering this, we present Chromatin Interaction Predictor (ChIPr), a suite of regression models based on deep neural networks (DNN), random forest, and gradient boosting. ChIPr predicts cohesin-mediated chromatin interaction strength between any two given loci in the genome. It was trained mainly using two top-tiered ENCODE cell lines GM12878 and K562 and tested on other cell lines (H1 and HepG2). Comprehensive tests showed that ChIPr predictions correlated well with the original ChIA-PET data at the peak-level resolution and at bin sizes 25kb and 5kb. In addition, ChIPr accurately captured most of the cell-type-dependent loops identified by both, ChIA-PET and Hi-C data, respectively. Extensive feature testing highlighted genomic distance and RAD21 (a cohesin component) ChIP-Seq signals as the most important inputs contributing to ChIPr, in determining the chromatin interaction strength. A standard ChIPr model requires three experimental inputs: ChIP-Seq signals for RAD21, H3K27ac (enhancer/active chromatin mark), and H3K27me3 (inactive chromatin mark). While a reduced model requires a single experiment input: ChIP-Seq signals for RAD21 and performs equally well. Additionally, the use of ChIPr was further extended to predict the contact maps for several prostate cancer (PCa) cell lines whose ChIA-PET data were not available on the ENCODE portal, for instance, RWPE1 and VCaP. We generated the necessary ChIP-seq data for these cell lines to use as input to the model to further examine the changes in the 3D interactions of PCa driver genes and also their isolated neighbourhoods and decipher useful biological insights about the gene regulatory networks involved. Integrative analysis revealed novel insights into the role of CTCF motif, its orientation, and CTCF binding on the prevalence and strength of cohesin-mediated chromatin interactions (~50% - 80% of the RAD21 interactions were enriched with CTCF, < 15% of the interactions had no CTCF ChIP-seq binding in both of the two peaks). Lastly, we also noticed a subset of RAD21 interactions with CTCF binding in only one or none of the two anchor peaks, which were significantly enriched with enhancers. Taken together, this study outlines the general features of genome folding and opens new avenues to analyze spatial genome organization in specimens with limited cell numbers. Citation Format: Khyati Chandratre. Accurate prediction of cohesin-mediated 3D genome organization using 2D chromatin features [abstract]. In: Proceedings of the AACR Special Conference: Advances in Prostate Cancer Research; 2023 Mar 15-18; Denver, Colorado. Philadelphia (PA): AACR; Cancer Res 2023;83(11 Suppl):Abstract nr A054.

Conserved noncoding sequences correlate with distant gene contacts in Arabidopsis and Brassica.

Physical contact between genes distant on chromosomes is a potentially important way for genes to coordinate their expressions. To investigate the potential importance of distant contacts, we performed high-throughput chromatin conformation capture (Hi-C) experiments on leaf nuclei isolated from Brassica rapa and Brassica oleracea. We then combined our results with published Hi-C data from Arabidopsis thaliana. We found that distant genes come into physical contact and do so preferentially between the proximal promoter of one gene and the downstream region of another gene. Genes with higher numbers of conserved noncoding sequences (CNSs) nearby were more likely to have contact with distant genes. With more CNSs came higher numbers of transcription factor binding sites and more histone modifications associated with the activity. In addition, for the genes we studied, distant contacting genes with CNSs were more likely to be transcriptionally coordinated. These observations suggest that CNSs may enrich active histone modifications and recruit transcription factors, correlating with distant contacts to ensure coordinated expression. This study advances our knowledge of gene contacts and provides insights into the relationship between CNSs and distant gene contacts in plants.