A new application of the phase-field method for understanding the mechanisms of nuclear architecture reorganization

Specific features of nuclear architecture are important for the functional organization of the nucleus, and chromatin consists of two forms, heterochromatin and euchromatin. Conventional nuclear architecture is observed when heterochromatin is enriched at nuclear periphery, and it represents the primary structure in the majority of eukaryotic cells, including the rod cells of diurnal mammals. In contrast to this, inverted nuclear architecture is observed when the heterochromatin is distributed at the center of the nucleus, which occurs in the rod cells of nocturnal mammals. The inverted architecture found in the rod cells of the adult mouse is formed through the reorganization of conventional architecture during terminal differentiation. Although a previous experimental approach has demonstrated the relationship between these two nuclear architecture types at the molecular level, the mechanisms underlying long-range reorganization processes remain unknown. The details of nuclear structures and their spatial and temporal dynamics remain to be elucidated. Therefore, a comprehensive approach, using mathematical modeling, is required, in order to address these questions. Here, we propose a new mathematical approach to the understanding of nuclear architecture dynamics using the phase-field method. We successfully recreated the process of nuclear architecture reorganization, and showed that it is robustly induced by physical features, independent of a specific genotype. Our study demonstrates the potential of phase-field method application in the life science fields. Electronic supplementary material The online version of this article (doi:10.1007/s00285-016-1031-3) contains supplementary material, which is available to authorized users.


