ChIA-PET2: a versatile and flexible pipeline for ChIA-PET data analysis

A systematic analyses of different bioinformatics pipelines for genomic data and its impact on deep learning models for chromatin loop prediction

Genomic data analysis has witnessed a surge in complexity and volume, primarily driven by the advent of high-throughput technologies. In particular, studying chromatin loops and structures has become pivotal in understanding gene regulation and genome organization. This systematic investigation explores the realm of specialized bioinformatics pipelines designed specifically for the analysis of chromatin loops and structures. Our investigation incorporates two protein (CTCF and Cohesin) factor-specific loop interaction datasets from six distinct pipelines, amassing a comprehensive collection of 36 diverse datasets. Through a meticulous review of existing literature, we offer a holistic perspective on the methodologies, tools and algorithms underpinning the analysis of this multifaceted genomic feature. We illuminate the vast array of approaches deployed, encompassing pivotal aspects such as data preparation pipeline, preprocessing, statistical features and modelling techniques. Beyond this, we rigorously assess the strengths and limitations inherent in these bioinformatics pipelines, shedding light on the interplay between data quality and the performance of deep learning models, ultimately advancing our comprehension of genomic intricacies.

Decoding the genomic landscape of chromatin-associated biomolecular condensates

Biomolecular condensates play a significant role in chromatin activities, primarily by concentrating and compartmentalizing proteins and/or nucleic acids. However, their genomic landscapes and compositions remain largely unexplored due to a lack of dedicated computational tools for systematic identification in vivo. To address this, we develop CondSigDetector, a computational framework designed to detect condensate-like chromatin-associated protein co-occupancy signatures (CondSigs), to predict genomic loci and component proteins of distinct chromatin-associated biomolecular condensates. Applying this framework to mouse embryonic stem cells (mESC) and human K562 cells enable us to depict the high-resolution genomic landscape of chromatin-associated biomolecular condensates, and uncover both known and potentially unknown biomolecular condensates. Multi-omics analysis and experimental validation further verify the condensation properties of CondSigs. Additionally, our investigation sheds light on the impact of chromatin-associated biomolecular condensates on chromatin activities. Collectively, CondSigDetector provides an approach to decode the genomic landscape of chromatin-associated condensates, facilitating a deeper understanding of their biological functions and underlying mechanisms in cells.

CTCF mutation at R567 causes developmental disorders via 3D genome rearrangement and abnormal neurodevelopment

The three-dimensional genome structure organized by CTCF is required for development. Clinically identified mutations in CTCF have been linked to adverse developmental outcomes. Nevertheless, the underlying mechanism remains elusive. In this investigation, we explore the regulatory roles of a clinically relevant R567W point mutation, located within the 11th zinc finger of CTCF, by introducing this mutation into both murine models and human embryonic stem cell-derived cortical organoid models. Mice with homozygous CTCFR567W mutation exhibit growth impediments, resulting in postnatal mortality, and deviations in brain, heart, and lung development at the pathological and single-cell transcriptome levels. This mutation induces premature stem-like cell exhaustion, accelerates the maturation of GABAergic neurons, and disrupts neurodevelopmental and synaptic pathways. Additionally, it specifically hinders CTCF binding to peripheral motifs upstream to the core consensus site, causing alterations in local chromatin structure and gene expression, particularly at the clustered protocadherin locus. Comparative analysis using human cortical organoids mirrors the consequences induced by this mutation. In summary, this study elucidates the influence of the CTCFR567W mutation on human neurodevelopmental disorders, paving the way for potential therapeutic interventions.

A continuum of zinc finger transcription factor retention on native chromatin underlies dynamic genome organization

Transcription factor (TF) residence on chromatin translates into quantitative transcriptional or structural outcomes on genome. Commonly used formaldehyde crosslinking fixes TF-DNA interactions cumulatively and compromises the measured occupancy level. Here we mapped the occupancy level of global or individual zinc finger TFs like CTCF and MAZ, in the form of highly resolved footprints, on native chromatin. By incorporating reinforcing perturbation conditions, we established S-score, a quantitative metric to proxy the continuum of CTCF or MAZ retention across different motifs on native chromatin. The native chromatin-retained CTCF sites harbor sequence features within CTCF motifs better explained by S-score than the metrics obtained from other crosslinking or native assays. CTCF retention on native chromatin correlates with local SUMOylation level, and anti-correlates with transcriptional activity. The S-score successfully delineates the otherwise-masked differential stability of chromatin structures mediated by CTCF, or by MAZ independent of CTCF. Overall, our study established a paradigm continuum of TF retention across binding sites on native chromatin, explaining the dynamic genome organization.

Statistical modeling and significance estimation of multi-way chromatin contacts with HyperloopFinder

Recent advances in chromatin conformation capture technologies, such as SPRITE and Pore-C, have enabled the detection of simultaneous contacts among multiple chromatin loci. This has made it possible to investigate the cooperative transcriptional regulation involving multiple genes and regulatory elements at the resolution of a single molecule. However, these technologies are unavoidably subject to the random polymer looping effect and technical biases, making it challenging to distinguish genuine regulatory relationships directly from random polymer interactions. Here, we present HyperloopFinder, a method for identifying regulatory multi-way chromatin contacts (hyperloops) by jointly modeling the random polymer looping effect and technical biases to estimate the statistical significance of multi-way contacts. The results show that our model can accurately estimate the expected interaction frequency of multi-way contacts based on the distance distribution of pairwise contacts, revealing that most multi-way contacts can be formed by randomly linking the pairwise contacts adjacent to each other. Moreover, we observed the spatial colocalization of the interaction sites of hyperloops from image-based data. Our results also revealed that hyperloops can function as scaffolds for the cooperation among multiple genes and regulatory elements. In summary, our work contributes novel insights into higher-order chromatin structures and functions and has the potential to enhance our understanding of transcriptional regulation and other cellular processes.

TEAD2 initiates ground-state pluripotency by mediating chromatin looping

The transition of mouse embryonic stem cells (ESCs) between serum/LIF and 2i(MEK and GSK3 kinase inhibitor)/LIF culture conditions serves as a valuable model for exploring the mechanisms underlying ground and confused pluripotent states. Regulatory networks comprising core and ancillary pluripotency factors drive the gene expression programs defining stable naïve pluripotency. In our study, we systematically screened factors essential for ESC pluripotency, identifying TEAD2 as an ancillary factor maintaining ground-state pluripotency in 2i/LIF ESCs and facilitating the transition from serum/LIF to 2i/LIF ESCs. TEAD2 exhibits increased binding to chromatin in 2i/LIF ESCs, targeting active chromatin regions to regulate the expression of 2i-specific genes. In addition, TEAD2 facilitates the expression of 2i-specific genes by mediating enhancer-promoter interactions during the serum/LIF to 2i/LIF transition. Notably, deletion of Tead2 results in reduction of a specific set of enhancer-promoter interactions without significantly affecting binding of chromatin architecture proteins, CCCTC-binding factor (CTCF), and Yin Yang 1 (YY1). In summary, our findings highlight a novel prominent role of TEAD2 in orchestrating higher-order chromatin structures of 2i-specific genes to sustain ground-state pluripotency.

MMCT-Loop: a mix model-based pipeline for calling targeted 3D chromatin loops

