Polymer Modeling Predicts Chromosome Reorganization in Senescence

Combinatorics and topological weights of chromatin loop networks

Polymer physics models suggest that chromatin spontaneously folds into loop networks with transcription units (TUs), such as enhancers and promoters, as anchors. Here we use combinatoric arguments to enumerate the emergent chromatin loop networks, both in the case where TUs are labeled and where they are unlabeled. We then combine these mathematical results with those of computer simulations aimed at finding the inter-TU energy required to form a target loop network. We show that different topologies are vastly different in terms of both their combinatorial weight and energy of formation. We explain the latter result qualitatively by computing the topological weight of a given network-i.e., its partition function in statistical mechanics language-in the approximation where excluded volume interactions are neglected. Our results show that networks featuring local loops are statistically more likely with respect to networks including more nonlocal contacts. We suggest our classification of loop networks, together with our estimate of the combinatorial and topological weight of each network, will be relevant to catalog three-dimensional structures of chromatin fibers around eukaryotic genes, and to estimate their relative frequency in both simulations and experiments.

The nexus of nuclear envelope dynamics, circular economy and cancer cell pathophysiology

The nuclear envelope (NE) is a critical component in maintaining the function and structure of the eukaryotic nucleus. The NE and lamina are disassembled during each cell cycle to enable an open mitosis. Nuclear architecture construction and deconstruction is a prime example of a circular economy, as it fulfills a highly efficient recycling program bound to continuous assessment of the quality and functionality of the building blocks. Alterations in the nuclear dynamics and lamina structure have emerged as important contributors to both oncogenic transformation and cancer progression. However, the knowledge of the NE breakdown and reassembly is still limited to a fraction of participating proteins and complexes. As cancer cells contain highly diverse nuclei in terms of DNA content, but also in terms of nuclear number, size, and shape, it is of great interest to understand the intricate relationship between these nuclear features in cancer cell pathophysiology. In this review, we provide insights into how those NE dynamics are regulated, and how lamina destabilization processes may alter the NE circular economy. Moreover, we expand the knowledge of the lamina-associated domain region by using strategic algorithms, including Artificial Intelligence, to infer protein associations, assess their function and location, and predict cancer-type specificity with implications for the future of cancer diagnosis, prognosis and treatment. Using this approach we identified NUP98 and MECP2 as potential proteins that exhibit upregulation in Acute Myeloid Leukemia (LAML) patients with implications for early diagnosis.

Mechanism of phase condensation for chromosome architecture and function

Chromosomal phase separation is involved in a broad spectrum of chromosome organization and functional processes. Nonetheless, the intricacy of this process has left its molecular mechanism unclear. Here, we introduce the principles governing phase separation and its connections to physiological roles in this context. Our primary focus is contrasting two phase separation mechanisms: self-association-induced phase separation (SIPS) and bridging-induced phase separation (BIPS). We provide a comprehensive discussion of the distinct features characterizing these mechanisms and offer illustrative examples that suggest their broad applicability. With a detailed understanding of these mechanisms, we explore their associations with nucleosomes and chromosomal biological functions. This comprehensive review contributes to the exploration of uncharted territory in the intricate interplay between chromosome architecture and function.

Choreography of lamina-associated domains: structure meets dynamics

Lamina-associated domains are large regions of heterochromatin positioned at the nuclear periphery. These domains have been implicated in gene repression, especially in the context of development. In mammals, LAD organization is dependent on nuclear lamins, inner nuclear membrane proteins, and chromatin state. In addition, chromatin readers and modifier proteins have been implicated in this organization, potentially serving as molecular tethers that interact with both nuclear envelope proteins and chromatin. More recent studies have focused on teasing apart the rules that govern dynamic LAD organization and how LAD organization, in turn, relates to gene regulation and overall 3D genome organization. This review highlights recent studies in mammalian cells uncovering factors that instruct the choreography of LAD organization, re-organization, and dynamics at the nuclear lamina, including LAD dynamics in interphase and through mitotic exit, when LAD organization is re-established, as well as intra-LAD subdomain variations.

Transcription modulates chromatin dynamics and locus configuration sampling

In living cells, the 3D structure of gene loci is dynamic, but this is not revealed by 3C and FISH experiments in fixed samples, leaving a notable gap in our understanding. To overcome these limitations, we applied the highly predictive heteromorphic polymer (HiP-HoP) model to determine chromatin fiber mobility at the Pax6 locus in three mouse cell lines with different transcription states. While transcriptional activity minimally affects movement of 40-kbp regions, we observed that motion of smaller 1-kbp regions depends strongly on local disruption to chromatin fiber structure marked by H3K27 acetylation. This also substantially influenced locus configuration dynamics by modulating protein-mediated promoter-enhancer loops. Importantly, these simulations indicate that chromatin dynamics are sufficiently fast to sample all possible locus conformations within minutes, generating wide dynamic variability within single cells. This combination of simulation and experimental validation provides insight into how transcriptional activity influences chromatin structure and gene dynamics.

Dynamic 3D genome reorganization during senescence: defining cell states through chromatin

