Fine-scale chromatin interaction maps reveal the cis-regulatory landscape of human lincRNA genes

Evolutionary analysis of gene ages across TADs associates chromatin topology with whole-genome duplications

Topologically associated domains (TADs) are interaction subnetworks of chromosomal regions in 3D genomes. TAD boundaries frequently coincide with genome breaks while boundary deletion is under negative selection, suggesting that TADs may facilitate genome rearrangements and evolution. We show that genes co-localize by evolutionary age in humans and mice, resulting in TADs having different proportions of younger and older genes. We observe a major transition in the age co-localization patterns between the genes born during vertebrate whole-genome duplications (WGDs) or before and those born afterward. We also find that genes recently duplicated in primates and rodents are more frequently essential when they are located in old-enriched TADs and interact with genes that last duplicated during the WGD. Therefore, the evolutionary relevance of recent genes may increase when located in TADs with established regulatory networks. Our data suggest that TADs could play a role in organizing ancestral functions and evolutionary novelty.

Hi-C, a chromatin 3D structure technique advancing the functional genomics of immune cells

The functional performance of immune cells relies on a complex transcriptional regulatory network. The three-dimensional structure of chromatin can affect chromatin status and gene expression patterns, and plays an important regulatory role in gene transcription. Currently available techniques for studying chromatin spatial structure include chromatin conformation capture techniques and their derivatives, chromatin accessibility sequencing techniques, and others. Additionally, the recently emerged deep learning technology can be utilized as a tool to enhance the analysis of data. In this review, we elucidate the definition and significance of the three-dimensional chromatin structure, summarize the technologies available for studying it, and describe the research progress on the chromatin spatial structure of dendritic cells, macrophages, T cells, B cells, and neutrophils.

Enhancer mutations modulate the severity of chemotherapy-induced myelosuppression

Non-small cell lung cancer is often diagnosed at advanced stages, and many patients are still treated with classical chemotherapy. The unselective nature of chemotherapy often results in severe myelosuppression. Previous studies showed that protein-coding mutations could not fully explain the predisposition to myelosuppression. Here, we investigate the possible role of enhancer mutations in myelosuppression susceptibility. We produced transcriptome and promoter-interaction maps (using HiCap) of three blood stem-like cell lines treated with carboplatin or gemcitabine. Taking advantage of publicly available enhancer datasets, we validated HiCap results in silico and in living cells using epigenetic CRISPR technology. We also developed a network approach for interactome analysis and detection of differentially interacting genes. Differential interaction analysis provided additional information on relevant genes and pathways for myelosuppression compared with differential gene expression analysis at the bulk level. Moreover, we showed that enhancers of differentially interacting genes are highly enriched for variants associated with differing levels of myelosuppression. Altogether, our work represents a prominent example of integrative transcriptome and gene regulatory datasets analysis for the functional annotation of noncoding mutations.

Cephalopod-omics: Emerging Fields and Technologies in Cephalopod Biology

Few animal groups can claim the level of wonder that cephalopods instill in the minds of researchers and the general public. Much of cephalopod biology, however, remains unexplored: the largest invertebrate brain, difficult husbandry conditions, and complex (meta-)genomes, among many other things, have hindered progress in addressing key questions. However, recent technological advancements in sequencing, imaging, and genetic manipulation have opened new avenues for exploring the biology of these extraordinary animals. The cephalopod molecular biology community is thus experiencing a large influx of researchers, emerging from different fields, accelerating the pace of research in this clade. In the first post-pandemic event at the Cephalopod International Advisory Council (CIAC) conference in April 2022, over 40 participants from all over the world met and discussed key challenges and perspectives for current cephalopod molecular biology and evolution. Our particular focus was on the fields of comparative and regulatory genomics, gene manipulation, single-cell transcriptomics, metagenomics, and microbial interactions. This article is a result of this joint effort, summarizing the latest insights from these emerging fields, their bottlenecks, and potential solutions. The article highlights the interdisciplinary nature of the cephalopod-omics community and provides an emphasis on continuous consolidation of efforts and collaboration in this rapidly evolving field.

3D genome organization and epigenetic regulation in autoimmune diseases

Three-dimensional (3D) genomics is an emerging field of research that investigates the relationship between gene regulatory function and the spatial structure of chromatin. Chromatin folding can be studied using chromosome conformation capture (3C) technology and 3C-based derivative sequencing technologies, including chromosome conformation capture-on-chip (4C), chromosome conformation capture carbon copy (5C), and high-throughput chromosome conformation capture (Hi-C), which allow scientists to capture 3D conformations from a single site to the entire genome. A comprehensive analysis of the relationships between various regulatory components and gene function also requires the integration of multi-omics data such as genomics, transcriptomics, and epigenomics. 3D genome folding is involved in immune cell differentiation, activation, and dysfunction and participates in a wide range of diseases, including autoimmune diseases. We describe hierarchical 3D chromatin organization in this review and conclude with characteristics of C-techniques and multi-omics applications of the 3D genome. In addition, we describe the relationship between 3D genome structure and the differentiation and maturation of immune cells and address how changes in chromosome folding contribute to autoimmune diseases.

Determining chromatin architecture with Micro Capture-C

Micro Capture-C (MCC) is a chromatin conformation capture (3C) method for visualizing reproducible three-dimensional contacts of specified regions of the genome at base pair resolution. These methods are an established family of techniques that use proximity ligation to assay the topology of chromatin. MCC can generate data at substantially higher resolution than previous techniques through multiple refinements of the 3C method. Using a sequence agnostic nuclease, the maintenance of cellular integrity and full sequencing of the ligation junctions, MCC achieves subnucleosomal levels of resolution, which can be used to reveal transcription factor binding sites analogous to DNAse I footprinting. Gene dense regions, close-range enhancer-promoter contacts, individual enhancers within super-enhancers and multiple other types of loci or regulatory regions that were previously challenging to assay with conventional 3C techniques, are readily observed using MCC. MCC requires training in common molecular biology techniques and bioinformatics to perform the experiment and analyze the data. The protocol can be expected to be completed in a 3 week timeframe for experienced molecular biologists.

Characterization of Gene Regulatory Elements in Human Fetal Cortical Development: Enhancing Our Understanding of Neurodevelopmental Disorders and Evolution