Protein-specific Chromatin Conformation Capture (3C)-based technologies have become essential for identifying distal genomic interactions with critical roles in gene regulation. The standard techniques include Chromatin Interaction Analysis by Paired-End Tag (ChIA-PET), in situ Hi-C followed by chromatin immunoprecipitation (HiChIP) also known as PLAC-seq. To identify chromatin interactions from these data, a variety of computational methods have emerged. Although these state-of-art methods address many issues with loop calling, only few methods can fit different data types simultaneously, and the accuracy as well as the efficiency these approaches remains limited. Here we have generated a pipeline, MMCT-Loop, which ensures the accurate identification of strong loops as well as dynamic or weak loops through a mixed model. MMCT-Loop outperforms existing methods in accuracy, and the detected loops show higher activation functionality. To highlight the utility of MMCT-Loop, we applied it to conformational data derived from neural stem cell (NSCs) and uncovered several previously unidentified regulatory regions for key master regulators of stem cell identity. MMCT-Loop is an accurate and efficient loop caller for targeted conformation capture data, which supports raw data or pre-processed valid pairs as input, the output interactions are formatted and easily uploaded to a genome browser for visualization.

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.

ChIPr: accurate prediction of cohesin-mediated 3D genome organization from 2D chromatin features

The three-dimensional genome organization influences diverse nuclear processes. Here we present Chromatin Interaction Predictor (ChIPr), a suite of regression models based on deep neural networks, random forest, and gradient boosting to predict cohesin-mediated chromatin interaction strength between any two loci in the genome. The predictions of ChIPr correlate well with ChIA-PET data in four cell lines. The standard ChIPr model requires three experimental inputs: ChIP-Seq signals for RAD21, H3K27ac, and H3K27me3 but works well with just RAD21 signal. Integrative analysis reveals novel insights into the role of CTCF motif, its orientation, and CTCF binding on cohesin-mediated chromatin interactions.

Protocol to capture transcription factor-mediated 3D chromatin interactions using affinity tag-based BL-HiChIP

Pioneer transcription factors (TFs) can directly establish higher-order chromatin interactions to instruct gene transcription. Here, we present a protocol for capturing TF-mediated 3D chromatin interactions using affinity tag-based bridge linker (BL)-Hi-chromatin immunoprecipitation (HiChIP). We describe steps for constructing FLAG-tagged TF, performing BL-HiChIP, and preparing the library. We then detail procedures for sequencing, data analysis, and quality control. This protocol has potential applications in 3D chromatin analysis centered on any specific TF in any type of cells without the need of optimal antibodies. For complete details on the use and execution of this protocol, please refer to Ren et al. (2022).1.

Domains Rearranged Methylase 2 maintains DNA methylation at large DNA hypomethylated shores and long-range chromatin interactions in rice

DNA methylation plays an important role in gene regulation and genomic stability. However, large DNA hypomethylated regions known as DNA methylation valleys (DMVs) or canyons have also been suggested to serve unique regulatory functions, largely unknown in rice (Oryza sativa). Here, we describe the DMVs in rice seedlings, which were highly enriched with developmental and transcription regulatory genes. Further detailed analysis indicated that grand DMVs (gDMVs) might be derived from nuclear integrants of organelle DNA (NORGs). Furthermore, Domains Rearranged Methylase 2 (OsDRM2) maintained DNA methylation at short DMV (sDMV) shores. Epigenetic maps indicated that sDMVs were marked with H3K4me3 and/or H3K27me3, although the loss of DNA methylation had a negligible effect on histone modification within these regions. In addition, we constructed H3K27me3-associated interaction maps for homozygous T-DNA insertion mutant of the gene (osdrm2) and wild type (WT). From a global perspective, most (90%) compartments were stable between osdrm2 and WT plants. At a high resolution, we observed a dramatic loss of long-range chromatin loops in osdrm2, which suffered an extensive loss of non-CG (CHG and CHH, H = A, T, or C) methylation. From another viewpoint, the loss of non-CG methylation at sDMV shores in osdrm2 could disrupt H3K27me3-mediated chromatin interaction networks. Overall, our results demonstrated that DMVs are a key genomic feature in rice and are precisely regulated by epigenetic modifications, including DNA methylation and histone modifications. OsDRM2 maintained DNA methylation at sDMV shores, while OsDRM2 deficiency strongly affected three-dimensional (3D) genome architectures.

Integrated 3D genome, epigenome and transcriptome analyses reveal transcriptional coordination of circadian rhythm in rice

Photoperiods integrate with the circadian clock to coordinate gene expression rhythms and thus ensure plant fitness to the environment. Genome-wide characterization and comparison of rhythmic genes under different light conditions revealed delayed phase under constant darkness (DD) and reduced amplitude under constant light (LL) in rice. Interestingly, ChIP-seq and RNA-seq profiling of rhythmic genes exhibit synchronous circadian oscillation in H3K9ac modifications at their loci and long non-coding RNAs (lncRNAs) expression at proximal loci. To investigate how gene expression rhythm is regulated in rice, we profiled the open chromatin regions and transcription factor (TF) footprints by time-series ATAC-seq. Although open chromatin regions did not show circadian change, a significant number of TFs were identified to rhythmically associate with chromatin and drive gene expression in a time-dependent manner. Further transcriptional regulatory networks mapping uncovered significant correlation between core clock genes and transcription factors involved in light/temperature signaling. In situ Hi-C of ZT8-specific expressed genes displayed highly connected chromatin association at the same time, whereas this ZT8 chromatin connection network dissociates at ZT20, suggesting the circadian control of gene expression by dynamic spatial chromatin conformation. These findings together implicate the existence of a synchronization mechanism between circadian H3K9ac modifications, chromatin association of TF and gene expression, and provides insights into circadian dynamics of spatial chromatin conformation that associate with gene expression rhythms.

Inferring CTCF-binding patterns and anchored loops across human tissues and cell types

CCCTC-binding factor (CTCF) is a transcription regulator with a complex role in gene regulation. The recognition and effects of CTCF on DNA sequences, chromosome barriers, and enhancer blocking are not well understood. Existing computational tools struggle to assess the regulatory potential of CTCF-binding sites and their impact on chromatin loop formation. Here we have developed a deep-learning model, DeepAnchor, to accurately characterize CTCF binding using high-resolution genomic/epigenomic features. This has revealed distinct chromatin and sequence patterns for CTCF-mediated insulation and looping. An optimized implementation of a previous loop model based on DeepAnchor score excels in predicting CTCF-anchored loops. We have established a compendium of CTCF-anchored loops across 52 human tissue/cell types, and this suggests that genomic disruption of these loops could be a general mechanism of disease pathogenesis. These computational models and resources can help investigate how CTCF-mediated cis-regulatory elements shape context-specific gene regulation in cell development and disease progression.

3D organization of regulatory elements for transcriptional regulation in Arabidopsis

BACKGROUND

Although spatial organization of compartments and topologically associating domains at large scale is relatively well studied, the spatial organization of regulatory elements at fine scale is poorly understood in plants.

RESULTS

