3D Chromatin Architecture Research Articles

3D chromatin architecture, BRD4, and Mediator have distinct roles in regulating genome-wide transcriptional bursting and gene network.

Discontinuous transcription is evolutionarily conserved and a fundamental feature of gene regulation; yet, the exact mechanisms underlying transcriptional bursting are unresolved. Analyses of bursting transcriptome-wide have focused on the role of cis-regulatory elements, but other factors that regulate this process remain elusive. We applied mathematical modeling to single-cell RNA sequencing data to infer bursting dynamics transcriptome-wide under multiple conditions to identify possible molecular mechanisms. We found that Mediator complex subunit 26 (MED26) primarily regulates frequency, MYC regulates burst size, while cohesin and Bromodomain-containing protein 4 (BRD4) can modulate both. Despite comparable effects on RNA levels among these perturbations, acute depletion of MED26 had the most profound impact on the entire gene regulatory network, acting downstream of chromatin spatial architecture and without affecting TATA box-binding protein (TBP) recruitment. These results indicate that later steps in the initiation of transcriptional bursts are primary nodes for integrating gene networks in single cells.

Abstract We139: 3D Chromatin Architectures and Transcription Regulation in Diabetic Endothelial Dysfunction

Cardiovascular disease is the major cause of mortality in people with type 2 diabetes, primarily due to endothelial dysfunction. Recent studies have suggested that transcriptional reprogramming plays a central role in endothelial dysfunction, but its regulation in the context of diabetes remains poorly understood. In this study, we utilized 3D chromatin mapping analysis (H3K27ac HiChIP) to investigate the transcriptional regulation in aortic endothelial cells (db/db mice) and HUVEC cells treated with high glucose (HUVEC high-glu ). We found that these cells were characterized by a massive loss of enhancer-promoter loops, which was associated with transcriptional downregulation. Through transcription factor motif scan, ChIP-seq and ATAC-seq analysis, we uncovered that the lost loop anchors are demarcated by reduced JUNB (transcription factor) binding while their chromatin accessibility remains unchanged. Moreover, we demonstrated that overexpression of JUNB significantly restored enhancer-promoter loops lost in HUVEC high-glu and ameliorated endothelial dysfunctions both in vitro (cell scratch, Transwell), ex vivo (aorta relaxation/contraction) and in vivo (hindlimb ischemia) assays. We further showed that JUNB protein (but not RNA) expression is reduced in aortic endothelial cells (db/db), HUVEC high-glu and primary endothelial cells from diabetic patients. Through RIME (Rapid immunoprecipitation mass spectrometry of endogenous protein) assay, we identified that RBBP6 (RB Binding Protein 6, Ubiquitin Ligase) is a JUNB interacting protein in HUVEC high-glu . We further showed that overexpression of RBBP6 downregulated JUNB expression while proteasome inhibitor MG132 treatment largely abolished the inhibitory effect of RBBP6 on JUNB protein levels in normal HUVEC cells, suggesting that RPPB6 promoted JUNB degradation through ubiquitination. Notably, knockdown of RBBP6 in HUVEC high-glu significantly upregulated JUNB expression, enhancer-promoter loopings and improved endothelial dysfunction. Overall, our findings provide novel 3D chromatin insights into diabetic endothelial dysfunction and illustrate a model whereby RBBP6-mediated JUNB degradation leads to global loss of enhancer-promoter looping and transcriptional inhibition in diabetic endothelial dysfunction.

Evolutionary patterns and functional effects of 3D chromatin structures in butterflies with extensive genome rearrangements

Chromosome rearrangements may distort 3D chromatin architectures and thus change gene regulation, yet how 3D chromatin structures evolve in insects is largely unknown. Here, we obtain chromosome-level genomes for four butterfly species, Graphium cloanthus, Graphium sarpedon, Graphium eurypylus with 2n = 30, 40, and 60, respectively, and Papilio bianor with 2n = 60. Together with large-scale Hi-C data, we find that inter-chromosome rearrangements very rarely disrupted the pre-existing 3D chromatin structure of ancestral chromosomes. However, some intra-chromosome rearrangements changed 3D chromatin structures compared to the ancestral configuration. We find that new TADs and subTADs have emerged across the rearrangement sites where their adjacent compartments exhibit uniform types. Two intra-chromosome rearrangements altered Rel and lft regulation, potentially contributing to wing patterning differentiation and host plant choice. Notably, butterflies exhibited chromatin loops between Hox gene cluster ANT-C and BX-C, unlike Drosophila. Our CRISPR-Cas9 experiments in butterflies confirm that knocking out the CTCF binding site of the loops in BX-C affected the phenotypes regulated by Antp in ANT-C, resulting in legless larva. Our results reveal evolutionary patterns of insect 3D chromatin structures and provide evidence that 3D chromatin structure changes can play important roles in the evolution of traits.