The neocortex is the region that most distinguishes human brain from other mammals and primates [Annu Rev Genet. 2021 Nov;55(1):555-81]. Studying the development of human cortex is important in understanding the evolutionary changes occurring in humans relative to other primates, as well as in elucidating mechanisms underlying neurodevelopmental disorders. Cortical development is a highly regulated process, spatially and temporally coordinated by expression of essential transcriptional factors in response to signaling pathways [Neuron. 2019 Sep;103(6):980-1004]. Enhancers are the most well-understood cis-acting, non-protein-coding regulatory elements that regulate gene expression [Nat Rev Genet. 2014 Apr;15(4):272-86]. Importantly, given the conservation of both DNA sequence and molecular function of the majority of proteins across mammals [Genome Res. 2003 Dec;13(12):2507-18], enhancers [Science. 2015 Mar;347(6226):1155-9], which are far more divergent at the sequence level, likely account for the phenotypes that distinguish the human brain by changing the regulation of gene expression. In this review, we will revisit the conceptual framework of gene regulation during human brain development, as well as the evolution of technologies to study transcriptional regulation, with recent advances in genome biology that open a window allowing us to systematically characterize cis-regulatory elements in developing human brain [Hum Mol Genet. 2022 Oct;31(R1):R84-96]. We provide an update on work to characterize the suite of all enhancers in the developing human brain and the implications for understanding neuropsychiatric disorders. Finally, we discuss emerging therapeutic ideas that utilize our emerging knowledge of enhancer function.

Predicting enhancer-promoter interaction based on epigenomic signals

Introduction: The physical interactions between enhancers and promoters are often involved in gene transcriptional regulation. High tissue-specific enhancer-promoter interactions (EPIs) are responsible for the differential expression of genes. Experimental methods are time-consuming and labor-intensive in measuring EPIs. An alternative approach, machine learning, has been widely used to predict EPIs. However, most existing machine learning methods require a large number of functional genomic and epigenomic features as input, which limits the application to different cell lines. Methods: In this paper, we developed a random forest model, HARD (H3K27ac, ATAC-seq, RAD21, and Distance), to predict EPI using only four types of features. Results: Independent tests on a benchmark dataset showed that HARD outperforms other models with the fewest features. Discussion: Our results revealed that chromatin accessibility and the binding of cohesin are important for cell-line-specific EPIs. Furthermore, we trained the HARD model in the GM12878 cell line and performed testing in the HeLa cell line. The cross-cell-lines prediction also performs well, suggesting it has the potential to be applied to other cell lines.

3D genomics and its applications in precision medicine

Three-dimensional (3D) genomics is an emerging discipline that studies the three-dimensional structure of chromatin and the three-dimensional and functions of genomes. It mainly studies the three-dimensional conformation and functional regulation of intranuclear genomes, such as DNA replication, DNA recombination, genome folding, gene expression regulation, transcription factor regulation mechanism, and the maintenance of three-dimensional conformation of genomes. Self-chromosomal conformation capture (3C) technology has been developed, and 3D genomics and related fields have developed rapidly. In addition, chromatin interaction analysis techniques developed by 3C technologies, such as paired-end tag sequencing (ChIA-PET) and whole-genome chromosome conformation capture (Hi-C), enable scientists to further study the relationship between chromatin conformation and gene regulation in different species. Thus, the spatial conformation of plant, animal, and microbial genomes, transcriptional regulation mechanisms, interaction patterns of chromosomes, and the formation mechanism of spatiotemporal specificity of genomes are revealed. With the help of new experimental technologies, the identification of key genes and signal pathways related to life activities and diseases is sustaining the rapid development of life science, agriculture, and medicine. In this paper, the concept and development of 3D genomics and its application in agricultural science, life science, and medicine are introduced, which provides a theoretical basis for the study of biological life processes.

Analysis of Chromatin Interaction and Accessibility by Trac-Looping

Spatial organization of the genome modulates pivotal biological processes. The emerging new technologies have provided novel insights into genome structure and its role in regulating cell activities. To examine the genome-wide chromatin interactions at accessible chromatin regions, we developed a DNA transposase-mediated analysis of chromatin looping (Trac-looping) method for simultaneously detecting chromatin interactions and chromatin accessibility. Here, we describe a detailed protocol of generating Trac-looping libraries.

ChIATAC is an efficient strategy for multi-omics mapping of 3D epigenomes from low-cell inputs

Connecting genes to their cis-regulatory elements has been enabled by genome-wide mapping of chromatin interactions using proximity ligation in ChIA-PET, Hi-C, and their derivatives. However, these methods require millions of input cells for high-quality data and thus are unsuitable for many studies when only limited cells are available. Conversely, epigenomic profiling via transposase digestion in ATAC-seq requires only hundreds to thousands of cells to robustly map open chromatin associated with transcription activity, but it cannot directly connect active genes to their distal enhancers. Here, we combine proximity ligation in ChIA-PET and transposase accessibility in ATAC-seq into ChIATAC to efficiently map interactions between open chromatin loci in low numbers of input cells. We validate ChIATAC in Drosophila cells and optimize it for mapping 3D epigenomes in human cells robustly. Applying ChIATAC to primary human T cells, we reveal mechanisms that topologically regulate transcriptional programs during T cell activation.

A Comparison of Topologically Associating Domain Callers Based on Hi-C Data

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.

When 3D genome technology meets viral infection, including SARS-CoV-2

Mammalian chromosomes undergo varying degrees of compression to form three-dimensional genome structures. These three-dimensional structures undergo dynamic and precise chromatin interactions to achieve precise spatial and temporal regulation of gene expression. Most eukaryotic DNA viruses can invade their genomes into the nucleus. However, it is still poorly understood how the viral genome is precisely positioned after entering the host cell nucleus to find the most suitable location and whether it can specifically interact with the host genome to hijack the host transcriptional factories or even integrate into the host genome to complete its transcription and replication rapidly. Chromosome conformation capture technology can reveal long-range chromatin interactions between different chromosomal sites in the nucleus, potentially providing a reference for viral DNA-host chromatin interactions. This review summarized the research progress on the three-dimensional interaction between virus and host genome and the impact of virus integration into the host genome on gene transcription regulation, aiming to provide new insights into chromatin interaction and viral gene transcription regulation, laying the foundation for the treatment of infectious diseases.