Introduction
In eukaryotes, the genome, where the genetic information is stored in the DNA molecule, shows a hierarchical structure, and this information is integrated into the chromosome of a cell nucleus. Although DNA is a long molecule that can be compacted, forming a highly condensed chromatin structure in the nucleus, the transcription of DNA represents a dynamic process. Studies showed that DNA represents a part of an ordered, folded structure in the cell nucleus, and that the formation of this structure is most likely tightly regulated. Chromatin fibers are formed when DNA molecule is wrapped around the histones, and transcriptionally inactivated and condensed DNA region is known as heterochromatin, while the more transcriptionally active and less condensed region is called euchromatin. Chromatin fibers consist of alternating euchromatin and heterochromatin structures that can interact with each other depending on the alterations in the cellular processes. In the interphase nuclei, the chromatin fibers of each chromosome are highly compartmentalized, and none of the domain structures interact, and this is called a chromosome territory (Chubb et al. 2002;Cremer and Cremer 2010).
The positioning of heterochromatin and euchromatin is related to gene expression. In the nucleus, heterochromatin and euchromatin are spatially segregated, which contributes to the organization of nuclear function. The development of fluorescence imaging and electron microscopy revealed the spatial segregation of chromatin types into distinct subnuclear compartments , and that heterochromatic clusters are not randomly distributed, but are enriched at the nuclear periphery and around the nucleoli . This is called conventional architecture, and it represents a nearly universal nuclear structure, found in the majority of eukaryotic cells. In contrast to this, certain nuclei exhibit an inverted architecture, where heterochromatin is located at the center of the nucleus and euchromatin is enriched at the periphery ( Fig. 1a) (Solovei et al. 2009). Solovei et al. (2009) demonstrated that these different types of nuclear architecture are associated with different mammalian lifestyles (e.g., diurnal versus nocturnal), and are determined by the epigenetic changes. For example, the nuclear architecture of rod photoreceptor cells in the retina of diurnal mammals is typically conventional, while that in nocturnal animals is typically inverted (Solovei et al. 2013). The inverted nuclear architecture of the mouse retina rod cells is determined by the transformation of the conventional architecture, as shown in Fig. 1b. This process is accompanied by the Fig. 1 The organization of the rod cell nucleus during postnatal development (P0-P28) and in a 9-monthold mouse (9 m). a Left panel FISH, using probes for L1-rich heterochromatin (red), euchromatin (green), and chromocenters (blue). Right panel euchromatin distribution of histone H3K4me3; histone modifications (green), nucleoli (blue, arrowheads), and nuclear counterstaining (red). b Reorganization of the nuclear architecture. The nuclei of mouse retinas dissected at birth (P0), 6, 14, 21, and 28 days after birth (P6, P14, P21, P28, respectively), and at 9 months (9 m). When a mouse is born, the rod cell has conventional architecture. However, a few days later, the heterochromatin domains close to nuclear envelope move and fuse with neighboring heterochromatin domains, which consequently leads to the inverted architecture with a single cluster of heterochromatin domain. c Distribution of chromosome subregions in the nuclei with the conventional (left) and inverted (right) architecture: euchromatin (green, white), heterochromatin (red, gray), chromocenters (blue), nucleolus (yellow). The chromosomes in the nucleus with the conventional architecture relocate during the reorganization and form the inverted architecture. The diagrams in a-c) and the descriptions are adapted from Solovei et al. (2009) (colour figure online) relocation of chromosomes from their position observed in the conventional nuclear organization, like slices of pizza, and the creation of a single heterochromatin cluster (hetero-cluster) in the inverted nuclear architecture, as shown by 2D imaging (Fig. 1c). At the time of birth, the nuclei of the rod cells in mice exhibit conventional nuclear architecture. However, the distribution of chromosomes and heterochromatic clusters changes slowly, and the inverted architecture is formed during terminal differentiation. The relationships between the different types of nuclear architecture and nuclear functions are not clear. However, a previous study suggested that different types of nuclear architecture may result in the different rates of light collection efficiency, and that the inverted architecture is more suitable for this process, in comparison with the conventional nuclear architecture (Solovei et al. 2009).
Detailed analyses of this reorganization process showed that the conventional architecture is reorganized through the transformation of the nuclear shape from elliptical to circular, which is accompanied by a decrease in nuclear volume by approximately 40 % (Fig. 1b) (Solovei et al. 2009). Furthermore, structural differences between the conventional and inverted architectures can be attributed to the activity of the nuclear envelope proteins, lamin B receptor (LBR) and lamin A (Lmna), which sequentially tether peripheral heterochromatin (Solovei et al. 2013). The conventional architec-ture is associated with LBR or lamin A/C expression, while the expression of these molecules was not found in the cells with the inverted architecture. However, specific mechanisms mediating these events are still poorly understood, and the factors involved in the association between nuclear architecture changes and other events, such as the alterations in nuclear size or shape, are unknown. It remains unclear how the same inverted structure is consistently created in the nuclei of rod cells. Some theoretical studies have described the spatial organization of nuclear chromatin (Finan et al. 2011;Heermann 2011;Awazu 2014;Ganai et al. 2014). The approach used in these studies is based on a model in which chromatin is represented as loops or strings, and the heterochromatin and euchromatin states are described by the differences in entropic forces or mobility at the protein structure level. These studies showed that chromatin fibers can create chromosome territories, and longrange structures composed of heterochromatin and euchromatin domains within the nucleus. However, these studies were not able to explain the specific mechanisms responsible for the formation of the conventional and inverted types of architecture. String models (Finan et al. 2011;Awazu 2014;Ganai et al. 2014) have shown that the heterochromatic domain is distributed at the nuclear periphery and not at the center of the nucleus, when no relationship between the nuclear envelope and heterochromatin is assumed, which is in contrast to the previously observed inverted architecture (Solovei et al. 2009(Solovei et al. , 2013. This indicates that the string model may not be appropriate for the determination of chromatin structure details using macroscopic models. Further macroscopic descriptions should be included in the model in order to integrate the long-range dynamics of chromatin and its molecular characteristics, and facilitate the elucidation of the key mechanisms involved in the reorganization process. Therefore, we chose an approach that involves a higher macroscopic description by capturing chromatin as a domain, and subsequently describes the overall dynamics of chromatin based on the variations in the domain structures, since the previously described string models have shown that chromatin fibers occupy discrete territories, and that the movement of each chromatin fiber is confined within a domain. The phasefield method has been applied to a wide range of problems related to the complex dynamics of the domain interface, especially in the materials science (Provatas and Elder 2010;Takagi and Yamanaka 2012). Recently, it was used for the numerical simulation of vesicles and their bio-mechanical properties (Maitre et al. 2009;Wang and Du 2008), and was applied in cell dynamics investigations, because it simulates the complex domains of higher dimensions, such as cell shape (Akiyama et al. 2010;Shao et al. 2010). A novel method, using multi-phase-fields, applied to the studies of cell and tissue dynamics, has been proposed as well, leading to the investigations of cell division, cell adhesion, and cell sorting in higher dimensions (Nonomura 2012). Here, we chose to use the multi-phase-field method proposed by Nonomura (2012), and we developed a novel application of this method for the analysis of chromatin dynamics. Capturing chromatin within a domain can yield the information about the distribution of heterochromatin domains, as shown by the image analysis data (Solovei et al. 2009). This allowed us to focus on the elucidation of the mechanisms underlying the reorganization process in nocturnal mammals.
Two contrasting models have been used to describe chromosome territories, in which the chromatin in different chromosomes is either separated by an interchromatin compartment, or not, and in the latter model, it is able to expand into the surrounding territories (Cremer and Cremer 2010). Recently, this latter model has been supported by the findings that chromatin represents a dynamic structure, which can diffuse (Chubb et al. 2002;Gasser 2002). Chromosomes have also been found to intermingle in the interphase nuclei of human cells, and the intermingling pattern was shown to be altered depending on transcriptional activity level and chromosome condensation (Branco and Pombo 2006). Although the intermingling phase is crucial for nuclear functions, the mechanisms underlying the effects of intermingling in chromosome territories on the chromatin structures are unknown. In this study, we suggest a novel role of the intermingling phase, influencing the spatial structure of long-range chromatin distributions.
Here, we first formulate a mathematical model using the phase-field method, and show that the two types of nuclear architecture can be successfully recreated. Following this, we demonstrate how nuclear structures can be dynamically altered depending on the physical features. Afterward, we reveal that the long-range distribution of euchromatin and heterochromatin can be influenced by the size and shape of the nucleus, and finally, we discuss how the degree of intermingling between chromosome territories or between heterochromatin and euchromatin domains (He-Eu domains) plays an important role in the determination of the nuclear architecture.

Model formulation
Chromatin consists of alternating euchromatin and heterochromatin regions, and it is folded into chromosome territories, in which heterochromatin and euchromatin are separated. We first simplified chromatin features, using a domain composed of euchromatin and heterochromatin regions (Fig. 2a) and applied the multi-phase-field method (Nonomura 2012). We defined the three domains as follows: nucleus, chromosomes, and heterochromatin, using a continuous function called the phase-field, which describes two stable states, the present state (1) and the absent state (0) (Fig. 2b). We defined heterochromatin without distinguishing between L1-rich heterochromatin and chromocenters, in order to simplify the model further. In the definition of the three domains, the complement of the heterochromatin region in each chromosome domain represents the euchromatin domain. Therefore, an increase in the heterochromatin implies a decrease in the euchromatin in each chromosome domain. Note that we formulated a mathematical model for two dimensions of eight chromosome territories in order to compare our results directly with the results obtained previously by image analysis (Solovei et al. 2009).

Phase-field model of nuclear architecture
We designed the model using energy functions defined by the nucleus, φ 0 (x, t), chromosome territories, φ m (x, t) (1 ≤ m ≤ N ), and heterochromatin, ψ(x, t), where x ∈ in R n and t > 0. For the simulations, is given by L x × L y where L x and L y are horizontal and perpendicular lengths, respectively (black square shown in Fig. 2b). N represents the total number of chromosomes in the nucleus and N = 8.  The concept underlying the phase-field method. Blue region represents the nuclear domain and is described by the phase-field φ 0 . Green regions are depicted by the phase-fields φ 1 · · · φ N , and they represent chromosome territories, while red regions, indicated by phase-field ψ, represent heterochromatin. The thickness of the intermingling region between chromosome territories is defined by δ φ , and the thickness of the intermingling region between heterochromatin and euchromatin is defined by δ φ (colour figure online) First, we defined the basal free energy functions for chromosome territories and heterochromatin, according to the following equation: where φ , ψ > 0 are gradient coefficients. Note that we defined the phase-field of the nucleus, φ 0 directly, using a function satisfying φ 0 = 0 and φ 0 = 1 for the interior and the exterior regions of the nucleus (Fig. 2b), with a sufficiently small thickness of the interface. Next, we defined spatial range restriction for the chromatin domains, in order to avoid biologically unfeasible events. We made three assumptions based on previous observations and results: (S 1 ) Each chromosome is restricted in its nuclear domain (Alberts et al. 2002), (S 2 ) No heterochromatin region can escape from the given chromosome domain (Alberts et al. 2002), (S 3 ) The chromosome domains are preferentially separated from each other and constitute separate territories (Cremer and Cremer 2010).
The assumptions, S 1 -S 3 are described as where β 0 , β ψ , and β φ are positive constants and indicate the intensities of domain for more details about h). E 0 and E 1 are fundamental formulations describing nuclear and chromatin domains.
In the third step, we defined chromosome and heterochromatin volumes, so that they can reflect the changes in chromatin territories. Three assumptions were made: (R 1 ) The nuclear space is fully occupied by chromosomes, (R 2 ) The chromosome can be contracted or expanded to a given volume, (R 3 ) Heterochromatin is converted from/to euchromatin within each chromosome. Hara et al. (2013) showed that the chromosome condensation is affected by the number of chromosomes per nuclear size, and the reduction of nuclear size leads to the reduction of the size of the condensed chromosome. Therefore, we made the assumption R 1 , and that the total volume of chromosomes is related to the nuclear volume. Chromosomes are condensed before a cell division and expand after cell division. However, it is unclear how the volume of the chromosome is regulated, and we simply assumed R 2 , based on the observations. R 3 is assumed based on the changes of the DNA molecule transcriptional states during the differentiation. R 1 -R 3 are formulated using the equation where α 0 , α V , α v > 0 are the energy intensity constants of each volume. In the first term of E 2 , [1 − h(φ 0 )]dx corresponds to the nuclear volume, because h(φ 0 = 0) = 0 and h(φ 0 = 1) = 1. R 1 indicates that the total volume of chromosomes has a minimal energy when its equals nuclear volume. This leads to the full occupancy of the nucleus. In the R 2 and R 3 equations, is the volume of the heterochromatin in an m-th chromosome at time t, and they are given by are the functions of the target volumes of condensed chromosomes and heterochromatin, or those that extend to and have a minimal energy at the target volume.
In the final step, we include LBR and lamin A/C, which have a role in the interactions between heterochromatin and the nuclear envelope (Solovei et al. 2013), in our model. We expressed this through the intensity of affinity between nuclear function, φ 0 , and heterochromatin region function, ψ, as follows: where γ is a constant that determines the intensity of the affinity. Note that γ > 0 reflects the preference of heterochromatin to remain at the nuclear periphery, which is interpreted as the expression of LBR and lamin A/C in the nuclear envelope and their tethering of heterochromatin to the nuclear periphery. In contrast to this, γ = 0 represents the absence of LBR and lamin A/C , which leads to the lack of heterochromatin and nuclear envelope interactions. The total energy of chromatin dynamics is given as (1) Afterward, we determined the functional derivatives of Eq. (1) with respect to φ m (1 ≤ m ≤ N ) and ψ, which drives the time for the system to evolve, satisfying where μ is a positive constant that represents the mobility of each phase-field. The calculations of the previous equations generated the reaction-diffusion model as follows: where A m and B are given as the functions of , t)). The reaction-diffusion system (2) describes the interface of chromosome territories and heterochromatin domains that changes their dynamics, depending on A m and B, respectively. We non-dimensionalized the model (2) for time scale, T , and spatial scale, L, incorporating t = Tt and x = L Q x into the model (2). The mobility constant, μ, is defined by and set We then obtained the following system by removing the tilde: In the model, we assumed that the size and shape of the nucleus change independently of chromatin states, and the phase-field φ 0 is set as an independent function, describing the nucleus independently of the other phase-field functions.

