Organization and function of the 3D genome

Chromatin phase separation and nuclear shape fluctuations are correlated in a polymer model of the nucleus

Abnormal cell nuclear shapes are hallmarks of diseases, including progeria, muscular dystrophy, and many cancers. Experiments have shown that disruption of heterochromatin and increases in euchromatin lead to nuclear deformations, such as blebs and ruptures. However, the physical mechanisms through which chromatin governs nuclear shape are poorly understood. To investigate how heterochromatin and euchromatin might govern nuclear morphology, we studied chromatin microphase separation in a composite coarse-grained polymer and elastic shell simulation model. By varying chromatin density, heterochromatin composition, and heterochromatin-lamina interactions, we show how the chromatin phase organization may perturb nuclear shape. Increasing chromatin density stabilizes the lamina against large fluctuations. However, increasing heterochromatin levels or heterochromatin-lamina interactions enhances nuclear shape fluctuations by a "wetting"-like interaction. In contrast, fluctuations are insensitive to heterochromatin's internal structure. Our simulations suggest that peripheral heterochromatin accumulation could perturb nuclear morphology, while nuclear shape stabilization likely occurs through mechanisms other than chromatin microphase organization.

An atlas of the tomato epigenome reveals that KRYPTONITE shapes TAD-like boundaries through the control of H3K9ac distribution

In recent years, the exploration of genome three-dimensional (3D) conformation has yielded profound insights into the regulation of gene expression and cellular functions in both animals and plants. While animals exhibit a characteristic genome topology defined by topologically associating domains (TADs), plants display similar features with a more diverse conformation across species. Employing advanced high-throughput sequencing and microscopy techniques, we investigated the landscape of 26 histone modifications and RNA polymerase II distribution in tomato (Solanum lycopersicum). Our study unveiled a rich and nuanced epigenetic landscape, shedding light on distinct chromatin states associated with heterochromatin formation and gene silencing. Moreover, we elucidated the intricate interplay between these chromatin states and the overall topology of the genome. Employing a genetic approach, we delved into the role of the histone modification H3K9ac in genome topology. Notably, our investigation revealed that the ectopic deposition of this chromatin mark triggered a reorganization of the 3D chromatin structure, defining different TAD-like borders. Our work emphasizes the critical role of H3K9ac in shaping the topology of the tomato genome, providing valuable insights into the epigenetic landscape of this agriculturally significant crop species.

Regulation of 3D genome organization during T cell activation

Within the three-dimensional (3D) nuclear space, the genome organizes into a series of orderly structures that impose important influences on gene regulation. T lymphocytes, crucial players in adaptive immune responses, undergo intricate transcriptional remodeling upon activation, leading to differentiation into specific effector and memory T cell subsets. Recent evidence suggests that T cell activation is accompanied by dynamic changes in genome architecture at multiple levels, providing a unique biological context to explore the functional relevance and molecular mechanisms of 3D genome organization. Here, we summarize recent advances that link the reorganization of genome architecture to the remodeling of transcriptional programs and conversion of cell fates during T cell activation and differentiation. We further discuss how various chromatin architecture regulators, including CCCTC-binding factor and several transcription factors, collectively modulate the genome architecture during this process.

Enhancing Hi-C contact matrices for loop detection with Capricorn: a multiview diffusion model

MOTIVATION

High-resolution Hi-C contact matrices reveal the detailed three-dimensional architecture of the genome, but high-coverage experimental Hi-C data are expensive to generate. Simultaneously, chromatin structure analyses struggle with extremely sparse contact matrices. To address this problem, computational methods to enhance low-coverage contact matrices have been developed, but existing methods are largely based on resolution enhancement methods for natural images and hence often employ models that do not distinguish between biologically meaningful contacts, such as loops and other stochastic contacts.

RESULTS

We present Capricorn, a machine learning model for Hi-C resolution enhancement that incorporates small-scale chromatin features as additional views of the input Hi-C contact matrix and leverages a diffusion probability model backbone to generate a high-coverage matrix. We show that Capricorn outperforms the state of the art in a cross-cell-line setting, improving on existing methods by 17% in mean squared error and 26% in F1 score for chromatin loop identification from the generated high-coverage data. We also demonstrate that Capricorn performs well in the cross-chromosome setting and cross-chromosome, cross-cell-line setting, improving the downstream loop F1 score by 14% relative to existing methods. We further show that our multiview idea can also be used to improve several existing methods, HiCARN and HiCNN, indicating the wide applicability of this approach. Finally, we use DNA sequence to validate discovered loops and find that the fraction of CTCF-supported loops from Capricorn is similar to those identified from the high-coverage data. Capricorn is a powerful Hi-C resolution enhancement method that enables scientists to find chromatin features that cannot be identified in the low-coverage contact matrix.

AVAILABILITY AND IMPLEMENTATION

Implementation of Capricorn and source code for reproducing all figures in this paper are available at https://github.com/CHNFTQ/Capricorn.

Proximal telomeric decompaction due to telomere shortening drives FOXC1-dependent myocardial senescence

Telomeres, TTAGGGn DNA repeat sequences located at the ends of eukaryotic chromosomes, play a pivotal role in aging and are targets of DNA damage response. Although we and others have demonstrated presence of short telomeres in genetic cardiomyopathic and heart failure cardiomyocytes, little is known about the role of telomere lengths in cardiomyocyte. Here, we demonstrate that in heart failure patient cardiomyocytes, telomeres are shortened compared to healthy controls. We generated isogenic human induced pluripotent stem cell derived cardiomyocytes (hiPSC-CMs) with short telomeres (sTL-CMs) and normal telomeres (nTL-CMs) as model. Compared to nTL-CMs, short telomeres result in cardiac dysfunction and expression of senescent markers. Using Hi-C and RNASeq, we observe that short telomeres induced TAD insulation decrease near telomeric ends and this correlated with a transcription upregulation in sTL-CMs. FOXC1, a key transcription factor involved in early cardiogenesis, was upregulated in sTL-CMs and its protein levels were negatively correlated with telomere lengths in heart failure patients. Overexpression of FOXC1 induced hiPSC-CM aging, mitochondrial and contractile dysfunction; knockdown of FOXC1 rescued these phenotypes. Overall, the work presented demonstrate that increased chromatin accessibility due to telomere shortening resulted in the induction of FOXC1-dependent expression network responsible for contractile dysfunction and myocardial senescence.

Extensive long-range polycomb interactions and weak compartmentalization are hallmarks of human neuronal 3D genome

Chromatin architecture regulates gene expression and shapes cellular identity, particularly in neuronal cells. Specifically, polycomb group (PcG) proteins enable establishment and maintenance of neuronal cell type by reorganizing chromatin into repressive domains that limit the expression of fate-determining genes and sustain distinct gene expression patterns in neurons. Here, we map the 3D genome architecture in neuronal and non-neuronal cells isolated from the Wernicke's area of four human brains and comprehensively analyze neuron-specific aspects of chromatin organization. We find that genome segregation into active and inactive compartments is greatly reduced in neurons compared to other brain cells. Furthermore, neuronal Hi-C maps reveal strong long-range interactions, forming a specific network of PcG-mediated contacts in neurons that is nearly absent in other brain cells. These interacting loci contain developmental transcription factors with repressed expression in neurons and other mature brain cells. But only in neurons, they are rich in bivalent promoters occupied by H3K4me3 histone modification together with H3K27me3, which points to a possible functional role of PcG contacts in neurons. Importantly, other layers of chromatin organization also exhibit a distinct structure in neurons, characterized by an increase in short-range interactions and a decrease in long-range ones.