The evolution and diversification of oakleaf butterflies

Oakleaf butterflies in the genus Kallima have a polymorphic wing phenotype, enabling these insects to masquerade as dead leaves. This iconic example of protective resemblance provides an interesting evolutionary paradigm that can be employed to study biodiversity. We integrated multi-omic data analyses and functional validation to infer the evolutionary history of Kallima species and investigate the genetic basis of their variable leaf wing patterns. We find that Kallima butterflies diversified in the eastern Himalayas and dispersed to East and Southeast Asia. Moreover, we find that leaf wing polymorphism is controlled by the wing patterning gene cortex, which has been maintained in Kallima by long-term balancing selection. Our results provide macroevolutionary and microevolutionary insights into a model species originating from a mountain ecosystem.

scHiCEmbed: Bin-Specific Embeddings of Single-Cell Hi-C Data Using Graph Auto-Encoders

Most publicly accessible single-cell Hi-C data are sparse and cannot reach a higher resolution. Therefore, learning latent representations (bin-specific embeddings) of sparse single-cell Hi-C matrices would provide us with a novel way of mining valuable information hidden in the limited number of single-cell Hi-C contacts. We present scHiCEmbed, an unsupervised computational method for learning bin-specific embeddings of single-cell Hi-C data, and the computational system is applied to the tasks of 3D structure reconstruction of whole genomes and detection of topologically associating domains (TAD). The only input of scHiCEmbed is a raw or scHiCluster-imputed single-cell Hi-C matrix. The main process of scHiCEmbed is to embed each node/bin in a higher dimensional space using graph auto-encoders. The learned n-by-3 bin-specific embedding/latent matrix is considered the final reconstructed 3D genome structure. For TAD detection, we use constrained hierarchical clustering on the latent matrix to classify bins: S_Dbw is used to determine the optimal number of clusters, and each cluster is considered as one potential TAD. Our reconstructed 3D structures for individual chromatins at different cell stages reveal the expanding process of chromatins during the cell cycle. We observe that the TADs called from single-cell Hi-C data are not shared across individual cells and that the TAD boundaries called from raw or imputed single-cell Hi-C are significantly different from those called from bulk Hi-C, confirming the cell-to-cell variability in terms of TAD definitions. The source code for scHiCEmbed is publicly available, and the URL can be found in the conclusion section.

A comparison of topologically associating domain callers over mammals at high resolution

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.

Regulation of LncRNAs in Melanoma and Their Functional Roles in the Metastatic Process

Long non-coding RNAs (lncRNAs) are key regulators of numerous intracellular processes leading to tumorigenesis. They are frequently deregulated in cancer, functioning as oncogenes or tumor suppressors. As they act through multiple mechanisms, it is not surprising that they may exert dual functions in the same tumor. In melanoma, a highly invasive and metastatic tumor with the propensity to rapidly develop drug resistance, lncRNAs play different roles in: (i) guiding the phenotype switch and leading to metastasis formation; (ii) predicting the response of melanoma patients to immunotherapy; (iii) triggering adaptive responses to therapy and acquisition of drug resistance phenotypes. In this review we summarize the most recent findings on the lncRNAs involved in melanoma growth and spreading to distant sites, focusing on their role as biomarkers for disease diagnosis and patient prognosis, or targets for novel therapeutic approaches.

Capture-C: a modular and flexible approach for high-resolution chromosome conformation capture

Chromosome conformation capture (3C) methods measure the spatial proximity between DNA elements in the cell nucleus. Many methods have been developed to sample 3C material, including the Capture-C family of protocols. Capture-C methods use oligonucleotides to enrich for interactions of interest from sequencing-ready 3C libraries. This approach is modular and has been adapted and optimized to work for sampling of disperse DNA elements (NuTi Capture-C), including from low cell inputs (LI Capture-C), as well as to generate Hi-C like maps for specific regions of interest (Tiled-C) and to interrogate multiway interactions (Tri-C). We present the design, experimental protocol and analysis pipeline for NuTi Capture-C in addition to the variations for generation of LI Capture-C, Tiled-C and Tri-C data. The entire procedure can be performed in 3 weeks and requires standard molecular biology skills and equipment, access to a next-generation sequencing platform, and basic bioinformatic skills. Implemented with other sequencing technologies, these methods can be used to identify regulatory interactions and to compare the structural organization of the genome in different cell types and genetic models.

Accessible Region Conformation Capture (ARC-C) gives high-resolution insights into genome architecture and regulation

Nuclear organization and chromatin interactions are important for genome function, yet determining chromatin connections at high resolution remains a major challenge. To address this, we developed Accessible Region Conformation Capture (ARC-C), which profiles interactions between regulatory elements genome-wide without a capture step. Applied to Caenorhabditis elegans, ARC-C identifies approximately 15,000 significant interactions between regulatory elements at 500-bp resolution. Of 105 TFs or chromatin regulators tested, we find that the binding sites of 60 are enriched for interacting with each other, making them candidates for mediating interactions. These include cohesin and condensin II. Applying ARC-C to a mutant of transcription factor BLMP-1 detected changes in interactions between its targets. ARC-C simultaneously profiles domain-level architecture, and we observe that C. elegans chromatin domains defined by either active or repressive modifications form topologically associating domains (TADs) that interact with A/B (active/inactive) compartment-like structure. Furthermore, we discover that inactive compartment interactions are dependent on H3K9 methylation. ARC-C is a powerful new tool to interrogate genome architecture and regulatory interactions at high resolution.

Inducible lncRNA transgenic mice reveal continual role of HOTAIR in promoting breast cancer metastasis

