Insulated Neighborhoods: Structural and Functional Units of Mammalian Gene Control

Dynamic RNA polymerase II occupancy drives differentiation of the intestine under the direction of HNF4

Terminal differentiation requires massive restructuring of the transcriptome. During intestinal differentiation, the expression patterns of nearly 4,000 genes are altered as cells transition from progenitor cells in crypts to differentiated cells in villi. We identify dynamic occupancy of RNA polymerase II (Pol II) to gene promoters as the primary driver of transcriptomic shifts during intestinal differentiation in vivo. Changes in enhancer-promoter looping interactions accompany dynamic Pol II occupancy and are dependent upon HNF4, a pro-differentiation transcription factor. Using genetic loss-of-function, chromatin immunoprecipitation sequencing (ChIP-seq), and immunoprecipitation (IP) mass spectrometry, we demonstrate that HNF4 collaborates with chromatin remodelers and loop-stabilizing proteins and facilitates Pol II occupancy at hundreds of genes pivotal to differentiation. We also explore alternate mechanisms that drive differentiation gene expression and find that pause-release of Pol II and post-transcriptional mRNA stability regulate smaller subsets of differentially expressed genes. These studies provide insights into the mechanisms of differentiation in renewing adult tissue.

Effects of super-enhancers in cancer metastasis: mechanisms and therapeutic targets

Metastasis remains the principal cause of cancer-related lethality despite advancements in cancer treatment. Dysfunctional epigenetic alterations are crucial in the metastatic cascade. Among these, super-enhancers (SEs), emerging as new epigenetic regulators, consist of large clusters of regulatory elements that drive the high-level expression of genes essential for the oncogenic process, upon which cancer cells develop a profound dependency. These SE-driven oncogenes play an important role in regulating various facets of metastasis, including the promotion of tumor proliferation in primary and distal metastatic organs, facilitating cellular migration and invasion into the vasculature, triggering epithelial-mesenchymal transition, enhancing cancer stem cell-like properties, circumventing immune detection, and adapting to the heterogeneity of metastatic niches. This heavy reliance on SE-mediated transcription delineates a vulnerable target for therapeutic intervention in cancer cells. In this article, we review current insights into the characteristics, identification methodologies, formation, and activation mechanisms of SEs. We also elaborate the oncogenic roles and regulatory functions of SEs in the context of cancer metastasis. Ultimately, we discuss the potential of SEs as novel therapeutic targets and their implications in clinical oncology, offering insights into future directions for innovative cancer treatment strategies.

SATB1, senescence and senescence-related diseases

Aging leads to an accumulation of cellular mutations and damage, increasing the risk of senescence, apoptosis, and malignant transformation. Cellular senescence, which is pivotal in aging, acts as both a guard against cellular transformation and as a check against cancer progression. It is marked by stable cell cycle arrest, widespread macromolecular changes, a pro-inflammatory profile, and altered gene expression. However, it remains to be determined whether these differing subsets of senescent cells result from unique intrinsic programs or are influenced by their environmental contexts. Multiple transcription regulators and chromatin modifiers contribute to these alterations. Special AT-rich sequence-binding protein 1 (SATB1) stands out as a crucial regulator in this process, orchestrating gene expression by structuring chromatin into loop domains and anchoring DNA elements. This review provides an overview of cellular senescence and delves into the role of SATB1 in senescence-related diseases. It highlights SATB1's potential in developing antiaging and anticancer strategies, potentially contributing to improved quality of life and addressing aging-related diseases.

Tet-mediated DNA methylation dynamics affect chromosome organization

DNA Methylation is a significant epigenetic modification that can modulate chromosome states, but its role in orchestrating chromosome organization has not been well elucidated. Here we systematically assessed the effects of DNA Methylation on chromosome organization with a multi-omics strategy to capture DNA Methylation and high-order chromosome interaction simultaneously on mouse embryonic stem cells with DNA methylation dioxygenase Tet triple knock-out (Tet-TKO). Globally, upon Tet-TKO, we observed weakened compartmentalization, corresponding to decreased methylation differences between CpG island (CGI) rich and poor domains. Tet-TKO could also induce hypermethylation for the CTCF binding peaks in TAD boundaries and chromatin loop anchors. Accordingly, CTCF peak generally weakened upon Tet-TKO, which results in weakened TAD structure and depletion of long-range chromatin loops. Genes that lost enhancer-promoter looping upon Tet-TKO showed DNA hypermethylation in their gene bodies, which may compensate for the disruption of gene expression. We also observed distinct effects of Tet1 and Tet2 on chromatin organization and increased DNA methylation correlation on spatially interacted fragments upon Tet inactivation. Our work showed the broad effects of Tet inactivation and DNA methylation dynamics on chromosome organization.

Structural basis for the preservation of a subset of topologically associating domains in interphase chromosomes upon cohesin depletion

Compartment formation in interphase chromosomes is a result of spatial segregation between euchromatin and heterochromatin on a few megabase pairs (Mbp) scale. On the sub-Mbp scales, topologically associating domains (TADs) appear as interacting domains along the diagonal in the ensemble averaged Hi-C contact map. Hi-C experiments showed that most of the TADs vanish upon deleting cohesin, while the compartment structure is maintained, and perhaps even enhanced. However, closer inspection of the data reveals that a non-negligible fraction of TADs is preserved (P-TADs) after cohesin loss. Imaging experiments show that, at the single-cell level, TAD-like structures are present even without cohesin. To provide a structural basis for these findings, we first used polymer simulations to show that certain TADs with epigenetic switches across their boundaries survive after depletion of loops. More importantly, the three-dimensional structures show that many of the P-TADs have sharp physical boundaries. Informed by the simulations, we analyzed the Hi-C maps (with and without cohesin) in mouse liver and human colorectal carcinoma cell lines, which affirmed that epigenetic switches and physical boundaries (calculated using the predicted 3D structures using the data-driven HIPPS method that uses Hi-C as the input) explain the origin of the P-TADs. Single-cell structures display TAD-like features in the absence of cohesin that are remarkably similar to the findings in imaging experiments. Some P-TADs, with physical boundaries, are relevant to the retention of enhancer-promoter/promoter-promoter interactions. Overall, our study shows that preservation of a subset of TADs upon removing cohesin is a robust phenomenon that is valid across multiple cell lines.

Cohesin mediated loop extrusion from active enhancers form chromatin jets in C. elegans

