The 3D Genome as Moderator of Chromosomal Communication

Genome-wide analysis of the biophysical properties of chromatin and nuclear proteins in living cells with Hi-D

To understand the dynamic nature of the genome, the localization and rearrangement of DNA and DNA-binding proteins must be analyzed across the entire nucleus of single living cells. Recently, we developed a computational light microscopy technique, called high-resolution diffusion (Hi-D) mapping, which can accurately detect, classify and map diffusion dynamics and biophysical parameters such as the diffusion constant, the anomalous exponent, drift velocity and model physical diffusion from the data at a high spatial resolution across the genome in living cells. Hi-D combines dense optical flow to detect and track local chromatin and nuclear protein motion genome-wide and Bayesian inference to characterize this local movement at nanoscale resolution. Here we present the Python implementation of Hi-D, with an option for parallelizing the calculations to run on multicore central processing units (CPUs). The functionality of Hi-D is presented to the users via user-friendly documented Python notebooks. Hi-D reduces the analysis time to less than 1 h using a multicore CPU with a single compute node. We also present different applications of Hi-D for live-imaging of DNA, histone H2B and RNA polymerase II sequences acquired with spinning disk confocal and super-resolution structured illumination microscopy.

OpenNucleome for high-resolution nuclear structural and dynamical modeling

The intricate structural organization of the human nucleus is fundamental to cellular function and gene regulation. Recent advancements in experimental techniques, including high-throughput sequencing and microscopy, have provided valuable insights into nuclear organization. Computational modeling has played significant roles in interpreting experimental observations by reconstructing high-resolution structural ensembles and uncovering organization principles. However, the absence of standardized modeling tools poses challenges for furthering nuclear investigations. We present OpenNucleome-an open-source software designed for conducting GPU-accelerated molecular dynamics simulations of the human nucleus. OpenNucleome offers particle-based representations of chromosomes at a resolution of 100 KB, encompassing nuclear lamina, nucleoli, and speckles. This software furnishes highly accurate structural models of nuclear architecture, affording the means for dynamic simulations of condensate formation, fusion, and exploration of non-equilibrium effects. We applied OpenNucleome to uncover the mechanisms driving the emergence of 'fixed points' within the nucleus-signifying genomic loci robustly anchored in proximity to specific nuclear bodies for functional purposes. This anchoring remains resilient even amidst significant fluctuations in chromosome radial positions and nuclear shapes within individual cells. Our findings lend support to a nuclear zoning model that elucidates genome functionality. We anticipate OpenNucleome to serve as a valuable tool for nuclear investigations, streamlining mechanistic explorations and enhancing the interpretation of experimental observations.

Evolutional heterochromatin condensation delineates chromocenter formation and retrotransposon silencing in plants

Heterochromatic condensates (chromocenters) are critical for maintaining the silencing of heterochromatin. It is therefore puzzling that the presence of chromocenters is variable across plant species. Here we reveal that variations in the plant heterochromatin protein ADCP1 confer a diversity in chromocenter formation via phase separation. ADCP1 physically interacts with the high mobility group protein HMGA to form a complex and mediates heterochromatin condensation by multivalent interactions. The loss of intrinsically disordered regions (IDRs) in ADCP1 homologues during evolution has led to the absence of prominent chromocenter formation in various plant species, and introduction of IDR-containing ADCP1 with HMGA promotes heterochromatin condensation and retrotransposon silencing. Moreover, plants in the Cucurbitaceae group have evolved an IDR-containing chimaera of ADCP1 and HMGA, which remarkably enables formation of chromocenters. Together, our work uncovers a coevolved mechanism of phase separation in packing heterochromatin and silencing retrotransposons.

Enhancers on the edge - how the nuclear envelope controls gene regulatory elements

Precise temporal and sequential control of gene expression during development and in response to environmental stimuli requires tight regulation of the physical contact between gene regulatory elements and promoters. Current models describing how the genome folds in 3D space to establish these interactions often ignore the role of the most stable structural nuclear feature - the nuclear envelope. While contributions of 3D folding within/between topologically associated domains (TADs) have been extensively described, mechanical contributions from the nuclear envelope can impact enhancer-promoter interactions both directly and indirectly through influencing intra/inter-TAD interactions. Importantly, these nuclear envelope contributions clearly link this mechanism to development and, when defective, to human disease. Here, we discuss evidence for nuclear envelope regulation of tissue-specific enhancer-promoter pairings, potential mechanisms for this regulation, exciting recent findings that other regulatory elements such as microRNAs and long noncoding RNAs are under nuclear envelope regulation, the possible involvement of condensates, and how disruption of this regulation can lead to disease.

SATB1, senescence and senescence-related diseases

Aging leads to an accumulation of cellular mutations and damage, increasing the risk of senescence, apoptosis, and malignant transformation. Cellular senescence, which is pivotal in aging, acts as both a guard against cellular transformation and as a check against cancer progression. It is marked by stable cell cycle arrest, widespread macromolecular changes, a pro-inflammatory profile, and altered gene expression. However, it remains to be determined whether these differing subsets of senescent cells result from unique intrinsic programs or are influenced by their environmental contexts. Multiple transcription regulators and chromatin modifiers contribute to these alterations. Special AT-rich sequence-binding protein 1 (SATB1) stands out as a crucial regulator in this process, orchestrating gene expression by structuring chromatin into loop domains and anchoring DNA elements. This review provides an overview of cellular senescence and delves into the role of SATB1 in senescence-related diseases. It highlights SATB1's potential in developing antiaging and anticancer strategies, potentially contributing to improved quality of life and addressing aging-related diseases.

Unsupervised representation learning of chromatin images identifies changes in cell state and tissue organization in DCIS