Cellular senescence, a cell state characterized by growth arrest and insensitivity to growth stimulatory hormones, is accompanied by a massive change in chromatin organization. Senescence can be induced by a range of physiological signals and pathological stresses and was originally thought to be an irreversible state, implicated in normal development, wound healing, tumor suppression and aging. Recently cellular senescence was shown to be reversible in some cases, with exit being triggered by the modulation of the cell's transcriptional program by the four Yamanaka factors, the suppression of p53 or H3K9me3, PDK1, and/or depletion of AP-1. Coincident with senescence reversal are changes in chromatin organization, most notably the loss of senescence-associated heterochromatin foci (SAHF) found in oncogene-induced senescence. In addition to fixed-cell imaging, chromatin conformation capture and multi-omics have been used to examine chromatin reorganization at different spatial resolutions during senescence. They identify determinants of SAHF formation and other key features that differentiate distinct types of senescence. Not surprisingly, multiple factors, including the time of induction, the type of stress experienced, and the type of cell involved, influence the global reorganization of chromatin in senescence. Here we discuss how changes in the three-dimensional organization of the genome contribute to the regulation of transcription at different stages of senescence. In particular, the distinct contributions of heterochromatin- and lamina-mediated interactions, changes in gene expression, and other cellular control mechanisms are discussed. We propose that high-resolution temporal and spatial analyses of the chromatin landscape during senescence will identify early markers of the different senescence states to help guide clinical diagnosis.

Beyond SAHF: An integrative view of chromatin compartmentalization during senescence

Cellular senescence, a persistent form of cell cycle arrest, has been linked to the formation of heterochromatic foci, accompanied by additional concentric epigenetic layers. However, senescence is a highly heterogeneous phenotype, and the formation of these structures is context dependent. Recent developments in the understanding of the high-order chromatin organization have opened new avenues for contextualizing the nuclear and chromatin phenotypes of senescence. Oncogene-induced senescence displays prominent foci and typically exhibits increased chromatin compartmentalization, based on the chromosome conformation assays, as marked by increased transcompaction and segregation of the heterochromatin and euchromatin. However, other types of senescence (e.g., replicative senescence) exhibit comparatively lower levels of compartmentalization. Thus, a more integrative view of the global rearrangement of the chromatin architecture that occurs during senescence is emerging, with potential functional implications for the heterogeneity of the senescence phenotype.

Strong interactions between highly dynamic lamina-associated domains and the nuclear envelope stabilize the 3D architecture of Drosophila interphase chromatin

BACKGROUND

Interactions among topologically associating domains (TADs), and between the nuclear envelope (NE) and lamina-associated domains (LADs) are expected to shape various aspects of three-dimensional (3D) chromatin structure and dynamics; however, relevant genome-wide experiments that may provide statistically significant conclusions remain difficult.

RESULTS

We have developed a coarse-grained dynamical model of D. melanogaster nuclei at TAD resolution that explicitly accounts for four distinct epigenetic classes of TADs and LAD-NE interactions. The model is parameterized to reproduce the experimental Hi-C map of the wild type (WT) nuclei; it describes time evolution of the chromatin over the G1 phase of the interphase. The simulations include an ensemble of nuclei, corresponding to the experimentally observed set of several possible mutual arrangements of chromosomal arms. The model is validated against multiple structural features of chromatin from several different experiments not used in model development. Predicted positioning of all LADs at the NE is highly dynamic-the same LAD can attach, detach and move far away from the NE multiple times during interphase. The probabilities of LADs to be in contact with the NE vary by an order of magnitude, despite all having the same affinity to the NE in the model. These probabilities are mostly determined by a highly variable local linear density of LADs along the genome, which also has the same strong effect on the predicted positioning of individual TADs -- higher probability of a TAD to be near NE is largely determined by a higher linear density of LADs surrounding this TAD. The distribution of LADs along the chromosome chains plays a notable role in maintaining a non-random average global structure of chromatin. Relatively high affinity of LADs to the NE in the WT nuclei substantially reduces sensitivity of the global radial chromatin distribution to variations in the strength of TAD-TAD interactions compared to the lamin depleted nuclei, where a small (0.5 kT) increase of cross-type TAD-TAD interactions doubles the chromatin density in the central nucleus region.

CONCLUSIONS

A dynamical model of the entire fruit fly genome makes multiple genome-wide predictions of biological interest. The distribution of LADs along the chromatin chains affects their probabilities to be in contact with the NE and radial positioning of highly mobile TADs, playing a notable role in creating a non-random average global structure of the chromatin. We conjecture that an important role of attractive LAD-NE interactions is to stabilize global chromatin structure against inevitable cell-to-cell variations in TAD-TAD interactions.

Mesoscale, long-time mixing of chromosomes and its connection to polymer dynamics

Chromosomes are arranged in distinct territories within the nucleus of animal cells. Recent experiments have shown that these territories overlap at their edges, suggesting partial mixing during interphase. Experiments that knock-down of condensin II proteins during interphase indicate increased chromosome mixing, which demonstrates control of the mixing. In this study, we use a generic polymer simulation to quantify the dynamics of chromosome mixing over time. We introduce the chromosome mixing index, which quantifies the mixing of distinct chromosomes in the nucleus. We find that the chromosome mixing index in a small confinement volume (as a model of the nucleus), increases as a power-law of the time, with the scaling exponent varying non-monotonically with self-interaction and volume fraction. By comparing the chromosome mixing index with both monomer subdiffusion due to (non-topological) intermingling of chromosomes as well as even slower reptation, we show that for relatively large volume fractions, the scaling exponent of the chromosome mixing index is related to Rouse dynamics for relatively weak chromosome attractions and to reptation for strong attractions. In addition, we extend our model to more realistically account for the situation of the Drosophila chromosome by including the heterogeneity of the polymers and their lengths to account for microphase separation of euchromatin and heterochromatin and their interactions with the nuclear lamina. We find that the interaction with the lamina further impedes chromosome mixing.