In mammals, cohesin and CTCF organize the 3D genome into topologically associated domains (TADs) to regulate communication between cis-regulatory elements. However, many organisms, including S. cerevisiae , C. elegans , and A. thaliana lack CTCF. Here, we use C. elegans as a model to investigate the function of cohesin in 3D genome organization in an animal without CTCF. We use auxin-inducible degradation to acutely deplete SMC-3 or its negative regulator WAPL-1 from somatic cells. Using Hi-C data, we identify a cohesin-dependent 3D genome organization feature called chromatin jets (aka fountains). These are population average reflections of DNA loops that are ∼20-40 kb in scale and often cover a few transcribed genes. The jets emerge from NIPBL occupied segments, and the trajectory of the jets coincides with cohesin binding. Cohesin translocation from jet origins depends on a fully intact complex and is extended upon WAPL-1 depletion. Hi-C results support the idea that cohesin is preferentially loaded at NIPBL occupied sites and loop extrudes in an effectively two-sided manner. The location of putative loading sites coincide with active enhancers and the strength of chromatin jet pattern correlates with transcription. Hi-C analyses upon WAPL-1 depletion reveal unequal loop extrusion processivity on each side and stalling due to cohesin molecules colliding. Compared to mammalian systems, average processivity of C. elegans cohesin is ∼10-fold shorter and NIPBL binding does not depend on cohesin. We conclude that the processivity of cohesin scales with genome size; and regardless of CTCF presence, preferential loading of cohesin at enhancers is a conserved mechanism of genome organization that regulates the interaction of gene regulatory elements in 3D.

TP63-TRIM29 axis regulates enhancer methylation and chromosomal instability in prostate cancer

BACKGROUND

Prostate adenocarcinoma (PRAD) is the second leading cause of cancer-related deaths in men. High variability in DNA methylation and a high rate of large genomic rearrangements are often observed in PRAD.

RESULTS

To investigate the reasons for such high variance, we integrated DNA methylation, RNA-seq, and copy number alterations datasets from The Cancer Genome Atlas (TCGA), focusing on PRAD, and employed weighted gene co-expression network analysis (WGCNA). Our results show that only single cluster of co-expressed genes is associated with genomic and epigenomic instability. Within this cluster, TP63 and TRIM29 are key transcription regulators and are downregulated in PRAD. We discovered that TP63 regulates the level of enhancer methylation in prostate basal epithelial cells. TRIM29 forms a complex with TP63 and together regulates the expression of genes specific to the prostate basal epithelium. In addition, TRIM29 binds DNA repair proteins and prevents the formation of the TMPRSS2:ERG gene fusion typically observed in PRAD.

CONCLUSION

Our study demonstrates that TRIM29 and TP63 are important regulators in maintaining the identity of the basal epithelium under physiological conditions. Furthermore, we uncover the role of TRIM29 in PRAD development.

3D organization of enhancers in MuSCs

Skeletal muscle stem cells (MuSCs), also known as satellite cells, are essential for muscle growth and injury induced regeneration. In healthy adult muscle, MuSCs remain in a quiescent state located in a specialized niche beneath the basal lamina. Upon injury, these dormant MuSCs can quickly activate to re-enter the cell cycle and differentiate into new myofibers, while a subset undergoes self-renewal and returns to quiescence to restore the stem cell pool. The myogenic lineage progression is intricately controlled by complex intrinsic and extrinsic cues and coupled with dynamic transcriptional programs. In transcriptional regulation, enhancers are key regulatory elements controlling spatiotemporal gene expression through physical contacting promoters of target genes. The three-dimensional (3D) chromatin architecture is known to orchestrate the establishment of proper enhancer-promoter interactions throughout development and aging. However, studies dissecting the 3D organization of enhancers in MuSCs are just emerging. Here, we provide an overview of the general properties of enhancers and newly developed methods for assessing their activity. In particular, we summarize recent discoveries regarding the 3D rewiring of enhancers during MuSC specification, lineage progression as well as aging.

HDAC inhibitors as pharmacological treatment for Duchenne muscular dystrophy: a discovery journey from bench to patients

Earlier evidence that targeting the balance between histone acetyltransferases (HATs) and deacetylases (HDACs), through exposure to HDAC inhibitors (HDACis), could enhance skeletal myogenesis, prompted interest in using HDACis to promote muscle regeneration. Further identification of constitutive HDAC activation in dystrophin-deficient muscles, caused by dysregulated nitric oxide (NO) signaling, provided the rationale for HDACi-based therapeutic interventions for Duchenne muscular dystrophy (DMD). In this review, we describe the molecular, preclinical, and clinical evidence supporting the efficacy of HDACis in countering disease progression by targeting pathogenic networks of gene expression in multiple muscle-resident cell types of patients with DMD. Given that givinostat is paving the way for HDACi-based interventions in DMD, next-generation HDACis with optimized therapeutic profiles and efficacy could be also explored for synergistic combinations with other therapeutic strategies.

DeepCBS: shedding light on the impact of mutations occurring at CTCF binding sites

CTCF-mediated chromatin loops create insulated neighborhoods that constrain promoter-enhancer interactions, serving as a unit of gene regulation. Disruption of the CTCF binding sites (CBS) will lead to the destruction of insulated neighborhoods, which in turn can cause dysregulation of the contained genes. In a recent study, it is found that CTCF/cohesin binding sites are a major mutational hotspot in the cancer genome. Mutations can affect CTCF binding, causing the disruption of insulated neighborhoods. And our analysis reveals a significant enrichment of well-known proto-oncogenes in insulated neighborhoods with mutations specifically occurring in anchor regions. It can be assumed that some mutations disrupt CTCF binding, leading to the disruption of insulated neighborhoods and subsequent activation of proto-oncogenes within these insulated neighborhoods. To explore the consequences of such mutations, we develop DeepCBS, a computational tool capable of analyzing mutations at CTCF binding sites, predicting their influence on insulated neighborhoods, and investigating the potential activation of proto-oncogenes. Futhermore, DeepCBS is applied to somatic mutation data of liver cancer. As a result, 87 mutations that disrupt CTCF binding sites are identified, which leads to the identification of 237 disrupted insulated neighborhoods containing a total of 135 genes. Integrative analysis of gene expression differences in liver cancer further highlights three genes: ARHGEF39, UBE2C and DQX1. Among them, ARHGEF39 and UBE2C have been reported in the literature as potential oncogenes involved in the development of liver cancer. The results indicate that DQX1 may be a potential oncogene in liver cancer and may contribute to tumor immune escape. In conclusion, DeepCBS is a promising method to analyze impacts of mutations occurring at CTCF binding sites on the insulator function of CTCF, with potential extensions to shed light on the effects of mutations on other functions of CTCF.

Topological reorganization and functional alteration of distinct genomic components in gallbladder cancer

Altered three-dimensional architecture of chromatin influences various genomic regulators and subsequent gene expression in human cancer. However, knowledge of the topological rearrangement of genomic hierarchical layers in cancer is largely limited. Here, by taking advantage of in situ Hi-C, RNA-sequencing, and chromatin immunoprecipitation sequencing (ChIP-seq), we investigated structural reorganization and functional changes in chromosomal compartments, topologically associated domains (TADs), and CCCTC binding factor (CTCF)-mediated loops in gallbladder cancer (GBC) tissues and cell lines. We observed that the chromosomal compartment A/B switch was correlated with CTCF binding levels and gene expression changes. Increased inter-TAD interactions with weaker TAD boundaries were identified in cancer cell lines relative to normal controls. Furthermore, the chromatin short loops and cancer unique loops associated with chromatin remodeling and epithelial-mesenchymal transition activation were enriched in cancer compared with their control counterparts. Cancer-specific enhancer-promoter loops, which contain multiple transcription factor binding motifs, acted as a central element to regulate aberrant gene expression. Depletion of individual enhancers in each loop anchor that connects with promoters led to the inhibition of their corresponding gene expressions. Collectively, our data offer the landscape of hierarchical layers of cancer genome and functional alterations that contribute to the development of GBC.