HOTAIR is a 2.2-kb long noncoding RNA (lncRNA) whose dysregulation has been linked to oncogenesis, defects in pattern formation during early development, and irregularities during the process of epithelial-to-mesenchymal transition (EMT). However, the oncogenic transformation determined by HOTAIR in vivo and its impact on chromatin dynamics are incompletely understood. Here, we generate a transgenic mouse model with doxycycline-inducible expression of human HOTAIR in the context of the MMTV-PyMT breast cancer-prone background to systematically interrogate the cellular mechanisms by which human HOTAIR lncRNA acts to promote breast cancer progression. We show that sustained high levels of HOTAIR over time increased breast metastatic capacity and invasiveness in breast cancer cells, promoting migration and subsequent metastasis to the lung. Subsequent withdrawal of HOTAIR overexpression reverted the metastatic phenotype, indicating oncogenic lncRNA addiction. Furthermore, HOTAIR overexpression altered both the cellular transcriptome and chromatin accessibility landscape of multiple metastasis-associated genes and promoted EMT. These alterations are abrogated within several cell cycles after HOTAIR expression is reverted to basal levels, indicating an erasable lncRNA-associated epigenetic memory. These results suggest that a continual role for HOTAIR in programming a metastatic gene regulatory program. Targeting HOTAIR lncRNA may potentially serve as a therapeutic strategy to ameliorate breast cancer progression.

Hi-C Data Formats

Processing, storing, and visualizing high-resolution Hi-C data required development of efficient data formats. A sparse matrix format saving only nonzero values has become the norm. A "zoomable" matrix style also became popular, storing multiple resolutions in a single file for interactive visualization. This chapter discusses the latest matrix file formats such as .hic and .mcool, and other intermediate formats including SAM/BAM and random-accessible contact lists.

Understanding 3D Genome Organization and Its Effect on Transcriptional Gene Regulation Under Environmental Stress in Plant: A Chromatin Perspective

The genome of a eukaryotic organism is comprised of a supra-molecular complex of chromatin fibers and intricately folded three-dimensional (3D) structures. Chromosomal interactions and topological changes in response to the developmental and/or environmental stimuli affect gene expression. Chromatin architecture plays important roles in DNA replication, gene expression, and genome integrity. Higher-order chromatin organizations like chromosome territories (CTs), A/B compartments, topologically associating domains (TADs), and chromatin loops vary among cells, tissues, and species depending on the developmental stage and/or environmental conditions (4D genomics). Every chromosome occupies a separate territory in the interphase nucleus and forms the top layer of hierarchical structure (CTs) in most of the eukaryotes. While the A and B compartments are associated with active (euchromatic) and inactive (heterochromatic) chromatin, respectively, having well-defined genomic/epigenomic features, TADs are the structural units of chromatin. Chromatin architecture like TADs as well as the local interactions between promoter and regulatory elements correlates with the chromatin activity, which alters during environmental stresses due to relocalization of the architectural proteins. Moreover, chromatin looping brings the gene and regulatory elements in close proximity for interactions. The intricate relationship between nucleotide sequence and chromatin architecture requires a more comprehensive understanding to unravel the genome organization and genetic plasticity. During the last decade, advances in chromatin conformation capture techniques for unravelling 3D genome organizations have improved our understanding of genome biology. However, the recent advances, such as Hi-C and ChIA-PET, have substantially increased the resolution, throughput as well our interest in analysing genome organizations. The present review provides an overview of the historical and contemporary perspectives of chromosome conformation capture technologies, their applications in functional genomics, and the constraints in predicting 3D genome organization. We also discuss the future perspectives of understanding high-order chromatin organizations in deciphering transcriptional regulation of gene expression under environmental stress (4D genomics). These might help design the climate-smart crop to meet the ever-growing demands of food, feed, and fodder.

Interrogating cell type-specific cooperation of transcriptional regulators in 3D chromatin

Context-specific activities of transcription regulators (TRs) in the nucleus modulate spatiotemporal gene expression precisely. Using the largest ChIP-seq data and chromatin loops in the human K562 cell line, we initially interrogated TR cooperation in 3D chromatin via a graphical model and revealed many known and novel TRs manipulating context-specific pathways. To explore TR cooperation across broad tissue/cell types, we systematically leveraged large-scale open chromatin profiles, computational footprinting, and high-resolution chromatin interactions to investigate tissue/cell type-specific TR cooperation. We first delineated a landscape of TR cooperation across 40 human tissue/cell types. Network modularity analyses uncovered the commonality and specificity of TR cooperation in different conditions. We also demonstrated that TR cooperation information can better interpret the disease-causal variants identified by genome-wide association studies and recapitulate cell states during neural development. Our study characterizes shared and unique patterns of TR cooperation associated with the cell type specificity of gene regulation in 3D chromatin.

•

Computational inference of transcriptional regulator (TR) cooperation in 3D chromatin

•

A landscape of 3D TR cooperation across 40 human tissue/cell types

•

TR cooperation can better interpret the disease-causal variants identified by GWAS

•

Cooperation of certain TRs shapes context-specific gene regulation in cell development

Identifying and characterizing lincRNA genomic clusters reveals its cooperative functions in human cancer

Background

Emerging evidence has revealed that some long intergenic non-coding RNAs (lincRNAs) are likely to form clusters on the same chromosome, and lincRNA genomic clusters might play critical roles in the pathophysiological mechanism. However, the comprehensive investigation of lincRNA clustering is rarely studied, particularly the characterization of their functional significance across different cancer types.

Methods

In this study, we firstly constructed a computational method basing a sliding window approach for systematically identifying lincRNA genomic clusters. We then dissected these lincRNA genomic clusters to identify common characteristics in cooperative expression, conservation among divergent species, targeted miRNAs, and CNV frequency. Next, we performed comprehensive analyses in differentially-expressed patterns and overall survival outcomes for patients from The Cancer Genome Atlas (TCGA) and The Genotype-Tissue Expression (GTEx) across multiple cancer types. Finally, we explored the underlying mechanisms of lincRNA genomic clusters by functional enrichment analysis, pathway analysis, and drug-target interaction.

Results

We identified lincRNA genomic clusters according to the algorithm. Clustering lincRNAs tended to be co-expressed, highly conserved, targeted by more miRNAs, and with similar deletion and duplication frequency, suggesting that lincRNA genomic clusters may exert their effects by acting in combination. We further systematically explored conserved and cancer-specific lincRNA genomic clusters, indicating they were involved in some important mechanisms of disease occurrence through diverse approaches. Furthermore, lincRNA genomic clusters can serve as biomarkers with potential clinical significance and involve in specific pathological processes in the development of cancer. Moreover, a lincRNA genomic cluster named Cluster127 in DLK1-DIO3 imprinted locus was discovered, which contained MEG3, MEG8, MEG9, MIR381HG, LINC02285, AL132709.5, and AL132709.1. Further analysis indicated that Cluster127 may have the potential for predicting prognosis in cancer and could play their roles by participating in the regulation of PI3K-AKT signaling pathway.