Ductal carcinoma in situ (DCIS) is a pre-invasive tumor that can progress to invasive breast cancer, a leading cause of cancer death. We generate a large-scale tissue microarray dataset of chromatin images, from 560 samples from 122 female patients in 3 disease stages and 11 phenotypic categories. Using representation learning on chromatin images alone, without multiplexed staining or high-throughput sequencing, we identify eight morphological cell states and tissue features marking DCIS. All cell states are observed in all disease stages with different proportions, indicating that cell states enriched in invasive cancer exist in small fractions in normal breast tissue. Tissue-level analysis reveals significant changes in the spatial organization of cell states across disease stages, which is predictive of disease stage and phenotypic category. Taken together, we show that chromatin imaging represents a powerful measure of cell state and disease stage of DCIS, providing a simple and effective tumor biomarker.

Polymer physics models reveal structural folding features of single-molecule gene chromatin conformations

Here, we employ polymer physics models of chromatin to investigate the 3D folding of a 2Mb wide genomic region encompassing the human LTN1 gene, a crucial DNA locus involved in key cellular functions. Through extensive Molecular Dynamics simulations, we reconstruct in-silico the ensemble of single-molecule LTN1 3D structures, which we benchmark against recent in-situ Hi-C 2.0 data. The model-derived single molecules are then used to predict structural folding features at the single-cell level, providing testable predictions for super-resolution microscopy experiments.

Enhancer selectivity in space and time: from enhancer-promoter interactions to promoter activation

The primary regulators of metazoan gene expression are enhancers, originally functionally defined as DNA sequences that can activate transcription at promoters in an orientation-independent and distance-independent manner. Despite being crucial for gene regulation in animals, what mechanisms underlie enhancer selectivity for promoters, and more fundamentally, how enhancers interact with promoters and activate transcription, remain poorly understood. In this Review, we first discuss current models of enhancer-promoter interactions in space and time and how enhancers affect transcription activation. Next, we discuss different mechanisms that mediate enhancer selectivity, including repression, biochemical compatibility and regulation of 3D genome structure. Through 3D polymer simulations, we illustrate how the ability of 3D genome folding mechanisms to mediate enhancer selectivity strongly varies for different enhancer-promoter interaction mechanisms. Finally, we discuss how recent technical advances may provide new insights into mechanisms of enhancer-promoter interactions and how technical biases in methods such as Hi-C and Micro-C and imaging techniques may affect their interpretation.

From Nucleosomes to Compartments: Physicochemical Interactions Underlying Chromatin Organization

Chromatin organization plays a critical role in cellular function by regulating access to genetic information. However, understanding chromatin folding is challenging due to its complex, multiscale nature. Significant progress has been made in studying in vitro systems, uncovering the structure of individual nucleosomes and their arrays, and elucidating the role of physicochemical forces in stabilizing these structures. Additionally, remarkable advancements have been achieved in characterizing chromatin organization in vivo, particularly at the whole-chromosome level, revealing important features such as chromatin loops, topologically associating domains, and nuclear compartments. However, bridging the gap between in vitro and in vivo studies remains challenging. The resemblance between in vitro and in vivo chromatin conformations and the relevance of internucleosomal interactions for chromatin folding in vivo are subjects of debate. This article reviews experimental and computational studies conducted at various length scales, highlighting the significance of intrinsic interactions between nucleosomes and their roles in chromatin folding in vivo.

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.

Differential 3D genome architecture and imprinted gene expression: cause or consequence?

Imprinted genes provide an attractive paradigm to unravel links between transcription and genome architecture. The parental allele-specific expression of these essential genes - which are clustered in chromosomal domains - is mediated by parental methylation imprints at key regulatory DNA sequences. Recent chromatin conformation capture (3C)-based studies show differential organization of topologically associating domains between the parental chromosomes at imprinted domains, in embryonic stem and differentiated cells. At several imprinted domains, differentially methylated regions show allelic binding of the insulator protein CTCF, and linked focal retention of cohesin, at the non-methylated allele only. This generates differential patterns of chromatin looping between the parental chromosomes, already in the early embryo, and thereby facilitates the allelic gene expression. Recent research evokes also the opposite scenario, in which allelic transcription contributes to the differential genome organization, similarly as reported for imprinted X chromosome inactivation. This may occur through epigenetic effects on CTCF binding, through structural effects of RNA Polymerase II, or through imprinted long non-coding RNAs that have chromatin repressive functions. The emerging picture is that epigenetically-controlled differential genome architecture precedes and facilitates imprinted gene expression during development, and that at some domains, conversely, the mono-allelic gene expression also influences genome architecture.

Dynamic RNA polymerase II occupancy drives differentiation of the intestine under the direction of HNF4

Terminal differentiation requires massive restructuring of the transcriptome. During intestinal differentiation, the expression patterns of nearly 4,000 genes are altered as cells transition from progenitor cells in crypts to differentiated cells in villi. We identify dynamic occupancy of RNA polymerase II (Pol II) to gene promoters as the primary driver of transcriptomic shifts during intestinal differentiation in vivo. Changes in enhancer-promoter looping interactions accompany dynamic Pol II occupancy and are dependent upon HNF4, a pro-differentiation transcription factor. Using genetic loss-of-function, chromatin immunoprecipitation sequencing (ChIP-seq), and immunoprecipitation (IP) mass spectrometry, we demonstrate that HNF4 collaborates with chromatin remodelers and loop-stabilizing proteins and facilitates Pol II occupancy at hundreds of genes pivotal to differentiation. We also explore alternate mechanisms that drive differentiation gene expression and find that pause-release of Pol II and post-transcriptional mRNA stability regulate smaller subsets of differentially expressed genes. These studies provide insights into the mechanisms of differentiation in renewing adult tissue.

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.

A high-throughput synthetic biology approach for studying combinatorial chromatin-based transcriptional regulation

