The insulator binding protein CTCF positions 20 nucleosomes around its binding sites across the human genome

Nucleosome dynamics render heterochromatin accessible in living human cells

The eukaryotic genome is packaged into chromatin, which is composed of a nucleosomal filament that coils up to form more compact structures. Chromatin exists in two main forms: euchromatin, which is relatively decondensed and enriched in transcriptionally active genes, and heterochromatin, which is condensed and transcriptionally repressed. It is widely accepted that chromatin architecture modulates DNA accessibility, restricting the access of sequence-specific, gene-regulatory, transcription factors to the genome. However, the evidence for this model derives primarily from experiments with isolated nuclei, in which chromatin remodeling has ceased, resulting in a static chromatin structure. Here, using a DNA methyltransferase to measure accessibility in vivo, we show that both euchromatin and heterochromatin are fully accessible in living human cells, whereas centromeric α-satellite chromatin is partly inaccessible. We conclude that all nucleosomes in euchromatin and heterochromatin are highly dynamic in living cells, except for nucleosomes in centromeric chromatin.

High-resolution CTCF footprinting reveals impact of chromatin state on cohesin extrusion

Cohesin-mediated DNA loop extrusion enables gene regulation by distal enhancers through the establishment of chromosome structure and long-range enhancer-promoter interactions. The best characterized cohesin-related structures, such as topologically associating domains (TADs) anchored at convergent CTCF binding sites, represent static conformations. Consequently, loop extrusion dynamics remain poorly understood. To better characterize static and dynamically extruding chromatin loop structures, we use MNase-based 3D genome assays to simultaneously determine CTCF and cohesin localization as well as the 3D contacts they mediate. Here we present CTCF Analyzer (with) Multinomial Estimation (CAMEL), a tool that identifies CTCF footprints at near base-pair resolution in CTCF MNase HiChiP. We also use Region Capture Micro-C to identify a CTCF-adjacent footprint that is attributed to cohesin occupancy. We leverage this substantial advance in resolution to determine that the fully extruded (CTCF-CTCF loop) state is rare genome-wide with locus-specific variation from ~1-10%. We further investigate the impact of chromatin state on loop extrusion dynamics and find that active regulatory elements impede cohesin extrusion. These findings support a model of topological regulation whereby the transient, partially extruded state facilitates enhancer-promoter contacts that can regulate transcription.

BINDER achieves accurate identification of hierarchical TADs by comprehensively characterizing consensus TAD boundaries

As crucial chromatin structures, hierarchical TADs play important roles in epigenetic organization, transcriptional activity, gene regulation, and cell differentiation. Currently, it remains a highly challenging task to accurately identify hierarchical TADs in a computational manner. The key bottleneck for existing TAD callers lies in the difficulty in the prediction of precise TAD boundaries. We solve this problem by introducing a novel algorithm, called BINDER, which conducts a boundary consensus approach, and then precisely locate hierarchical TAD boundaries by developing a multifaceted boundary characterization strategy. In comparison with other leading TAD callers, BINDER shows significant improvement in identifying hierarchical TADs and exhibits the strongest robustness with ultrasparse data, which supports the importance of boundary identification in calling hierarchical TADs. Applying BINDER to experimental data and mouse hematopoietic cases, we find that the hierarchical TADs identified by BINDER show strong biological relevance in their epigenetic organization, transcriptional activity, DNA motifs, and coregulation during cellular differentiation. BINDER discovers differences in the enrichment of two specific transcription factors, CHD1 and CHD2, at TAD boundaries with different hierarchies. It also observes variations in the gene expression of TADs with different hierarchies during cellular differentiation.

Beyond genomic weaving: molecular roles for CTCF outside cohesin loop extrusion

CCCTC-binding factor (CTCF) is a key regulator of 3D genome organization and transcriptional activity. Beyond its well-characterized role in facilitating cohesin-mediated loop extrusion, CTCF exhibits several cohesin-independent activities relevant to chromatin structure and various nuclear processes. These functions include patterning of nucleosome arrangement and chromatin accessibility through interactions with ATP-dependent chromatin remodelers. In addition to influencing transcription, DNA replication, and DNA repair in ways that are separable from its role in loop extrusion, CTCF also interacts with RNA and contributes to RNA splicing and condensation of transcriptional activators. Here, we review recent insight into cohesin-independent activities of CTCF, highlighting its multifaceted roles in chromatin biology and transcriptional regulation.

Interpreting the CTCF-mediated sequence grammar of genome folding with AkitaV2

Interphase mammalian genomes are folded in 3D with complex locus-specific patterns that impact gene regulation. CTCF (CCCTC-binding factor) is a key architectural protein that binds specific DNA sites, halts cohesin-mediated loop extrusion, and enables long-range chromatin interactions. There are hundreds of thousands of annotated CTCF-binding sites in mammalian genomes; disruptions of some result in distinct phenotypes, while others have no visible effect. Despite their importance, the determinants of which CTCF sites are necessary for genome folding and gene regulation remain unclear. Here, we update and utilize Akita, a convolutional neural network model, to extract the sequence preferences and grammar of CTCF contributing to genome folding. Our analyses of individual CTCF sites reveal four predictions: (i) only a small fraction of genomic sites are impactful; (ii) impact is highly dependent on sequences flanking the core CTCF binding motif; (iii) core and flanking nucleotides contribute largely additively to the overall impact of a site; (iv) sites created as combinations of different core and flanking sequences have impacts proportional to the product of their average impacts, i.e. they are broadly compatible. Our analysis of collections of CTCF sites make two predictions for multi-motif grammar: (i) insulation strength depends on the number of CTCF sites within a cluster, and (ii) pattern formation is governed by the orientation and spacing of these sites, rather than any inherent specialization of the CTCF motifs themselves. In sum, we present a framework for using neural network models to probe the sequences instructing genome folding and provide a number of predictions to guide future experimental inquiries.

Interphase chromatin biophysics and mechanics: new perspectives and open questions

The physical organization and properties of chromatin within the interphase nucleus are intimately linked to a wide range of functional DNA-based processes. In this context, interphase chromatin mechanics - that is, how chromatin, physically, responds to forces - is gaining increasing attention. Recent methodological advances for probing the force-response of chromatin in cellulo open new avenues for research, as well as new questions. This review discusses emerging views from these approaches and others, including recent in vitro single-molecule studies of cohesin and condensin motor activities, providing insights into physical and material aspects of chromatin, its plasticity in the context of functional processes, its nonequilibrium or 'active matter' properties, and the importance of factors such as chromatin fiber tension and stiffness. This growing field offers exciting opportunities to better understand the interplay between interphase chromosome structure, dynamics, mechanics, and functions.

Comprehensive evaluation of the impact of whole-genome bisulfite sequencing (WGBS) on the fragmentomic characteristics of plasma cell-free DNA

BACKGROUND

Cell-free DNA (cfDNA) is non-randomly fragmented in human body fluids. Analyzing such fragmentation patterns of cfDNA holds great promise for liquid biopsy. Whole-genome bisulfite sequencing (WGBS) is widely used for cfDNA methylation profiling. However, its applicability for studying fragmentomic characteristics remains largely unexplored.