Intertwining roles of R-loops and G-quadruplexes in DNA repair, transcription and genome organization

R-loops are three-stranded nucleic acid structures that are abundant and widespread across the genome and that have important physiological roles in many nuclear processes. Their accumulation is observed in cancers and neurodegenerative disorders. Recent studies have implicated a function for R-loops and G-quadruplex (G4) structures, which can form on the displaced single strand of R-loops, in three-dimensional genome organization in both physiological and pathological contexts. Here we discuss the interconnected functions of DNA:RNA hybrids and G4s within R-loops, their impact on DNA repair and gene regulatory networks, and their emerging roles in genome organization during development and disease.

Reprogramming of 3D genome structure underlying HSPC development in zebrafish

BACKGROUND

Development of hematopoietic stem and progenitor cells (HSPC) is a multi-staged complex process that conserved between zebrafish and mammals. Understanding the mechanism underlying HSPC development is a holy grail of hematopoietic biology, which is helpful for HSPC clinical application. Chromatin conformation plays important roles in transcriptional regulation and cell fate decision; however, its dynamic and role in HSPC development is poorly investigated.

METHODS

We performed chromatin structure and multi-omics dissection across different stages of HSPC developmental trajectory in zebrafish for the first time, including Hi-C, RNA-seq, ATAC-seq, H3K4me3 and H3K27ac ChIP-seq.

RESULTS

The chromatin organization of zebrafish HSPC resemble mammalian cells with similar hierarchical structure. We revealed the multi-scale reorganization of chromatin structure and its influence on transcriptional regulation and transition of cell fate during HSPC development. Nascent HSPC is featured by loose conformation with obscure structure at all layers. Notably, PU.1 was identified as a potential factor mediating formation of promoter-involved loops and regulating gene expression of HSPC.

CONCLUSIONS

Our results provided a global view of chromatin structure dynamics associated with development of zebrafish HSPC and discovered key transcription factors involved in HSPC chromatin interactions, which will provide new insights into the epigenetic regulatory mechanisms underlying vertebrate HSPC fate decision.

Dynamic control of gene expression by ISGF3 and IRF1 during IFNβ and IFNγ signaling

Type I interferons (IFN-I, including IFNβ) and IFNγ produce overlapping, yet clearly distinct immunological activities. Recent data show that the distinctness of global transcriptional responses to the two IFN types is not apparent when comparing their immediate effects. By analyzing nascent transcripts induced by IFN-I or IFNγ over a period of 48 h, we now show that the distinctiveness of the transcriptomes emerges over time and is based on differential employment of the ISGF3 complex as well as of the second-tier transcription factor IRF1. The distinct transcriptional properties of ISGF3 and IRF1 correspond with a largely diverse nuclear protein interactome. Mechanistically, we describe the specific input of ISGF3 and IRF1 into enhancer activation and the regulation of chromatin accessibility at interferon-stimulated genes (ISG). We further report differences between the IFN types in altering RNA polymerase II pausing at ISG 5' ends. Our data provide insight how transcriptional regulators create immunological identities of IFN-I and IFNγ.

Multiscale Bayesian simulations reveal functional chromatin condensation of gene loci

Chromatin, the complex assembly of DNA and associated proteins, plays a pivotal role in orchestrating various genomic functions. To aid our understanding of the principles underlying chromatin organization, we introduce Hi-C metainference, a Bayesian approach that integrates Hi-C contact frequencies into multiscale prior models of chromatin. This approach combines both bottom-up (the physics-based prior) and top-down (the data-driven posterior) strategies to characterize the 3D organization of a target genomic locus. We first demonstrate the capability of this method to accurately reconstruct the structural ensemble and the dynamics of a system from contact information. We then apply the approach to investigate the Sox2, Pou5f1, and Nanog loci of mouse embryonic stem cells using a bottom-up chromatin model at 1 kb resolution. We observe that the studied loci are conformationally heterogeneous and organized as crumpled globules, favoring contacts between distant enhancers and promoters. Using nucleosome-resolution simulations, we then reveal how the Nanog gene is functionally organized across the multiple scales of chromatin. At the local level, we identify diverse tetranucleosome folding motifs with a characteristic distribution along the genome, predominantly open at cis-regulatory elements and compact in between. At the larger scale, we find that enhancer-promoter contacts are driven by the transient condensation of chromatin into compact domains stabilized by extensive internucleosome interactions. Overall, this work highlights the condensed, but dynamic nature of chromatin in vivo, contributing to a deeper understanding of gene structure-function relationships.

Evolutionary insights into the organization of chromatin structure and landscape of transcriptional regulation in plants

Development of complex traits necessitates the functioning and coordination of intricate regulatory networks involving multiple genes. Understanding 3D chromatin structure can facilitate insight into the regulation of gene expression by regulatory elements. This potential, of visualizing the role of chromatin organization in the evolution and function of regulatory elements, remains largely unexplored. Here, we describe new perspectives that arise from the dual considerations of sequence variation of regulatory elements and chromatin structure, with a special focus on whole-genome doubling or polyploidy. We underscore the significance of hierarchical chromatin organization in gene regulation during evolution. In addition, we describe strategies for exploring chromatin organization in future investigations of regulatory evolution in plants, enabling insights into the evolutionary influence of regulatory elements on gene expression and, hence, phenotypes.

Multiscale 3D genome organization underlies duck fatty liver with no adipose inflammation or serious injury

Non-alcoholic fatty liver disease (NAFLD) is the most common form of chronic liver disease. Little is known about how gene expression and chromatin structure are regulated in NAFLD due to lack of suitable model. Ducks naturally develop fatty liver similar to serious human non-alcoholic fatty liver (NAFL) without adipose inflammation and liver fibrosis, thus serves as a good model for investigating molecular mechanisms of adipose metabolism and anti-inflammation. Here, we constructed a NAFLD model without adipose inflammation and liver fibrosis in ducks. By performing dynamic pathological and transcriptomic analyses, we identified critical genes involving in regulation of the NF-κB and MHCII signaling, which usually lead to adipose inflammation and liver fibrosis. We further generated dynamic three-dimensional chromatin maps during liver fatty formation and recovery. This showed that ducks enlarged hepatocyte cell nuclei to reduce inter-chromosomal interaction, decompress chromatin structure, and alter strength of intra-TAD and loop interactions during fatty liver formation. These changes partially contributed to the tight control the NF-κB and the MHCII signaling. Our analysis uncovers duck chromatin reorganization might be advantageous to maintain liver regenerative capacity and reduce adipose inflammation. These findings shed light on new strategies for NAFLD control.

Chromatin organization in myelodysplastic syndrome

Disordered chromatin organization has emerged as a new aspect of the pathogenesis of myelodysplastic syndrome (MDS). Characterized by lineage dysplasia and a high transformation rate to acute myeloid leukemia (AML), the genetic determinant of MDS is thought to be the main driver of the disease's progression. Among the recurrently mutated pathways, alterations in chromatin organization, such as the cohesin complex, have a profound impact on hematopoietic stem cell (HSC) function and lineage commitment. The cohesin complex is a ring-like structure comprised of structural maintenance of chromosomes (SMC), RAD21, and STAG proteins that involve three-dimensional (3D) genome organization via loop extrusion in mammalian cells. The partial loss of the functional cohesin ring leads to altered chromatin accessibility specific to key hematopoietic transcription factors, which is thought to be the molecular mechanism of cohesin dysfunction. Currently, there are no specific targeting agents for cohesin mutant MDS/AML. Potential therapeutic strategies have been proposed based on the current understanding of cohesin mutant leukemogenesis. Here, we will review the recent advances in investigation and targeting approaches against cohesin mutant MDS/AML.