The construction of synthetic gene circuits requires the rational combination of multiple regulatory components, but predicting their behavior can be challenging due to poorly understood component interactions and unexpected emergent behaviors. In eukaryotes, chromatin regulators (CRs) are essential regulatory components that orchestrate gene expression. Here, we develop a screening platform to investigate the impact of CR pairs on transcriptional activity in yeast. We construct a combinatorial library consisting of over 1,900 CR pairs and use a high-throughput workflow to characterize the impact of CR co-recruitment on gene expression. We recapitulate known interactions and discover several instances of CR pairs with emergent behaviors. We also demonstrate that supervised machine learning models trained with low-dimensional amino acid embeddings accurately predict the impact of CR co-recruitment on transcriptional activity. This work introduces a scalable platform and machine learning approach that can be used to study how networks of regulatory components impact gene expression.

DNA packaging by molecular motors: from bacteriophage to human chromosomes

Dense packaging of genomic DNA is crucial for organismal survival, as DNA length always far exceeds the dimensions of the cells that contain it. Organisms, therefore, use sophisticated machineries to package their genomes. These systems range across kingdoms from a single ultra-powerful rotary motor that spools the DNA into a bacteriophage head, to hundreds of thousands of relatively weak molecular motors that coordinate the compaction of mitotic chromosomes in eukaryotic cells. Recent technological advances, such as DNA proximity-based sequencing approaches, polymer modelling and in vitro reconstitution of DNA loop extrusion, have shed light on the biological mechanisms driving DNA organization in different systems. Here, we discuss DNA packaging in bacteriophage, bacteria and eukaryotic cells, which, despite their extreme variation in size, structure and genomic content, all rely on the action of molecular motors to package their genomes.

Microglia homeostasis mediated by epigenetic ARID1A regulates neural progenitor cells response and leads to autism-like behaviors

Microglia are resident macrophages of the central nervous system that selectively emerge in embryonic cortical proliferative zones and regulate neurogenesis by altering molecular and phenotypic states. Despite their important roles in inflammatory phagocytosis and neurodegenerative diseases, microglial homeostasis during early brain development has not been fully elucidated. Here, we demonstrate a notable interplay between microglial homeostasis and neural progenitor cell signal transduction during embryonic neurogenesis. ARID1A, an epigenetic subunit of the SWI/SNF chromatin-remodeling complex, disrupts genome-wide H3K9me3 occupancy in microglia and changes the epigenetic chromatin landscape of regulatory elements that influence the switching of microglial states. Perturbation of microglial homeostasis impairs the release of PRG3, which regulates neural progenitor cell self-renewal and differentiation during embryonic development. Furthermore, the loss of microglia-driven PRG3 alters the downstream cascade of the Wnt/β-catenin signaling pathway through its interaction with the neural progenitor receptor LRP6, which leads to misplaced regulation in neuronal development and causes autism-like behaviors at later stages. Thus, during early fetal brain development, microglia progress toward a more homeostatic competent phenotype, which might render neural progenitor cells respond to environmental cross-talk perturbations.

Single-Macromolecule Studies of Eukaryotic Genomic Maintenance

Genomes are self-organized and self-maintained as long, complex macromolecules of chromatin. The inherent heterogeneity, stochasticity, phase separation, and chromatin dynamics of genome operation make it challenging to study genomes using ensemble methods. Various single-molecule force-, fluorescent-, and sequencing-based techniques rooted in different disciplines have been developed to fill critical gaps in the capabilities of bulk measurements, each providing unique, otherwise inaccessible, insights into the structure and maintenance of the genome. Capable of capturing molecular-level details about the organization, conformational changes, and packaging of genetic material, as well as processive and stochastic movements of maintenance factors, a single-molecule toolbox provides an excellent opportunity for collaborative research to understand how genetic material functions in health and malfunctions in disease. In this review, we discuss novel insights brought to genomic sciences by single-molecule techniques and their potential to continue to revolutionize the field-one molecule at a time.

Chromatin compaction by Polycomb group proteins revisited

The chromatin compaction activity of Polycomb group proteins has traditionally been considered essential for transcriptional repression. However, there is very little information on how Polycomb group proteins compact chromatin at the molecular level and no causal link between the compactness of chromatin and transcriptional repression. Recently, a more complete picture of Polycomb-dependent chromatin architecture has started to emerge, owing to advanced methods for imaging and chromosome conformation capture. Discoveries into Polycomb-driven phase separation add another layer of complexity. Recent observations generally imply that Polycomb group proteins modulate chromatin structure at multiple scales to reduce its dynamics and segregate it from active domains. Hence, it is reasonable to hypothesise that Polycomb group proteins maintain the energetically favourable state of compacted chromatin, rather than actively compact it.

Phase separation and transcriptional regulation in cancer development

Liquid-liquid phase separation, a novel biochemical phenomenon, has been increasingly studied for its medical applications. It underlies the formation of membrane-less organelles and is involved in many cellular and biological processes. During transcriptional regulation, dynamic condensates are formed through interactions between transcriptional elements, such as transcription factors, coactivators, and mediators. Cancer is a disease characterized by uncontrolled cell proliferation, but the precise mechanisms underlying tumorigenesis often remain to be elucidated. Emerging evidence has linked abnormal transcriptional condensates to several diseases, especially cancer, implying that phase separation plays an important role in tumorigenesis. Condensates formed by phase separation may have an effect on gene transcription in tumors. In the present review, we focus on the correlation between phase separation and transcriptional regulation, as well as how this phenomenon contributes to cancer development.

Identifying Protein Interactions with Viral DNA Genomes during Virus Infection

Viruses exploit the host cell machinery to enable infection and propagation. This review discusses the complex landscape of DNA virus-host interactions, focusing primarily on herpesviruses and adenoviruses, which replicate in the nucleus of infected cells, and vaccinia virus, which replicates in the cytoplasm. We discuss experimental approaches used to discover and validate interactions of host proteins with viral genomes and how these interactions impact processes that occur during infection, including the host DNA damage response and viral genome replication, repair, and transcription. We highlight the current state of knowledge regarding virus-host protein interactions and also outline emerging areas and future directions for research.

Potential roles of inter-chromosomal interactions in cell fate determination