METHODS

We performed paired WGBS and whole-genome sequencing (WGS) on 66 peripheral plasma samples from 58 pregnant women. Then, we systematically compared the fragmentation patterns of cell-free nuclear DNA and mitochondrial DNA (mtDNA) sequenced from these two approaches. Additionally, we evaluated the extent of the size shortening in fetal-derived cfDNA and estimated the fetal DNA fraction in maternal plasma using both sequencing methods.

RESULTS

Compared to WGS samples, WGBS samples demonstrated a significantly lower genome coverage and higher GC content in cfDNA. They also showed a significant decrease in the size of cell-free nuclear DNA, along with alterations in the end motif pattern that were specifically associated with CpG and "CC" sites. While there was a slight shift in the inferred nucleosome footprint from cfDNA coverages in WGBS samples, the cfDNA coverage patterns in CTCF and TSS regions remained highly consistent between these two sequencing methods. Both methods accurately reflected gene expression levels through their TSS coverages. Additionally, WGBS samples exhibited an increased abundance and longer length of mtDNA in plasma. Furthermore, we observed the size shortening of fetal cfDNA in plasma consistently, with a highly correlated fetal DNA fraction inferred by cfDNA coverage between WGBS and WGS samples (r = 0.996). However, the estimated fetal cfDNA fraction in WGBS samples was approximately 7 % lower than in WGS samples.

CONCLUSIONS

We confirmed that WGBS can introduce artificial breakages to cfDNA, leading to altered fragmentomic patterns in both nuclear and mitochondrial DNA. However, WGBS cfDNA remains suitable for analyzing certain cfDNA fragmentomic characteristics, such as coverage in genome regulation regions and the essential characteristics of fetal DNA in maternal plasma.

KAS-ATAC reveals the genome-wide single-stranded accessible chromatin landscape of the human genome

Gene regulation in most eukaryotes involves two fundamental processes: alterations in genome packaging by nucleosomes, with active cis-regulatory elements (CREs) generally characterized by open-chromatin configuration, and transcriptional activation. Mapping these physical properties and biochemical activities, through profiling chromatin accessibility and active transcription, is a key tool for understanding the logic and mechanisms of transcription and its regulation. However, the relationship between these two states has not been accessible to simultaneous measurement. To this end, we developed KAS-ATAC, a combination of the kethoxal-assisted ssDNA sequencing (KAS-seq) and assay for transposase-accessible chromatin using sequencing (ATAC-seq) methods for mapping single-stranded DNA (and thus active transcription) and chromatin accessibility, respectively, enabling the genome-wide identification of DNA fragments that are simultaneously accessible and contain ssDNA. We use KAS-ATAC to evaluate levels of active transcription over different CRE classes, to estimate absolute levels of transcribed accessible DNA over CREs, to map nucleosomal configurations associated with RNA polymerase activities, and to assess transcription factor association with transcribed DNA through transcription factor binding site (TFBS) footprinting. We observe lower levels of transcription over distal enhancers compared with promoters and distinct nucleosomal configurations around transcription initiation sites associated with active transcription. We find that most TFs associate equally with transcribed and nontranscribed DNA, but a few factors specifically do not exhibit footprints over ssDNA-containing fragments. We anticipate KAS-ATAC to continue to derive useful insights into chromatin organization and transcriptional regulation in other contexts in the future.

A genome wide code to define cell-type specific CTCF binding and chromatin organization

CTCF-mediated chromatin folding plays a key role in gene regulation, however the mechanisms underlying cell type-specific control are not fully elucidated. Comprehensive analyses reveal that CTCF binding stability and cohesin overlap in mice and humans is regulated by species specific differences in accessibility, the presence of CTCF extended Up and Downstream binding site sequences and motifs corresponding to expressed TFs enriched at most bound sites. CTCF repositions nucleosomes asymmetrically, with linker spacing altered by accessibility, while cohesin-mediated nucleosome phasing is affected by surrounding motifs that similarly impact chromatin insulation. Importantly, cell type-specific transcriptional programs determine motif enrichment at CTCF bound sites, reflecting the stabilizing/destabilizing effect of individual TFs. These studies identify mechanisms underlying cell type-specific CTCF profiles, linked to local and long-range chromatin organization.

Neuronal Progenitors Suffer Genotoxic Stress in the Drosophila Clock Mutant per0

The physiological role and the molecular architecture of the circadian clock in fully developed organisms are well established. Yet, we have a limited understanding of the function of the clock during ontogenesis. We have used a null mutant (per0) of the clock gene period (per) in Drosophila melanogaster to ask whether PER may play a role during normal brain development. In third-instar larvae, we have observed that the absence of functional per results in increased genotoxic stress compared to wild-type controls. We have detected increased double-strand DNA breaks in the central nervous system and chromosome aberrations in dividing neuronal precursor cells. We have demonstrated that reactive oxygen species (ROS) are causal to the genotoxic effect and that expression of PER in glia is necessary and sufficient to suppress such a phenotype. Finally, we have shown that the absence of PER may result in less condensed chromatin, which contributes to DNA damage.

Dissecting the cis-regulatory syntax of transcription initiation with deep learning

Despite extensive characterization of mammalian Pol II transcription, the DNA sequence determinants of transcription initiation at a third of human promoters and most enhancers remain poorly understood. We trained and interpreted a neural network called ProCapNet that accurately models base-resolution initiation profiles from PRO-cap experiments using local DNA sequence. ProCapNet learns sequence motifs with distinct effects on initiation rates and TSS positioning and uncovers context-specific cryptic initiator elements intertwined within other TF motifs. ProCapNet annotates predictive motifs in nearly all actively transcribed regulatory elements across multiple cell-lines, revealing a shared cis-regulatory logic across promoters and enhancers and a highly epistatic sequence syntax of cooperative and competitive motif interactions. ProCapNet models of steady-state RAMPAGE profiles distill initiation signals on par with models trained directly on PRO-cap profiles. ProCapNet learns a largely cell-type-agnostic cis-regulatory code of initiation complementing sequence drivers of cell-type-specific chromatin state critical for accurate prediction of cell-type-specific transcription initiation.

GHT-SELEX demonstrates unexpectedly high intrinsic sequence specificity and complex DNA binding of many human transcription factors

A long-standing challenge in human regulatory genomics is that transcription factor (TF) DNA-binding motifs are short and degenerate, while the genome is large. Motif scans therefore produce many false-positive binding site predictions. By surveying 179 TFs across 25 families using >1,500 cyclic in vitro selection experiments with fragmented, naked, and unmodified genomic DNA - a method we term GHT-SELEX (Genomic HT-SELEX) - we find that many human TFs possess much higher sequence specificity than anticipated. Moreover, genomic binding regions from GHT-SELEX are often surprisingly similar to those obtained in vivo (i.e. ChIP-seq peaks). We find that comparable specificity can also be obtained from motif scans, but performance is highly dependent on derivation and use of the motifs, including accounting for multiple local matches in the scans. We also observe alternative engagement of multiple DNA-binding domains within the same protein: long C2H2 zinc finger proteins often utilize modular DNA recognition, engaging different subsets of their DNA binding domain (DBD) arrays to recognize multiple types of distinct target sites, frequently evolving via internal duplication and divergence of one or more DBDs. Thus, contrary to conventional wisdom, it is common for TFs to possess sufficient intrinsic specificity to independently delineate cellular targets.