Here we perform high-resolution chromatin interaction analysis using paired-end tag sequencing approach. We map chromatin interactions tethered with RNA polymerase II and associated with heterochromatic, transcriptionally active, and Polycomb-repressive histone modifications in Arabidopsis. Analysis of the regulatory repertoire shows that distal active cis-regulatory elements are linked to their target genes through long-range chromatin interactions with increased expression of the target genes, while poised cis-regulatory elements are linked to their target genes through long-range chromatin interactions with depressed expression of the target genes. Furthermore, we demonstrate that transcription factor MYC2 is critical for chromatin spatial organization, and propose that MYC2 occupancy and MYC2-mediated chromatin interactions coordinately facilitate transcription within the framework of 3D chromatin architecture. Analysis of functionally related gene-defined chromatin connectivity networks reveals that genes implicated in flowering-time control are functionally compartmentalized into separate subdomains via their spatial activity in the leaf or shoot apical meristem, linking active mark- or Polycomb-repressive mark-associated chromatin conformation to coordinated gene expression.

CONCLUSION

The results reveal that the regulation of gene transcription in Arabidopsis is not only by linear juxtaposition, but also by long-range chromatin interactions. Our study uncovers the fine scale genome organization of Arabidopsis and the potential roles of such organization in orchestrating transcription and development.

Chromatin-associated RNA Dictates the ecDNA Interactome in the Nucleus

Extrachromosomal DNA (ecDNA) promotes cancer by driving copy number heterogeneity and amplifying oncogenes along with functional enhancers. More recent studies suggest two additional mechanisms for further enhancing their oncogenic potential, one via forming ecDNA hubs to augment oncogene expression 1 and the other through acting as portable enhancers to trans-activate target genes 2. However, it has remained entirely elusive about how ecDNA explores the three-dimensional space of the nucleus and whether different ecDNA have distinct interacting mechanisms. Here, by profiling the DNA-DNA and DNA-RNA interactomes in tumor cells harboring different types of ecDNAs in comparison with similarly amplified homogenously staining regions (HSRs) in the chromosome, we show that specific ecDNA interactome is dictated by ecDNA-borne nascent RNA. We demonstrate that the ecDNA co-amplifying PVT1 and MYC utilize nascent noncoding PVT1 transcripts to mediate specific trans-activation of both ecDNA and chromosomal genes. In contrast, the ecDNA amplifying EGFR is weak in this property because of more efficient splicing to remove chromatin-associated nascent RNA. These findings reveal a noncoding RNA-orchestrated program hijacked by cancer cells to enhance the functional impact of amplified oncogenes and associated regulatory elements.

Multiomic analysis of cohesin reveals that ZBTB transcription factors contribute to chromatin interactions

One bottleneck in understanding the principles of 3D chromatin structures is caused by the paucity of known regulators. Cohesin is essential for 3D chromatin organization, and its interacting partners are candidate regulators. Here, we performed proteomic profiling of the cohesin in chromatin and identified transcription factors, RNA-binding proteins and chromatin regulators associated with cohesin. Acute protein degradation followed by time-series genomic binding quantitation and BAT Hi-C analysis were conducted, and the results showed that the transcription factor ZBTB21 contributes to cohesin chromatin binding, 3D chromatin interactions and transcriptional repression. Strikingly, multiomic analyses revealed that the other four ZBTB factors interacted with cohesin, and double degradation of ZBTB21 and ZBTB7B led to a further decrease in cohesin chromatin occupancy. We propose that multiple ZBTB transcription factors orchestrate the chromatin binding of cohesin to regulate chromatin interactions, and we provide a catalog of many additional proteins associated with cohesin that warrant further investigation.

Anchoring super-enhancer-driven oncogenic lncRNAs for anti-tumor therapy in hepatocellular carcinoma

Super-enhancer (SE) plays a vital role in the determination of cell identity and fate. Up-regulated expression of coding genes is frequently associated with SE. However, the transcription dysregulation driven by SE, from the viewpoint of long non-coding RNA (lncRNA), remains unclear. Here, SE-associated lncRNAs in HCC are comprehensively outlined for the first time. This study integrally screens and identifies several novel SE-associated lncRNAs that are highly abundant and sensitive to JQ1. Especially, HSAL3 is identified as an uncharacterized SE-driven oncogenic lncRNA, which is activated by transcription factors HCFC1 and HSF1 via its super-enhancer. HSAL3 interference negatively regulates NOTCH signaling, implying the potential mechanism of its tumor-promoting role. The expression of HSAL3 is increased in HCC samples, and higher HSAL3 expression indicates an inferior overall survival of HCC patients. Furthermore, siHSAL3 loaded nanoparticles exert anti-tumor effect on HCC in vitro and in vivo. In conclusion, this is the first comprehensive survey of SE-associated lncRNAs in HCC. HSAL3 is a novel SE-driven oncogenic lncRNA, and siHSAL3 loaded nanoparticles are therapeutic candidates for HCC. This work sheds lights on the merit of anchoring SE-driven oncogenic lncRNAs for HCC treatment.

Mapping Multi-factor-mediated Chromatin Interactions to Assess Dysregulation of Lung Cancer-related Genes

Studies on the lung cancer genome are indispensable for developing a cure for lung cancer. Whole-genome resequencing, genome-wide association studies, and transcriptome sequencing have greatly improved our understanding of the cancer genome. However, dysregulation of long-range chromatin interactions in lung cancer remains poorly described. To better understand the three-dimensional (3D) genomic interaction features of the lung cancer genome, we used the A549 cell line as a model system and generated high-resolution chromatin interactions associated with RNA polymerase II (RNAPII), CCCTC-binding factor (CTCF), enhancer of zeste homolog 2 (EZH2), and histone 3 lysine 27 trimethylation (H3K27me3) using long-read chromatin interaction analysis by paired-end tag sequencing (ChIA-PET). Analysis showed that EZH2/H3K27me3-mediated interactions further repressed target genes, either through loops or domains, and their distributions along the genome were distinct from and complementary to those associated with RNAPII. Cancer-related genes were highly enriched with chromatin interactions, and chromatin interactions specific to the A549 cell line were associated with oncogenes and tumor suppressor genes, such as additional repressive interactions on FOXO4 and promoter-promoter interactions between NF1 and RNF135. Knockout of an anchor associated with chromatin interactions reversed the dysregulation of cancer-related genes, suggesting that chromatin interactions are essential for proper expression of lung cancer-related genes. These findings demonstrate the 3D landscape and gene regulatory relationships of the lung cancer genome.

Haplotype mapping of H3K27me3-associated chromatin interactions defines topological regulation of gene silencing in rice

Histone modification H3K27me3 is an important chromatin mark that plays vital roles in repressing expression of developmental genes. Here, we construct high-resolution 3D genome maps using long-read chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) and characterize H3K27me3-associated chromatin interactions in an elite rice hybrid, Shanyou 63. We find that many H3K27me3-marked regions may function as silencer-like regulatory elements. The silencer-like elements can come into proximity with distal target genes via forming chromatin loops in 3D space of the nuclei, regulating gene silencing and plant traits. Natural and induced deletion of silencers upregulate expression of distal connected genes. Furthermore, we identify extensive allele-specific chromatin loops. We find that genetic variations alter allelic chromatin topology, thus modulating allelic gene imprinting in rice hybrids. In conclusion, the characterization of silencer-like regulatory elements and haplotype-resolved chromatin interaction maps provide insights into the understanding of molecular mechanisms underlying allelic gene silencing and plant trait controlling.