Mammalian genomic DNA is packed in a small nucleus, and its folding and organization in the nucleus are critical for gene regulation and cell fate determination. In interphase, chromosomes are compartmentalized into certain nuclear spaces and territories that are considered incompatible with each other. The regulation of gene expression is influenced by the epigenetic characteristics of topologically associated domains and A/B compartments within chromosomes (intrachromosomal). Previously, interactions among chromosomes detected via chromosome conformation capture-based methods were considered noise or artificial errors. However, recent studies based on newly developed ligation-independent methods have shown that inter-chromosomal interactions play important roles in gene regulation. This review summarizes the recent understanding of spatial genomic organization in mammalian interphase nuclei and discusses the potential mechanisms that determine cell identity. In addition, this review highlights the potential role of inter-chromosomal interactions in early mouse development.

Individual transcription factors modulate both the micromovement of chromatin and its long-range structure

The control of eukaryotic gene expression is intimately connected to highly dynamic chromatin structures. Gene regulation relies on activator and repressor transcription factors (TFs) that induce local chromatin opening and closing. However, it is unclear how nucleus-wide chromatin organization responds dynamically to the activity of specific TFs. Here, we examined how two TFs with opposite effects on local chromatin accessibility modulate chromatin dynamics nucleus-wide. We combine high-resolution diffusion mapping and dense flow reconstruction and correlation in living cells to obtain an imaging-based, nanometer-scale analysis of local diffusion processes and long-range coordinated movements of both chromatin and TFs. We show that the expression of either an individual transcriptional activator (CDX2) or repressor (SIX6) with large numbers of binding sites increases chromatin mobility nucleus-wide, yet they induce opposite coherent chromatin motions at the micron scale. Hi-C analysis of higher-order chromatin structures shows that induction of the pioneer factor CDX2 leads both to changes in local chromatin interactions and the distribution of A and B compartments, thus relating the micromovement of chromatin with changes in compartmental structures. Given that inhibition of transcription initiation and elongation by RNA Pol II has a partial impact on the global chromatin dynamics induced by CDX2, we suggest that CDX2 overexpression alters chromatin structure dynamics both dependently and independently of transcription. Our biophysical analysis shows that sequence-specific TFs can influence chromatin structure on multiple architectural levels, arguing that local chromatin changes brought by TFs alter long-range chromatin mobility and its organization.

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.

Histone H1.0 couples cellular mechanical behaviors to chromatin structure

Tuning of genome structure and function is accomplished by chromatin-binding proteins, which determine the transcriptome and phenotype of the cell. Here we investigate how communication between extracellular stress and chromatin structure may regulate cellular mechanical behaviors. We demonstrate that histone H1.0, which compacts nucleosomes into higher-order chromatin fibers, controls genome organization and cellular stress response. We show that histone H1.0 has privileged expression in fibroblasts across tissue types and that its expression is necessary and sufficient to induce myofibroblast activation. Depletion of histone H1.0 prevents cytokine-induced fibroblast contraction, proliferation and migration via inhibition of a transcriptome comprising extracellular matrix, cytoskeletal and contractile genes, through a process that involves locus-specific H3K27 acetylation. Transient depletion of histone H1.0 in vivo prevents fibrosis in cardiac muscle. These findings identify an unexpected role of linker histones to orchestrate cellular mechanical behaviors, directly coupling force generation, nuclear organization and gene transcription.

Effect of loops on the mean-square displacement of Rouse-model chromatin

Chromatin polymer dynamics are commonly described using the classical Rouse model. The subsequent discovery, however, of intermediate-scale chromatin organization known as topologically associating domains (TADs) in experimental Hi-C contact maps for chromosomes across the tree of life, together with the success of loop extrusion factor (LEF) model in explaining TAD formation, motivates efforts to understand the effect of loops and loop extrusion on chromatin dynamics. This paper seeks to fulfill this need by combining LEF-model simulations with extended Rouse-model polymer simulations to investigate the dynamics of chromatin with loops and dynamic loop extrusion. We show that loops significantly suppress the averaged mean-square displacement (MSD) of a gene locus, consistent with recent experiments that track fluorescently labeled chromatin loci. We also find that loops reduce the MSD's stretching exponent from the classical Rouse-model value of 1/2 to a loop-density-dependent value in the 0.45-0.40 range. Remarkably, stretching exponent values in this range have also been observed in recent experiments [Weber et al., Phys. Rev. Lett. 104, 238102 (2010)0031-900710.1103/PhysRevLett.104.238102; Bailey et al., Mol. Biol. Cell 34, ar78 (2023)1059-152410.1091/mbc.E23-04-0119]. We also show that the dynamics of loop extrusion itself negligibly affects chromatin mobility. By studying static "rosette" loop configurations, we also demonstrate that chromatin MSDs and stretching exponents depend on the location of the locus in question relative to the position of the loops and on the local friction environment.

OpenNucleome for high resolution nuclear structural and dynamical modeling

The intricate structural organization of the human nucleus is fundamental to cellular function and gene regulation. Recent advancements in experimental techniques, including high-throughput sequencing and microscopy, have provided valuable insights into nuclear organization. Computational modeling has played significant roles in interpreting experimental observations by reconstructing high-resolution structural ensembles and uncovering organization principles. However, the absence of standardized modeling tools poses challenges for furthering nuclear investigations. We present OpenNucleome-an open-source software designed for conducting GPU-accelerated molecular dynamics simulations of the human nucleus. OpenNucleome offers particle-based representations of chromosomes at a resolution of 100 KB, encompassing nuclear lamina, nucleoli, and speckles. This software furnishes highly accurate structural models of nuclear architecture, affording the means for dynamic simulations of condensate formation, fusion, and exploration of non-equilibrium effects. We applied OpenNucleome to uncover the mechanisms driving the emergence of "fixed points" within the nucleus-signifying genomic loci robustly anchored in proximity to specific nuclear bodies for functional purposes. This anchoring remains resilient even amidst significant fluctuations in chromosome radial positions and nuclear shapes within individual cells. Our findings lend support to a nuclear zoning model that elucidates genome functionality. We anticipate OpenNucleome to serve as a valuable tool for nuclear investigations, streamlining mechanistic explorations and enhancing the interpretation of experimental observations.