Targeting the 3D genome by anthracyclines for chemotherapeutic effects

The chromatins are folded into three-dimensional (3D) structures inside cells, which coordinates the regulation of gene transcription by the non-coding regulatory elements. Aberrant chromatin 3D folding has been shown in many diseases, such as acute myeloid leukemia (AML), and may contribute to tumorigenesis. The anthracycline topoisomerase II inhibitors can induce histone eviction and DNA damage. We performed genome-wide high-resolution mapping of the chemotherapeutic effects of various clinically used anthracycline drugs. ATAC-seq was used to profile the histone eviction effects of different anthracyclines. TOP2A ChIP-seq was used to profile the potential DNA damage regions. Integrated analyses show that different anthracyclines have distinct target selectivity on epigenomic regions, based on their respective ATAC-seq and ChIP-seq profiles. We identified the underlying molecular mechanism that unique anthracycline variants selectively target chromatin looping anchors via disrupting CTCF binding, suggesting an additional potential therapeutic effect on the 3D genome. We further performed Hi-C experiments, and data from K562 cells treated with the selective anthracycline drugs indicate that the 3D chromatin organization is disrupted. Furthermore, AML patients receiving anthracycline drugs showed altered chromatin structures around potential looping anchors, which linked to distinct clinical outcomes. Our data indicate that anthracyclines are potent and selective epigenomic targeting drugs and can target the 3D genome for anticancer therapy, which could be used for personalized medicine to treat tumors with aberrant 3D chromatin structures.

The characteristics of CTCF binding sequences contribute to enhancer blocking activity

While the elements encoding enhancers and promoters have been relatively well studied, the full spectrum of insulator elements which bind the CCCTC binding factor (CTCF), is relatively poorly characterized. This is partly due to the genomic context of CTCF sites greatly influencing their roles and activity. Here we have developed an experimental system to determine the ability of minimal, consistently sized, individual CTCF elements to interpose between enhancers and promoters and thereby reduce gene expression during differentiation. Importantly, each element is tested in the identical location thereby minimising the effect of genomic context. We found no correlation between the ability of CTCF elements to block enhancer-promoter activity with the degree of evolutionary conservation; their resemblance to the consensus core sequences; or the number of CTCF core motifs harboured in the element. Nevertheless, we have shown that the strongest enhancer-promoter blockers include a previously described bound element lying upstream of the CTCF core motif. In addition, we found other uncharacterised DNaseI footprints located close to the core motif that may affect function. We have developed an assay of CTCF sequences which will enable researchers to sub-classify individual CTCF elements in a uniform and unbiased way.

A novel approach for in vivo DNA footprinting using short double-stranded cell-free DNA from plasma

Here, we present a method for enrichment of double-stranded cfDNA with an average length of ∼40 bp from cfDNA for high-throughput DNA sequencing. This class of cfDNA is enriched at gene promoters and binding sites of transcription factors or structural DNA-binding proteins, so that a genome-wide DNA footprint is directly captured from liquid biopsies. In short double-stranded cfDNA from healthy individuals, we find significant enrichment of 203 transcription factor motifs. Additionally, short double-stranded cfDNA signals at specific genomic regions correlate negatively with DNA methylation, positively with H3K4me3 histone modifications and gene transcription. The diagnostic potential of short double-stranded cell-free DNA (cfDNA) in blood plasma has not yet been recognized. When comparing short double-stranded cfDNA from patient samples of pancreatic ductal adenocarcinoma with colorectal carcinoma or septic with postoperative controls, we identify 136 and 241 differentially enriched loci, respectively. Using these differentially enriched loci, the disease types can be clearly distinguished by principal component analysis, demonstrating the diagnostic potential of short double-stranded cfDNA signals as a new class of biomarkers for liquid biopsies.

Semiconservative transmission of DNA N 6-adenine methylation in a unicellular eukaryote

Although DNA N 6-adenine methylation (6mA) is best known in prokaryotes, its presence in eukaryotes has recently generated great interest. Biochemical and genetic evidence supports that AMT1, an MT-A70 family methyltransferase (MTase), is crucial for 6mA deposition in unicellular eukaryotes. Nonetheless, the 6mA transmission mechanism remains to be elucidated. Taking advantage of single-molecule real-time circular consensus sequencing (SMRT CCS), here we provide definitive evidence for semiconservative transmission of 6mA in Tetrahymena thermophila In wild-type (WT) cells, 6mA occurs at the self-complementary ApT dinucleotide, mostly in full methylation (full-6mApT); after DNA replication, hemi-methylation (hemi-6mApT) is transiently present on the parental strand, opposite to the daughter strand readily labeled by 5-bromo-2'-deoxyuridine (BrdU). In ΔAMT1 cells, 6mA predominantly occurs as hemi-6mApT. Hemi-to-full conversion in WT cells is fast, robust, and processive, whereas de novo methylation in ΔAMT1 cells is slow and sporadic. In Tetrahymena, regularly spaced 6mA clusters coincide with the linker DNA of nucleosomes arrayed in the gene body. Importantly, in vitro methylation of human chromatin by the reconstituted AMT1 complex recapitulates preferential targeting of hemi-6mApT sites in linker DNA, supporting AMT1's intrinsic and autonomous role in maintenance methylation. We conclude that 6mA is transmitted by a semiconservative mechanism: full-6mApT is split by DNA replication into hemi-6mApT, which is restored to full-6mApT by AMT1-dependent maintenance methylation. Our study dissects AMT1-dependent maintenance methylation and AMT1-independent de novo methylation, reveals a 6mA transmission pathway with a striking similarity to 5-methylcytosine (5mC) transmission at the CpG dinucleotide, and establishes 6mA as a bona fide eukaryotic epigenetic mark.

FACT maintains chromatin architecture and thereby stimulates RNA polymerase II pausing during transcription in vivo

Facilitates chromatin transcription (FACT) is a histone chaperone that supports transcription through chromatin in vitro, but its functional roles in vivo remain unclear. Here, we analyze the in vivo functions of FACT with the use of multi-omics analysis after rapid FACT depletion from human cells. We show that FACT depletion destabilizes chromatin and leads to transcriptional defects, including defective promoter-proximal pausing and elongation, and increased premature termination of RNA polymerase II. Unexpectedly, our analysis revealed that promoter-proximal pausing depends not only on the negative elongation factor (NELF) but also on the +1 nucleosome, which is maintained by FACT.

DNA methylation patterns of transcription factor binding regions characterize their functional and evolutionary contexts

BACKGROUND