From compartments to loops: understanding the unique chromatin organization in neuronal cells

The three-dimensional organization of the genome plays a central role in the regulation of cellular functions, particularly in the human brain. This review explores the intricacies of chromatin organization, highlighting the distinct structural patterns observed between neuronal and non-neuronal brain cells. We integrate findings from recent studies to elucidate the characteristics of various levels of chromatin organization, from differential compartmentalization and topologically associating domains (TADs) to chromatin loop formation. By defining the unique chromatin landscapes of neuronal and non-neuronal brain cells, these distinct structures contribute to the regulation of gene expression specific to each cell type. In particular, we discuss potential functional implications of unique neuronal chromatin organization characteristics, such as weaker compartmentalization, neuron-specific TAD boundaries enriched with active histone marks, and an increased number of chromatin loops. Additionally, we explore the role of Polycomb group (PcG) proteins in shaping cell-type-specific chromatin patterns. This review further emphasizes the impact of variations in chromatin architecture between neuronal and non-neuronal cells on brain development and the onset of neurological disorders. It highlights the need for further research to elucidate the details of chromatin organization in the human brain in order to unravel the complexities of brain function and the genetic mechanisms underlying neurological disorders. This research will help bridge a significant gap in our comprehension of the interplay between chromatin structure and cell functions.

C2c: Predicting Micro-C from Hi-C

MOTIVATION

High-resolution Hi-C data, capable of detecting chromatin features below the level of Topologically Associating Domains (TADs), significantly enhance our understanding of gene regulation. Micro-C, a variant of Hi-C incorporating a micrococcal nuclease (MNase) digestion step to examine interactions between nucleosome pairs, has been developed to overcome the resolution limitations of Hi-C. However, Micro-C experiments pose greater technical challenges compared to Hi-C, owing to the need for precise MNase digestion control and higher-resolution sequencing. Therefore, developing computational methods to derive Micro-C data from existing Hi-C datasets could lead to better usage of a large amount of existing Hi-C data in the scientific community and cost savings.

RESULTS

We developed C2c ("high" or upper case C to "micro" or lower case c), a computational tool based on a residual neural network to learn the mapping between Hi-C and Micro-C contact matrices and then predict Micro-C contact matrices based on Hi-C contact matrices. Our evaluation results show that the predicted Micro-C contact matrices reveal more chromatin loops than the input Hi-C contact matrices, and more of the loops detected from predicted Micro-C match the promoter-enhancer interactions. Furthermore, we found that the mutual loops from real and predicted Micro-C better match the ChIA-PET data compared to Hi-C and real Micro-C loops, and the predicted Micro-C leads to more TAD-boundaries detected compared to the Hi-C data. The website URL of C2c can be found in the Data Availability Statement.

Statistical modeling and significance estimation of multi-way chromatin contacts with HyperloopFinder

Recent advances in chromatin conformation capture technologies, such as SPRITE and Pore-C, have enabled the detection of simultaneous contacts among multiple chromatin loci. This has made it possible to investigate the cooperative transcriptional regulation involving multiple genes and regulatory elements at the resolution of a single molecule. However, these technologies are unavoidably subject to the random polymer looping effect and technical biases, making it challenging to distinguish genuine regulatory relationships directly from random polymer interactions. Here, we present HyperloopFinder, a method for identifying regulatory multi-way chromatin contacts (hyperloops) by jointly modeling the random polymer looping effect and technical biases to estimate the statistical significance of multi-way contacts. The results show that our model can accurately estimate the expected interaction frequency of multi-way contacts based on the distance distribution of pairwise contacts, revealing that most multi-way contacts can be formed by randomly linking the pairwise contacts adjacent to each other. Moreover, we observed the spatial colocalization of the interaction sites of hyperloops from image-based data. Our results also revealed that hyperloops can function as scaffolds for the cooperation among multiple genes and regulatory elements. In summary, our work contributes novel insights into higher-order chromatin structures and functions and has the potential to enhance our understanding of transcriptional regulation and other cellular processes.

Drought response revealed by chromatin organization variation and transcriptional regulation in cotton

BACKGROUND

Cotton is a major world cash crop and an important source of natural fiber, oil, and protein. Drought stress is becoming a restrictive factor affecting cotton production. To facilitate the development of drought-tolerant cotton varieties, it is necessary to study the molecular mechanism of drought stress response by exploring key drought-resistant genes and related regulatory factors.

RESULTS

In this study, two cotton varieties, ZY007 (drought-sensitive) and ZY168 (drought-tolerant), showing obvious phenotypic differences under drought stress, were selected. A total of 25,898 drought-induced genes were identified, exhibiting significant enrichment in pathways related to plant stress responses. Under drought induction, At subgenome expression bias was observed at the whole-genome level, which may be due to stronger inhibition of Dt subgenome expression. A gene co-expression module that was significantly associated with drought resistance was identified. About 90% of topologically associating domain (TAD) boundaries were stable, and 6613 TAD variation events were identified between the two varieties under drought. We identified 92 genes in ZY007 and 98 in ZY168 related to chromatin 3D structural variation and induced by drought stress. These genes are closely linked to the cotton response to drought stress through canonical hormone-responsive pathways, modulation of kinase and phosphatase activities, facilitation of calcium ion transport, and other related molecular mechanisms.

CONCLUSIONS

These results lay a foundation for elucidating the molecular mechanism of the cotton drought response and provide important regulatory locus and gene resources for the future molecular breeding of drought-resistant cotton varieties.

CD-Loop: a chromatin loop detection method based on the diffusion model

Motivation

In recent years, there have been significant advances in various chromatin conformation capture techniques, and annotating the topological structure from Hi-C contact maps has become crucial for studying the three-dimensional structure of chromosomes. However, the structure and function of chromatin loops are highly dynamic and diverse, influenced by multiple factors. Therefore, obtaining the three-dimensional structure of the genome remains a challenging task. Among many chromatin loop prediction methods, it is difficult to fully extract features from the contact map and make accurate predictions at low sequencing depths.

Results

In this study, we put forward a deep learning framework based on the diffusion model called CD-Loop for predicting accurate chromatin loops. First, by pre-training the input data, we obtain prior probabilities for predicting the classification of the Hi-C contact map. Then, by combining the denoising process based on the diffusion model and the prior probability obtained by pre-training, candidate loops were predicted from the input Hi-C contact map. Finally, CD-Loop uses a density-based clustering algorithm to cluster the candidate chromatin loops and predict the final chromatin loops. We compared CD-Loop with the currently popular methods, such as Peakachu, Chromosight, and Mustache, and found that in different cell types, species, and sequencing depths, CD-Loop outperforms other methods in loop annotation. We conclude that CD-Loop can accurately predict chromatin loops and reveal cell-type specificity. The code is available at https://github.com/wangyang199897/CD-Loop.

An RNA-dependent and phase-separated active subnuclear compartment safeguards repressive chromatin domains