Prediction of CTCF loop anchor based on machine learning

Introduction: Various activities in biological cells are affected by three-dimensional genome structure. The insulators play an important role in the organization of higher-order structure. CTCF is a representative of mammalian insulators, which can produce barriers to prevent the continuous extrusion of chromatin loop. As a multifunctional protein, CTCF has tens of thousands of binding sites in the genome, but only a portion of them can be used as anchors of chromatin loops. It is still unclear how cells select the anchor in the process of chromatin looping.

Methods: In this paper, a comparative analysis is performed to investigate the sequence preference and binding strength of anchor and non-anchor CTCF binding sites. Furthermore, a machine learning model based on the CTCF binding intensity and DNA sequence is proposed to predict which CTCF sites can form chromatin loop anchors.

Results: The accuracy of the machine learning model that we constructed for predicting the anchor of the chromatin loop mediated by CTCF reached 0.8646. And we find that the formation of loop anchor is mainly influenced by the CTCF binding strength and binding pattern (which can be interpreted as the binding of different zinc fingers).

Discussion: In conclusion, our results suggest that The CTCF core motif and it’s flanking sequence may be responsible for the binding specificity. This work contributes to understanding the mechanism of loop anchor selection and provides a reference for the prediction of CTCF-mediated chromatin loops.

An atlas of gene regulatory networks for memory CD4 + T cells in youth and old age

Aging profoundly affects immune-system function, promoting susceptibility to pathogens, cancers and chronic inflammation. We previously identified a population of IL-10-producing, T follicular helper-like cells (" Tfh10 "), linked to suppressed vaccine responses in aged mice. Here, we integrate single-cell ( sc )RNA-seq, scATAC-seq and genome-scale modeling to characterize Tfh10 - and the full CD4 + memory T cell ( CD4 + TM ) compartment - in young and old mice. We identified 13 CD4 + TM populations, which we validated through cross-comparison to prior scRNA-seq studies. We built gene regulatory networks ( GRNs ) that predict transcription-factor control of gene expression in each T-cell population and how these circuits change with age. Through integration with pan-cell aging atlases, we identified intercellular-signaling networks driving age-dependent changes in CD4 + TM. Our atlas of finely resolved CD4 + TM subsets, GRNs and cell-cell communication networks is a comprehensive resource of predicted regulatory mechanisms operative in memory T cells, presenting new opportunities to improve immune responses in the elderly.

An autoimmune pleiotropic SNP modulates IRF5 alternative promoter usage through ZBTB3-mediated chromatin looping

Genetic sharing is extensively observed for autoimmune diseases, but the causal variants and their underlying molecular mechanisms remain largely unknown. Through systematic investigation of autoimmune disease pleiotropic loci, we found most of these shared genetic effects are transmitted from regulatory code. We used an evidence-based strategy to functionally prioritize causal pleiotropic variants and identify their target genes. A top-ranked pleiotropic variant, rs4728142, yielded many lines of evidence as being causal. Mechanistically, the rs4728142-containing region interacts with the IRF5 alternative promoter in an allele-specific manner and orchestrates its upstream enhancer to regulate IRF5 alternative promoter usage through chromatin looping. A putative structural regulator, ZBTB3, mediates the allele-specific loop to promote IRF5-short transcript expression at the rs4728142 risk allele, resulting in IRF5 overactivation and M1 macrophage polarization. Together, our findings establish a causal mechanism between the regulatory variant and fine-scale molecular phenotype underlying the dysfunction of pleiotropic genes in human autoimmunity.

MINE is a method for detecting spatial density of regulatory chromatin interactions based on a multi-modal network

Chromatin interactions play essential roles in chromatin conformation and gene expression. However, few tools exist to analyze the spatial density of regulatory chromatin interactions (SD-RCI). Here, we present the multi-modal network (MINE) toolkit, including MINE-Loop, MINE-Density, and MINE-Viewer. The MINE-Loop network aims to enhance the detection of RCIs, MINE-Density quantifies the SD--RCI, and MINE-Viewer facilitates 3D visualization of the density of chromatin interactions and participating regulatory factors (e.g., transcription factors). We applied MINE to investigate the relationship between the SD-RCI and chromatin volume change in HeLa cells before and after liquid-liquid phase separation. Changes in SD-RCI before and after treating the HeLa cells with 1,6-hexanediol suggest that changes in chromatin organization was related to the degree of activation or repression of genes. Together, the MINE toolkit enables quantitative studies on different aspects of chromatin conformation and regulatory activity.

Prostate Cancer Transcriptomic Regulation by the Interplay of Germline Risk Alleles, Somatic Mutations, and 3D Genomic Architecture

Prostate cancer is one of the most heritable human cancers. Genome-wide association studies have identified at least 185 prostate cancer germline risk alleles, most noncoding. We used integrative three-dimensional (3D) spatial genomics to identify the chromatin interaction targets of 45 prostate cancer risk alleles, 31 of which were associated with the transcriptional regulation of target genes in 565 localized prostate tumors. To supplement these 31, we verified transcriptional targets for 56 additional risk alleles using linear proximity and linkage disequilibrium analysis in localized prostate tumors. Some individual risk alleles influenced multiple target genes; others specifically influenced only distal genes while leaving proximal ones unaffected. Several risk alleles exhibited widespread germline-somatic interactions in transcriptional regulation, having different effects in tumors with loss of PTEN or RB1 relative to those without. These data clarify functional prostate cancer risk alleles in large linkage blocks and outline a strategy to model multidimensional transcriptional regulation.

SIGNIFICANCE

Many prostate cancer germline risk alleles are enriched in the noncoding regions of the genome and are hypothesized to regulate transcription. We present a 3D genomics framework to unravel risk SNP function and describe the widespread germline-somatic interplay in transcription control. This article is highlighted in the In This Issue feature, p. 2711.

Upregulation of Superenhancer-Driven LncRNA FASRL by USF1 Promotes De Novo Fatty Acid Biosynthesis to Exacerbate Hepatocellular Carcinoma

Superenhancers drive abnormal gene expression in tumors and promote malignancy. However, the relationship between superenhancer-associated long noncoding RNA (lncRNA) and abnormal metabolism is unknown. This study identifies a novel lncRNA, fatty acid synthesis-related lncRNA (FASRL), whose expression is driven by upstream stimulatory factor 1 (USF1) through its superenhancer. FASRL promotes hepatocellular carcinoma (HCC) cell proliferation in vitro and in vivo. Furthermore, FASRL binds to acetyl-CoA carboxylase 1 (ACACA), a fatty acid synthesis rate-limiting enzyme, increasing fatty acid synthesis via the fatty acid metabolism pathway. Moreover, the expression of FASRL, USF1, and ACACA is increased, and their high expression indicates a worse prognosis in HCC patients. In summary, USF1 drives FASRL transcription via a superenhancer. FASRL binding to ACACA increases fatty acid synthesis and lipid accumulation to mechanistically exacerbate HCC. FASRL may serve as a novel prognostic marker and treatment target in HCC.

CTCF acetylation at lysine 20 is required for the early cardiac mesoderm differentiation of embryonic stem cells