Conventional architecture
The conventional architecture is the primary structure in the majority of eukaryotic cells and it is evolutionary conserved from unicellular to multicellular organisms. Because the conventional architecture is formed by the chromosome uncoiling during the initial stage immediately after cell division, we assumed that the target volumes of chromosomes constant and that the rate of heterochromatin conversion at each chromosome is not likely to change significantly. Therefore, we setV

Inverted architecture
The inverted architecture is formed by the reorganization of conventional architecture during terminal differentiation, accompanied by a decrease in nuclear and chromosome volumes (Solovei et al. 2009). Therefore, the target volume of chromosome was set so thatV m (t) = rV m (r ∈ (0, 1)) in E 2 . In contrast to this, nuclear transcriptional activity is likely to be lower after the final cell division and during terminal differentiation, and therefore, the target volume of the m-th heterochromatin at time t in the m-th chromosome was set asv where ρ m (t) represents the rate of change from euchromatin to heterochromatin at the m-th chromosome. We assumed that the heterochromatin conversion rate increases, which is reflected by ρ m (t), upon choosing an increasing monotone function. We used a sigmoid function described by where ρ m (0) is v m (0)/V m (0), and α 1 , α 2 and t * are positive constants,ρ m is the increasing rate of v m /V m , and ρ m (0) +ρ m attains the maximal rate of v m /V m . The model for the inverted architecture is given by system (3) with

