Single-molecule tracking of Nanog and Oct4 in living mouse embryonic stem cells uncovers a feedback mechanism of pluripotency maintenance.

Nanog and Oct4 are core transcription factors that form part of a gene regulatory network to regulate hundreds of target genes for pluripotency maintenance in mouse embryonic stem cells (ESCs). To understand their function in the pluripotency maintenance, we visualised and quantified the dynamics of single molecules of Nanog and Oct4 in a mouse ESCs during pluripotency loss. Interestingly, Nanog interacted longer with its target loci upon reduced expression or at the onset of differentiation, suggesting a feedback mechanism to maintain the pluripotent state. The expression level and interaction time of Nanog and Oct4 correlate with their fluctuation and interaction frequency, respectively, which in turn depend on the ESC differentiation status. The DNA viscoelasticity near the Oct4 target locus remained flexible during differentiation, supporting its role either in chromatin opening or a preferred binding to uncondensed chromatin regions. Based on these results, we propose a new negative feedback mechanism for pluripotency maintenance via the DNA condensation state-dependent interplay of Nanog and Oct4.

Centromere/kinetochore is assembled through CENP-C oligomerization.

Kinetochore is an essential protein complex required for accurate chromosome segregation. The constitutive centromere-associated network (CCAN), a subcomplex of the kinetochore, associates with centromeric chromatin and provides a platform for the kinetochore assembly. The CCAN protein CENP-C is thought to be a central hub for the centromere/kinetochore organization. However, the role of CENP-C in CCAN assembly needs to be elucidated. Here, we demonstrate that both the CCAN-binding domain and the C-terminal region that includes the Cupin domain of CENP-C are necessary and sufficient for chicken CENP-C function. Structural and biochemical analyses reveal self-oligomerization of the Cupin domains of chicken and human CENP-C. We find that the CENP-C Cupin domain oligomerization is vital for CENP-C function, centromeric localization of CCAN, and centromeric chromatin organization. These results suggest that CENP-C facilitates the centromere/kinetochore assembly through its oligomerization.

Condensed but liquid-like domain organization of active chromatin regions in living human cells.

In eukaryotes, higher-order chromatin organization is spatiotemporally regulated as domains, for various cellular functions. However, their physical nature in living cells remains unclear (e.g., condensed domains or extended fiber loops; liquid-like or solid-like). Using novel approaches combining genomics, single-nucleosome imaging, and computational modeling, we investigated the physical organization and behavior of early DNA replicated regions in human cells, which correspond to Hi-C contact domains with active chromatin marks. Motion correlation analysis of two neighbor nucleosomes shows that nucleosomes form physically condensed domains with ~150-nm diameters, even in active chromatin regions. The mean-square displacement analysis between two neighbor nucleosomes demonstrates that nucleosomes behave like a liquid in the condensed domain on the ~150 nm/~0.5 s spatiotemporal scale, which facilitates chromatin accessibility. Beyond the micrometers/minutes scale, chromatin seems solid-like, which may contribute to maintaining genome integrity. Our study reveals the viscoelastic principle of the chromatin polymer; chromatin is locally dynamic and reactive but globally stable.

PHi-C2: interpreting Hi-C data as the dynamic 3D genome state.

SUMMARY: High-throughput chromosome conformation capture (Hi-C) is a widely used assay for studying the three-dimensional (3D) genome organization across the whole genome. Here, we present PHi-C2, a Python package supported by mathematical and biophysical polymer modeling that converts input Hi-C matrix data into the polymer model's dynamics, structural conformations and rheological features. The updated optimization algorithm for regenerating a highly similar Hi-C matrix provides a fast and accurate optimal solution compared to the previous version by eliminating the factors underlying the inefficiency of the optimization algorithm in the iterative optimization process. In addition, we have enabled a Google Colab workflow to run the algorithm, wherein users can easily change the parameters and check the results in the notebook. Overall, PHi-C2 represents a valuable tool for mining the dynamic 3D genome state embedded in Hi-C data. AVAILABILITY AND IMPLEMENTATION: PHi-C2 as the phic Python package is freely available under the GPL license and can be installed from the Python package index. The source code is available from GitHub at https://github.com/soyashinkai/PHi-C2. Moreover, users do not have to prepare a Python environment because PHi-C2 can run on Google Colab (https://bit.ly/3rlptGI). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Single-nucleosome imaging reveals steady-state motion of interphase chromatin in living human cells.