The CCCTC-binding factor (CTCF) protein and its modified forms regulate gene expression and genome organization. However, information on CTCF acetylation and its biological function is still lacking. Here, we show that CTCF can be acetylated at lysine 20 (CTCF-K20) by CREB-binding protein (CBP) and deacetylated by histone deacetylase 6 (HDAC6). CTCF-K20 is required for the CTCF interaction with CBP. A CTCF point mutation at lysine 20 had no effect on self-renewal but blocked the mesoderm differentiation of mouse embryonic stem cells (mESCs). The CTCF-K20 mutation reduced CTCF binding to the promoters and enhancers of genes associated with early cardiac mesoderm differentiation, resulting in diminished chromatin accessibility and decreased enhancer-promoter interactions, impairing gene expression. In summary, this study reveals the important roles of CTCF-K20 in regulating CTCF genomic functions and mESC differentiation into mesoderm.

DeepLoop robustly maps chromatin interactions from sparse allele-resolved or single-cell Hi-C data at kilobase resolution

Mapping chromatin loops from noisy Hi-C heatmaps remains a major challenge. Here we present DeepLoop, which performs rigorous bias correction followed by deep-learning-based signal enhancement for robust chromatin interaction mapping from low-depth Hi-C data. DeepLoop enables loop-resolution, single-cell Hi-C analysis. It also achieves a cross-platform convergence between different Hi-C protocols and micrococcal nuclease (micro-C). DeepLoop allowed us to map the genetic and epigenetic determinants of allele-specific chromatin interactions in the human genome. We nominate new loci with allele-specific interactions governed by imprinting or allelic DNA methylation. We also discovered that, in the inactivated X chromosome (X_i), local loops at the DXZ4 'megadomain' boundary escape X-inactivation but the FIRRE 'superloop' locus does not. Importantly, DeepLoop can pinpoint heterozygous single-nucleotide polymorphisms and large structure variants that cause allelic chromatin loops, many of which rewire enhancers with transcription consequences. Taken together, DeepLoop expands the use of Hi-C to provide loop-resolution insights into the genetics of the three-dimensional genome.

CTCF organizes inter-A compartment interactions through RYBP-dependent phase separation

Chromatin is spatially organized into three-dimensional structures at different levels including A/B compartments, topologically associating domains and loops. The canonical CTCF-mediated loop extrusion model can explain the formation of loops. However, the organization mechanisms underlying long-range chromatin interactions such as interactions between A-A compartments are still poorly understood. Here we show that different from the canonical loop extrusion model, RYBP-mediated phase separation of CTCF organizes inter-A compartment interactions. Based on this model, we designed and verified an induced CTCF phase separation system in embryonic stem cells (ESCs), which facilitated inter-A compartment interactions, improved self-renewal of ESCs and inhibited their differentiation toward neural progenitor cells. These findings support a novel and non-canonical role of CTCF in organizing long-range chromatin interactions via phase separation.

Enhancer-promoter interaction maps provide insights into skeletal muscle-related traits in pig genome

Background

Gene expression programs are intimately linked to the interplay of active cis regulatory elements mediated by chromatin contacts and associated RNAs. Genome-wide association studies (GWAS) have identified many variants in these regulatory elements that can contribute to phenotypic diversity. However, the functional interpretation of these variants remains nontrivial due to the lack of chromatin contact information or limited contact resolution. Furthermore, the distribution and role of chromatin-associated RNAs in gene expression and chromatin conformation remain poorly understood. To address this, we first present a comprehensive interaction map of nuclear dynamics of 3D chromatin-chromatin interactions (H3K27ac BL-HiChIP) and RNA-chromatin interactions (GRID-seq) to reveal genomic variants that contribute to complex skeletal muscle traits.

Results

In a genome-wide scan, we provide systematic fine mapping and gene prioritization from GWAS leading signals that underlie phenotypic variability of growth rate, meat quality, and carcass performance. A set of candidate functional variants and 54 target genes previously not detected were identified, with 71% of these candidate functional variants choosing to skip over their nearest gene to regulate the target gene in a long-range manner. The effects of three functional variants regulating KLF6 (related to days to 100 kg), MXRA8 (related to lean meat percentage), and TAF11 (related to loin muscle depth) were observed in two pig populations. Moreover, we find that this multi-omics interaction map consists of functional communities that are enriched in specific biological functions, and GWAS target genes can serve as core genes for exploring peripheral trait-relevant genes.

Conclusions

Our results provide a valuable resource of candidate functional variants for complex skeletal muscle-related traits and establish an integrated approach to complement existing 3D genomics by exploiting RNA-chromatin and chromatin-chromatin interactions for future association studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01322-2.

3D chromatin architecture and transcription regulation in cancer

Chromatin has distinct three-dimensional (3D) architectures important in key biological processes, such as cell cycle, replication, differentiation, and transcription regulation. In turn, aberrant 3D structures play a vital role in developing abnormalities and diseases such as cancer. This review discusses key 3D chromatin structures (topologically associating domain, lamina-associated domain, and enhancer-promoter interactions) and corresponding structural protein elements mediating 3D chromatin interactions [CCCTC-binding factor, polycomb group protein, cohesin, and Brother of the Regulator of Imprinted Sites (BORIS) protein] with a highlight of their associations with cancer. We also summarise the recent development of technologies and bioinformatics approaches to study the 3D chromatin interactions in gene expression regulation, including crosslinking and proximity ligation methods in the bulk cell population (ChIA-PET and HiChIP) or single-molecule resolution (ChIA-drop), and methods other than proximity ligation, such as GAM, SPRITE, and super-resolution microscopy techniques.

Comprehensive enhancer-target gene assignments improve gene set level interpretation of genome-wide regulatory data

BACKGROUND

Revealing the gene targets of distal regulatory elements is challenging yet critical for interpreting regulome data. Experiment-derived enhancer-gene links are restricted to a small set of enhancers and/or cell types, while the accuracy of genome-wide approaches remains elusive due to the lack of a systematic evaluation. We combined multiple spatial and in silico approaches for defining enhancer locations and linking them to their target genes aggregated across >500 cell types, generating 1860 human genome-wide distal enhancer-to-target gene definitions (EnTDefs). To evaluate performance, we used gene set enrichment (GSE) testing on 87 independent ENCODE ChIP-seq datasets of 34 transcription factors (TFs) and assessed concordance of results with known TF Gene Ontology annotations, and other benchmarks.

RESULTS

The top ranked 741 (40%) EnTDefs significantly outperform the common, naïve approach of linking distal regions to the nearest genes, and the top 10 EnTDefs perform well when applied to ChIP-seq data of other cell types. The GSE-based ranking of EnTDefs is highly concordant with ranking based on overlap with curated benchmarks of enhancer-gene interactions. Both our top general EnTDef and cell-type-specific EnTDefs significantly outperform seven independent computational and experiment-based enhancer-gene pair datasets. We show that using our top EnTDefs for GSE with either genome-wide DNA methylation or ATAC-seq data is able to better recapitulate the biological processes changed in gene expression data performed in parallel for the same experiment than our lower-ranked EnTDefs.

CONCLUSIONS

Our findings illustrate the power of our approach to provide genome-wide interpretation regardless of cell type.

YY1 safeguard multidimensional epigenetic landscape associated with extended pluripotency