Extent of the intermingling of territories
The phase-field approach assumes that the dynamics of the interface connecting the two states determines the dynamics of chromosome territories and heterochromatin. This implies that the thickness of the interface can be reinterpreted as the intermingling between chromosome territories or the intermingling between heterochromatin and euchromatin domains (He-Eu domains). We obtained the measure of intermingling directly by calculating the thickness of the interface. A novel translation for the phase-fields φ 1 , . . . , φ m , ψ in Eq.
(3) is defined by translating the phase-fields of chromatin as the relative DNA content (%). Therefore, φ m = 1 and ψ = 1 indicate the maximum density state of euchromatin and heterochromatin at the m-th chromosome, meaning, the highly condensed state. If we assume that the conditions S 2 and S 3 in E 1 are not too strong or too weak, we find that the phase-fields overlap in the interface region, which satisfies 0 < φ 1 , . . . , φ m , ψ < 1 for fixed β ψ and β φ . That is, the scale of intermingling of chromosome territories or He-Eu domains is determined by the interface thickness. Following this, the extent of intermingling of each phase-field, δ φ and δ ψ , can be explicitly calculated (Takagi and Yamanaka 2012), as follows: where λ is a constant that defines the interface region, so that λ ≤ φ 1 , . . . , φ m , ψ ≤ 1 − λ. Here, we mostly chose the value λ = 0.15. The scale of φ and ψ defines the scale of intermingling in our model.