Sister chromatid cohesion halts DNA loop expansion

Eukaryotic genomes are folded into DNA loops mediated by structural maintenance of chromosomes (SMC) complexes such as cohesin, condensin, and Smc5/6. This organization regulates different DNA-related processes along the cell cycle, such as transcription, recombination, segregation, and DNA repair. During the G2 stage, SMC-mediated DNA loops coexist with cohesin complexes involved in sister chromatid cohesion (SCC). However, the articulation between the establishment of SCC and the formation of SMC-mediated DNA loops along the chromatin remains unknown. Here, we show that SCC is indeed a barrier to cohesin-mediated DNA loop expansion along G2/M Saccharomyces cerevisiae chromosomes.

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.

Machine and deep learning methods for predicting 3D genome organization

Three-Dimensional (3D) chromatin interactions, such as enhancer-promoter interactions (EPIs), loops, Topologically Associating Domains (TADs), and A/B compartments play critical roles in a wide range of cellular processes by regulating gene expression. Recent development of chromatin conformation capture technologies has enabled genome-wide profiling of various 3D structures, even with single cells. However, current catalogs of 3D structures remain incomplete and unreliable due to differences in technology, tools, and low data resolution. Machine learning methods have emerged as an alternative to obtain missing 3D interactions and/or improve resolution. Such methods frequently use genome annotation data (ChIP-seq, DNAse-seq, etc.), DNA sequencing information (k-mers, Transcription Factor Binding Site (TFBS) motifs), and other genomic properties to learn the associations between genomic features and chromatin interactions. In this review, we discuss computational tools for predicting three types of 3D interactions (EPIs, chromatin interactions, TAD boundaries) and analyze their pros and cons. We also point out obstacles of computational prediction of 3D interactions and suggest future research directions.

Transcriptional changes are tightly coupled to chromatin reorganization during cellular aging

Human life expectancy is constantly increasing and aging has become a major risk factor for many diseases, although the underlying gene regulatory mechanisms are still unclear. Using transcriptomic and chromosomal conformation capture (Hi-C) data from human skin fibroblasts from individuals across different age groups, we identified a tight coupling between the changes in co-regulation and co-localization of genes. We obtained transcription factors, cofactors, and chromatin regulators that could drive the cellular aging process by developing a time-course prize-collecting Steiner tree algorithm. In particular, by combining RNA-Seq data from different age groups and protein-protein interaction data we determined the key transcription regulators and gene regulatory changes at different life stage transitions. We then mapped these transcription regulators to the 3D reorganization of chromatin in young and old skin fibroblasts. Collectively, we identified key transcription regulators whose target genes are spatially rearranged and correlate with changes in their expression, thereby providing potential targets for reverting cellular aging.

Enhancer-promoter interactions can form independently of genomic distance and be functional across TAD boundaries

Topologically Associating Domains (TADs) have been suggested to facilitate and constrain enhancer-promoter interactions. However, the role of TAD boundaries in effectively restricting these interactions remains unclear. Here, we show that a significant proportion of enhancer-promoter interactions are established across TAD boundaries in Drosophila embryos, but that developmental genes are strikingly enriched in intra- but not inter-TAD interactions. We pursued this observation using the twist locus, a master regulator of mesoderm development, and systematically relocated one of its enhancers to various genomic locations. While this developmental gene can establish inter-TAD interactions with its enhancer, the functionality of these interactions remains limited, highlighting the existence of topological constraints. Furthermore, contrary to intra-TAD interactions, the formation of inter-TAD enhancer-promoter interactions is not solely driven by genomic distance, with distal interactions sometimes favored over proximal ones. These observations suggest that other general mechanisms must exist to establish and maintain specific enhancer-promoter interactions across large distances.

KIT mutations and expression: current knowledge and new insights for overcoming IM resistance in GIST

Gastrointestinal stromal tumor (GIST) is the most common sarcoma located in gastrointestinal tract and derived from the interstitial cell of Cajal (ICC) lineage. Both ICC and GIST cells highly rely on KIT signal pathway. Clinically, about 80-90% of treatment-naive GIST patients harbor primary KIT mutations, and special KIT-targeted TKI, imatinib (IM) showing dramatic efficacy but resistance invariably occur, 90% of them was due to the second resistance mutations emerging within the KIT gene. Although there are multiple variants of KIT mutant which did not show complete uniform biologic characteristics, most of them have high KIT expression level. Notably, the high expression level of KIT gene is not correlated to its gene amplification. Recently, accumulating evidences strongly indicated that the gene coding, epigenetic regulation, and pre- or post- protein translation of KIT mutants in GIST were quite different from that of wild type (WT) KIT. In this review, we elucidate the biologic mechanism of KIT variants and update the underlying mechanism of the expression of KIT gene, which are exclusively regulated in GIST, providing a promising yet evidence-based therapeutic landscape and possible target for the conquer of IM resistance. Video Abstract.

Power-law behavior of transcriptional bursting regulated by enhancer-promoter communication

Revealing how transcriptional bursting kinetics are genomically encoded is challenging because genome structures are stochastic at the organization level and are suggestively linked to gene transcription. To address this challenge, we develop a generic theoretical framework that integrates chromatin dynamics, enhancer-promoter (E-P) communication, and gene-state switching to study transcriptional bursting. The theory predicts that power law can be a general rule to quantitatively describe bursting modulations by E-P spatial communication. Specifically, burst frequency and burst size are up-regulated by E-P communication strength, following power laws with positive exponents. Analysis of the scaling exponents further reveals that burst frequency is preferentially regulated. Bursting kinetics are down-regulated by E-P genomic distance with negative power-law exponents, and this negative modulation desensitizes at large distances. The mutual information between burst frequency (or burst size) and E-P spatial distance further reveals essential characteristics of the information transfer from E-P communication to transcriptional bursting kinetics. These findings, which are in agreement with experimental observations, not only reveal fundamental principles of E-P communication in transcriptional bursting but also are essential for understanding cellular decision-making.