Although extended pluripotent stem cells (EPSCs) have the potential to form both embryonic and extraembryonic lineages, how their transcriptional regulatory mechanism differs from that of embryonic stem cells (ESCs) remains unclear. Here, we discovered that YY1 binds to specific open chromatin regions in EPSCs. Yy1 depletion in EPSCs leads to a gene expression pattern more similar to that of ESCs than control EPSCs. Moreover, Yy1 depletion triggers a series of epigenetic crosstalk activities, including changes in DNA methylation, histone modifications and high-order chromatin structures. Yy1 depletion in EPSCs disrupts the enhancer-promoter (EP) interactions of EPSC-specific genes, including Dnmt3l. Yy1 loss results in DNA hypomethylation and dramatically reduces the enrichment of H3K4me3 and H3K27ac on the promoters of EPSC-specific genes by upregulating the expression of Kdm5c and Hdac6 through facilitating the formation of CCCTC-binding factor (CTCF)-mediated EP interactions surrounding their loci. Furthermore, single-cell RNA sequencing (scRNA-seq) experiments revealed that YY1 is required for the derivation of extraembryonic endoderm (XEN)-like cells from EPSCs in vitro. Together, this study reveals that YY1 functions as a key regulator of multidimensional epigenetic crosstalk associated with extended pluripotency.

CRISPR-mediated multiplexed live cell imaging of nonrepetitive genomic loci with one guide RNA per locus

Three-dimensional (3D) structures of the genome are dynamic, heterogeneous and functionally important. Live cell imaging has become the leading method for chromatin dynamics tracking. However, existing CRISPR- and TALE-based genomic labeling techniques have been hampered by laborious protocols and are ineffective in labeling non-repetitive sequences. Here, we report a versatile CRISPR/Casilio-based imaging method that allows for a nonrepetitive genomic locus to be labeled using one guide RNA. We construct Casilio dual-color probes to visualize the dynamic interactions of DNA elements in single live cells in the presence or absence of the cohesin subunit RAD21. Using a three-color palette, we track the dynamic 3D locations of multiple reference points along a chromatin loop. Casilio imaging reveals intercellular heterogeneity and interallelic asynchrony in chromatin interaction dynamics, underscoring the importance of studying genome structures in 4D.

CTCF functions as an insulator for somatic genes and a chromatin remodeler for pluripotency genes during reprogramming

CTCF mediates chromatin insulation and long-distance enhancer-promoter (EP) interactions; however, little is known about how these regulatory functions are partitioned among target genes in key biological processes. Here, we show that Ctcf expression is progressively increased during induced pluripotency. In this process, CTCF first functions as a chromatin insulator responsible for direct silencing of the somatic gene expression program and, interestingly, elevated Ctcf expression next ensures chromatin accessibility and contributes to increased EP interactions for a fraction of pluripotency-associated genes. Therefore, CTCF functions in a context-specific manner to modulate the 3D genome to enable cellular reprogramming. We further discover that these context-specific CTCF functions also enlist SMARCA5, an imitation switch (ISWI) chromatin remodeler, together rewiring the epigenome to facilitate cell-fate switch. These findings reveal the dual functions of CTCF in conjunction with a key chromatin remodeler to drive reprogramming toward pluripotency.

Transcription Elongation Machinery Is a Druggable Dependency and Potentiates Immunotherapy in Glioblastoma Stem Cells

Glioblastoma (GBM) is the most lethal primary brain cancer characterized by therapeutic resistance, which is promoted by GBM stem cells (GSC). Here, we interrogated gene expression and whole-genome CRISPR/Cas9 screening in a large panel of patient-derived GSCs, differentiated GBM cells (DGC), and neural stem cells (NSC) to identify master regulators of GSC stemness, revealing an essential transcription state with increased RNA polymerase II-mediated transcription. The YY1 and transcriptional CDK9 complex was essential for GSC survival and maintenance in vitro and in vivo. YY1 interacted with CDK9 to regulate transcription elongation in GSCs. Genetic or pharmacologic targeting of the YY1-CDK9 complex elicited RNA m⁶A modification-dependent interferon responses, reduced regulatory T-cell infiltration, and augmented efficacy of immune checkpoint therapy in GBM. Collectively, these results suggest that YY1-CDK9 transcription elongation complex defines a targetable cell state with active transcription, suppressed interferon responses, and immunotherapy resistance in GBM. SIGNIFICANCE: Effective strategies to rewire immunosuppressive microenvironment and enhance immunotherapy response are still lacking in GBM. YY1-driven transcriptional elongation machinery represents a druggable target to activate interferon response and enhance anti-PD-1 response through regulating the m⁶A modification program, linking epigenetic regulation to immunomodulatory function in GBM.This article is highlighted in the In This Issue feature, p. 275.

Bacon: a comprehensive computational benchmarking framework for evaluating targeted chromatin conformation capture-specific methodologies

Chromatin conformation capture (3C)-based technologies have enabled the accurate detection of topological genomic interactions, and the adoption of ChIP techniques to 3C-based protocols makes it possible to identify long-range interactions. To analyze these large and complex datasets, computational methods are undergoing rapid and expansive evolution. Thus, a thorough evaluation of these analytical pipelines is necessary to identify which commonly used algorithms and processing pipelines need to be improved. Here we present a comprehensive benchmark framework, Bacon, to evaluate the performance of several computational methods. Finally, we provide practical recommendations for users working with HiChIP and/or ChIA-PET analyses.

Diurnal RNAPII-tethered chromatin interactions are associated with rhythmic gene expression in rice

BACKGROUND

The daily cycling of plant physiological processes is speculated to arise from the coordinated rhythms of gene expression. However, the dynamics of diurnal 3D genome architecture and their potential functions underlying the rhythmic gene expression remain unclear.

RESULTS

Here, we reveal the genome-wide rhythmic occupancy of RNA polymerase II (RNAPII), which precedes mRNA accumulation by approximately 2 h. Rhythmic RNAPII binding dynamically correlates with RNAPII-mediated chromatin architecture remodeling at the genomic level of chromatin interactions, spatial clusters, and chromatin connectivity maps, which are associated with the circadian rhythm of gene expression. Rhythmically expressed genes within the same peak phases of expression are preferentially tethered by RNAPII for coordinated transcription. RNAPII-associated chromatin spatial clusters (CSCs) show high plasticity during the circadian cycle, and rhythmically expressed genes in the morning phase and non-rhythmically expressed genes in the evening phase tend to be enriched in RNAPII-associated CSCs to orchestrate expression. Core circadian clock genes are associated with RNAPII-mediated highly connected chromatin connectivity networks in the morning in contrast to the scattered, sporadic spatial chromatin connectivity in the evening; this indicates that they are transcribed within physical proximity to each other during the AM circadian window and are located in discrete "transcriptional factory" foci in the evening, linking chromatin architecture to coordinated transcription outputs.

CONCLUSION

Our findings uncover fundamental diurnal genome folding principles in plants and reveal a distinct higher-order chromosome organization that is crucial for coordinating diurnal dynamics of transcriptional regulation.

Assessing Specific Networks of Chromatin Interactions with HiChIP

