Peak-valley-peak pattern of histone modifications delineates active regulatory elements and their directionality

Data-driven projections of candidate enhancer-activating SNPs in immune regulation

BACKGROUND

Millions of single nucleotide polymorphisms (SNPs) have been identified in humans, but the functionality of almost all SNPs remains unclear. While current research focuses primarily on SNPs altering one amino acid to another one, the majority of SNPs are located in intergenic spaces. Some of these SNPs can be found in candidate cis-regulatory elements (CREs) such as promoters and enhancers, potentially destroying or creating DNA-binding motifs for transcription factors (TFs) and, hence, deregulating the expression of nearby genes. These aspects are understudied due to the sheer number of SNPs and TF binding motifs, making it challenging to identify SNPs that yield phenotypic changes or altered gene expression.

RESULTS

We developed a data-driven computational protocol to prioritize high-potential SNPs informed from former knowledge for experimental validation. We evaluated the protocol by investigating SNPs in CREs in the Janus kinase (JAK) - Signal Transducer and Activator of Transcription (-STAT) signaling pathway, which is activated by a plethora of cytokines and crucial in controlling immune responses and has been implicated in diseases like cancer, autoimmune disorders, and responses to viral infections. The protocol involves scanning the entire human genome (hg38) to pinpoint DNA sequences that deviate by only one nucleotide from the canonical binding sites (TTCnnnGAA) for STAT TFs. We narrowed down from an initial pool of 3,301,512 SNPs across 17,039,967 nearly complete STAT motifs and identified six potential gain-of-function SNPs in regions likely to influence regulation within the JAK-STAT pathway. This selection was guided by publicly available open chromatin and gene expression data and further refined by filtering for proximity to immune response genes and conservation between the mouse and human genomes.

CONCLUSION

Our findings highlight the value of combining genomic, epigenomic, and cross-species conservation data to effectively narrow down millions of SNPs to a smaller number with a high potential to induce interferon regulation of nearby genes. These SNPs can finally be reviewed manually, laying the groundwork for a more focused and efficient exploration of regulatory SNPs in an experimental setting.

Regulation of transcription patterns, poly(ADP-ribose), and RNA-DNA hybrids by the ATM protein kinase

The ataxia telangiectasia mutated (ATM) protein kinase is a master regulator of the DNA damage response and also an important sensor of oxidative stress. Analysis of gene expression in ataxia-telangiectasia (A-T) patient brain tissue shows that large-scale transcriptional changes occur in patient cerebellum that correlate with the expression level and guanine-cytosine (GC) content of transcribed genes. In human neuron-like cells in culture, we map locations of poly(ADP-ribose) and RNA-DNA hybrid accumulation genome-wide with ATM inhibition and find that these marks also coincide with high transcription levels, active transcription histone marks, and high GC content. Antioxidant treatment reverses the accumulation of R-loops in transcribed regions, consistent with the central role of reactive oxygen species in promoting these lesions. Based on these results, we postulate that transcription-associated lesions accumulate in ATM-deficient cells and that the single-strand breaks and PARylation at these sites ultimately generate changes in transcription that compromise cerebellum function and lead to neurodegeneration over time in A-T patients.

Regulation of transcription patterns, poly-ADP-ribose, and RNA-DNA hybrids by the ATM protein kinase

The ATM protein kinase is a master regulator of the DNA damage response and also an important sensor of oxidative stress. Analysis of gene expression in Ataxia-telangiectasia patient brain tissue shows that large-scale transcriptional changes occur in patient cerebellum that correlate with expression level and GC content of transcribed genes. In human neuron-like cells in culture we map locations of poly-ADP-ribose and RNA-DNA hybrid accumulation genome-wide with ATM inhibition and find that these marks also coincide with high transcription levels, active transcription histone marks, and high GC content. Antioxidant treatment reverses the accumulation of R-loops in transcribed regions, consistent with the central role of ROS in promoting these lesions. Based on these results we postulate that transcription-associated lesions accumulate in ATM-deficient cells and that the single-strand breaks and PARylation at these sites ultimately generate changes in transcription that compromise cerebellum function and lead to neurodegeneration over time in A-T patients.

Epigenetic signature of human immune aging in the GESTALT study