Boundary stacking interactions enable cross-TAD enhancer-promoter communication during limb development

Although promoters and their enhancers are frequently contained within a topologically associating domain (TAD), some developmentally important genes have their promoter and enhancers within different TADs. Hypotheses about molecular mechanisms enabling cross-TAD interactions remain to be assessed. To test these hypotheses, we used optical reconstruction of chromatin architecture to characterize the conformations of the Pitx1 locus on single chromosomes in developing mouse limbs. Our data support a model in which neighboring boundaries are stacked as a result of loop extrusion, bringing boundary-proximal cis-elements into contact. This stacking interaction also contributes to the appearance of architectural stripes in the population average maps. Through molecular dynamics simulations, we found that increasing boundary strengths facilitates the formation of the stacked boundary conformation, counter-intuitively facilitating border bypass. This work provides a revised view of the TAD borders' function, both facilitating and preventing cis-regulatory interactions, and introduces a framework to distinguish border-crossing from border-respecting enhancer-promoter pairs.

Chromatin structure from high resolution microscopy: Scaling laws and microphase separation

Recent advances in experimental fluorescence microscopy allow high accuracy determination (resolution of 50 nm) of the three-dimensional physical location of multiple (up to ∼10^{2}) tagged regions of the chromosome. We investigate publicly available microscopy data for two loci of the human Chr21 obtained from multiplexed fluorescence in situ hybridization (FISH) methods for different cell lines and treatments. Inspired by polymer physics models, our analysis centers around distance distributions between different tags with the aim being to unravel the chromatin conformational arrangements. We show that for any specific genomic site, there are (at least) two different conformational arrangements of chromatin, implying coexisting distinct topologies which we refer to as phase α and phase β. These two phases show different scaling behaviors: the former is consistent with a crumpled globule, while the latter indicates a confined, but more extended conformation, such as a looped domain. The identification of these distinct phases sheds light on the coexistence of multiple chromatin topologies and provides insights into the effects of cellular context and/or treatments on chromatin structure.

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.

3D genome structural variations play important roles in regulating seed oil content of Brassica napus

Dissecting the complex regulatory mechanism of seed oil content (SOC) is one of the main research goals in Brassica napus. Increasing evidence suggests that genome architecture is linked to multiple biological functions. However, the effect of genome architecture on SOC regulation remains unclear. Here, we used high-throughput chromatin conformation capture to characterize differences in the three-dimensional (3D) landscape of genome architecture of seeds from two B. napus lines, N53-2 (with high SOC) and Ken-C8 (with low SOC). Bioinformatics analysis demonstrated that differentially accessible regions and differentially expressed genes between N53-2 and Ken-C8 were preferentially enriched in regions with quantitative trait loci (QTLs)/associated genomic regions (AGRs) for SOC. A multi-omics analysis demonstrated that expression of SOC-related genes was tightly correlated with genome structural variations in QTLs/AGRs of B. napus. The candidate gene BnaA09g48250D, which showed structural variation in a QTL/AGR on chrA09, was identified by fine-mapping of a KN double-haploid population derived from hybridization of N53-2 and Ken-C8. Overexpression and knockout of BnaA09g48250D led to significant increases and decreases in SOC, respectively, in the transgenic lines. Taken together, our results reveal the 3D genome architecture of B. napus seeds and the roles of genome structural variations in SOC regulation, enriching our understanding of the molecular mechanisms of SOC regulation from the perspective of spatial chromatin structure.

Replicating Chromosomes in Whole-Cell Models of Bacteria

Computational models of cells cannot be considered complete unless they include the most fundamental process of life, the replication of genetic material. In a recent study, we presented a computational framework to model systems of replicating bacterial chromosomes as polymers at 10 bp resolution with Brownian dynamics. This approach was used to investigate changes in chromosome organization during replication and extend the applicability of an existing whole-cell model (WCM) for a genetically minimal bacterium, JCVI-syn3A, to the entire cell cycle. To achieve cell-scale chromosome structures that are realistic, we modeled the chromosome as a self-avoiding homopolymer with bending and torsional stiffnesses that capture the essential mechanical properties of dsDNA in Syn3A. Additionally, the polymer interacts with ribosomes distributed according to cryo-electron tomograms of Syn3A. The polymer model was further augmented by computational models of loop extrusion by structural maintenance of chromosomes (SMC) protein complexes and topoisomerase action, and the modeling and analysis of multi-fork replication states.

A chromatin-regulated biphasic circuit coordinates IL-1β-mediated inflammation

Inflammation is characterized by a biphasic cycle consisting initially of a proinflammatory phase that is subsequently resolved by anti-inflammatory processes. Interleukin-1β (IL-1β) is a master regulator of proinflammation and is encoded within the same topologically associating domain (TAD) as IL-37, which is an anti-inflammatory cytokine that opposes the function of IL-1β. Within this TAD, we identified a long noncoding RNA called AMANZI, which negatively regulates IL-1β expression and trained immunity through the induction of IL37 transcription. We found that the activation of IL37 occurs through the formation of a dynamic long-range chromatin contact that leads to the temporal delay of anti-inflammatory responses. The common variant rs16944 present in AMANZI augments this regulatory circuit, predisposing individuals to enhanced proinflammation or immunosuppression. Our work illuminates a chromatin-mediated biphasic circuit coordinating expression of IL-1β and IL-37, thereby regulating two functionally opposed states of inflammation from within a single TAD.

BRWD1 orchestrates small pre-B cell chromatin topology by converting static to dynamic cohesin