The introduction of chromosome conformation capture (3C)-based technologies coupled with next-generation sequencing have significantly advanced our understanding of how the genetic material is organized within the eukaryotic nucleus. Three-dimensional (3D) genomic organization occurs at hierarchical levels, ranging from chromosome territories and subnuclear compartments to smaller self-associated domains and fine-scale chromatin interactions. The latter can be further categorized into different subtypes, such as structural or regulatory, based either on their presumed functionality and/or the factors that mediate their formation. Various enrichment strategies coupled with 3C-based technologies have been developed to prospectively isolate and quantify chromatin interactions around regions occupied by specific proteins or marks of interest. These approaches not only enable high-resolution characterization of the selected chromatin contacts at a cost-effective manner, but also offer important biological insights into their organizational principles and regulatory function. In this chapter, we will focus on the recently developed HiChIP technology with an emphasis on the discovery of putative active enhancers and promoter interactions in cell types of interest. We will describe the specific steps for designing, performing and analyzing successful HiChIP experiments as well as important limitations and considerations.

Quantifying the phase separation property of chromatin-associated proteins under physiological conditions using an anti-1,6-hexanediol index

BACKGROUND

Liquid-liquid phase separation (LLPS) is an important organizing principle for biomolecular condensation and chromosome compartmentalization. However, while many proteins have been reported to undergo LLPS, quantitative and global analysis of chromatin LLPS property remains absent.

RESULTS

Here, by combining chromatin-associated protein pull-down, quantitative proteomics and 1,6-hexanediol (1,6-HD) treatment, we develop Hi-MS and define an anti-1,6-HD index of chromatin-associated proteins (AICAP) to quantify 1,6-HD sensitivity of chromatin-associated proteins under physiological conditions. Compared with known physicochemical properties involved in phase separation, we find that proteins with lower AICAP are associated with higher content of disordered regions, higher hydrophobic residue preference, higher mobility and higher predicted LLPS potential. We also construct BL-Hi-C libraries following 1,6-HD treatment to study the sensitivity of chromatin conformation to 1,6-HD treatment. We find that the active chromatin and high-order structures, as well as the proteins enriched in corresponding regions, are more sensitive to 1,6-HD treatment.

CONCLUSIONS

Our work provides a global quantitative measurement of LLPS properties of chromatin-associated proteins and higher-order chromatin structure. Hi-MS and AICAP data provide an experimental tool and quantitative resources valuable for future studies of biomolecular condensates.

CCIP: Predicting CTCF-mediated chromatin loops with transitivity

Motivation

CTCF-mediated chromatin loops underlie the formation of topological associating domains and serve as the structural basis for transcriptional regulation. However, the formation mechanism of these loops remains unclear, and the genome-wide mapping of these loops is costly and difficult. Motivated by the recent studies on the formation mechanism of CTCF-mediated loops, we studied the possibility of making use of transitivity-related information of interacting CTCF anchors to predict CTCF loops computationally. In this context, transitivity arises when two CTCF anchors interact with the same third anchor by the loop extrusion mechanism and bring themselves close to each other spatially to form an indirect loop.

Results

To determine whether transitivity is informative for predicting CTCF loops and to obtain an accurate and low-cost predicting method, we proposed a two-stage random-forest-based machine learning method, CTCF-mediated Chromatin Interaction Prediction (CCIP), to predict CTCF-mediated chromatin loops. Our two-stage learning approach makes it possible for us to train a prediction model by taking advantage of transitivity-related information as well as functional genome data and genomic data. Experimental studies showed that our method predicts CTCF-mediated loops more accurately than other methods and that transitivity, when used as a properly defined attribute, is informative for predicting CTCF loops. Furthermore, we found that transitivity explains the formation of tandem CTCF loops and facilitates enhancer–promoter interactions. Our work contributes to the understanding of the formation mechanism and function of CTCF-mediated chromatin loops.

Availability and implementation

The source code of CCIP can be accessed at: https://github.com/GaoLabXDU/CCIP.

Supplementary information

Supplementary data are available at Bioinformatics online.

3D genome alterations associated with dysregulated HOXA13 expression in high-risk T-lineage acute lymphoblastic leukemia

3D genome alternations can dysregulate gene expression by rewiring enhancer-promoter interactions and lead to diseases. We report integrated analyses of 3D genome alterations and differential gene expressions in 18 newly diagnosed T-lineage acute lymphoblastic leukemia (T-ALL) patients and 4 healthy controls. 3D genome organizations at the levels of compartment, topologically associated domains and loop could hierarchically classify different subtypes of T-ALL according to T cell differentiation trajectory, similar to gene expressions-based classification. Thirty-four previously unrecognized translocations and 44 translocation-mediated neo-loops are mapped by Hi-C analysis. We find that neo-loops formed in the non-coding region of the genome could potentially regulate ectopic expressions of TLX3, TAL2 and HOXA transcription factors via enhancer hijacking. Importantly, both translocation-mediated neo-loops and NUP98-related fusions are associated with HOXA13 ectopic expressions. Patients with HOXA11-A13 expressions, but not other genes in the HOXA cluster, have immature immunophenotype and poor outcomes. Here, we highlight the potentially important roles of 3D genome alterations in the etiology and prognosis of T-ALL.

The non-coding genome of T-ALL has not been extensively studied. Here, the authors conduct RNA-seq, ATAC-seq and Hi-C seq analyses and find that T-ALL associated neo-loops may regulate key transcription factors including HOXA13; the aberrant expression of which is associated with poor prognosis.

Phase separation of OCT4 controls TAD reorganization to promote cell fate transitions

Topological-associated domains (TADs) are thought to be relatively stable across cell types, although some TAD reorganization has been observed during cellular differentiation. However, little is known about the mechanisms through which TAD reorganization affects cell fate or how master transcription factors affect TAD structures during cell fate transitions. Here, we show extensive TAD reorganization during somatic cell reprogramming, which is correlated with gene transcription and changes in cellular identity. Manipulating TAD reorganization promotes reprogramming, and the dynamics of concentrated chromatin loops in OCT4 phase separated condensates contribute to TAD reorganization. Disrupting OCT4 phase separation attenuates TAD reorganization and reprogramming, which can be rescued by fusing an intrinsically disordered region (IDR) to OCT4. We developed an approach termed TAD reorganization-based multiomics analysis (TADMAN), which identified reprogramming regulators. Together, these findings elucidate a role and mechanism of TAD reorganization, regulated by OCT4 phase separation, in cellular reprogramming.

Loop competition and extrusion model predicts CTCF interaction specificity

Three-dimensional chromatin looping interactions play an important role in constraining enhancer–promoter interactions and mediating transcriptional gene regulation. CTCF is thought to play a critical role in the formation of these loops, but the specificity of which CTCF binding events form loops and which do not is difficult to predict. Loops often have convergent CTCF binding site motif orientation, but this constraint alone is only weakly predictive of genome-wide interaction data. Here we present an easily interpretable and simple mathematical model of CTCF mediated loop formation which is consistent with Cohesin extrusion and can predict ChIA-PET CTCF looping interaction measurements with high accuracy. Competition between overlapping loops is a critical determinant of loop specificity. We show that this model is consistent with observed chromatin interaction frequency changes induced by CTCF binding site deletion, inversion, and mutation, and is also consistent with observed constraints on validated enhancer–promoter interactions.

Boundaries of topologically associated domains in genomes are marked by CTCF and cohesin binding. Here the authors predict CTCF interaction specificity by building a simple mathematical model with features including loop competition and extrusion.