Phase transitions in heterochromatin organization

Heterochromatin formation has been proposed to involve phase transitions on the level of the three-dimensional folding of heterochromatin regions and the liquid-liquid demixing of heterochromatin proteins. Here, I outline the hallmarks of such transitions and the current challenges to detect them in living cells. I further discuss the abundance and properties of prominent heterochromatin proteins and relate them to their potential role in driving phase transitions. Recent data from mouse fibroblasts indicate that pericentric heterochromatin is organized via a reordering transition on the level of heterochromatin regions that does not necessarily involve liquid-liquid demixing of heterochromatin proteins. Evaluating key hallmarks of the different candidate phase transition mechanisms across cell types and species will be needed to complete the current picture.

Compartmentalization with nuclear landmarks yields random, yet precise, genome organization

The 3D organization of eukaryotic genomes plays an important role in genome function. While significant progress has been made in deciphering the folding mechanisms of individual chromosomes, the principles of the dynamic large-scale spatial arrangement of all chromosomes inside the nucleus are poorly understood. We use polymer simulations to model the diploid human genome compartmentalization relative to nuclear bodies such as nuclear lamina, nucleoli, and speckles. We show that a self-organization process based on a cophase separation between chromosomes and nuclear bodies can capture various features of genome organization, including the formation of chromosome territories, phase separation of A/B compartments, and the liquid property of nuclear bodies. The simulated 3D structures quantitatively reproduce both sequencing-based genomic mapping and imaging assays that probe chromatin interaction with nuclear bodies. Importantly, our model captures the heterogeneous distribution of chromosome positioning across cells while simultaneously producing well-defined distances between active chromatin and nuclear speckles. Such heterogeneity and preciseness of genome organization can coexist due to the nonspecificity of phase separation and the slow chromosome dynamics. Together, our work reveals that the cophase separation provides a robust mechanism for us to produce functionally important 3D contacts without requiring thermodynamic equilibration that can be difficult to achieve.

4D epigenomics: deciphering the coupling between genome folding and epigenomic regulation with biophysical modeling

Recent experimental observations suggest a strong coupling between the 3D nuclear chromosome organization and epigenomics. However, the mechanistic and functional bases of such interplay remain elusive. In this review, we describe how biophysical modeling has been instrumental in characterizing how genome folding may impact the formation of epigenomic domains and, conversely, how epigenomic marks may affect chromosome conformation. Finally, we discuss how this mutual feedback loop between chromatin organization and epigenome regulation, via the formation of physicochemical nanoreactors, may represent a key functional role of 3D compartmentalization in the assembly and maintenance of stable - but yet plastic - epigenomic landscapes.

Predicting Adhesive Free Energies of Polymer-Surface Interactions with Machine Learning

Polymer-surface interactions are crucial to many biological processes and industrial applications. Here we propose a machine learning method to connect a model polymer's sequence with its adhesion to decorated surfaces. We simulate the adhesive free energies of 20000 unique coarse-grained one-dimensional polymer sequences interacting with functionalized surfaces and build support vector regression models that demonstrate inexpensive and reliable prediction of the adhesive free energy as a function of sequence. Our work highlights the promising integration of coarse-grained simulation with data-driven machine learning methods for the design of functional polymers and represents an important step toward linking polymer compositions with polymer-surface interactions.

Characterizing the variation in chromosome structure ensembles in the context of the nuclear microenvironment

Inside the nucleus, chromosomes are subjected to direct physical interaction between different components, active forces, and thermal noise, leading to the formation of an ensemble of three-dimensional structures. However, it is still not well understood to what extent and how the structural ensemble varies from one chromosome region or cell-type to another. We designed a statistical analysis technique and applied it to single-cell chromosome imaging data to reveal the heterogeneity of individual chromosome structures. By analyzing the resulting structural landscape, we find that the largest dynamic variation is the overall radius of gyration of the chromatin region, followed by domain reorganization within the region. By comparing different human cell-lines and experimental perturbation data using this statistical analysis technique and a network-based similarity quantification approach, we identify both cell-type and condition-specific features of the structural landscapes. We identify a relationship between epigenetic state and the properties of chromosome structure fluctuation and validate this relationship through polymer simulations. Overall, our study suggests that the types of variation in a chromosome structure ensemble are cell-type as well as region-specific and can be attributed to constraints placed on the structure by factors such as variation in epigenetic state.

Recent work has revealed principles of how chromosomes are folded and structured inside the human nucleus. It is now even possible to microscopically trace the path of chromosomes in 3D in individual cells. With this data, we can start to examine how much variation exists in chromosome structure and what biological factors may restrict or enhance this variation. Are chromosomes stuck in just a few possible positions or do they move around more freely, sampling many configurations? Here, we use a mathematical approach to compare chromosome structure variation in different cell types, at different locations along the genome, and when key structural proteins are removed. Through these comparisons and dynamic simulations of chromosome behavior, we identify factors that may constrain or promote variation in chromosome structure.

Phase Separation and Correlated Motions in Motorized Genome