YY1 Binding to Regulatory Elements That Lack Enhancer Activity Promotes Locus Folding and Gene Activation

Enhancers activate their cognate promoters over huge distances but how enhancer/promoter interactions become established is not completely understood. There is strong evidence that cohesin-mediated loop extrusion is involved but this does not appear to be a universal mechanism. Here, we identify an element within the mouse immunoglobulin lambda (Igλ) light chain locus, HSCλ1, that has characteristics of active regulatory elements but lacks intrinsic enhancer or promoter activity. Remarkably, knock-out of the YY1 binding site from HSCλ1 reduces Igλ transcription significantly and disrupts enhancer/promoter interactions, even though these elements are >10 kb from HSCλ1. Genome-wide analyses of mouse embryonic stem cells identified 2671 similar YY1-bound, putative genome organizing elements that lie within CTCF/cohesin loop boundaries but that lack intrinsic enhancer activity. We suggest that such elements play a fundamental role in locus folding and in facilitating enhancer/promoter interactions.

Dynamic RNA Polymerase II Recruitment Drives Differentiation of the Intestine under the direction of HNF4

Terminal differentiation requires a massive restructuring of the transcriptome. During intestinal differentiation, the expression patterns of nearly 4000 genes are altered as cells transition from progenitor cells in crypts to differentiated cells in villi. We identified dynamic recruitment of RNA Polymerase II (Pol II) to gene promoters as the primary driver of transcriptomic shifts during intestinal differentiation in vivo. Changes in enhancer-promoter looping interactions accompany dynamic Pol II recruitment and are dependent upon HNF4, a pro-differentiation transcription factor. Using genetic loss-of- function, ChIP-seq and IP mass spectrometry, we demonstrate that HNF4 collaborates with chromatin remodelers and loop-stabilizing proteins and facilitates Pol II recruitment at hundreds of genes pivotal to differentiation. We also explore alternate mechanisms which drive differentiation gene expression and find pause-release of Pol II and post- transcriptional mRNA stability regulate smaller subsets of differentially expressed genes. These studies provide insights into the mechanisms of differentiation in a renewing adult tissue.

Temporal analyses reveal a pivotal role for sense and antisense enhancer RNAs in coordinate immunoglobulin lambda locus activation

Transcription enhancers are essential activators of V(D)J recombination that orchestrate non-coding transcription through complementary, unrearranged gene segments. How transcription is coordinately increased at spatially distinct promoters, however, remains poorly understood. Using the murine immunoglobulin lambda (Igλ) locus as model, we find that three enhancer-like elements in the 3' Igλ domain, Eλ3-1, HSCλ1 and HSE-1, show strikingly similar transcription factor binding dynamics and close spatial proximity, suggesting that they form an active enhancer hub. Temporal analyses show coordinate recruitment of complementary V and J gene segments to this hub, with comparable transcription factor binding dynamics to that at enhancers. We find further that E2A, p300, Mediator and Integrator bind to enhancers as early events, whereas YY1 recruitment and eRNA synthesis occur later, corresponding to transcription activation. Remarkably, the interplay between sense and antisense enhancer RNA is central to both active enhancer hub formation and coordinate Igλ transcription: Antisense Eλ3-1 eRNA represses Igλ activation whereas temporal analyses demonstrate that accumulating levels of sense eRNA boost YY1 recruitment to stabilise enhancer hub/promoter interactions and lead to coordinate transcription activation. These studies therefore demonstrate for the first time a critical role for threshold levels of sense versus antisense eRNA in locus activation.

MYC function and regulation in physiological perspective

MYC, a key member of the Myc-proto-oncogene family, is a universal transcription amplifier that regulates almost every physiological process in a cell including cell cycle, proliferation, metabolism, differentiation, and apoptosis. MYC interacts with several cofactors, chromatin modifiers, and regulators to direct gene expression. MYC levels are tightly regulated, and deregulation of MYC has been associated with numerous diseases including cancer. Understanding the comprehensive biology of MYC under physiological conditions is an utmost necessity to demark biological functions of MYC from its pathological functions. Here we review the recent advances in biological mechanisms, functions, and regulation of MYC. We also emphasize the role of MYC as a global transcription amplifier.

Regulation of CTCF loop formation during pancreatic cell differentiation

Transcription reprogramming during cell differentiation involves targeting enhancers to genes responsible for establishment of cell fates. To understand the contribution of CTCF-mediated chromatin organization to cell lineage commitment, we analyzed 3D chromatin architecture during the differentiation of human embryonic stem cells into pancreatic islet organoids. We find that CTCF loops are formed and disassembled at different stages of the differentiation process by either recruitment of CTCF to new anchor sites or use of pre-existing sites not previously involved in loop formation. Recruitment of CTCF to new sites in the genome involves demethylation of H3K9me3 to H3K9me2, demethylation of DNA, recruitment of pioneer factors, and positioning of nucleosomes flanking the new CTCF sites. Existing CTCF sites not involved in loop formation become functional loop anchors via the establishment of new cohesin loading sites containing NIPBL and YY1 at sites between the new anchors. In both cases, formation of new CTCF loops leads to strengthening of enhancer promoter interactions and increased transcription of genes adjacent to loop anchors. These results suggest an important role for CTCF and cohesin in controlling gene expression during cell differentiation.

Defining super-enhancers by highly ranked histone H4 multi-acetylation levels identifies transcription factors associated with glioblastoma stem-like properties

BACKGROUND

Super-enhancers (SEs), which activate genes involved in cell-type specificity, have mainly been defined as genomic regions with top-ranked enrichment(s) of histone H3 with acetylated K27 (H3K27ac) and/or transcription coactivator(s) including a bromodomain and extra-terminal domain (BET) family protein, BRD4. However, BRD4 preferentially binds to multi-acetylated histone H4, typically with acetylated K5 and K8 (H4K5acK8ac), leading us to hypothesize that SEs should be defined by high H4K5acK8ac enrichment at least as well as by that of H3K27ac.

RESULTS

Here, we conducted genome-wide profiling of H4K5acK8ac and H3K27ac, BRD4 binding, and the transcriptome by using a BET inhibitor, JQ1, in three human glial cell lines. When SEs were defined as having the top ranks for H4K5acK8ac or H3K27ac signal, 43% of H4K5acK8ac-ranked SEs were distinct from H3K27ac-ranked SEs in a glioblastoma stem-like cell (GSC) line. CRISPR-Cas9-mediated deletion of the H4K5acK8ac-preferred SEs associated with MYCN and NFIC decreased the stem-like properties in GSCs.

CONCLUSIONS

Collectively, our data highlights H4K5acK8ac's utility for identifying genes regulating cell-type specificity.

Multiple lineage-specific epigenetic landscapes at the antigen receptor loci

