Search

Histone modifications play a crucial role in regulating chromatin architecture and gene expression. Here we develop a multiscale model for incorporating methylation in our nucleosome-resolution physics-based chromatin model to investigate the mechanisms by which H3K9 and H3K27 trimethylation (H3K9me3 and H3K27me3) influence chromatin structure and gene regulation. We apply three types of energy terms for this purpose: short-range potentials are derived from all-atom molecular dynamics simulations of wildtype and methylated chromatosomes, which revealed subtle local changes; medium-range potentials are derived by incorporating contacts between HP1 and nucleosomes modified by H3K9me3, to incorporate experimental results of enhanced contacts for short chromatin fibers (12 nucleosomes); for long-range interactions we identify H3K9me3- and H3K27me3-associated contacts based on Hi-C maps with a machine learning approach. These combined multiscale effects can model methylation as a first approximation in our mesoscale chromatin model, and applications to gene systems offer new insights into the epigenetic regulation of genomes mediated by H3K9me3 and H3K27me3., (© 2025 Author(s). Published under an exclusive license by AIP Publishing.)

DNA is heterogeneously packaged into chromatin, which is further organized into topologically associating domains (TADs) with sharp boundaries. These boundary locations are critical for genome regulation. Here, we explore how the distribution of DNA-packing density across chromatin affects the TAD boundary locations. We develop a polymer-physics-based model that utilizes DNA accessibility data to parameterize DNA-packing density along chromosomes, treating them as heteropolymers, and simulates the stochastic folding of these heteropolymers within a nucleus to yield a conformation ensemble. Such an ensemble reproduces a subset (over 60%) of TAD boundaries across the human genome, as confirmed by Hi-C data. Additionally, it reproduces the spatial distance matrices of 2-Mb genomic regions provided by FISH experiments. Furthermore, our model suggests that utilizing DNA accessibility data alone as input is sufficient to predict the emergence and disappearance of key TADs during early T cell differentiation. We show that stochastic folding of heteropolymers in a confined space can replicate both the prevalence of chromatin domain structures and the cell-to-cell variation in domain boundary positions observed in single-cell experiments. Furthermore, regions of lower DNA-packing density preferentially form domain boundaries, and this preference drives the emergence of TAD boundaries observed in ensemble-averaged Hi-C maps. The enrichment of TAD boundaries at CTCF binding sites can be attributed to the influence of CTCF binding on local DNA-packing density in our model. Collectively, our findings establish a strong link between TAD boundaries and regions of lower DNA-packing density, providing insights into the mechanisms underlying TAD formation., Competing Interests: Competing interests statement:The authors declare no competing interest.

Glioblastoma multiforme (GBM) progression is associated with changes in matrix stiffness, and different regions of the tumor niche exhibit distinct stiffnesses. Using elastic hydrogels, previous work has demonstrated that matrix stiffness modulates GBM behavior and drug responses. However, brain tissue is viscoelastic, and how stiffness impacts the GBM invasive phenotype and response to therapy within a viscoelastic niche remains largely unclear. Here, we report a three-dimensional (3D) viscoelastic GBM hydrogel system that models the stiffness heterogeneity present within the tumor niche. We find that GBM cells exhibit enhanced migratory ability, proliferation, and resistance to radiation in soft matrices, mimicking the tumor core and perifocal margins. Conversely, GBM cells remain confined and demonstrate increased resistance to chemotherapy in stiff matrices mimicking edematous tumor regions. We identify that stiffness-induced changes in the GBM phenotype are regulated by nuclear mechanosensing and chromatin condensation. Pharmacologically decondensing the chromatin significantly impedes GBM migration and overcomes stiffness-induced chemoresistance and radioresistance. Our findings highlight that stiffness regulates aggressive GBM behavior in viscoelastic matrices through mechanotransduction processes. Finally, we reveal the critical role of chromatin condensation in mediating GBM migration and therapy resistance, offering a potential new therapeutic target to improve GBM treatment outcomes.

The investigation of chromatin organization in single cells holds great promise for identifying causal relationships between genome structure and function. However, analysis of single-molecule data is hampered by extreme yet inherent heterogeneity, making it challenging to determine the contributions of individual chromatin fibers to bulk trends. To address this challenge, we propose ChromaFactor, a novel computational approach based on non-negative matrix factorization that deconvolves single-molecule chromatin organization datasets into their most salient primary components. ChromaFactor provides the ability to identify trends accounting for the maximum variance in the dataset while simultaneously describing the contribution of individual molecules to each component. Applying our approach to two single-molecule imaging datasets across different genomic scales, we find that these primary components demonstrate significant correlation with key functional phenotypes, including active transcription, enhancer-promoter distance, and genomic compartment. Also, we find that some bulk trends exist at the single-cell level, but only in a small fraction of cells, suggesting that critical changes in genome organization may be driven by specific rare subpopulations rather than occurring uniformly across all cells. ChromaFactor offers a robust tool for understanding the complex interplay between chromatin structure and function on individual DNA molecules, pinpointing which subpopulations drive functional changes and fostering new insights into cellular heterogeneity and its implications for bulk genomic phenomena., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2025 Gunsalus et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Proteins are subject to aging in the form of spontaneous, nonenzymatic post-translational modifications (PTMs). One such PTM is the formation of the β-linked isomer l-isoaspartic acid (isoAsp) from aspartic acid (Asp) or asparagine residues, which tends to occur in long-lived proteins. Histones can exhibit half-lives on the order of 100 days, and unsurprisingly, isoAsp formation has been observed in nearly every histone family. Delineating the molecular consequences of isoAsp formation in histones is challenging due to the multitude of processes that occur on such time scales. To isolate the effects of a specific isoAsp modification thus necessitates precise in vitro characterization with well-defined substrates. Here, we adapt a protein semisynthesis approach to generate full-length variants of histone H4 in which the canonical Asp at position 24 is replaced by its isoAsp isomer (H4isoD24). This variant was incorporated into chromatin templates, and the resulting constructs were used to interrogate key parameters of chromatin integrity and maintenance in vitro : compaction, nucleosome remodeling, and methylation of H4 lysine 20 (H4K20). Remarkably, despite its disruptive changes to the backbone's spacing and direction, isoD24 did not dramatically disrupt Mg 2+ -mediated chromatin self-association or nucleosome repositioning by the remodeler Chd1. In contrast, H4isoD24 significantly inhibited both Set8- and Suv4-20h1-catalyzed methylation at H4K20. These results suggest that H4isoD24 gives rise to a complex reorganization of the chromatin functional landscape, in which macroscopic processes show robustness and local mechanisms exhibit vulnerability to the presence of this mark.