DNA methylation is an important epigenetic modification which has numerous roles in modulating genome function. Its levels are spatially correlated across the genome, typically high in repressed regions but low in transcription factor (TF) binding sites and active regulatory regions. However, the mechanisms establishing genome-wide and TF binding site methylation patterns are still unclear.

RESULTS

Here we use a comparative approach to investigate the association of DNA methylation to TF binding evolution in mammals. Specifically, we experimentally profile DNA methylation and combine this with published occupancy profiles of five distinct TFs (CTCF, CEBPA, HNF4A, ONECUT1, FOXA1) in the liver of five mammalian species (human, macaque, mouse, rat, dog). TF binding sites are lowly methylated, but they often also have intermediate methylation levels. Furthermore, biding sites are influenced by the methylation status of CpGs in their wider binding regions even when CpGs are absent from the core binding motif. Employing a classification and clustering approach, we extract distinct and species-conserved patterns of DNA methylation levels at TF binding regions. CEBPA, HNF4A, ONECUT1, and FOXA1 share the same methylation patterns, while CTCF's differ. These patterns characterize alternative functions and chromatin landscapes of TF-bound regions. Leveraging our phylogenetic framework, we find DNA methylation gain upon evolutionary loss of TF occupancy, indicating coordinated evolution. Furthermore, each methylation pattern has its own evolutionary trajectory reflecting its genomic contexts.

CONCLUSIONS

Our epigenomic analyses indicate a role for DNA methylation in TF binding changes across species including that specific DNA methylation profiles characterize TF binding and are associated with their regulatory activity, chromatin contexts, and evolutionary trajectories.

DNA fragility at topologically associated domain boundaries is promoted by alternative DNA secondary structure and topoisomerase II activity

CCCTC-binding factor (CTCF) binding sites are hotspots of genome instability. Although many factors have been associated with CTCF binding site fragility, no study has integrated all fragility-related factors to understand the mechanism(s) of how they work together. Using an unbiased, genome-wide approach, we found that DNA double-strand breaks (DSBs) are enriched at strong, but not weak, CTCF binding sites in five human cell types. Energetically favorable alternative DNA secondary structures underlie strong CTCF binding sites. These structures coincided with the location of topoisomerase II (TOP2) cleavage complex, suggesting that DNA secondary structure acts as a recognition sequence for TOP2 binding and cleavage at CTCF binding sites. Furthermore, CTCF knockdown significantly increased DSBs at strong CTCF binding sites and at CTCF sites that are located at topologically associated domain (TAD) boundaries. TAD boundary-associated CTCF sites that lost CTCF upon knockdown displayed increased DSBs when compared to the gained sites, and those lost sites are overrepresented with G-quadruplexes, suggesting that the structures act as boundary insulators in the absence of CTCF, and contribute to increased DSBs. These results model how alternative DNA secondary structures facilitate recruitment of TOP2 to CTCF binding sites, providing mechanistic insight into DNA fragility at CTCF binding sites.

Genome organization across scales: mechanistic insights from in vitro reconstitution studies

Eukaryotic genomes are compacted and organized into distinct three-dimensional (3D) structures, which range from small-scale nucleosome arrays to large-scale chromatin domains. These chromatin structures play an important role in the regulation of transcription and other nuclear processes. The molecular mechanisms that drive the formation of chromatin structures across scales and the relationship between chromatin structure and function remain incompletely understood. Because the processes involved are complex and interconnected, it is often challenging to dissect the underlying principles in the nuclear environment. Therefore, in vitro reconstitution systems provide a valuable approach to gain insight into the molecular mechanisms by which chromatin structures are formed and to determine the cause-consequence relationships between the processes involved. In this review, we give an overview of in vitro approaches that have been used to study chromatin structures across scales and how they have increased our understanding of the formation and function of these structures. We start by discussing in vitro studies that have given insight into the mechanisms of nucleosome positioning. Next, we discuss recent efforts to reconstitute larger-scale chromatin domains and loops and the resulting insights into the principles of genome organization. We conclude with an outlook on potential future applications of chromatin reconstitution systems and how they may contribute to answering open questions concerning chromatin architecture.

Apparent RNA bridging between PRC2 and chromatin is an artifact of non-specific chromatin precipitation upon RNA degradation

Polycomb repressive complex 2 (PRC2) modifies chromatin to maintain repression of genes specific for other cell lineages. In vitro, RNA inhibits PRC2 activity, but the effect of RNA on PRC2 in cells is less clear, with studies concluding that RNA either antagonizes or promotes PRC2 chromatin association. The addition of RNase A to chromatin immunoprecipitation reactions has been reported to reduce detection of PRC2 target sites, suggesting the existence of RNA bridges connecting PRC2 to chromatin. Here, we show that the apparent loss of PRC2 chromatin association after RNase A treatment is due to non-specific chromatin precipitation. RNA degradation precipitates chromatin out of solution, thereby masking enrichment of specific DNA sequences in chromatin immunoprecipitation reactions. Maintaining chromatin solubility by the addition of poly-L-glutamic acid rescues detection of PRC2 chromatin occupancy upon RNA degradation. These findings undermine support for the model that RNA bridges PRC2 and chromatin in cells.

In vitro reconstitution of chromatin domains shows a role for nucleosome positioning in 3D genome organization

Eukaryotic genomes are organized into chromatin domains. The molecular mechanisms driving the formation of these domains are difficult to dissect in vivo and remain poorly understood. Here we reconstitute Saccharomyces cerevisiae chromatin in vitro and determine its 3D organization at subnucleosome resolution by micrococcal nuclease-based chromosome conformation capture and molecular dynamics simulations. We show that regularly spaced and phased nucleosome arrays form chromatin domains in vitro that resemble domains in vivo. This demonstrates that neither loop extrusion nor transcription is required for basic domain formation in yeast. In addition, we find that the boundaries of reconstituted domains correspond to nucleosome-free regions and that insulation strength scales with their width. Finally, we show that domain compaction depends on nucleosome linker length, with longer linkers forming more compact structures. Together, our results demonstrate that regular nucleosome positioning is important for the formation of chromatin domains and provide a proof-of-principle for bottom-up 3D genome studies.

Cell-free DNA fragmentomics and second malignant neoplasm risk in patients with PTEN hamartoma tumor syndrome

Individuals with PTEN hamartoma tumor syndrome (PHTS) harbor pathogenic germline PTEN variants that confer a significantly increased lifetime risk of various organ-specific cancers including second primary malignant neoplasms (SMNs). Currently, there are no reliable biomarkers that can predict individual-level cancer risk. Despite the highly promising value of cell-free DNA (cfDNA) as a biomarker for underlying sporadic cancers, the utility of cfDNA in individuals with known cancer-associated germline variants and subclinical cancers remains poorly understood. We perform ultra-low-pass whole-genome sequencing (ULP-WGS) of cfDNA from plasma samples from patients with PHTS and cancer as well as those without cancer. Analysis of cfDNA reveals that patients with PHTS and SMNs have distinct cfDNA size distribution, aberrant genome-wide fragmentation, and differential fragment end motif frequencies. Our work provides evidence that cfDNA profiles may be used as a marker for SMN risk in patients with PHTS.