The human genome is arranged in the cell nucleus nonrandomly, and phase separation has been proposed as an important driving force for genome organization. However, the cell nucleus is an active system, and the contribution of nonequilibrium activities to phase separation and genome structure and dynamics remains to be explored. We simulated the genome using an energy function parametrized with chromosome conformation capture (Hi-C) data with the presence of active, nondirectional forces that break the detailed balance. We found that active forces that may arise from transcription and chromatin remodeling can dramatically impact the spatial localization of heterochromatin. When applied to euchromatin, active forces can drive heterochromatin to the nuclear envelope and compete with passive interactions among heterochromatin that tend to pull them in opposite directions. Furthermore, active forces induce long-range spatial correlations among genomic loci beyond single chromosome territories. We further showed that the impact of active forces could be understood from the effective temperature defined as the fluctuation-dissipation ratio. Our study suggests that nonequilibrium activities can significantly impact genome structure and dynamics, producing unexpected collective phenomena.

Restructuring of Lamina-Associated Domains in Senescence and Cancer

Induction of cellular senescence or cancer is associated with a reshaping of the nuclear envelope and a broad reorganization of heterochromatin. At the periphery of mammalian nuclei, heterochromatin is stabilized at the nuclear lamina via lamina-associated domains (LADs). Alterations in the composition of the nuclear lamina during senescence lead to a loss of peripheral heterochromatin, repositioning of LADs, and changes in epigenetic states of LADs. Cancer initiation and progression are also accompanied by a massive reprogramming of the epigenome, particularly in domains coinciding with LADs. Here, we review recent knowledge on alterations in chromatin organization and in the epigenome that affect LADs and related genomic domains in senescence and cancer.

Deciphering aging at three-dimensional genomic resolution

Aging is characterized by progressive functional declines at the organismal, organic, and cellular levels and increased susceptibility to aging-related diseases. Epigenetic alteration is a hallmark of aging, senescent cells show epigenomic changes at multiple scales, such as 3D genome reorganization, alterations of histone modifications and chromatin accessibility, and DNA hypomethylation. Chromosome conformation capture (3C)-based technologies have enabled the generation of key information on genomic reorganizations during senescence. A comprehensive understanding of epigenomic alterations during aging will yield important insights into the underlying epigenetic mechanism for aging regulation, the identification of aging-related biomarkers, and the development of potential aging intervention targets.

Biology and Model Predictions of the Dynamics and Heterogeneity of Chromatin-Nuclear Lamina Interactions

Associations of chromatin with the nuclear lamina, at the nuclear periphery, help shape the genome in 3 dimensions. The genomic landscape of lamina-associated domains (LADs) is well characterized, but much remains unknown on the physical and mechanistic properties of chromatin conformation at the nuclear lamina. Computational models of chromatin folding at, and interactions with, a surface representing the nuclear lamina are emerging in attempts to characterize these properties and predict chromatin behavior at the lamina in health and disease. Here, we highlight the heterogeneous nature of the nuclear lamina and LADs, outline the main 3-dimensional chromatin structural modeling methods, review applications of modeling chromatin-lamina interactions and discuss biological insights inferred from these models in normal and disease states. Lastly, we address perspectives on future developments in modeling chromatin interactions with the nuclear lamina.

Large-scale chromatin reorganization reactivates placenta-specific genes that drive cellular aging

Nuclear deformation, a hallmark frequently observed in senescent cells, is presumed to be associated with the erosion of chromatin organization at the nuclear periphery. However, how such gradual changes in higher-order genome organization impinge on local epigenetic modifications to drive cellular mechanisms of aging has remained enigmatic. Here, through large-scale epigenomic analyses of isogenic young, senescent, and progeroid human mesenchymal progenitor cells (hMPCs), we delineate a hierarchy of integrated structural state changes that manifest as heterochromatin loss in repressive compartments, euchromatin weakening in active compartments, switching in interfacing topological compartments, and increasing epigenetic entropy. We found that the epigenetic de-repression unlocks the expression of pregnancy-specific beta-1 glycoprotein (PSG) genes that exacerbate hMPC aging and serve as potential aging biomarkers. Our analyses provide a rich resource for uncovering the principles of epigenomic landscape organization and its changes in cellular aging and for identifying aging drivers and intervention targets with a genome-topology-based mechanism.

Modeling the 3D Genome Using Hi-C and Nuclear Lamin-Genome Contacts

The three-dimensional (3D) organization of the genome is shaped by interactions with multiple structures within the nucleus, affecting gene expression outcomes. Technological breakthroughs in recent years have generated vast data reflecting various aspects of nuclear genome architecture in space and time. Integrating these datasets into comprehensive 3D genome models can reveal new insights into genome structure and regulation in normal and disease states. In this chapter, we provide a step-by-step guide on how to generate publication-ready integrated 3D genome models from (raw) Hi-C data and from lamin-genome (LAD) contacts.

The Physical Behavior of Interphase Chromosomes: Polymer Theory and Coarse-Grain Computer Simulations

Fluorescence in situ hybridization and chromosome conformation capture methods point to the same conclusion: that chromosomes appear to the external observer as compact structures with a highly nonrandom three-dimensional organization. In this work, we recapitulate the efforts made by us and other groups to rationalize this behavior in terms of the mathematical language and tools of polymer physics. After a brief introduction dedicated to some crucial experiments dissecting the structure of interphase chromosomes, we discuss at a nonspecialistic level some fundamental aspects of theoretical and numerical polymer physics. Then, we inglobe biological and polymer aspects into a polymer model for interphase chromosomes which moves from the observation that mutual topological constraints, such as those typically present between polymer chains in ordinary melts, induce slow chain dynamics and "constraint" chromosomes to resemble double-folded randomly branched polymer conformations. By explicitly turning these ideas into a multi-scale numerical algorithm which is described here in full details, we can design accurate model polymer conformations for interphase chromosomes and offer them for systematic comparison to experiments. The review is concluded by discussing the limitations of our approach and pointing to promising perspectives for future work.