Barrier-to-autointegration factor (BAF) associates with mitotic chromosomes and promotes nuclear envelope assembly by recruiting proteins, such as Lamins, required for the reconstruction of the nuclear envelope and lamina. BAF also mediates chromatin anchoring to the nuclear lamina via Lamin A/C. However, the mechanism by which BAF and Lamin A/C bind chromatin and affect the chromatin organization remains elusive. Here we report the cryo-electron microscopy structures of BAF-Lamin A/C-nucleosome complexes. We find that the BAF dimer complexed with the Lamin A/C IgF domain occupies the nucleosomal dyad position, forming a tripartite nucleosomal DNA binding structure. We also show that the Lamin A/C Lys486 and His506 residues, which are reportedly mutated in lipodystrophy patients, directly contact the DNA at the nucleosomal dyad. Excess BAF-Lamin A/C complexes symmetrically bind other nucleosomal DNA sites and connect two BAF-Lamin A/C-nucleosome complexes. Although the linker histone H1 competes with BAF-Lamin A/C binding at the nucleosomal dyad region, the two BAF-Lamin A/C molecules still bridge two nucleosomes. These findings provide insights into the mechanism by which BAF, Lamin A/C, and/or histone H1 bind nucleosomes and influence chromatin organization within the nucleus., Competing Interests: Competing interests: The authors declare no competing interests., (© 2025. The Author(s).)