Parameters and numerical methods
The estimation of all parameters of the model from experimental data is difficult, because the reports describing chromatin dynamics are not numerous. We typically performed simulations with a non-dimensionalized system and then verified representative parameters by estimating dimensional scales through the comparisons of qualitative dynamics and previously obtained temporal data about the reorganization process (Solovei et al. 2009). In system (3), we defined the mobility of phase-fields as μ = T −1 . That is, we can determine μ by estimating the time scale, T . For this, we directly compared qualitative dynamics of a non-dimensionalized system (3) with the previously reported experimental data (Solovei et al. 2009) (Fig. 1). We were then able to estimate T from t = Tt, which consequently allowed the determination of μ. With representative parameters set, T was estimated to be 5 h, and we obtained μ = 1/5 (h −1 ).
Additionally, we were able to estimate the spatial scale directly from the experimental images obtained previously by Solovei et al. (2009), in which the size of the nucleus of the rod cell in two dimensions was determined to be approximately 4-5 µm (x-axis) and 6-8 µm (y-axis) in P0 cells. We chose 5 µm as the x-axis diameter and and 8 µm as the y-axis diameter of the elliptic nucleus, and used L x = 6 µ m for the x-axis and L y = 9 µm for the y-axis in space for numerical simulations. In the non-dimensional system, we used 1.2 × 1.8 square, so that L = 5 µm allows L x = 6 µm for the x-axis and L y = 9 µm in the dimensional system.
With T , μ and L known, we are able to confirm whether ε 2 φ = 2 φ L −2 and ε 2 ψ = 2 ψ L −2 fall within a reasonable range of parameters in the dimensional system (3). First, we scaled the dimensionless parameters that we used for simulations with T , μ and L, and obtained the dimensional values shown in Table 1. μ 2 φ and μ 2 ψ from Eq. (3) are considered the diffusion coefficients of euchromatin and heterochromatin. Chromatin mobility in mammalian nuclei is reported not to exceed 0.4 µm for time periods of over 1 h (Abney et al. 1997), and therefore, the diffusion coefficients of chromatins can be estimated on a scale of less than 0.16 µm 2 /h. The parameter values in Table 1 show that the diffusion rates of chromatin we used are μ 2 φ ∈ [0.9 × 10 −3 , 3.7 × 10 −3 µm 2 /h] for euchromatin and μ 2 ψ ∈ [0.98 × 10 −3 , 3.92 × 10 −3 µm 2 /h] for heterochromatin, since both of these are sufficiently smaller than 0.16 µm 2 /h.  (7). This scale is also considered to be in a feasible range, because the length of chromatin fibers is thought to be approximately 300 nm in eukaryotic genomes (Rosa and Everaers 2008).
The representative dimensional/nondimensional parameters are presented in Table  1, except for α 0 , α V , α v , β 0 , β φ , β ψ , which are the same in both non-dimensional and dimensional systems, so that α 0 = 25/6, α V = 10/6, α v = 20/3, β 0 = 5/3, β φ = 1, and β ψ = 2/3. We solved the Eq. (3) using an explicit numerical algorithm written in C programing language. In order to solve the Laplacian terms in the model, we used the standard method by applying twice the central difference operator (Morton and Mayers 1994). Therefore, the numerical simulation does not blow up if (diffusion constant)×[(timestep)/(gridsize of x-axis) 2 + (timestep)/(gridsize of y-axis) 2 ]≤ 1/2, with sufficiently small time steps. We took a time step of 6 × 10 −4 for the grid size 6 × 10 −3 of x-and y-axis. The ratio of minimal interface width to the chosen grid size was approximately 26.33, and the grid size was sufficiently small not to influence the movement of interfaces.