Machine learning enables pan-cancer identification of mutational hotspots at persistent CTCF binding sites.

CCCTC-binding factor (CTCF) is an insulator protein that binds to a highly conserved DNA motif and facilitates regulation of three-dimensional (3D) nuclear architecture and transcription. CTCF binding sites (CTCF-BSs) reside in non-coding DNA and are frequently mutated in cancer. Our previous study identified a small subclass of CTCF-BSs that are resistant to CTCF knock down, termed persistent CTCF binding sites (P-CTCF-BSs). P-CTCF-BSs show high binding conservation and potentially regulate cell-type constitutive 3D chromatin architecture. Here, using ICGC sequencing data we made the striking observation that P-CTCF-BSs display a highly elevated mutation rate in breast and prostate cancer when compared to all CTCF-BSs. To address whether P-CTCF-BS mutations are also enriched in other cell-types, we developed CTCF-INSITE-a tool utilising machine learning to predict persistence based on genetic and epigenetic features of experimentally-determined P-CTCF-BSs. Notably, predicted P-CTCF-BSs also show a significantly elevated mutational burden in all 12 cancer-types tested. Enrichment was even stronger for P-CTCF-BS mutations with predicted functional impact to CTCF binding and chromatin looping. Using in vitro binding assays we validated that P-CTCF-BS cancer mutations, predicted to be disruptive, indeed reduced CTCF binding. Together this study reveals a new subclass of cancer specific CTCF-BS DNA mutations and provides insights into their importance in genome organization in a pan-cancer setting.

3D Chromatin Alteration by Disrupting β-Catenin/CBP Interaction Is Enriched with Insulin Signaling in Pancreatic Cancer.

The therapeutic potential of targeting the β-catenin/CBP interaction has been demonstrated in a variety of preclinical tumor models with a small molecule inhibitor, ICG-001, characterized as a β-catenin/CBP antagonist. Despite the high binding specificity of ICG-001 for the N-terminus of CBP, this β-catenin/CBP antagonist exhibits pleiotropic effects. Our recent studies found global changes in three-dimensional (3D) chromatin architecture in response to disruption of the β-catenin/CBP interaction in pancreatic cancer cells. However, an understanding of how the functional crosstalk between the antagonist and the β-catenin/CBP interaction affects changes in 3D chromatin architecture and, thereby, gene expression and downstream effects remains to be elucidated. Here, we perform Hi-C analyses on canonical and patient-derived pancreatic cancer cells before and after treatment with ICG-001. In addition to global alteration of 3D chromatin domains, we unexpectedly identify insulin signaling genes enriched in the altered chromatin domains. We further demonstrate that the chromatin loops associated with insulin signaling genes are significantly weakened after ICG-001 treatment. We finally elicit the deletion of a looping of IRS1-a key insulin signaling gene-significantly impeding pancreatic cancer cell growth, indicating that looping-mediated insulin signaling might act as an oncogenic pathway to promote pancreatic cancer progression. Our work shows that targeting aberrant insulin chromatin looping in pancreatic cancer might provide a therapeutic benefit.

Mariner: explore the Hi-Cs.

3D chromatin structure plays an important role in regulating gene expression and alterations to this structure can result in developmental abnormalities and disease. While genomic approaches like Hi-C and Micro-C can provide valuable insights in 3D chromatin architecture, the resulting datasets are extremely large and difficult to manipulate. Here, we present mariner, a rapid and memory efficient tool to extract, aggregate, and plot data from Hi-C matrices within the R/Bioconductor environment. Mariner simplifies the process of querying and extracting contacts from multiple Hi-C files using a parallel and block-processing approach. Modular functions allow complete workflow customization for advanced users, yet all-in-one functions are available for running the most common types of analyses. Finally, tight integration with existing Bioconductor infrastructure enables complete analysis and visualization of Hi-C data in R. Available on GitHub at https://github.com/EricSDavis/mariner and on Bioconductor at https://www.bioconductor.org/packages/release/bioc/html/mariner.html.

Three-dimensional chromatin landscapes in somatotroph tumour.