The nucleus is composed of functionally distinct membraneless compartments that undergo phase separation (PS). However, whether different subnuclear compartments are connected remains elusive. We identified a type of nuclear body with PS features composed of BAZ2A that associates with active chromatin. BAZ2A bodies depend on RNA transcription and BAZ2A non-disordered RNA-binding TAM domain. Although BAZ2A and H3K27me3 occupancies anticorrelate in the linear genome, in the nuclear space, BAZ2A bodies contact H3K27me3 bodies. BAZ2A-body disruption promotes BAZ2A invasion into H3K27me3 domains, causing H3K27me3-body loss and gene upregulation. Weak BAZ2A-RNA interactions, such as with nascent transcripts, promote BAZ2A bodies, whereas the strong binder long non-coding RNA (lncRNA) Malat1 impairs them while mediating BAZ2A association to chromatin at nuclear speckles. In addition to unraveling a direct connection between nuclear active and repressive compartments through PS mechanisms, the results also showed that the strength of RNA-protein interactions regulates this process, contributing to nuclear organization and the regulation of chromatin and gene expression.

Unveiling Alterations of Epigenetic Modifications and Chromatin Architecture Leading to Lipid Metabolic Reprogramming during the Evolutionary Trastuzumab Adaptation of HER2-Positive Breast Cancer

Secondary trastuzumab resistance represents an evolutionary adaptation of HER2-positive breast cancer during anti-HER2 treatment. Most current studies have tended to prioritize HER2 and its associated signaling pathways, often overlooking broader but seemingly less relevant cellular processes, along with their associated genetic and epigenetic mechanisms. Here, transcriptome data is not only characterized but also examined epigenomic and 3D genome architecture information in both trastuzumab-sensitive and secondary-resistant breast cancer cells. The findings reveal that the global metabolic reprogramming associated with trastuzumab resistance may stem from genome-wide alterations in both histone modifications and chromatin structure. Specifically, the transcriptional activities of key genes involved in lipid metabolism appear to be regulated by variant promoter H3K27me3 and H3K4me3 modifications, as well as promoter-enhancer interactions. These discoveries offer valuable insights into how cancer cells adapt to anti-tumor drugs and have the potential to impact future diagnostic and treatment strategies.

Three-dimensional chromatin landscapes in somatotroph tumour

BACKGROUND

The three-dimensional (3D) genome architecture plays a critical role inregulating gene expression. However, the specific alterations in thisarchitecture within somatotroph tumors and their implications for gene expression remain largely unexplored.

METHODS

We employed Hi-C and RNA-seq analyses to compare the 3D genomic structures of somatotroph tumors with normal pituitary tissue. This comprehensive approachenabled the characterization of A/B compartments, topologically associateddomains (TADs), and chromatin loops, integrating these with gene expression patterns.

RESULTS

We observed a decrease in both the frequency of chromosomal interactions andthe size of TADs in tumor tissue compared to normal tissue. Conversely, the number of TADs and chromatin loops was found to be increased in tumors. Integrated analysis of Hi-C and RNA-seq data demonstrated that changes inhigher-order chromat in structure were associated with alterations in gene expression. Specifically, genes in A compartments showed higher density and increased expression relative to those in B compartments. Moreover, the weakand enhanced insulation boundaries were identified, and the associated genes were enriched in the Wnt/β-Catenin signaling pathway. We identified the gainedand lost loops in tumor and integrated these differences with transcriptional changes to examine the functional relevance of the identified loops. Notably, we observed an enhanced insulation boundary and a greater number of loops in the TCF7L2 gene region within tumors, which was accompanied by an upregulation of TCF7L2 expression. Subsequently, TCF7L2 expression was confirmed through qRT-PCR, and upregulated TCF7L2 prompted cell proliferation and growth hormone (GH) secretion in vitro.

CONCLUSION

Our results provide comprehensive 3D chromatin architecture maps of somatotroph tumors and offer a valuable resource for furthering the understanding of the underlying biology and mechanisms of gene expression regulation.

Genome structural dynamics: insights from Gaussian network analysis of Hi-C data

Characterization of the spatiotemporal properties of the chromatin is essential to gaining insights into the physical bases of gene co-expression, transcriptional regulation and epigenetic modifications. The Gaussian network model (GNM) has proven in recent work to serve as a useful tool for modeling chromatin structural dynamics, using as input high-throughput chromosome conformation capture data. We focus here on the exploration of the collective dynamics of chromosomal structures at hierarchical levels of resolution, from single gene loci to topologically associating domains or entire chromosomes. The GNM permits us to identify long-range interactions between gene loci, shedding light on the role of cross-correlations between distal regions of the chromosomes in regulating gene expression. Notably, GNM analysis performed across diverse cell lines highlights the conservation of the global/cooperative movements of the chromatin across different types of cells. Variations driven by localized couplings between genomic loci, on the other hand, underlie cell differentiation, underscoring the significance of the four-dimensional properties of the genome in defining cellular identity. Finally, we demonstrate the close relation between the cell type-dependent mobility profiles of gene loci and their gene expression patterns, providing a clear demonstration of the role of chromosomal 4D features in defining cell-specific differential expression of genes.

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

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

Droplet Hi-C for Fast and Scalable Profiling of Chromatin Architecture in Single Cells

Comprehensive analysis of chromatin architecture is crucial for understanding the gene regulatory programs during development and in disease pathogenesis, yet current methods often inadequately address the unique challenges presented by analysis of heterogeneous tissue samples. Here, we introduce Droplet Hi-C, which employs a commercial microfluidic device for high-throughput, single-cell chromatin conformation profiling in droplets. Using Droplet Hi-C, we mapped the chromatin architecture at single-cell resolution from the mouse cortex and analyzed gene regulatory programs in major cortical cell types. Additionally, we used this technique to detect copy number variation (CNV), structural variations (SVs) and extrachromosomal DNA (ecDNA) in cancer cells, revealing clonal dynamics and other oncogenic events during treatment. We further refined this technique to allow for joint profiling of chromatin architecture and transcriptome in single cells, facilitating a more comprehensive exploration of the links between chromatin architecture and gene expression in both normal tissues and tumors. Thus, Droplet Hi-C not only addresses critical gaps in chromatin analysis of heterogeneous tissues but also emerges as a versatile tool enhancing our understanding of gene regulation in health and disease.

Analysis of long-range chromatin contacts, compartments and looping between mouse embryonic stem cells, lens epithelium and lens fibers

BACKGROUND

Nuclear organization of interphase chromosomes involves individual chromosome territories, "open" and "closed" chromatin compartments, topologically associated domains (TADs) and chromatin loops. The DNA- and RNA-binding transcription factor CTCF together with the cohesin complex serve as major organizers of chromatin architecture. Cellular differentiation is driven by temporally and spatially coordinated gene expression that requires chromatin changes of individual loci of various complexities. Lens differentiation represents an advantageous system to probe transcriptional mechanisms underlying tissue-specific gene expression including high transcriptional outputs of individual crystallin genes until the mature lens fiber cells degrade their nuclei.

RESULTS

Chromatin organization between mouse embryonic stem (ES) cells, newborn (P0.5) lens epithelium and fiber cells were analyzed using Hi-C. Localization of CTCF in both lens chromatins was determined by ChIP-seq and compared with ES cells. Quantitative analyses show major differences between number and size of TADs and chromatin loop size between these three cell types. In depth analyses show similarities between lens samples exemplified by overlaps between compartments A and B. Lens epithelium-specific CTCF peaks are found in mostly methylated genomic regions while lens fiber-specific and shared peaks occur mostly within unmethylated DNA regions. Major differences in TADs and loops are illustrated at the ~ 500 kb Pax6 locus, encoding the critical lens regulatory transcription factor and within a larger ~ 15 Mb WAGR locus, containing Pax6 and other loci linked to human congenital diseases. Lens and ES cell Hi-C data (TADs and loops) together with ATAC-seq, CTCF, H3K27ac, H3K27me3 and ENCODE cis-regulatory sites are shown in detail for the Pax6, Sox1 and Hif1a loci, multiple crystallin genes and other important loci required for lens morphogenesis. The majority of crystallin loci are marked by unexpectedly high CTCF-binding across their transcribed regions.