Age-associated DNA methylation in blood cells convey information on health status. However, the mechanisms that drive these changes in circulating cells and their relationships to gene regulation are unknown. We identified age-associated DNA methylation sites in six purified blood-borne immune cell types (naive B, naive CD4+ and CD8+ T cells, granulocytes, monocytes, and NK cells) collected from healthy individuals interspersed over a wide age range. Of the thousands of age-associated sites, only 350 sites were differentially methylated in the same direction in all cell types and validated in an independent longitudinal cohort. Genes close to age-associated hypomethylated sites were enriched for collagen biosynthesis and complement cascade pathways, while genes close to hypermethylated sites mapped to neuronal pathways. In silico analyses showed that in most cell types, the age-associated hypo- and hypermethylated sites were enriched for ARNT (HIF1β) and REST transcription factor (TF) motifs, respectively, which are both master regulators of hypoxia response. To conclude, despite spatial heterogeneity, there is a commonality in the putative regulatory role with respect to TF motifs and histone modifications at and around these sites. These features suggest that DNA methylation changes in healthy aging may be adaptive responses to fluctuations of oxygen availability.

Methyl-CpG binding domain 2 (Mbd2) is an epigenetic regulator of autism-risk genes and cognition

The Methyl-CpG-Binding Domain Protein family has been implicated in neurodevelopmental disorders. The Methyl-CpG-binding domain 2 (Mbd2) binds methylated DNA and was shown to play an important role in cancer and immunity. Some evidence linked this protein to neurodevelopment. However, its exact role in neurodevelopment and brain function is mostly unknown. Here we show that Mbd2-deficiency in mice (Mbd2-/-) results in deficits in cognitive, social and emotional functions. Mbd2 binds regulatory DNA regions of neuronal genes in the hippocampus and loss of Mbd2 alters the expression of hundreds of genes with a robust down-regulation of neuronal gene pathways. Further, a genome-wide DNA methylation analysis found an altered DNA methylation pattern in regulatory DNA regions of neuronal genes in Mbd2-/- mice. Differentially expressed genes significantly overlap with gene-expression changes observed in brains of Autism Spectrum Disorder (ASD) individuals. Notably, downregulated genes are significantly enriched for human ortholog ASD risk genes. Observed hippocampal morphological abnormalities were similar to those found in individuals with ASD and ASD rodent models. Hippocampal Mbd2 knockdown partially recapitulates the behavioral phenotypes observed in Mbd2-/- mice. These findings suggest that Mbd2 is a novel epigenetic regulator of genes that are associated with ASD in humans. Mbd2 loss causes behavioral alterations that resemble those found in ASD individuals.

TF-Prioritizer: a Java pipeline to prioritize condition-specific transcription factors

BACKGROUND

Eukaryotic gene expression is controlled by cis-regulatory elements (CREs), including promoters and enhancers, which are bound by transcription factors (TFs). Differential expression of TFs and their binding affinity at putative CREs determine tissue- and developmental-specific transcriptional activity. Consolidating genomic datasets can offer further insights into the accessibility of CREs, TF activity, and, thus, gene regulation. However, the integration and analysis of multimodal datasets are hampered by considerable technical challenges. While methods for highlighting differential TF activity from combined chromatin state data (e.g., chromatin immunoprecipitation [ChIP], ATAC, or DNase sequencing) and RNA sequencing data exist, they do not offer convenient usability, have limited support for large-scale data processing, and provide only minimal functionality for visually interpreting results.

RESULTS

We developed TF-Prioritizer, an automated pipeline that prioritizes condition-specific TFs from multimodal data and generates an interactive web report. We demonstrated its potential by identifying known TFs along with their target genes, as well as previously unreported TFs active in lactating mouse mammary glands. Additionally, we studied a variety of ENCODE datasets for cell lines K562 and MCF-7, including 12 histone modification ChIP sequencing as well as ATAC and DNase sequencing datasets, where we observe and discuss assay-specific differences.

CONCLUSION

TF-Prioritizer accepts ATAC, DNase, or ChIP sequencing and RNA sequencing data as input and identifies TFs with differential activity, thus offering an understanding of genome-wide gene regulation, potential pathogenesis, and therapeutic targets in biomedical research.

Differential regulation of alternate promoter regions in Sox17 during endodermal and vascular endothelial development