Conclusions

Clarification of the lincRNA genomic clusters specific roles in human cancers could be beneficial for understanding the molecular pathogenesis of different cancer types.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-021-03179-5.

The macro and micro of chromosome conformation capture

The 3D organization of the genome facilitates gene regulation, replication, and repair, making it a key feature of genomic function and one that remains to be properly understood. Over the past two decades, a variety of chromosome conformation capture (3C) methods have delineated genome folding from megabase-scale compartments and topologically associating domains (TADs) down to kilobase-scale enhancer-promoter interactions. Understanding the functional role of each layer of genome organization is a gateway to understanding cell state, development, and disease. Here, we discuss the evolution of 3C-based technologies for mapping 3D genome organization. We focus on genomics methods and provide a historical account of the development from 3C to Hi-C. We also discuss ChIP-based techniques that focus on 3D genome organization mediated by specific proteins, capture-based methods that focus on particular regions or regulatory elements, 3C-orthogonal methods that do not rely on restriction digestion and proximity ligation, and methods for mapping the DNA-RNA and RNA-RNA interactomes. We consider the biological discoveries that have come from these methods, examine the mechanistic contributions of CTCF, cohesin, and loop extrusion to genomic folding, and detail the 3D genome field's current understanding of nuclear architecture. Finally, we give special consideration to Micro-C as an emerging frontier in chromosome conformation capture and discuss recent Micro-C findings uncovering fine-scale chromatin organization in unprecedented detail. This article is categorized under: Gene Expression and Transcriptional Hierarchies > Regulatory Mechanisms Gene Expression and Transcriptional Hierarchies > Gene Networks and Genomics.

Comprehensive characterization genetic regulation and chromatin landscape of enhancer-associated long non-coding RNAs and their implication in human cancer