CONCLUSIONS

Our study has generated the first data on 3-dimensional (3D) nuclear organization in lens epithelium and lens fibers and directly compared these data with ES cells. These findings generate novel insights into lens-specific transcriptional gene control, open new research avenues to study transcriptional condensates in lens fiber cells, and enable studies of non-coding genetic variants linked to cataract and other lens and ocular abnormalities.

Integrative analysis of transcriptomic and epigenomic data reveals distinct patterns for developmental and housekeeping gene regulation

BACKGROUND

Regulation of transcription is central to the emergence of new cell types during development, and it often involves activation of genes via proximal and distal regulatory regions. The activity of regulatory elements is determined by transcription factors (TFs) and epigenetic marks, but despite extensive mapping of such patterns, the extraction of regulatory principles remains challenging.

RESULTS

Here we study differentially and similarly expressed genes along with their associated epigenomic profiles, chromatin accessibility and DNA methylation, during lineage specification at gastrulation in mice. Comparison of the three lineages allows us to identify genomic and epigenomic features that distinguish the two classes of genes. We show that differentially expressed genes are primarily regulated by distal elements, while similarly expressed genes are controlled by proximal housekeeping regulatory programs. Differentially expressed genes are relatively isolated within topologically associated domains, while similarly expressed genes tend to be located in gene clusters. Transcription of differentially expressed genes is associated with differentially open chromatin at distal elements including enhancers, while that of similarly expressed genes is associated with ubiquitously accessible chromatin at promoters.

CONCLUSION

Based on these associations of (linearly) distal genes' transcription start sites (TSSs) and putative enhancers for developmental genes, our findings allow us to link putative enhancers to their target promoters and to infer lineage-specific repertoires of putative driver transcription factors, within which we define subgroups of pioneers and co-operators.

Cohesin regulation and roles in chromosome structure and function

Chromosome structure regulates DNA-templated processes such as transcription of genes. Dynamic changes to chromosome structure occur during development and in disease contexts. The cohesin complex is a molecular motor that regulates chromosome structure by generating DNA loops that bring two distal genomic sites into close spatial proximity. There are many open questions regarding the formation and dissolution of DNA loops, as well as the role(s) of DNA loops in regulating transcription of the interphase genome. This review focuses on recent discoveries that provide molecular insights into the role of cohesin and chromosome structure in gene transcription during development and disease.

A systematic analyses of different bioinformatics pipelines for genomic data and its impact on deep learning models for chromatin loop prediction

Genomic data analysis has witnessed a surge in complexity and volume, primarily driven by the advent of high-throughput technologies. In particular, studying chromatin loops and structures has become pivotal in understanding gene regulation and genome organization. This systematic investigation explores the realm of specialized bioinformatics pipelines designed specifically for the analysis of chromatin loops and structures. Our investigation incorporates two protein (CTCF and Cohesin) factor-specific loop interaction datasets from six distinct pipelines, amassing a comprehensive collection of 36 diverse datasets. Through a meticulous review of existing literature, we offer a holistic perspective on the methodologies, tools and algorithms underpinning the analysis of this multifaceted genomic feature. We illuminate the vast array of approaches deployed, encompassing pivotal aspects such as data preparation pipeline, preprocessing, statistical features and modelling techniques. Beyond this, we rigorously assess the strengths and limitations inherent in these bioinformatics pipelines, shedding light on the interplay between data quality and the performance of deep learning models, ultimately advancing our comprehension of genomic intricacies.

Select EZH2 inhibitors enhance viral mimicry effects of DNMT inhibition through a mechanism involving NFAT:AP-1 signaling

DNA methyltransferase inhibitor (DNMTi) efficacy in solid tumors is limited. Colon cancer cells exposed to DNMTi accumulate lysine-27 trimethylation on histone H3 (H3K27me3). We propose this Enhancer of Zeste Homolog 2 (EZH2)-dependent repressive modification limits DNMTi efficacy. Here, we show that low-dose DNMTi treatment sensitizes colon cancer cells to select EZH2 inhibitors (EZH2is). Integrative epigenomic analysis reveals that DNMTi-induced H3K27me3 accumulates at genomic regions poised with EZH2. Notably, combined EZH2i and DNMTi alters the epigenomic landscape to transcriptionally up-regulate the calcium-induced nuclear factor of activated T cells (NFAT):activating protein 1 (AP-1) signaling pathway. Blocking this pathway limits transcriptional activating effects of these drugs, including transposable element and innate immune response gene expression involved in viral defense. Analysis of primary human colon cancer specimens reveals positive correlations between DNMTi-, innate immune response-, and calcium signaling-associated transcription profiles. Collectively, we show that compensatory EZH2 activity limits DNMTi efficacy in colon cancer and link NFAT:AP-1 signaling to epigenetic therapy-induced viral mimicry.

Post-GWAS functional analyses of CNTNAP5 suggests its role in glaucomatous neurodegeneration

Primary angle closure glaucoma (PACG) affects more than 20 million people worldwide, with an increased prevalence in south-east Asia. In a prior haplotype-based GWAS, we identified a novel CNTNAP5 genic region, significantly associated with PACG. In the current study, we have extended our perception of CNTNAP5 involvement in glaucomatous neurodegeneration in a zebrafish model, through investigating phenotypic consequences pertinent to retinal degeneration upon knockdown of cntnap5 by translation-blocking morpholinos. While cntnap5 knockdown was successfully validated using an antibody, immunofluorescence followed by western blot analyses in cntnap5-morphant (MO) zebrafish revealed increased expression of acetylated tubulin indicative of perturbed cytoarchitecture of retinal layers. Moreover, significant loss of Nissl substance is observed in the neuro-retinal layers of cntnap5-MO zebrafish eye, indicating neurodegeneration. Additionally, in spontaneous movement behavioural analysis, cntnap5-MO zebrafish have a significantly lower average distance traversed in light phase compared to mismatch-controls, whereas no significant difference was observed in the dark phase, corroborating with vision loss in the cntnap5-MO zebrafish. This study provides the first direct functional evidence of a putative role of CNTNAP5 in visual neurodegeneration.

Identifying TAD-like domains on single-cell Hi-C data by graph embedding and changepoint detection

MOTIVATION

Topologically associating domains (TADs) are fundamental building blocks of 3D genome. TAD-like domains in single cells are regarded as the underlying genesis of TADs discovered in bulk cells. Understanding the organization of TAD-like domains helps to get deeper insights into their regulatory functions. Unfortunately, it remains a challenge to identify TAD-like domains on single-cell Hi-C data due to its ultra-sparsity.

RESULTS

We propose scKTLD, an in silico tool for the identification of TAD-like domains on single-cell Hi-C data. It takes Hi-C contact matrix as the adjacency matrix for a graph, embeds the graph structures into a low-dimensional space with the help of sparse matrix factorization followed by spectral propagation, and the TAD-like domains can be identified using a kernel-based changepoint detection in the embedding space. The results tell that our scKTLD is superior to the other methods on the sparse contact matrices, including downsampled bulk Hi-C data as well as simulated and experimental single-cell Hi-C data. Besides, we demonstrated the conservation of TAD-like domain boundaries at single-cell level apart from heterogeneity within and across cell types, and found that the boundaries with higher frequency across single cells are more enriched for architectural proteins and chromatin marks, and they preferentially occur at TAD boundaries in bulk cells, especially at those with higher hierarchical levels.

