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.

Introduction

Genomic DNA is packed and folded three-dimensionally in the cell nuclei. In the nuclei of eukaryotic cells, the nucleosome is a basic unit consisting of an approximately 147-bp DNA wrapped around core histones [1]. Recent experimental evidences suggest that the nucleosome is irregularly folded without the 30-nm chromatin fiber [2-7]. On the other hand, at the scale of the whole nucleus, interphase chromosomes occupy distinct chromosome territories [8]. This highly organized chromosome structure allows for effective regulation of various genome functions.

By virtue of recent developments of chromosome conformation capture (3C) techniques, the genome-wide chromosome organization has been revealed by detecting the physical contact frequencies between pairs of genomic loci [9]. More recently, 3C derivatives, Hi-C and 5C profiles demonstrated that metazoan genomes are partitioned into submegabase-sized chromatin domains (CDs) including topologically associating domains (TADs) [10-12]. TADs are considered to be a regulatory and structural unit of the genome [13]; genome loci located in the same TAD are associated with each other, whereas genomic interactions are sharply depleted between adjacent domains. For even single-cell Hi-C, individual chromosomes maintain domain organization [14]. Furthermore, kilobase-resolution in situ Hi-C maps identified not only small contact domains but also CTCF-mediated loop domains [15, 16].

In contrast, dynamic aspects of chromatin have been shown by live-cell imaging experiments [17-24]. In particular, single-nucleosome imaging in living mammalian cells has revealed local nucleosome fluctuations caused by the thermal random force [25-27]. The mean-squared displacement (MSD) of dynamic nucleosome movement clearly shows subdiffusive motion, where Dapp is the apparent diffusion coefficient with dimension m2/s. This means that nucleosome movement must be affected by restrictions from some factors but thermal noise. Therefore, there must be a way that the dynamic aspect is consistent with aspects of the genome organization. A theory is required to relate the dynamic aspects described by Dapp and β to the structural features of CDs. To date, the subdiffusive exponent β has been considered to depend on the folding structure of nucleosome fibers [28] and the viscoelasticity of the thermal environment [29, 30].

The fractal nature of chromatin architecture as well as nucleus environment has been revealed recently [9, 31, 32]. The topological structure of CDs can be described by use of the fractal manner. Here, we propose a polymer model for a CD, whose conformational state is assumed to be expressed by the fractal dimension df in a viscoelastic medium with the exponent 0 < α < 1. Although not only the strings and binders switch model [33] but also the block copolymer model [34] can explain aspects of chromatin folding and chromosome architecture in Hi-C experiment datasets, in our model we abstract information on the conformational states of CDs and interpret their dynamic features by using size scaling according to the fractal dimensions. Accordingly, the analytical form of the MSD of nucleosomes in CDs can be derived in terms of polymer physics. As a result, the structural information of CDs, such as the size and conformational state expressed by the fractal dimension, can be derived from the MSD data of dynamic nucleosomes.

Results

Polymer model

CDs characterized by fractal dimensions

To construct a model of CDs, we assumed that a nucleosome fiber is represented as a polymer bead chain and forms a CD with size scaling, (Fig 1A), where N is the number of nucleosome beads in the CD, and 〈⋅〉CD represents the average for all nucleosome beads within the CD at thermal equilibrium. In polymer physics, the exponent 1/df corresponds to the size exponent ν [35, 36]. A nucleosome fiber in a CD not only has the excluded volume as a physical polymer, but also forms chromatin loops for transcriptional regulation [15, 16, 37]. Therefore, nucleosome fibers can interact with each other within the same CD through both attractive and repulsive interactions. Here, we assume that the effective conformational state of CDs is phenomenologically represented by the fractal dimension. Note that the states with df = 1, 2, and 3 correspond to a straight line, the ideal chain [35], and the fractal globule [9, 28, 38], respectively (Fig 1A).

Fig 1

Schematic illustration of our polymer model for CDs.

(A) A nucleosome fiber is represented as a polymer bead chain and forms a CD. The size scaling of CDs is expressed as , where the fractal dimension represents the effective conformational state of CDs: df = 1, 2, and 3 correspond to a straight line, the ideal chain, and the fractal globule, respectively. (B) The viscoelasticity of the medium, where the movement of particles shows the subdiffusive FBM 〈Δ(t)2〉∼t, is described using the friction coefficient with memory, γ(t)∼t−. When a nucleosome with coordinates (t) dynamically fluctuates in the viscoelastic medium, the movement of nucleosomes in CDs shows subdiffusion: 〈Δ(t)2〉∼t.