Lymphocyte development consists of sequential and mutually exclusive cell states of proliferative selection and antigen receptor gene recombination. Transitions between each state require large, coordinated changes in epigenetic landscapes and transcriptional programs. How this occurs remains unclear. Here we demonstrate that in small pre-B cells, the lineage and stage-specific epigenetic reader bromodomain and WD repeat-containing protein 1 (BRWD1) reorders three-dimensional chromatin topology to affect the transition between proliferative and gene recombination molecular programs. BRWD1 regulated the switch between poised and active enhancers interacting with promoters, and coordinated this switch with Igk locus contraction. BRWD1 did so by converting chromatin-bound static to dynamic cohesin competent to mediate long-range looping. ATP-depletion revealed cohesin conversion to be the main energetic mechanism dictating dynamic chromatin looping. Our findings provide a new mechanism of cohesin regulation and reveal how cohesin function can be dictated by lineage contextual mechanisms to facilitate specific cell fate transitions.

Single Nucleotide Polymorphisms (SNPs) in the Shadows: Uncovering their Function in Non-Coding Region of Esophageal Cancer

Esophageal cancer is a complex disease influenced by genetic and environmental factors. Single nucleotide polymorphisms (SNPs) in non-coding regions of the genome have emerged as crucial contributors to esophageal cancer susceptibility. This review provides a comprehensive overview of the role of SNPs in non-coding regions and their association with esophageal cancer. The accumulation of SNPs in the genome has been implicated in esophageal cancer risk. Various studies have identified specific locations in the genome where SNPs are more likely to occur, suggesting a location-specific response. Chromatin conformational studies have shed light on the localization of SNPs and their impact on gene transcription, posttranscriptional modifications, gene expression regulation, and histone modification. Furthermore, miRNA-related SNPs have been found to play a significant role in esophageal squamous cell carcinoma (ESCC). These SNPs can affect miRNA binding sites, thereby altering target gene regulation and contributing to ESCC development. Additionally, the risk of ESCC has been linked to base excision repair, suggesting that SNPs in this pathway may influence disease susceptibility. Somatic DNA segment alterations and modified expression quantitative trait loci (eQTL) have also been associated with ESCC. These alterations can lead to disrupted gene expression and cellular processes, ultimately contributing to cancer development and progression. Moreover, SNPs have been found to be associated with the long non-coding RNA HOTAIR, which plays a crucial role in ESCC pathogenesis. This review concludes with a discussion of the current and future perspectives in the field of SNPs in non-coding regions and their relevance to esophageal cancer. Understanding the functional implications of these SNPs may lead to the identification of novel therapeutic targets and the development of personalized approaches for esophageal cancer prevention and treatment.

Direct Reprogramming of Mouse Fibroblasts to Osteoblast-like Cells Using Runx2/Dlx5 Factors on Engineered Stiff Hydrogels

Direct reprogramming of somatic cells into functional cells still faces major limitations in terms of efficiency and achieving functional maturity of the reprogramed cells. While different approaches have been developed commonly based on exploiting biochemical signals, introducing appropriate mechanical cues that stimulate the reprogramming process is rarely reported. In this study, collagen-coated polyacrylamide (PAM) hydrogels with stiffness close to that of collagenous bone (40 kPa) were adopted to augment the direct reprogramming process of mouse fibroblasts to osteoblastic-like cells. The results suggested that culturing cells on a hydrogel substrate enhanced the overexpression of osteogenic transcription factors using nonviral vectors and improved the yield of osteoblast-like cells. Particularly, a synergistic effect on achieving osteogenic functionality has been observed for the mechanical cues and overexpression of transcriptional factors, leading to enhanced osteogenic transformation and production of bone mineral matrix. Animal experiments suggested that reprogramed cells generated on matrix hydrogels accelerated bone regeneration and stimulated ectopic osteogenesis. Mechanism analysis suggested the critical involvement of actomyosin contraction and mechanical signal-mediated pathways like the RhoA-ROCK pathway, leading to a synergistic effect on the key transcriptional processes, including chromatin remodeling, nuclear translocation, and epigenetic transition. This study provides insights into the mechanical cue-enhanced direct reprogramming and cell therapy.

An alternative NURF complex sustains acute myeloid leukemia by regulating the accessibility of insulator regions

Efficient treatment of acute myeloid leukemia (AML) patients remains a challenge despite recent therapeutic advances. Here, using a CRISPRi screen targeting chromatin factors, we identified the nucleosome-remodeling factor (NURF) subunit BPTF as an essential regulator of AML cell survival. We demonstrate that BPTF forms an alternative NURF chromatin remodeling complex with SMARCA5 and BAP18, which regulates the accessibility of a large set of insulator regions in leukemic cells. This ensures efficient CTCF binding and boundary formation between topologically associated domains that is essential for maintaining the leukemic transcriptional programs. We also demonstrate that the well-studied PHD2-BROMO chromatin reader domains of BPTF, while contributing to complex recruitment to chromatin, are dispensable for leukemic cell growth. Taken together, our results uncover how the alternative NURF complex contributes to leukemia and provide a rationale for its targeting in AML.

Adhesome Receptor Clustering is Accompanied by the Colocalization of the Associated Genes in the Cell Nucleus

Proteins on the cell membrane cluster to respond to extracellular signals; for example, adhesion proteins cluster to enhance extracellular matrix sensing; or T-cell receptors cluster to enhance antigen sensing. Importantly, the maturation of such receptor clusters requires transcriptional control to adapt and reinforce the extracellular signal sensing. However, it has been unclear how such efficient clustering mechanisms are encoded at the level of the genes that code for these receptor proteins. Using the adhesome as an example, we show that genes that code for adhesome receptor proteins are spatially co-localized and co-regulated within the cell nucleus. Towards this, we use Hi-C maps combined with RNA-seq data of adherent cells to map the correspondence between adhesome receptor proteins and their associated genes. Interestingly, we find that the transcription factors that regulate these genes are also co-localized with the adhesome gene loci, thereby potentially facilitating a transcriptional reinforcement of the extracellular matrix sensing machinery. Collectively, our results highlight an important layer of transcriptional control of cellular signal sensing.