AVAILABILITY AND IMPLEMENTATION

scKTLD is freely available at https://github.com/lhqxinghun/scKTLD.

SpecLoop predicts cell type-specific chromatin loop via transcription factor cooperation

Cell-type-Specific Chromatin Loops (CSCLs) are crucial for gene regulation and cell fate determination. However, the mechanisms governing their establishment remain elusive. Here, we present SpecLoop, a network regularization-based machine learning framework, to investigate the role of transcription factors (TFs) cooperation in CSCL formation. SpecLoop integrates multi-omics data, including gene expression, chromatin accessibility, sequence, protein-protein interaction, and TF binding motif data, to predict CSCLs and identify TF cooperations. Using high resolution Hi-C data as the gold standard, SpecLoop accurately predicts CSCL in GM12878, IMR90, HeLa-S3, K562, HUVEC, HMEC, and NHEK seven cell types, with the AUROC values ranging from 0.8645 to 0.9852 and AUPR values ranging from 0.8654 to 0.9734. Notably SpecLoop demonstrates improved accuracy in predicting long-distance CSCLs and identifies TF complexes with strong predictive ability. Our study systematically explores the TFs and TF pairs associated with CSCL through effective integration of diverse omics data. SpecLoop is freely available at https://github.com/AMSSwanglab/SpecLoop.

The impact of DNA methylation on CTCF-mediated 3D genome organization

Cytosine DNA methylation is a highly conserved epigenetic mark in eukaryotes. Although the role of DNA methylation at gene promoters and repetitive elements has been extensively studied, the function of DNA methylation in other genomic contexts remains less clear. In the nucleus of mammalian cells, the genome is spatially organized at different levels, and strongly influences myriad genomic processes. There are a number of factors that regulate the three-dimensional (3D) organization of the genome, with the CTCF insulator protein being among the most well-characterized. Pertinently, CTCF binding has been reported as being DNA methylation-sensitive in certain contexts, perhaps most notably in the process of genomic imprinting. Therefore, it stands to reason that DNA methylation may play a broader role in the regulation of chromatin architecture. Here we summarize the current understanding that is relevant to both the mammalian DNA methylation and chromatin architecture fields and attempt to assess the extent to which DNA methylation impacts the folding of the genome. The focus is in early embryonic development and cellular transitions when the epigenome is in flux, but we also describe insights from pathological contexts, such as cancer, in which the epigenome and 3D genome organization are misregulated.

Non-coding DNA variants for risk in lupus

Systemic Lupus Erythematosus (SLE) is a multifactorial autoimmune disease that arises from a dynamic interplay between genetics and environmental triggers. The advent of sophisticated genomics technology has catalyzed a shift in our understanding of disease etiology, spotlighting the pivotal role of non-coding DNA variants in SLE pathogenesis. In this review, we present a comprehensive examination of the non-coding variants associated with SLE, shedding light on their role in influencing disease risk and progression. We discuss the latest methodological advancements that have been instrumental in the identification and functional characterization of these genomic elements, with a special focus on the transformative power of CRISPR-based gene-editing technologies. Additionally, the review probes into the therapeutic opportunities that arise from modulating non-coding regions associated with SLE. Through an exploration of the complex network of non-coding DNA, this review aspires to decode the genetic puzzle of SLE and set the stage for groundbreaking gene-based therapeutic interventions and the advancement of precision medicine strategies tailored to SLE management.

Alteration of chromosome structure impacts gene expressions implicated in pancreatic ductal adenocarcinoma cells

BACKGROUND

Pancreatic ductal adenocarcinoma (PDAC) is a lethal malignancy with a five-year survival rate of approximately 10%. Genetic mutations are pivotal drivers in PDAC pathogenesis, but recent investigations also revealed the involvement of non-genetic alterations in the disease development. In this study, we undertook a multi-omics approach, encompassing ATAC-seq, RNA-seq, ChIP-seq, and Hi-C methodologies, to dissect gene expression alterations arising from changes in chromosome accessibility and chromatin three-dimensional interactions in PDAC.

RESULTS

Our findings indicate that chromosomal structural alterations can lead to abnormal expressions on key genes during PDAC development. Notably, overexpression of oncogenes FGFR2, FOXA2, CYP2R1, and CPOX can be attributed to the augmentation of promoter accessibility, coupled with long-range interactions with distal elements. Additionally, our findings indicate that chromosomal structural alterations caused by genomic instability can lead to abnormal expressions in PDACs. As an example, by analyzing chromosomal changes, we identified a putative oncogenic gene, LPAR1, which shows upregulated expression in both PDAC cell lines and clinical samples. The overexpression is correlated with alterations in LPAR1-associated 3D genome structure and chromatin state. We further demonstrated that high LPAR1 activity is required for enhanced PDAC cell migration in vitro.

CONCLUSIONS

Collectively, our findings reveal that the chromosomal conformational alterations, in addition to the well-known genetic mutations, are critical for PDAC tumorigenesis.

Posterior inference of Hi-C contact frequency through sampling

Hi-C is one of the most widely used approaches to study three-dimensional genome conformations. Contacts captured by a Hi-C experiment are represented in a contact frequency matrix. Due to the limited sequencing depth and other factors, Hi-C contact frequency matrices are only approximations of the true interaction frequencies and are further reported without any quantification of uncertainty. Hence, downstream analyses based on Hi-C contact maps (e.g., TAD and loop annotation) are themselves point estimations. Here, we present the Hi-C interaction frequency sampler (HiCSampler) that reliably infers the posterior distribution of the interaction frequency for a given Hi-C contact map by exploiting dependencies between neighboring loci. Posterior predictive checks demonstrate that HiCSampler can infer highly predictive chromosomal interaction frequency. Summary statistics calculated by HiCSampler provide a measurement of the uncertainty for Hi-C experiments, and samples inferred by HiCSampler are ready for use by most downstream analysis tools off the shelf and permit uncertainty measurements in these analyses without modifications.

Integrating Dynamic 3D Chromatin Architecture and Gene Expression Alterations Reveal Heterosis in Brassica rapa

Heterosis plays a significant role in enhancing variety, boosting yield, and raising economic value in crops, but the molecular mechanism is still unclear. We analyzed the transcriptomes and 3D genomes of a hybrid (F1) and its parents (w30 and 082). The analysis of the expression revealed a total of 485 specially expressed genes (SEGs), 173 differentially expressed genes (DEGs) above the parental expression level, more actively expressed genes, and up-regulated DEGs in the F1. Further study revealed that the DEGs detected in the F1 and its parents were mainly involved in the response to auxin, plant hormone signal transduction, DNA metabolic process, purine metabolism, starch, and sucrose metabolism, which suggested that these biological processes may play a crucial role in the heterosis of Brassica rapa. The analysis of 3D genome data revealed that hybrid F1 plants tend to contain more transcriptionally active A chromatin compartments after hybridization. Supplementaryly, the F1 had a smaller TAD (topologically associated domain) genome length, but the number was the highest, and the expression change in activated TAD was higher than that of repressed TAD. More specific TAD boundaries were detected between the parents and F1. Subsequently, 140 DEGs with genomic structural variants were selected as potential candidate genes. We found two DEGs with consistent expression changes in A/B compartments and TADs. Our findings suggested that genomic structural variants, such as TADs and A/B chromatin compartments, may affect gene expression and contribute to heterosis in Brassica rapa. This study provides further insight into the molecular mechanism of heterosis in Brassica rapa.

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