Simulation results
When presenting the simulation results, nucleus was depicted in blue, chromosome territories in green, and heterochromatin in red (Fig. 2b). The complementary region of heterochromatin in each green domain is implicitly euchromatin. We first generated the conventional architecture and then show the inverted architecture through the regeneration of the reorganization process. Following this, we explored the mechanisms underlying the reorganization process and the pattern formation of chromatin.

Conventional architecture
The highly condensed chromosome in the early stage of mitosis or meiosis becomes uncoiled after cell division, and the chromatin in the mouse rod cells has the conventional architecture at the time of birth. In order to regenerate the conventional architecture, we initiated simulations under the initial condition, as shown in Fig. 3, and compared the results obtained by changing the parameters that control the intensity of domain territories and the affinity between heterochromatin and the nuclear envelope. A high intensity of domain territory indicates that chromosome territories are strongly confined. Zero affinity indicates that there are no interactions between heterochromatin and the nuclear envelope, but a positive affinity value represents a tethering of heterochromatin to LBR or lamin A/C on the nuclear envelope.
The results are shown in Fig. 3a, b. We first demonstrate that the conventional architecture is obtained with the positive affinity. In contrast to this, when there is no affinity, heterochromatin domains are confined completely to the center of each chromosome if the domain separation is strong, or they are fused with adjacent heterochromatin if the domain separation is weak. The obtained results indicate that positive affinity is a required condition for the localization of heterochromatin at the nuclear envelope. In particular, 6 days after the cell division, the difference in territory intensities results in slightly different types of conventional architecture. In Fig. 3a, it is shown that heterochromatin accumulates at the territories between chromosomes instead of in the region of the nuclear envelope. Therefore, small domes are generated between A Strong territory intensity B Weak territory intensity  Fig. 3b, however, heterochromatin is shown to be distributed almost homogeneously along the nuclear envelope.
The division of rod cells in the mouse ceases 5 days after birth (Solovei et al. 2009), and, the conventional architecture of P0 in Fig. 1 is likely last less than 5 days. We chose to use data from day 3 as the initial condition for inverted architecture simulations, since it is not sensitive to the parameter of territory intensity. Our simulation results indicate that the affinity between heterochromatin and the nuclear envelope is crucial for the formation of the conventional architecture, which confirms that the expression of LBR and lamin A/C is indispensable for the generation of conventional nuclear organization, as previously shown (Solovei et al. 2013).

Mechanism of single hetero-cluster inverted architecture generation
The successful regeneration of the reorganization is based on the following five assumptions: (1) The size of the nucleus decreases.
(2) The shape of the nucleus changes from elliptical to circular.
(3) The rate of conversion of heterochromatin to euchromatin increases. (4) There is a lack of affinity between the nuclear envelope and heterochromatin. (5) The nucleus is fully occupied by chromosomes.
We evaluated these five conditions sequentially, in order to define which conditions are indispensable for the induction of the single hetero-cluster inverted architecture generation. First, we fixed the size of the nucleus, in order to determine whether a decreased nuclear size and circular shape are indispensable. Surprisingly, the obtained data, presented in Fig. 5a indicated that the increase in the heterochromatin conversion rate is sufficient for the formation of the inverted architecture. This demonstrated that the change to a circular nuclear shape is not crucial for this process. However, when the increase in the heterochromatin conversion rate is insufficient, the inverted

Fixed Nucleus Decreasing Nucleus
Circle Ellipse