Predicting genome organisation and function with mechanistic modelling

Fitting-free mechanistic models based on polymer simulations predict chromatin folding in 3D by focussing on the underlying biophysical mechanisms. This class of models has been increasingly used in conjunction with experiments to study the spatial organisation of eukaryotic chromosomes. Feedback from experiments to models leads to successive model refinement and has previously led to the discovery of new principles for genome organisation. Here, we review the basis of mechanistic polymer simulations, explain some of the more recent approaches and the contexts in which they have been useful to explain chromosome biology, and speculate on how they might be used in the future.

The interplay of chromatin phase separation and lamina interactions in nuclear organization

The genetic material of eukaryotes is segregated into transcriptionally active euchromatin and silent heterochromatin compartments. The spatial arrangement of chromatin compartments evolves over the course of cellular life in a process that remains poorly understood. The latest nuclear imaging experiments reveal a number of dynamical signatures of chromatin that are reminiscent of active multiphase liquids. This includes the observations of viscoelastic response, coherent motions, Ostwald ripening, and coalescence of chromatin compartments. There is also growing evidence that liquid-liquid phase separation of protein and nucleic acid components is the underlying mechanism for the dynamical behavior of chromatin. To dissect the organizational and dynamical implications of chromatin's liquid behavior, we have devised a phenomenological field-theoretic model of the nucleus as a multiphase condensate of liquid chromatin types. Employing the liquid chromatin model of the Drosophila nucleus, we have carried out an extensive set of simulations with an objective to shed light on the dynamics and chromatin patterning observed in the latest nuclear imaging experiments. Our simulations reveal the emergence of experimentally detected mesoscale chromatin channels and spheroidal droplets which arise from the dynamic interplay of chromatin type to type interactions and intermingling of chromosomal territories. We also quantitatively reproduce coherent motions of chromatin domains observed in displacement correlation spectroscopy measurements which are explained within the framework of our model by phase separation of chromatin types operating within constrained intrachromosomal and interchromosomal boundaries. Finally, we illuminate the role of heterochromatin-lamina interactions in the nuclear organization by showing that these interactions enhance the mobility of euchromatin and indirectly introduce correlated motions of heterochromatin droplets.

HMGs as rheostats of chromosomal structure and cell proliferation

High mobility group proteins (HMGs) are the most abundant nuclear proteins next to histones and are robustly expressed across tissues and organs. HMGs can uniquely bend or bind distorted DNA, and are central to such processes as transcription, recombination, and DNA repair. However, their dynamic association with chromatin renders capturing HMGs on chromosomes challenging. Recent work has changed this and now implicates these factors in spatial genome organization. Here, I revisit older and review recent literature to describe how HMGs rewire spatial chromatin interactions to sustain homeostasis or promote cellular aging. I propose a 'rheostat' model to explain how HMG-box proteins (HMGBs), and to some extent HMG A proteins (HMGAs), may control cellular aging and, likely, cancer progression.

Reduction of lamin B receptor levels by miR-340-5p disrupts chromatin, promotes cell senescence and enhances senolysis

A major stress response influenced by microRNAs (miRNAs) is senescence, a state of indefinite growth arrest triggered by sublethal cell damage. Here, through bioinformatic analysis and experimental validation, we identified miR-340-5p as a novel miRNA that foments cellular senescence. miR-340-5p was highly abundant in diverse senescence models, and miR-340-5p overexpression in proliferating cells rendered them senescent. Among the target mRNAs, miR-340-5p prominently reduced the levels of LBR mRNA, encoding lamin B receptor (LBR). Loss of LBR by ectopic overexpression of miR-340-5p derepressed heterochromatin in lamina-associated domains, promoting the expression of DNA repetitive elements characteristic of senescence. Importantly, overexpressing miR-340-5p enhanced cellular sensitivity to senolytic compounds, while antagonization of miR-340-5p reduced senescent cell markers and engendered resistance to senolytic-induced cell death. We propose that miR-340-5p can be exploited for removing senescent cells to restore tissue homeostasis and mitigate damage by senescent cells in pathologies of human aging.

Multiscale modeling of genome organization with maximum entropy optimization

Three-dimensional (3D) organization of the human genome plays an essential role in all DNA-templated processes, including gene transcription, gene regulation, and DNA replication. Computational modeling can be an effective way of building high-resolution genome structures and improving our understanding of these molecular processes. However, it faces significant challenges as the human genome consists of over 6 × 109 base pairs, a system size that exceeds the capacity of traditional modeling approaches. In this perspective, we review the progress that has been made in modeling the human genome. Coarse-grained models parameterized to reproduce experimental data via the maximum entropy optimization algorithm serve as effective means to study genome organization at various length scales. They have provided insight into the principles of whole-genome organization and enabled de novo predictions of chromosome structures from epigenetic modifications. Applications of these models at a near-atomistic resolution further revealed physicochemical interactions that drive the phase separation of disordered proteins and dictate chromatin stability in situ. We conclude with an outlook on the opportunities and challenges in studying chromosome dynamics.

The loss of heterochromatin is associated with multiscale three-dimensional genome reorganization and aberrant transcription during cellular senescence