Antigen receptors (AgRs) expressed on B and T cells provide the adaptive immune system with ability to detect numerous foreign antigens. Epigenetic features of B cell receptor (BCR) and T cell receptor (TCR) genes were previously studied in lymphocytes, but little is known about their epigenetic features in other cells. Here, we explored histone modifications and transcription markers at the BCR and TCR loci in lymphocytes (pro-B, DP T cells, and mature CD4+ T cells), compared to embryonic stem (ES) cells and neurons. In B cells, the BCR loci exhibited active histone modifications and transcriptional markers indicative of active loci. Similar results were observed at the TCR loci in T cells. All loci were largely inactive in neurons. Surprisingly, in ES cells all AgR loci displayed a high degree of active histone modifications and markers of active transcription. Locations of these active histone modifications in ES cells were largely distinct from those in pro-B cells, and co-localized at numerous binding locations for transcription factors Oct4, Sox2, and Nanog. ES and pro-B cells also showed distinct binding patterns for the ubiquitous transcription factor YY1 and chromatin remodeler Brg1. On the contrary, there were many overlapping CCCTC-binding factor (CTCF) binding patterns when comparing ES cells, pro-B cells, and neurons. Our study identifies epigenetic features in ES cells and lymphocytes that may be related to ES cell pluripotency and lymphocyte tissue-specific activation at the AgR loci.

CTCF and R-loops are boundaries of cohesin-mediated DNA looping

Cohesin and CCCTC-binding factor (CTCF) are key regulatory proteins of three-dimensional (3D) genome organization. Cohesin extrudes DNA loops that are anchored by CTCF in a polar orientation. Here, we present direct evidence that CTCF binding polarity controls cohesin-mediated DNA looping. Using single-molecule imaging, we demonstrate that a critical N-terminal motif of CTCF blocks cohesin translocation and DNA looping. The cryo-EM structure of the cohesin-CTCF complex reveals that this CTCF motif ahead of zinc fingers can only reach its binding site on the STAG1 cohesin subunit when the N terminus of CTCF faces cohesin. Remarkably, a C-terminally oriented CTCF accelerates DNA compaction by cohesin. DNA-bound Cas9 and Cas12a ribonucleoproteins are also polar cohesin barriers, indicating that stalling may be intrinsic to cohesin itself. Finally, we show that RNA-DNA hybrids (R-loops) block cohesin-mediated DNA compaction in vitro and are enriched with cohesin subunits in vivo, likely forming TAD boundaries.

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.

CapsNetYY1: identifying YY1-mediated chromatin loops based on a capsule network architecture

BACKGROUND

Previous studies have identified that chromosome structure plays a very important role in gene control. The transcription factor Yin Yang 1 (YY1), a multifunctional DNA binding protein, could form a dimer to mediate chromatin loops and active enhancer-promoter interactions. The deletion of YY1 or point mutations at the YY1 binding sites significantly inhibit the enhancer-promoter interactions and affect gene expression. To date, only a few computational methods are available for identifying YY1-mediated chromatin loops.

RESULTS

We proposed a novel model named CapsNetYY1, which was based on capsule network architecture to identify whether a pair of YY1 motifs can form a chromatin loop. Firstly, we encode the DNA sequence using one-hot encoding method. Secondly, multi-scale convolution layer is used to extract local features of the sequence, and bidirectional gated recurrent unit is used to learn the features across time steps. Finally, capsule networks (convolution capsule layer and digital capsule layer) used to extract higher level features and recognize YY1-mediated chromatin loops. Compared with DeepYY1, the only prediction for YY1-mediated chromatin loops, our model CapsNetYY1 achieved the better performance on the independent datasets (AUC [Formula: see text]).

CONCLUSION

The results indicate that CapsNetYY1 is an excellent method for identifying YY1-mediated chromatin loops. We believe that the CapsNetYY1 method will be used for predictive classification of other DNA sequences.

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.

Systematic mapping and modeling of 3D enhancer-promoter interactions in early mouse embryonic lineages reveal regulatory principles that determine the levels and cell-type specificity of gene expression

Mammalian embryogenesis commences with two pivotal and binary cell fate decisions that give rise to three essential lineages, the trophectoderm (TE), the epiblast (EPI) and the primitive endoderm (PrE). Although key signaling pathways and transcription factors that control these early embryonic decisions have been identified, the non-coding regulatory elements via which transcriptional regulators enact these fates remain understudied. To address this gap, we have characterized, at a genome-wide scale, enhancer activity and 3D connectivity in embryo-derived stem cell lines that represent each of the early developmental fates. We observed extensive enhancer remodeling and fine-scale 3D chromatin rewiring among the three lineages, which strongly associate with transcriptional changes, although there are distinct groups of genes that are irresponsive to topological changes. In each lineage, a high degree of connectivity or "hubness" positively correlates with levels of gene expression and enriches for cell-type specific and essential genes. Genes within 3D hubs also show a significantly stronger probability of coregulation across lineages, compared to genes in linear proximity or within the same contact domains. By incorporating 3D chromatin features, we build a novel predictive model for transcriptional regulation (3D-HiChAT), which outperformed models that use only 1D promoter or proximal variables in predicting levels and cell-type specificity of gene expression. Using 3D-HiChAT, we performed genome-wide in silico perturbations to nominate candidate functional enhancers and hubs in each cell lineage, and with CRISPRi experiments we validated several novel enhancers that control expression of one or more genes in their respective lineages. Our study comprehensively identifies 3D regulatory hubs associated with the earliest mammalian lineages and describes their relationship to gene expression and cell identity, providing a framework to understand lineage-specific transcriptional behaviors.

Hi-TrAC detects active sub-TADs and reveals internal organizations of super-enhancers

The spatial folding of eukaryotic genome plays a key role in genome function. We report here that our recently developed method, Hi-TrAC, which specializes in detecting chromatin loops among accessible genomic regions, can detect active sub-TADs with a median size of 100 kb, most of which harbor one or two cell specifically expressed genes and regulatory elements such as super-enhancers organized into nested interaction domains. These active sub-TADs are characterized by highly enriched histone mark H3K4me1 and chromatin-binding proteins, including Cohesin complex. Deletion of selected sub-TAD boundaries have different impacts, such as decreased chromatin interaction and gene expression within the sub-TADs or compromised insulation between the sub-TADs, depending on the specific chromatin environment. We show that knocking down core subunit of the Cohesin complex using shRNAs in human cells or decreasing the H3K4me1 modification by deleting the H3K4 methyltransferase Mll4 gene in mouse Th17 cells disrupted the sub-TADs structure. Our data also suggest that super-enhancers exist as an equilibrium globule structure, while inaccessible chromatin regions exist as a fractal globule structure. In summary, Hi-TrAC serves as a highly sensitive and inexpensive approach to study dynamic changes of active sub-TADs, providing more explicit insights into delicate genome structures and functions.

Integrating Omics and CRISPR Technology for Identification and Verification of Genomic Safe Harbor Loci in the Chicken Genome

BACKGROUND

One of the most prominent questions in the field of transgenesis is 'Where in the genome to integrate a transgene?'. Escape from epigenetic silencing and promoter shutdown of the transgene needs reliable genomic safe harbor (GSH) loci. Advances in genome engineering technologies combined with multi-omics bioinformatics data have enabled rational evaluation of GSH loci in the host genome. Currently, no validated GSH loci have been evaluated in the chicken genome.

