Skip to main content

Thank you for visiting You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Chromatin network retards nucleoli coalescence


Nuclear bodies are membraneless condensates that may form via liquid-liquid phase separation. The viscoelastic chromatin network could impact their stability and may hold the key for understanding experimental observations that defy predictions of classical theories. However, quantitative studies on the role of the chromatin network in phase separation have remained challenging. Using a diploid human genome model parameterized with chromosome conformation capture (Hi-C) data, we study the thermodynamics and kinetics of nucleoli formation. Dynamical simulations predict the formation of multiple droplets for nucleolar particles that experience specific interactions with nucleolus-associated domains (NADs). Coarsening dynamics, surface tension, and coalescence kinetics of the simulated droplets are all in quantitative agreement with experimental measurements for nucleoli. Free energy calculations further support that a two-droplet state, often observed for nucleoli in somatic cells, is metastable and separated from the single-droplet state with an entropic barrier. Our study suggests that nucleoli-chromatin interactions facilitate droplets’ nucleation but hinder their coarsening due to the coupled motion between droplets and the chromatin network: as droplets coalesce, the chromatin network becomes increasingly constrained. Therefore, the chromatin network supports a nucleation and arrest mechanism to stabilize the multi-droplet state for nucleoli and possibly for other nuclear bodies.


Nuclear bodies are pervasive in eukaryotic cells and play a diverse set of functions1, including RNA metabolism, transcriptional regulation2, genome organization3, etc. They are membraneless structures that mainly consist of protein and RNA molecules4. Their lack of a lipid-rich barrier allows rapid exchange of components with the nucleoplasm in responses to environmental cues and stress signaling5. Nuclear bodies also effectively increase the local concentration of enzymes involved in particular functions to facilitate more efficient cellular reactions6.

Increasing evidence supports that nuclear bodies function as biomolecular condensates formed via liquid-liquid phase separation (LLPS)4,7,8,9. They exhibit round morphologies and dynamic fluid properties10,11. Two nuclear bodies can fuse into larger condensates following growth kinetics with similar scaling behavior as that observed for simple liquids12,13,14. In addition, their assembly was shown to be concentration-dependent, and the coarsening and growth dynamics can be quantitatively modeled with classical theories of phase separation15. At the molecular level, detailed mechanistic models for LLPS are beginning to emerge as well. In particular, low complexity domains and intrinsically disordered regions are enriched in many of the proteins associated with nuclear bodies16,17,18,19,20. These features enable non-specific, multivalent interactions that drive the formation of dynamical condensates.

However, several observations of nuclear bodies appear to defy predictions from classical nucleation and phase separation theories. In particular, these theories predict the thermodynamic equilibrium to consist of a single condensate that minimizes the surface energy21,22. On the other hand, multiple nucleoli (~2–5) can stably coexist in the same nucleus11,12,23,24,25,26, as can paraspeckles27 and nuclear speckles28. The exact number of nuclear bodies is sensitive to various factors, including cell volume29 and nuclear lamina composition30,31. It has been proposed that non-equilibrium activities can dynamically alter protein-protein interactions to stabilize the multi-droplet state15,32,33,34,35. In addition, the chromatin network may suppress droplet coarsening through mechanical frustration as well36,37,38.

Nuclear bodies and chromatin are also known to form attractive interactions27, further complicating phase separation inside the nucleus beyond mechanical stress. For example, the upstream binding factor (UBF), which is a DNA binding protein and a key component of nucleoli, is known to recognize ribosomal DNA (rDNA) repeats to seed the rapid formation of nucleoli after cell division39. Correspondingly, rDNA and other chromosome segments, which are collectively noted as nucleolus-association domains (NADs)9,40,41, can be seen inside and adjacent to nucleoli42. Paraspeckles27 and speckles28,43 have been found in spatial proximity with chromatin as well. In addition to proteins, nuclear bodies can harbor non-coding RNA to contact chromatin either by recruiting intermediate protein molecules or by forming RNA-DNA duplex or triplex44. Since chromatin forms a viscoelastic network spanning the nucleus13,45,46,47,48,49, its interactions with nuclear bodies could impact the thermodynamics and kinetics of phase separation.

We carry out molecular dynamics simulations to investigate nucleoli formation with a computational model that explicitly considers nucleolar particles, the chromatin network, and the interactions between the two. We represent the chromatin network using a diploid human genome model that provides explicit polymer configurations for individual chromosomes. Interactions within and among chromosomes are optimized based on chromosome conformation capture (Hi-C) experiments to ensure in vivo relevance. The simulated dynamical process of droplet growth, coarsening, and coalescence are in quantitative agreement with experimental measurements. Importantly, our simulations predict the formation of multiple droplets, much like the coexistence of several nucleoli seen in the nucleus. We show that a two-droplet state is metastable and separated from the single-droplet state with an entropic barrier. The barrier arises from the chromatin network, which becomes more constrained upon droplet coalescence. Nucleolar particle-chromatin interactions link the motion between the chromatin network and the droplets, and stronger interactions are shown to produce more droplets. Our study provides insight into the critical role of the chromatin network on the formation of nucleoli and nuclear bodies in general.


Phase separation with chromatin network leads to multiple droplets

Leveraging a recently introduced computational model for the diploid human genome50,51, we studied the impact of the chromatin network on nucleoli formation. We modeled the genome at the one megabase (Mb) resolution as a collection of 46 chromosomes inside spherical confinement (Fig. 1). Each chromosome is represented as a string of beads, which can be assigned as one of three types, A, B, or C. A and B correspond to the two compartment types that contribute to the checkerboard patterns typically seen in Hi-C contact maps52,53, and C marks centromeric regions. Interactions among the beads were optimized to reproduce various average contact probabilities determined from Hi-C experiments for GM12878 cells using the maximum entropy optimization algorithm54,55 (see Supplementary Material). Because of the non-equilibrium nature of the system, these experimentally derived interactions and temperature represent effective approximations to the steady-state distribution of genome organization (see “Methods”). While the model does not explicitly include histone modifications, transcription factors, or molecular motors, reproducing the contact probabilities between genomic segments measured in situ effectively allows it to account for the contribution of these factors to genome organization.

Fig. 1: Overview of the computational model used for studying phase separation inside the nucleus.
figure 1

The diploid human genome model represents each one of the 46 chromosomes as a string of beads confined in the nuclear envelope. Each bead is identified as compartment type A, B, or centromeric region C. Nucleolar particles share favorable interactions to promote phase separation and bind specifically with nucleolus-associated domains (NADs).