ERG orchestrates chromatin interactions to drive prostate cell fate reprogramming

Although cancer is commonly perceived as a disease of dedifferentiation, the hallmark of early-stage prostate cancer is paradoxically the loss of more plastic basal cells and the abnormal proliferation of more differentiated secretory luminal cells. However, the mechanism of prostate cancer proluminal differentiation is largely unknown. Through integrating analysis of the transcription factors (TFs) from 806 human prostate cancers, we found that ERG was highly correlated with prostate cancer luminal subtyping. ERG overexpression in luminal epithelial cells inhibited those cells' normal plasticity to transdifferentiate into a basal lineage, and ERG superseded PTEN loss, which favored basal differentiation. ERG KO disrupted prostate cell luminal differentiation, whereas AR KO had no such effects. Trp63 is a known master regulator of the prostate basal lineage. Through analysis of 3D chromatin architecture, we found that ERG bound and inhibited the enhancer activity and chromatin looping of a Trp63 distal enhancer, thereby silencing its gene expression. Specific deletion of the distal ERG binding site resulted in the loss of ERG-mediated inhibition of basal differentiation. Thus, ERG, in its fundamental role in lineage differentiation in prostate cancer initiation, orchestrated chromatin interactions and regulated prostate cell lineage toward a proluminal program.

3 '-5 ' crosstalk contributes to transcriptional bursting

BACKGROUND

Transcription in mammalian cells is a complex stochastic process involving shuttling of polymerase between genes and phase-separated liquid condensates. It occurs in bursts, which results in vastly different numbers of an mRNA species in isogenic cell populations. Several factors contributing to transcriptional bursting have been identified, usually classified as intrinsic, in other words local to single genes, or extrinsic, relating to the macroscopic state of the cell. However, some possible contributors have not been explored yet. Here, we focus on processes at the 3 ' and 5 ' ends of a gene that enable reinitiation of transcription upon termination.

RESULTS

Using Bayesian methodology, we measure the transcriptional bursting in inducible transgenes, showing that perturbation of polymerase shuttling typically reduces burst size, increases burst frequency, and thus limits transcriptional noise. Analysis based on paired-end tag sequencing (PolII ChIA-PET) suggests that this effect is genome wide. The observed noise patterns are also reproduced by a generative model that captures major characteristics of the polymerase flux between the ends of a gene and a phase-separated compartment.

CONCLUSIONS

Interactions between the 3 ' and 5 ' ends of a gene, which facilitate polymerase recycling, are major contributors to transcriptional noise.

The landscape of RNA polymerase II-associated chromatin interactions in prostate cancer

Transcriptional dysregulation is a hallmark of prostate cancer (PCa). We mapped the RNA polymerase II-associated (RNA Pol II-associated) chromatin interactions in normal prostate cells and PCa cells. We discovered thousands of enhancer-promoter, enhancer-enhancer, as well as promoter-promoter chromatin interactions. These transcriptional hubs operate within the framework set by structural proteins - CTCF and cohesins - and are regulated by the cooperative action of master transcription factors, such as the androgen receptor (AR) and FOXA1. By combining analyses from metastatic castration-resistant PCa (mCRPC) specimens, we show that AR locus amplification contributes to the transcriptional upregulation of the AR gene by increasing the total number of chromatin interaction modules comprising the AR gene and its distal enhancer. We deconvoluted the transcription control modules of several PCa genes, notably the biomarker KLK3, lineage-restricted genes (KRT8, KRT18, HOXB13, FOXA1, ZBTB16), the drug target EZH2, and the oncogene MYC. By integrating clinical PCa data, we defined a germline-somatic interplay between the PCa risk allele rs684232 and the somatically acquired TMPRSS2-ERG gene fusion in the transcriptional regulation of multiple target genes - VPS53, FAM57A, and GEMIN4. Our studies implicate changes in genome organization as a critical determinant of aberrant transcriptional regulation in PCa.

Integrative Epigenomic Analysis of Transcriptional Regulation of Human CircRNAs

Circular RNAs (circRNAs) are evolutionarily conserved and abundant non-coding RNAs whose functions and regulatory mechanisms remain largely unknown. Here, we identify and characterize an epigenomically distinct group of circRNAs (TAH-circRNAs), which are transcribed to a higher level than their host genes. By integrative analysis of cistromic and transcriptomic data, we find that compared with other circRNAs, TAH-circRNAs are expressed more abundantly and have more transcription factors (TFs) binding sites and lower DNA methylation levels. Concordantly, TAH-circRNAs are enriched in open and active chromatin regions. Importantly, ChIA-PET results showed that 23–52% of transcription start sites (TSSs) of TAH-circRNAs have direct interactions with cis-regulatory regions, strongly suggesting their independent transcriptional regulation from host genes. In addition, we characterize molecular features of super-enhancer-driven circRNAs in cancer biology. Together, this study comprehensively analyzes epigenomic characteristics of circRNAs and identifies a distinct group of TAH-circRNAs that are independently transcribed via enhancers and super-enhancers by TFs. These findings substantially advance our understanding of the regulatory mechanism of circRNAs and may have important implications for future investigations of this class of non-coding RNAs.

ChIAMM: A Mixture Model for Statistical Analysis of Long-Range Chromatin Interactions From ChIA-PET Experiments

Chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) is an important experimental method for detecting specific protein-mediated chromatin loops genome-wide at high resolution. Here, we proposed a new statistical approach with a mixture model, chromatin interaction analysis using mixture model (ChIAMM), to detect significant chromatin interactions from ChIA-PET data. The statistical model is cast into a Bayesian framework to consider more systematic biases: the genomic distance, local enrichment, mappability, and GC content. Using different ChIA-PET datasets, we evaluated the performance of ChIAMM and compared it with the existing methods, including ChIA-PET Tool, ChiaSig, Mango, ChIA-PET2, and ChIAPoP. The result showed that the new approach performed better than most top existing methods in detecting significant chromatin interactions in ChIA-PET experiments.

3D landscape of Hepatitis B virus interactions with human chromatins

Hepatitis B viral (HBV) DNAs, including covalently closed circular DNA (cccDNA) and integrated HBV DNA forms, are considered to be primary contributors to the development and progression of HBV-associated liver diseases. However, it remains largely unclear how HBV DNAs communicate with human chromatin. Here we employed a highly sensitive technology, 3C-high-throughput genome-wide translocation sequencing (3C-HTGTS), to globally identify HBV DNA-host DNA contacts in cellular models of HBV infection. HBV DNA does not randomly position in host genome but instead preferentially establishes contacts with the host DNA at active chromatin regions. HBV DNA-host DNA contacts are significantly enriched at H3K4me1-marked regions modified by KMT2C/D; this histone modification is also observed in the HBV cccDNA mini-chromosome and strongly influences HBV transcription. On the other hand, chromatin loop formed by integrated HBV DNA with host genomic DNA was found in transcriptionally active regions. Furthermore, HBV infection influences host gene expression accompanied with HBV DNA-host DNA contacts. Our study provides a 3D landscape of spatial organization of cccDNA and integrated HBV DNA within the human genome, which lays the foundation for a better understanding of the mechanisms how HBV involves in liver disease development and progression.