Long non-coding RNAs (lncRNAs) that emanate from enhancer regions (defined as enhancer-associated lncRNAs, or elncRNAs) are emerging as critical regulators in disease progression. However, their biological characteristics and clinical relevance have not been fully portrayed. Here, based on the traditional expression quantitative loci (eQTL) and our optimized residual eQTL method, we comprehensively described the genetic effect on elncRNA expression in more than 300 lymphoblastoid cell lines. Meanwhile, a chromatin atlas of elncRNAs relative to the genetic regulation state was depicted. By applying the maximum likelihood estimate method, we successfully identified causal elncRNAs for protein-coding gene expression reprogramming and showed their associated single nucleotide polymorphisms (SNPs) favor binding of transcription factors. Further epigenome analysis revealed two immune-associated elncRNAs AL662844.4 and LINC01215 possess high levels of H3K27ac and H3K4me1 in human cancer. Besides, pan-cancer analysis of 3D genome, transcriptome, and regulatome data showed they potentially regulate tumor-immune cell interaction through affecting MHC class I genes and CD47, respectively. Moreover, our study showed there exist associations between elncRNA and patient survival. Finally, we made a user-friendly web interface available for exploring the regulatory relationship of SNP-elncRNA-protein-coding gene triplets (http://bio-bigdata.hrbmu.edu.cn/elncVarReg). Our study provides critical mechanistic insights for elncRNA function and illustrates their implications in human cancer.

Seeing the forest through the trees: prioritising potentially functional interactions from Hi-C

Eukaryotic genomes are highly organised within the nucleus of a cell, allowing widely dispersed regulatory elements such as enhancers to interact with gene promoters through physical contacts in three-dimensional space. Recent chromosome conformation capture methodologies such as Hi-C have enabled the analysis of interacting regions of the genome providing a valuable insight into the three-dimensional organisation of the chromatin in the nucleus, including chromosome compartmentalisation and gene expression. Complicating the analysis of Hi-C data, however, is the massive amount of identified interactions, many of which do not directly drive gene function, thus hindering the identification of potentially biologically functional 3D interactions. In this review, we collate and examine the downstream analysis of Hi-C data with particular focus on methods that prioritise potentially functional interactions. We classify three groups of approaches: structural-based discovery methods, e.g. A/B compartments and topologically associated domains, detection of statistically significant chromatin interactions, and the use of epigenomic data integration to narrow down useful interaction information. Careful use of these three approaches is crucial to successfully identifying potentially functional interactions within the genome.

Chromatin interaction neural network (ChINN): a machine learning-based method for predicting chromatin interactions from DNA sequences

Chromatin interactions play important roles in regulating gene expression. However, the availability of genome-wide chromatin interaction data is limited. We develop a computational method, chromatin interaction neural network (ChINN), to predict chromatin interactions between open chromatin regions using only DNA sequences. ChINN predicts CTCF- and RNA polymerase II-associated and Hi-C chromatin interactions. ChINN shows good across-sample performances and captures various sequence features for chromatin interaction prediction. We apply ChINN to 6 chronic lymphocytic leukemia (CLL) patient samples and a published cohort of 84 CLL open chromatin samples. Our results demonstrate extensive heterogeneity in chromatin interactions among CLL patient samples.

Understanding 3D genome organization by multidisciplinary methods

Understanding how chromatin is folded in the nucleus is fundamental to understanding its function. Although 3D genome organization has been historically difficult to study owing to a lack of relevant methodologies, major technological breakthroughs in genome-wide mapping of chromatin contacts and advances in imaging technologies in the twenty-first century considerably improved our understanding of chromosome conformation and nuclear architecture. In this Review, we discuss methods of 3D genome organization analysis, including sequencing-based techniques, such as Hi-C and its derivatives, Micro-C, DamID and others; microscopy-based techniques, such as super-resolution imaging coupled with fluorescence in situ hybridization (FISH), multiplex FISH, in situ genome sequencing and live microscopy methods; and computational and modelling approaches. We describe the most commonly used techniques and their contribution to our current knowledge of nuclear architecture and, finally, we provide a perspective on up-and-coming methods that open possibilities for future major discoveries.

Three-dimensional genome organization in the central nervous system, implications for neuropsychological disorders

Chromosomes in eukaryotic cell nuclei are highly compacted and finely organized into hierarchical three-dimensional (3D) configuration. In recent years, scientists have gained deeper understandings of 3D genome structures and revealed novel evidence linking 3D genome organization to various important cell events on the molecular level. Most importantly, alteration of 3D genome architecture has emerged as an intriguing higher order mechanism that connects disease-related genetic variants in multiple heterogenous and polygenic neuropsychological disorders, delivering novel insights into the etiology. In this review, we provide a brief overview of the hierarchical structures of 3D genome and two proposed regulatory models, loop extrusion and phase separation. We then focus on recent Hi-C data in the central nervous system and discuss 3D genome alterations during normal brain development and in mature neurons. Most importantly, we make a comprehensive review on current knowledge and discuss the role of 3D genome in multiple neuropsychological disorders, including schizophrenia, repeat expansion disorders, 22q11 deletion syndrome, and others.

EnHiC: learning fine-resolution Hi-C contact maps using a generative adversarial framework

Motivation

The high-throughput chromosome conformation capture (Hi-C) technique has enabled genome-wide mapping of chromatin interactions. However, high-resolution Hi-C data requires costly, deep sequencing; therefore, it has only been achieved for a limited number of cell types. Machine learning models based on neural networks have been developed as a remedy to this problem.

Results

In this work, we propose a novel method, EnHiC, for predicting high-resolution Hi-C matrices from low-resolution input data based on a generative adversarial network (GAN) framework. Inspired by non-negative matrix factorization, our model fully exploits the unique properties of Hi-C matrices and extracts rank-1 features from multi-scale low-resolution matrices to enhance the resolution. Using three human Hi-C datasets, we demonstrated that EnHiC accurately and reliably enhanced the resolution of Hi-C matrices and outperformed other GAN-based models. Moreover, EnHiC-predicted high-resolution matrices facilitated the accurate detection of topologically associated domains and fine-scale chromatin interactions.

Availability and implementation

EnHiC is publicly available at https://github.com/wmalab/EnHiC.

Supplementary information

Supplementary data are available at Bioinformatics online.

Long noncoding RNAs in cancer metastasis

Metastasis is a major contributor to cancer-associated deaths. It is characterized by a multistep process that occurs through the acquisition of molecular and phenotypic changes enabling cancer cells from a primary tumour to disseminate and colonize at distant organ sites. Over the past decade, the discovery and characterization of long noncoding RNAs (lncRNAs) have revealed the diversity of their regulatory roles, including key contributions throughout the metastatic cascade. Here, we review how lncRNAs promote metastasis by functioning in discrete pro-metastatic steps including the epithelial-mesenchymal transition, invasion and migration and organotrophic colonization, and by influencing the metastatic tumour microenvironment, often by interacting within ribonucleoprotein complexes or directly with other nucleic acid entities. We discuss well-characterized lncRNAs with in vivo phenotypes and highlight mechanistic commonalities such as convergence with the TGFβ-ZEB1/ZEB2 axis or the nuclear factor-κB pathway, in addition to lncRNAs with controversial mechanisms and the influence of methodologies on mechanistic interpretation. Furthermore, some lncRNAs can help identify tumours with increased metastatic risk and spur novel therapeutic strategies, with several lncRNAs having shown potential as novel targets for antisense oligonucleotide therapy in animal models. In addition to well-characterized examples of lncRNAs functioning in metastasis, we discuss controversies and ongoing challenges in lncRNA biology. Finally, we present areas for future study for this rapidly evolving field.

Plant 3D genomics: the exploration and application of chromatin organization

Eukaryotic genomes are highly folded for packing into higher-order chromatin structures in the nucleus. With the emergence of state-of-the-art chromosome conformation capture methods and microscopic imaging techniques, the spatial organization of chromatin and its functional implications have been interrogated. Our knowledge of 3D chromatin organization in plants has improved dramatically in the past few years, building on the early advances in animal systems. Here, we review recent advances in 3D genome mapping approaches, our understanding of the sophisticated organization of spatial structures, and the application of 3D genomic principles in plants. We also discuss directions for future developments in 3D genomics in plants.

Advances in higher-order chromatin architecture: the move towards 4D genome

In eukaryotes, the genome is hierarchically packed inside the nucleus, which facilitates physical contact between cis-regulatory elements (CREs), such as enhancers and promoters. Accumulating evidence highlights the critical role of higher-order chromatin structure in precise regulation of spatiotemporal gene expression under diverse biological contexts including lineage commitment and cell activation by external stimulus. Genomics and imaging-based technologies, such as Hi-C and DNA fluorescence in situ hybridization (FISH), have revealed the key principles of genome folding, while newly developed tools focus on improvement in resolution, throughput and modality at single-cell and population levels, and challenge the knowledge obtained through conventional approaches. In this review, we discuss recent advances in our understanding of principles of higher-order chromosome conformation and technologies to investigate 4D chromatin interactions.

Hi-CO: 3D genome structure analysis with nucleosome resolution

The nucleosome is the basic organizational unit of the genome. The folding structure of nucleosomes is closely related to genome functions, and has been reported to be in dynamic interplay with binding of various nuclear proteins to genomic loci. Here, we describe our high-throughput chromosome conformation capture with nucleosome orientation (Hi-CO) technology to derive 3D nucleosome positions with their orientations at every genomic locus in the nucleus. This technology consists of an experimental procedure for nucleosome proximity analysis and a computational procedure for 3D modeling. The experimental procedure is based on an improved method of high-throughput chromosome conformation capture (Hi-C) analysis. Whereas conventional Hi-C allows spatial proximity analysis among genomic loci with 1-10 kbp resolution, our Hi-CO allows proximity analysis among DNA entry or exit points at every nucleosome locus. This analysis is realized by carrying out ligations among the entry/exit points in every nucleosome in a micrococcal-nuclease-fragmented genome, and by quantifying frequencies of ligation products with next-generation sequencing. Our protocol has enabled this analysis by cleanly excluding unwanted non-ligation products that are abundant owing to the frequent genome fragmentation by micrococcal nuclease. The computational procedure is based on simulated annealing-molecular dynamics, which allows determination of optimized 3D positions and orientations of every nucleosome that satisfies the proximity ligation data sufficiently well. Typically, examination of the Saccharomyces cerevisiae genome with 130 million sequencing reads facilitates analysis of a total of 66,360 nucleosome loci with 6.8 nm resolution. The technique requires 2-3 weeks for sequencing library preparation and 2 weeks for simulation.

covNorm: An R package for coverage based normalization of Hi-C and capture Hi-C data

Hi-C and capture Hi-C have greatly advanced our understanding of the principles of higher-order chromatin structure. In line with the evolution of the Hi-C protocols, there is a demand for an advanced computational method that can be applied to the various forms of Hi-C protocols and effectively remove innate biases. To resolve this issue, we developed an implicit normalization method named “covNorm” and implemented it as an R package. The proposed method can perform a complete procedure of data processing for Hi-C and its variants. Starting from the negative binomial model-based normalization for DNA fragment coverages, removal of genomic distance-dependent background and calling of the significant interactions can be applied sequentially. The performance evaluation of covNorm showed enhanced or similar reproducibility in terms of HiC-spector score, correlation of compartment A/B profiles, and detection of reproducible significant long-range chromatin contacts compared to baseline methods in the benchmark datasets. The developed method is powerful in terms of effective normalization of Hi-C and capture Hi-C data, detection of long-range chromatin contacts, and readily extendibility to the other derivative Hi-C protocols. The covNorm R package is freely available at GitHub: https://github.com/kaistcbfg/covNormRpkg.

The Toolbox for Untangling Chromosome Architecture in Immune Cells

The code of life is not only encrypted in the sequence of DNA but also in the way it is organized into chromosomes. Chromosome architecture is gradually being recognized as an important player in regulating cell activities (e.g., controlling spatiotemporal gene expression). In the past decade, the toolbox for elucidating genome structure has been expanding, providing an opportunity to explore this under charted territory. In this review, we will introduce the recent advancements in approaches for mapping spatial organization of the genome, emphasizing applications of these techniques to immune cells, and trying to bridge chromosome structure with immune cell activities.

Application of Hi-C and other omics data analysis in human cancer and cell differentiation research

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.

A cookbook for DNase Hi-C

Background

The Hi-C technique is widely employed to study the 3-dimensional chromatin architecture and to assemble genomes. The conventional in situ Hi-C protocol employs restriction enzymes to digest chromatin, which results in nonuniform genomic coverage. Using sequence-agnostic restriction enzymes, such as DNAse I, could help to overcome this limitation.

Results

In this study, we compare different DNAse Hi-C protocols and identify the critical steps that significantly affect the efficiency of the protocol. In particular, we show that the SDS quenching strategy strongly affects subsequent chromatin digestion. The presence of biotinylated oligonucleotide adapters may lead to ligase reaction by-products, which can be avoided by rational design of the adapter sequences. Moreover, the use of nucleotide-exchange enzymes for biotin fill-in enables simultaneous labelling and repair of DNA ends, similar to the conventional Hi-C protocol. These improvements simplify the protocol, making it less expensive and time-consuming.

Conclusions

We propose a new robust protocol for the preparation of DNAse Hi-C libraries from cultured human cells and blood samples supplemented with experimental controls and computational tools for the evaluation of library quality.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-021-00389-5.

Molecular and computational approaches to map regulatory elements in 3D chromatin structure

Epigenetic marks do not change the sequence of DNA but affect gene expression in a cell-type specific manner by altering the activities of regulatory elements. Development of new molecular biology assays, sequencing technologies, and computational approaches enables us to profile the human epigenome in three-dimensional structure genome-wide. Here we describe various molecular biology techniques and bioinformatic tools that have been developed to measure the activities of regulatory elements and their chromatin interactions. Moreover, we list currently available three-dimensional epigenomic data sets that are generated in various human cell types and tissues to assist in the design and analysis of research projects.

ASHIC: hierarchical Bayesian modeling of diploid chromatin contacts and structures

The recently developed Hi-C technique has been widely applied to map genome-wide chromatin interactions. However, current methods for analyzing diploid Hi-C data cannot fully distinguish between homologous chromosomes. Consequently, the existing diploid Hi-C analyses are based on sparse and inaccurate allele-specific contact matrices, which might lead to incorrect modeling of diploid genome architecture. Here we present ASHIC, a hierarchical Bayesian framework to model allele-specific chromatin organizations in diploid genomes. We developed two models under the Bayesian framework: the Poisson-multinomial (ASHIC-PM) model and the zero-inflated Poisson-multinomial (ASHIC-ZIPM) model. The proposed ASHIC methods impute allele-specific contact maps from diploid Hi-C data and simultaneously infer allelic 3D structures. Through simulation studies, we demonstrated that ASHIC methods outperformed existing approaches, especially under low coverage and low SNP density conditions. Additionally, in the analyses of diploid Hi-C datasets in mouse and human, our ASHIC-ZIPM method produced fine-resolution diploid chromatin maps and 3D structures and provided insights into the allelic chromatin organizations and functions. To summarize, our work provides a statistically rigorous framework for investigating fine-scale allele-specific chromatin conformations. The ASHIC software is publicly available at https://github.com/wmalab/ASHIC.

Long non-coding RNAs in prostate tumorigenesis and therapy (Review)

Prostate cancer (PCa) is one of the most frequently diagnosed malignancy. Although there have been many advances in PCa diagnosis and therapy, the concrete mechanism remains unknown. Long non-coding RNAs (lncRNAs) are novel biomarkers associated with PCa, and their dysregulated expression is closely associated with risk stratification, diagnosis and carcinogenesis. Accumulating evidence has suggested that lncRNAs play important roles in prostate tumorigenesis through relevant pathways, such as androgen receptor interaction and PI3K/Akt. The present review systematically summarized the potential clinical utility of lncRNAs and provided a novel guide for their function in PCa.

Trans-ancestral dissection of urate- and gout-associated major loci SLC2A9 and ABCG2 reveals primate-specific regulatory effects

Gout is a complex inflammatory arthritis affecting ~20% of people with an elevated serum urate level (hyperuricemia). Gout and hyperuricemia are essentially specific to humans and other higher primates, with varied prevalence across ancestral groups. SLC2A9 and ABCG2 are major loci associated with both urate and gout in multiple ancestral groups. However, fine mapping has been challenging due to extensive linkage disequilibrium underlying the associated regions. We used trans-ancestral fine mapping integrated with primate-specific genomic information to address this challenge. Trans-ancestral meta-analyses of GWAS cohorts of either European (EUR) or East Asian (EAS) ancestry resulted in single-variant resolution mappings for SLC2A9 (rs3775948 for urate and rs4697701 for gout) and ABCG2 (rs2622621 for gout). Tests of colocalization of variants in both urate and gout suggested existence of a shared candidate causal variant for SLC2A9 only in EUR and for ABCG2 only in EAS. The fine-mapped gout variant rs4697701 was within an ancient enhancer, whereas rs2622621 was within a primate-specific transposable element, both supported by functional evidence from the Roadmap Epigenomics project in human primary tissues relevant to urate and gout. Additional primate-specific elements were found near both loci and those adjacent to SLC2A9 overlapped with known statistical epistatic interactions associated with urate as well as multiple super-enhancers identified in urate-relevant tissues. We conclude that by leveraging ancestral differences trans-ancestral fine mapping has identified ancestral and functional variants for SLC2A9 or ABCG2 with primate-specific regulatory effects on urate and gout.

Renal Cell Tumors: Uncovering the Biomarker Potential of ncRNAs

Renal cell tumors (RCT) remain as one of the most common and lethal urological tumors worldwide. Discrimination between (1) benign and malignant disease, (2) indolent and aggressive tumors, and (3) patient responsiveness to a specific therapy is of major clinical importance, allowing for a more efficient patient management. Nonetheless, currently available tools provide limited information and novel strategies are needed. Over the years, a putative role of non-coding RNAs (ncRNAs) as disease biomarkers has gained relevance and is now one of the most prolific fields in biological sciences. Herein, we extensively sought the most significant reports on ncRNAs as potential RCTs' diagnostic, prognostic, predictive, and monitoring biomarkers. We could conclude that ncRNAs, either alone or in combination with currently used clinical and pathological parameters, might represent key elements to improve patient management, potentiating the implementation of precision medicine. Nevertheless, most ncRNA biomarkers require large-scale validation studies, prior to clinical implementation.

Advances in optical mapping for genomic research

Recent advances in optical mapping have allowed the construction of improved genome assemblies with greater contiguity. Optical mapping also enables genome comparison and identification of large-scale structural variations. Association of these large-scale genomic features with biological functions is an important goal in plant and animal breeding and in medical research. Optical mapping has also been used in microbiology and still plays an important role in strain typing and epidemiological studies. Here, we review the development of optical mapping in recent decades to illustrate its importance in genomic research. We detail its applications and algorithms to show its specific advantages. Finally, we discuss the challenges required to facilitate the optimization of optical mapping and improve its future development and application.

Exonuclease combinations reduce noises in 3D genomics technologies

Chromosome conformation-capture technologies are widely used in 3D genomics; however, experimentally, such methods have high-noise limitations and, therefore, require significant bioinformatics efforts to extract reliable distal interactions. Miscellaneous undesired linear DNAs, present during proximity-ligation, represent a main noise source, which needs to be minimized or eliminated. In this study, different exonuclease combinations were tested to remove linear DNA fragments from a circularized DNA preparation. This method efficiently removed linear DNAs, raised the proportion of annulation and increased the valid-pairs ratio from ∼40% to ∼80% for enhanced interaction detection in standard Hi-C. This strategy is applicable for development of various 3D genomics technologies, or optimization of Hi-C sequencing efficiency.

HiCoP, a simple and robust method for detecting interactions of regulatory regions

Background

Chromatin physical interactions provide essential information for understanding the regulation of cis-elements like enhancers, promoters, and insulators in cell development and differentiation. The Hi-C assay is a technique detecting chromatin structures of the whole genome, but not sensitive to interactions of regulatory elements. Several methods, like HiChIP, DNase-C, and OCEAN-C, have been developed for enriching interactions of regulatory regions, but all of them have some shortcomings. New simple, efficient, and robust methods are still in need for detecting interactions of regulatory regions.

Results

We developed a new, simple, and robust assay called CoP (Column Purified chromatin) for profiling of open chromatin regions by directly purifying fragmentized crosslinked chromatin with a DNA purification column. The accessible chromatin regions, including active enhancers, promoters, and insulators, were significantly enriched in CoP chromatin. The CoP-seq assay can efficiently detect open chromatin regions, especially active promoters, with a high signal-to-noise ratio. We integrated the CoP-seq and Hi-C technique (HiCoP) to detect interactions of accessible chromatin regions, which represent active cis-regulatory elements in cells. We observed that the HiCoP captured the peaks in the promoters-associated enhancer regions. HiCoP detected more promoter–enhancer (P–E), promoter–promoter (P–P), and enhancer–enhancer (E–E) interactions within 20 kb–5 Mb than Hi-C. Most of the loops identified by HiCoP are associated with the expressed genes.

Conclusion

CoP assay can efficiently enrich open chromatin regions. When CoP assay was integrated with Hi-C assay, it provides a simple, robust, alternative technique for profiling accessible chromatin regions and chromatin conformation simultaneously.

Architectural proteins for the formation and maintenance of the 3D genome

Eukaryotic genomes are densely packaged into hierarchical three-dimensional (3D) structures that contain information about gene regulation and many other biological processes. With the development of imaging and sequencing-based technologies, 3D genome studies have revealed that the high-order chromatin structure is composed of hierarchical levels, including chromosome territories, A/B compartments, topologically associated domains, and chromatin loops. However, how this chromatin architecture is formed and maintained is not completely clear. In this review, we introduce experimental methods to investigate the 3D genome, review major architectural proteins that regulate 3D chromatin organization in mammalian cells, such as CTCF (CCCTC-binding factor), cohesin, lamins, and transcription factors, and discuss relevant mechanisms such as phase separation.

Technologies to study spatial genome organization: beyond 3C

The way that chromatin is organized in three-dimensional nuclear space is now acknowledged as a factor critical for the major cell processes, like transcription, replication and cell division. Researchers have been armed with new molecular and imaging technologies to study this structure-to-function link of genomes, spearheaded by the introduction of the 'chromosome conformation capture' technology more than a decade ago. However, this technology is not without shortcomings, and novel variants and orthogonal approaches are being developed to overcome these. As a result, the field of nuclear organization is constantly fueled by methods of increasing resolution and/or throughput that strive to eliminate systematic biases and increase precision. In this review, we attempt to highlight the most recent advances in technology that promise to provide novel insights on how chromosomes fold and function.