RESULTS

Here, we analyzed and experimentally examined two GSH loci in the genome of chicken cells. To this end, putative GSH loci including chicken HIPP-like (cHIPP; between DRG1 and EIF4ENIF1 genes) and chicken ROSA-like (cROSA; upstream of the THUMPD3 gene) were predicted using multi-omics bioinformatics data. Then, the durable expression of the transgene was validated by experimental characterization of continuously-cultured isogenous cell clones harboring DsRed2-ΔCMV-EGFP cassette in the predicted loci. The weakened form of the CMV promoter (ΔCMV) allowed the precise evaluation of GSH loci in a locus-dependent manner compared to the full-length CMV promoter.

CONCLUSIONS

cHIPP and cROSA loci introduced in this study can be reliably exploited for consistent bio-manufacturing of recombinant proteins in the genetically-engineered chickens. Also, results showed that the genomic context dictates the expression of transgene controlled by ΔCMV in GSH loci.

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.

Regulatory architecture of housekeeping genes is driven by promoter assemblies

Genes that are key to cell identity are generally regulated by cell-type-specific enhancer elements bound by transcription factors, some of which facilitate looping to distant gene promoters. In contrast, genes that encode housekeeping functions, whose regulation is essential for normal cell metabolism and growth, generally lack interactions with distal enhancers. We find that Ronin (Thap11) assembles multiple promoters of housekeeping and metabolic genes to regulate gene expression. This behavior is analogous to how enhancers are brought together with promoters to regulate cell identity genes. Thus, Ronin-dependent promoter assemblies provide a mechanism to explain why housekeeping genes can forgo distal enhancer elements and why Ronin is important for cellular metabolism and growth control. We propose that clustering of regulatory elements is a mechanism common to cell identity and housekeeping genes but is accomplished by different factors binding distinct control elements to establish enhancer-promoter or promoter-promoter interactions, respectively.

Loop stacking organizes genome folding from TADs to chromosomes

Although population-level analyses revealed significant roles for CTCF and cohesin in mammalian genome organization, their contributions at the single-cell level remain incompletely understood. Here, we used a super-resolution microscopy approach to measure the effects of removal of CTCF or cohesin in mouse embryonic stem cells. Single-chromosome traces revealed cohesin-dependent loops, frequently stacked at their loop anchors forming multi-way contacts (hubs), bridging across TAD boundaries. Despite these bridging interactions, chromatin in intervening TADs was not intermixed, remaining separated in distinct loops around the hub. At the multi-TAD scale, steric effects from loop stacking insulated local chromatin from ultra-long range (>4 Mb) contacts. Upon cohesin removal, the chromosomes were more disordered and increased cell-cell variability in gene expression. Our data revise the TAD-centric understanding of CTCF and cohesin and provide a multi-scale, structural picture of how they organize the genome on the single-cell level through distinct contributions to loop stacking.

Novel evidence that the ABO blood group shapes erythropoiesis and results in higher hematocrit for blood group B carriers

The ABO blood group (BG) system is of great importance for blood transfusion and organ transplantation. Since the same transcription factors (TFs) and microRNAs (miRNAs) govern the expression of ABO BG antigens and regulate erythropoiesis, we hypothesized functional connections between both processes. We found significantly higher hemoglobin and hematocrit values in BG B blood donors compared to BG A. Furthermore, we observed that erythropoiesis in BG B hematopoietic stem/progenitor cells (HSPCs) was accelerated compared to BG A HSPCs. Specifically, BG B HSPCs yielded more lineage-specific progenitors in a shorter time (B: 31.3 ± 2.2% vs. A: 22.5 ± 3.0%). Moreover, non-BG A individuals exhibited more terminally differentiated RBCs with higher enucleation rates containing more hemoglobin compared to BG A. Additionally, we detected increased levels of miRNA-215-5p and -182-5p and decreased expression of their target TFs RUNX1 and HES-1 mRNAs in erythroid BG B precursor cells compared to BG A. This highlights the important roles of these factors for the disappearance of differentiation-specific glycan antigens and the appearance of cancer-specific glycan antigens. Our work contributes to a deeper understanding of erythropoiesis gene regulatory networks and identifies its interference with BG-specific gene expression regulations particularly in diseases, where ABO BGs determine treatment susceptibility and disease progression.

3D enhancer-promoter interactions and multi-connected hubs: Organizational principles and functional roles

The spatiotemporal control of gene expression is dependent on the activity of cis-acting regulatory sequences, called enhancers, which regulate target genes over variable genomic distances and, often, by skipping intermediate promoters, suggesting mechanisms that control enhancer-promoter communication. Recent genomics and imaging technologies have revealed highly complex enhancer-promoter interaction networks, whereas advanced functional studies have started interrogating the forces behind the physical and functional communication among multiple enhancers and promoters. In this review, we first summarize our current understanding of the factors involved in enhancer-promoter communication, with a particular focus on recent papers that have revealed new layers of complexities to old questions. In the second part of the review, we focus on a subset of highly connected enhancer-promoter "hubs" and discuss their potential functions in signal integration and gene regulation, as well as the putative factors that might determine their dynamics and assembly.

SATB1 regulates 3D genome architecture in T cells by constraining chromatin interactions surrounding CTCF-binding sites

Special AT-rich sequence binding protein 1 (SATB1) has long been proposed to act as a global chromatin loop organizer in T cells. However, the exact functions of SATB1 in spatial genome organization remain elusive. Here we show that the depletion of SATB1 in human and murine T cells leads to transcriptional dysregulation for genes involved in T cell activation, as well as alterations of 3D genome architecture at multiple levels, including compartments, topologically associating domains, and loops. Importantly, SATB1 extensively colocalizes with CTCF throughout the genome. Depletion of SATB1 leads to increased chromatin contacts among and across the SATB1/CTCF co-occupied sites, thereby affecting the transcription of critical regulators of T cell activation. The loss of SATB1 does not affect CTCF occupancy but significantly reduces the retention of CTCF in the nuclear matrix. Collectively, our data show that SATB1 contributes to 3D genome organization by constraining chromatin topology surrounding CTCF-binding sites.

Epigenetic regulation by ASXL1 in myeloid malignancies

Myeloid malignancies are clonal hematopoietic disorders that are comprised of a spectrum of genetically heterogeneous disorders, including myelodysplastic syndromes (MDS), myeloproliferative neoplasms (MPN), chronic myelomonocytic leukemia (CMML), and acute myeloid leukemia (AML). Myeloid malignancies are characterized by excessive proliferation, abnormal self-renewal, and/or differentiation defects of hematopoietic stem cells (HSCs) and myeloid progenitor cells hematopoietic stem/progenitor cells (HSPCs). Myeloid malignancies can be caused by genetic and epigenetic alterations that provoke key cellular functions, such as self-renewal, proliferation, biased lineage commitment, and differentiation. Advances in next-generation sequencing led to the identification of multiple mutations in myeloid neoplasms, and many new gene mutations were identified as key factors in driving the pathogenesis of myeloid malignancies. The polycomb protein ASXL1 was identified to be frequently mutated in all forms of myeloid malignancies, with mutational frequencies of 20%, 43%, 10%, and 20% in MDS, CMML, MPN, and AML, respectively. Significantly, ASXL1 mutations are associated with a poor prognosis in all forms of myeloid malignancies. The fact that ASXL1 mutations are associated with poor prognosis in patients with CMML, MDS, and AML, points to the possibility that ASXL1 mutation is a key factor in the development of myeloid malignancies. This review summarizes the recent advances in understanding myeloid malignancies with a specific focus on ASXL1 mutations.