Fig. 5
Mechanism leading to the generation of a single hetero-cluster inverted architecture. The patterns were generated from P0. a Nuclear size and shape are fixed under conditions (3)-(5), defined in the text. The parameters used are the same as in Fig. 4, except for ρ m (t) = 0. b The heterochromatin conversion rate is fixed for a fixed-sized nucleus, and circular and elliptical nuclei that decrease in size (80 % cut-off) under conditions (4)-(5), defined in the text. The parameters used are the same as in Fig. 4, except for (ρ 1 ,ρ 2 ,ρ 3 ,ρ 4 ,ρ 5 ,ρ 6 ,ρ 7 ,ρ 8 ) = (0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0). c The affinity between the nuclear envelope and heterochromatin is present for a fixed-sized nucleus and circular and elliptical nuclei that decrease in size (80 % cut-off) under conditions (3) and (5), defined in the text. The parameters used are the same as in Fig. 4, except for γ = 0.022. (d) The unoccupied space in the nucleus is considered a fixed nucleus and circular and elliptical nuclei that decrease in size (40 % cut-off) under conditions (3) and (4), defined in the text. The final volumes of euchromatin and heterochromatin decrease to approximately 49 % compared with P0 volumes, in each panel. The parameters used are the same as in Fig. 4 architecture formation, with a single heterochromatin cluster, cannot be achieved, even when the decrease in the size of the nucleus is sufficiently small, regardless of its shape (Fig. 5b). When the heterochromatin conversion rate is constant, it is possible for the heterochromatin domain to fuse, so that the number of clusters decreases. However, a single cluster of heterochromatin cannot form without a sufficient increase in heterochromatin conversion rate. Next, we assessed the effect of the affinity between the nuclear envelope and heterochromatin. The data presented in Fig. 5b confirm that the absence of affinity is indispensable for the formation of inverted architecture, regardless of nuclear size and shape. When the affinity between the nuclear envelope and heterochromatin exists, it remains at the nuclear periphery and the inverted architecture is never formed. Figure 5d shows that the inter-chromatin compartment and the space between the chromosome and the nuclear envelope must be completely occupied, in order to induce the formation of a single heterochromatin cluster. When the total volume occupied by chromosomes is insufficient to fill the nucleus, chromosome territories are separated more strictly and heterochromatin is unable to fuse, regardless of the nuclear shape. Therefore, chromosome volumes must decrease to a size no smaller than that of the nuclear volume, in order for the single hetero-cluster inverted architecture to form.
A decrease in nuclear size and nuclear deformation are not indispensable for the reorganization process, but the increase in the rate of conversion of heterochromatin to euchromatin, the absence of affinity, and the absence of unoccupied space in the nucleus are crucial for the induction of the transformation from the conventional to the inverted architecture.

Nuclear size and deformation effects on nuclear patterns
In the previous section, we have found that the deformation of nuclear size and shape were not indispensable for the formation of the single hetero-cluster inverted architecture. However, it was reported that the size of the rod cell nuclei differs by approximately 40 % between P0 and 9 m and the nuclear shape changes from an elliptical to a circular (Solovei et al. 2009). Therefore, we investigated how these two features influence chromatin dynamics during the reorganization.
First, we found that the nuclear size and shape influence the time scale of reorganization ( Fig. 6A1 and B1). The formation of the same type of inverted patterns in the fixed-domain case required almost twice as long to complete, compared with the formation in a 40 % smaller nucleus. Similarly, the formation in a circular nucleus lasts longer than in an elliptical nucleus, because the distance that two central chromosomes need to cover in order to reach the nuclear periphery is longer in larger and circular nuclei compared with that in the smaller and elliptical nuclei. The speed of the movement of the interface between chromatin territories is similar.
The nuclear shape significantly influences chromosomal rearrangement (Fig. 6B3). Although the two geometries minimally alter the difference between heterochromatin patterns, the position of chromosomes in an elliptical nucleus changes more dynamically compared with that in a circular nucleus.

Intermingling of chromosome and He-Eu domains during the formation of a single hetero-cluster nuclear architecture
In order to identify the role of intermingling on the reorganization of nuclear architecture, we focused on the extent of the intermingling between chromosome territories (δ φ ) and the extent of intermingling between a euchromatin and heterochromatin domains (δ ψ ), which were directly calculated in our model. We investigated how the inverted architecture is changed by the extent of intermingling (Fig. 7). First, the extent of intermingling exerts a significant effect on the terminal pattern of inverted architecture, and the two regions of intermingling, δ φ and δ ψ , play different roles. In the case of He-Eu domains (Fig. 7a), the larger intermingling region is observed prior to the generation of a single heterochromatin cluster. In contrast to this, the intermingling between chromosome territories shows the opposite effects (Fig. 7c), i.e., the smaller intermingling region is observed prior to generation of a single heterochromatin cluster. These observations indicate that there is a mini-
Next, we focused on the relationship between the extent of intermingling and the nuclear size required for the formation of the inverted architecture (Fig. 7b, d). We identified a minimum thickness of He-Eu domains necessary for the formation of the single hetero-cluster inverted architecture, which is presented in Fig. 7b and a maximum thickness of chromosome territories necessary for the formation of the single hetero-cluster inverted architecture, presented in Fig. 7d. The obtained data suggest that the extent of intermingling of He-Eu domains needs to be expanded for the formation of the single hetero-cluster inverted architecture, when the nucleus is large. In contrast to this, the extent of intermingling of chromosome territories needs to be reduced when nucleus is large, and the single hetero-cluster inverted architecture is more likely to form when the chromosomes occupy strictly defined territories.