Sox17 gene expression is essential for both endothelial and endodermal cell differentiation. To better understand the genetic basis for the expression of multiple Sox17 mRNA forms, we identified and performed CRISPR/Cas9 mutagenesis of two evolutionarily conserved promoter regions (CRs). The deletion of the upstream and endothelial cell-specific CR1 caused only a modest increase in lympho-vasculogenesis likely via reduced Notch signaling downstream of SOX17. In contrast, the deletion of the downstream CR2 region, which functions in both endothelial and endodermal cells, impairs both vascular and endodermal development causing death by embryonic day 12.5. Analyses of 3D chromatin looping, transcription factor binding, histone modification, and chromatin accessibility data at the Sox17 locus and surrounding region further support differential regulation of the two promoters during the development.

Identification of putative enhancer-like elements predicts regulatory networks active in planarian adult stem cells

Planarians have become an established model system to study regeneration and stem cells, but the regulatory elements in the genome remain almost entirely undescribed. Here, by integrating epigenetic and expression data we use multiple sources of evidence to predict enhancer elements active in the adult stem cell populations that drive regeneration. We have used ChIP-seq data to identify genomic regions with histone modifications consistent with enhancer activity, and ATAC-seq data to identify accessible chromatin. Overlapping these signals allowed for the identification of a set of high-confidence candidate enhancers predicted to be active in planarian adult stem cells. These enhancers are enriched for predicted transcription factor (TF) binding sites for TFs and TF families expressed in planarian adult stem cells. Footprinting analyses provided further evidence that these potential TF binding sites are likely to be occupied in adult stem cells. We integrated these analyses to build testable hypotheses for the regulatory function of TFs in stem cells, both with respect to how pluripotency might be regulated, and to how lineage differentiation programs are controlled. We found that our predicted GRNs were independently supported by existing TF RNAi/RNA-seq datasets, providing further evidence that our work predicts active enhancers that regulate adult stem cells and regenerative mechanisms.

Myod1 and GR coordinate myofiber-specific transcriptional enhancers

Skeletal muscle is a dynamic tissue the size of which can be remodeled through the concerted actions of various cues. Here, we investigated the skeletal muscle transcriptional program and identified key tissue-specific regulatory genetic elements. Our results show that Myod1 is bound to numerous skeletal muscle enhancers in collaboration with the glucocorticoid receptor (GR) to control gene expression. Remarkably, transcriptional activation controlled by these factors occurs through direct contacts with the promoter region of target genes, via the CpG-bound transcription factor Nrf1, and the formation of Ctcf-anchored chromatin loops, in a myofiber-specific manner. Moreover, we demonstrate that GR negatively controls muscle mass and strength in mice by down-regulating anabolic pathways. Taken together, our data establish Myod1, GR and Nrf1 as key players of muscle-specific enhancer-promoter communication that orchestrate myofiber size regulation.

Novel ChIP-seq simulating program with superior versatility: isChIP

Chromatin immunoprecipitation followed by next-generation sequencing (ChIP-seq) is recognized as an extremely powerful tool to study the interaction of numerous transcription factors and other chromatin-associated proteins with DNA. The core problem in the optimization of ChIP-seq protocol and the following computational data analysis is that a 'true' pattern of binding events for a given protein factor is unknown. Computer simulation of the ChIP-seq process based on 'a-priory known binding template' can contribute to a drastically reduce the number of wet lab experiments and finally help achieve radical optimization of the entire processing pipeline. We present a newly developed ChIP-sequencing simulation algorithm implemented in the novel software, in silico ChIP-seq (isChIP). We demonstrate that isChIP closely approximates real ChIP-seq protocols and is able to model data similar to those obtained from experimental sequencing. We validated isChIP using publicly available datasets generated for well-characterized transcription factors Oct4 and Sox2. Although the novel software is compatible with the Illumina protocols by default, it can also successfully perform simulations with a number of alternative sequencing platforms such as Roche454, Ion Torrent and SOLiD as well as model ChIP -Exo. The versatility of isChIP was demonstrated through modelling a wide range of binding events, including those of transcription factors and chromatin modifiers. We also performed a comparative analysis against a few existing ChIP-seq simulators and showed the fundamental superiority of our model. Due to its ability to utilize known binding templates, isChIP can potentially be employed to help investigators choose the most appropriate analytical software through benchmarking of available ChIP-seq programs and optimize the experimental parameters of ChIP-seq protocol. isChIP software is freely available at https://github.com/fnaumenko/isChIP.

Characterization of histone modification patterns and prediction of novel promoters using functional principal component analysis