BACKGROUND

Three-dimensional (3D) chromatin architecture frequently altered in cancer. However, its changes during the pathogenesis of hepatocellular carcinoma (HCC) remained elusive.

METHODS

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

CONCLUSION

This study demonstrates the crucial role of chromosomal structure variation in HCC oncogenesis and potential novel biomarkers of HCC, laying a foundation for cancer precision medicine development.

Exploring noncoding variants in genetic diseases: from detection to functional insights

Previous studies on genetic diseases predominantly focused on protein-coding variations, overlooking the vast noncoding regions in the human genome. The development of high-throughput sequencing technologies and functional genomics tools has enabled the systematic identification of functional noncoding variants. These variants can impact gene expression, regulation, and chromatin conformation, thereby contributing to disease pathogenesis. Understanding the mechanisms that underlie the impact of noncoding variants on genetic diseases is indispensable for the development of precisely targeted therapies and the implementation of personalized medicine strategies. The intricacies of noncoding regions introduce a multitude of challenges and research opportunities. In this review, we introduce a spectrum of noncoding variants involved in genetic diseases, along with research strategies and advanced technologies for their precise identification and in-depth understanding of the complexity of the noncoding genome. We will delve into the research challenges and propose potential solutions for unraveling the genetic basis of rare and complex diseases.

Explicit ion modeling predicts physicochemical interactions for chromatin organization

Molecular mechanisms that dictate chromatin organization in vivo are under active investigation, and the extent to which intrinsic interactions contribute to this process remains debatable. A central quantity for evaluating their contribution is the strength of nucleosome-nucleosome binding, which previous experiments have estimated to range from 2 to 14 kBT. We introduce an explicit ion model to dramatically enhance the accuracy of residue-level coarse-grained modeling approaches across a wide range of ionic concentrations. This model allows for de novo predictions of chromatin organization and remains computationally efficient, enabling large-scale conformational sampling for free energy calculations. It reproduces the energetics of protein-DNA binding and unwinding of single nucleosomal DNA, and resolves the differential impact of mono- and divalent ions on chromatin conformations. Moreover, we showed that the model can reconcile various experiments on quantifying nucleosomal interactions, providing an explanation for the large discrepancy between existing estimations. We predict the interaction strength at physiological conditions to be 9 kBT, a value that is nonetheless sensitive to DNA linker length and the presence of linker histones. Our study strongly supports the contribution of physicochemical interactions to the phase behavior of chromatin aggregates and chromatin organization inside the nucleus.

Chromatin phase separation and nuclear shape fluctuations are correlated in a polymer model of the nucleus

Abnormalities in the shapes of mammalian cell nuclei are hallmarks of a variety of diseases, including progeria, muscular dystrophy, and various cancers. Experiments have shown that there is a causal relationship between chromatin organization and nuclear morphology. Decreases in heterochromatin levels, perturbations to heterochromatin organization, and increases in euchromatin levels all lead to misshapen nuclei, which exhibit deformations, such as nuclear blebs and nuclear ruptures. However, the polymer physical mechanisms of how chromatin governs nuclear shape and integrity are poorly understood. To investigate how heterochromatin and euchromatin, which are thought to microphase separate in vivo , govern nuclear morphology, we implemented a composite coarse-grained polymer and elastic shell model. By varying chromatin volume fraction (density), heterochromatin levels and structure, and heterochromatin-lamina interactions, we show how the spatial organization of chromatin polymer phases within the nucleus could perturb nuclear shape in some scenarios. Increasing the volume fraction of chromatin in the cell nucleus stabilizes the nuclear lamina against large fluctuations. However, surprisingly, we find that increasing heterochromatin levels or heterochromatin-lamina interactions enhances nuclear shape fluctuations in our simulations by a "wetting"-like interaction. In contrast, shape fluctuations are largely insensitive to the internal structure of the heterochromatin, such as the presence or absence of chromatin-chromatin crosslinks. Therefore, our simulations suggest that heterochromatin accumulation at the nuclear periphery could perturb nuclear morphology in a nucleus or nuclear region that is sufficiently soft, while stabilization of the nucleus via heterochromatin likely occurs through mechanisms other than chromatin microphase organization.

The Dynamics of Histone Modifications during Mammalian Zygotic Genome Activation

Mammalian fertilization initiates the reprogramming of oocytes and sperm, forming a totipotent zygote. During this intricate process, the zygotic genome undergoes a maternal-to-zygotic transition (MZT) and subsequent zygotic genome activation (ZGA), marking the initiation of transcriptional control and gene expression post-fertilization. Histone modifications are pivotal in shaping cellular identity and gene expression in many mammals. Recent advances in chromatin analysis have enabled detailed explorations of histone modifications during ZGA. This review delves into conserved and unique regulatory strategies, providing essential insights into the dynamic changes in histone modifications and their variants during ZGA in mammals. The objective is to explore recent advancements in leading mechanisms related to histone modifications governing this embryonic development phase in depth. These considerations will be useful for informing future therapeutic approaches that target epigenetic regulation in diverse biological contexts. It will also contribute to the extensive areas of evolutionary and developmental biology and possibly lay the foundation for future research and discussion on this seminal topic.

scENCORE: leveraging single-cell epigenetic data to predict chromatin conformation using graph embedding

Dynamic compartmentalization of eukaryotic DNA into active and repressed states enables diverse transcriptional programs to arise from a single genetic blueprint, whereas its dysregulation can be strongly linked to a broad spectrum of diseases. While single-cell Hi-C experiments allow for chromosome conformation profiling across many cells, they are still expensive and not widely available for most labs. Here, we propose an alternate approach, scENCORE, to computationally reconstruct chromatin compartments from the more affordable and widely accessible single-cell epigenetic data. First, scENCORE constructs a long-range epigenetic correlation graph to mimic chromatin interaction frequencies, where nodes and edges represent genome bins and their correlations. Then, it learns the node embeddings to cluster genome regions into A/B compartments and aligns different graphs to quantify chromatin conformation changes across conditions. Benchmarking using cell-type-matched Hi-C experiments demonstrates that scENCORE can robustly reconstruct A/B compartments in a cell-type-specific manner. Furthermore, our chromatin confirmation switching studies highlight substantial compartment-switching events that may introduce substantial regulatory and transcriptional changes in psychiatric disease. In summary, scENCORE allows accurate and cost-effective A/B compartment reconstruction to delineate higher-order chromatin structure heterogeneity in complex tissues.

Cancer Evolution: A Multifaceted Affair

Cancer cells adapt and survive through the acquisition and selection of molecular modifications. This process defines cancer evolution. Building on a theoretical framework based on heritable genetic changes has provided insights into the mechanisms supporting cancer evolution. However, cancer hallmarks also emerge via heritable nongenetic mechanisms, including epigenetic and chromatin topological changes, and interactions between tumor cells and the tumor microenvironment. Recent findings on tumor evolutionary mechanisms draw a multifaceted picture where heterogeneous forces interact and influence each other while shaping tumor progression. A comprehensive characterization of the cancer evolutionary toolkit is required to improve personalized medicine and biomarker discovery.

SIGNIFICANCE