The three-dimensional (3D) genome architecture plays a critical role inregulating gene expression. However, the specific alterations in thisarchitecture within somatotroph tumors and their implications for gene expression remain largely unexplored. We employed Hi-C and RNA-seq analyses to compare the 3D genomic structures of somatotroph tumors with normal pituitary tissue. This comprehensive approachenabled the characterization of A/B compartments, topologically associateddomains (TADs), and chromatin loops, integrating these with gene expression patterns. We observed a decrease in both the frequency of chromosomal interactions andthe size of TADs in tumor tissue compared to normal tissue. Conversely, the number of TADs and chromatin loops was found to be increased in tumors. Integrated analysis of Hi-C and RNA-seq data demonstrated that changes inhigher-order chromat in structure were associated with alterations in gene expression. Specifically, genes in A compartments showed higher density and increased expression relative to those in B compartments. Moreover, the weakand enhanced insulation boundaries were identified, and the associated genes were enriched in the Wnt/β-Catenin signaling pathway. We identified the gainedand lost loops in tumor and integrated these differences with transcriptional changes to examine the functional relevance of the identified loops. Notably, we observed an enhanced insulation boundary and a greater number of loops in the TCF7L2 gene region within tumors, which was accompanied by an upregulation of TCF7L2 expression. Subsequently, TCF7L2 expression was confirmed through qRT-PCR, and upregulated TCF7L2 prompted cell proliferation and growth hormone (GH) secretion in vitro. Our results provide comprehensive 3D chromatin architecture maps of somatotroph tumors and offer a valuable resource for furthering the understanding of the underlying biology and mechanisms of gene expression regulation.

A CTCF-binding site in the Mdm1-Il22-Ifng locus shapes cytokine expression profiles and plays a critical role in early Th1 cell fate specification

Cytokine expression during Tcell differentiation is a highly regulated process that involves long-range promoter-enhancer and CTCF-CTCF contacts at cytokine loci. Here, we investigated the impact of dynamic chromatin loop formation within the topologically associating domain (TAD) in regulating the expression of interferon gamma (IFN-γ) and interleukin-22 (IL-22); these cytokine loci are closely located in the genome and are associated with complex enhancer landscapes, which are selectively active in type 1 and type 3 lymphocytes. In situ Hi-C analyses revealed inducible TADs that insulated Ifng and Il22 enhancers during Th1 cell differentiation. Targeted deletion of a 17bp boundary motif of these TADs imbalanced Th1- and Th17-associated immunity, both invitro and invivo, upon Toxoplasma gondii infection. In contrast, this boundary element was dispensable for cytokine regulation in natural killer cells. Our findings suggest that precise cytokine regulation relies on lineage- and developmental stage-specific interactions of 3D chromatin architectures and enhancer landscapes.

Building Haplotype-Resolved 3D Genome Maps of Chicken Skeletal Muscle.

Haplotype-resolved 3D chromatin architecture related to allelic differences in avian skeletal muscle development has not been addressed so far, although chicken husbandry for meat consumption has been prevalent feature of cultures on every continent for more than thousands of years. Here, high-resolution Hi-C diploid maps (1.2-kb maximum resolution) are generated for skeletal muscle tissues in chicken across three developmental stages (embryonic day 15 to day 30 post-hatching). The sequence features governing spatial arrangement of chromosomes and characterize homolog pairing in the nucleus, are identified. Multi-scale characterization of chromatin reorganization between stages from myogenesis in the fetus to myofiber hypertrophy after hatching show concordant changes in transcriptional regulation by relevant signaling pathways. Further interrogation of parent-of-origin-specific chromatin conformation supported that genomic imprinting is absent in birds. This study also reveals promoter-enhancer interaction (PEI) differences between broiler and layer haplotypes in skeletal muscle development-related genes are related to genetic variation between breeds, however, only a minority of breed-specific variations likely contribute to phenotypic divergence in skeletal muscle potentially via allelic PEI rewiring. Beyond defining the haplotype-specific 3D chromatin architecture in chicken, this study provides a rich resource for investigating allelic regulatory divergence among chicken breeds.

The sirtuin-associated human senescence program converges on the activation of placenta-specific gene PAPPA

Sirtuins are pro-longevity genes with chromatin modulation potential, but how these properties are connected is not well understood. Here, we generated a panel of isogeneic human stem cell lines with SIRT1-SIRT7 knockouts and found that any sirtuin deficiency leads to accelerated cellular senescence. Through large-scale epigenomic analyses, we show how sirtuin deficiency alters genome organization and that genomic regions sensitive to sirtuin deficiency are preferentially enriched in active enhancers, thereby promoting interactions within topologically associated domains and the formation of de novo enhancer-promoter loops. In all sirtuin-deficient human stem cell lines, we found that chromatin contacts are rewired to promote aberrant activation of the placenta-specific gene PAPPA, which controls the pro-senescence effects associated with sirtuin deficiency and serves as a potential aging biomarker. Based on our survey of the 3D chromatin architecture, we established connections between sirtuins and potential target genes, thereby informing the development of strategies for aging interventions.