Characterization of distinct histone methylation and acetylation binding patterns in promoters and prediction of novel regulatory regions remains an important area of genomic research, as it is hypothesized that distinct chromatin signatures may specify unique genomic functions. However, methods that have been proposed in the literature are either descriptive in nature or are fully parametric and hence more restrictive in pattern discovery. In this article, we propose a two-step non-parametric statistical inference procedure to characterize unique histone modification patterns and apply it to analyzing the binding patterns of four histone marks, H3K4me2, H3K4me3, H3K9ac, and H4K20me1, in human B-lymphoblastoid cells. In the first step, we used a functional principal component analysis method to represent the concatenated binding patterns of these four histone marks around the transcription start sites as smooth curves. In the second step, we clustered these curves to reveal several unique classes of binding patterns. These uncovered patterns were used in turn to scan the whole-genome to predict novel and alternative promoters. Our analyses show that there are three distinct promoter binding patterns of active genes. Further, 19654 regions not within known gene promoters were found to overlap with human ESTs, CpG islands, or common SNPs, indicative of their potential role in gene regulation, including being potential novel promoter regions.

Ancestrally Duplicated Conserved Noncoding Element Suggests Dual Regulatory Roles of HOTAIR in cis and trans

HOTAIR was proposed to regulate either HoxD cluster genes in trans or HoxC cluster genes in cis, a mechanism that remains unclear. We have identified a 32-nucleotide conserved noncoding element (CNE) as HOTAIR ancient sequence that likely originated at the root of vertebrate. The second round of whole-genome duplication resulted in one copy of the CNE within HOTAIR and another copy embedded in noncoding transcript of HOXD11. Paralogous CNEs underwent compensatory mutations, exhibit sequence complementarity with respect to transcripts directionality, and have high affinity in vitro. The HOTAIR CNE resembled a poised enhancer in stem cells and an active enhancer in HOTAIR-expressing cells. HOTAIR expression is positively correlated with HOXC11 in cis and negatively correlated with HOXD11 in trans. We propose a dual modality of HOTAIR regulation where transcription of HOTAIR and its embedded enhancer regulates HOXC11 in cis and sequence complementarity between paralogous CNEs suggests HOXD11 regulation in trans.

Specific subfamilies of transposable elements contribute to different domains of T lymphocyte enhancers

Transposable elements (TEs) compose nearly half of mammalian genomes and provide building blocks for cis-regulatory elements. Using high-throughput sequencing, we show that 84 TE subfamilies are overrepresented, and distributed in a lineage-specific fashion in core and boundary domains of CD8⁺ T cell enhancers. Endogenous retroviruses are most significantly enriched in core domains with accessible chromatin, and bear recognition motifs for immune-related transcription factors. In contrast, short interspersed elements (SINEs) are preferentially overrepresented in nucleosome-containing boundaries. A substantial proportion of these SINEs harbor a high density of the enhancer-specific histone mark H3K4me1 and carry sequences that match enhancer boundary nucleotide composition. Motifs with regulatory features are better preserved within enhancer-enriched TE copies compared to their subfamily equivalents located in gene deserts. TE-rich and TE-poor enhancers associate with both shared and unique gene groups and are enriched in overlapping functions related to lymphocyte and leukocyte biology. The majority of T cell enhancers are shared with other immune lineages and are accessible in common hematopoietic progenitors. A higher proportion of immune tissue-specific enhancers are TE-rich compared to enhancers specific to other tissues, correlating with higher TE occurrence in immune gene-associated genomic regions. Our results suggest that during evolution, TEs abundant in these regions and carrying motifs potentially beneficial for enhancer architecture and immune functions were particularly frequently incorporated by evolving enhancers. Their putative selection and regulatory cooption may have accelerated the evolution of immune regulatory networks.

Enhancer and Transcription Factor Dynamics during Myeloid Differentiation Reveal an Early Differentiation Block in Cebpa null Progenitors

Transcription factors PU.1 and CEBPA are required for the proper coordination of enhancer activity during granulocytic-monocytic (GM) lineage differentiation to form myeloid cells. However, precisely how these factors control the chronology of enhancer establishment during differentiation is not known. Through integrated analyses of enhancer dynamics, transcription factor binding, and proximal gene expression during successive stages of murine GM-lineage differentiation, we unravel the distinct kinetics by which PU.1 and CEBPA coordinate GM enhancer activity. We find no evidence of a pioneering function of PU.1 during late GM-lineage differentiation. Instead, we delineate a set of enhancers that gain accessibility in a CEBPA-dependent manner, suggesting a pioneering function of CEBPA. Analyses of Cebpa null bone marrow demonstrate that CEBPA controls PU.1 levels and, unexpectedly, that the loss of CEBPA results in an early differentiation block. Taken together, our data provide insights into how PU.1 and CEBPA functionally interact to drive GM-lineage differentiation.