Tumor evolution is fueled by multiple enabling mechanisms. Importantly, genetic instability, epigenetic reprogramming, and interactions with the tumor microenvironment are neither alternative nor independent evolutionary mechanisms. As demonstrated by findings highlighted in this perspective, experimental and theoretical approaches must account for multiple evolutionary mechanisms and their interactions to ultimately understand, predict, and steer tumor evolution.

Genome binding properties of Zic transcription factors underlie their changing functions during neuronal maturation

Background

The Zic family of transcription factors (TFs) promote both proliferation and maturation of cerebellar granule neurons (CGNs), raising the question of how a single, constitutively expressed TF family can support distinct developmental processes. Here we use an integrative experimental and bioinformatic approach to discover the regulatory relationship between Zic TF binding and changing programs of gene transcription during CGN differentiation.

Results

We first established a bioinformatic pipeline to integrate Zic ChIP-seq data from the developing mouse cerebellum with other genomic datasets from the same tissue. In newborn CGNs, Zic TF binding predominates at active enhancers that are co-bound by developmentally-regulated TFs including Atoh1, whereas in mature CGNs, Zic TF binding consolidates toward promoters where it co-localizes with activity-regulated TFs. We then performed CUT&RUN-seq in differentiating CGNs to define both the time course of developmental shifts in Zic TF binding and their relationship to gene expression. Mapping Zic TF binding sites to genes using chromatin looping, we identified the set of Zic target genes that have altered expression in RNA-seq from Zic1 or Zic2 knockdown CGNs.

Conclusion

Our data show that Zic TFs are required for both induction and repression of distinct, developmentally regulated target genes through a mechanism that is largely independent of changes in Zic TF binding. We suggest that the differential collaboration of Zic TFs with other TF families underlies the shift in their biological functions across CGN development.

Transcriptional reprogramming at the intersection of the heat shock response and proteostasis

Cellular homeostasis is constantly challenged by a myriad of extrinsic and intrinsic stressors. To mitigate the stress-induced damage, cells activate transient survival programs. The heat shock response (HSR) is an evolutionarily well-conserved survival program that is activated in response to proteotoxic stress. The HSR encompasses a dual regulation of transcription, characterized by rapid activation of genes encoding molecular chaperones and concomitant global attenuation of non-chaperone genes. Recent genome-wide approaches have delineated the molecular depth of stress-induced transcriptional reprogramming. The dramatic rewiring of gene and enhancer networks is driven by key transcription factors, including heat shock factors (HSFs), that together with chromatin-modifying enzymes remodel the 3D chromatin architecture, determining the selection of either gene activation or repression. Here, we highlight the current advancements of molecular mechanisms driving transcriptional reprogramming during acute heat stress. We also discuss the emerging implications of HSF-mediated stress signaling in the context of physiological and pathological conditions.

Deciphering the Role of Rapidly Evolving Conserved Elements in Primate Brain Development and Exploring Their Potential Involvement in Alzheimer's Disease

Although previous studies have identified human-specific accelerated regions as playing a key role in the recent evolution of the human brain, the characteristics and cellular functions of rapidly evolving conserved elements (RECEs) in ancestral primate lineages remain largely unexplored. Here, based on large-scale primate genome assemblies, we identify 888 RECEs that have been highly conserved in primates that exhibit significantly accelerated substitution rates in the ancestor of the Simiiformes. This primate lineage exhibits remarkable morphological innovations, including an expanded brain mass. Integrative multiomic analyses reveal that RECEs harbor sequences with potential cis-regulatory functions that are activated in the adult human brain. Importantly, genes linked to RECEs exhibit pronounced expression trajectories in the adult brain relative to the fetal stage. Furthermore, we observed an increase in the chromatin accessibility of RECEs in oligodendrocytes from individuals with Alzheimer's disease (AD) compared to that of a control group, indicating that these RECEs may contribute to brain aging and AD. Our findings serve to expand our knowledge of the genetic underpinnings of brain function during primate evolution.

Mapping nucleosome-resolution chromatin organization and enhancer-promoter loops in plants using Micro-C-XL

Although chromatin organizations in plants have been dissected at the scales of compartments and topologically associating domain (TAD)-like domains, there remains a gap in resolving fine-scale structures. Here, we use Micro-C-XL, a high-throughput chromosome conformation capture (Hi-C)-based technology that involves micrococcal nuclease (instead of restriction enzymes) and long cross-linkers, to dissect single nucleosome-resolution chromatin organization in Arabidopsis. Insulation analysis reveals more than 14,000 boundaries, which mostly include chromatin accessibility, epigenetic modifications, and transcription factors. Micro-C-XL reveals associations between RNA Pols and local chromatin organizations, suggesting that gene transcription substantially contributes to the establishment of local chromatin domains. By perturbing Pol II both genetically and chemically at the gene level, we confirm its function in regulating chromatin organization. Visible loops and stripes are assigned to super-enhancers and their targeted genes, thus providing direct insights for the identification and mechanistic analysis of distal CREs and their working modes in plants. We further investigate possible factors regulating these chromatin loops. Subsequently, we expand Micro-C-XL to soybean and rice. In summary, we use Micro-C-XL for analyses of plants, which reveal fine-scale chromatin organization and enhancer-promoter loops and provide insights regarding three-dimensional genomes in plants.

Modeling the compaction of bacterial chromosomes by biomolecular crowding and the cross-linking protein H-NS

Cells orchestrate the action of various molecules toward organizing their chromosomes. Using a coarse-grained computational model, we study the compaction of bacterial chromosomes by the cross-linking protein H-NS and cellular crowders. In this work, H-NS, modeled as a mobile "binder," can bind to a chromosome-like polymer with a characteristic binding energy. The simulation results reported here clarify the relative role of biomolecular crowding and H-NS in condensing a bacterial chromosome in a quantitative manner. In particular, they shed light on the nature and degree of crowder and H-NS synergetics: while the presence of crowders enhances H-NS binding to a chromosome-like polymer, the presence of H-NS makes crowding effects more efficient, suggesting two-way synergetics in chain compaction. Also, the results show how crowding effects promote clustering of bound H-NS. For a sufficiently large concentration of H-NS, the cluster size increases with the volume fraction of crowders.

The evolving views of hematopoiesis: from embryo to adulthood and from in vivo to in vitro

The hematopoietic system composed of hematopoietic stem and progenitor cells (HSPCs) and their differentiated lineages serves as an ideal model to uncover generic principles of cell fate transitions. From gastrulation onwards, there successively emerge primitive hematopoiesis (that produces specialized hematopoietic cells), pro-definitive hematopoiesis (that produces lineage-restricted progenitor cells), and definitive hematopoiesis (that produces multipotent HSPCs). These nascent lineages develop in several transient hematopoietic sites and finally colonize into lifelong hematopoietic sites. The development and maintenance of hematopoietic lineages are orchestrated by cell-intrinsic gene regulatory networks and cell-extrinsic microenvironmental cues. Owing to the progressive methodology (e.g., high-throughput lineage tracing and single-cell functional and omics analyses), our understanding of the developmental origin of hematopoietic lineages and functional properties of certain hematopoietic organs has been updated; meanwhile, new paradigms to characterize rare cell types, cell heterogeneity and its causes, and comprehensive regulatory landscapes have been provided. Here, we review the evolving views of HSPC biology during developmental and postnatal hematopoiesis. Moreover, we discuss recent advances in the in vitro induction and expansion of HSPCs, with a focus on the implications for clinical applications.