Nucleosome fiber fluctuation in viscoelastic medium with memory

The subdiffusive motion of tracer particles in living cells, 〈[(t) − (0)]2〉∼t, has been observed [29, 39–41]. There are several physical models for generating subdiffusion, including: (i) the generalized Langevin equation (GLE), which is consistent with fractional Brownian motion (FBM) [42-45], and (ii) the continuous-time random walk [46]. Since some experiments have shown that the movement of chromosomal loci displays the FBM [23, 29], here, we adopt the former model to describe the friction effect with memory in the viscoelastic medium [39, 47, 48] that satisfies the fluctuation-dissipation relation (FDR) [35, 49]: the GLE , with friction coefficient with memory of γ(t) ∼ t−, generates the subdiffusive FBM. The thermal random force (t) satisfies the FDR 〈g(t)gλ(t′)〉 = kBTγ(t − t′)δ, where k is the Boltzmann constant, T is the temperature of the environment, and the suffixes κ and λ represent x, y and z.

Here, we focus on the concrete description of our polymer model. A CD is assumed to be formed by N + 1 nucleosome beads at positions {0, 1, ⋯, } (Fig 1B), and adjacent beads are connected via a harmonic spring so that the effective bond length is beff, and long-range interactions exist such that the phenomenological size scaling of the CDs is proportional to . Moreover, as mentioned above, the friction effect between each nucleosome and the viscoelastic medium is assumed to be described by the friction coefficient with memory [30, 44, 45, 47], where the dimension of the coefficient γ is kg/s2−, and the Laplace transform of γ(t) has a simple form γ s (see Eq S19 in S1 Text). In the continuous limit [35], the Langevin equation of nucleosomes is described as where the long-range interaction force ((n,t) including attractive and repulsive interactions results in the size scaling and the thermal random force (n, t) satisfies the FDR: 〈g(n, t)gλ(m, t′)〉 = kB Tγ(t − t′)δ(n − m)δ. Our model for df = 2 formally corresponds to the classical Rouse model in the viscoelastic medium [30], where the force ((n, t) apparently vanishes. Hence, the additional long-range interaction force generating the scaling (Eq 4) has an important role in our model, and enables us to calculate the MSD analytically. Here, we do not take into account the hydrodynamic interactions between nucleosomes, which are discussed in Discussion and S1 Text, Section II.

Analytical calculation shows that the MSD of nucleosomes within fractal CDs is subdiffusive in viscoelastic medium

A standard approach for treating Eq 3 is to use the normal coordinates for p = 0, 1, 2, ⋯; however, the nonlinearity of the long-range interaction makes it difficult to deal with the equation in this manner. Therefore, to simplify the analysis, firstly, we assume that nucleosome fluctuations within the CD reach thermal equilibrium after the relaxation time τ, which is explicitly described below (Eqs 11 and 12). Second, we use an approximation to transform the nonlinear Langevin equation (Eq 3) into a linear equation by averaging under thermal equilibrium with respect to the normal coordinates The term in the left hand side and the second term in the right hand side (RHS) are straightforwardly derived according to the normal coordinates, in which satisfies 〈(t)〉 = 0 and the FDR (see S1 Text, Section IA). Instead of the linearity of Eq 5, the parameter k implicitly includes the nonlinear effect such as the long-range interactions, and is determined by the variance of over the thermal relaxation time [30] (see S1 Text, Section IB): Finally, to calculate the thermal average , the effective size scaling (Eq 4) generated by the long-range interactions is used. The asymptotic form for large p is calculated as follows (see S1 Text, Section IC): A is a dimensionless constant depending on the fractal dimension: . We shall refer to the above approximation as the linearization approximation, which is on the same level of the approximation as the preaveraging approximation in terms of polymer physics [35, 50]. From this point forward, to avoid complicated expressions caused by this asymptotic form, we regard the asymptotic sign ‘≃’ as equality.

Next, let us consider the MSD of nucleosomes in CDs. Since the inverse transform of normal coordinates is and the correlation between different modes vanishes, the MSD of the n-th nucleosome, ϕ(n, t)≡〈[(n, t) − (n, 0)]2〉, is expressed as where the correlation function is defined as C(t)≡〈(t) ⋅ (0)〉. Multiplying Eq 5 by (0) and averaging with 〈(t) ⋅ (0)〉 = 〈(t)〉⋅〈(0)〉 = 0, we can derive that the correlation function for p ≥ 1 satisfies The first term for p = 0 in the RHS of Eq 8 corresponds to the MSD of the center of the CD, and the motion obeys and the FDR . According to the fluctuation-dissipation theorem [49], the motion of the center of mass is subdiffusive with exponent α (see S1 Text, Section IE): where represents the relaxation time of nucleosome fluctuations in the CD.

On the other hand, the second term in the RHS of Eq 8 describes the fluctuations of many modes inside the CD. Using the Laplace transformation and the thermal equilibrium initial state, the solution of Eq 9 can be derived as follows (see S1 Text, Section ID): where E(x) is the Mittag-Leffler function. According to the polymer physics [35] for t ≪ τ, ϕ(n, t) is dominated by terms with large p. Moreover, since the MSD in our experiment (Fig 2E) is calculated by averaging the nucleosome trajectories at various positions in CDs, the term can be replaced by the average 1/2. Therefore, according to the asymptotic form of the Mittag-Leffler function, E(−x) ≃ exp[−x/Γ(1 + α)] for x ≪ 1, and the conversion of the sum into the integral, we obtain for t ≪ τ where is a dimensionless constant (see S1 Text, Section IF). Thus, in our model, subdiffusive motion of single nucleosomes is a typical feature, assuming both fractal CDs and viscoelastic medium.

Fig 2

Single-nucleosome imaging and analysis.

(A) Single-nucleosome image of a human HeLa cell nucleus expressing H2B-PA-mCherry. Each dot represents single nucleosome. (B) Evidence that each dot represents single-nucleosome molecule. Each H2B-PA-mCherry dot shows single-step photobleaching. The vertical axis represents the fluorescence intensity of each H2B-PA-mCherry dot. The horizontal axis is the tracking time series (each photobleaching point is set as time 0; the average and the standard deviation at each time point were calculated for 50 dots.). Due to the clear single-step photobleaching profile of the H2B-PA-mCherry dots, each dot shows a single H2B-PA-mCherry molecule in a single nucleosome. (C) A scheme for nuclear interior (Top) and periphery (Bottom) imaging. Focal plane (red) in the living cells is shown. See also S1 Fig. (D) Representative trajectories of fluorescently labeled single nucleosome (50 ms per frame). (E) Plots of the MSD at the interior and periphery regions. These fit well with the MSD curves using Eq 1.

Nucleosome movement is much greater in the nuclear interior than at the nuclear periphery

In order to apply our model to living human cells, single-particle imaging of nucleosomes was performed by observation of PA-mCherry labels [51] attached to histone H2B in human HeLa cells (Fig 2A). The clear single-step photobleaching profile of the H2B-PA-mCherry dots shows a single H2B-PA-mCherry molecule in a single nucleosome (Fig 2B). We tracked approximately 40,000 dots representing single nucleosomes (S1 Table). Fig 2D shows representative trajectories of the dynamic nucleosome movement in single cells.

Here, to evaluate the state of CDs according to their position in the nucleus, we focused on the nuclear interior and periphery (or surface) (Fig 2C and S1 Fig), and calculated the MSD. The nuclear periphery is a heterochromatin-rich region, which presumably shows much less active transcription than the interior. The plots of the MSD at each region, in time interval t up to 0.5 s, are shown in Fig 2E (normal scale) and S2 Fig (log-log scale) (also see S1 Table). The MSD at the interior is higher than that at the periphery. This result implies that nucleosome movement within CDs in the euchromatin-rich interior region is higher than that in the heterochromatin-rich periphery region.

As we analytically derived the subdiffusive MSD (Eq 13), the experimental result clearly shows subdiffusion of single-nucleosomes: using Eq 1, the plots fit well with the MSD curves 0.018 t0.44 μm2 and 0.013 t0.39 μm2 for the interior and the periphery, respectively.

MSD is lower at the nuclear periphery than the interior, indicating that heterochromatin-rich CDs are more compact

Comparing Eqs 1 and 13, β and Dapp are calculated as where . It turns out that these values contain statistical information of the CD structures, 〈R〉CD and df. Since β and Dapp can be determined by the fitting in our experiments, we can therefore estimate 〈R〉CD and df, inversely.

The lower MSD at the periphery than at the interior, Dapp,periphery < Dapp,interior and βperiphery < βinterior, reflects the fact that the CDs near the periphery are in a more compact conformational state and are smaller in size than those at the interior: df,periphery > df,interior and 〈R〉CD,periphery < 〈R〉CD,interior. This property is consistent with the conventional distribution of heterochromatin: the CDs in the heterochromatin-rich nuclear periphery are more compact than those in the euchromatin-rich interior [52].

Discussion

To estimate the structural information of CDs through solving Eqs 14 and 15 inversely, the values of N, α, and γ in mammalian living cell nuclei are required. The average size of TADs was determined to be 880 kb from mouse embryonic stem cells (mESCs), with a range of 100 kb to 5 Mb [10]. Here, we assume a CD size of 1 Mb, which corresponds to 〈N〉CD = 5000 nucleosomes. To the best of our knowledge, few studies have estimated the friction effect in viscoelastic cell nuclei. Therefore, we use the value of the diffusion coefficient of enhanced green fluorescent protein (EGFP)-monomer around interphase chromatin, DEGFP = 20.6 μm2/s [25], measured by fluorescence correlation spectroscopy, in which α is assumed to be 1. In general, as a result of the FDR in a viscoelastic medium with α, the diffusion coefficient of a diffusive particle for one degree of freedom is kBT/[Γ(1 + α) ⋅ γ] (see Eq. S34 in S1 Text). Since the contribution of Γ(1 + α) is within the range 1 ≤ 1/Γ(1 + α)<1.13 for 0 < α ≤ 1, the friction coefficient of EGFP in the nucleus can be approximately regarded as the diffusion coefficient as γ = kBT/DEGFP. The hydrodynamic radius of a nucleosome bead with an H2B-PA-mCherry is assumed to be approximately quadruple for the EGFP. This means that the friction effect is also 4 times larger [48]. Accordingly, we use γ = 4kBT/DEGFP. Finally, the structural information of CDs is estimated by calculating

β could be measured in our experiment, although the value of α could not be determined simultaneously. Hence, Eq 16 represents the relationship between α and df, as shown in Fig 3A. Under this constrained condition, according to Eqs 16 and 17, the values of the structural information within the nuclear interior and periphery regions are calculated and mapped as a function of α (Fig 3B). Since fluorescence correlation spectroscopy measurements of GFP have shown that the value of α is close to 0.79 in not HeLa but NRK nuclei [31], as an example, we summarize the estimated values for α = 0.8 and α = 0.9 in Table 1. The exponent β = 0.4 for the fractal globule model [28] corresponds to the value for df = 3 and α = 1 in Eq 14. Furthermore, our previous results have shown smaller exponents β = 0.37 and 0.31 for interphase chromatin and mitotic chromosome, respectively [25]. Unless considering the case of 0 < α < 1, this smaller exponent cannot be explained. Note that α has only minor effects on C (see S3 Fig).

Fig 3

Structural information of CDs provided by single-nucleosome dynamics.

(A) The fractal dimension df within the nuclear interior region (β = 0.44) and the periphery region (β = 0.39) for α, according to Eq 16. (B) df and the size 〈R〉CD of CDs, and (C) df and the relaxation time τ of nucleosomes in CDs within the nuclear interior region (colored circle) and the periphery region (colored square) calculated for various α values, according to Eqs 16, 17 and 18. (D) 〈R〉CD and the τ within the nuclear interior region (colored circle) and the periphery region (colored square) calculated for α = 0.9 and various 〈N〉CD values, corresponding to the range of 200 kb to 4 Mb, according to Eqs 17 and 18.

Table 1

Estimated values of the fractal dimension df, the size 〈R〉CD, and the relaxation time τ for α = 0.8 and α = 0.9.

Region	β	Dapp (μm 2/sβ)	α = 0.8			α = 0.9
			df	〈R〉CD (nm)	τdf,α (s)	df	〈R〉CD (nm)	τdf,α (s)
Interior	0.44	0.018	1.64	358	4.65	2.09	268	2.69
Periphery	0.39	0.013	2.10	191	1.31	2.62	166	1.01

The exponent β and the apparent diffusion coefficient Dapp are obtained by fitting of the MSD results at the interior and periphery regions (Fig 2C). In the calculation, we used the following values: 〈N〉CD = 5000 nucleosomes, γ = 4kBT/DEGFP and DEGFP = 20.6 μm2/s.

The relaxation time of nucleosomes in CDs is calculated as and is mapped as a function of α and df (Fig 3C). The short relaxation time (∼ s) means that the thermal equilibrium, which is the precondition for the linearization approximation, were fulfilled in our experiments. In measurements of long-term single-nucleosome movements, the MSD is expected to show a transition toward movement of the center of CDs with the exponent α (Eq 10). This would enable estimating α, df, 〈R〉CD, and τ without requiring the use of the assumptive values described above, such as 〈N〉CD and DEGFP. The long-term (≫τ) imaging of chromatin dynamics in mammalian nuclei might reveal this transition motion [19, 21, 24].

As mentioned at the beginning of this section, the measured TAD size of mESCs is in the range of 100 kb to 5 Mb. Fig 3D shows the relationship between 〈R〉CD and τ for α = 0.9 as a function of 〈N〉CD, corresponding to the range of 200 kb to 4 Mb, according to Eqs 17 and 18. The relaxation time within several tens of seconds is consistent with the assumption of the linearization approximation as mentioned above. Moreover, the estimated CD size within 100–500 nm is also consistent with observed radius for chromatin domains as detected by super-resolution imaging [53].

The critical assumption of the linearization approximation is that nonlinear and complicated long-range interactions can be replaced by the mean-field fluctuation near thermal equilibrium within a fractal CD. Our result, that the estimated CD size is about 100–500 nm and the relaxation time is at most a few seconds, implies that the condition of the approximation is fulfilled. This kind of approximations has been discussed well in polymer physics [35], where the results of the approximation are not much different from those of more sophisticated calculations including the renormalization group theory. Furthermore, we have already reported that thermal fluctuation plays a dominant role in chromatin dynamics within CDs during a few seconds of observation for each fluorescent nucleosome [25-27]. On the other hand, successive ATP-dependent active processes on chromatin might play an important role in chromosome folding during mitosis [54]. In such a case, we cannot apply our theory due to the non-equilibrium nature.

In addition, non-equilibrium fluctuations driven by ATP-dependent cell activities affect chromatin dynamics [17, 20]. In order to directly take the effect into account, we have to add non-equilibrium fluctuation noise term to Eq 3. For a polymer with specific non-equilibrium fluctuations, where the correlation of the added athermal noise exponentially decays, a theoretical result has been obtained [55]. However, properties of non-equilibrium fluctuations on chromatin are unclear. It was suggested that a decrease of the chromatin persistence length of a CD may occur due to ATP-driven nucleosome remodeling [56]. We might consider an effect of ATP-driven remodeling on changes of the domain size and the fractal dimension.

Recent high-resolution Hi-C and ChIA-PET data have shown that architecture proteins such as CTCF and cohesin play important roles in CD organization at the boundaries [15, 16, 57]. However, because of experimental limitations including specific nucleosome labeling and microscopy resolution, we could not distinguish the nucleosome movements between the domain center and boundaries. Furthermore, different epigenetic states [53] including posttranslational modifications [58] affect the spatial organization of chromatin domains. These effects seem to regulate nucleosome-nucleosome interactions. If we could directly observe the nucleosome movements depending on the interactions in vivo using a novel labeling technology, we would be able to extend our framework; which will be a challenging issue in the future.

Here, we considered a locally clustered polymer with effective size scaling (Eq 4) in the absence of hydrodynamic interactions (HIs) as a model of CDs. The inverse proportion of k to N, except for the contribution from , in Eq 6 reflects the lack of HIs in our model; that is, the hydrodynamic field goes through nucleosome beads without interactions. The hydrodynamic effect of surface monomers in a polymer blob on the exponent β has been argued in [28]. Applying their discussion to our results, Eq 16 changes into df = c(2α/β − 2), where the coefficient c is within the range 1 ≤ c < 1.09. The effect is expected to be small. One can also consider a polymer model including HIs, which would affect the mobility matrix and work cooperatively within a polymer blob [35, 50]. In such a situation, the HI cancels out the effect of the size scaling described by the fractal dimension df: β = α ⋅ 2/3, and β does not depend on df (see S1 Text, Section II).

Conclusion

Our results indicate that our proposed model serves as a strong method for extracting the structural information of CDs from observations of dynamic nucleosome movement. Super-resolution microscopy techniques can be used to elucidate the spatial size of CDs according to different epigenetic states [53]. On the other hand, development of an effective imaging technique to reveal the fractal dimensions remains a challenge for the future. The conformational state of CDs characterized by the fractal dimension must be associated with the accessibility of transcription factors, depending on the physical size of those factors [59]. Beyond the pioneer computational work of analyzing interphase chromosomes based on the chromatin fibers [60], further development of not only a large-scale chromosome model based on the results of a genome-wide association study [61] but also restraint-based three-dimensional modeling of genomes [62] is expected to provide novel insight and open the door toward further discovery on the relationship between dynamic genome organization and stochastic gene expression.

Materials and Methods

Cell isolation and culture

To observe single nucleosomes and analyze their local dynamics in living human cells, histone H2B was fused with photoactivatable (PA)-red fluorescent protein (mCherry) [51] and expressed in HeLa cells as described previously [25]. The cell lines expressing H2B-PA-mCherry at a very low level were isolated. The cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% fetal bovine serum (FBS) (vol/vol) at 37°C in 5% CO2 (vol/vol). The cells were plated 24–48 h before the experiment onto Iwaki glass bottom dishes treated with poly-lysine. Before the experiment, the medium was replaced by DMEM F-12 (non phenol red) with 15% FBS. The cells were then set on the microscope stage kept in a custom-built 37°C microscope incubator enclosure with 5% CO2 (vol/vol) delivery throughout the experiment.

Microscopy

For single-nucleosome imaging, an oblique illumination microscope was used to illuminate a limited thin area within the cell (Nikon laser TIRF microscope system Ti with sapphire 564-nm laser). In general, PA-mCherry exhibits red fluorescence only after activation by a 405-nm laser [51]. However, we unexpectedly found that a relatively small number (∼100/time frame/nucleus) of H2B-PA-mCherry molecules were continuously and stochastically activated even without UV laser stimulation. Fig 2A shows a typical single-nucleosome image of a living HeLa cell. Due to the clear single-step photobleaching profile of the H2B-PA-mCherry dots, each dot in the nucleus represents a single H2B-PA-mCherry in a single nucleosome (Fig 2B). Nucleosome signals were recorded in the interphase chromatin of the nuclear interior and periphery in living HeLa cells at a frame rate of ca. 50 ms/frame. Note that the two different focal planes for the nuclear interior and periphery (Fig 2C) were precisely ensured by nuclear surface labeling with Nup107 (a nuclear pore component)-Venus (a bright yellow fluorescent protein) [63] (see S1 Fig).

Tracking and data analysis

Local nucleosome fluctuation was observed (ca. 60 nm movement/50 ms), presumably caused by Brownian motion. The free MATLAB software u-track [64] was used for single-nucleosome tracking. The dots were fitted to an assumed Gaussian point spread function to determine the precise center of the signals with higher resolution. Finally, we obtained data set of two-dimensional M trajectories , where the suffix j ∈ {1, ⋯, M} represents the sample number for the tracked time-interval [0, t]; t ≡ i × 50 ms. Several representative trajectories of fluorescently tagged single nucleosomes are shown in Fig 2D (bar = 100 nm).

According to observed regions, we calculated the ensemble-averaged MSD of single nucleosomes: . Here, in order to obtain the three-dimensional value, we multiplied the two-dimensional value by 3/2 on the assumption of isotropy. Plots of the MSDs of single nucleosomes in interphase chromatin at the nuclear interior (10 cells) and the nuclear periphery (10 cells) from 0 to 0.5 s are shown in Fig 2E. The plots for single nucleosomes were fitted with the subdiffusion model (Eq 1) using R-software. The standard error of the mean (SEM), which is the standard deviation of the sampling distribution of the mean, for MSD(t) was sufficiently small. The number of trajectories M and the SEM of MSD(t) are summarized in S1 Table.

Further details on derivations of the theoretical results and remarks on the hydrodynamic effect for the model.

(PDF)

Click here for additional data file.

A schematic representation for nuclear interior (Top left) and periphery (Top right) imaging.

Illumination laser (green) and focal plane (red) in the living cells are shown. Note that the two different focal planes were precisely verified by nuclear surface labeling with Nup107 (a nuclear pore component)-Venus (a bright yellow fluorescent protein) [63]. The nuclear rim signals (Bottom left) and dot signals in ellipse shape (Bottom right) show the middle layer of nucleoplasm and the nuclear surface, respectively. Bar shows 5 μm.

(TIF)

Click here for additional data file.

Plots of the MSD (Fig 2E) on log-log scale.

(EPS)

Click here for additional data file.

The function C of the fractal dimension df for α = 0.6, 0.7, 0.8, 0.9, and 1.0.

α has only a slight effect on C.

(EPS)

Click here for additional data file.

The number of tracked trajectories M and the standard error of the mean (SEM) of the MSD at the nuclear interior region and the periphery region.

The measurements at each region were performed using 10 cells.

(PDF)

Click here for additional data file.