EpiSAFARI: sensitive detection of valleys in epigenetic signals for enhancing annotations of functional elements

MOTIVATION

Functional genomics experiments generate genomewide signal profiles that are dense information sources for annotating the regulatory elements. These profiles measure epigenetic activity at the nucleotide resolution and they exhibit distinctive patterns as they fluctuate along the genome. Most notable of these patterns are the valley patterns that are prevalently observed in assays such as ChIP Sequencing and bisulfite sequencing. The genomic positions of valleys pinpoint locations of cis-regulatory elements such as enhancers and insulators. Systematic identification of the valleys provides novel information for delineating the annotation of regulatory elements. Nevertheless, the valleys are not reported by majority of the analysis pipelines.

RESULTS

We describe EpiSAFARI, a computational method for sensitive detection of valleys from diverse types of epigenetic profiles. EpiSAFARI employs a novel smoothing method for decreasing noise in signal profiles and accounts for technical factors such as sparse signals, mappability and nucleotide content. In performance comparisons, EpiSAFARI performs favorably in terms of accuracy. The histone modification valleys detected by EpiSAFARI exhibit high conservation, transcription factor binding and they are enriched in nascent transcription. In addition, the large clusters of histone valleys are found to be enriched at the promoters of the developmentally associated genes. Differential histone valleys exhibit concordance with differential DNase signal at cell line specific valleys. DNA methylation valleys exhibit elevated conservation and high transcription factor binding. Specifically, we observed enriched binding of transcription factors associated with chromatin structure around methyl-valleys.

AVAILABILITY AND IMPLEMENTATION

EpiSAFARI is publicly available at https://github.com/harmancilab/EpiSAFARI.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The DNA methylation landscape of enhancers in the guinea pig hippocampus

AIM

To determine the state of methylation of DNA molecules in the guinea pig hippocampus that are associated with either poised or active enhancers.

METHODS

We used sequential chromatin immunoprecipitation-bisulfite-sequencing with an antibody to H3K4me1 to map the state of methylation of DNA that is found within enhancers. Actively transcribing transcription start sites were mapped by chromatin immunoprecipitation-sequencing with an antibody to RNApolII-PS5. Total DNA methylation was mapped using reduced representation bisulfite sequencing.

RESULTS

DNA that overlaps with H3K4me1 binding regions in the genome is heavily methylated. However, DNA molecules that are found in H3K4me1 chromatin are hypomethylated, while DNA found in enhancers that are associated with active transcription is further demethylated. Differential methylation in enhancers is spotted in single CGs, bimodal and corresponds to transcription factor binding sites.

CONCLUSION

Our study delineates the DNA methylation status of H3K4 me1-bound regions in the hippocampus in active and inactive genes.

Class 3 semaphorins are transcriptionally regulated by 1,25(OH)2D3 in osteoblasts

The vitamin D endocrine system is essential for calcium metabolism and skeletal integrity. 1,25-dihydroxyvitamin D₃ [1,25(OH)₂D₃] regulates bone mineral homeostasis and acts directly on osteoblasts. In the present study we characterized the transcriptional regulation of the class 3 semaphorin (Sema3) gene family by 1,25(OH)₂D₃ in osteoblastic cells. Class 3 semaphorins are secreted proteins that regulate cell growth, morphology and migration, and were recently shown to be involved in bone homeostasis. In ST2, MC3T3-E1 and primary calvarial osteoblast cell cultures we found that all members of the Sema3 gene family were expressed, and that Sema3e and Sema3f were the most strongly induced 1,25(OH)₂D₃ target genes among the studied cell types. In addition, transcription of Sema3b and Sema3c was upregulated, whereas Sema3d and Sema3g was downregulated by 1,25(OH)₂D₃ in different osteoblastic cells. Chromatin immunoprecipitation analysis linked to DNA sequencing (ChIP-seq analysis) revealed the presence of the vitamin D receptor at multiple genomic loci in the proximity of Sema3 genes, demonstrating that the genes are primary 1,25(OH)₂D₃ targets. Furthermore, we showed that recombinant SEMA3E and SEMA3F protein were able to inhibit osteoblast proliferation. However, recombinant SEMA3s did not affect ST2 cell migration. The expression of class 3 semaphorins in osteoblasts together with their regulation by 1,25(OH)₂D₃ suggests that these genes, involved in the regulation of bone homeostasis, are additional mediators for 1,25(OH)₂D₃ signaling in osteoblasts.