Discussion
We studied two types of nuclear architecture and the process of reorganization of the photoreceptor rod cell nucleus from the conventional to the formation of inverted architecture, using mathematical modeling, with the phase-field method. Our analyses demonstrate that an increase in the rate of conversion of heterochromatin to euchromatin and the loss of affinity of the nuclear envelope for heterochromatin in the absence of both LBR and lamin A/C expression are indispensable for this reorganization. Furthermore, the extent of intermingling between chromosome territories strongly influences the formation of the single hetero-cluster inverted architecture, related to a specific nuclear size. These findings suggest that the transformation of the chromatin from euchromatin to heterochromatin state induces initially the long-range movement of the chromosome territories. Additionally, the physical features of chromosome territories or nucleus play a crucial role in the determination of chromatin dynamics.
It was shown previously that nuclear size and shape are associated with mammalian lifestyle (Solovei et al. 2009). In nocturnal mammals, a smaller and circular nucleus with a single hetero-cluster inverted architecture reduces light scattering. Here, we analyzed the direct influence of nuclear size and shape on nuclear architecture. Our results demonstrate that the size and shape of the nucleus are crucial for the determination of inverted architecture pattern, and the temporal scale of the reorganization process. This is the first theoretical study to consider the effect of nuclear deformation on chromatin reorganization.
We explored the effects of chromosome intermingling on nuclear architecture and showed that the intermingling of chromosome territories and He-Eu domains can play opposite roles in the creation of the single hetero-cluster inverted architecture, based on the nuclear size. This suggests that nuclear size can play an important role in chromatin distribution, independently of, or together with, the intermingling dynamics.
Previous studies (Awazu 2014; Finan et al. 2011;Ganai et al. 2014;Heermann 2011), employing the microscopic modeling approach to chromatin dynamics, have shown that a phase separation between heterochromatin and euchromatin exist, which leads to the formation of long-range clusters in heterochromatin. Therefore, the conversion of euchromatin to heterochromatin, and heterochromatin mobility can be investigated by analyzing heterochromatin domain dynamics. Since our model is based on a macroscopic approach, we were able to capture the overall chromatin dynamics in the nucleus, even when the detailed features of chromatin at the fiber level were not known. Many previous studies have successfully applied the macroscopic approach to gain understanding of micro-scale problems. For example, the phase separation problems in polymer solution systems, such as block copolymer dynamics, have been solved by the macroscopic approach, which does not depend on the microscopic details of polymers. This confirms that the macroscopic approach may represent a very useful tool for the facilitation of the understanding of soft matter systems (Fredrickson et al. 2002;Pinna and Zvelindovsky 2012;Yamada et al. 2004).
The details of chromatin structures and how chromatin dynamics is regulated spatially and temporally has not been completely elucidated, and our mathematical model, using a top-down approach, may help understand these structures and processes. The findings presented here suggest a potential mechanism underlying the reorganization of nuclear architecture. Although the properties of chromatin are simplified in our model, and we have made assumptions about several features of the reorganization process, the obtained results suggest that understanding the primary physical features of the nucleus may be sufficient to allow the understanding of the core mechanism of nuclear reorganization. Our study represents the first step toward the understanding of chromatin dynamics by incorporating the information about the molecular properties of chromatin with the macro level chromatin dynamics data. Furthermore, our model demonstrates the applicability of phase-field method in life science investigations. Therefore, W (φ) describes a double-well potential which has local minimums at φ = 0 and 1 if | f s − f l | < 1/12. The constants f s and f l correspond the free energy densities for the phase described by φ = 1 and φ = 0 By taking the functional derivative of E basal with respect to φ, we can induce a time evolution model of φ driven by ∂φ/∂t = −μδ E basal /δφ, as follows: where μ is a positive constant. Eq. (8) provides a smooth front solution connecting the region φ = 1 and φ = 0 so that the direction of the front movement depends on sgn( f s − f l ). Formulating the terms of f s and f l using the function h(φ) is the core of the phase-field modeling. This model gives a traveling wave solution connecting φ = 0 and φ = 1. In fact, the depth of the well of the function W (φ) is controlled by the constants, f s and f l , as shown in Fig. 8. If f l < f s ( f l > f s ), φ = 0 (φ = 1) region expands. In the case of f s = f l , this implies that the phase interface is stopped.