Compartmentalization with nuclear landmarks yields random, yet precise, genome organization

The 3D organization of eukaryotic genomes plays an important role in genome function. While significant progress has been made in deciphering the folding mechanisms of individual chromosomes, the principles of the dynamic large-scale spatial arrangement of all chromosomes inside the nucleus are poorly understood. We use polymer simulations to model the diploid human genome compartmentalization relative to nuclear bodies such as nuclear lamina, nucleoli, and speckles. We show that a self-organization process based on a cophase separation between chromosomes and nuclear bodies can capture various features of genome organization, including the formation of chromosome territories, phase separation of A/B compartments, and the liquid property of nuclear bodies. The simulated 3D structures quantitatively reproduce both sequencing-based genomic mapping and imaging assays that probe chromatin interaction with nuclear bodies. Importantly, our model captures the heterogeneous distribution of chromosome positioning across cells while simultaneously producing well-defined distances between active chromatin and nuclear speckles. Such heterogeneity and preciseness of genome organization can coexist due to the nonspecificity of phase separation and the slow chromosome dynamics. Together, our work reveals that the cophase separation provides a robust mechanism for us to produce functionally important 3D contacts without requiring thermodynamic equilibration that can be difficult to achieve.

CTCF mediates CD8+ effector differentiation through dynamic redistribution and genomic reorganization

Differentiation of effector CD8+ T cells is instructed by stably and dynamically expressed transcription regulators. Here we show that naive-to-effector differentiation was accompanied by dynamic CTCF redistribution and extensive chromatin architectural changes. Upon CD8+ T cell activation, CTCF acquired de novo binding sites and anchored novel chromatin interactions, and these changes were associated with increased chromatin accessibility and elevated expression of cytotoxic program genes including Tbx21, Ifng, and Klrg1. CTCF was also evicted from its ex-binding sites in naive state, with concomitantly reduced chromatin interactions in effector cells, as observed at memory precursor-associated genes including Il7r, Sell, and Tcf7. Genetic ablation of CTCF indeed diminished cytotoxic gene expression, but paradoxically elevated expression of memory precursor genes. Comparative Hi-C analysis revealed that key memory precursor genes were harbored within insulated neighborhoods demarcated by constitutive CTCF binding, and their induction was likely due to disrupted CTCF-dependent insulation. CTCF thus promotes cytotoxic effector differentiation by integrating local chromatin accessibility control and higher-order genomic reorganization.

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.

A DNA tumor virus globally reprograms host 3D genome architecture to achieve immortal growth

Epstein-Barr virus (EBV) immortalization of resting B lymphocytes (RBLs) to lymphoblastoid cell lines (LCLs) models human DNA tumor virus oncogenesis. RBL and LCL chromatin interaction maps are compared to identify the spatial and temporal genome architectural changes during EBV B cell transformation. EBV induces global genome reorganization where contact domains frequently merge or subdivide during transformation. Repressed B compartments in RBLs frequently switch to active A compartments in LCLs. LCLs gain 40% new contact domain boundaries. Newly gained LCL boundaries have strong CTCF binding at their borders while in RBLs, the same sites have much less CTCF binding. Some LCL CTCF sites also have EBV nuclear antigen (EBNA) leader protein EBNALP binding. LCLs have more local interactions than RBLs at LCL dependency factors and super-enhancer targets. RNA Pol II HiChIP and FISH of RBL and LCL further validate the Hi-C results. EBNA3A inactivation globally alters LCL genome interactions. EBNA3A inactivation reduces CTCF and RAD21 DNA binding. EBNA3C inactivation rewires the looping at the CDKN2A/B and AICDA loci. Disruption of a CTCF site at AICDA locus increases AICDA expression. These data suggest that EBV controls lymphocyte growth by globally reorganizing host genome architecture to facilitate the expression of key oncogenes.

Species-specific regulation of XIST by the JPX/FTX orthologs

X chromosome inactivation (XCI) is an essential process, yet it initiates with remarkable diversity in various mammalian species. XIST, the main trigger of XCI, is controlled in the mouse by an interplay of lncRNA genes (LRGs), some of which evolved concomitantly to XIST and have orthologues across all placental mammals. Here, we addressed the functional conservation of human orthologues of two such LRGs, FTX and JPX. By combining analysis of single-cell RNA-seq data from early human embryogenesis with various functional assays in matched human and mouse pluripotent stem- or differentiated post-XCI cells, we demonstrate major functional differences for these orthologues between species, independently of primary sequence conservation. While the function of FTX is not conserved in humans, JPX stands as a major regulator of XIST expression in both species. However, we show that different entities of JPX control the production of XIST at various steps depending on the species. Altogether, our study highlights the functional versatility of LRGs across evolution, and reveals that functional conservation of orthologous LRGs may involve diversified mechanisms of action. These findings represent a striking example of how the evolvability of LRGs can provide adaptative flexibility to constrained gene regulatory networks.

Network Approaches to Uncover Pathogenesis and Therapeutic Targets of Inflammatory Bowel Diseases

Inflammatory bowel diseases (IBD) are currently recognized to involve chronic intestinal inflammation in genetically susceptible individuals. Patients with IBD mainly develop gastrointestinal inflammation, but it is sometimes accompanied by extraintestinal manifestations such as arthritis, erythema nodosum, episcleritis, pyoderma gangrenosum, uveitis, and primary sclerosing cholangitis. These clinical aspects imply the importance of interorgan networks in IBD. In the gastrointestinal tract, immune cells are influenced by multiple local environmental factors including microbiota, dietary environment, and intercellular networks, which further alter molecular networks in immune cells. Therefore, deciphering networks at interorgan, intercellular, and intracellular levels should help to obtain a comprehensive understanding of IBD. This review focuses on the intestinal immune system, which governs the physiological and pathological functions of the digestive system in harmony with the other organs.

Auxin-inducible degron 2 system deciphers functions of CTCF domains in transcriptional regulation

BACKGROUND

CTCF is a well-established chromatin architectural protein that also plays various roles in transcriptional regulation. While CTCF biology has been extensively studied, how the domains of CTCF function to regulate transcription remains unknown. Additionally, the original auxin-inducible degron 1 (AID1) system has limitations in investigating the function of CTCF.

RESULTS