CIPHER: a flexible and extensive workflow platform for integrative next-generation sequencing data analysis and genomic regulatory element prediction

BACKGROUND

Next-generation sequencing (NGS) approaches are commonly used to identify key regulatory networks that drive transcriptional programs. Although these technologies are frequently used in biological studies, NGS data analysis remains a challenging, time-consuming, and often irreproducible process. Therefore, there is a need for a comprehensive and flexible workflow platform that can accelerate data processing and analysis so more time can be spent on functional studies.

RESULTS

We have developed an integrative, stand-alone workflow platform, named CIPHER, for the systematic analysis of several commonly used NGS datasets including ChIP-seq, RNA-seq, MNase-seq, DNase-seq, GRO-seq, and ATAC-seq data. CIPHER implements various open source software packages, in-house scripts, and Docker containers to analyze and process single-ended and pair-ended datasets. CIPHER's pipelines conduct extensive quality and contamination control checks, as well as comprehensive downstream analysis. A typical CIPHER workflow includes: (1) raw sequence evaluation, (2) read trimming and adapter removal, (3) read mapping and quality filtering, (4) visualization track generation, and (5) extensive quality control assessment. Furthermore, CIPHER conducts downstream analysis such as: narrow and broad peak calling, peak annotation, and motif identification for ChIP-seq, differential gene expression analysis for RNA-seq, nucleosome positioning for MNase-seq, DNase hypersensitive site mapping, site annotation and motif identification for DNase-seq, analysis of nascent transcription from Global-Run On (GRO-seq) data, and characterization of chromatin accessibility from ATAC-seq datasets. In addition, CIPHER contains an "analysis" mode that completes complex bioinformatics tasks such as enhancer discovery and provides functions to integrate various datasets together.

CONCLUSIONS

Using public and simulated data, we demonstrate that CIPHER is an efficient and comprehensive workflow platform that can analyze several NGS datasets commonly used in genome biology studies. Additionally, CIPHER's integrative "analysis" mode allows researchers to elicit important biological information from the combined dataset analysis.

Initial high-resolution microscopic mapping of active and inactive regulatory sequences proves non-random 3D arrangements in chromatin domain clusters

BACKGROUND

The association of active transcription regulatory elements (TREs) with DNAse I hypersensitivity (DHS[+]) and an 'open' local chromatin configuration has long been known. However, the 3D topography of TREs within the nuclear landscape of individual cells in relation to their active or inactive status has remained elusive. Here, we explored the 3D nuclear topography of active and inactive TREs in the context of a recently proposed model for a functionally defined nuclear architecture, where an active and an inactive nuclear compartment (ANC-INC) form two spatially co-aligned and functionally interacting networks.

RESULTS

Using 3D structured illumination microscopy, we performed 3D FISH with differently labeled DNA probe sets targeting either sites with DHS[+], apparently active TREs, or DHS[-] sites harboring inactive TREs. Using an in-house image analysis tool, DNA targets were quantitatively mapped on chromatin compaction shaped 3D nuclear landscapes. Our analyses present evidence for a radial 3D organization of chromatin domain clusters (CDCs) with layers of increasing chromatin compaction from the periphery to the CDC core. Segments harboring active TREs are significantly enriched at the decondensed periphery of CDCs with loops penetrating into interchromatin compartment channels, constituting the ANC. In contrast, segments lacking active TREs (DHS[-]) are enriched toward the compacted interior of CDCs (INC).

CONCLUSIONS

Our results add further evidence in support of the ANC-INC network model. The different 3D topographies of DHS[+] and DHS[-] sites suggest positional changes of TREs between the ANC and INC depending on their functional state, which might provide additional protection against an inappropriate activation. Our finding of a structural organization of CDCs based on radially arranged layers of different chromatin compaction levels indicates a complex higher-order chromatin organization beyond a dichotomic classification of chromatin into an 'open,' active and 'closed,' inactive state.