We introduced additional coarse-grained particles to model molecules that make up the nucleoli. The size and number of these particles were chosen based on the concentration of nucleolar proteins and the volume fraction of nucleoli (see “Methods”). Nucleolar particles share favorable interactions with each other and with nucleolus-associated domains (NADs), which are chromatin regions strongly associated with nucleoli41. A non-specific interaction term was also introduced between nucleolar particles and non-NAD chromatin regions. While our results are relatively robust with respect to the strength of these interactions (Supplementary Fig. 1), we found that, with the chosen values, the surface tension of simulated droplets compares favorably to experimental values for nucleoli (see Supplementary Material).

Starting from an equilibrated genome structure and randomly distributed nucleolar particles (see “Methods”), we carried out a total of twelve independent molecular dynamics simulations. These simulations lasted for 20 million steps, much longer than the relaxation timescale of chromosome conformations (Supplementary Fig. 2). In all but one case, the nucleolar particles aggregated into multiple droplets that persisted to the end of the simulations (Fig. 2b). This result contrasts with simulations performed without the chromatin network, where nucleolar particles always condense into a single droplet (Supplementary Fig. 3). Notably, the emergence of the multi-droplet state is insensitive to the configuration used to initialize the simulations (Supplementary Fig. 4), the interactions between chromosomes (Supplementary Figs. 5,6), and the resolution of the genome model (Supplementary Fig. 7). We found that nucleolar particles forming the droplets undergo dynamical exchange with the surrounding nucleoplasm while maintaining the droplet size (~1 μm) on timescales of several tens of minutes (Supplementary Fig. 8). Dynamical exchange of materials has been observed in fluorescence recovery after photobleaching (FRAP) experiments56,57 and directly supports the liquid-like property of the droplets. The droplets are preferentially localized at the interior of the nucleus (Supplementary Fig. 9A). Because of their close association with these droplets, NADs are closer to the nuclear center than other heterochromatin regions as well58. However, not all NADs bind to the droplets, and a significant fraction of them localize towards the nuclear envelope (Supplementary Fig. 9B). Two classes of NADs that vary in nuclear localization have indeed been observed in prior studies59.

Fig. 2: Multiple droplets form in dynamical simulations of the nucleus model.
figure 2

a Representative initial (left) and final (right) configurations obtained from simulations, with nucleolar particles in yellow, NADs in black, and the rest of the genome in gray. b Probability distribution of the number of droplets observed at the end of simulation trajectories. c Representative configuration of the genome that illustrates the formation of chromosome territories. d Radial distributions of the different chromatin types that support their phase separation and preferential nuclear localization. An example genome configuration is shown as the inset, with the three types colored in red (A), blue (B), and orange (C), respectively. e Simulated radial chromosome positions correlate strongly with experimental values63. Error bars correspond to the standard deviation of the 12 mean values estimated using individual simulation trajectories. Homologous chromosomes were averaged together. RN is the radius of the nucleus used in polymer simulations.

Analyzing the simulated genome structures, we found that the model with nucleolar particles reproduces global features of the genome organization, including the formation of chromosome territories60, the compartmentalization of heterochromatin/euchromatin61, and the clustering of centromeric regions62, as shown in Fig. 2c, d. The simulated chromosome radial positions agree well with experimental values63, and the Pearson’s correlation coefficient between the two is 0.89 (Fig. 2e). Because of its coarse resolution, the model will inevitably miss certain features of genome organizations, including the formation of chromatin loops and topologically associating domains53. While these structural motifs at fine scales are crucial for an accurate representation of the genome organization, they are less significant for exploring the mechanisms of nuclear body formation, at least at a qualitative level (see Supplementary Fig. 7).

An entropic barrier hinders droplet coalescence

The emergence of a long-lived state with multiple droplets, while in contrast with predictions of the classical nucleation theory21, is consistent with the coexistence of multiple nucleoli in the nucleus11,13,24,25,26. To better understand the stability of the multi-droplet state, we computed the free energy profile as a function of the radius of gyration (Rg) for a two-droplet system (see “Methods”). Rg is defined as \(\sqrt{\frac{1}{N}\mathop{\sum }\nolimits_{i = 1}^{N}{\left|{{{{{{{{\bf{r}}}}}}}}}_{i}-{{{{{{{{\bf{r}}}}}}}}}_{{{\mbox{com}}}}\right|}^{2}}\), where ri is the coordinate of the ith nucleolar particle and rcom corresponds to the center of mass. The summation includes all N nucleolar particles in either one of the two droplets. As the size of individual droplets remains stable, changes in Rg will be mainly driven by variations in the distance of the two droplets. However, different from a simple center-of-mass distance, which becomes ill-defined if the lists of nucleolar particles participating in droplet formation are not updated on the fly, Rg is relatively invariant with respect to the flux of particles between the two droplets (Supplementary Fig. 10). Simulations of the full system with both chromosomes and nucleolar particles were used to compute the free energy profile. Therefore, the impact of the chromatin network was accounted for implicitly even though it was not included in the definition of Rg. Umbrella sampling and temperature replica exchange were used to enhance conformational exploration.

As shown in Fig. 3a, the free energy profile exhibits two basins. While the basin at Rg ≈ 1.5 μm corresponds to the two droplet state observed before, an additional minimum with all nucleolar particles participating in a single droplet appears at smaller values (≈0.75 μm). The two basins are separated from each other with a transition state at Rg ≈ 1.13 μm. Representative configurations at the transition state show a dumbbell shape with the establishment of a thin bridge between the two droplets. Consistent with predictions of the classical nucleation theory, the one droplet state remains as the global minimum. However, the merging of the droplets is kinetically constrained due to the presence of a barrier that is ~7 kBT in height. The barrier height is much larger than the error bars (0.5 kBT) estimated via block averaging (see “Methods”), supporting its statistical significance.

Fig. 3: An entropic barrier hinders droplet coalescence and drives the metastability of a two-droplet state.
figure 3

a Free energy profile as a function of the radius of gyration (Rg) that effectively measures the distance between two droplets. Energetic (black) and entropic (red) contributions to the free energy profile before (c) and after (b) the coalescence barrier. Error bars were calculated as standard deviations of the mean using block averaging by dividing the simulation trajectories into five blocks of equal length.

To reveal the nature of the barrier, we decomposed the free energy into entropic and energetic contributions. Using the free energy profiles at different temperatures (Supplementary Fig. 11), we computed the entropy change along the collective variable with a finite difference method64,65. As shown in Fig. 3b, c, contributions from the entropy (–TΔS) continue to increase as Rg decreases, and the droplets coalesce. While restricting the motion of the two droplets to smaller distances is naturally unfavorable, the entropic penalty is intensified here due to the increasingly restricted motion of chromosomes as well. As mentioned before, nucleolar particles form extensive contacts with NADs via specific interactions, and such contacts enforce correlative motions between droplets and chromosomes. The potential energy, on the other hand, favors droplet merging and decreases continuously along the collective variable due to the increase in nucleolar particle-nucleolar particle and nucleolar particle-chromatin contacts. Therefore, the transition state arises due to the presence of the chromatin network, and dissolving the polymeric topology of chromosomes indeed removes the barrier (Supplementary Fig. 3C).