Cell-type differential targeting of SETDB1 prevents aberrant CTCF binding, chromatin looping, and cis-regulatory interactions

SETDB1 is an essential histone methyltransferase that deposits histone H3 lysine 9 trimethylation (H3K9me3) to transcriptionally repress genes and repetitive elements. The function of differential H3K9me3 enrichment between cell-types remains unclear. Here, we demonstrate mutual exclusivity of H3K9me3 and CTCF across mouse tissues from different developmental timepoints. We analyze SETDB1 depleted cells and discover that H3K9me3 prevents aberrant CTCF binding independently of DNA methylation and H3K9me2. Such sites are enriched with SINE B2 retrotransposons. Moreover, analysis of higher-order genome architecture reveals that large chromatin structures including topologically associated domains and subnuclear compartments, remain intact in SETDB1 depleted cells. However, chromatin loops and local 3D interactions are disrupted, leading to transcriptional changes by modifying pre-existing chromatin landscapes. Specific genes with altered expression show differential interactions with dysregulated cis-regulatory elements. Collectively, we find that cell-type specific targets of SETDB1 maintain cellular identities by modulating CTCF binding, which shape nuclear architecture and transcriptomic networks.

Corticosteroid-induced chromatin loop dynamics at the FKBP5 gene

FKBP5 is a 115-kb-long glucocorticoid-inducible gene implicated in psychiatric disorders. To investigate the complexities of chromatin interaction frequencies at the FKBP5 topologically associated domain (TAD), we deployed 15 one-to-all chromatin capture viewpoints near gene promoters, enhancers, introns, and CTCF-loop anchors. This revealed a "one-TAD-one-gene" structure encompassing the FKBP5 promoter and its enhancers. The FKBP5 promoter and its two glucocorticoid-stimulated enhancers roam the entire TAD while displaying subtle cell type-specific interactomes. The FKBP5 TAD consists of two nested CTCF loops that are coordinated by one CTCF site in the eighth intron of FKBP5 and another beyond its polyadenylation site, 61 kb further. Loop extension correlates with transcription increases through the intronic CTCF site. This is efficiently compensated for, since the short loop is restored even under high transcription regimes. The boundaries of the FKBP5 TAD consist of divergent CTCF site patterns, harbor multiple smaller genes, and are resilient to glucocorticoid stimulation. Interestingly, both FKBP5 TAD boundaries harbor H3K27me3-marked heterochromatin blocks that may reinforce them. We propose that cis-acting genetic and epigenetic polymorphisms underlying FKBP5 expression variation are likely to reside within a 240-kb region that consists of the FKBP5 TAD, its left sub-TAD, and both its boundaries.

High-resolution CTCF footprinting reveals impact of chromatin state on cohesin extrusion dynamics

DNA looping is vital for establishing many enhancer-promoter interactions. While CTCF is known to anchor many cohesin-mediated loops, the looped chromatin fiber appears to predominantly exist in a poorly characterized actively extruding state. To better characterize extruding chromatin loop structures, we used CTCF MNase HiChIP data to determine both CTCF binding at high resolution and 3D contact information. Here we present FactorFinder, a tool that identifies CTCF binding sites at near base-pair resolution. We leverage this substantial advance in resolution to determine that the fully extruded (CTCF-CTCF) state is rare genome-wide with locus-specific variation from ~1-10%. We further investigate the impact of chromatin state on loop extrusion dynamics, and find that active enhancers and RNA Pol II impede cohesin extrusion, facilitating an enrichment of enhancer-promoter contacts in the partially extruded loop state. We propose a model of topological regulation whereby the transient, partially extruded states play active roles in transcription.

An integrative single-cell multi-omics profiling of human pancreatic islets identifies T1D associated genes and regulatory signals

Genome wide association studies (GWAS) have identified over 100 signals associated with type 1 diabetes (T1D). However, translating any given T1D GWAS signal into mechanistic insights, including putative causal variants and the context (cell type and cell state) in which they function, has been limited. Here, we present a comprehensive multi-omic integrative analysis of single-cell/nucleus resolution profiles of gene expression and chromatin accessibility in healthy and autoantibody+ (AAB+) human islets, as well as islets under multiple T1D stimulatory conditions. We broadly nominate effector cell types for all T1D GWAS signals. We further nominated higher-resolution contexts, including effector cell types, regulatory elements, and genes for three independent T1D risk variants acting through islet cells within the pancreas at the DLK1/MEG3, RASGRP1, and TOX loci. Subsequently, we created isogenic gene knockouts DLK1-/-, RASGRP1-/-, and TOX-/-, and the corresponding regulatory region knockout, RASGRP1Δ, and DLK1Δ hESCs. Loss of RASGRP1 or DLK1, as well as knockout of the regulatory region of RASGRP1 or DLK1, increased β cell apoptosis. Additionally, pancreatic β cells derived from isogenic hESCs carrying the risk allele of rs3783355A/A exhibited increased β cell death. Finally, RNA-seq and ATAC-seq identified five genes upregulated in both RASGRP1-/- and DLK1-/- β-like cells, four of which are associated with T1D. Together, this work reports an integrative approach for combining single cell multi-omics, GWAS, and isogenic hESC-derived β-like cells to prioritize the T1D associated signals and their underlying context-specific cell types, genes, SNPs, and regulatory elements, to illuminate biological functions and molecular mechanisms.

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.

G4access identifies G-quadruplexes and their associations with open chromatin and imprinting control regions

Metazoan promoters are enriched in secondary DNA structure-forming motifs, such as G-quadruplexes (G4s). Here we describe 'G4access', an approach to isolate and sequence G4s associated with open chromatin via nuclease digestion. G4access is antibody- and crosslinking-independent and enriches for computationally predicted G4s (pG4s), most of which are confirmed in vitro. Using G4access in human and mouse cells, we identify cell-type-specific G4 enrichment correlated with nucleosome exclusion and promoter transcription. G4access allows measurement of variations in G4 repertoire usage following G4 ligand treatment, HDAC and G4 helicases inhibitors. Applying G4access to cells from reciprocal hybrid mouse crosses suggests a role for G4s in the control of active imprinting regions. Consistently, we also observed that G4access peaks are unmethylated, while methylation at pG4s correlates with nucleosome repositioning on DNA. Overall, our study provides a new tool for studying G4s in cellular dynamics and highlights their association with open chromatin, transcription and their antagonism to DNA methylation.

Histone exchange sensors reveal variant specific dynamics in mouse embryonic stem cells