Heterochromatin remodeling is critical for various cell processes. In particular, the "loss of heterochromatin" phenotype in cellular senescence is associated with the process of aging and age-related disorders. Although biological processes of senescent cells, including senescence-associated heterochromatin foci (SAHF) formation, chromosome compaction, and redistribution of key proteins, have been closely associated with high-order chromatin structure, the relationship between the high-order chromatin reorganization and the loss of heterochromatin phenotype during senescence has not been fully understood. By using senescent and deep senescent fibroblasts induced by DNA damage harboring the "loss of heterochromatin" phenotype, we observed progressive 3D reorganization of heterochromatin during senescence. Facultative and constitutive heterochromatin marked by H3K27me3 and H3K9me3, respectively, show different alterations. Facultative heterochromatin tends to switch from the repressive B-compartment to the active A-compartment, whereas constitutive heterochromatin shows no significant changes at the compartment level but enhanced interactions between themselves. Both types of heterochromatin show increased chromatin accessibility and gene expression leakage during senescence. Furthermore, increased chromatin accessibility in potential CTCF binding sites accompanies the establishment of novel loops in constitutive heterochromatin. Finally, we also observed aberrant expression of repetitive elements, including LTR (long terminal repeat) and satellite classes. Overall, facultative and constitutive heterochromatin show both similar and distinct multiscale alterations in the 3D map, chromatin accessibility, and gene expression leakage. This study provides an epigenomic map of heterochromatin reorganization during senescence.

Principles of 3D compartmentalization of the human genome

Chromatin is organized in the nucleus via CTCF loops and compartmental domains. Here, we compare different cell types to identify distinct paradigms of compartmental domain formation in human tissues. We identify and quantify compartmental forces correlated with histone modifications characteristic of transcriptional activity and previously underappreciated roles for distinct compartmental domains correlated with the presence of H3K27me3 and H3K9me3, respectively. We present a computer simulation model capable of predicting compartmental organization based on the biochemical characteristics of independent chromatin features. Using this model, we show that the underlying forces responsible for compartmental domain formation in human cells are conserved and that the diverse compartmentalization patterns seen across cell types are due to differences in chromatin features. We extend these findings to Drosophila to suggest that the same principles are at work beyond humans. These results offer mechanistic insights into the fundamental forces driving the 3D organization of the genome.

Using high-resolution Hi-C data and computer simulations, Nichols and Corces show that compartments arise as a consequence of interactions among proteins that correlate with the presence of H3K27ac, H3K27me3, and H3K9me3, suggesting that human cells contain at least three distinct compartments. The same principles apply to other organisms.

Interphase Chromosomes in Replicative Senescence: Chromosome Positioning as a Senescence Biomarker and the Lack of Nuclear Motor-Driven Chromosome Repositioning in Senescent Cells

This study demonstrates, and confirms, that chromosome territory positioning is altered in primary senescent human dermal fibroblasts (HDFs). The chromosome territory positioning pattern is very similar to that found in HDFs made quiescent either by serum starvation or confluence; but not completely. A few chromosomes are found in different locations. One chromosome in particular stands out, chromosome 10, which is located in an intermediate location in young proliferating HDFs, but is found at the nuclear periphery in quiescent cells and in an opposing location of the nuclear interior in senescent HDFs. We have previously demonstrated that individual chromosome territories can be actively and rapidly relocated, with 15 min, after removal of serum from the culture media. These chromosome relocations require nuclear motor activity through the presence of nuclear myosin 1β (NM1β). We now also demonstrate rapid chromosome movement in HDFs after heat-shock at 42°C. Others have shown that heat shock genes are actively relocated using nuclear motor protein activity via actin or NM1β (Khanna et al., 2014; Pradhan et al., 2020). However, this current study reveals, that in senescent HDFs, chromosomes can no longer be relocated to expected nuclear locations upon these two types of stimuli. This coincides with a entirely different organisation and distribution of NM1β within senescent HDFs.

Mesoscale phase separation of chromatin in the nucleus

Intact-organism imaging of Drosophila larvae reveals and quantifies chromatin-aqueous phase separation. The chromatin can be organized near the lamina layer of the nuclear envelope, conventionally fill the nucleus, be organized centrally, or as a wetting droplet. These transitions are controlled by changes in nuclear volume and the interaction of chromatin with the lamina (part of the nuclear envelope) at the nuclear periphery. Using a simple polymeric model that includes the key features of chromatin self-attraction and its binding to the lamina, we demonstrate theoretically that it is the competition of these two effects that determines the mode of chromatin distribution. The qualitative trends as well as the composition profiles obtained in our simulations compare well with the observed intact-organism imaging and quantification. Since the simulations contain only a small number of physical variables we can identify the generic mechanisms underlying the changes in the observed phase separations.

The functional impact of nuclear reorganization in cellular senescence

Cellular senescence is the irreversible cell cycle arrest in response to DNA damage. Because senescent cells accumulate with age and contribute to chronic inflammation, they are promising therapeutic targets for healthspan extension. The senescent phenotype can vary depending on cell type and on the specific insults that induce senescence. This variability is also reflected in the extensive remodeling of the genome organization within the nucleus of senescent cells. Here, we give an overview of the nuclear changes that occur in different forms of senescence, including changes to chromatin state and composition and to the three-dimensional organization of the genome, as well as alterations to the nuclear envelope and to the accessibility of repetitive genomic regions. Many of these changes are shared across all forms of senescence, implicating nuclear organization as a fundamental driver of the senescent state and of how senescent cells interact with the surrounding tissue.

A DNA Sequence Based Polymer Model for Chromatin Folding

The recent development of sequencing technology and imaging methods has provided an unprecedented understanding of the inter-phase chromatin folding in mammalian nuclei. It was found that chromatin folds into topological-associated domains (TADs) of hundreds of kilo base pairs (kbps), and is further divided into spatially segregated compartments (A and B). The compartment B tends to be located near to the periphery or the nuclear center and interacts with other domains of compartments B, while compartment A tends to be located between compartment B and interacts inside the domains. These spatial domains are found to highly correlate with the mosaic CpG island (CGI) density. High CGI density corresponds to compartments A and small TADs, and vice versa. The variation of contact probability as a function of sequential distance roughly follows a power-law decay. Different chromosomes tend to segregate to occupy different chromosome territories. A model that can integrate these properties at multiple length scales and match many aspects is highly desired. Here, we report a DNA-sequence based coarse-grained block copolymer model that considers different interactions between blocks of different CGI density, interactions of TAD formation, as well as interactions between chromatin and the nuclear envelope. This model captures the various single-chromosome properties and partially reproduces the formation of chromosome territories.

Predicting Genome Architecture: Challenges and Solutions

Genome architecture plays a pivotal role in gene regulation. The use of high-throughput methods for chromatin profiling and 3-D interaction mapping provide rich experimental data sets describing genome organization and dynamics. These data challenge development of new models and algorithms connecting genome architecture with epigenetic marks. In this review, we describe how chromatin architecture could be reconstructed from epigenetic data using biophysical or statistical approaches. We discuss the applicability and limitations of these methods for understanding the mechanisms of chromatin organization. We also highlight the emergence of new predictive approaches for scoring effects of structural variations in human cells.

Physical constraints in polymer modeling of chromatin associations with the nuclear periphery at kilobase scale

Interactions of chromatin with the nuclear lamina imposes a radial genome distribution important for nuclear functions. How physical properties of chromatin affect these interactions is unclear. We used polymer simulations to model how physical parameters of chromatin affect its interaction with the lamina. Impact of polymer stiffness is greater than stretching on its configurations at the lamina; these are manifested as trains describing extended interactions, and loops describing desorbed regions . Conferring an attraction potential leads to persistent interaction and adsorption-desorption regimes manifested by fluctuations between trains and loops. These are modulated by polymer stiffness and stretching, with a dominant impact of stiffness on resulting structural configurations. We infer that flexible euchromatin is more prone to stochastic interactions with lamins than rigid heterochromatin characterizing constitutive LADs. Our models provide insights on the physical properties of chromatin as a polymer which affect the dynamics and patterns of interactions with the nuclear lamina.

Reflections on the organization and the physical state of chromatin in eukaryotic cells

In recent years, our perception of chromatin structure and organization in the cell nucleus has changed in fundamental ways. The 30 nm chromatin fiber has lost its status as an essential in vivo structure. Hi-C and related biochemical methods, advanced electron and super-resolved fluorescence microscopy, together with concepts from soft matter physics, have revolutionized the field. A comprehensive understanding of the structural and functional interactions that regulate cell cycle and cell type specific nuclear functions appears within reach, but it requires the integration of top-down and bottom-up approachs. In this review, I present an update on nuclear architecture studies with an emphasis on organization and the controversy regarding the physical state of chromatin in cells.

Heterochromatin as an Important Driver of Genome Organization

Heterochromatin is a constituent of eukaryotic genomes with functions spanning from gene expression silencing to constraining DNA replication and repair. Inside the nucleus, heterochromatin segregates spatially from euchromatin and is localized preferentially toward the nuclear periphery and surrounding the nucleolus. Despite being an abundant nuclear compartment, little is known about how heterochromatin regulates and participates in the mechanisms driving genome organization. Here, we review pioneer and recent evidence that explores the functional role of heterochromatin in the formation of distinct chromatin compartments and how failure of the molecular mechanisms forming heterochromatin leads to disarray of genome conformation and disease.

Genome folding and refolding in differentiation and cellular senescence

The spatial conformation of chromatin within the confines of eukaryotic cell nuclei is now acknowledged as a decisive epigenetic mechanism for the modulation of such cellular functions as gene expression regulation, DNA replication or DNA damage repair. Of course, these processes are tightly regulated during organismal development and markedly affected by cellular ageing. Thus, the question that arises is to what extent does folding or refolding of the genome in three-dimensional space underlie the progression of development or ageing? Herein, we discuss recent experimental and modelling evidence to address this question and revisit how these seemingly different processed might represent two sides of the same coin.

Multiplexed imaging of nucleome architectures in single cells of mammalian tissue

The three-dimensional architecture of the genome affects genomic functions. Multiple genome architectures at different length scales, including chromatin loops, domains, compartments, and lamina- and nucleolus-associated regions, have been discovered. However, how these structures are arranged in the same cell and how they are mutually correlated in different cell types in mammalian tissue are largely unknown. Here, we develop Multiplexed Imaging of Nucleome Architectures that measures multiscale chromatin folding, copy numbers of numerous RNA species, and associations of numerous genomic regions with nuclear lamina, nucleoli and surface of chromosomes in the same, single cells. We apply this method in mouse fetal liver, and identify de novo cell-type-specific chromatin architectures associated with gene expression, as well as cell-type-independent principles of chromatin organization. Polymer simulation shows that both intra-chromosomal self-associating interactions and extra-chromosomal interactions are necessary to establish the observed organization. Our results illustrate a multi-faceted picture and physical principles of chromatin organization.