While Rg suffices for monitoring the progression of droplet coalescence, whether it serves as a “good” reaction coordinate or not requires further investigation. In particular, a good reaction coordinate should provide insight into the bottleneck that limits the reaction. Furthermore, trajectories initialized from the identified transition state will have an equal probability of committing to the reactant or product, or the so-called committor probability adopts a value of 0.566,67. Otherwise, the transition state may bear little relevance to the reaction, and mechanisms derived from it can be misleading. To evaluate the significance of the transition state for droplet coalescence, we carried out additional simulations starting from random configurations with Rg values at 1.13 μm. For each configuration, we initialized ten independent 200,000-step-long simulations with randomized velocities. We then counted the number of simulations that end up with the two-droplet state versus the single-droplet state. As shown in Fig. 4, among all these simulations, 56% of them led to the single-droplet state while the rest 44% ended up in the two-droplet state. These results strongly support the usefulness of Rg for studying coalescence and the relevance of the identified transition state for mechanistic interpretation. Since the chromatin network was not included for defining the reaction coordinate and can vary significantly at a given value for Rg (Supplementary Fig. 12), it may play a secondary or passive role in mediating coalescence.

Fig. 4: Dynamical characterization of droplet coalescence.
figure 4

a Trajectories initialized from transition state configurations share approximately equal probability of committing to the two and single droplet state. The number of trajectories that ended up in the single droplet state is shown as red blocks in the bottom panel. Example configurations of the single-droplet (left), transition (middle), and two-droplet (right) state are shown on top. b Time evolution of the radius of the neck region between two droplets averaged over twelve independent simulations.

We further isolated trajectories initialized from the transition state that led to the single-droplet state and computed the evolution of the neck radius as a function of time. As shown in Supplementary Fig. 13, the neck region was identified as the minimum of the density profile of nucleolar particles along the principal axis with the largest variance. As shown in Fig. 4b, by plotting the normalized neck radius (2R(t)/R0) with respect to the time, we obtained a power-law relationship with exponent 0.51. R0 is the average radius of the droplets before fusion. This exponent agrees with the experimental value determined for nucleoli13 and suggests that droplet coalescence proceeds in the low Reynolds number regime dominated by viscous effects from the outer fluid, i.e., the nucleoplasm68,69.

Chromatin network gives rise to slow coarsening dynamics

The thermodynamic analysis presented in Fig. 3 suggests that the chromatin network acts much like entropic springs to hinder the coalescence of droplets70. As the droplets move close to each other, they pull on the network and restrict the accessible polymer configurations. Restoring forces from the network to maximize configurational entropy counter droplet merging and give rise to the barrier. A similar mechanism could potentially impact the coarsening dynamics and the pathway leading to the formation of multiple droplets. To better understand the role of nucleolar particle-chromatin interactions in the overall phase separation process, we analyzed the dynamical trajectories at the onset of cluster formation.

We first monitored the time evolution of the number of clusters formed along the dynamical trajectories shown in Fig. 2. The clusters were identified as high-density regions of nucleolar particles across the entire nucleus using the DBSCAN (Density-Based Spatial Clustering of Applications with Noise) algorithm71 (see Supplementary Materials). A typical trajectory is shown in Fig. 5a, and starts with zero clusters due to the random distribution of nucleolar particles in the initial configuration. The sudden appearance of nine clusters at time ~ 14 min suggests that nucleation can occur at multiple sites almost simultaneously. As time proceeds, the droplets began to merge or evaporate, and the trajectory eventually stabilizes to the two-droplet state.

Fig. 5: Coarsening dynamics and pathways for phase separation with the chromatin network.
figure 5

a Time evolution of the number of clusters observed along a dynamical trajectory. b Detailed characterization of the two switching events that follow the Brownian motion-induced coalescence (b1) and the diffusion-limited Ostwald ripening (b2) path. The average pair-wise distance for each cluster remains relatively constant along the first path but increases significantly along the second one due to particle evaporation. c Power-law scaling of the average cluster size as a function of time.

We further recorded the average size of all clusters measured by their mean Rg. The increase of Rg as a function of the time follows a power-law scaling with Rgtβ. The exponent β = 0.1 differs from the phase separation process without the chromatin network. As shown in Supplementary Fig. 14, for simulations performed without the network, the initial cluster number is smaller and the exponent for cluster size growth is close to the theoretical value of 1/372,73.

To provide insight into the appearance of the abnormal exponent β, we characterized the switching events that led to the decrease in cluster numbers. Specifically, we tracked the average pair-wise distance \(\overline{{r}_{ij}(t)}\) of nucleolar particles within each cluster. In some cases, \(\overline{{r}_{ij}(t)}\) for all clusters remain relatively constant and the switching occurs through the Brownian motion-induced coalescence (BMC) pathway (Fig. 5b1). In other cases, we observe a significant increase in \(\overline{{r}_{ij}(t)}\) before the switching occurs (Fig. 5b2). This increase indicates an “evaporation” of the cluster following the diffusion-limited Ostwald ripening (DOR) path. We then identified all the switching events (42 in total) from the dynamical simulations and found that ~76% of them follow the BMC path. The rest of the switching events proceed via the DOR path and often involve smaller clusters to reduce the penalty associated with cluster disassembly (Supplementary Fig. 15).

The dominance of the BMC pathway explains the dramatic slowdown of the coarsening dynamics. In particular, the scaling exponent of 1/3 was predicted based on a normal diffusion model in which distances between droplets scale linearly in time, i.e., x2(t) Dt74. However, as shown in Supplementary Fig. 16, most of the clusters exhibit sub-diffusion and x2(t) Dt1/2. Assuming that the average size of droplets r(t) is proportional to their mean distance, we arrive at r(t) t1/6 for the observed abnormal diffusion. The exponent now is closer to the value shown in Fig. 5c. The sub-diffusive motion arises from both the elastic stress produced by the viscoelastic chromatin network and the physical tethering of droplets to the chromosomes. In addition to the sub-diffusive motion, the chromatin network could further reduce the exponent and slow down the Brownian diffusion dominated coarsening dynamics by hindering droplet coalescence through entropic barriers similar to that shown in Fig. 3.

We note that the abnormal diffusion and slower coarsening kinetics have been directly observed by Lee et al. as well when monitoring the coarsening dynamics of model condensates based on intrinsically disordered protein regions in the nucleus36. In particular, they revealed a coarsening exponent of ~0.12, which is close to the value shown in Fig. 5c. The scaling exponent for nucleolar coarsening in vivo is also in good agreement with the simulated value when considering short time kinetics before 5 min15.

Nucleoli-chromatin interactions promote cluster nucleation