We employ an improved auxin-inducible degron technology, AID2, to facilitate the study of acute depletion of CTCF while overcoming the limitations of the previous AID system. As previously observed through the AID1 system and steady-state RNA analysis, the new AID2 system combined with SLAM-seq confirms that CTCF depletion leads to modest nascent and steady-state transcript changes. A CTCF domain sgRNA library screening identifies the zinc finger (ZF) domain as the region within CTCF with the most functional relevance, including ZFs 1 and 10. Removal of ZFs 1 and 10 reveals genomic regions that independently require these ZFs for DNA binding and transcriptional regulation. Notably, loci regulated by either ZF1 or ZF10 exhibit unique CTCF binding motifs specific to each ZF.

CONCLUSIONS

By extensively comparing the AID1 and AID2 systems for CTCF degradation in SEM cells, we confirm that AID2 degradation is superior for achieving miniAID-tagged protein degradation without the limitations of the AID1 system. The model we create that combines AID2 depletion of CTCF with exogenous overexpression of CTCF mutants allows us to demonstrate how peripheral ZFs intricately orchestrate transcriptional regulation in a cellular context for the first time.

Аpplication of massive parallel reporter analysis in biotechnology and medicine

The development and functioning of an organism relies on tissue-specific gene programs. Genome regulatory elements play a key role in the regulation of such programs, and disruptions in their function can lead to the development of various pathologies, including cancers, malformations and autoimmune diseases. The emergence of high-throughput genomic studies has led to massively parallel reporter analysis (MPRA) methods, which allow the functional verification and identification of regulatory elements on a genome-wide scale. Initially MPRA was used as a tool to investigate fundamental aspects of epigenetics, but the approach also has great potential for clinical and practical biotechnology. Currently, MPRA is used for validation of clinically significant mutations, identification of tissue-specific regulatory elements, search for the most promising loci for transgene integration, and is an indispensable tool for creating highly efficient expression systems, the range of application of which extends from approaches for protein development and design of next-generation therapeutic antibody superproducers to gene therapy. In this review, the main principles and areas of practical application of high-throughput reporter assays will be discussed.

The Interplay of Transcription and Genome Topology Programs T Cell Development and Differentiation

T cells are essential for mounting defense against various pathogens and malignantly transformed cells. Thymic development and peripheral T cell differentiation are highly orchestrated biological processes that require precise gene regulation. Higher-order genome organization on multiple scales, in the form of chromatin loops, topologically associating domains and compartments, provides pivotal control of T cell gene expression. CTCF and the cohesin machinery are ubiquitously expressed architectural proteins responsible for establishing chromatin structures. Recent studies indicate that transcription factors, such as T lineage-defining Tcf1 and TCR-induced Batf, may have intrinsic ability and/or engage CTCF to shape chromatin architecture. In this article, we summarize current knowledge on the dynamic changes in genome topology that underlie normal or leukemic T cell development, CD4+ helper T cell differentiation, and CD8+ cytotoxic T cell functions. The knowledge lays a solid foundation for elucidating the causative link of spatial chromatin configuration to transcriptional and functional output in T cells.

HNF1B-driven three-dimensional chromatin structure for molecular classification in pancreatic cancers

The molecular subtypes of pancreatic cancer (PC), either classical/progenitor-like or basal/squamous-like, are currently a major topic of research because of their direct association with clinical outcomes. Some transcription factors (TFs) have been reported to be associated with these subtypes. However, the mechanisms by which these molecular signatures of PCs are established remain unknown. Epigenetic regulatory processes, supported by dynamic changes in the chromatin structure, are essential for transcriptional profiles. Previously, we reported the importance of open chromatin profiles in the biological features and transcriptional status of PCs. Here, we aimed to analyze the relationships between three-dimensional (3D) genome structures and the molecular subtypes of human PCs using Hi-C analysis. We observed a correlation of the specific elements of 3D genome modules, including compartments, topologically associating domains, and enhancer-promoter loops, with the expression of related genes. We focused on HNF1B, a TF that is implicated in the progenitor subtype. Forced expression of HNF1B in squamous-type PC organoids induced the upregulation and downregulation of genes associated with progenitor and squamous subtypes, respectively. Long-range genomic interactions induced by HNF1B were accompanied by compartment modulation and H3K27ac redistribution. We also found that these HNF1B-induced changes in subtype-related gene expression required an intrinsically disordered region, suggesting a possible involvement of phase separation in compartment modulation. Thus, mapping of 3D structural changes induced by TFs, such as HNF1B, may become a useful resource for further understanding the molecular features of PCs.

3D-Epigenomic Regulation of Gene Transcription in Hepatocellular Carcinoma

The fundamental cause of transcription dysregulation in hepatocellular carcinoma (HCC) remains elusive. To investigate the underlying mechanisms, comprehensive 3D-epigenomic analyses are performed in cellular models of THLE2 (a normal hepatocytes cell line) and HepG2 (a hepatocellular carcinoma cell line) using integrative approaches for chromatin topology, genomic and epigenomic variation, and transcriptional output. Comparing the 3D-epigenomes in THLE2 and HepG2 reveal that most HCC-associated genes are organized in complex chromatin interactions mediated by RNA polymerase II (RNAPII). Incorporation of genome-wide association studies (GWAS) data enables the identification of non-coding genetic variants that are enriched in distal enhancers connecting to the promoters of HCC-associated genes via long-range chromatin interactions, highlighting their functional roles. Interestingly, CTCF binding and looping proximal to HCC-associated genes appear to form chromatin architectures that overarch RNAPII-mediated chromatin interactions. It is further demonstrated that epigenetic variants by DNA hypomethylation at a subset of CTCF motifs proximal to HCC-associated genes can modify chromatin topological configuration, which in turn alter RNAPII-mediated chromatin interactions and lead to dysregulation of transcription. Together, the 3D-epigenomic analyses provide novel insights of multifaceted interplays involving genetics, epigenetics, and chromatin topology in HCC cells.

Comparison of CRISPR-Cas9-mediated megabase-scale genome deletion methods in mouse embryonic stem cells

The genome contains large functional units ranging in size from hundreds of kilobases to megabases, such as gene clusters and topologically associating domains. To analyse these large functional units, the technique of deleting the entire functional unit is effective. However, deletion of such large regions is less efficient than conventional genome editing, especially in cultured cells, and a method that can ensure success is anticipated. Here, we compared methods to delete the 2.5-Mb Krüppel-associated box zinc finger protein (KRAB-ZFP) gene cluster in mouse embryonic stem cells using CRISPR-Cas9. Three methods were used: first, deletion by non-homologous end joining (NHEJ); second, homology-directed repair (HDR) using a single-stranded oligodeoxynucleotide (ssODN); and third, HDR employing targeting vectors with a selectable marker and 1-kb homology arms. NHEJ-mediated deletion was achieved in 9% of the transfected cells. Inversion was also detected at similar efficiency. The deletion frequency of NHEJ and HDR was found to be comparable when the ssODN was transfected. Deletion frequency was highest when targeting vectors were introduced, with deletions occurring in 31-63% of the drug-resistant clones. Biallelic deletion was observed when targeting vectors were used. This study will serve as a benchmark for the introduction of large deletions into the genome.

Integrating extrusion complex-associated pattern to predict cell type-specific long-range chromatin loops

The chromatin loop plays a critical role in the study of gene expression and disease. Supervised learning-based algorithms to predict the chromatin loops require large priori information to satisfy the model construction, while the prediction sensitivity of unsupervised learning-based algorithms is still unsatisfactory. Therefore, we propose an unsupervised algorithm, Ecomap-loop. It takes advantage of extrusion complex-associated patterns, including CTCF, RAD21, and SMC enrichments, as well as the orientation distribution of CTCF motif of loops to build feature matrices; then the eigen decomposition model is employed to obtain the cell type-specific loops. We compare the performance of Ecomap-loop with the state-of-the-art unsupervised algorithm using Hi-C, ChIA-PET, expression quantitative trait locus (eQTL), and CRISPR interference (CRISPRi) screen data; the results show that Ecomap-loop achieves the best in four cell types. In addition, the functional analysis reveals the ability of Ecomap-loop to predict active functionality-related and cell type-specific loops.

Systematic discovery and functional dissection of enhancers needed for cancer cell fitness and proliferation

A scarcity of functionally validated enhancers in the human genome presents a significant hurdle to understanding how these cis-regulatory elements contribute to human diseases. We carry out highly multiplexed CRISPR-based perturbation and sequencing to identify enhancers required for cell proliferation and fitness in 10 human cancer cell lines. Our results suggest that the cell fitness enhancers, unlike their target genes, display high cell-type specificity of chromatin features. They typically adopt a modular structure, comprised of activating elements enriched for motifs of oncogenic transcription factors, surrounded by repressive elements enriched for motifs recognized by transcription factors with tumor suppressor functions. We further identify cell fitness enhancers that are selectively accessible in clinical tumor samples, and the levels of chromatin accessibility are associated with patient survival. These results reveal functional enhancers across multiple cancer cell lines, characterize their context-dependent chromatin organization, and yield insights into altered transcription programs in cancer cells.

MyoD-Induced Trans-Differentiation: A Paradigm for Dissecting the Molecular Mechanisms of Cell Commitment, Differentiation and Reprogramming

The discovery of the skeletal muscle-specific transcription factor MyoD represents a milestone in the field of transcriptional regulation during differentiation and cell-fate reprogramming. MyoD was the first tissue-specific factor found capable of converting non-muscle somatic cells into skeletal muscle cells. A unique feature of MyoD, with respect to other lineage-specific factors able to drive trans-differentiation processes, is its ability to dramatically change the cell fate even when expressed alone. The present review will outline the molecular strategies by which MyoD reprograms the transcriptional regulation of the cell of origin during the myogenic conversion, focusing on the activation and coordination of a complex network of co-factors and epigenetic mechanisms. Some molecular roadblocks, found to restrain MyoD-dependent trans-differentiation, and the possible ways for overcoming these barriers, will also be discussed. Indeed, they are of critical importance not only to expand our knowledge of basic muscle biology but also to improve the generation skeletal muscle cells for translational research.

High-resolution genome topology of human retina uncovers super enhancer-promoter interactions at tissue-specific and multifactorial disease loci

Chromatin organization and enhancer-promoter contacts establish unique spatiotemporal gene expression patterns in distinct cell types. Non-coding genetic variants can influence cellular phenotypes by modifying higher-order transcriptional hubs and consequently gene expression. To elucidate genomic regulation in human retina, we mapped chromatin contacts at high resolution and integrated with super-enhancers (SEs), histone marks, binding of CTCF and select transcription factors. We show that topologically associated domains (TADs) with central SEs exhibit stronger insulation and augmented contact with retinal genes relative to TADs with edge SEs. Merging genome-wide expression quantitative trait loci (eQTLs) with topology map reveals physical links between 100 eQTLs and corresponding eGenes associated with retinal neurodegeneration. Additionally, we uncover candidate genes for susceptibility variants linked to age-related macular degeneration and glaucoma. Our study of high-resolution genomic architecture of human retina provides insights into genetic control of tissue-specific functions, suggests paradigms for missing heritability, and enables the dissection of common blinding disease phenotypes.

Analysis of chromatin data supports a role for CD14+ monocytes/macrophages in mediating genetic risk for juvenile idiopathic arthritis

Introduction

Genome wide association studies (GWAS) have identified multiple regions that confer genetic risk for the polyarticular/oligoarticular forms of juvenile idiopathic arthritis (JIA). However, genome-wide scans do not identify the cells impacted by genetic polymorphisms on the risk haplotypes or the genes impacted by those variants. We have shown that genetic variants driving JIA risk are likely to affect both innate and adaptive immune functions. We provide additional evidence that JIA risk variants impact innate immunity.

Materials and methods

We queried publicly available H3K4me1/H3K27ac ChIP-seq data in CD14+ monocytes to determine whether the linkage disequilibrium (LD) blocks incorporating the SNPs that tag JIA risk loci showed enrichment for these epigenetic marks. We also queried monocyte/macrophage GROseq data, a functional readout of active enhancers. We defined the topologically associated domains (TADs) encompassing enhancers on the risk haplotypes and identified genes within those TADs expressed in monocytes. We performed ontology analyses of these genes to identify cellular processes that may be impacted by these variants. We also used whole blood RNAseq data from the Genotype-Tissue Expression (GTEx) data base to determine whether SNPs lying within monocyte GROseq peaks influence plausible target genes within the TADs encompassing the JIA risk haplotypes.

Results

The LD blocks encompassing the JIA genetic risk regions were enriched for H3K4me1/H3K27ac ChIPseq peaks (p=0.00021 and p=0.022) when compared to genome background. Eleven and sixteen JIA were enriched for resting and activated macrophage GROseq peaks, respectively risk regions (p=0.04385 and p=0.00004). We identified 321 expressed genes within the TADs encompassing the JIA haplotypes in human monocytes. Ontological analysis of these genes showed enrichment for multiple immune functions. Finally, we found that SNPs lying within the GROseq peaks are strongly associated with expression levels of plausible target genes in human whole blood.

Conclusions

These findings support the idea that both innate and adaptive immunity are impacted by JIA genetic risk variants.

The transcriptional coactivator RUVBL2 regulates Pol II clustering with diverse transcription factors

RNA polymerase II (Pol II) apparatuses are compartmentalized into transcriptional clusters. Whether protein factors control these clusters remains unknown. In this study, we find that the ATPase-associated with diverse cellular activities (AAA + ) ATPase RUVBL2 co-occupies promoters with Pol II and various transcription factors. RUVBL2 interacts with unphosphorylated Pol II in chromatin to promote RPB1 carboxy-terminal domain (CTD) clustering and transcription initiation. Rapid depletion of RUVBL2 leads to a decrease in the number of Pol II clusters and inhibits nascent RNA synthesis, and tethering RUVBL2 to an active promoter enhances Pol II clustering at the promoter. We also identify target genes that are directly linked to the RUVBL2-Pol II axis. Many of these genes are hallmarks of cancers and encode proteins with diverse cellular functions. Our results demonstrate an emerging activity for RUVBL2 in regulating Pol II cluster formation in the nucleus.

RNA polymerase II (Pol II) transcription factories play a central role in gene expression and 3D chromatin organization. Here, the authors demonstrate that RUVBL2 directly regulates Pol II clustering at active gene promoters.