Single-cell combinatorial indexing (sci) methods have addressed major limitations of throughput and cost for many single-cell modalities. With the incorporation of linear amplification and three-level barcoding in our suite of methods called sci-L3, we further addressed the limitations of uniformity in single-cell genome amplification. Here, we build on the generalizability of sci-L3 by extending it to template strand sequencing (sci-L3-Strand-seq), genome conformation capture (sci-L3-Hi-C), and the joint profiling of RNA and chromatin accessibility (sci-L3-RNA/ATAC). We demonstrate the ease of adapting sci-L3 to these new modalities by only requiring a single-step modification of the original protocol. As a proof of principle, we show our ability to detect sister chromatid exchanges, genome compartmentalization, and cell state-specific features in thousands of single cells. We anticipate sci-L3 to be compatible with additional modalities, including DNA methylation (sci-MET) and chromatin-associated factors (CUT&Tag), and ultimately enable a multi-omics readout of them., (© The Author(s) 2025. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Topologically associating domains (TADs) are essential components of three-dimensional (3D) genome organization and significantly influence gene transcription regulation. However, accurately identifying TADs from sparse chromatin contact maps and exploring the structural and functional elements within TADs remain challenging. To this end, we develop TADGATE, a graph attention auto-encoder that can generate imputed maps from sparse Hi-C contact maps while adaptively preserving or enhancing the underlying topological structures, thereby facilitating TAD identification. TADGATE captures specific attention patterns with two types of units within TADs and demonstrates TAD organization relates to chromatin compartmentalization with diverse biological properties. We identify many structural and functional elements within TADs, with their abundance reflecting the overall properties of these domains. We applied TADGATE to sparse and noisy Hi-C contact maps from 21 human tissues or cell lines. That improved the clarity of TAD structures, allowing us to investigate conserved and cell-type-specific boundaries and uncover cell-type-specific transcriptional regulatory mechanisms associated with topological domains. We also demonstrated TADGATE's capability to fill in sparse single-cell Hi-C contact maps and identify TAD-like domains within them, revealing the specific domain boundaries with distinct heterogeneity and the shared backbone boundaries characterized by strong CTCF enrichment and high gene expression levels., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

We analyze multiplexed fluorescence in situ hybridization (m-FISH) data for human and mouse cell lines. The m-FISH technique uses fluorescently-labeled single-stranded probes which hybridize to specific chromosomal regions, thereby allowing the measurement of the spatial positions of up to ∼100 tagged sites for several thousands of interphase chromosomes. Our analysis focuses on a wide range of different cell lines and two distinct organisms and provides a unified picture of chromatin structure for scales ranging from 5 kb (kilobases) up to 2 Mb (megabases), thus covering a genomic region of almost three orders of magnitude. Confirming recent analysis [Remini et al., Phys. Rev. E 109, 024408 (2024)], we show that there are two characteristic arrangements of chromatin referred to as phase α (crumpled globule) and phase β (looped domain) and discuss the physical properties of these phases. We show that a simple heterogeneous random walk model captures the main behavior observed in experiments and brings considerable insights into chromosomal structure., (© 2025 Author(s). Published under an exclusive license by AIP Publishing.)

Chromatin loop formation plays a crucial role in 3D genome interactions, with misfolding potentially leading to irregular gene expression and various diseases. While experimental tools such as Hi-C have advanced our understanding of genome interactions, the biophysical principles underlying chromatin loop formation remain elusive. This review examines computational approaches to chromatin folding, focusing on polymer models that elucidate chromatin loop mechanics. We discuss three key models: (1) the multi-loop-subcompartment model, which investigates the structural effects of loops on chromatin conformation; (2) the strings and binders switch model, capturing thermodynamic chromatin aggregation; and (3) the loop extrusion model, revealing the role of structural maintenance of chromosome complexes. In addition, we explore advanced models that address chromatin clustering heterogeneity in biological processes and disease progression. The review concludes with an outlook on open questions and current trends in chromatin loop formation and genome interactions, emphasizing the physical and computational challenges in the field., (© 2025 Author(s). Published under an exclusive license by AIP Publishing.)

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., Competing Interests: D.R.K. is an employee of Calico Life Sciences, LLC. All other authors have declared no competing interests exist., (Copyright: © 2025 Smaruj et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Motivation: Techniques such as high-throughput chromosome conformation capture (Hi-C) have provided a wealth of information on nucleus organization and genome important for understanding gene expression regulation. Genome-Wide Association Studies have identified numerous loci associated with complex traits. Expression quantitative trait loci (eQTL) studies have further linked the genetic variants to alteration in expression levels of associated target genes across individuals. However, the functional roles of many eQTLs in noncoding regions remain unclear. Current joint analyses of Hi-C and eQTLs data lack advanced computational tools, limiting what can be learned from these data., Results: We developed a computational method for simultaneous analysis of Hi-C and eQTL data, capable of identifying a small set of nonrandom interactions from all Hi-C interactions. Using these nonrandom interactions, we reconstructed large ensembles (×105) of high-resolution single-cell 3D chromatin conformations with thorough sampling, accurately replicating Hi-C measurements. Our results revealed many-body interactions in chromatin conformation at the single-cell level within eQTL loci, providing a detailed view of how 3D chromatin structures form the physical foundation for gene regulation, including how genetic variants of eQTLs affect the expression of associated eGenes. Furthermore, our method can deconvolve chromatin heterogeneity and investigate the spatial associations of eQTLs and eGenes at subpopulation level, revealing their regulatory impacts on gene expression. Together, ensemble modeling of thoroughly sampled single-cell chromatin conformations combined with eQTL data, helps decipher how 3D chromatin structures provide the physical basis for gene regulation, expression control, and aid in understanding the overall structure-function relationships of genome organization., Availability and Implementation: It is available at https://github.com/uic-liang-lab/3DChromFolding-eQTL-Loci., (© The Author(s) 2025. Published by Oxford University Press.)

The interphase chromatin structure is extremely complex, precise and dynamic. Experimental methods can only show the frequency of interaction of the various parts of the chromatin. Therefore, it is extremely important to develop theoretical methods to predict the chromatin structure. In this publication, we implemented an extended version of the SBS model described by Barbieri et al. and created the ChroMC program that is easy to use and freely available (https://github.com/regulomics/chroMC) to other users. We also describe the necessary factors for the effective modeling of the chromatin structure in Drosophila melanogaster . We compared results of chromatin structure predictions using two methods: Monte Carlo and Molecular Dynamic. Our simulations suggest that incorporating black, non-reactive chromatin is necessary for successful prediction of chromatin structure, while the loop extrusion model with a long range attraction potential or Lennard-Jones (with local attraction force) as well as using Hi-C data as input are not essential for the basic structure reconstruction. We also proposed a new way to calculate the similarity of the properties of contact maps including the calculation of local similarity.Communicated by Ramaswamy H. Sarma.

Identifying phase-separated structures remains challenging, and effective intervention methods are currently lacking 1 . Here we screened for phase-separated proteins in breast tumour cells and identified forkhead (FKH) box protein M1 (FOXM1) as the most prominent candidate. Oncogenic FOXM1 underwent liquid-liquid phase separation (LLPS) with FKH consensus DNA element, and compartmentalized the transcription apparatus in the nucleus, thereby sustaining chromatin accessibility and super-enhancer landscapes crucial for tumour metastatic outgrowth. Screening an epigenetics compound library identified AMPK agonists as suppressors of FOXM1 condensation. AMPK phosphorylated FOXM1 in the intrinsically disordered region (IDR), perturbing condensates, reducing oncogenic transcription, accumulating double-stranded DNA to stimulate innate immune responses, and endowing discrete FOXM1 with the ability to activate immunogenicity-related gene expressions. By developing a genetic code-expansion orthogonal system, we demonstrated that a phosphoryl moiety at a specific IDR1 site causes electrostatic repulsion, thereby abolishing FOXM1 LLPS and aggregation. A peptide targeting IDR1 and carrying the AMPK-phosphorylated residue was designed to disrupt FOXM1 LLPS and was shown to inhibit tumour malignancy, rescue tumour immunogenicity and improve tumour immunotherapy. Together, these findings provide novel and in-depth insights on function and mechanism of FOXM1 and develop methodologies that hold promising implications in clinics., Competing Interests: Competing interests: The authors have filed a patent application on this work (202411045018.3)., (© 2025. The Author(s), under exclusive licence to Springer Nature Limited.)

Cells can deform their local niche in three dimensions via whole-cell movements such as spreading, migration or volume expansion. These behaviours, occurring over hours to days, influence long-term cell fates including differentiation. Here we report a whole-cell movement that occurs in sliding hydrogels at the minutes timescale, termed cell tumbling, characterized by three-dimensional cell dynamics and hydrogel deformation elicited by heightened seconds-to-minutes-scale cytoskeletal and nuclear activity. Studies inhibiting or promoting the cell tumbling of mesenchymal stem cells show that this behaviour enhances differentiation into chondrocytes. Further, it is associated with a decrease in global chromatin accessibility, which is required for enhanced differentiation. Cell tumbling also occurs during differentiation into other lineages and its differentiation-enhancing effects are validated in various hydrogel platforms. Our results establish that cell tumbling is an additional regulator of stem cell differentiation, mediated by rapid niche deformation and nuclear mechanotransduction., Competing Interests: Competing interests: The authors declare no competing interests., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

Whilst DNA encodes our genetic blueprint as individual nucleobases, as well as epigenetic annotations in the form of biochemical marks, it also carries an extra layer of topological information -, the local over or underwinding of the double helix, known as DNA supercoiling. Supercoiling is a fundamental property of DNA that can be viewed as "topological epigenetics": it stores energy and structural information, and is tightly linked to fundamental processes; however, its quantification and study, by experiments and modelling alike, is challenging. We review experimental and simulation techniques to study supercoiling and its partition into twist and writhe, especially in the context of chromatin. We then discuss the dynamics of transcription-driven supercoiling in vitro and in vivo, and of supercoiling propagation along mammalian genomes. We finally provide evidence from the literature and potential mechanisms linking this ethereal topological mark to gene expression and chromosome instabilities in genetic diseases and cancer., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Cis-regulatory elements (CREs) control gene expression and are dynamic in their structure and function, reflecting changes in the composition of diverse effector proteins over time 1 . However, methods for measuring the organization of effector proteins at CREs across the genome are limited, hampering efforts to connect CRE structure to their function in cell fate and disease. Here we developed PRINT, a computational method that identifies footprints of DNA-protein interactions from bulk and single-cell chromatin accessibility data across multiple scales of protein size. Using these multiscale footprints, we created the seq2PRINT framework, which uses deep learning to allow precise inference of transcription factor and nucleosome binding and interprets regulatory logic at CREs. Applying seq2PRINT to single-cell chromatin accessibility data from human bone marrow, we observe sequential establishment and widening of CREs centred on pioneer factors across haematopoiesis. We further discover age-associated alterations in the structure of CREs in murine haematopoietic stem cells, including widespread reduction of nucleosome footprints and gain of de novo identified Ets composite motifs. Collectively, we establish a method for obtaining rich insights into DNA-binding protein dynamics from chromatin accessibility data, and reveal the architecture of regulatory elements across differentiation and ageing., Competing Interests: Competing interests: J.D.B. holds patents related to ATAC–seq, is a consultant for the Treehouse Family Foundation, and is a SAB member of Camp4 and seqWell. J.D.B. and S.M. hold a patent based on SHARE–seq. M.A.H. holds patents related to CRISPR-mediated interference and activation, and is a consultant for Akuous, Inc., DEM Biopharma and Gordian Biotechnology. A.J.W. is a scientific advisor for Frequency Therapeutics and Kate Therapeutics, is also a cofounder and scientific advisory board member and holds private equity in Elevian, Inc., a company that aims to develop medicines to restore regenerative capacity. Elevian also provides sponsored research to the Wagers laboratory. The other authors declare no competing interests., (© 2025. The Author(s).)

Cellular identity requires the concerted action of multiple transcription factors (TFs) bound together to enhancers of cell-type-specific genes. Despite TFs recognizing specific DNA motifs within accessible chromatin, this information is insufficient to explain how TFs select enhancers 1 . Here we compared four different TF combinations that induce different cell states, analysing TF genome occupancy, chromatin accessibility, nucleosome positioning and 3D genome organization at the nucleosome resolution. We show that motif recognition on mononucleosomes can decipher only the individual binding of TFs. When bound together, TFs act cooperatively or competitively to target nucleosome arrays with defined 3D organization, displaying motifs in particular patterns. In one combination, motif directionality funnels TF combinatorial binding along chromatin loops, before infiltrating laterally to adjacent enhancers. In other combinations, TFs assemble on motif-dense and highly interconnected loop junctions, and subsequently translocate to nearby lineage-specific sites. We propose a guided-search model in which motif grammar on nucleosome fibres acts as signpost elements, directing TF combinatorial binding to enhancers., Competing Interests: Competing interests: The authors declare no competing interests., (© 2024. The Author(s).)

Spatial organization of chromatin is essential for cellular functioning. However, the precise mechanisms governing sequence-dependent positioning of nucleosomes on DNA remain unknown in detail. Existing algorithms, considering the sequence-dependent deformability of DNA and its interactions with the histone globular domains, predict rotational setting of only 65% of human nucleosomes mapped in vivo. To uncover additional factors responsible for the nucleosome positioning, we analyzed potential involvement of the histone N-tails in this process. To this aim, we reconstituted the H2A/H4 N-tailless nucleosomes on human BRCA1 DNA (∼100 kb) and compared their positions and sequences with those of the wild-type nucleosomes. We found that removal of the histone N-tails promoted displacement of the predominant positions of nucleosomes, accompanied by redistribution of the AT-rich and GC-rich motifs in nucleosome sequences. Importantly, most of these sequence changes occurred at superhelical locations (SHLs) ±4, ±1, and ± 2, where the H2A and H4 N-tails interact with the DNA minor grooves. Furthermore, a substantial number of H4-tailless nucleosomes exhibit rotational settings opposite to that of the wild-type nucleosomes, the effect known to change the topological properties of chromatin fiber. Thus, the histone N-tails are operative in the selection of nucleosome positions, which may have wide-ranging implications for epigenetic modulation of chromatin states., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interest with the contents of this article., (Published by Elsevier Inc.)

During transcription, RNA polymerase II traverses through chromatin, and posttranslational modifications including histone methylations mark regions of active transcription. Histone protein H3 lysine 36 trimethylation (H3K36me3), which is established by the histone methyltransferase SET domain containing 2 (SETD2), suppresses cryptic transcription, regulates splicing, and serves as a binding site for transcription elongation factors. The mechanism by which the transcription machinery coordinates the deposition of H3K36me3 is not well understood. Here we provide cryo-electron microscopy structures of mammalian RNA polymerase II-DSIF-SPT6-PAF1c-TFIIS-IWS1-SETD2-nucleosome elongation complexes, revealing that the transcription machinery regulates H3K36me3 deposition by SETD2 on downstream and upstream nucleosomes. SPT6 binds the exposed H2A-H2B dimer during transcription, and the SPT6 death-like domain mediates an interaction with SETD2 bound to a nucleosome upstream of RNA polymerase II.

Breakthroughs in high-throughput sequencing and microscopic imaging technologies have revealed that chromatin structures vary considerably between cells of the same type. However, a thorough characterization of this heterogeneity remains elusive due to the labor-intensive and time-consuming nature of these experiments. To address these challenges, we introduce ChromoGen, a generative model based on state-of-the-art artificial intelligence techniques that efficiently predicts three-dimensional, single-cell chromatin conformations de novo with both region and cell type specificity. These generated conformations accurately reproduce experimental results at both the single-cell and population levels. Moreover, ChromoGen successfully transfers to cell types excluded from the training data using just DNA sequence and widely available DNase-seq data, thus providing access to chromatin structures in myriad cell types. These achievements come at a remarkably low computational cost. Therefore, ChromoGen enables the systematic investigation of single-cell chromatin organization, its heterogeneity, and its relationship to sequencing data, all while remaining economical.

Many transcription factors (TFs) have been shown to bind to super-enhancers, forming transcriptional condensates to activate transcription in various cellular systems. However, the genomic and epigenomic determinants of phase-separated transcriptional condensate formation remain poorly understood. Questions regarding which TFs tend to associate with transcriptional condensates and what factors influence their association are largely unanswered. Here we systematically analyzed 571 DNA sequence motifs across the human genome and 6650 TF binding profiles across different cell types to identify the molecular features contributing to the formation of transcriptional condensates. We found that the genomic distributions of sequence motifs for different TFs exhibit distinct clustering tendencies. Notably, TF motifs with a high genomic clustering tendency are significantly associated with super-enhancers. TF binding profiles showing a high genomic clustering tendency are further enriched at cell-type-specific super-enhancers. TFs with a high binding clustering tendency also possess high liquid-liquid phase separation abilities. Compared to nonclustered TF binding, densely clustered TF binding sites are more enriched at cell-type-specific super-enhancers with higher chromatin accessibility, elevated chromatin interaction and stronger association with cancer outcomes. Our results indicate that the clustered genomic binding patterns and the phase separation properties of TFs collectively contribute to the formation of transcriptional condensates., (© The Author(s) 2025. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Replication stress induced by exposure to extrinsic agents can lead to DNA breaks at common fragile sites, which are regions in the genome known to be prone to structural instability. The γH2A.X chromatin immunoprecipitation (ChIP) assay serves as a powerful tool in genotoxicity studies, as γH2A.X phosphorylation is a well-established marker for DNA double-strand breaks. Traditional γH2A.X ChIP assays, however, are often labor-intensive and involve multiple, time-consuming steps. In this study, we present a simplified yet effective method that combines subcellular fractionation with native ChIP to isolate γH2A.X-associated complexes. This approach is particularly suitable for analyzing γH2A.X-chromatin interactions with enhanced specificity and efficiency. Using subcellular fractionation, chromatin-unbound materials are effectively removed, resulting in a purified chromatin fraction. Subsequent micrococcal nuclease (MNase) digestion under mild conditions allows chromatin fragmentation while preserving physiological interactions between γH2A.X and its associated protein complexes. This preservation is essential for studying native interaction partners involved in DNA damage response pathways. This optimized native ChIP protocol substantially reduces the time and labor associated with conventional γH2A.X ChIP assays. The streamlined procedure not only simplifies the workflow but also yields highly reproducible results, making it particularly advantageous in settings where high-throughput processing of multiple samples is required. This method has broad applicability in studies focused on genome stability, DNA repair, and chromatin biology, where accurate and efficient detection of DNA damage sites is critical. By employing optimized protocols and streamlined steps, this method enables the detection of DNA damage at fragile sites with improved sensitivity and minimal sample handling, making it a valuable tool for studies on genome stability and DNA damage response.

Single-cell RNA-seq methods can be used to delineate cell types and states at unprecedented resolution but do little to explain why certain genes are expressed. Single-cell ATAC-seq and multiome (ATAC + RNA) have emerged to give a complementary view of the cell state. It is however unclear what additional information can be extracted from ATAC-seq data besides transcription factor binding sites. Here, we show that ATAC-seq telomere-like reads counter-inituively cannot be used to infer telomere length, as they mostly originate from the subtelomere, but can be used as a biomarker for chromatin condensation. Using long-read sequencing, we further show that modern hyperactive Tn5 does not duplicate 9 bp of its target sequence, contrary to common belief. We provide a new tool, Telomemore, which can quantify nonaligning subtelomeric reads. By analyzing several public datasets and generating new multiome fibroblast and B-cell atlases, we show how this new readout can aid single-cell data interpretation. We show how drivers of condensation processes can be inferred, and how it complements common RNA-seq-based cell cycle inference, which fails for monocytes. Telomemore-based analysis of the condensation state is thus a valuable complement to the single-cell analysis toolbox., (© The Author(s) 2025. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Eukaryotic DNA is packaged in the cell nucleus into chromatin, composed of arrays of DNA-histone protein octamer complexes, the nucleosomes. Over the past decade, it has become clear that chromatin structure in vivo is not a hierarchy of well-organized folded nucleosome fibers but displays considerable conformational variability and heterogeneity. In vitro and in vivo studies, as well as computational modeling, have revealed that attractive nucleosome-nucleosome interaction with an essential role of nucleosome stacking defines chromatin compaction. The internal structure of compacted nucleosome arrays is regulated by the flexible and dynamic histone N-terminal tails. Since DNA is a highly negatively charged polyelectrolyte, electrostatic forces make a decisive contribution to chromatin formation and require the histones, particularly histone tails, to carry a significant positive charge. This also results in an essential role of mobile cations of the cytoplasm (K+, Na+, Mg2+) in regulating electrostatic interactions. Building on a previously successfully established bottom-up coarse-grained (CG) nucleosome model, we have developed a CG nucleosome array (chromatin fiber) model with the explicit presence of mobile ions and studied its conformational variability as a function of Na+ and Mg2+ ion concentration. With progressively elevated ion concentrations, we identified four main conformational states of nucleosome arrays characterized as extended, flexible, nucleosome-clutched, and globular fibers., (© 2025 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution-NonCommercial 4.0 International (CC BY-NC) license (https://creativecommons.org/licenses/by-nc/4.0/).)

In single cells, variably sized nanoscale chromatin structures are observed, but it is unknown whether these form a cohesive framework that regulates RNA transcription. Here, we demonstrate that the human genome is an emergent, self-assembling, reinforcement learning system. Conformationally defined heterogeneous, nanoscopic packing domains form by the interplay of transcription, nucleosome remodeling, and loop extrusion. We show that packing domains are not topologically associated domains. Instead, packing domains exist across a structure-function life cycle that couples heterochromatin and transcription in situ, explaining how heterochromatin enzyme inhibition can produce a paradoxical decrease in transcription by destabilizing domain cores. Applied to development and aging, we show the pairing of heterochromatin and transcription at myogenic genes that could be disrupted by nuclear swelling. In sum, packing domains represent a foundation to explore the interactions of chromatin and transcription at the single-cell level in human health.

Complex deep learning models trained on very large datasets have become key enabling tools for current research in natural language processing and computer vision. By providing pre-trained models that can be fine-tuned for specific applications, they enable researchers to create accurate models with minimal effort and computational resources. Large scale genomics deep learning models come in two flavors: the first are large language models of DNA sequences trained in a self-supervised fashion, similar to the corresponding natural language models; the second are supervised learning models that leverage large scale genomics datasets from ENCODE and other sources. We argue that these models are the equivalent of foundation models in natural language processing in their utility, as they encode within them chromatin state in its different aspects, providing useful representations that allow quick deployment of accurate models of gene regulation. We demonstrate this premise by leveraging the recently created Sei model to develop simple, interpretable models of intron retention, and demonstrate their advantage over models based on the DNA language model DNABERT-2. Our work also demonstrates the impact of chromatin state on the regulation of intron retention. Using representations learned by Sei, our model is able to discover the involvement of transcription factors and chromatin marks in regulating intron retention, providing better accuracy than a recently published custom model developed for this purpose., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2025 Daoud, Ben-Hur. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

In the last years, it has been proved that some viruses are able to re-structure chromatin organization and alter the epigenomic landscape of the host genome. In addition, they are able to affect the physical mechanisms shaping chromatin 3D structure, with a consequent impact on gene activity. Here, we investigate with polymer physics genome re-organization of the host genome upon SARS-CoV-2 viral infection and how it can impact structural variability within the population of single-cell chromatin configurations. Using published Hi-C data and molecular dynamics simulations, we build ensembles of 3D configurations representing single-cell chromatin conformations in control and SARS-CoV-2 infected conditions. We focus on genomic length scales of TADs and consider, as a case study, models of real loci containing DDX58 and IL6 genes, belonging, respectively, to the antiviral interferon response and pro-inflammatory genes. Clustering analysis applied to the ensemble of polymer configurations reveals a generally increased variability and a more heterogeneous population of 3D structures in infected conditions. This points toward a scenario in which viral infection leads to a loss of chromatin structural specificity with, likely, a consequent impact on the correct regulation of host cell genes., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

Investigating the relationship, particularly the lead-lag effect, between time series is a common question across various disciplines, especially when uncovering biological processes. However, analyzing time series presents several challenges. Firstly, due to technical reasons, the time points at which observations are made are not at uniform intervals. Secondly, some lead-lag effects are transient, necessitating time-lag estimation based on a limited number of time points. Thirdly, external factors also impact these time series, requiring a similarity metric to assess the lead-lag relationship. To counter these issues, we introduce a model grounded in the Gaussian process, affording the flexibility to estimate lead-lag effects for irregular time series. In addition, our method outputs dissimilarity scores, thereby broadening its applications to include tasks such as ranking or clustering multiple pairwise time series when considering their strength of lead-lag effects with external factors. Crucially, we offer a series of theoretical proofs to substantiate the validity of our proposed kernels and the identifiability of kernel parameters. Our model demonstrates advances in various simulations and real-world applications, particularly in the study of dynamic chromatin interactions, compared to other leading methods., (© The Author(s) 2025. Published by Oxford University Press on behalf of The International Biometric Society.)

The three-dimensional (3D) organization of the human genome plays a crucial role in gene regulation. EXPloration of Regulatory Epigenome with Spatial and Sequence Observations (EXPRESSO) is a novel multi-omics database for exploration and visualization of multi-layered 3D genomic features across 46 different human tissues. Integrating 1360 3D genomic datasets (Hi-C, HiChIP, ChIA-PET) and 842 1D genomic and transcriptomic datasets (ChIP-seq, ATAC-seq, RNA-seq) from the same biosample, EXPRESSO provides a comprehensive resource for studying the interplay between 3D genome architecture and transcription regulation. This database offers diverse 3D genomic feature types (compartments, contact matrix, contact domains, stripes as diagonal lines extending from a genomic locus in contact matrix, chromatin loops, etc.) and user-friendly interface for both data exploration and download. Other key features include REpresentational State Transfer application programming interfaces for programmatic access, advanced visualization tools for 3D genomic features and web-based applications that correlate 3D genomic features with gene expression and epigenomic modifications. By providing extensive datasets and tools, EXPRESSO aims to deepen our understanding of 3D genomic architecture and its implications for human health and disease, serving as a vital resource for the research community. EXPRESSO is freely available at https://expresso.sustech.edu.cn., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Gene regulation is inherently multiscale, but scale-adaptive machine learning methods that fully exploit this property in single-nucleus accessibility data are still lacking. Here, we develop ChromatinHD, a pair of scale-adaptive models that uses the raw accessibility data, without peak-calling or windows, to link regions to gene expression and determine differentially accessible chromatin. We show how ChromatinHD consistently outperforms existing peak and window-based approaches and find that this is due to a large number of uniquely captured, functional accessibility changes within and outside of putative cis-regulatory regions. Furthermore, ChromatinHD can delineate collaborating regulatory regions, including their preferential genomic conformations, that drive gene expression. Finally, our models also use changes in ATAC-seq fragment lengths to identify dense binding of transcription factors, a feature not captured by footprinting methods. Altogether, ChromatinHD, available at https://chromatinhd.org , is a suite of computational tools that enables a data-driven understanding of chromatin accessibility at various scales and how it relates to gene expression., Competing Interests: Competing interests: The authors declare no competing interests., (© 2024. The Author(s).)

The intricate interplay between lamins and chromatin underpins the structural integrity and functional organization of the eukaryotic nucleus. Lamins, type V intermediate filament proteins, form a robust meshwork beneath the inner nuclear membrane that is crucial for sustaining nuclear architecture through interactions with lamin-associated domains (LADs). LADs are predominantly heterochromatic regions in which compacted chromatin is enriched at the nuclear periphery, interacting with lamins and lamin-associated proteins. Disruptions of these interactions are implicated in a spectrum of diseases, including laminopathies, cancer, and age-related pathologies, highlighting the importance of lamin-LAD interactions. Thus, a detailed understanding of lamin-chromatin interactions may provide new insights into chromatin organization and shed light on the mechanism behind certain disease states. Here, we discuss the current state of knowledge of lamin-chromatin interactions from a biochemical and structural point of view., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

In eukaryotic cells, transcription by RNA polymerase II occurs in the context of chromatin, requiring the transcription machinery to navigate through nucleosomes as it traverses gene bodies. Recent advances in structural biology have provided unprecedented insights into the mechanisms underlying transcription elongation. This review presents a structural perspective on transcription through chromatin, focusing on the latest findings from high-resolution structures of transcribing RNA polymerase II-nucleosome complexes. I discuss how RNA polymerase II, in concert with elongation factors such as SPT4/5, SPT6, ELOF1, and the PAF1 complex, engages with and transcribes through nucleosomes. The review examines the stepwise unwrapping of nucleosomal DNA as polymerase advances, the roles of elongation factors in facilitating this process, and the mechanisms of nucleosome retention and transfer during transcription. This structural perspective provides a foundation for understanding the intricate interplay between the transcription machinery and chromatin, offering insights into how cells balance the need for genetic accessibility with the maintenance of genome stability and epigenetic regulation., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Recent analyses revealed the essential function of chromatin structure in maintaining and regulating genomic information. Advancements in microscopy, nuclear structure observation techniques, and the development of methods utilizing next-generation sequencers (NGSs) have significantly progressed these discoveries. Methods utilizing NGS enable genome-wide analysis, which is challenging with microscopy, and have elucidated concepts of important chromatin structures such as a loop structure, a domain structure called topologically associating domains (TADs), and compartments. In this chapter, I introduce chromatin interaction techniques using NGS and outline the principles and features of each method., (© 2025. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

Hi-C methods reveal 3D genome features but lack correspondence to dynamic chromatin behavior. PHi-C2, Python software, addresses this gap by transforming Hi-C data into polymer models. After the optimization algorithm, it enables us to calculate 3D conformations and conduct dynamic simulations, providing insights into chromatin dynamics, including the mean-squared displacement and rheological properties. This chapter introduces PHi-C2 usage, offering a tutorial for comprehensive 4D genome analysis., (© 2025. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

DNA supercoiling is a biophysical feature of the double helix with a pivotal role in biological processes. However, understanding of DNA supercoiling in the chromatin remains limited. Here, we developed azide-trimethylpsoralen sequencing (ATMP-seq), a DNA supercoiling assay offering quantitative accuracy while minimizing genomic bias and background noise. Using ATMP-seq, we directly visualized transcription-dependent negative and positive twin-supercoiled domains around genes and mapped kilobase-resolution DNA supercoiling throughout the human genome. Remarkably, we discovered megabase-scale supercoiling domains (SDs) across all chromosomes that are modulated mainly by topoisomerases I and IIβ. Transcription activities, but not the consequent supercoiling accumulation in the local region, contribute to SD formation, indicating the long-range propagation of transcription-generated supercoiling. Genome-wide SDs colocalize with A/B compartments in both human and Drosophila cells but are distinct from topologically associating domains (TADs), with negative supercoiling accumulation at TAD boundaries. Furthermore, genome-wide DNA supercoiling varies between cell states and types and regulates human gene expression, underscoring the importance of supercoiling dynamics in chromatin regulation and function., Competing Interests: Competing interests: The authors declare no competing interests., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Hi-C is a popular ligation-based technique to detect 3D physical chromosome structure within the nucleus using cross-linking and next-generation sequencing. As an unbiased genome-wide assay based on chromosome conformation capture, it provides rich insights into chromosome structure, dynamic chromosome folding and interactions, and the regulatory state of a cell. Bioinformatics analyses of Hi-C data require dedicated protocols as most genome alignment tools assume that both paired-end reads will map to the same chromosome, resulting in large two-dimensional matrices as processed data. Here, we outline the necessary steps to generate high-quality aligned Hi-C data by separately mapping each read while correcting for biases from restriction enzyme digests. We introduce our own custom open-source pipeline, which enables users to select an aligner of their choosing with high accuracy and performance. This enables users to generate high-resolution datasets with fast turnaround and fewer unmapped reads. Finally, we discuss recent innovations in experimental techniques, bioinformatics techniques, and their applications in clinical testing for diagnostics., (© 2025. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

Chromatin modeling enables the characterization of chromatin architecture at a resolution so far unachievable with experimental techniques. Polymer models fill our knowledge gap on a wide range of structures, from chromatin loops to nuclear compartments. Many physical properties already known for polymers can thus explain the dynamics of chromatin. With molecular simulations, it is possible to probe an ensemble of conformations, which attest to the variability observed in individual cells and the general behavior of a population of cells. In this review, we describe universal characteristics of polymers that chromatin carries. We introduce how these characteristics can be assessed with polymer simulations while also addressing specific aspects of chromatin and its environment. Finally, we give examples of plant chromatin models that, despite their paucity, augur well for the future of polymer simulations to plant chromosome biology., (© 2025. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

Histones are positively charged proteins found in the chromatin of eukaryotic cells. They regulate gene expression and are required for the organization and packaging of DNA within the nucleus. Histones are extremely conserved, allowing for transcription, replication, and repair. This review delves into their complex structure and function in DNA assembly, their role in nucleosome assembly, and the higher-order chromatin structures they generate. We look at the five different types of histone proteins: H1, H2A, H2B, H3, H4, and their variations. These histones bind with DNA to produce nucleosomes, the basic units of chromatin that are essential for compacting DNA and controlling its accessibility. Their dynamic control of chromatin accessibility has important implications for genomic stability and cellular activities. We elucidate regulatory mechanisms in both normal and pathological situations by investigating their structural features, diverse interaction mechanisms, and chromatin impact. In addition, we discuss the functions of histone post-translational modifications (PTMs) and their significance in various disorders. These alterations, which include methylation, acetylation, phosphorylation, and ubiquitination, are crucial in regulating histone function and chromatin dynamics. We specifically describe and explore the role of changed histones in the evolution of cancer, neurological disorders, sepsis, autoimmune illnesses, and inflammatory conditions. This comprehensive review emphasizes histone's critical role in genomic integrity and their potential as therapeutic targets in various diseases., (Copyright © 2025. Published by Elsevier Inc.)

The formation of biomolecular condensates in vitro and in vivo has become an increasingly important subject of studies. One particular area of interest is the phase separation of chromatin in the nucleus. However, the interplay of condensed chromatin and chromatin-binding enzymes has barely been studied as of now. Here, we show an optical tweezer-based assay that uses controlled fusion of two condensates to probe the effect of enzymes, such as chromatin remodeling motor proteins, on the fluidity of condensates. The assay provides a powerful tool that enables the study of if and how chromatin-enzyme interactions alter biophysical properties in these dense chromatin condensates., (© 2025. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

Polymer modeling has been playing an increasingly important role in complementing 3D genome experiments, both to aid their interpretation and to reveal the underlying molecular mechanisms. This chapter illustrates an application of Hi-C metainference, a Bayesian approach to explore the 3D organization of a target genomic region by integrating experimental contact frequencies into a prior model of chromatin. The method reconstructs the conformational ensemble of the target locus by combining molecular dynamics simulation and Monte Carlo sampling from the posterior probability distribution given the data. Using prior chromatin models at both 1 kb and nucleosome resolution, we apply this approach to a 30 kb locus of mouse embryonic stem cells consisting of two well-defined domains linking several gene promoters together. Retaining the advantages of both physics-based and data-driven strategies, Hi-C metainference can provide an experimentally consistent representation of the system while at the same time retaining molecular details necessary to derive physical insights., (© 2025. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

In organisms ranging from vertebrates to plants, major components of centromeres are rapidly evolving repeat sequences, such as tandem repeats (TRs) and transposable elements (TEs), which harbour centromere-specific histone H3 (CENH3) 1,2 . Complete centromere structures recently determined in human and Arabidopsis suggest frequent integration and purging of retrotransposons within the TR regions of centromeres 3-5 . Despite the high impact of 'centrophilic' retrotransposons on the paradox of rapid centromere evolution, the mechanisms involved in centromere targeting remain poorly understood in any organism. Here we show that both Ty3 and Ty1 long terminal repeat retrotransposons rapidly turnover within the centromeric TRs of Arabidopsis species. We demonstrate that the Ty1/Copia element Tal1 (Transposon of Arabidopsis lyrata 1) integrates de novo into regions occupied by CENH3 in Arabidopsis thaliana, and that ectopic expansion of the CENH3 region results in spread of Tal1 integration regions. The integration spectra of chimeric TEs reveal the key structural variations responsible for contrasting chromatin-targeting specificities to centromeres versus gene-rich regions, which have recurrently converted during the evolution of these TEs. Our findings show the impact of centromeric chromatin on TE-mediated rapid centromere evolution, with relevance across eukaryotic genomes., Competing Interests: Competing interests: The authors declare no competing interests., (© 2025. The Author(s).)

Pioneer factors are transcription factors capable of binding to nucleosomal DNA, initiating chromatin opening, and facilitating gene expression. By overcoming nucleosomes, pioneer factors enable cellular reprogramming, tissue-specific gene expression, and genome response to external stimuli. Here we discuss the recent literature on how pioneer factors modulate chromatin architecture at multiple levels, from local chromatin accessibility to large-scale genome organization, including chromatin compartments, topologically associating domains, and enhancer-promoter looping. Understanding the mechanisms by which pioneer factors modulate chromatin organization dynamics is key to understand their broader impact on gene expression regulation., Competing Interests: Declaration of competing interest None., (Copyright © 2025 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Similarly to life in our planet, where thousands of species inhabit the same ecosystem, the cell nucleus hosts different essential processes that share the same territory, making the interaction between them unavoidable. DNA replication and transcription are essential processes that copy and decode the information contained in our genomes, sharing -and competing for- the same chromatin template. Both activities are executed by large macromolecular machines with similar requirements to access the DNA, remodel the nucleosomes ahead of them and reassemble the chromatin make-up behind. Mechanistically, both processes cannot simultaneously act on the same DNA sequence, but emerging evidence shows that they frequently interact. Here we revise recent data on how transcription and replication occur in chromatin highlighting the symbiotic relationship between both processes, which might help explain how their activities contribute to shape the structure and function of the mammalian genome., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have influenced the work reported in this paper., (Copyright © 2025 The Authors. Published by Elsevier Ltd.. All rights reserved.)

The CRISPR-associated endonuclease Cas9 derived from prokaryotes is used as a genome editing, which targets specific genomic loci by single guide RNAs (sgRNAs). The eukaryotes, the target of genome editing, store their genome DNA in chromatin, in which the nucleosome is a basic unit. Despite previous structural analyses focusing on Cas9 cleaving free DNA, structural insights into Cas9 targeting of DNA within nucleosomes are limited, leading to uncertainties in understanding how Cas9 operates in the eukaryotic genome. In the present study, we perform native-polyacrylamide gel electrophoresis (PAGE)  analyses and find that Cas9 targets the linker DNA and the entry-exit DNA region of the nucleosome but not the DNA tightly wrapped around the histone octamer. We further determine cryo-electron microscopy (cryo-EM) structure of the Cas9-sgRNA-nucleosome ternary complex that targets linker DNA in nucleosomes. The structure suggests interactions between Cas9 and nucleosomes at multiple sites. Mutants that reduce the interaction between nucleosomal DNA and Cas9 improve nucleosomal DNA cleavage activity in vitro, although inhibition by the interaction between Cas9 and nucleosomes is limited in vivo. These findings will contribute to the development of novel genome editing tools in chromatin., Competing Interests: Competing interests: O.N. is a co-founder, board member and scientific advisor for Curreio. The remaining authors declare no competing interests., (© 2024. The Author(s).)

Sequencing-based Hi-C technology has been widely used to study the three-dimensional structure of chromatin. More recently, the development of single-cell Hi-C technology has enabled the study of chromatin structural variations between individual cells. However, single-cell Hi-C data are often highly sparse, necessitating the use of imputation algorithms to address insufficient sampling. Current methods encounter challenges such as significant discrepancies from real structural features, limited reproducibility, slower computational speeds, or reliance on large amounts of training data, which hinder their broader applicability. In this study, we improved the previously published CTG (Hi-C To Geometry) algorithm to introduce the single-cell CTG (scCTG) algorithm, which combines convolution and diffusion processes to yield the spatial distance matrix for various types of single-cell chromatin structure data. scCTG algorithm shows a good performance in terms of computational efficiency, robustness, and correlation with physical spatial distances. The scCTG algorithm can be applied to effectively identify compartments and insulation strength for each locus, providing deeper insights into the relationship between chromatin structure and gene expression at the single-cell level., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

The dynamic organization of chromatin plays an essential role in the regulation of genetic activity, interconverting between open and compact forms at the global level. The mechanisms underlying these large-scale changes remain a topic of widespread interest. The simulations of nucleosome-decorated DNA reported herein reveal profound effects of the nucleosome itself on overall chromatin properties. Models that capture the long-range communication between proteins on nucleosome-decorated DNA chains incorporate DNA pathways different from those that were previously proposed based on ultracentrifugation and chemical cross-linking data. New quantitative biochemical assays measuring the rates of communication between interacting proteins bound to a promoter and an enhancer at the ends of saturated, precisely positioned, nucleosome-decorated DNA chains reveal a chromatin architecture with a three-nucleosome repeat, a model inconsistent with the two-start configurations deduced from earlier physical studies. Accompanying computations uncover small differences in the twisting of successive base pairs that seemingly give rise to the observed global properties. These data suggest that the novel state of chromatin determined under physiological conditions differs from that deduced under standard physical conditions, likely reflecting the different salt conditions used in the two types of experiments. This novel chromatin state may be important for a number of DNA transactions that occur in the cell nucleus., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

We demonstrated a novel approach for protein-protein interaction (PPI) profiling of histone H3 using intracellular photocatalytic-proximity labeling (iPPL). This approach identified that the combination of acriflavine as a photocatalyst and 1-methyl-4-arylurazol (MAUra) as a protein labeling agent was the most efficient strategy to proceed the protein proximity labeling reaction. Furthermore, the identification of the labeled amino acids in histone H3 interacting proteins, histone lysine N -methyltransferase EZH2, showed that the amino acid in EZH2 within a few nanometers from histone H3 is labeled by iPPL. This restricted labeling radius allows for more-focused PPI profiling, compared to conventional proximity labeling methods.

The scarcity of zebrafish-specific motif databases presents a challenge to the analysis of transcription factor (TF) motif within zebrafish single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) data, thus hindering the identification of regulatory elements throughout zebrafish embryonic development. Here, we provide a protocol to analyze single-nucleus chromatin accessibility dataset during zebrafish early embryogenesis. We describe steps for fragment file retrieval, sample integration, quality control, Latent Semantic Indexing (LSI) clustering, and peak calling via ArchR. Crucially, we detail procedures for zebrafish-specific motif database construction, motif enrichment, and TF footprinting analysis. For complete details on the use and execution of this protocol, please refer to Lin et al. 1 ., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

Here, we present a protocol for using d-rG4-seq, a technique for mapping RNA G-quadruplex (rG4) for chromatin-bound RNA. We describe steps for identifying in vivo rG4 structures based on differential sensitivity of rG4 to dimethyl sulfate (DMS) modification, folding in the presence of monovalent cations, K+ versus Li+, and reverse transcriptase (RT) readthrough when folded. We then detail procedures for isolating RNA from the chromatin-bound fractions to enrich for epigenetic regulators and comparing in vitro versus in vivo profiles. For complete details on the use and execution of this protocol, please refer to Lee et al. 1 ., Competing Interests: Declaration of interests J.T.L. is a cofounder of Translate Bio and Fulcrum, an advisor to Skyhawk Therapeutics, and a Non-Executive Director of the GSK., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)