Our results indicate that, while nucleolar particle-chromatin interactions increase the overall stability of the single-droplet state, they retard the coarsening kinetics by giving rise to sub-diffusion and entropic barriers. To more directly probe their impact on droplet coalescence, we recomputed the free energy profile at stronger (ϵ = 2.0 kBT) and weaker (ϵ = 1.6 kBT) nucleolar particle-chromatin interactions. As shown in Fig. 6a, while the stability of the merged state varied significantly, the transition and the two-droplet state are much less affected. These simulations again support the entropic origin of the free energy barrier.

Fig. 6: Nucleoli-chromatin interactions promote cluster nucleation.
figure 6

a The free energy profiles of droplet coalescence at various nucleolar particle-chromatin interactions. The result for ϵ = 1.8 is identical to the one presented in Fig. 3. Error bars were calculated as standard deviation of the mean using block averaging by dividing the simulation trajectories into five blocks of equal length. b Number of droplets formed at the end of 12 independent dynamical simulations at various strengths of nucleolar particle-chromatin interactions. Error bars correspond to standard deviations of results from independent trajectories. c The maximum number of clusters observed in dynamical trajectories performed at various nucleolar particle-chromatin interactions. d The fraction of cluster switching events following the Brownian motion-induced coalescence (BMC, blue) and the diffusion-limited Ostwald ripening (DOR, red) at various nucleolar particle-chromatin interactions.

We performed additional long-timescale simulations following the same protocol as those shown in Fig. 2. Starting from the same initial configurations, we varied nucleolar particle-chromatin interactions from 1.4 kBT to 2.8 kBT and carried out 12 independent 20-million-step long simulations for each interaction strength. We then computed the number of droplets formed at the end of these simulations. As shown in Fig. 6b, the number of droplets increases with the interaction strength. Stronger interactions facilitate the nucleation of nucleolar clusters on the chromatin, as can be seen from the increase of cluster numbers at the onset of phase separation (Fig. 6c). Cluster coarsening, again, primarily follow the BMC pathway, though the ratio of DOR increases as well due to the instability of the nucleated clusters with a smaller number of nucleolar particles (Fig. 6d). Since merging of these clusters along the BMC pathway is hindered by entropic barriers, more nucleation naturally leads to increased droplet numbers at the end of the simulations. These results further support the role of the chromatin network in hindering droplet coalescence.


While nuclear bodies have been shown to exhibit liquid properties, the role of the surrounding environment on their formation and stability is less known. We modeled the dynamical process of phase separation that drives nucleoli formation in the presence of the chromatin network. The diploid genome model parameterized from Hi-C data was used to describe the interactions within and among chromosomes accurately. Simulations carried out with coarse-grained nucleolar particles succeeded in producing multiple droplets with dynamical behaviors comparable to nucleoli. We found that nucleoli-chromatin interactions facilitate the nucleation of condensates and retard their coalescence, stabilizing the multiple-droplet state.

Despite our best effort, the coarsened representations for the genome and nucleoli are bound to be approximate, and certain features of the biological system would inevitably be missed. While it may struggle at quantitative predictions, the model serves as a valuable tool for mechanistic explorations. In particular, since the nucleation and arrest mechanism arises from the generic polymer nature of the chromatin network, it should be insensitive to many of the model assumptions. Indeed, as shown in Supplementary Fig. 17, the multi-droplet state was also observed using a genome model parameterized with Hi-C data generated from clinical tissue samples with colorectal cancer51. The tumor Hi-C data and genome organization differ significantly from the one presented in the main text. The stability of the multi-droplet state is also preserved in additional perturbations that directly vary the resolution of the genome model (Supplementary Fig. 7), chromosome-chromosome interactions (Supplementary Fig. 5,6), nucleolar particle size (Supplementary Fig. 18), and nucleolar particle-chromatin interactions (Supplementary Fig. 1).

Specific interactions between nucleolar particles and chromatin are crucial for the nucleation step and the emergence of the multi-droplet state (Supplementary Fig. 1). As such interactions may be present for other nuclear bodies, a similar mechanism could contribute to their formation as well. For example, several recent studies support the close contact between a subset of chromatin and speckles with distances at ~100 nm or less using techniques based on high throughput sequencing75 or imaging76. Furthermore, these speckle-associated chromatin regions are largely conserved across cell-types77, supporting the presence of non-random mechanisms for their maintenance. Close contacts with chromatin are likely stabilized by RNA molecules that are known to localize at speckle periphery78,79. Non-coding RNAs have also been found in paraspeckles and PML bodies to mediate their interactions with chromatin27.

The mechanism uncovered here differs from the suppression of droplet growth in an elastic matrix discussed in several recent studies36,37,80,81. In contrast to our focus on the viscoelastic nature of the chromatin network and its attractive interactions with phase-separating agents, these studies were concerned with cross-linked gels and repulsive interactions. In the presence of permanent cross-linking, droplet-induced swelling of the polymer network can give rise to compressive stresses that arrest the coarsening dynamics completely. Therefore, the thermodynamic equilibrium becomes the monodisperse state with many droplets whose size is largely determined by the network mesh. Because of their favorable interactions in our model, there is substantial wetting of the chromatin network by nucleolar particles. The droplets incur minimal stress on the network (Supplementary Fig. 20), and the multi-droplet state remains metastable (Fig. 3). We note that the two mechanisms are not necessarily mutually exclusive and could contribute to the formation of different condensates inside the nucleus since the dynamical properties of chromatin are known to be timescale dependent82.

Finally, it’s worth mentioning that the nucleus is a non-equilibrium system, and active processes could contribute to the stability of the multi-droplet state as well34,35. For example, enzymes such as kinase could add post-translational modifications to condensate proteins and regulate their ability in engaging multivalent interactions. Explicitly modeling the active processes within the simulation framework outlined here may be necessary to account for the complementary mechanisms and provide a more comprehensive understanding of in vivo phase separation and nuclear body formation.


Details on the setup and simulations of the nucleus model

We included both the genome and nucleolar proteins to simulate phase separation inside the nucleus. Following the same setup and interactions as in ref. 50, we modeled all 46 chromosomes of the diploid human genome as beads on a string at the 1 MB resolution. In addition, a total of 500 particles were introduced to represent nucleolar proteins. See below for details on estimations for the particle size and number.

In addition to their attractive self-interactions, nucleolar particles bind with chromatin via specific and non-specific interactions. All three types of interactions were modeled with the cut and shifted Lennard-Jones potential

$${U}_{{{{{{{{\rm{LJ}}}}}}}}}({{{{{{{\bf{r}}}}}}}})=4\epsilon \left[{\left(\frac{{\sigma }_{l}}{{r}_{ij}}\right)}^{12}-{\left(\frac{{\sigma }_{l}}{{r}_{ij}}\right)}^{6}\right]+{E}_{{{{{{{{\rm{cut}}}}}}}}}$$