Eviction of histones from nucleosomes and their exchange with newly synthesized or alternative variants is a central epigenetic determinant. Here, we define the genome-wide occupancy and exchange pattern of canonical and non-canonical histone variants in mouse embryonic stem cells by genetically encoded exchange sensors. While exchange of all measured variants scales with transcription, we describe variant-specific associations with transcription elongation and Polycomb binding. We found considerable exchange of H3.1 and H2B variants in heterochromatin and repeat elements, contrasting the occupancy and little exchange of H3.3 in these regions. This unexpected association between H3.3 occupancy and exchange of canonical variants is also evident in active promoters and enhancers, and further validated by reduced H3.1 dynamics following depletion of H3.3-specific chaperone, HIRA. Finally, analyzing transgenic mice harboring H3.1 or H3.3 sensors demonstrates the vast potential of this system for studying histone exchange and its impact on gene expression regulation in vivo.

Harnessing male germline epigenomics for the genetic improvement in cattle

Sperm is essential for successful artificial insemination in dairy cattle, and its quality can be influenced by both epigenetic modification and epigenetic inheritance. The bovine germline differentiation is characterized by epigenetic reprogramming, while intergenerational and transgenerational epigenetic inheritance can influence the offspring's development through the transmission of epigenetic features to the offspring via the germline. Therefore, the selection of bulls with superior sperm quality for the production and fertility traits requires a better understanding of the epigenetic mechanism and more accurate identifications of epigenetic biomarkers. We have comprehensively reviewed the current progress in the studies of bovine sperm epigenome in terms of both resources and biological discovery in order to provide perspectives on how to harness this valuable information for genetic improvement in the cattle breeding industry.

Single-Molecule Mapping of Chromatin Accessibility Using NOMe-seq/dSMF

The bulk of gene expression regulation in most organisms is accomplished through the action of transcription factors (TFs) on cis-regulatory elements (CREs). In eukaryotes, these CREs are generally characterized by nucleosomal depletion and thus higher physical accessibility of DNA. Many methods exploit this property to map regions of high average accessibility, and thus putative active CREs, in bulk. However, these techniques do not provide information about coordinated patterns of accessibility along the same DNA molecule, nor do they map the absolute levels of occupancy/accessibility. SMF (Single-Molecule Footprinting) fills these gaps by leveraging recombinant DNA cytosine methyltransferases (MTase) to mark accessible locations on individual DNA molecules. In this chapter, we discuss current methods and important considerations for performing SMF experiments.

Loss of the Ash2l subunit of histone H3K4 methyltransferase complexes reduces chromatin accessibility at promoters

Changes in gene expression programs are intimately linked to cell fate decisions. Post-translational modifications of core histones contribute to control gene expression. Methylation of lysine 4 of histone H3 (H3K4) correlates with active promoters and gene transcription. This modification is catalyzed by KMT2 methyltransferases, which require interaction with 4 core subunits, WDR5, RBBP5, ASH2L and DPY30, for catalytic activity. Ash2l is necessary for organismal development and for tissue homeostasis. In mouse embryo fibroblasts (MEFs), Ash2l loss results in gene repression, provoking a senescence phenotype. We now find that upon knockout of Ash2l both H3K4 mono- and tri-methylation (H3K4me1 and me3, respectively) were deregulated. In particular, loss of H3K4me3 at promoters correlated with gene repression, especially at CpG island promoters. Ash2l loss resulted in increased loading of histone H3 and reduced chromatin accessibility at promoters, accompanied by an increase of repressing and a decrease of activating histone marks. Moreover, we observed altered binding of CTCF upon Ash2l loss. Lost and gained binding was noticed at promoter-associated and intergenic sites, respectively. Thus, Ash2l loss and reduction of H3K4me3 correlate with altered chromatin accessibility and transcription factor binding. These findings contribute to a more detailed understanding of mechanistic consequences of H3K4me3 loss and associated repression of gene transcription and thus of the observed cellular consequences.

From Nucleus to Organs: Insights of Aryl Hydrocarbon Receptor Molecular Mechanisms

The aryl hydrocarbon receptor (AHR) is a markedly established regulator of a plethora of cellular and molecular processes. Its initial role in the detoxification of xenobiotic compounds has been partially overshadowed by its involvement in homeostatic and organ physiology processes. In fact, the discovery of its ability to bind specific target regulatory sequences has allowed for the understanding of how AHR modulates such processes. Thereby, AHR presents functions in transcriptional regulation, chromatin architecture modifications and participation in different key signaling pathways. Interestingly, such fields of influence end up affecting organ and tissue homeostasis, including regenerative response both to endogenous and exogenous stimuli. Therefore, from classical spheres such as canonical transcriptional regulation in embryonic development, cell migration, differentiation or tumor progression to modern approaches in epigenetics, senescence, immune system or microbiome, this review covers all aspects derived from the balance between regulation/deregulation of AHR and its physio-pathological consequences.

CTCF shapes chromatin structure and gene expression in health and disease

CCCTC-binding factor (CTCF) is an eleven zinc finger (ZF), multivalent transcriptional regulator, that recognizes numerous motifs thanks to the deployment of distinct combinations of its ZFs. The great majority of the ~50,000 genomic locations bound by the CTCF protein in a given cell type is intergenic, and a fraction of these sites overlaps with transcriptional enhancers. Furthermore, a proportion of the regions bound by CTCF intersect genes and promoters. This suggests multiple ways in which CTCF may impact gene expression. At promoters, CTCF can directly affect transcription. At more distal sites, CTCF may orchestrate interactions between regulatory elements and help separate eu- and heterochromatic areas in the genome, exerting a chromatin barrier function. In this review, we outline how CTCF contributes to the regulation of the three-dimensional structure of chromatin and the formation of chromatin domains. We discuss how CTCF binding and architectural functions are regulated. We examine the literature implicating CTCF in controlling gene expression in development and disease both by acting as an insulator and a factor facilitating regulatory elements to efficiently interact with each other in the nuclear space.

Transcription factor-nucleosome dynamics from plasma cfDNA identifies ER-driven states in breast cancer

Genome-wide binding profiles of estrogen receptor (ER) and FOXA1 reflect cancer state in ER+ breast cancer. However, routine profiling of tumor transcription factor (TF) binding is impractical in the clinic. Here, we show that plasma cell-free DNA (cfDNA) contains high-resolution ER and FOXA1 tumor binding profiles for breast cancer. Enrichment of TF footprints in plasma reflects the binding strength of the TF in originating tissue. We defined pure in vivo tumor TF signatures in plasma using ER+ breast cancer xenografts, which can distinguish xenografts with distinct ER states. Furthermore, state-specific ER-binding signatures can partition human breast tumors into groups with significantly different ER expression and mortality. Last, TF footprints in human plasma samples can identify the presence of ER+ breast cancer. Thus, plasma TF footprints enable minimally invasive mapping of the regulatory landscape of breast cancer in humans and open vast possibilities for clinical applications across multiple tumor types.

NURF301 contributes to gypsy chromatin insulator-mediated nuclear organization

Chromatin insulators are DNA-protein complexes that can prevent the spread of repressive chromatin and block communication between enhancers and promoters to regulate gene expression. In Drosophila, the gypsy chromatin insulator complex consists of three core proteins: CP190, Su(Hw), and Mod(mdg4)67.2. These factors concentrate at nuclear foci termed insulator bodies, and changes in insulator body localization have been observed in mutants defective for insulator function. Here, we identified NURF301/E(bx), a nucleosome remodeling factor, as a novel regulator of gypsy insulator body localization through a high-throughput RNAi imaging screen. NURF301 promotes gypsy-dependent insulator barrier activity and physically interacts with gypsy insulator proteins. Using ChIP-seq, we found that NURF301 co-localizes with insulator proteins genome-wide, and NURF301 promotes chromatin association of Su(Hw) and CP190 at gypsy insulator binding sites. These effects correlate with NURF301-dependent nucleosome repositioning. At the same time, CP190 and Su(Hw) both facilitate recruitment of NURF301 to chromatin. Finally, Oligopaint FISH combined with immunofluorescence revealed that NURF301 promotes 3D contact between insulator bodies and gypsy insulator DNA binding sites, and NURF301 is required for proper nuclear positioning of gypsy binding sites. Our data provide new insights into how a nucleosome remodeling factor and insulator proteins cooperatively contribute to nuclear organization.

Generating specificity in genome regulation through transcription factor sensitivity to chromatin

Cell type-specific gene expression relies on transcription factors (TFs) binding DNA sequence motifs embedded in chromatin. Understanding how motifs are accessed in chromatin is crucial to comprehend differential transcriptional responses and the phenotypic impact of sequence variation. Chromatin obstacles to TF binding range from DNA methylation to restriction of DNA access by nucleosomes depending on their position, composition and modification. In vivo and in vitro approaches now enable the study of TF binding in chromatin at unprecedented resolution. Emerging insights suggest that TFs vary in their ability to navigate chromatin states. However, it remains challenging to link binding and transcriptional outcomes to molecular characteristics of TFs or the local chromatin substrate. Here, we discuss our current understanding of how TFs access DNA in chromatin and novel techniques and directions towards a better understanding of this critical step in genome regulation.

Multiomic atlas with functional stratification and developmental dynamics of zebrafish cis-regulatory elements