Reorganization of 3D genome architecture provides insights into pathogenesis of early fatty liver disease in laying hens

BackgroundFatty liver disease causes huge economic losses in the poultry industry due to its high occurrence and lethality rate. Three-dimensional (3D) chromatin architecture takes part in disease processing by regulating transcriptional reprogramming. The study is carried out to investigate the alterations of hepatic 3D genome and H3K27ac profiling in early fatty liver (FLS) and reveal their effect on hepatic transcriptional reprogramming in laying hens.ResultsResults show that FLS model is constructed with obvious phenotypes including hepatic visible lipid deposition as well as higher total triglyceride and cholesterol in serum. A/B compartment switching, topologically associating domain (TAD) and chromatin loop changes are identified by high-throughput/resolution chromosome conformation capture (HiC) technology. Targeted genes of these alternations in hepatic 3D genome organization significantly enrich pathways related to lipid metabolism and hepatic damage. H3K27ac differential peaks and differential expression genes (DEGs) identified through RNA-seq analysis are also enriched in these pathways. Notably, certain DEGs are found to correspond with changes in 3D chromatin structure and H3K27ac binding in their promoters. DNA motif analysis reveals that candidate transcription factors are implicated in regulating transcriptional reprogramming. Furthermore, disturbed folate metabolism is observed, as evidenced by lower folate levels and altered enzyme expression.ConclusionOur findings establish a link between transcriptional reprogramming changes and 3D chromatin structure variations during early FLS formation, which provides candidate transcription factors and folate as targets for FLS prevention or treatment.

First person – Mallika Bhattacharya

ABSTRACT First Person is a series of interviews with the first authors of a selection of papers published in Journal of Cell Science, helping researchers promote themselves alongside their papers. Mallika Bhattacharya is first author on ‘ Steroid hormone signaling synchronizes cell migration machinery, adhesion and polarity to direct collective movement’, published in JCS. Mallika conducted the research described in this article while a PhD student in Dr Michelle Starz-Gaiano's lab at University of Maryland, Baltimore, USA. She is now a postdoc in the lab of Dr Elissa Lei at NIDDK, Bethesda, USA. She is excited by basic biology questions in the field of developmental biology and is currently investigating the regulation of 3D chromatin architecture during development.

Haplotype-resolved 3D chromatin architecture of the hybrid pig.

In diploid mammals, allele-specific three-dimensional (3D) genome architecture may lead to imbalanced gene expression. Through ultradeep in situ Hi-C sequencing of three representative somatic tissues (liver, skeletal muscle, and brain) from hybrid pigs generated by reciprocal crosses of phenotypically and physiologically divergent Berkshire and Tibetan pigs, we uncover extensive chromatin reorganization between homologous chromosomes across multiple scales. Haplotype-based interrogation of multi-omic data revealed the tissue dependence of 3D chromatin conformation, suggesting that parent-of-origin-specific conformation may drive gene imprinting. We quantify the effects of genetic variations and histone modifications on allelic differences of long-range promoter-enhancer contacts, which likely contribute to the phenotypic differences between the parental pig breeds. We also observe the fine structure of somatically paired homologous chromosomes in the pig genome, which has a functional implication genome-wide. This work illustrates how allele-specific chromatin architecture facilitates concomitant shifts in allele-biased gene expression, as well as the possible consequential phenotypic changes in mammals.

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-promoterinteractions 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.

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

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

Foxp3 orchestrates reorganization of chromatin architecture to establish regulatory T cell identity

Chromatin conformation reorganization is emerging as an important layer of regulation for gene expression and lineage specification. Yet, how lineage-specific transcription factors contribute to the establishment of cell type-specific 3D chromatin architecture in the immune cells remains unclear, especially for the late stages of T cell subset differentiation and maturation. Regulatory T cells (Treg) are mainly generated in the thymus as a subpopulation of T cells specializing in suppressing excessive immune responses. Here, by comprehensively mapping 3D chromatin organization during Treg cell differentiation, we show that Treg-specific chromatin structures were progressively established during its lineage specification, and highly associated with Treg signature gene expression. Additionally, the binding sites of Foxp3, a Treg lineage specifying transcription factor, were highly enriched at Treg-specific chromatin loop anchors. Further comparison of the chromatin interactions between wide-type Tregs versus Treg cells from Foxp3 knock-in/knockout or newly-generated Foxp3 domain-swap mutant mouse revealed that Foxp3 was essential for the establishment of Treg-specific 3D chromatin architecture, although it was not dependent on the formation of the Foxp3 domain-swapped dimer. These results highlighted an underappreciated role of Foxp3 in modulating Treg-specific 3D chromatin structure formation.