The three-dimensional architecture of the genome affects genomic functions. Here, the authors developed Multiplexed Imaging of Nucleome Architectures to measure multiscale chromatin folding, RNA profiles, and associations of numerous genomic regions with nuclear lamina and nucleoli in the same, single cells in heterogeneous tissue.

Mechanistic modeling of chromatin folding to understand function

Experimental approaches have been applied to address questions in understanding three-dimensional chromatin organization and function. As datasets increase in size and complexity, it becomes a challenge to reach a mechanistic interpretation of experimental results. Polymer simulations and mechanistic modeling have been applied to explain experimental observations and their links to different aspects of genome function. Here we provide a guide for biologists, explaining different simulation approaches and the contexts in which they have been used.

A Lamin-Associated Chromatin Model for Chromosome Organization

We propose a simple model for chromatin organization based on the interaction of the chromatin fibers with lamin proteins along the nuclear membrane. Lamin proteins are known to be a major factor that influences chromatin organization and hence gene expression in the cells. We provide a quantitative understanding of lamin-associated chromatin organization in a crowded macromolecular environment by systematically varying the heteropolymer segment distribution and the strength of the lamin-chromatin attractive interaction. Our minimal polymer model reproduces the formation of lamin-associated-domains and provides an in silico tool for quantifying domain length distributions for different distributions of heteropolymer segments. We show that a Gaussian distribution of heteropolymer segments, coupled with strong lamin-chromatin interactions, can qualitatively reproduce observed length distributions of lamin-associated-domains. Further, lamin-mediated interaction can enhance the formation of chromosome territories as well as the organization of chromatin into tightly packed heterochromatin and the loosely packed gene-rich euchromatin regions.

Lamina-associated domains: peripheral matters and internal affairs

At the nuclear periphery, associations of chromatin with the nuclear lamina through lamina-associated domains (LADs) aid functional organization of the genome. We review the organization of LADs and provide evidence of LAD heterogeneity from cell ensemble and single-cell data. LADs are typically repressive environments in the genome; nonetheless, we discuss findings of lamin interactions with regulatory elements of active genes, and the role lamins may play in genome regulation. We address the relationship between LADs and other genome organizers, and the involvement of LADs in laminopathies. The current data lay the basis for future studies on the significance of lamin-chromatin interactions in health and disease.

4D Genome Rewiring during Oncogene-Induced and Replicative Senescence

To understand the role of the extensive senescence-associated 3D genome reorganization, we generated genome-wide chromatin interaction maps, epigenome, replication-timing, whole-genome bisulfite sequencing, and gene expression profiles from cells entering replicative senescence (RS) or upon oncogene-induced senescence (OIS). We identify senescence-associated heterochromatin domains (SAHDs). Differential intra- versus inter-SAHD interactions lead to the formation of senescence-associated heterochromatin foci (SAHFs) in OIS but not in RS. This OIS-specific configuration brings active genes located in genomic regions adjacent to SAHDs in close spatial proximity and favors their expression. We also identify DNMT1 as a factor that induces SAHFs by promoting HMGA2 expression. Upon DNMT1 depletion, OIS cells transition to a 3D genome conformation akin to that of cells in replicative senescence. These data show how multi-omics and imaging can identify critical features of RS and OIS and discover determinants of acute senescence and SAHF formation.

•

Deep multi-omics characterization of replicative and oncogene-induced senescence

•

Senescence-associated heterochromatin domains (SAHDs) form SAHFs via 3D changes

•

DNMT1 is required for SAHF formation via regulation of HMGA2 expression

•

SAHF formation leads to expression of SAHF-adjacent genes via 3D chromatin contacts

Integrating transposable elements in the 3D genome

Chromosome organisation is increasingly recognised as an essential component of genome regulation, cell fate and cell health. Within the realm of transposable elements (TEs) however, the spatial information of how genomes are folded is still only rarely integrated in experimental studies or accounted for in modelling. Whilst polymer physics is recognised as an important tool to understand the mechanisms of genome folding, in this commentary we discuss its potential applicability to aspects of TE biology. Based on recent works on the relationship between genome organisation and TE integration, we argue that existing polymer models may be extended to create a predictive framework for the study of TE integration patterns. We suggest that these models may offer orthogonal and generic insights into the integration profiles (or "topography") of TEs across organisms. In addition, we provide simple polymer physics arguments and preliminary molecular dynamics simulations of TEs inserting into heterogeneously flexible polymers. By considering this simple model, we show how polymer folding and local flexibility may generically affect TE integration patterns. The preliminary discussion reported in this commentary is aimed to lay the foundations for a large-scale analysis of TE integration dynamics and topography as a function of the three-dimensional host genome.

Polymer models for the mechanisms of chromatin 3D folding: review and perspective

In this perspective paper, classical physical models for mammalian interphase chromatin folding are reviewed.

Nonequilibrium Theory of Epigenomic Microphase Separation in the Cell Nucleus

Understanding the spatial organization of the genome in the cell nucleus is one of the current grand challenges in biophysics. Certain biochemical-or epigenetic-marks that are deposited along the genome are thought to play an important, yet poorly understood, role in determining genome organization and cell identity. The physical principles underlying the interplay between epigenetic dynamics and genome folding remain elusive. Here we propose and study a theory that assumes a coupling between epigenetic mark and genome densities, and which can be applied at the scale of the whole nucleus. We show that equilibrium models are not compatible with experiments and a qualitative agreement is recovered by accounting for nonequilibrium processes that can stabilize microphase separated epigenomic domains. We finally discuss the potential biophysical origin of these terms.