Dynamic chromatin behavior plays a critical role in various genome functions. However, it remains unclear how chromatin behavior changes during interphase, where the nucleus enlarges and genomic DNA doubles. While the previously reported chromatin movements varied during interphase when measured using a minute or longer time scale, we unveil that local chromatin motion captured by single-nucleosome imaging/tracking on a second time scale remained steady throughout G1, S, and G2 phases in live human cells. This motion mode appeared to change beyond this time scale. A defined genomic region also behaved similarly. Combined with Brownian dynamics modeling, our results suggest that this steady-state chromatin motion was mainly driven by thermal fluctuations. Steady-state motion temporarily increased following a DNA damage response. Our findings support the viscoelastic properties of chromatin. We propose that the observed steady-state chromatin motion allows cells to conduct housekeeping functions, such as transcription and DNA replication, under similar environments during interphase.

Toward understanding the dynamic state of 3D genome.

The three-dimensional (3D) genome organization and its role in biological activities have been investigated for over a decade in the field of cell biology. Recent studies using live-imaging and polymer simulation have suggested that the higher-order chromatin structures are dynamic; the stochastic fluctuations of nucleosomes and genomic loci cannot be captured by bulk-based chromosome conformation capture techniques (Hi-C). In this review, we focus on the physical nature of the 3D genome architecture. We first describe how to decode bulk Hi-C data with polymer modeling. We then introduce our recently developed PHi-C method, a computational tool for modeling the fluctuations of the 3D genome organization in the presence of stochastic thermal noise. We also present another new method that analyzes the dynamic rheology property (represented as microrheology spectra) as a measure of the flexibility and rigidity of genomic regions over time. By applying these methods to real Hi-C data, we highlighted a temporal hierarchy embedded in the 3D genome organization; chromatin interaction boundaries are more rigid than the boundary interior, while functional domains emerge as dynamic fluctuations within a particular time interval. Our methods may bridge the gap between live-cell imaging and Hi-C data and elucidate the nature of the dynamic 3D genome organization.

Microrheology for Hi-C Data Reveals the Spectrum of the Dynamic 3D Genome Organization.

The one-dimensional information of genomic DNA is hierarchically packed inside the eukaryotic cell nucleus and organized in a three-dimensional (3D) space. Genome-wide chromosome conformation capture (Hi-C) methods have uncovered the 3D genome organization and revealed multiscale chromatin domains of compartments and topologically associating domains (TADs). Moreover, single-nucleosome live-cell imaging experiments have revealed the dynamic organization of chromatin domains caused by stochastic thermal fluctuations. However, the mechanism underlying the dynamic regulation of such hierarchical and structural chromatin units within the microscale thermal medium remains unclear. Microrheology is a way to measure dynamic viscoelastic properties coupling between thermal microenvironment and mechanical response. Here, we propose a new, to our knowledge, microrheology for Hi-C data to analyze the dynamic compliance property as a measure of rigidness and flexibility of genomic regions along with the time evolution. Our method allows the conversion of an Hi-C matrix into the spectrum of the dynamic rheological property along the genomic coordinate of a single chromosome. To demonstrate the power of the technique, we analyzed Hi-C data during the neural differentiation of mouse embryonic stem cells. We found that TAD boundaries behave as more rigid nodes than the intra-TAD regions. The spectrum clearly shows the dynamic viscoelasticity of chromatin domain formation at different timescales. Furthermore, we characterized the appearance of synchronous and liquid-like intercompartment interactions in differentiated cells. Together, our microrheology data derived from Hi-C data provide physical insights into the dynamics of the 3D genome organization.

PHi-C: deciphering Hi-C data into polymer dynamics.

Genomes are spatiotemporally organized within the cell nucleus. Genome-wide chromosome conformation capture (Hi-C) technologies have uncovered the 3D genome organization. Furthermore, live-cell imaging experiments have revealed that genomes are functional in 4D. Although computational modeling methods can convert 2D Hi-C data into population-averaged static 3D genome models, exploring 4D genome nature based on 2D Hi-C data remains lacking. Here, we describe a 4D simulation method, PHi-C (polymer dynamics deciphered from Hi-C data), that depicts 4D genome features from 2D Hi-C data by polymer modeling. PHi-C allows users to interpret 2D Hi-C data as physical interaction parameters within single chromosomes. The physical interaction parameters can then be used in the simulations and analyses to demonstrate dynamic characteristics of genomic loci and chromosomes as observed in live-cell imaging experiments. PHi-C is available at https://github.com/soyashinkai/PHi-C.

Spindle pole body movement is affected by glucose and ammonium chloride in fission yeast.

The complexity of chromatin dynamics is orchestrated by several active processes. In fission yeast, the centromeres are clustered around the spindle pole body (SPB) and oscillate in a microtubule- and adenosine triphosphate (ATP)-dependent manner. However, whether and how SPB oscillation are affected by different environmental conditions remain poorly understood. In this study, we quantitated movements of the SPB component, which colocalizes with the centromere in fission yeast. We found that SPB movement was significantly reduced at low glucose concentrations. Movement of the SPB was also affected by the presence of ammonium chloride. Power spectral analysis revealed that periodic movement of the SPB is disrupted by low glucose concentrations. Measurement of ATP levels in living cells by quantitative single-cell imaging suggests that ATP levels are not the only determinant of SPB movement. Our results provide novel insight into how SPB movement is regulated by cellular energy status and additional factors such as the medium nutritional composition.

Bridging the dynamics and organization of chromatin domains by mathematical modeling.

The genome is 3-dimensionally organized in the cell, and the mammalian genome DNA is partitioned into submegabase-sized chromatin domains. Genome functions are regulated within and across the domains according to their organization, whereas the chromatin itself is highly dynamic. However, the details of such dynamic organization of chromatin domains in living cells remain unclear. To unify chromatin dynamics and organization, we recently demonstrated that structural information of chromatin domains in living human cells can be extracted from analyses of the subdiffusive nucleosome movement using mathematical modeling. Our mathematical analysis suggested that as the chromatin domain becomes smaller and more compact, nucleosome movement becomes increasingly restricted. Here, we show the implication of these results for bridging the gap between chromatin dynamics and organization, and provide physical insight into chromatin domains as efficient units to conduct genome functions in the thermal noisy environment of the cell.

Dynamic Nucleosome Movement Provides Structural Information of Topological Chromatin Domains in Living Human Cells.

The mammalian genome is organized into submegabase-sized chromatin domains (CDs) including topologically associating domains, which have been identified using chromosome conformation capture-based methods. Single-nucleosome imaging in living mammalian cells has revealed subdiffusively dynamic nucleosome movement. It is unclear how single nucleosomes within CDs fluctuate and how the CD structure reflects the nucleosome movement. Here, we present a polymer model wherein CDs are characterized by fractal dimensions and the nucleosome fibers fluctuate in a viscoelastic medium with memory. We analytically show that the mean-squared displacement (MSD) of nucleosome fluctuations within CDs is subdiffusive. The diffusion coefficient and the subdiffusive exponent depend on the structural information of CDs. This analytical result enabled us to extract information from the single-nucleosome imaging data for HeLa cells. Our observation that the MSD is lower at the nuclear periphery region than the interior region indicates that CDs in the heterochromatin-rich nuclear periphery region are more compact than those in the euchromatin-rich interior region with respect to the fractal dimensions as well as the size. Finally, we evaluated that the average size of CDs is in the range of 100-500 nm and that the relaxation time of nucleosome movement within CDs is a few seconds. Our results provide physical and dynamic insights into the genome architecture in living cells.

Generalized Lyapunov exponent as a unified characterization of dynamical instabilities.

The Lyapunov exponent characterizes an exponential growth rate of the difference of nearby orbits. A positive Lyapunov exponent (exponential dynamical instability) is a manifestation of chaos. Here, we propose the Lyapunov pair, which is based on the generalized Lyapunov exponent, as a unified characterization of nonexponential and exponential dynamical instabilities in one-dimensional maps. Chaos is classified into three different types, i.e., superexponential, exponential, and subexponential chaos. Using one-dimensional maps, we demonstrate superexponential and subexponential chaos and quantify the dynamical instabilities by the Lyapunov pair. In subexponential chaos, we show superweak chaos, which means that the growth of the difference of nearby orbits is slower than a stretched exponential growth. The scaling of the growth is analytically studied by a recently developed theory of a continuous accumulation process, which is related to infinite ergodic theory.