Three-dimensional chromatin landscapes in hepatocellular carcinoma associated with hepatitis B virus.

Three-dimensional (3D) chromatin architecture frequently altered in cancer. However, its changes during the pathogenesis of hepatocellular carcinoma (HCC) remained elusive. Hi-C and RNA-seq were applied to study the 3D chromatin landscapes and gene expression of HCC and ANHT. Hi-C Pro was used to generate genome-wide raw interaction matrices, which were normalized via iterative correction (ICE). Moreover, the chromosomes were divided into different compartments according to the first principal component (E1). Furthermore, topologically associated domains (TADs) were visualized via WashU Epigenome Browser. Furthermore, differential expression analysis of ANHT and HCC was performed using the DESeq2 R package. Additionally, dysregulated genes associated with 3D genome architecture altered were confirmed using TCGA, qRT-PCR, immunohistochemistry (IHC), etc. RESULTS: First, the intrachromosomal interactions of chr1, chr2, chr5, and chr11 were significantly different, and the interchromosomal interactions of chr4-chr10, chr13-chr21, chr15-chr22, and chr16-chr19 are remarkably different between ANHT and HCC, which resulted in the up-regulation of TP53I3 and ZNF738 and the down-regulation of APOC3 and APOA5 in HCC. Second, 49 compartment regions on 18 chromosomes have significantly switched (A-B or B-A) during HCC tumorigenesis, contributing to up-regulation of RAP2A. Finally, a tumor-specific TAD boundary located on chr5: 6271000-6478000 and enhancer hijacking were identified in HCC tissues, potentially associated with the elevated expression of MED10, whose expression were associated with poor prognosis of HCC patients. This study demonstrates the crucial role of chromosomal structure variation in HCC oncogenesis and potential novel biomarkers of HCC, laying a foundation for cancer precision medicine development.

Evolutionary insights into 3D genome organization and epigenetic landscape of Vigna mungo.

Eukaryotic genomes show an intricate three-dimensional (3D) organization within the nucleus that regulates multiple biological processes including gene expression. Contrary to animals, understanding of 3D genome organization in plants remains at a nascent stage. Here, we investigate the evolution of 3D chromatin architecture in legumes. By using cutting-edge PacBio, Illumina, and Hi-C contact reads, we report a gap-free, chromosome-scale reference genome assembly of Vigna mungo, an important minor legume cultivated in Southeast Asia. We spatially resolved V. mungo chromosomes into euchromatic, transcriptionally active A compartment and heterochromatic, transcriptionally-dormant B compartment. We report the presence of TAD-like-regions throughout the diagonal of the HiC matrix that resembled transcriptional quiescent centers based on their genomic and epigenomic features. We observed high syntenic breakpoints but also high coverage of syntenic sequences and conserved blocks in boundary regions than in the TAD-like region domains. Our findings present unprecedented evolutionary insights into spatial 3D genome organization and epigenetic patterns and their interaction within the V. mungo genome. This will aid future genomics and epigenomics research and breeding programs of V. mungo.

Involvement of CCCTC-binding factor in epigenetic regulation of cancer.

A major global health burden continues to be borne by the complex and multifaceted disease of cancer. Epigenetic changes, which are essential for the emergence and spread of cancer, have drawn a huge amount of attention recently. The CCCTC-binding factor (CTCF), which takes part in a wide range of cellular processes including genomic imprinting, X chromosome inactivation, 3D chromatin architecture, local modifications of histone, and RNA polymerase II-mediated gene transcription, stands out among the diverse array of epigenetic regulators. CTCF not only functions as an architectural protein but also modulates DNA methylation and histone modifications. Epigenetic regulation of cancer has already been the focus of plenty of studies. Understanding the role of CTCF in the cancer epigenetic landscape may lead to the development of novel targeted therapeutic strategies for cancer. CTCF has already earned its status as a tumor suppressor gene by acting like a homeostatic regulator of genome integrity and function. Moreover, CTCF has a direct effect on many important transcriptional regulators that control the cell cycle, apoptosis, senescence, and differentiation. As we learn more about CTCF-mediated epigenetic modifications and transcriptional regulations, the possibility of utilizing CTCF as a diagnostic marker and therapeutic target for cancer will also increase. Thus, the current review intends to promote personalized and precision-based therapeutics for cancer patients by shedding light on the complex interplay between CTCF and epigenetic processes.

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.