Cell-type-specific 3D-genome organization and transcription regulation in the brain

3D organization of the genome plays a critical role in regulating gene expression. However, it remains unclear how chromatin organization differs among different cell types in the brain. Here we used genome-scale DNA and RNA imaging to investigate 3D-genome organization in transcriptionally distinct cell types in the primary motor cortex of the mouse brain. We uncovered a wide spectrum of differences in the nuclear architecture and 3D-genome organization among different cell types, ranging from the physical size of the cell nucleus to the active-inactive chromatin compartmentalization and radial positioning of chromatin loci within the nucleus. These cell-type-dependent variations in nuclear architecture and chromatin organization exhibited strong correlation with both total transcriptional activity of the cell and transcriptional regulation of cell-type-specific marker genes. Moreover, we found that the methylated-DNA-binding protein MeCP2 regulates transcription in a divergent manner, depending on the nuclear radial positions of chromatin loci, through modulating active-inactive chromatin compartmentalization.

Transposon wave remodeled the epigenomic landscape in the rapid evolution of X-Chromosome dosage compensation

Sex chromosome dosage compensation is a model to understand the coordinated evolution of transcription; however, the advanced age of the sex chromosomes in model systems makes it difficult to study how the complex regulatory mechanisms underlying chromosome-wide dosage compensation can evolve. The sex chromosomes of Poecilia picta have undergone recent and rapid divergence, resulting in widespread gene loss on the male Y, coupled with complete X Chromosome dosage compensation, the first case reported in a fish. The recent de novo origin of dosage compensation presents a unique opportunity to understand the genetic and evolutionary basis of coordinated chromosomal gene regulation. By combining a new chromosome-level assembly of P. picta with whole-genome bisulfite sequencing and RNA-seq data, we determine that the YY1 transcription factor (YY1) DNA binding motif is associated with male-specific hypomethylated regions on the X, but not the autosomes. These YY1 motifs are the result of a recent and rapid repetitive element expansion on the P. picta X Chromosome, which is absent in closely related species that lack dosage compensation. Taken together, our results present compelling support that a disruptive wave of repetitive element insertions carrying YY1 motifs resulted in the remodeling of the X Chromosome epigenomic landscape and the rapid de novo origin of a dosage compensation system.

4D nucleome equation predicts gene expression controlled by long-range enhancer-promoter interaction

Recent experimental evidence strongly supports that three-dimensional (3D) long-range enhancer-promoter (E-P) interactions have important influences on gene-expression dynamics, but it is unclear how the interaction information is translated into gene expression over time (4D). To address this question, we developed a general theoretical framework (named as a 4D nucleome equation), which integrates E-P interactions on chromatin and biochemical reactions of gene transcription. With this equation, we first present the distribution of mRNA counts as a function of the E-P genomic distance and then reveal a power-law scaling of the expression level in this distance. Interestingly, we find that long-range E-P interactions can induce bimodal and trimodal mRNA distributions. The 4D nucleome equation also allows for model selection and parameter inference. When this equation is applied to the mouse embryonic stem cell smRNA-FISH data and the E-P genomic-distance data, the predicted E-P contact probability and mRNA distribution are in good agreement with experimental results. Further statistical inference indicates that the E-P interactions prefer to modulate the mRNA level by controlling promoter activation and transcription initiation rates. Our model and results provide quantitative insights into both spatiotemporal gene-expression determinants (i.e., long-range E-P interactions) and cellular fates during development.

Lineage specific 3D genome structure in the adult human brain and neurodevelopmental changes in the chromatin interactome

The human brain is a complex organ comprised of distinct cell types, and the contribution of the 3D genome to lineage specific gene expression remains poorly understood. To decipher cell type specific genome architecture, and characterize fine scale changes in the chromatin interactome across neural development, we compared the 3D genome of the human fetal cortical plate to that of neurons and glia isolated from the adult prefrontal cortex. We found that neurons have weaker genome compartmentalization compared to glia, but stronger TADs, which emerge during fetal development. Furthermore, relative to glia, the neuronal genome shifts more strongly towards repressive compartments. Neurons have differential TAD boundaries that are proximal to active promoters involved in neurodevelopmental processes. CRISPRi on CNTNAP2 in hIPSC-derived neurons reveals that transcriptional inactivation correlates with loss of insulation at the differential boundary. Finally, re-wiring of chromatin loops during neural development is associated with transcriptional and functional changes. Importantly, differential loops in the fetal cortex are associated with autism GWAS loci, suggesting a neuropsychiatric disease mechanism affecting the chromatin interactome. Furthermore, neural development involves gaining enhancer-promoter loops that upregulate genes that control synaptic activity. Altogether, our study provides multi-scale insights on the 3D genome in the human brain.

Dynamic RNA Polymerase II Recruitment Drives Differentiation of the Intestine under the direction of HNF4

Terminal differentiation requires a massive restructuring of the transcriptome. During intestinal differentiation, the expression patterns of nearly 4000 genes are altered as cells transition from progenitor cells in crypts to differentiated cells in villi. We identified dynamic recruitment of RNA Polymerase II (Pol II) to gene promoters as the primary driver of transcriptomic shifts during intestinal differentiation in vivo. Changes in enhancer-promoter looping interactions accompany dynamic Pol II recruitment and are dependent upon HNF4, a pro-differentiation transcription factor. Using genetic loss-of- function, ChIP-seq and IP mass spectrometry, we demonstrate that HNF4 collaborates with chromatin remodelers and loop-stabilizing proteins and facilitates Pol II recruitment at hundreds of genes pivotal to differentiation. We also explore alternate mechanisms which drive differentiation gene expression and find pause-release of Pol II and post- transcriptional mRNA stability regulate smaller subsets of differentially expressed genes. These studies provide insights into the mechanisms of differentiation in a renewing adult tissue.