Zebrafish, a popular organism for studying embryonic development and for modeling human diseases, has so far lacked a systematic functional annotation program akin to those in other animal models. To address this, we formed the international DANIO-CODE consortium and created a central repository to store and process zebrafish developmental functional genomic data. Our data coordination center ( https://danio-code.zfin.org ) combines a total of 1,802 sets of unpublished and re-analyzed published genomic data, which we used to improve existing annotations and show its utility in experimental design. We identified over 140,000 cis-regulatory elements throughout development, including classes with distinct features dependent on their activity in time and space. We delineated the distinct distance topology and chromatin features between regulatory elements active during zygotic genome activation and those active during organogenesis. Finally, we matched regulatory elements and epigenomic landscapes between zebrafish and mouse and predicted functional relationships between them beyond sequence similarity, thus extending the utility of zebrafish developmental genomics to mammals.

Mammalian PERIOD2 regulates H2A.Z incorporation in chromatin to orchestrate circadian negative feedback

Mammalian circadian oscillators are built on a feedback loop in which the activity of the transcription factor CLOCK-BMAL1 is repressed by the PER-CRY complex. Here, we show that murine Per^-/- fibroblasts display aberrant nucleosome occupancy around transcription start sites (TSSs) and at promoter-proximal and distal CTCF sites due to impaired histone H2A.Z deposition. Knocking out H2A.Z mimicked the Per null chromatin state and disrupted cellular rhythms. We found that endogenous mPER2 complexes retained CTCF as well as the specific H2A.Z-deposition chaperone YL1-a component of the ATP-dependent remodeler SRCAP and p400-TIP60 complex. While depleting YL1 or mutating chaperone-binding sites on H2A.Z lengthened the circadian period, H2A.Z deletion abrogated BMAL1 chromatin recruitment and promoted its proteasomal degradation. We propose that a PER2-mediated H2A.Z deposition pathway (1) compacts CLOCK-BMAL1 binding sites to establish negative feedback, (2) organizes circadian chromatin landscapes using CTCF and (3) bookmarks genomic loci for BMAL1 binding to impinge on the positive arm of the subsequent cycle.

Alcohol induced increases in sperm Histone H3 lysine 4 trimethylation correlate with increased placental CTCF occupancy and altered developmental programming

Using a mouse model, studies by our group reveal that paternal preconception alcohol intake affects offspring fetal-placental growth, with long-lasting consequences on adult metabolism. Here, we tested the hypothesis that chronic preconception male alcohol exposure impacts histone enrichment in sperm and that these changes are associated with altered developmental programming in the placenta. Using chromatin immunoprecipitation, we find alcohol-induced increases in sperm histone H3 lysine 4 trimethylation (H3K4me3) that map to promoters and presumptive enhancer regions enriched in genes driving neurogenesis and craniofacial development. Given the colocalization of H3K4me3 with the chromatin binding factor CTCF across both sperm and embryos, we next examined CTCF localization in the placenta. We find global changes in CTCF binding within placentae derived from the male offspring of alcohol-exposed sires. Furthermore, altered CTCF localization correlates with dysregulated gene expression across multiple gene clusters; however, these transcriptional changes only occur in male offspring. Finally, we identified a correlation between genomic regions exhibiting alcohol-induced increases in sperm H3K4me3 and increased CTCF binding in male placentae. Collectively, our analysis demonstrates that the chromatin landscape of sperm is sensitive to chronic alcohol exposure and that a subset of these affected regions exhibits increased placental CTCF enrichment.

CTCF: a misguided jack-of-all-trades in cancer cells

The emergence and progression of cancers is accompanied by a dysregulation of transcriptional programs. The three-dimensional (3D) organization of the human genome has emerged as an important multi-level mediator of gene transcription and regulation. In cancer cells, this organization can be restructured, providing a framework for the deregulation of gene activity. The CTCF protein, initially identified as the product from a tumor suppressor gene, is a jack-of-all-trades for the formation of 3D genome organization in normal cells. Here, we summarize how CTCF is involved in the multi-level organization of the human genome and we discuss emerging insights into how perturbed CTCF function and DNA binding causes the activation of oncogenes in cancer cells, mostly through a process of enhancer hijacking. Moreover, we highlight non-canonical functions of CTCF that can be relevant for the emergence of cancers as well. Finally, we provide guidelines for the computational identification of perturbed CTCF binding and reorganized 3D genome structure in cancer cells.

Cancer cells use self-inflicted DNA breaks to evade growth limits imposed by genotoxic stress

Genotoxic therapy such as radiation serves as a frontline cancer treatment, yet acquired resistance that leads to tumor reoccurrence is frequent. We found that cancer cells maintain viability during irradiation by reversibly increasing genome-wide DNA breaks, thereby limiting premature mitotic progression. We identify caspase-activated DNase (CAD) as the nuclease inflicting these de novo DNA lesions at defined loci, which are in proximity to chromatin-modifying CCCTC-binding factor (CTCF) sites. CAD nuclease activity is governed through phosphorylation by DNA damage response kinases, independent of caspase activity. In turn, loss of CAD activity impairs cell fate decisions, rendering cancer cells vulnerable to radiation-induced DNA double-strand breaks. Our observations highlight a cancer-selective survival adaptation, whereby tumor cells deploy regulated DNA breaks to delimit the detrimental effects of therapy-evoked DNA damage.

Aryl Hydrocarbon Receptor: From Homeostasis to Tumor Progression

Transcription factor aryl hydrocarbon receptor (AHR) has emerged as one of the main regulators involved both in different homeostatic cell functions and tumor progression. Being a member of the family of basic-helix-loop-helix (bHLH) transcriptional regulators, this intracellular receptor has become a key member in differentiation, pluripotency, chromatin dynamics and cell reprogramming processes, with plenty of new targets identified in the last decade. Besides this role in tissue homeostasis, one enthralling feature of AHR is its capacity of acting as an oncogene or tumor suppressor depending on the specific organ, tissue and cell type. Together with its well-known modulation of cell adhesion and migration in a cell-type specific manner in epithelial-mesenchymal transition (EMT), this duality has also contributed to the arise of its clinical interest, highlighting a new potential as therapeutic tool, diagnosis and prognosis marker. Therefore, a deregulation of AHR-controlled pathways may have a causal role in contributing to physiological and homeostatic failures, tumor progression and dissemination. With that firmly in mind, this review will address the remarkable capability of AHR to exert a different function influenced by the phenotype of the target cell and its potential consequences.

Overview of the Complex Relationship between Epigenetics Markers, CTG Repeat Instability and Symptoms in Myotonic Dystrophy Type 1

Among the trinucleotide repeat disorders, myotonic dystrophy type 1 (DM1) is one of the most complex neuromuscular diseases caused by an unstable CTG repeat expansion in the DMPK gene. DM1 patients exhibit high variability in the dynamics of CTG repeat instability and in the manifestations and progression of the disease. The largest expanded alleles are generally associated with the earliest and most severe clinical form. However, CTG repeat length alone is not sufficient to predict disease severity and progression, suggesting the involvement of other factors. Several data support the role of epigenetic alterations in clinical and genetic variability. By highlighting epigenetic alterations in DM1, this review provides a new avenue on how these changes can serve as biomarkers to predict clinical features and the mutation behavior.

Prediction and comparative analysis of CTCF binding sites based on a first principle approach

We calculated the patterns for the CCCTC transcription factor (CTCF) binding sites across many genomes on a first principle approach. The validation of the first principle method was done on the human as well as on the mouse genome. The predicted human CTCF binding sites are consistent with the consensus sequence, ChIP-seq data for the K562 cell, nucleosome positions for IMR90 cell as well as the CTCF binding sites in the mouse HOXA gene. The analysis of Homo sapiens, Mus musculus, Sus scrofa, Capra hircus and Drosophila melanogaster whole genomes shows: binding sites are organized in cluster-like groups, where two consecutive sites obey a power-law with coefficient ranging from to 0.3292 0.0068 to 0.5409 0.0064; the distance between these groups varies from 18.08 0.52kbp to 42.1 2.0kbp. The genome of Aedes aegypti does not show a power law, but 19.9% of binding sites are 144 4 and 287 5bp distant of each other. We run negative tests, confirming the under-representation of CTCF binding sites in Caenorhabditis elegans, Plasmodium falciparum and Arabidopsis thaliana complete genomes.

DNA Double Strand Break Repair and Its Control by Nucleosome Remodeling

DNA double strand breaks (DSBs) are repaired in eukaryotes by one of several cellular mechanisms. The decision-making process controlling DSB repair takes place at the step of DNA end resection, the nucleolytic processing of DNA ends, which generates single-stranded DNA overhangs. Dependent on the length of the overhang, a corresponding DSB repair mechanism is engaged. Interestingly, nucleosomes-the fundamental unit of chromatin-influence the activity of resection nucleases and nucleosome remodelers have emerged as key regulators of DSB repair. Nucleosome remodelers share a common enzymatic mechanism, but for global genome organization specific remodelers have been shown to exert distinct activities. Specifically, different remodelers have been found to slide and evict, position or edit nucleosomes. It is an open question whether the same remodelers exert the same function also in the context of DSBs. Here, we will review recent advances in our understanding of nucleosome remodelers at DSBs: to what extent nucleosome sliding, eviction, positioning and editing can be observed at DSBs and how these activities affect the DSB repair decision.

Single-Molecule Multikilobase-Scale Profiling of Chromatin Accessibility Using m6A-SMAC-Seq and m6A-CpG-GpC-SMAC-Seq

A hallmark feature of active cis-regulatory elements (CREs) in eukaryotes is their nucleosomal depletion and, accordingly, higher accessibility to enzymatic treatment. This property has been the basis of a number of sequencing-based assays for genome-wide identification and tracking the activity of CREs across different biological conditions, such as DNAse-seq, ATAC-seq , NOMeseq, and others. However, the fragmentation of DNA inherent to many of these assays and the limited read length of short-read sequencing platforms have so far not allowed the simultaneous measurement of the chromatin accessibility state of CREs located distally from each other. The combination of labeling accessible DNA with DNA modifications and nanopore sequencing has made it possible to develop such assays. Here, we provide a detailed protocol for carrying out the SMAC-seq assay (Single-Molecule long-read Accessible Chromatin mapping sequencing), in its m6A-SMAC-seq and m6A-CpG-GpC-SMAC-seq variants, together with methods for data processing and analysis, and discuss key experimental and analytical considerations for working with SMAC-seq datasets.

DeNOPA: decoding nucleosome positions sensitively with sparse ATAC-seq data

As the basal bricks, the dynamics and arrangement of nucleosomes orchestrate the higher architecture of chromatin in a fundamental way, thereby affecting almost all nuclear biology processes. Thanks to its rather simple protocol, assay for transposase-accessible chromatin using sequencing (ATAC)-seq has been rapidly adopted as a major tool for chromatin-accessible profiling at both bulk and single-cell levels; however, to picture the arrangement of nucleosomes per se remains a challenge with ATAC-seq. In the present work, we introduce a novel ATAC-seq analysis toolkit, named decoding nucleosome organization profile based on ATAC-seq data (deNOPA), to predict nucleosome positions. Assessments showed that deNOPA outperformed state-of-the-art tools with ultra-sparse ATAC-seq data, e.g. no more than 0.5 fragment per base pair. The remarkable performance of deNOPA was fueled by the short fragment reads, which compose nearly half of sequenced reads in the ATAC-seq libraries and are commonly discarded by state-of-the-art nucleosome positioning tools. However, we found that the short fragment reads enrich information on nucleosome positions and that the linker regions were predicted by reads from both short and long fragments using Gaussian smoothing. Last, using deNOPA, we showed that the dynamics of nucleosome organization may not directly couple with chromatin accessibility in the cis-regulatory regions when human cells respond to heat shock stimulation. Our deNOPA provides a powerful tool with which to analyze the dynamics of chromatin at nucleosome position level with ultra-sparse ATAC-seq data.