for rij < rc and zero otherwise, where rc = 2.0 σ. \({E}_{{{{{{{{\rm{cut}}}}}}}}}=-4\epsilon [{(\frac{{\sigma }_{l}}{{r}_{c}})}^{12}-{(\frac{{\sigma }_{l}}{{r}_{c}})}^{6}]\). ϵ = 1.8, 1.8, and 1.0 kBT, and σl = 0.5, 0.75 and 0.75 σ for nucleolar particle-nucleolar particle, nucleolar particle-NAD and nucleolar particle-non-NAD interactions. The values for σl were chosen based on the size of nucleolar particles (σp = 0.5σ) and chromatin beads (σ) with arithmetic averaging. NADs were identified using the high-resolution sequencing data generated from ref. 41. We processed the raw nucleolar-to-genomic ratios to generate signal data at the 1 Mb resolution. Only genomic loci with signals higher than a threshold value (15.0) were labeled as NADs and homologous chromosomes share identical NADs.

We used the software package LAMMPS83 to perform molecular dynamics (MD) simulations in reduced units. Constant-temperature (T = 1.0 in reduced unit) simulations were carried out via the Langevin dynamics with a damping coefficient γ = 10.0 τB and a simulation time step of dt = 0.008 τB, where τB is the Brownian time unit. Configurations were recorded at every 2000 simulation steps for analysis. The initial configuration of MD simulations was built as follows. We first obtained chromosome conformations from the end configuration of a 20-million-step long simulation of the genome-only model carried out in our previous study50. Next, 500 nucleolar particles were placed with random positions inside the spherical confinement. We further relaxed the system to avoid steric clashes by performing 400,000 step MD simulations. ϵ in Eq. (1) was set to 1.0 kBT for both specific and non-specific interactions to prevent cluster formation during the equilibration period. The last configuration of this equilibration simulation was then used to initialize the phase separation simulations. Analyses of simulation trajectories were carried out with in-house Python and MATLAB scripts.

We mapped the reduced units in MD simulations to physical units for direct comparison with experimental measurements. The length scale unit was determined by assuming a typical nucleus size with a radius of 5 μm. Since the spherical confinement mimicking the nucleus in our model adopts a radius of 19.7 σ, we have σ = 0.254 μm for the reduced length unit. We further determined the reduced time unit with the expression \({\tau }_{B}=\frac{30\pi \eta {\sigma }^{3}}{{k}_{{{{B}}}}T}\), which was obtained by matching the diffusion coefficient from simulations with that in the nucleus (see Supplementary Material). Using a value of the nucleoplasmic viscosity as η = 10−2 Pas, we estimated τB = 3.6 s.

We note that the energy unit (kBT) in our model should be viewed as an effective temperature instead of strictly the biological value (310 K). Since the nucleus is a non-equilibrium system with constant perturbations from molecular motors and chemical reactions, the ensemble of genome organizations collected over a population of cells as in Hi-C experiments is unlikely to be in thermodynamic equilibrium. However, non-equilibrium distributions can often be well approximated by renormalized Boltzmann distributions with effective potentials and temperatures84. It is these effective quantities that we inferred from Hi-C contact maps to describe the interactions among chromosomes.

Estimating the size and number of nucleolar particles

The number of nucleolar particles used in simulations was selected based on the concentration of a representative nucleolar protein NPM1 as follows. For nucleoli of 2RNu in diameter, and a protein concentration of 1 μM85, the number of nucleolar particles can be estimated as \({N}_{{{{{{{{\rm{p}}}}}}}}}=\frac{4\pi }{3}\cdot {N}_{A}\cdot {R}_{{{{{{{{\rm{Nu}}}}}}}}}^{3}\cdot 1\,\mu\)M, where NA is the Avogadro constant. We used 500 particles in simulations, and the corresponding RNu = 0.58 μm matches well with the observed size of nucleoli11,13.

The size of the nucleolar particles can be determined assuming a space-filling model for the nucleoli as

$$\frac{\frac{4\pi }{3}\cdot {\left(\frac{{2}^{\frac{1}{6}}\cdot {\sigma }_{{{{P}}}}}{2}\right)}^{3}\cdot {N}_{{{{P}}}}}{\frac{4\pi }{3}\cdot {R}_{{{N}}\,}^{3}}={\left(\frac{{R}_{{{{Nu}}}}}{{R}_{{{{N}}}}}\right)}^{3}.$$

Using NP = 500, RN = 19.7, σ = 5 μm and RNu = 0.58 μm, we have the diameter of nucleolar particles σp = 0.5 σ. 21/6σp is the equilibrium distance between neighboring nucleolar particles in the Lennard–Jones potential.

We note that the above estimation is crude since many additional proteins and RNA molecules are present inside the nucleoli. Therefore, the nucleolar particles should be viewed as molecular aggregates with size σP = 0.5σ–0.12 μm, rather than a single protein molecule. Given the size of a typical protein as 5–10 nm86, the number of molecules within a single coarse-grained particle can be on the order of 103. This number, while large, is, in fact, in the same order as the amount of distinct molecules that make up the nucleoli87.

The molecular aggregate interpretation of nucleolar particles is also consistent with our use of ϵ = 1.8 kBT for specific interactions. This value appears much smaller than the strength of typical non-bonded interactions, such as hydrogen bonds. ϵ should be viewed as “free energy” that accounts for the averaging over the heterogeneous interaction pattern of different molecules and the entropic penalty arising from confining individual molecules together. Furthermore, as stated above, the effective temperature in our model is presumably higher than 310 K, resulting in more significant values for ϵ.

Simulation details for free energy calculations

We computed the free energy profile for coalescence using umbrella sampling and temperature replica exchange with 16 umbrella windows88,89. We defined the collective variable as \({R}_{g}=\sqrt{\frac{1}{N}\mathop{\sum }\nolimits_{i = 1}^{N}{\left|{{{{{{{{\bf{r}}}}}}}}}_{i}-{{{{{{{{\bf{r}}}}}}}}}_{{{\mbox{com}}}}\right|}^{2}}\). ri is the Cartesian coordinate of the i-th nuclear particles found inside one of the two droplets and rcom is the center of mass. Indices of nucleolar particles in the two droplets were identified using the DBSCAN algorithm (see Supplementary Materials). A harmonic potential \(\frac{K}{2}\cdot {\left({R}_{g}-{R}_{c}\right)}^{2}\) was introduced in each umbrella window to facilitate the sampling of configurations at targeted distances. We chose the center of these windows, Rc, to be evenly spaced between 2.5 σ and 6.0 σ with an increment of 0.25 σ, except for the first window whose Rc = 2.0 σ. The spring constant K was chosen as follows:

$$K=\left\{\begin{array}{ll}100.0&2.00\le {R}_{0} \, < \, 3.50\\ 150.0&3.50\le {R}_{0}\, < \, 4.00\\ 200.0&{R}_{0}\ge 4.00\hfill\end{array}\right.$$

Eight replicas were used within each umbrella window with temperatures ranging from 1.00 to 1.14 with 0.02 increments. Exchanges between these replicas were attempted at every 100-time steps. As shown in Supplementary Fig. 21, our choice of the temperature grid allows frequent exchange among replicas.

These simulations were initialized from a typical two-droplet configuration recorded at the end of a dynamical simulation and lasted twelve million steps. Configurations were recorded at every 400 steps. To compute the error bars and evaluate simulation convergence, we divided the simulation trajectories into five consecutive blocks with equal length. Free energy profiles were then calculated using only data collected from each block with the weighted histogram method (WHAM), and error bars were determined as the standard deviation of the five independent estimations.

Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data availability

The data that support this study are available from the corresponding author upon reasonable request. Simulation trajectories generated in this study have been deposited in the Zenodo database ( Source data are provided with this paper.

Code availability

The computer code used in this study is available at


  1. Strom, A. R. & Brangwynne, C. P. The liquid nucleome–phase transitions in the nucleus at a glance. J. Cell Sci. 132, jcs235093 (2019).

    CAS  PubMed  PubMed Central  Google Scholar 

  2. Henninger, J. E. et al. RNA-mediated feedback control of transcriptional condensates. Cell 184, 207–225 (2021).

    CAS  PubMed  Google Scholar 

  3. Dundr, M. Nuclear bodies: multifunctional companions of the genome. Curr. Opin. Cell Biol. 24, 415–422 (2012).

    CAS  PubMed  PubMed Central  Google Scholar 

  4. Banani, S. F., Lee, H. O., Hyman, A. A. & Rosen, M. K. Biomolecular condensates: organizers of cellular biochemistry. Nat. Rev. Mol. Cell Biol. 18, 285–298 (2017).

    CAS  PubMed  PubMed Central  Google Scholar 

  5. Brangwynne, C. P. et al. Germline p granules are liquid droplets that localize by controlled dissolution/condensation. Science 324, 1729–1732 (2009).

    ADS  CAS  PubMed  Google Scholar 

  6. Mao, Y. S., Zhang, B. & Spector, D. L. Biogenesis and function of nuclear bodies. Trends Genet. 27, 295–306 (2011).

    CAS  PubMed  PubMed Central  Google Scholar 

  7. Boeynaems, S. et al. Protein phase separation: a new phase in cell biology. Trends Cell Biol. 28, 420–435 (2018).

    CAS  PubMed  PubMed Central  Google Scholar 

  8. Shin, Y. & Brangwynne, C. P. Liquid phase condensation in cell physiology and disease. Science 357, eaaf4382 (2017).

    PubMed  Google Scholar 

  9. Lafontaine, D. L., Riback, J. A., Bascetin, R. & Brangwynne, C. P. The nucleolus as a multiphase liquid condensate. Nat. Rev. Mol. Cell Biol. 1–18 (2020).

  10. Courchaine, E. M., Lu, A. & Neugebauer, K. M. Droplet organelles? EMBO J. 35, 1603–1612 (2016).

    CAS  PubMed  PubMed Central  Google Scholar 

  11. Caragine, C. M., Haley, S. C. & Zidovska, A. Nucleolar dynamics and interactions with nucleoplasm in living cells. Elife 8, e47533 (2019).

    CAS  PubMed  PubMed Central  Google Scholar 

  12. Brangwynne, C. P., Mitchison, T. J. & Hyman, A. A. Active liquid-like behavior of nucleoli determines their size and shape in Xenopus laevis oocytes. Proc. Natl. Acad. Sci. 108, 4334–4339 (2011).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  13. Caragine, C. M., Haley, S. C. & Zidovska, A. Surface fluctuations and coalescence of nucleolar droplets in the human cell nucleus. Phys. Rev. Lett. 121, 148101 (2018).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  14. Kim, J., Han, K. Y., Khanna, N., Ha, T. & Belmont, A. S. Nuclear speckle fusion via long-range directional motion regulates speckle morphology after transcriptional inhibition. J. Cell Sci. 132, jcs226563 (2019).

    CAS  PubMed  PubMed Central  Google Scholar 

  15. Berry, J., Weber, S. C., Vaidya, N., Haataja, M. & Brangwynne, C. P. Rna transcription modulates phase transition-driven nuclear body assembly. Proc. Natl. Acad. Sci. 112, E5237–E5245 (2015).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  16. Larson, A. G. et al. Liquid droplet formation by hp1α suggests a role for phase separation in heterochromatin. Nature 547, 236–240 (2017).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  17. Strom, A. R. et al. Phase separation drives heterochromatin domain formation. Nature 547, 241–245 (2017).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  18. Hnisz, D., Shrinivas, K., Young, R. A., Chakraborty, A. K. & Sharp, P. A. A phase separation model for transcriptional control. Cell 169, 13–23 (2017).

    CAS  PubMed  PubMed Central  Google Scholar 

  19. Chong, S. et al. Imaging dynamic and selective low-complexity domain interactions that control gene transcription. Science 361, eaar2555 (2018).

    PubMed  PubMed Central  Google Scholar 

  20. Gibson, B. A. et al. Organization of chromatin by intrinsic and regulated phase separation. Cell 179, 470–484 (2019).

    CAS  PubMed  PubMed Central  Google Scholar 

  21. Kashchiev, D. Nucleation (Elsevier, 2000).

  22. McGaughey, A. & Ward, C. Droplet stability in a finite system: consideration of the solid–vapor interface. J. Appl. Phys. 93, 3619–3626 (2003).

    ADS  CAS  Google Scholar 

  23. Farley, K. I., Surovtseva, Y., Merkel, J. & Baserga, S. J. Determinants of mammalian nucleolar architecture. Chromosoma 124, 323–331 (2015).

    CAS  PubMed  PubMed Central  Google Scholar 

  24. Heride, C. et al. Distance between homologous chromosomes results from chromosome positioning constraints. J. Cell Sci. 123, 4063–4075 (2010).

    CAS  PubMed  Google Scholar 

  25. Royo, F. et al. Spatial link between nucleoli and expression of the zac1 gene. Chromosoma 118, 711 (2009).

    CAS  PubMed  PubMed Central  Google Scholar 

  26. Kind, J. et al. Single-cell dynamics of genome-nuclear lamina interactions. Cell 153, 178–192 (2013).

    CAS  PubMed  Google Scholar 

  27. Grosch, M., Ittermann, S., Shaposhnikov, D. & Drukker, M. Chromatin-associated membraneless organelles in regulation of cellular differentiation. Stem Cell Rep. 15, 1220-1232 (2020).

    CAS  Google Scholar 

  28. Spector, D. L. & Lamond, A. I. Nuclear speckles. Cold Spring Harb. Perspect. Biol. 3, a000646 (2011).

    PubMed  PubMed Central  Google Scholar 

  29. Weber, S. C. & Brangwynne, C. P. Inverse size scaling of the nucleolus by a concentration-dependent phase transition. Curr. Biol. 25, 641–646 (2015).

    CAS  PubMed  PubMed Central  Google Scholar 

  30. Dillinger, S., Straub, T. & Nemeth, A. Nucleolus association of chromosomal domains is largely maintained in cellular senescence despite massive nuclear reorganisation. PLoS One 12, e0178821 (2017).

    PubMed  PubMed Central  Google Scholar 

  31. Martin, C. et al. Lamin b1 maintains the functional plasticity of nucleoli. J. Cell Sci. 122, 1551–1562 (2009).

    CAS  PubMed  PubMed Central  Google Scholar 

  32. Zwicker, D., Hyman, A. A. & Juelicher, F. Suppression of Ostwald ripening in active emulsions. Phys. Rev. E 92, 012317 (2015).

    ADS  Google Scholar 

  33. Falahati, H. & Wieschaus, E. Independent active and thermodynamic processes govern the nucleolus assembly in vivo. Proc. Natl. Acad. Sci. 114, 1335–1340 (2017).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  34. Wurtz, J. D. & Lee, C. F. Chemical-reaction-controlled phase separated drops: formation, size selection, and coarsening. Phys. Rev. Lett. 120, 078102 (2018).

    ADS  CAS  PubMed  Google Scholar 

  35. Söding, J., Zwicker, D., Sohrabi-Jahromi, S., Boehning, M. & Kirschbaum, J. Mechanisms for active regulation of biomolecular condensates. Trends Cell Biol. 30, 4–14 (2020).

    PubMed  Google Scholar 

  36. Lee, D. S., Wingreen, N. S. & Brangwynne, C. P. Chromatin mechanics dictates subdiffusion and coarsening dynamics of embedded condensates. Nature Physics. 17, 531–538 (2021).

    ADS  CAS  Google Scholar 

  37. Zhang, Y., Lee, D. S., Meir, Y., Brangwynne, C. P. & Wingreen, N. S. Mechanical frustration of phase separation in the cell nucleus by chromatin. Phys. Rev. Lett. 126, 258102 (2021).

    ADS  CAS  PubMed  Google Scholar 

  38. Ronceray, P., Mao, S., Košmrlj, A. & Haataja, M. P. Liquid demixing in elastic networks: cavitation, permeation, or size selection? arXivpreprint (2021).

  39. McStay, B. Nucleolar organizer regions: genomic ‘dark matter requiring illumination. Genes Dev. 30, 1598–1610 (2016).

    CAS  PubMed  PubMed Central  Google Scholar 

  40. Németh, A. & Längst, G. Genome organization in and around the nucleolus. Trends Genet. 27, 149–156 (2011).

    PubMed  Google Scholar 

  41. van Koningsbruggen, S. et al. High-resolution whole-genome sequencing reveals that specific chromatin domains from most human chromosomes associate with nucleoli. Mol. Biol. Cell 21, 3735–3748 (2010).

    PubMed  PubMed Central  Google Scholar 

  42. Boisvert, F.-M., van Koningsbruggen, S., Navascués, J. & Lamond, A. I. The multifunctional nucleolus. Nat. Rev. Mol. Cell Biol. 8, 574–585 (2007).

    CAS  PubMed  Google Scholar 

  43. Galganski, L., Urbanek, M. O. & Krzyzosiak, W. J. Nuclear speckles: molecular organization, biological function and role in disease. Nucleic Acids Res. 45, 10350–10368 (2017).

    CAS  PubMed  PubMed Central  Google Scholar 

  44. Khelifi, G. & Hussein, S. M. A new view of genome organization through RNA directed interactions. Front. Cell Dev. Biol. 8, 517 (2020).

    PubMed  PubMed Central  Google Scholar 

  45. Bronstein, I. et al. Transient anomalous diffusion of telomeres in the nucleus of mammalian cells. Phys. Rev. Lett. 103, 018102 (2009).

    ADS  CAS  PubMed  Google Scholar 

  46. Weber, S. C., Theriot, J. A. & Spakowitz, A. J. Subdiffusive motion of a polymer composed of subdiffusive monomers. Phys. Rev. E 82, 011913 (2010).

    ADS  Google Scholar 

  47. Kepten, E., Bronshtein, I. & Garini, Y. Ergodicity convergence test suggests telomere motion obeys fractional dynamics. Phys. Rev. E 83, 041919 (2011).

    ADS  CAS  Google Scholar 

  48. Burnecki, K. et al. Universal algorithm for identification of fractional Brownian motion. A case of telomere subdiffusion. Biophys. J. 103, 1839–1847 (2012).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  49. Lucas, J. S., Zhang, Y., Dudko, O. K. & Murre, C. 3D trajectories adopted by coding and regulatory DNA elements: first-passage times for genomic interactions. Cell 158, 339–352 (2014).

    CAS  PubMed  PubMed Central  Google Scholar 

  50. Qi, Y. et al. Data-driven polymer model for mechanistic exploration of diploid genome organization. Biophys. J. 119, 1905–1916 (2020).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  51. Johnstone, S. E. et al. Large-scale topological changes restrain malignant progression in colorectal cancer. Cell 182, 1474–1489 (2020).

    CAS  PubMed  PubMed Central  Google Scholar 

  52. Lieberman-Aiden, E. et al. Comprehensive mapping of long-range interactions reveals folding principles of the human genome. Science 326, 289–293 (2009).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  53. Rao, S. S. et al. A 3D map of the human genome at kilobase resolution reveals principles of chromatin looping. Cell 159, 1665–1680 (2014).

    CAS  PubMed  PubMed Central  Google Scholar 

  54. Zhang, B. & Wolynes, P. G. Topology, structures, and energy landscapes of human chromosomes. Proc. Natl. Acad. Sci. 112, 6062–6067 (2015).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  55. Qi, Y. & Zhang, B. Predicting three-dimensional genome organization with chromatin states. PLoS Comput. Biol. 15, e1007024 (2019).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  56. Phair, R. D. & Misteli, T. High mobility of proteins in the mammalian cell nucleus. Nature 404, 604–609 (2000).

    ADS  CAS  PubMed  Google Scholar 

  57. Chen, D. & Huang, S. Nucleolar components involved in ribosome biogenesis cycle between the nucleolus and nucleoplasm in interphase cells. J. Cell Biol. 153, 169–176 (2001).

    CAS  PubMed  PubMed Central  Google Scholar 

  58. Stevens, T. J. et al. 3D structures of individual mammalian genomes studied by single-cell Hi-C. Nature 544, 59–64 (2017).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  59. Vertii, A. et al. Two contrasting classes of nucleolus-associated domains in mouse fibroblast heterochromatin. Genome Res. 29, 1235–1249 (2019).

    CAS  PubMed  PubMed Central  Google Scholar 

  60. Cremer, T. & Cremer, M. Chromosome territories. Cold Spring Harb. Perspect. Biol. 2, a003889–a003889 (2010).

    PubMed  PubMed Central  Google Scholar 

  61. Bickmore, W. A. & van Steensel, B. Genome architecture: domain organization of interphase chromosomes. Cell 152, 1270–1284 (2013).

    CAS  PubMed  Google Scholar 

  62. Tjong, H. et al. Population-based 3D genome structure analysis reveals driving forces in spatial genome organization. Proc. Natl. Acad. Sci. 113, E1663–E1672 (2016).

    CAS  PubMed  PubMed Central  Google Scholar 

  63. Boyle, S. et al. The spatial organization of human chromosomes within the nuclei of normal and emerin-mutant cells. Hum. Mol. Genet. 10, 211–9 (2001).

    CAS  PubMed  Google Scholar 

  64. van Gunsteren, W. F., Daura, X. & Mark, A. E. Computation of free energy. Helvetica Chim. Acta 85, 3113–3129 (2002).

    Google Scholar 

  65. Parsons, T. & Zhang, B. Critical role of histone tail entropy in nucleosome unwinding. J. Chem. Phys. 150, 185103 (2019).

    ADS  PubMed  Google Scholar 

  66. Du, R., Pande, V. S., Grosberg, A. Y., Tanaka, T. & Shakhnovich, E. S. On the transition coordinate for protein folding. J. Chem. Phys. 108, 334–350 (1998).

    ADS  CAS  Google Scholar 

  67. Peters, B. Reaction coordinates and mechanistic hypothesis tests. Annu. Rev. Phys. Chem. 67, 669–690 (2016).

    ADS  CAS  PubMed  Google Scholar 

  68. Eggers, J., Lister, J. R. & Stone, H. A. Coalescence of liquid drops. J. Fluid Mech. 401, 293–310 (1999).

    ADS  MathSciNet  CAS  MATH  Google Scholar 

  69. Paulsen, J. D., Burton, J. C. & Nagel, S. R. Viscous to inertial crossover in liquid drop coalescence. Phys. Rev. Lett. 106, 114501 (2011).

    ADS  PubMed  Google Scholar 

  70. Flory, P. J. Principles of Polymer Chemistry (Cornell University Press, 1953).

  71. Ester, M. et al. A density-based algorithm for discovering clusters in large spatial databases with noise. In KDD, vol. 96, 226–231 (1996).

  72. Onuki, A. Phase Transition Dynamics (Cambridge University Press, 2002).

  73. Ratke, L. & Voorhees, P. W.Growth and Coarsening: Ostwald Ripening in Material Processing (Springer Science & Business Media, 2013).

  74. Stanich, C. A. et al. Coarsening dynamics of domains in lipid membranes. Biophys. J. 105, 444–454 (2013).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  75. Chen, Y. et al. Mapping 3D genome organization relative to nuclear compartments using TSA-Seq as a cytological ruler. J. Cell Biol. 217, 4025–4048 (2018).

    CAS  PubMed  PubMed Central  Google Scholar 

  76. Su, J.-H., Zheng, P., Kinrot, S. S., Bintu, B. & Zhuang, X. Genome-scale imaging of the 3D organization and transcriptional activity of chromatin. Cell 182, 1641–1659 (2020).

    CAS  PubMed  PubMed Central  Google Scholar 

  77. Zhang, L. et al. TSA-seq reveals a largely conserved genome organization relative to nuclear speckles with small position changes tightly correlated with gene expression changes. Genome Res. 31, 251–264 (2021).

    PubMed Central  Google Scholar 

  78. Fei, J. et al. Quantitative analysis of multilayer organization of proteins and RNA in nuclear speckles at super resolution. J. Cell Sci. 130, 4180–4192 (2017).

    CAS  PubMed  PubMed Central  Google Scholar 

  79. Huang, K. et al. Physical and data structure of 3D genome. Sci. Adv. 6, eaay4055 (2020).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  80. Style, R. W. et al. Liquid-liquid phase separation in an elastic network. Phys. Rev. X 8, 011028 (2018).

    CAS  Google Scholar 

  81. Rosowski, K. A. et al. Elastic ripening and inhibition of liquid-liquid phase separation. Nat. Phys. 16, 422–425 (2020).

    CAS  PubMed  PubMed Central  Google Scholar 

  82. Zidovska, A. Chromatin: liquid or solid? Cell 183, 1737–1739 (2020).

    CAS  PubMed  Google Scholar 

  83. Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 117, 1–19 (1995).

    ADS  CAS  MATH  Google Scholar 

  84. Jiang, Z. & Zhang, B. Theory of active chromatin remodeling. Phys. Rev. Lett. 123, 208102 (2019).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  85. Zhu, L. et al. Controlling the material properties and rRNA processing function of the nucleolus using light. Proc. Natl Acad. Sci. 116, 17330–17335 (2019).

    CAS  PubMed  PubMed Central  Google Scholar 

  86. Lee, H. H. et al. Crystal structure of human nucleophosmin-core reveals plasticity of the pentamer–pentamer interface. Proteins: Struct. Funct. Bioinforma 69, 672–678 (2007).

    CAS  Google Scholar 

  87. Scherl, A. et al. Functional proteomic analysis of human nucleolus. Mol. Biol. Cell 13, 4100–4109 (2002).

    CAS  PubMed  PubMed Central  Google Scholar 

  88. Sugita, Y. & Okamoto, Y. Replica-exchange molecular dynamics method for protein folding. Chem. Phys. Lett. 314, 141–151 (1999).

    ADS  CAS  Google Scholar 

  89. Souaille, M. & Roux, B. Extension to the weighted histogram analysis method: combining umbrella sampling with free energy calculations. Computer Phys. Commun. 135, 40–57 (2001).

    ADS  CAS  MATH  Google Scholar 

Download references


We thank Dr. Xingcheng Lin, Dr. Kartik Kamat, and Dr. Xinqiang Ding for helpful discussions. This work was supported by the National Institutes of Health (Grant R35GM133580 to B.Z.) and by the Gordon and Betty Moore Foundation (Grant GBMF9162.15 to B.Z.).

Author information

Authors and Affiliations



Y.Q. and B.Z. designed and performed the research, analyzed the data, and wrote the paper.

Corresponding author

Correspondence to Bin Zhang.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review information Nature Communications thanks the anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Qi, Y., Zhang, B. Chromatin network retards nucleoli coalescence. Nat Commun 12, 6824 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.


Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing