Skip to main content

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

Abstract

Background

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

Results

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

Conclusions

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

Background

Interphase chromosomes are intricately folded and packaged inside the cell nuclei where they interact with nuclear bodies and remodeling factors [1,2,3,4,5]. This complex 3D organization of the genome plays a crucial role in regulation of gene expression, DNA replication, and DNA double-strand break repair [6,7,8,9,10,11,12,13]. Chromosomes are spatially partitioned into sub-megabase topologically associating domains (TADs), which have emerged as the fundamental structural and functional units of chromatin organization in the eukaryotic cell nuclei, from fruit fly to human [1, 2, 14,15,16,17]. Interactions among TADs contribute to higher level, hierarchical genome organization [18]: multiple TADs within each chromatin compartment and multiple compartments within each chromosome territory. Boundaries of TADs are conserved among cell types and sometimes across species [14, 16, 19,20,21,22,23,24,25]; this feature makes TAD a natural structural unit of chromatin useful for modeling purposes and reasoning [26, 27], including in this work.

TADs interact among themselves forming dynamic, relatively short-living 3D-structures, which segregate chromatin into mutually excluded active (type A) and inactive (type B) compartments [18, 28]. Averaged over time of the entire interphase, the effect of TAD–TAD interactions is seen in experimental Hi-C maps course-grained to TAD resolution [1, 14, 26, 29]. However, such Hi-C maps do not reveal how mobile are the 3D positions of each individual TAD. Live-cell imaging studies of interphase nuclei in different organisms showed that some chromatin loci are mobile [30,31,32,33,34,35]. In addition, studies suggested that constrains on the motion of chromatin loci can result from their interaction with structural elements such as nucleolus [36], or nuclear lamina (NL) which lines the inner surface of the nuclear envelope (NE) [37, 38], or the nuclear pore complexes [31, 39,40,41,42]. However, genome-wide analysis of nuclear dynamics of all “freely floating” and “tethered” chromatin loci has not been performed. DamID experiments have identified hundreds of genome regions that can “anchor” to the NL [37, 43,44,45,46]. These genome regions that have specific affinity to the NL are called lamina-associated domains (LADs) [47,48,49,50].

Positioning of LADs in the nuclei of mammalian cells [34, 48, 49, 51, 52] was shown to be stochastic: averaged across cells, a certain % of LADs is bound to the NE in a given cell, but binding of specific LADs varies from cell to cell. Single-cell DamID studies demonstrate [34, 51] that LADs identified in a human cell population may in fact be located either at the nuclear periphery (30% of LADs) or in the nuclear interior in individual cells [34]. LADs that contacted the NL were shown to move during interphase [34], but this movement is confined to a relatively narrow 1 \(\mu\)m layer next to the NE. These observations raise several questions [50]. Is positioning of LADs at the NE only stochastic between cells, i.e., LADs have mostly stable positions within a cell, but highly variable positions between cells, or is it truly dynamic/mobile in each cell, i.e. individual LADs can change their positions relative to the NE significantly within a cell during interphase? If LADs are mobile, on what time-scale? Also, do different LADs have different probabilities of being attached to the lamina?

States of transcription activation and repression are also linked to the positioning of chromatin with respect to the NE. LADs are typically repressive chromatin environment [43, 44, 46, 50]. This observation supports the notion that the nuclear periphery is generally occupied by inactive chromatin [53]. TADs and LADs are similar in size: each is approximately 1 Mb in the human genome [14, 44] and approximately 90–100 kb in the D. melanogaster genome [1, 15, 47]. There are 1169 TADs and 412 LADs in the D. melanogaster genome [1, 47]. The role of LAD–NE interactions and their interplay with compartmentalization of TADs has been explored in Drosophila [35, 54] and mammals [52, 55,56,57,58,59,60]. For example, depletion of lamins from Drosophila S2 cells leads to chromatin compaction and reduction in spatial segregation of the chromatin into active and inactive compartments [35]. Simulations of multiple copies of mouse chromosomes 1 and 2 [55] have disproved the role of lamina as the main driver of compartmentalization [48]. At the same time, attractions between heterochromatic TADs emerged as the main force of compartmentalization, while LAD–NE interactions are crucial for controlling the global spatial morphology of the nucleus [54, 55, 57]. However, a number of questions remains unanswered, related to sensitivity of the general organization principles with respect to global loss of chromatin–lamina interactions, which may occur in disease or senescence [52, 61,62,63,64]. Also, the strength of TAD–TAD interactions can naturally vary during life of organisms. For example, the chromatin compartmentalization is weaker in embryonic cells and stronger in adult cells of Anopheles mosquitoes [65]. Given that attractive TAD–TAD interactions play a major role in compartmentalization, will the global chromatin architecture, such as its radial distribution, become more sensitive or less sensitive to the variations in these interactions upon lamin depletion? Will the cell-to-cell variability in the 3D chromatin organization increase or decrease upon disruption of LAD–NE interactions?

Understanding the principles and factors leading to the formation of non-random 3D genome organization and, ultimately, to the structure–function relationships in chromatin is prerequisite to understanding cell physiology [4, 55, 66]. However, relevant experiments that may provide statistically significant conclusions remain difficult. Computational models that faithfully reproduce available experimental data are indispensable for the 3D genome reconstruction problem (3D-GRP) based on DNA-proximity ligation data such as 5C, Hi-C and Pore-C, and they can generate valuable predictions and guide experiment [26, 52, 55, 56, 67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92]. A particularly strong feature of computational models is that they can answer questions that may be very hard to address experimentally [93]. For example, tracing the movement of a given TAD or a few chromosomal loci upon detachment from the NE caused by a lamin depletion is possible experimentally [35, 94], but making statistically significant statements based on t1000s of TADs in the nucleus would be extremely laborious. Another example is obtaining time-resolved Hi-C maps for a synchronized cell population: possible, but difficult experimentally, with only a handful of studies so far [95,96,97,98,99] Models can generate testable hypotheses, the most promising of which can be checked experimentally.

Since simulating an entire eukaryotic genome on biologically meaningful time-scales at fully atomistic resolution is still far out of reach [100], current practical models accept various levels of approximation determined by the balance of research goals and computational feasibility. A large class of recent computational models, which aim to understanding factors affecting 3D chromatin structure, employ the “beads-on-a-string” coarse-graining approach [52, 60, 68, 77, 82, 83, 90, 101, 102], borrowed from polymer physics models [103, 104]. For complex organisms, such as mammals, various further approximations are often made to reach the desired temporal or/and spatial resolution, such as considering only a small subset of the dozens of the original chromosomes [55] or even only one chromosome [52, 56, 59]. In that respect, modeling “simpler” nuclei [71, 90], in particular of higher eukaryotes, e.g., that of the well-studied fruit fly, with its only five major chromosome arms, offers a computational advantage that may translate into an ability to model the entire nucleus [26, 77, 83], which, in turn, may help to answer questions otherwise difficult to address. Consensus conclusions about chromatin organization that emerge from using models of different types applied to substantially different organisms are of value, as these conclusions hint at conservation of general principles across species.

Our work contributes to developing a qualitative and quantitative understanding of various aspects of 3D structure and dynamics of chromatin in D. melanogaster nuclei by constructing and employing a computational model to make various testable predictions. The model is trained to reproduce contact probabilities between TADs from the experimental Hi-C data [1, 26] and the average fraction of LADs near the NE [43], and is validated on multiple structural features of chromatin from several different experiments. It allows us to make multiple, biologically relevant predictions followed by a discussion about their possible biological significance. We have used the model to investigate sensitivity of the spatial organization of the interphase chromatin to the variation of the interactions among chromatin domains, clarifying the role of LAD–NE interactions in stabilization of the local and global chromatin structure. The model makes a genome-wide prediction of a highly dynamic nature of LAD positioning in fruit fly interphase nuclei—most of the LADs within a single nucleus attach to and detach from the NE multiple times, moving far away from the nuclear periphery on the time scale of the interphase. Previous experimental study [35] of LAD mobility in fruit fly was limited to only three LADs. Our model also predicts that the very different probabilities of individual LADs to be in contact with the NE are determined by the highly variable local linear density of LADs along the chromosomes which also affects the radial positioning of individual TAD, determining non-random average global structure of chromatin.

The dynamical model of the 3D chromatin organization in fruit fly proposed in this study offers a number of novel genome-wide insights, both biological and methodological, while also reinforcing robustness of some of the conclusions about chromatin organization previously made using models of mammalian nuclei. Combined with appropriate experimental data, our model can be used to make structure–function predictions on genome-wide level [105].

Methods

Background and rationale for choosing model building blocks at TAD level

Our model is coarse-grained at the resolution of individual TADs. Numerous chromatin contacts within each TAD make the intra-TAD interactions effectively stronger than the inter-TAD ones. This difference allows one to consider the intra-TAD interactions as the major factors responsible for the shapes of TADs, while neglecting the effects of the inter-TAD and TAD–NE interaction on these shapes. We use these considerations and omit the details of the internal structure of TADs, considering them as the smallest building blocks of our model and using the shape-independent inter-TAD and TAD–NE interactions.

In Drosophila melanogaster nuclei each TAD compartmentalizes on average \(\sim\)100 kb of chromatin [1]; there are four major epigenetic classes of TADs: Active, Null, PcG and HP1 [1, 106] (see Additional file 1: Fig. S1). Active TADs are defined by histone H3 modifications such as trimethylation of lysine 4 and 36 (H3K4me3 and H3K36me3). PcG TADs are enriched in Polycomb group proteins and histone mark H3K27me3. HP1 TADs are associated with classical heterochromatin marks such as H3K9me2 histone modification, heterochromatin proteins HP1 and Su(var)3-9. Null TADs lack known specific chromatin marks [1]. The differences in TAD–TAD interactions within each of these classes and between different classes are major factors responsible for the segregation of chromatin into active (type A) and inactive (type B) compartments [5, 55].

The four classes of TADs in fruit fly further refine the distinction between type A and type B compartments. Type A compartments are organized by Active class TADs and defined by early replication, proteins, and histone modifications involved in active transcription. Type B compartments are defined by late replication and modifications that silence genes. These compartments consist of transcriptionally silent TADs (PcG, HP1 and Null) and occupy a larger portion of the genome than type A compartments [1, 15, 107]

Even though our model is coarse-grained at TAD resolution and omits details of the internal structure of TADs, it provides a rather realistic description of the nuclear architecture by taking into account different mutual arrangements of chromosome arms (nucleus topologies) [108], different epigenetic classes of TADs [1, 106] and their interactions, and a proper distribution of LADs along chromosome chains [47].

Model: chromosome representation

We consider D. melanogaster female interphase nuclei with a diploid set of four chromosomes (2, 3, 4 and X). Hi-C data suggest [1] that D. melanogaster genome (excluding pericentromeric constitutive heterochromatin (HET) and centromeric (CEN) regions) is organized into 1169 TADs. Since the same TADs in homologous chromosomes are almost always in close proximity to each other [26, 30, 109,110,111], we represent each pair of homologous TADs by a single spherical bead. In addition to these 1169 TAD-beads, we introduce 4 beads representing CEN regions of the four chromosomes and 6 beads representing HET domains.

Using the “beads-on-a-string” model [103, 112], four chains of the beads represent four paired homologous chromosomes (Chr). In Chr 2 and 3, the L- and R-arms are connected via three-bead HET–CEN–HET structure. Chr 4 begins from a CEN–HET bead pair, and Chr X ends with an HET–CEN two-bead structure. The Chr chains are placed inside of a spherical boundary (see Fig. 1, top panel) which represents the nuclear envelope (NE). We also consider a nucleolus which is modeled as a spherical bead of 0.333 \(\mu\)m radius placed at half distance between the nucleus center and the NE [26].

Bead size and mass

The mass \(m_i\) of each bead, which represents a pair of homologous TADs, is related to the DNA sequence length \(L_i\) (in bp) in the corresponding TAD and consists of the DNA mass (\(\sim\)660 Da/bp) and the mass of proteins associated with a nucleosome (132500 Da/nucleosome) [113, 114]. Using a 200 bp value for the nucleosome repeat length [115,116,117] one can get the following: \(m_i = 2 L_i (132500/200 + 660)\).

We assume that the volume of each bead is proportional to the DNA length in the corresponding TAD, and we employ the modified TAD radii determined in [26] (“hard radii”), scaled by a factor 1.254031 to reflect the double volume (two homologous TADs) of our beads. The average size (diameter) of our beads representing TADs is 0.18 \(\mu\)m, ranging from 0.08 to 0.38 \(\mu\)m. Additional details can be found in Additional file 1.

Bead–bead interactions and bead types

Since our model is coarse-grained at TAD resolution, the interactions between beads reflect the interactions between different TADs, and depend only on the epigentetic class of TADs. These interactions are tuned to reproduce the experimental Hi-C map coarse-grained to TAD resolution [26]. See a more detailed description Additional file 1, and in the “Model development” section.

Interactions in 268 specific pairs of remote loci

See detailed description in the Additional file 1.

Arrangement of chromosome arms (nucleus topology) and preparation of their initial configurations

Fig. 1
figure 1

Top: A schematic showing the key elements of the model nucleus of D. melanogaster. Bottom panels: The four initial configurations (bottom panels) of the D. melanogaster chromosome arms [2 L (red), 2R (orange), 3 L (yellow), 3R (green), X (blue) and 4 (cyan)] that serve as the starting points of the simulations described here. The nucleolus is shown as a gray sphere (on bottom panels). The arms are fully “territorial” at the beginning of each simulation. The initial configurations correspond to the different mutual arrangements of the chromosome arms (nucleus topologies) experimentally determined by Hochstrasser et al. [108]. From left to right: CIS-X6S, CIS-X7N, TRANS-X3S and TRANS-X4N nucleus topologies. Here, the “CIS” configuration refers to chromosome arrangement in which the two “R” arms (orange and green) or two “L” arms (red and yellow) of autosomes are next to each other in 3D space, while “TRANS” arrangement is complementary to it

Several mutual arrangements of the chromosome chains that correspond to CIS and TRANS relative arrangements of the Chr 2 and 3 arms, and two most probable positions (“North”—“N” and “South”—“S”) of the Chr X have been observed experimentally [108], Fig. 1, bottom panels. The CIS arrangement refers to the one in which the R-arms (or the L-arms) of the autosomes are next to each other. The TRANS arrangement refers to the positions of the R- and L-arms of the autosomes next to each other. In the simulations, we use the following arrangements of the nucleus topology [108]: CIS-X6S, CIS-X7N, TRANS-X3S and TRANS-X4N. The CIS:TRANS arrangements ratio is assigned its experimental weight of 2, meaning that twice more replicas of CIS topologies are simulated. In addition, each of these 6 properly weighted arrangements has 3 replicas to account for 3 different nucleus sizes used. Thus, the entire simulated ensemble consists of 18 model nuclei.

To create initial configurations of the bead chains that represent Chr X, L- and R-arms of Chr 2 and 3, we generated them as separate linear strings of beads parallel to z-axis (the axis of centromere-to-telomere chromosome polarization). Each chain is initially separately restrained by two planes intersecting along the z-axis with the dihedral angle 72\(^{\circ }\). To collapse initially linear chain of beads into the nucleus sphere, a harmonic interaction between the nucleus center and each bead is applied, see Additional file 1: Eq. S1, with the effective equilibrium length 0.99 \(\mu\)m and harmonic spring parameter \(\kappa =0.5\) kT. A subsequent Langevin dynamics of the restrained chains over 10\(^5\) 1.36 ns time steps produced collapsed fractal-like Chr-arm configurations [112]. This harmonic interaction is removed from all subsequent simulation stages. The L- and R-arms of Chrs 2 and 3 are then connected through the CEN beads, Chr 4 is added at the “North” pole, and large restraining sphere is placed around the chromosomes. A set of short Langevin dynamics simulations with a decreasing radius of the restraining sphere brings the size of the model nuclei to their specified final values. The resulting configurations of all four chromosomes are presented in Fig. 1 and are used as the initial configurations for the simulations of the interphase chromatin dynamics.

Nuclear envelope and its interactions with chromatin chains

We model the NE, inner surface of which is covered by the nuclear lamina (NL), as a spherical boundary that restricts the motion of the chromosomes to the inside of the nucleus, and can attract LADs [37, 47,48,49]. Following [83], we map positions of 412 LADs (a median size of LADs is about 90 kb [47]) onto the chains of 1169 TADs. If TAD contains LAD, then the corresponding bead can attractively interact with the NL (NE). After the mapping, we have determined 350 TADs that contain LADs (we will call them L-TADs). Attractive interactions of LADs with the NL is due to the affinity of LADs to a meshwork of filaments formed by B-type lamins, with the average edge length smaller than 200 nm [38]. Furthermore, the LAD binding cites on the filaments are on the central \(\alpha\)-helical rod domains of B-Lamins [118], which are largely overlapping due to shifting of the lamin ends by 21–27 nm along the filaments [38], which is about 6 to 9 times smaller than the average bead size in our model. These considerations justify the coarse-graining of the LAD–NL interactions from a high-resolution, discrete molecular scale to the continuum scale in our model, approximating LAD–NL interactions as smooth, uniform interactions of L-TADs with a structure-less inner surface of the NE. We describe all the L-TAD–NE attractive interactions by the LJ-cos potential (see Eq. S2 Additional file 1) with a single well depth parameter \(\epsilon _L\). This is another simplification of the L-TAD–NE interactions (eliminating of the dependence of these interactions on the length of LAD in the L-TAD) which allows us to investigate what other factors can affect L-TAD–NE contact probabilities and radial distributions.

Fraction of LADs at the NE is calculated as the average (over the described above ensemble of 18 nuclei) fraction of L-TAD beads within (position of bead centers) 0.09 \(\mu\)m layer (average bead radius) at the NE.

To analyze the mobility of LADs we compute the probabilities of L-TADs (the centers of L-TAD bead) to be within a very thin, 0.2 \(\mu\)m layer immediately adjacent to the NE. The thickness of this layer roughly corresponds to the average diameter of beads in our model. When L-TAD is in this layer we define it as being in contact with the NE.

Interactions of the NE with TADs not containing LADs are described by the purely repulsive potential (see Eq. S3 in Additional file 1).

Lamin Mutant model

We refer to the lamin depleted or lamin knock-down nuclei as “Lamin mutant” nuclei. Our Lamin mutant nucleus model has the L-TAD–NE affinity parameter \(\epsilon _L\) reduced to essentially zero (0.1 kT, see a more detailed description in Additional file 1). All other interaction parameters are unchanged from the wild-type (WT) model. The Hi-C map for Lamin mutant model is calculated using long “11 h” trajectories (see “Matching simulation time to biological time” subsection in “Methods”.)

Experimental chromatin density profiles

The results of the simulations are compared with the experimental chromatin density data [119], briefly described in Additional file 1.

Simulation of chromatin dynamics

Langevin dynamics simulations have been performed using ESPResSo 3.3.1 package [120] by solving the Langevin equation of motion:

$$\begin{aligned} m_i \ddot{\textbf{r}_i} = -\nabla _i U - \gamma m_i \dot{\textbf{r}_i} + \textbf{L}_i(t). \end{aligned}$$
(1)

Here, \(\textbf{r}_i\) is the position of bead i with mass \(m_i\), U is the potential energy of the system. The last two terms describe the interaction with the solvent: a velocity dependent friction force, characterized by the parameter \(\gamma\), and a random Gaussian white noise force \(\textbf{L}_i\). The integration time step is \(t_{step}=0.01 \tau\), where \(\tau =\sigma _{min}\sqrt{m_{min}/\epsilon }\) is the LJ time scale [104, 121]. Using \(\epsilon =3\) kT we obtain \(\tau =136\) ns and, accordingly, \(t_{step}=1.36\) ns, which was used for all our simulations. For the simulations used to narrow down the range of the interactions parameters, we generated \(40\times 10^6\) time steps trajectories with \(\gamma = 1/\tau\). For the production simulations with the selected WT and Lamin mutant interaction parameter sets, we generated \(400\times 10^6\) time steps trajectories using \(\gamma = 0.01/\tau\). This much smaller \(\gamma\) allows one to effectively speed up the simulations (\(\sim 20\) times), which is one of the key benefits [122] of the implicit solvent approach used here. Also, in this regime of small friction the bead inertia and, accordingly, the difference in bead masses may become important.

Contact probability (Hi-C) map

The TAD–TAD contact probability map (Hi-C map) of the model is defined as the average of Hi-C maps (1169x1169) calculated over 18 trajectories of the ensemble of 18 model nuclei of 3 different sizes and 4 mutual arrangements of the chromosomes. (See a more detailed description in Additional file 1).

The final Hi-C maps for the developed WT and Lamin mutant models are calculated using long “11 h” trajectories (see the following subsection).

Matching simulation time to biological time

To relate the simulation timescale with the experiment [123], that is to map the simulation time onto real biological time, we compared [83] a diffusive motion of model beads with the experimental interphase chromatin diffusion, and used the match to estimate the scaling factor \(\lambda\) that converts the simulation time to real biological time. We calculated time dependencies of distance \(R_{i}(t)\) between bead i and the nucleus center for 9 randomly selected beads, which do not contain LADs. The dependence of the mean squared displacement (MSD) on the time interval \(\Delta t\), \(\langle \Delta R^2_i(\Delta t) \rangle = \langle [R_i(t+\Delta t) - R_i(t)]^2 \rangle\), was calculated for 3 different values of the Langevin friction parameter \(\gamma\) (\(0.01/\tau\), \(0.1/\tau\) and \(1/\tau\)). The averaging was performed over 9 selected beads along the 18 trajectories. We fit the first \(300 \cdot 10^3\) time steps of each of the three curves \(\langle \Delta R^2_i(\Delta t) \rangle\) (see Fig. S3 in Additional file 1) with the following equation for a sub-diffusive motion of chromosomal loci [123]:

$$\begin{aligned} \langle \Delta R^2_i(\Delta t) \rangle = 4D_{\text{app}} \left( \Delta t/\lambda \right) ^{0.39} , \end{aligned}$$
(2)

where \(D_{app}\) is the apparent diffusion coefficient, and 0.39 is the average experimental exponential scaling parameter. This equation (with \(\lambda =1\) and \(\Delta t\) in seconds) describes the experimentally observed diffusion of chromosomal loci over the time periods ranging from 1 to \(10^3\) s [123]. Very similar values for the exponential time scaling parameter have been obtained in the scaling theory (0.40) and in the accompanied dissipative particle dynamics simulations (0.38) for the thermal motion of monomers in the fractal globule state [124]. A slightly different value for the time scaling exponent for a single chromosomal locus MSD (0.52) has been recently reported in Ref. [125], consistent with the range seen in earlier experiments [123].

We obtain reasonable fits of the simulation derived MSDs with \(4D_{app}=0.061\) \(\mu\)m\(^2\), \(\Delta t\) in the number of time steps and \(\lambda =10^4\) s\(^{-1}\) for \(\gamma =0.01/\tau\), \(\lambda =2.8\cdot 10^4\) s\(^{-1}\) for \(\gamma =0.1/\tau\), and \(\lambda =20\cdot 10^4\) s\(^{-1}\) for \(\gamma =1/\tau\). These values of \(\lambda\) give us the number of simulation time steps that corresponds to 1 s of a real biological time. The scaling translates the \(40\times 10^6\) time step simulations with \(\gamma =1/\tau\) into 3 min of real nucleus time, and the long production \(400\times 10^6\) time step simulations with \(\gamma =0.01/\tau\) as corresponding to 667 min (11 h) of real time. The inevitable uncertainty in the experimental value of the scaling exponent used in Eq. 2 to map the simulation time onto real time suggests that one has to exercise caution in interpreting this mapping. We are nevertheless confident that a minute of the simulation time is still of the same order of magnitude as a minute of real biological time, not an hour or a second.

The Hi-C map calculated using long “11 h” trajectories for the WT parameter set has Pearson’s correlation coefficient with the experimental Hi-C map 0.954, close to 0.956 value obtained using “3 min” trajectories. These values suggest that reducing friction coefficient \(\gamma\) to 0.01 of its original \(1/\tau\) value employed at the stage of model parameters development, does not affect much of the local chromatin structure. A small deterioration of the original structure is expected due to time evolution and decay of the initial configurations.

Within the model, \(t=0\) corresponds to a point in the very beginning of G1 phase of the cell cycle, when all of the chromosomes are fully de-condensed.

Model development

Our model consists of 1179 “soft” beads that represent homologous pairs of TADs (1169 beads), pericentromeric constitutive heterochromatin domains (6 HET beads) and centromeric chromatin domains (4 CEN beads) (see “Methods”). The beads representing TADs are split into four types which correspond to four major epigenetic TAD classes: Active, Null, PcG and HP1 [1]. The beads are combined into four chains of homologous chromosomes (Chrs 2, 3, 4 and X) using “beads-on-a-string” model [68, 82, 101, 102]. Nucleolus is modeled as an additional constrained bead [26].

The nuclear envelope (NE), a spherical boundary surrounding the chromosome chains (see Fig. 1), constrains the motion of beads within the nucleus by TAD–NE repulsive interactions. At the same time, the NE, internal surface of which is lined by the nuclear lamina (NL)—a dense meshwork of protein filaments, can attract L-TADs (TADs that contain LADs) [43, 47]) representing LAD–NE interactions.

Compared to the previous fruit fly interphase chromatin model [83], and models developed for mammalian nuclei Refs. [52, 55, 56, 126], our model utilizes the attractive interactions between non-bonded beads. We assume that this presumably protein mediated TAD–TAD attraction [127, 128] is a TAD-class dependent [126], since different combinations of proteins are bound to different epigenetic classes of TADs [106] with their characteristic histone modifications.

Experimentally observed compartmentalization of chromatin [28] into active euchromatin (type A compartments) and more densely packed inactive heterochromatin (type B compartments) allows us to assume [55, 126] at least three types of TAD–TAD interactions: A–A, B–B and A–B. To allow the compartmentalization, some general relations (the so-called Flory–Huggins rule [129]) between these three types of TAD–TAD interactions have to be considered.

We also introduce a set of “specific” non-bonded TAD–TAD interactions for 268 TAD (bead) pairs which form “long-range” contacts with the increased probabilities [1]. These attractive interactions are described by a set of effective potentials with the well depths determined from the experimental contact probabilities for the corresponding TAD pairs (see “Methods” in Additional file 1).

Interaction parameters of the models

The attractive protein mediated interactions between TADs can spread over a several kT range of energies [55, 130].

To determine “optimum” bead–bead and L-TAD–NE interaction parameters of the model interphase nuclei we use simultaneously the following three major criteria (selection rules):

  1. (1)

    Maximum possible Pearson’s correlation coefficient between model derived TAD–TAD contact probability map (model Hi-C map) and the experimental WT Hi-C map [1], reduced to the TAD–TAD contact probabilities [26]);

  2. (2)

    The fraction of LADs (LAD containing beads, L-TADs) which are in contact with the NE matches the experiment for the WT nuclei: 25% [43]);

  3. (3)

    The commonly used condition for the chromatin compartmentalization is satisfied. This is the so-called Flory–Huggins rule [129]—the strength of the interactions between different (A and B) unit types should satisfy the phase separation criterion: interaction A–B < (A–A + B–B)/2.

Three stages of model development are used (see the description in Additional file 1, “Model development”).

The model Hi-C map derived from the nuclei dynamics using the final “optimum” WT set of the interaction parameters has Pearson’s correlation coefficient with the experimental Hi-C map 0.956. The contact probabilities in the experimental map include intrinsic experimental noise, which can affect the correlation with the model Hi-C map. Excluding the low-probability contacts (lower than 0.03), which are most affected by the noise component, effectively divides the experimental Hi-C map into two strata [131]; the corresponding Pearson’s coefficient between model and experiment increases to 0.967.

The “optimum” WT set of the interaction parameters is shown in the second column of the Table 1. We use this set for 10x longer simulations reported in the next section. The resulting chromatin density profiles will be compared with the available experimental data presented in [119]. We will also investigate how sensitive is the chromatin structure to the deviations in the parameters.

Results

Dynamic model of chromosomes in fruit-fly at TAD resolution

We have developed a coarse-grained “beads-on-a-string” model of D. melanogaster female interphase nuclei at TAD resolution (the average TAD size is \(\sim\)100 kb), where a bead represents a pair of homologous TADs in paired homologous chromosomes. Four main types of beads corresponding to four major epigenetic classes of TADs [1]—Active, Null, PcG and HP1—are used. The arrangements of the chromosome chains (nucleus topologies) in the model nuclei correspond to the experimentally observed CIS and TRANS mutual arrangements of the L- and R-chromosome arms and the position of X chromosome [108], see Fig. 1 in “Methods”. Unless otherwise specified, the results are averages over the ensemble of all different chromosome topologies and nucleus sizes (18 systems in total, see “Methods”).

Fig. 2
figure 2

The temporal evolution of the WT chromatin configurations (top row), and the model derived TAD–TAD contact probability (Hi-C) maps (bottom row). From left to right: The model chromatin configurations at t=0 min (starting configuration), and at time points corresponding to 30 min, 3 h and 11 h. The configurations shown correspond to TRANS-X3S nucleus topology (see Fig. 1), used as an example. The Hi-C maps (bottom row) are each averaged over 5 min (\(3 \cdot 10^6\) time steps) time intervals, and over the ensemble of 18 independent trajectories of the model nucleus, generated with the WT parameter set. The averaging includes all of the starting topologies, see “Methods

The main capabilities of the model are illustrated in Fig. 2, where we present the snapshots of the temporal evolution of the nucleus along the production trajectory, corresponding to about 11 h of biological time. The initial (\(t=0\)) and the following (at \(t=\) 30 min, 3 h and 11 h) full ensemble averaged Hi-C maps are shown as well. The snapshots and the maps show the expected small decay [83] of the perfect original chromosome territories and local structures: the probabilities of TAD–TAD contacts within each chromosome arm slightly decrease with time (by about 5–10%), although there are regions of contacts where these probabilities increase (see the difference maps in Fig. S5 in Additional file 1). The presented configurations and the maps demonstrate that reasonably distinct territories still exist 11 h into the interphase, providing the first “sanity check” of our model. The overall conclusion from the bottom panel of Fig. 2 (see also the difference maps in Fig. S5 in Additional file 1) is that the Hi-C map of a fruit fly nucleus at TAD resolution experiences only relatively small changes at the time-scale of the interphase. We refrain from a more detailed analysis of these changes here, as our focus is the role of the NE.

The model is trained to agree with the experimental Hi-C map, Fig. 3 (top panel) [1, 26], and lamin-DamID data [43] for the WT nuclei, and to satisfy the general polymer physics restrictions imposed on the strengths of the attractive interactions between beads of different types. The model derived TAD–TAD contact probability map, i.e., the model Hi-C map, Fig. 3 (bottom panel), has been calculated for the ensemble of “young” nuclei, corresponding to about 3 min of real biological time.

The selected WT parameter set provides a reasonably accurate reproduction of the experimental Hi-C map (Pearson’s correlation coefficient is 0.956). Some of the difference in the intensity of inter-arm TAD–TAD contacts between the experimental and model Hi-C maps can be explained by a noisy experimental data at this scale. We have also examined how the TAD–TAD contact probability decays as a function of the distance along the genome, within a single chromosome arm. Up to \(\sim\)1 Mb, the agreement with the experiment is quantitative, but between \(\sim\)1 and \(\sim\)20 Mb, the experimental contact probability decay is slower, indicating the presence of more compact structures at this length scale compared to the model prediction.

Fig. 3
figure 3

Top panel: Experimental Hi-C map for TAD–TAD contact probabilities in the WT D. melanogaster embryonic nuclei [26] (original data for embryonic nuclei from Ref. [1]). Bottom panel: Model derived TAD–TAD contact probability (Hi-C) map corresponding to young (0–3 min) nuclei. This Hi-C map reproduces many features expected from experiment: increased interactions within chromosomal arms, long-range chromatin contacts visible as bright spots, genome compartmentalization manifested as plaid-patterns of TAD–TAD contacts, and the Rabl-like configuration represented by interactions between chromosomal arms (e.g. 2 L and 2R) as “wings” stretched perpendicular to the main diagonal. Pearson’s correlation coefficient with the experimental Hi-C map is 0.956

The model is validated against multiple (ten) structural features of chromatin from several different experiments, which were not used for its training, see below. The model is capable of predicting the evolution of the 3D chromatin architecture (see Fig. 2) on time scales up to 11 h, which is enough to approximate the duration of the G1 phase of the fruit fly nuclei interphase [30, 132,133,134,135]. The Hi-C map calculated using “11 h” trajectories for the WT parameter set has Pearson’s correlation coefficient with the experimental [1] Hi-C map 0.954.

Fig. 4
figure 4

Radial chromatin density distributions in fruit fly nucelus. Top panel: Model WT nuclei (red curve, diamonds) and the Lamin mutant model nuclei (green curve, circles). The error bars are the standard deviations of the mean values for 18 nuclei. Bottom panel: Experimental WT nuclei (red curve, diamonds) and the Lamin mutant nuclei (green curve, circles) [Adapted from Supplementary Data, Figure S3 (Group 1) of Ref. [119]]

WT nucleus vs. Lamin mutant nucleus

To assess the role of the NE in organization of the chromatin architecture, a Lamin depleted nucleus model (simply referred to as “Lamin mutant”) has been created by reducing the LAD–NE affinity from relatively strong 4 kT WT value to essentially zero, see “Methods”. The rest of the model parameters were kept the same.

Two model derived WT and Lamin mutant radial TAD density distributions are shown in Fig. 4 (top panel). These distributions are in reasonable agreement with recently reported experimental chromatin density distributions [119] shown in Fig. 4 (bottom panel). In particular, the model faithfully reproduces non-trivial nuances of all of the experimental distributions reported in Ref. [119], particularly the “flatness” of the of density profile away from the NE for Lamin mutant (in experimentally seen distributions), as opposed to, e.g. continued growth of the density toward the center of the nucleus that can occur if the model parameters deviate from their optimal values, see below. This agreement provides a strong validation to the model, independent of the experimental data (WT Hi-C map) used for its training.

Comparing the WT and Lamin mutant chromatin distributions, Fig. 4, one can see a substantial shift of the chromatin density toward the NE in the WT nuclei. LAD–NE attractive interactions in the WT nuclei transform a more compact globule-like distribution in the Lamin mutants, with a very low chromatin density near the NE, into a more extended distribution, with a chromatin density peak near the NE and substantially reduced chromatin density in the central nucleus region. Similar effects of increased chromatin compaction and its movement away from the NE upon reduction of the LAD–NE attraction strength were previously observed in the experiments with Drosophila S2 cells [35], in the models of single human chromosomes [52, 56, 59] and in the model of several mouse chromosomes [55].

Comparing the model derived Hi-C maps for the Lamin mutant (Fig. 5) with that of the WT nuclei, our first conclusion is that the two are quite similar: the Pearson’s correlation coefficient between the two is 0.9989. Which means that, by and large, the significant global re-arrangements of the chromatin upon abrogation of the LAD–NE attraction, Fig. 4, have relatively little effect on the overall structure TAD–TAD contacts, as revealed by the Hi-C map. The relatively subtle differences between the two Hi-C maps are best revealed in the difference map, see Additional file 1: Fig. S8. One can see that the small (within 10%) changes in the contact probabilities are in agreement with the changes in the chromatin density distribution upon transition to the Lamin mutant seen in Fig. 4. A more compact and dense chromatin in the Lamin mutant has slightly higher TAD–TAD contact probabilities, both within the chromosome arms and between the arms.

Also, the specific TAD–TAD contacts produce more intensive spots on the Lamin mutant map (Fig. 5) compared to the WT model Hi-C map (see the difference map in Additional file 1: Fig. S8).

Fig. 5
figure 5

Model derived Lamin mutant TAD–TAD contact probability (Hi-C) map. Pearson’s correlation coefficient relative to the corresponding WT model Hi-C map is 0.9989, suggesting a high overall similarity between the WT and Lamin mutant contact maps

The model also reproduces key qualitative result from Ref. [35] that lamin depletion enhances interactions between active and inactive chromatin, impairing spatial segregation of active and inactive compartments. Our models provide a quantitative estimate for the increase of contacts between Null and Active TADs in the Lamin mutant relative to the WT nuclei (see the sums of the selected contact probabilities for each Null TAD with Active TADs in Fig. S11 in Additional file 1). The averaged relative increase is 22%. This number is a testable prediction; the predicted increase is yet another validation of the model, independent of the data sets used to build the model.

NE has a dual role, acting as an “attractive enclosure”

Comparing the model derived Lamin mutant and WT chromatin distributions (Fig. 4, top panel) one can see that the main consequence of the almost complete elimination of the LAD–NE affinity is a significant “global” change in the chromatin density profile. The peak of the WT chromatin density at about 1.9 \(\mu\)m from the nucleus center (at about 94% of the relative distance from the center) disappears, and the chromatin is shifted away from the NE. The density in the nearest to the NE layer of TADs drops 6 times compared to the WT value, while the density in the central region of the nucleus doubles, making the resulting chromatin distribution more compact. A similar general conclusion about the ability of LAD–NE interactions to re-distribute large blocks of chromatin was reached in the context of models of single human chromosomes interacting with the NE [52, 59]. We note that the respective parameter regimes that lead to this common conclusion, e.g., the ratios of TAD–TAD to LAD–NE interaction strengths, may be quite different from ours. On the other hand, it is reassuring that certain general principles of chromatin organization appear to be robust to model details.

What we have found unexpected is how quickly—on the minute time-scale—the more compact chromatin structure of the Lamin mutants deteriorates once the entire confining NE (not just Lamins) is removed completely, see Fig. S6 in Additional file 1. These combined results suggests a dual role for the NE. It is not about mere confinement of chromatin; but also it is not just “keeping interphase chromosomes slightly stretched” [35]. The NE acts as an “attractive enclosure”, which can both redistribute the chromatin, shifting most of it from the interior to the nuclear periphery due to the LAD–NE attraction in the WT nuclei, and, at the same time, the NE confines the chromatin, preventing it from decondensation in Lamin mutants, in the absence of the LAD–NE attraction. Without the NE, the chromatin will not shrink as it does in Lamin mutants. Possible biological implications of these findings are touched upon in “Discussion”.

3D chromatin architecture is sensitive to variations in the interactions of its key elements

To investigate how different interaction types affect the chromatin properties we have simulated the nuclei with the interaction parameters deviating from the WT and the corresponding Lamin mutant set values.

Effect of LAD–NE affinity variation

One of the question which clarifies the role of the NE is how the deviation of the LAD–NE affinity from its WT value (4 kT) affects the chromatin distribution? In Fig. S7 in the Additional file 1 (top panel), the radial distributions of TADs in the model at different levels of LAD–NE affinity (from 0.1 kT in Lamin mutant to 5 kT) are presented. One can see that the increase of LAD–NE affinity by 0.5 kT from its WT value leads to a substantial (\(\sim 2\)x) decrease of the TAD density in the central nucleus region and, at the same time, to a 20% increase of the chromatin density peak at the NE.

A more interesting behavior can be observed when LAD–NE affinity is decreased. The 0.5 kT decrease of the affinity results in a substantial reduction of the density peak value near the NE relative to the chromatin density in the central nucleus region (from 1.8 to 1.1) and its 0.2 \(\mu\)m shift from the NE. A further decrease of the affinity to 3 kT, which is only a 25% change, leads to a complete transformation of the chromatin density profile to the Lamin mutant-like profile: the density peak at the NE disappears, a substantial amount of the chromatin moves away from the NE, the interior density increases by 80% compared to its WT value. The summary of these effects is presented in Table 1. To the best of our knowledge such a sensitivity of chromatin distribution on the strength of LAD–NE interactions was not observed in the previous studies of mammalian nuclei [55].

Table 1 Wild-type (WT) model TAD interaction parameters and the effects of their variations on WT Chromatin (CHR) distribution and Chromosome arms contacts

The strong dependence of the chromatin density profile on the decrease of LAD–NE affinity [see Additional file 1: Fig. S7 (top panel)] suggests that both the number of LADs in the fruit fly genome and their affinity to the NE are likely “tuned” to be rather close to values corresponding to a transition between the Lamin mutant-like and WT-like chromatin distributions.

Effects of TAD–TAD interaction variation

We described in “Methods” how we selected the values of the TAD–TAD interaction parameters which produce a good agreement with the experimental Hi-C map [1, 26] and LAD distribution [43]. Here we will discuss what effects the deviations from those optimal parameters will have on the chromatin distribution and the TAD–TAD contact probabilities. The summary of these effects is presented in Table 1.

Active–Active TAD interactions: A large relative (400%), but small in the absolute value, increase (from 0.1 to 0.5 kT) does not lead to a noticeable change in the TAD radial distribution. Small variations of the intra-arm and inter-arm TAD–TAD contact probabilities are observed.

Null–Null TAD interactions: Decreasing the amplitude of these interactions from 1.5 kT (WT value) to 0.1 kT produces a small effect on the radial chromatin distribution—a 20% increase of the chromatin density in the central nucleus region, and a small decrease of the density near the peak at the nucleus periphery, Fig. S7 in Additional file 1 (bottom panel). Similar effect of partial chromatin redistribution toward the nuclear center upon reduction of the attraction between B-type beads has been observed in the single chromosome model of chromatin described in Ref. [52]. This small redistribution of the chromatin in our model results in a minor decrease of intra-arm and inter-arm TAD–TAD contact intensities. On the other hand, the increase of the Null–Null interactions by 1 kT (to 2.5 kT) leads to a substantial shift of the chromatin from the NE decreasing the fraction of L-TADs at the NE from 25 to 16% by transforming a more narrow main peak near the NE at 1.5\(-\)1.9 \(\mu\)m into a wider peak at 0.9\(-\)1.7 \(\mu\)m, Fig. S7 in Additional file 1 (bottom panel). Similar effect of the shift of chromatin density toward the central region upon increase of the mutual attraction between B-type beads, when some of them have affinity to the NE, has been observed in the one-chromosome chromatin model described in Ref. [59]. At the same time, the effect predicted in Ref. [52] is the opposite, at least at the same ratio of LAD–NE to B–B attractions as used in our model.

The shift of the chromatin density peak in our model when increasing the Null–Null attraction to 2.5 kT [see Fig. S7 in Additional file 1 (bottom panel, blue line)] is accompanied by a noticeable increase of intra-arm and inter-arm TAD–TAD contacts probabilities with a reduction of Pearson’s correlation coefficient with the experimental Hi-C map from 0.956 to 0.934. Since 228 of 492 Null TADs contain LADs (the most enriched in LADs class of TADs, see Fig. S1 in Additional file 1), these changes suggest a competition between the LAD–NE and Null–Null (TAD–TAD) interactions for the “optimal” chromatin structure. For comparison, only 54 of 494 Active TADs contain LADs, and the other two less numerous TAD classes contain even smaller numbers of L-TADs (50 of 131 PcG TADs and 18 of 52 HP1 TADs). This suggests that the competition between their TAD–TAD and LAD–NE interactions seems to be less important in the formation of the radial chromatin structure.

PcG–PcG TAD interactions: As in the case of Null–Null interactions, the decrease (from 1.5 to 0.1 kT) of the attraction between PcG TADs results in minor changes in the chromatin distribution and TAD–TAD contact probabilities. Due to the relative paucity of the PcG TADs we did not pursue further analysis of their selective influence on the chromatin structure.

Fig. 6
figure 6

Chromatin density distributions in the WT and Lamin mutant nuclei models with modified levels of cross-type TAD–TAD attractive interactions (parameter \(\epsilon _{AB}\)), both between A and B types of beads and between different types of B beads (e.g. Null–PcG interaction)

Cross-type TAD–TAD interactions in Lamin mutant vs. WT: Natural separation of TADs into active and inactive compartments suggests that cross-type TAD–TAD interactions would be detrimental to the cell. Does the WT level of LAD–NE affinity make chromatin density profile robust to changes in strength of cross-type TAD–TAD interactions? Decreasing the cross-type interactions from 0.5 kT (WT value) to 0.1 kT produces practically no change in the WT nuclei and minor changes in the Lamin mutant chromatin distribution (10–15% density decrease in the central regions due to small expansion of the chromatin toward the NE), Fig. 6. A small decrease of intra-arm and inter-arm TAD–TAD contacts is observed in both WT and Lamin mutant nuclei. On the other hand, increasing these interactions by 0.5 kT (to 1.0 kT) results in a substantial (80–100%) increase of the chromatin density in the central regions of the Lamin mutant due to further compaction of the chromatin in the globule-like state. The same change of the interactions in the WT nuclei leads to only a small increase of the density (about 20%) in the central regions, Fig. 6. A significant increase of intra-arm TAD–TAD contacts and a moderate increase of inter-arm TAD–TAD contacts are observed in both WT and Lamin mutant nuclei, as it can be seen on the Hi-C map differences (see Fig. S8 in Additional file 1).

The substantial chromatin density changes in the Lamin mutant and the small corresponding changes in the WT nuclei demonstrate that in the absence of the attractive LAD–NE interactions the chromatin distribution is very sensitive to the small changes in the cross-type TAD–TAD interactions. The reduction of this sensitivity in the presence of LAD–NE interactions suggests a stabilizing role of these interactions in maintaining native chromatin distribution in the WT nuclei and preventing cells from potentially negative effects of cross-type TAD–TAD interactions.

Probability of a TAD to be near the NE is determined by a highly variable local linear density of L-TADs (LADs) along the genome

Fig. 7
figure 7

Probabilities of TADs to be found within 0.4 \(\mu\)m layer adjacent to the NE (half the nuclear volume) for the WT (red) and Lamin mutant (green) model nuclei. Statistical error bars are smaller than symbol size. The WT TADs can be partitioned onto two major groups: the TADs with the probabilities greater than 0.5 and the TADs with the probabilities in the range 0.1\(-\)0.5. The TADs in the second group belong to the sections of chromosomes containing relatively long continuous stretches (10–16 TADs) without LADs. Null L-TADs #15, analyzed in [35] as cytological region 22A, is marked by yellow circles (for the WT and Lamin mutant nuclei). Null L-TADs #120, described in [35] as cytological region 36C, is marked by red triangles. PcG L-TAD #435, analyzed in [35] as cytological region 60D, is marked by blue squares. These three regions have been experimentally shown to detach from the NE in the Lamin mutant nuclei [35]. Our model agrees with the experiment for these three L-TAD regions, providing an additional validation for the model. (Also see model derived cumulative frequencies of radial positions of these regions in Fig. S12 in Additional file 1 for details)

To investigate the dynamic positioning of individual TADs within the nucleus relative to the NE, and the role of LAD–NE interactions in this positioning, we have partitioned the nucleus into two spherically symmetric compartments of equal volume: the central spherical region of 1.6 \(\mu\)m radius, and the spherical layer (0.4 \(\mu\)m thickness) adjacent to the NE. In the WT nuclei this 0.4 \(\mu\)m layer contains about 70% of all L-TADs. The probabilities of individual TADs to be in that adjacent to the NE layer for the WT and Lamin mutant nuclei are shown in Fig. 7.

One can see that in the WT nuclei the TADs are partitioned into two major groups: TADs that have high (0.5-\(-\)0.97) probability to be in the 0.4 \(\mu\)m layer near the NE (outer half the nuclear volume), and the ones that tend to stay away from the NE, with the probability in the range 0.1-\(-\)0.5 to be in that “near NE” layer. The TADs in the second group are in the sections of the chromosomes containing relatively long continuous stretches (10–16 TADs) without LADs. Quantitatively, one can characterize linear L-TAD density along a chromosome chain (L-TAD frequency of occurrence), \(f_L\), as the ratio of the numbers of L-TADs to the number of all TADs in a given chromosome section.

Fig. 8
figure 8

Probability of a TAD to be near NE is largely determined by the linear density of L-TADs surrounding this TAD. Shown are the probabilities of TADs/L-TADs to be found within 0.4 \(\mu\)m layer adjacent to the NE (upper trace), \(p_{NE}\), and the linear L-TAD density, \(f_L\), along the chromosome chain (lower trace). For visual clarity, only a fraction of 2 L chromosome is shown

We see a clear correlation between the probability (\(p_{NE}\)) of a TAD to be in the 0.4 \(\mu\)m “near NE” layer and the \(f_L\) value (see Fig. 8). For example, the first dip in the probability distribution for the 2 L-arm TADs from #70 to #72 (\(p_{NE}=0.13\), Fig. 7 and Fig. 8) is in the stretch of 17 TADs with only a single L-TAD (#78), and has the L-TAD density \(f_L = 0.05\). This density is 6 times lower that the genome-average value \(f_L = 0.30\). On the other hand, the regions in the first group (high \(p_{NE}\)) have a noticeably higher than average L-TAD density: in 2 L-arm–TADs from #156 to #166 (\(p_{NE} > 0.8\)) have the average \(f_L = 0.55\).

This stable separation of the average radial positions of different TADs suggests that the chromatin, despite being in a liquid-like state, has some “averaged global structure” determined by the irregular distribution of LADs (L-TADs) along the chromosomes.

In the Lamin mutant nuclei, the probabilities of TADs to stay closer to the NE (in the 0.4 \(\mu\)m “near NE” layer), presented in Fig. 7, are in the range 0.1-\(-\)0.4. The overall reduction of these probabilities compared to their WT values, reflecting the change of the TAD’s average radial positions, is in agreement with the detachment of the chromatin from the NE and its compactization in the Lamin mutant nuclei seen on the chromatin density distributions, Fig. 4.

Mobility of LADs. The nature of their stochastic distribution

It is well-known that the subset of LADs found at the NE differs substantially from cell to cell. Here we ask if LADs in an individual fruit fly nucleus may also be mobile, and if so, to what extent.

There are 412 LADs in D. melanogaster nucleus [47]; they are approximately evenly distributed along the three largest chromosomes: X, 2 and 3. As discussed in “Methods”, we have mapped these LADs onto the 1169 TADs [1] and have found LADs in 350 TADs (L-TADs). The WT value of the LAD–NE affinity (4 kT) in our model leads to the 25% fraction of L-TADs being, on average, in the nearest to the NE layer. Increasing or decreasing the LAD–NE affinity relative to its WT value in our simulations leads to an increase or a decrease of this fraction (see Fig. S4 in Additional file 1). This sensitivity (\(\sim\)8% per 0.5 kT change) suggests a dynamic balance between the number of L-TADs attached to the NE at any given moment and the rest of L-TADs.

LADs in individual nuclei are highly dynamic (mobile). We find that L-TADs can frequently attach to and detach from the NE during the interphase, and are not permanently anchored to the NE in any given cell. Visualization of motion of five randomly selected L-TADs in each chromosome shows that they are indeed highly mobile on time-scale of 20 min (see Fig. S10 in Additional file 1 and Additional files 2, 3, 4: Movie). One can see that, unlike larger LADs in relatively large human nuclei [34, 48, 49], fruit fly LADs can quickly move though a significant portion of the nuclear volume.

Fig. 9
figure 9

L-TADs are highly mobile within a single cell. Shown are the probabilities of L-TADs to be in contact with the NE. None of the L-TADs stay at the NE all the time, the maximum probability being only 0.85, suggesting a dynamic nature of all LAD–NE contacts. Large red/blue circles are the averages over 6 nuclei with 4 different nucleus topologies (4 CIS and 2 TRANS). Smaller black/green circles are the probabilities for a single CIS-X6S nucleus topology. All nuclei are of 2 \(\mu\)m radius. Statistical error bars are smaller than symbol size. Despite all 350 L-TADs in the model having the same affinity to the NE, there is a significant spread of their probabilities to be near the NE. A non-negligible number of L-TADs, 12%, have a very low (less than 0.2) probability to be at the NE

To confirm the suggestion about a highly dynamic nature of LADs on the genome-wide level, we have computed the probabilities of all the L-TADs to be in contact with the NE (see the definition in “Methods”), Fig. 9.

To distinguish between cell-to-cell stochasticity and LAD mobility within an individual nucleus, we compare the probabilities computed for the ensemble of 6 nuclei with the probabilities computed for a single nucleus.

Note that a hypothetical scenario in which L-TADs were static within each nucleus would yield a qualitatively different—binary—distribution of L-TADs for a single nucleus in Fig. 9 compared to the 6 nuclei average: it would look like L-TADs having only 0 or 1 probability to be in contact with the NE.

From our genome-wide model prediction for all of the L-TADs, the probabilities to be in contact with the NE are lower than 0.85. The fact that the probabilities computed for a single nucleus are noticeably less than 1.0 supports the hypothesis that during the interphase all L-TADs (LADs) bind to and unbind from the NE, within one and the same nucleus (see also Fig. S10 in Additional file 1 and Additional files 2, 3, 4: Movie). Even those L-TADs with the highest level of NE-binding probability are not permanently attached to the NE.

Probabilities of L-TADs to be in contact with the NE vary widely. Another interesting observation is a rather wide spread of the L-TAD contact probabilities seen in Fig. 9. Despite the fact that the model makes a simplification by considering all of the L-TAD–NE affinities to be the same (neglecting the variations in LAD sizes), it predicts very different levels of L-TAD–NE contact probabilities for different L-TADs.

Similar to a more general TAD positioning in a relatively wide layer close to the NE (see Fig. 7), we suggest that this large difference in the L-TAD binding—the contact probabilities vary from 0.07 to 0.84—is mostly related to a very different linear L-TAD density, \(f_L\), (shown in Fig. 8) in the chromosome regions surrounding these L-TADs. For example, the \(f_L\) varies from 0.08 for L-TAD #78 to 0.7 for L-TAD #1015, with the corresponding probabilities to be found at the NE being 0.08 and 0.84, respectively. Note that both \(f_L\) values are significantly different from the average linear L-TAD density for all four chromosomes, which is 0.30.

Fig. 10
figure 10

Examples of qualitatively different predicted radial distributions of several TADs, and their responses to lamin depletion. Top panel: Radial distribution of the PcG L-TAD #435 (cytological region 60D in Ref. [35]) in the WT and Lamin mutant nucleus models. Middle panel: Radial distribution of the Null L-TAD #120 (cytological region 36C in Ref. [35]) in the WT and Lamin mutant model nuclei. Bottom panel: Radial distribution of the Active TAD #22 in the WT and Lamin mutant model nuclei. The distributions of L-TADs (LADs) are bi-modal: the NE-bound mode is characterized by a high and very narrow density peak right at the NE, and the diffusive mode exhibits a noticeable LAD density in the nuclear interior. The bi-modality is consistent with the mobile nature of LADs

Bi-modal radial distribution of L-TADs in WT nuclei. To further explore the predicted dynamic nature of L-TAD binding to the NE, we have computed radial distributions of the positions of two L-TADs (PcG L-TAD #435 and Null L-TAD #120) in the WT and Lamin mutant nuclei, see Fig. 10 (top and middle panels). One can see that for the WT nuclei, where L-TADs have 4 kT affinity to the NE, the distributions are bi-modal: they have a bound mode, characterized by a large narrow density peak at the NE, and a diffusive mode with a noticeable density in the nucleus interior. In the case of the Null L-TAD #120 (Fig. 10, middle panel), the WT diffusive mode has a greater amplitude, close to the average single TAD density in the nucleus (\(\sim\) 0.03 \(\mu\)m\(^{-3}\)). In the case of the PcG L-TAD #435 (Fig. 10, top panel), most of the WT diffusive mode has a much smaller amplitude (\(\sim\) 0.01 \(\mu\)m\(^{-3}\)) with more of the TAD position density localized in the bound mode. These bi-modal radial L-TAD distributions demonstrate further, that in the WT nuclei LADs are not permanently bound to the NE, but are instead rather mobile, being able to attach and detach during the interphase. Similar bi-modal distributions of LADs have been observed in a polymer model of a 81 Mbp genomic region of human chromosome 5 [56].

In addition to L-TADs, the shift of the radial positioning toward the NE in the WT nuclei is also observed for most other TADs (see Fig. 7). A typical behavior of radial distributions of these TADs is demonstrated using as an example the distribution for the Active TAD #22, depicted in Fig. 10 (bottom panel). The wide distribution peak of this TAD in the Lamin mutant nuclei is shifted toward the NE in the WT nuclei, which is in agreement with different probabilities of this TAD to be within the 0.4 \(\mu\)m layer near the NE (half the nuclear volume) for the Lamin mutant and WT nuclei: 0.30 and 0.62, respectively.

An interesting exception related to this general chromatin density shift toward the NE in the WT nuclei (compared to the Lamin mutant nuclei) are small groups of TADs (e.g., TADs with IDs in the ranges 64–80 and 312–371 in 2 L and 2R chromosome arms) that demonstrate small shifts in the radial density distributions toward the NE in the Lamin mutant. We suggest that this opposite pattern is related to a very low linear L-TAD density in these selected regions of chromosomes and, accordingly, a very low averaged attraction to the NE of these regions, compared to the neighboring regions with higher L-TAD densities. Due to a competition among L-TADs for a limited space near the NE in the WT nuclei, these special groups of TADs with a very low averaged attraction to the NE are effectively pushed away from the NE. Consistent with the above explanation, the differences in the attraction towards the NE are eliminated in the Lamin mutant nuclei, leading to a much smaller spread of TAD probabilities to be near the NE. The remaining small variation—upturns in the green trace at the ends of each chromosome arm in Fig. 7—is consistent with a greater probability to be at the nuclear periphery for the centromeric and telomeric regions of chromosomes.

The distributions of radial TAD positioning discussed above further suggest that the TADs are highly dynamic in both WT and Lamin mutant interphase nuclei.

Discussion

This work has two types of outcomes: methodological advances and biological predictions. The latter are focused on the role of interactions between structural components of chromatin such as TADs, as well as their interactions with the NE, in determining the local structure and 3D global architecture of chromatin, and their stability.

Methodological advances

We have developed a novel coarse-grained “beads-on-a-string” model of chromatin of the entire D. melanogaster interphase nucleus at TAD resolution (\(\sim\)100 kb). One of the major considerations in choosing the TAD level of coarse-graining is that TADs are conserved, stable units of the fruit fly chromatin, interactions between which determine compartmentalization of the chromatin into active euchromatin and more densely packed inactive heterochromatin [1, 15, 107] and, together with the LAD–NE interactions, its global distribution.

The physics-based, as opposed to purely data-driven, approach taken here allows us to answer many “what if” questions hard to address experimentally, focusing on the role of the chromosome–NE interactions on the global and local chromatin structure and its stability, which are the subject of many recent experimental studies [35, 40, 55, 91, 119, 136,137,138].

Compared to many existing physics-based models of chromatin in higher eukaryotes, our approach has several methodological novelties; these make a tangible difference with respect to biologically relevant predictions our model can make.

First and foremost, we are simulating an entire biological system—an ensemble of nuclei, corresponding to the experimentally observed set of mutual spatial arrangements of chromosome arm (CIS and TRANS), properly weighted according to experiment. As we have shown, considering the entire biological system within a physics-based model proves important to reproducing experiment, specifically to obtain a better agreement with the experimental Hi-C map derived from a very large set of nuclei. A related novel aspect of the developed model is that we consider four distinct epigenetic classes of TADs as four corresponding types of beads, as well as two additional bead types—pericentromeric constitutive heterochromatin (HET) and centromeric regions (CEN) beads.

The second key methodological novelty of our modeling approach is that we are able to simulate temporal evolution of the structure of all chromosomes in the nucleus on the time-scale of the entire G1 phase of the interphase. That a computer simulation of an entire fruit fly nucleus at \(\sim\)100 kb resolution can reach the biological time-scales of 11 h or so is not entirely obvious. Here, we have successfully adapted to the field of chromatin simulations an approach developed in Ref. [122]—implicit solvation with low Langevin friction term. This approach was found successful in the study of protein folding and similar problems in structural biology that faced the same problem: biologically interesting time-scales are far out of reach of traditional simulation techniques often used in this field. While we can not map the resulting time-scales precisely onto real biological time, we are confident that what we call “minutes” in our simulations are not hours or seconds of real biological time (see Fig. S3 in Additional file 1).

To arrive at the key interaction parameters that define the behavior of the model, we have used three major criteria, or “rules”: (1) maximizing Pearson’s correlation coefficient between model derived TAD–TAD contact probability map and the experimental Hi-C map [1, 26]; (2) a match between the model derived and experimental set-averaged fractions of LADs contacting the NE [43]; (3) the conditions for the chromatin A/B compartmentalization—the Flory–Huggins phase separation criterion [129]. The criterion (2) is likely a novelty in the development of coarse-grained models of chromatin, which leads to important biological conclusions, see below.

The proposed model has been validated extensively against multiple experimental observations and trends, independent from those used in the model construction. A few examples are summarized below. The model derived TAD–TAD contact probability map reproduces all the main qualitative features of the experimental Hi-C maps obtained from Drosophila cells [1, 15]. These include increased interactions within chromosomal arms, long-range chromatin contacts visible as bright spots located off the main diagonal, genome compartmentalization manifested as plaid-patterns of TAD–TAD contacts, and the Rabl-like configuration represented by interactions between chromosome arms as “wings” stretched perpendicular to the main diagonal.

Further, the model reproduces the key qualitative results from Ref. [35] that lamin depletion enhances interactions (increases contact frequency) between active and inactive chromatin and leads to a chromatin compaction. Experimentally observed detachment of several cytological regions from the NE in lamin depleted nuclei [35] and their radial positioning are also faithfully reproduced by the model.

In another independent validation of the model, we have demonstrated that it reproduces, automatically, experimental chromatin density profiles [119] of both the WT and Lamin mutant nuclei, including some highly nuanced features such the “flatness” of the Lamin mutant density profile away from the NE, Fig. 4, as opposed to, e.g., continued growth of the density toward the center of the nucleus that can be seen when the model parameters deviate from their optimal values, Fig. 6. The fact that the model can reproduce this nuanced behavior is non-trivial, as small variations of the model parameters destroy the agreement (while still predicting the more trivial behavior of the chromatin moving away from the NE in the lamin depleted nuclei).

The overall behavior of TADs we describe in the manuscript, such as their radial positioning and TAD–NE contact probabilities, depends of the density of neighboring LADs along the chromatin chain and epigenetic type of chromatin in a TAD, rather than the TAD size. For this reason, we believe that our calculations of chromatin density distribution, chromatin–NE contact probabilities, and estimates of the LAD dynamics, are robust to the selection of the average TAD size. In the future, re-calibrating the model to a finer resolution [139] will unlikely invalidate our main conclusions made in this work.

In summary, the model reproduces about ten distinct features of Drosophila interphase chromatin observed in several independent experiments not used in model construction.

Biological predictions and speculations

Our first noteworthy conclusion is that the positioning of all LADs in D. melanogaster interphase nuclei is highly dynamic (mobile)—on the time scale of the G1 phase of the interphase -- the same LAD can attach, detach, move far away from the NE and then re-attach itself to the NE multiple times. This prediction is supported by multiple computational experiments. In particular, the analysis of the distributions of radial positions of single L-TADs: the distributions have two modes—NE-bound and diffusive. Consequently, none of the L-TADs spends all of the time at the NE. This conclusion goes beyond what is known from experiment for fruit fly nuclei: that LADs found at the NE differ from cell to cell. What we show is that, in any given cell nucleus, LADs are highly dynamic. We argue that this prediction is robust, as it is an inevitable consequence of the relatively low strength of the LAD–NE attraction. The specific value, 4 kT, of this attractive energy used by our model is not arbitrary, it is derived from the experimental fact that only a certain, limited fraction of LADs is on average bound to the NE [43]. A hypothetically much higher value of the LAD–NE affinity that would “glue” all or most LADs to the NE would be inconsistent with the experimental data used to construct the model.

It is worthwhile to compare our genome-wide predictions for LAD mobility in fruit fly with the corresponding experimental findings, which, to the best of our knowledge, are available for human nuclei [34]. In both cases, LADs are not static, but in human nuclei, LAD movements are confined [34] to a relatively narrow layer near the NE, while in fruit fly we see a relatively higher mobility overall, with many L-TADs traversing the nucleus, from the NE to the center. A more detailed analysis is warranted to quantify the similarities and differences between the nature of LAD mobility in these two organisms.

Observations of the dynamic nature of LADs in interphase nuclei raise a question of the effect of a LAD being in close proximity to the NE on the expression of the genes within that LAD. More specifically, does the expression level vary as the LAD moves between the periphery and the other parts of the nuclear interior? A study systematically tested mammalian promoters moved from their native LAD location to a more neutral chromatin environment and to a wide range of chromatin contexts inside LADs [140]. The study has demonstrated that it is the features encoded in the promoter sequence and variation in local chromatin composition that determine gene expression levels in LADs [140]. If the interplay between promoter sequence and local chromatin features is sufficient to determine the level of transcription inside LADs, then gene expression may be robust to the dynamic nature of LADs, at least in WT nuclei. Future genome-wide studies of the spatial–temporal transcription inside the nucleus may answer this question; combining experiment and computer modeling may be beneficial.

Related to the above conclusion about the dynamic nature of LAD binding to the NE is the prediction that, despite all of the L-TADs in our model having exactly the same affinity to the NE, the probability of L-TADs binding to the NE varies by an order of magnitude between L-TADs. We explain this variation by the corresponding pronounced variation, up to 9 times, of the local linear L-TAD density along the chromatin chains, in contrast to an earlier suggestion that it is the highly variable LAD–NE affinities of relatively large LADs in human cells that may be responsible for the differences in the frequency of LAD binding to the NE [51]. A potentially biologically relevant consequence of our finding is that the genetic/epigenetic features of a given TAD alone can not fully determine its fate with respect to probability of being found near the NE, even if the stochastic component of the positioning is eliminated by averaging over time and an ensemble of nuclei. The distribution of LADs along the genome strongly affects the average radial positioning of individual TADs, playing a notable role in maintaining a non-random average global structure of chromatin, within its overall liquid-like state.

We also find that the specific strength of the WT value of LAD–NE attraction puts the chromatin very near the “phase boundary”, separating two qualitatively different chromatin density distributions; a mere 12% (0.5 kT) decrease of the LAD–NE affinity strength changes the shape of the chromatin density profile appreciably, from the WT one to one that resembles the Lamin mutant density profile. Changing the LAD–NE affinity by 25% from its WT value (1 kT decrease) results in a drastic (60%) decrease in the fraction of L-TADs at the NE. One proposed biological consequence of being on the “phase boundary” is as follows. If we assume that the \(\sim\) 12% (0.5 kT) variation in LAD–NE affinity occurs naturally, then the high sensitivity of the chromatin structure to the strength of LAD–NE affinity might explain variability of chromatin architecture between nuclei of the same tissue and between different tissues. Indeed, some LADs are conserved between cell types, while others are more variable [141]. LADs that display less consistency between cells in a population tend to be specific to cells where genes, located in these LADs, are transcriptionally repressed [142]. Cell type-specific genes located in variable LADs are released from the NE upon cell type differentiation [49].

Another set of model predictions focuses on the potential role of LAD–NE interactions in the sensitivity, or lack thereof, of chromatin 3D architecture to other key interactions (TAD–TAD), which together create the delicate balance that determines the nuclear architecture. Recent studies [26, 52, 55, 83, 143], including this one, leave little room to debate the importance of LAD–NE interactions in genome organization. As in previous works, e.g., on mouse [55], agreement of the polymer model with experiment can only be achieved in a rather narrow window of parameters that determine TAD–TAD and LAD–NE interactions in fruit fly nucleus. Our model goes further, by allows us to differentiate between the four main types of TADs: we find that among transcriptionally repressed TAD types, Null–Null interactions have the strongest effect on the 3D chromatin architecture. Also, as previously reported, one must assume a relatively weak mutual attraction between Active-type TADs. The fact that very different models applied to very different organisms, from fruit fly to mammals [52, 55], arrive at several similar general conclusions regarding the role of the interplay of the interactions between chromatin units and the NE, speaks for a certain degree of conservation of chromatin organization across species. We would like to note, however, that in a resent computational study of mammalian nuclei [55], a much larger (compared to our work) relative variation of model LAD–NE affinities (up to 3 times) was found compatible with the chromatin distribution in the WT nuclei.

In contrast to the predicted strong effect of small changes in LAD–NE interactions on the radial chromatin distribution, this distribution is rather insensitive to even relatively large (30–100%) changes in the strength of the interactions between TADs. The changes in the interactions between the Null TADs have the most effect, but even that effect is significantly smaller than the changes in chromatin density profile resulting from similar relative changes in the LAD–NE affinity.

Critically, in contrast to WT nuclei, chromatin density distribution in the Lamin mutant is predicted to be sensitive to increase of the strength of cross-type TAD–TAD attractive interactions. This comparison suggests that another role of relatively strong LAD–NE interactions is in providing a stable global environment with a low sensitivity to small fluctuations in TAD–TAD interactions. We speculate that the low sensitivity of global chromatin architecture (radial chromatin density distribution) to TAD–TAD interactions hints at the possibility of a mechanism in which changes in TAD–TAD interactions can be more important for local regulation of gene transcription activity. We also make a testable prediction that many aspects of the chromatin architecture will be more variable among cells with fewer LAD–NE contacts and even more so in the Lamin-depleted cells. These may include TAD–TAD contacts, chromatin radial distribution, Rabl configuration, spatial segregation of chromosome territories. In fact, a greater variability of the chromatin radial distribution in proventriculus nuclei of the Drosophila Lamin mutants in comparison to WT nuclei can be observed by comparing three experimental groups ( [119]). Thus, a dramatic loss or dysfunction of lamins during aging or disease [49, 62,63,64, 144] may contribute to the increased disorder in gene expression due to greater variability of the global chromatin architecture predicted in this work.

Turning off LAD–NE interactions (simulating Lamin mutant) results in almost complete detachment of the chromatin from the NE and its compaction, causing a twofold increase of the chromatin density in the central region of the nucleus. At the same time, complete removal of the NE in the model results in chromatin decompaction and separation of the chromosome territories on very short time-scale, orders of magnitude shorter that the duration of the interphase. This result suggests that, without the NE, the interactions between TADs are not strong enough to keep the chromatin of the Lamin mutant in a globule-like form on biologically relevant time-scales of hours. At the onset of mitosis, the NE is disassembled and the entire genome is condensed into mitotic chromosomes [145]. In organisms with an open mitosis, such as Drosophila (except syncytial embryonic divisions), NE reformation occurs by recruitment of nuclear pore complexes and membrane components to the surface of the segregating chromosomes. The process begins in late anaphase with the binding of nuclear pore complex proteins to chromosomes and is completed with the recruitment and fusion of membranes during telophase [146]. Thus, the enclosing role of the NE is already established at the beginning of the interphase. This process could have evolved to prevent chromosomes from further unfolding and detaching from each other.

Taken together, these observations suggest a dual mechanical role of the NE in the WT nuclei: it is not simply a confinement or “stretcher” of the chromatin, but, rather, the NE acts as an “attractive enclosure”, which simultaneously expands and confines the chromatin, while stabilizing both its local and global 3D structure.

Contrary to the chromatin density profiles, which reflect the global 3D chromatin architecture, the predicted Hi-C maps are not very sensitive to the changes in the LAD–NE affinity, suggesting that local chromatin structure is determined mostly by the TAD–TAD interactions. Indeed, increasing the interactions strength between the most numerous strongly interacting Null TADs by 1 kT, from their 1.5 kT optimum WT value, leads to a substantial increase in the number of intra-arm and inter-arm contacts. A similar effect on Hi-C maps is predicted for the change in interactions between TADs of different types.

By and large, the predicted Lamin mutant Hi-C map looks rather similar to that of the WT one (Pearson’s correlation coefficient between these maps is 0.9989); this similarity, within each chromosome, was previously observed in lamin knockdown nuclei [35]. A more detailed analysis shows that most TAD–TAD contacts, including inter-chromosome contacts, are slightly enhanced in the Lamin mutant compared to the WT nucleus (see the difference map in Additional file 1, Fig. S8), consistent with the overall compaction of chromatin observed experimentally in LamA25 mutant [119]. Small areas where the contact frequency decreases in the Lamin mutant are predicted to be limited to close to diagonal intra-arm contacts. These findings are consistent with predictions made for polytene chromosomes in fruit fly [77].

The overall increase of TAD–TAD contacts and compaction of the chromatin in the Lamin mutant may bring TADs of different chromosomes closer to each other, facilitating inter-chromosome interactions, which is exactly what we observe in our simulations. In particular, the Lamin knockdown increases the chromatin density in a fraction of TADs enriched in active chromatin, and enhances interactions between active and inactive chromatin [35]. Using our model, we can go further and quantify some of these changes. From the model Hi-C maps, we have estimated the sums of contact probabilities of each Null TAD (inactive B-type TADs) with the Active TADs in the Lamin mutant and WT model nuclei (see Fig. S11 in Additional file 1). We predict a noticeable, 22% on average, increase of these active–inactive chromatin contacts in the Lamin depleted nuclei. The reduction of active–inactive chromatin contacts in WT suggests a stabilizing role of the LAD–NE interactions in maintaining native chromatin distribution and preventing cells from potentially detrimental effects of cross-type TAD–TAD interactions.

Taken together, our modeling data indicate that LAD–NE interactions play a diverse and prominent role in 3D genome organization.

Conclusions

We have developed a coarse-grained model of D. melanogaster interphase nuclei at TAD (\(\sim\)100 kb) resolution that describes dynamics and time evolution of all fruit fly chromosomes, and their interactions with the nuclear envelope (NE), on the time-scale of the entire interphase. The model takes into account different types of TAD–TAD interactions between different epigenetic classes of TADs, attractive interaction between LADs and the NE, and is tuned to reproduce the experimental Hi-C map and the fraction of LADs positioned at the NE. Several methodological novelties proved important to achieve good agreement with experiment, including explicitly accounting for different experimentally observed mutual spatial arrangements of the chromosome arms (nucleus topologies). The model has been validated against multiple distinct features of Drosophila interphase chromatin, not used in the fitting of its parameters.

We have used the model to explore, in detail, how several key characteristics of the chromatin 3D architecture, including the overall chromatin density distribution and Hi-C maps, are sensitive to the interaction strength between different classes of TADs, and between LADs and the NE. Some of our general conclusions agree with previous findings based on models of mammalian nuclei, which supports conservation of several general principles of chromatin organization across species.

Multiple genome-wide predictions have been made in this work. We predict a very dynamic nature of binding of LADs to the NE in D. melanogaster interphase nuclei. We also predict an increased sensitivity of global chromatin architecture to the fluctuations in TAD–TAD interactions in lamin depleted nuclei compared to the WT, where relatively strong LAD–NE interactions suppress this sensitivity. The proposed model predicts that radial positioning of TADs in the nuclei, including the probabilities of TADs to be in the high density layer at the NE, are largely determined by local linear (along the chromatin chain) densities of LADs around TADs, suggesting a significant role of LAD distribution in average 3D positioning of TADs.

We conjecture that one important role of the distribution of LADs along the chromosome chains and their attractive interactions with the NE is to create a non-random average global structure of chromatin and to protect its integrity and stability against inevitable cell-to-cell variations in TAD–TAD interactions. We also predict greater variability of the chromatin architecture due to loss or dysfunction of lamins, which may contribute to the increased disorder in gene expression during aging or disease.

Availability of data and materials

The modeling code ESPResSo 3.3.1 [120] used in this research is available at http://espressomd.org/. The software is free, open-source, published under the GNU General Public License (GPL3). Visualization of motion of randomly selected L-TADs in each chromosome on a time-scale of 20 min is available at https://drive.google.com/drive/folders/1Dwfe7hbCjARxnvhf7oP76sVaTRa0RE6Z and as supplementary movies in Additional files 2,3 and 4. Input files for ESPResSo 3.3.1 to reproduce the simulations, and corresponding representative snapshots of the resulting chromatin structures are available at: https://people.cs.vt.edu/~onufriev/CODES/DROSOPHILA_NUCLEUS.zip.

Abbreviations

3D:

Three-dimensional

CEN:

Centromeric

Chr:

Chromosome

HET:

Pericentromeric constitutive heterochromatin

LAD:

Lamina-associated domain

L-TAD:

TAD that contains LADs

LJ:

Lennard–Jones

MSD:

Mean squared displacement

NE:

Nuclear envelope

TAD:

Topologically associating domain

WT:

Wild-type

References

  1. Sexton T, Yaffe E, Kenigsberg E, Bantignies F, Leblanc B, Hoichman M, Parrinello H, Tanay A, Cavalli G. Three-dimensional folding and functional organization principles of the Drosophila genome. Cell. 2012;148(3):458–72.

    Article  CAS  PubMed  Google Scholar 

  2. Bickmore W, van Steensel B. Genome architecture: domain organization of interphase chromosomes. Cell. 2013;152(6):1270–84. https://doi.org/10.1016/j.cell.2013.02.001.

    Article  CAS  PubMed  Google Scholar 

  3. Bickmore WA. The spatial organization of the human genome. Ann Rev of Genomics Hum Genet. 2013;14(1):67–84. https://doi.org/10.1146/annurev-genom-091212-153515.

    Article  CAS  Google Scholar 

  4. Bonev B, Cavalli G. Organization and function of the 3D genome. Nat Rev Genet. 2016;17(11):661–78. https://doi.org/10.1038/nrg.2016.112.

    Article  CAS  PubMed  Google Scholar 

  5. Solovei I, Thanisch K, Feodorova Y. How to rule the nucleus: divide et impera. Curr Opin Cell Biol. 2016;40:47–59. https://doi.org/10.1016/j.ceb.2016.02.014.

    Article  CAS  PubMed  Google Scholar 

  6. Ciabrelli F, Cavalli G. Chromatin-driven behavior of topologically associating domains. J Mol Biol. 2015;427(3):608–25. https://doi.org/10.1016/j.jmb.2014.09.013.

    Article  CAS  PubMed  Google Scholar 

  7. Barton LJ, Soshnev AA, Geyer PK. Networking in the nucleus: a spotlight on LEM-domain proteins. Curr Opin Cell Biol. 2015;34:1–8. https://doi.org/10.1016/j.ceb.2015.03.005.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Lemaître C, Bickmore WA. Chromatin at the nuclear periphery and the regulation of genome functions. Histochem Cell Biol. 2015;144(2):111–22. https://doi.org/10.1007/s00418-015-1346-y.

    Article  CAS  PubMed  Google Scholar 

  9. Nguyen HQ, Bosco G. Gene positioning effects on expression in eukaryotes. Annu Rev Genet. 2015;49:627–46. https://doi.org/10.1146/annurev-genet-112414-055008.

    Article  CAS  PubMed  Google Scholar 

  10. Sexton T, Cavalli G. The role of chromosome domains in shaping the functional genome. Cell. 2015;160(6):1049–59. https://doi.org/10.1016/j.cell.2015.02.040.

    Article  CAS  PubMed  Google Scholar 

  11. Wang Y, Nagarajan M, Uhler C, Shivashankar GV. Orientation and repositioning of chromosomes correlate with cell geometry-dependent gene expression. Mol Biol Cell. 2017;28(14):1997–2009. https://doi.org/10.1091/mbc.e16-12-0825.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Diament A, Tuller T. Modeling three-dimensional genomic organization in evolution and pathogenesis. Semin Cell Dev Biol. 2019;90:78–93. https://doi.org/10.1016/j.semcdb.2018.07.008.

    Article  CAS  PubMed  Google Scholar 

  13. Hafner A, Boettiger A. The spatial organization of transcriptional control. Nat Rev Genet. 2023;24(1):53–68. https://doi.org/10.1038/s41576-022-00526-0.

    Article  CAS  PubMed  Google Scholar 

  14. Dixon JR, Selvaraj S, Yue F, Kim A, Li Y, Shen Y, Hu M, Liu JS, Ren B. Topological domains in mammalian genomes identified by analysis of chromatin interactions. Nature. 2012;485(7398):376–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Hou C, Li L, Qin ZS, Corces VG. Gene density, transcription, and insulators contribute to the partition of the Drosophila genome into physical domains. Mol Cell. 2012;48(3):471–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Acemel RD, Maeso I, Gómez-Skarmeta JL. Topologically associated domains: a successful scaffold for the evolution of gene regulation in animals. Wiley Interdiscip Rev Dev Biol. 2017;6(3):265.

    Article  Google Scholar 

  17. Szabo Q, Bantignies F, Cavalli G. Principles of genome folding into topologically associating domains. Sci Adv. 2019. https://doi.org/10.1126/sciadv.aaw1668.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Wang S, Su J-H, Beliveau BJ, Bintu B, Moffitt JR, Wu C-T, Zhuang X. Spatial organization of chromatin domains and compartments in single chromosomes. Science. 2016;353(6299):598–602.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Phillips-Cremins JE, Sauria ME, Sanyal A, Gerasimova TI, Lajoie BR, Bell JS, Ong C-T, Hookway TA, Guo C, Sun Y, et al. Architectural protein subclasses shape 3D organization of genomes during lineage commitment. Cell. 2013;153(6):1281–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Rao SS, Huntley MH, Durand NC, Stamenova EK, Bochkov ID, Robinson JT, Sanborn AL, Machol I, Omer AD, Lander ES, et al. A 3D map of the human genome at kilobase resolution reveals principles of chromatin looping. Cell. 2014;159(7):1665–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Dixon JR, Jung I, Selvaraj S, Shen Y, Antosiewicz-Bourget JE, Lee AY, Ye Z, Kim A, Rajagopal N, Xie W, et al. Chromatin architecture reorganization during stem cell differentiation. Nature. 2015;518(7539):331–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Smith EM, Lajoie BR, Jain G, Dekker J. Invariant tad boundaries constrain cell-type-specific looping interactions between promoters and distal elements around the CFTR locus. Am J Hum Genet. 2016;98(1):185–201. https://doi.org/10.1016/j.ajhg.2015.12.002.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Renschler G, Richard G, Valsecchi CIK, Toscano S, Arrigoni L, Ramírez F, Akhtar A. Hi-C guided assemblies reveal conserved regulatory topologies on X and autosomes despite extensive genome shuffling. Genes Dev. 2019;33(21–22):1591–612.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Torosin NS, Anand A, Golla TR, Cao W, Ellison CE. 3D genome evolution and reorganization in the Drosophila melanogaster species group. PLOS Genet. 2020;16(12):1–29. https://doi.org/10.1371/journal.pgen.1009229.

    Article  CAS  Google Scholar 

  25. McArthur E, Capra JA. Topologically associating domain boundaries that are stable across diverse cell types are evolutionarily constrained and enriched for heritability. Am J Hum Genet. 2021;108(2):269–83. https://doi.org/10.1016/j.ajhg.2021.01.001.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Li Q, Tjong H, Li X, Gong K, Zhou XJ, Chiolo I, Alber F. The three-dimensional genome organization of Drosophila melanogaster through data integration. Genome Biol. 2017;18(1):145–66.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Hua N, Tjong H, Shin H, Gong K, Zhou XJ, Alber F. Producing genome structure populations with the dynamic and automated PGS software. Nat Protoc. 2018;13(5):915–26. https://doi.org/10.1038/nprot.2018.008.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Lieberman-Aiden E, Van Berkum NL, Williams L, Imakaev M, Ragoczy T, Telling A, Amit I, Lajoie BR, Sabo PJ, Dorschner MO, et al. Comprehensive mapping of long-range interactions reveals folding principles of the human genome. Science. 2009;326(5950):289–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Nora EP, Lajoie BR, Schulz EG, Giorgetti L, Okamoto I, Servant N, Piolot T, van Berkum NL, Meisig J, Sedat J, et al. Spatial partitioning of the regulatory landscape of the X-inactivation centre. Nature. 2012;485(7398):381–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Csink AK, Henikoff S. Large-scale chromosomal movements during interphase progression in Drosophila. J Cell Biol. 1998;143(1):13–22. https://doi.org/10.1083/jcb.143.1.13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Chubb JR, Boyle S, Perry P, Bickmore WA. Chromatin motion is constrained by association with nuclear compartments in human cells. Curr Biol. 2002;12(6):439–45. https://doi.org/10.1016/s0960-9822(02)00695-4.

    Article  CAS  PubMed  Google Scholar 

  32. Spector DL. The dynamics of chromosome organization and gene regulation. Annu Rev Biochem. 2003;72(1):573–608. https://doi.org/10.1146/annurev.biochem.72.121801.161724.

    Article  CAS  PubMed  Google Scholar 

  33. Lanctôt C, Cheutin T, Cremer M, Cavalli G, Cremer T. Dynamic genome architecture in the nuclear space: regulation of gene expression in three dimensions. Nat Rev Genet. 2007;8(2):104–15. https://doi.org/10.1038/nrg2041.

    Article  CAS  PubMed  Google Scholar 

  34. Kind J, Pagie L, Ortabozkoyun H, Boyle S, de Vries SS, Janssen H, Amendola M, Nolen LD, Bickmore WA, van Steensel B. Single-cell dynamics of genome-nuclear lamina interactions. Cell. 2013;153(1):178–92. https://doi.org/10.1016/j.cell.2013.02.028.

    Article  CAS  PubMed  Google Scholar 

  35. Ulianov S, Doronin S, Khrameeva E, Kos P, Luzhin A, Starikov S, Galitsyna A, Nenasheva V, Ilyin A, Flyamer I, Mikhaleva E, Logacheva M, Gelfand M, Chertovich A, Gavrilov A, Razin S, Shevelyov Y. Nuclear lamina integrity is required for proper spatial organization of chromatin in Drosophila. Nat Commun. 2019;10(1):1176. https://doi.org/10.1038/s41467-019-09185-y.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Padeken J, Mendiburo MJ, Chlamydas S, Schwarz H-J, Kremmer E, Heun P. The nucleoplasmin homolog NLP mediates centromere clustering and anchoring to the nucleolus. Mol Cell. 2013;50(2):236–49. https://doi.org/10.1016/j.molcel.2013.03.002.

    Article  CAS  PubMed  Google Scholar 

  37. Prokocimer M, Davidovich M, Nissim-Rafinia M, Wiesel-Motiuk N, Bar DZ, Barkan R, Meshorer E, Gruenbaum Y. Nuclear lamins: key regulators of nuclear structure and activities. J Cell Mol Med. 2009;13(6):1059–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. de Leeuw R, Gruenbaum Y, Medalia O. Nuclear lamins: thin filaments with major functions. Trends Cell Biol. 2018;28(1):34–45. https://doi.org/10.1016/j.tcb.2017.08.004.

    Article  CAS  PubMed  Google Scholar 

  39. Marshall WF. Order and disorder in the nucleus. Curr Biol. 2002;12(5):185–92. https://doi.org/10.1016/S0960-9822(02)00724-8.

    Article  Google Scholar 

  40. Bronshtein I, Kepten E, Kanter I, Berezin S, Lindner M, Redwood AB, Mai S, Gonzalo S, Foisner R, Shav-Tal Y, Garini Y. Loss of lamin A function increases chromatin dynamics in the nuclear interior. Nat Commun. 2015;6(1):8044. https://doi.org/10.1038/ncomms9044.

    Article  CAS  PubMed  Google Scholar 

  41. Spichal M, Fabre E. The emerging role of the cytoskeleton in chromosome dynamics. Front Genet. 2017;8:60. https://doi.org/10.3389/fgene.2017.00060.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Caridi CP, D’Agostino C, Ryu T, Zapotoczny G, Delabaere L, Li X, Khodaverdian VY, Amaral N, Lin E, Rau AR, et al. Nuclear F-actin and myosins drive relocalization of heterochromatic breaks. Nature. 2018;559(7712):54–60. https://doi.org/10.1038/s41586-018-0242-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Pickersgill H, Kalverda B, De Wit E, Talhout W, Fornerod M, van Steensel B. Characterization of the Drosophila melanogaster genome at the nuclear lamina. Nat Genet. 2006;38(9):1005–14.

    Article  CAS  PubMed  Google Scholar 

  44. Guelen L, Pagie L, Brasset E, Meuleman W, Faza MB, Talhout W, Eussen BH, de Klein A, Wessels L, de Laat W, et al. Domain organization of human chromosomes revealed by mapping of nuclear lamina interactions. Nature. 2008;453(7197):948–51.

    Article  CAS  PubMed  Google Scholar 

  45. Peric-Hupkes D, Meuleman W, Pagie L, Bruggeman SW, Solovei I, Brugman W, Gräf S, Flicek P, Kerkhoven RM, van Lohuizen M, et al. Molecular maps of the reorganization of genome-nuclear lamina interactions during differentiation. Mol Cell. 2010;38(4):603–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Shevelyov YY, Ulianov SV. The nuclear lamina as an organizer of chromosome architecture. Cells. 2019;8(2):136. https://doi.org/10.3390/cells8020136.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. van Bemmel JG, Pagie L, Braunschweig U, Brugman W, Meuleman W, Kerkhoven RM, van Steensel B. The insulator protein su(hw) fine-tunes nuclear lamina interactions of the Drosophila genome. PloS ONE. 2010;5(11):15013.

    Article  Google Scholar 

  48. Steensel B, Belmont A. Lamina-associated domains: links with chromosome architecture, heterochromatin, and gene repression. Cell. 2017;169:780–91. https://doi.org/10.1016/j.cell.2017.04.022.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Briand N, Collas P. Lamina-associated domains: peripheral matters and internal affairs. Genome Biol. 2020;21:85.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Manzo SG, Dauban L, van Steensel B. Lamina-associated domains: tethers and looseners. Curr Opin Cell Biol. 2022;74:80–7. https://doi.org/10.1016/j.ceb.2022.01.004.

    Article  CAS  PubMed  Google Scholar 

  51. Kind J, Pagie L, de Vries SS, Nahidiazar L, Dey SS, Bienko M, Zhan Y, Lajoie B, de Graaf CA, Amendola M, Fudenberg G, Imakaev M, Mirny LA, Jalink K, Dekker J, van Oudenaarden A, van Steensel B. Genome-wide maps of nuclear lamina interactions in single human cells. Cell. 2015;163(1):134–47. https://doi.org/10.1016/j.cell.2015.08.040.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Chiang M, Michieletto D, Brackley CA, Rattanavirotkul N, Mohammed H, Marenduzzo D, Chandra T. Polymer modeling predicts chromosome reorganization in senescence. Cell Rep. 2019;28(12):3212–32236. https://doi.org/10.1016/j.celrep.2019.08.045.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Gonzalez-Sandoval A, Gasser SM. On tads and lads: spatial control over gene expression. Trends Genet. 2016;32(8):485–95.

    Article  CAS  PubMed  Google Scholar 

  54. Bajpai G, Amiad Pavlov D, Lorber D, Volk T, Safran S. Mesoscale phase separation of chromatin in the nucleus. eLife. 2021;10:63976. https://doi.org/10.7554/eLife.63976.

    Article  Google Scholar 

  55. Falk M, Feodorova Y, Naumova N, Imakaev M, Lajoie BR, Leonhardt H, Joffe B, Dekker J, Fudenberg G, Solovei I, Mirny LA. Heterochromatin drives compartmentalization of inverted and conventional nuclei. Nature. 2019;570(7761):395–9. https://doi.org/10.1038/s41586-019-1275-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Sati S, Bonev B, Szabo Q, Jost D, Bensadoun P, Serra F, Loubiere V, Papadopoulos GL, Rivera-Mulia J-C, Fritsch L, Bouret P, Castillo D, Gelpi JL, Orozco M, Vaillant C, Pellestor F, Bantignies F, Marti-Renom MA, Gilbert DM, Lemaitre J-M, Cavalli G. 4D genome rewiring during oncogene-induced and replicative senescence. Mol Cell. 2020;78(3):522–5389. https://doi.org/10.1016/j.molcel.2020.03.007.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Maji A, Ahmed JA, Roy S, Chakrabarti B, Mitra MK. A lamin-associated chromatin model for chromosome organization. Biophys J. 2020;118(12):3041–50. https://doi.org/10.1016/j.bpj.2020.05.014.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Laghmach R, Di Pierro M, Potoyan DA. Mesoscale liquid model of chromatin recapitulates nuclear order of eukaryotes. Biophys J. 2020;118(9):2130–40. https://doi.org/10.1016/j.bpj.2019.09.013.

    Article  CAS  PubMed  Google Scholar 

  59. Zhou R, Gao YQ. A DNA sequence based polymer model for chromatin folding. Int J Mol Sci. 2021. https://doi.org/10.3390/ijms22031328.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Brunet A, Destainville N, Collas P. Physical constraints in polymer modeling of chromatin associations with the nuclear periphery at kilobase scale. Nucleus. 2021;12(1):6–20. https://doi.org/10.1080/19491034.2020.1868105.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Lin W-Q, Ngian Z-K, Koh T-W, Ong C-T. Altered stability of nuclear lamin-B marks the onset of aging in male Drosophila. PLoS ONE. 2022;17(3):1–12. https://doi.org/10.1371/journal.pone.0265223.

    Article  CAS  Google Scholar 

  62. Nazer E. To be or not be (in the LAD): emerging roles of lamin proteins in transcriptional regulation. Biochem Soc Trans. 2022;50(2):1035–44. https://doi.org/10.1042/BST20210858.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Panatta E, Butera A, Celardo I, Leist M, Melino G, Amelio I. p53 regulates expression of nuclear envelope components in cancer cells. Biology Direct. 2022;17(12):38. https://doi.org/10.1186/s13062-022-00349-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Ayala-Guerrero L, Claudio-Galeana S, Furlan-Magaril M, Castro-Obregón S. Chromatin structure from development to ageing. In: Harris JR, Korolchuk VI, editors. Biochemistry and cell biology of ageing: part III biomedical science. Springer: New York; 2023. p. 7–51.

    Chapter  Google Scholar 

  65. Lukyanchikova V, Nuriddinov M, Belokopytova P, Taskina A, Liang J, Reijnders MJMF, Ruzzante L, Feron R, Waterhouse R, Wu Y, Mao C, Tu Z, Sharakhov I, Fishman V. Anopheles mosquitoes reveal new principles of 3D genome organization in insects. Nat Commun. 2022. https://doi.org/10.1038/s41467-022-29599-5.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Boettiger AN, Bintu B, Moffitt JR, Wang S, Beliveau BJ, Fudenberg G, Imakaev M, Mirny LA, Wu C-T, Zhuang X. Super-resolution imaging reveals distinct chromatin folding for different epigenetic states. Nature. 2016;529(7586):418–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Rosa A, Everaers R. Structure and dynamics of interphase chromosomes. PLoS Comput Biol. 2008;4(8):1000153. https://doi.org/10.1371/journal.pcbi.1000153.

    Article  CAS  Google Scholar 

  68. Barbieri M, Chotalia M, Fraser J, Lavitas L-M, Dostie J, Pombo A, Nicodemi M. Complexity of chromatin folding is captured by the strings and binders switch model. Proc Natl Acad Sci. 2012;109(40):16173–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Tjong H, Gong K, Chen L, Alber F. Physical tethering and volume exclusion determine higher-order genome organization in budding yeast. Genome Res. 2012;22(7):1295–305. https://doi.org/10.1101/gr.129437.111.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Tokuda N, Terada T, Sasai M. Dynamical modeling of three-dimensional genome organization in interphase budding yeast. Biophys J. 2012;102(2):296–304. https://doi.org/10.1016/j.bpj.2011.12.005.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Wong H, Marie-Nelly H, Herbert S, Carrivain P, Blanc H, Koszul R, Fabre E, Zimmer C. A predictive computational model of the dynamic 3D interphase yeast nucleus. Curr Biol. 2012;22(20):1881–90. https://doi.org/10.1016/j.cub.2012.07.069.

    Article  CAS  PubMed  Google Scholar 

  72. Kalhor R, Tjong H, Jayathilaka N, Alber F, Chen L. Genome architectures revealed by tethered chromosome conformation capture and population-based modeling. Nat Biotechnol. 2012;30(1):90–8.

    Article  CAS  Google Scholar 

  73. Jerabek H, Heermann D. How chromatin looping and nuclear envelope attachment affect genome organization in eukaryotic cell nuclei. In: Hancock R, Jeon KW, editors. New models of the cell nucleus: crowding, entropic forces, phase separation, and fractals., vol. 307. International review of cell and molecular biology. Cambridge: Academic Press; 2014. p. 351–81. https://doi.org/10.1016/B978-0-12-800046-5.00010-2.

    Chapter  Google Scholar 

  74. Bruinsma R, Grosberg A, Rabin Y, Zidovska A. Chromatin hydrodynamics. Biophys J. 2014;106(9):1871–81. https://doi.org/10.1016/j.bpj.2014.03.038.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  75. Sanborn A, Rao S, Huang S-C, Durand N, Huntley M, Jewett A, Bochkov I, Chinnappan D, Cutkosky A, Li J, Geeting K, Gnirke A, Melnikov A, McKenna D, Stamenova E, Lander E, Aiden E. Chromatin extrusion explains key features of loop and domain formation in wild-type and engineered genomes. Proc Natl Acad Sci. 2015;112(47):6456–65. https://doi.org/10.1073/pnas.1518552112.

    Article  CAS  Google Scholar 

  76. Zhang B, Wolynes PG. Topology, structures, and energy landscapes of human chromosomes. Proc Natl Acad Sci. 2015;112(19):6062–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Kinney N, Onufriev A, Sharakhov I. Quantified effects of chromosome-nuclear envelope attachments on 3D organization of chromosomes. Nucleus. 2015;6(3):212–24. https://doi.org/10.1080/19491034.2015.1056441.

    Article  PubMed  PubMed Central  Google Scholar 

  78. Tjong H, Li W, Kalhor R, Dai C, Hao S, Gong K, Zhou Y, Li H, Zhou XJ, Le Gros MA, et al. Population-based 3D genome structure analysis reveals driving forces in spatial genome organization. Proc Natl Acad Sci. 2016;113(12):1663–72.

    Article  Google Scholar 

  79. Di Pierro M, Cheng RR, Lieberman-Aiden E, Wolynes PG, Onuchic JN. De novo prediction of human chromosome structures: epigenetic marking patterns encode genome architecture. Proc Natl Acad Sci. 2017;114(46):12126–31.

    Article  PubMed  PubMed Central  Google Scholar 

  80. Serra F, Baù D, Goodstadt M, Castillo D, Filion GJ, Marti-Renom MA. Automatic analysis and 3D-modelling of Hi-C data using TADbit reveals structural features of the fly chromatin colors. PLoS Comput Biol. 2017;13(7):1005665.

    Article  Google Scholar 

  81. Gürsoy G, Xu Y, Kenter AL, Liang J. Computational construction of 3D chromatin ensembles and prediction of functional interactions of alpha-globin locus from 5C data. Nucleic Acids Res. 2017;45(20):11547–58. https://doi.org/10.1093/nar/gkx784.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. Haddad N, Jost D, Vaillant C. Perspectives: using polymer modeling to understand the formation and function of nuclear compartments. Chromosome Res. 2017;25(1):35–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  83. Kinney NA, Sharakhov IV, Onufriev AV. Chromosome-nuclear envelope attachments affect interphase chromosome territories and entanglement. Epigenet Chromatin. 2018;11:3.

    Article  Google Scholar 

  84. MacPherson Q, Beltran B, Spakowitz AJ. Bottom-up modeling of chromatin segregation due to epigenetic modifications. Proc Natl Acad Sci. 2018;115(50):12739–44. https://doi.org/10.1073/pnas.1812268115.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Shi G, Liu L, Hyeon C, Thirumalai D. Interphase human chromosome exhibits out of equilibrium glassy dynamics. Nat Commun. 2018;9(1):3161. https://doi.org/10.1038/s41467-018-05606-6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  86. Ghosh SK, Jost D. How epigenome drives chromatin folding and dynamics, insights from efficient coarse-grained models of chromosomes. PLoS Comput Biol. 2018;14(5):1006159.

    Article  Google Scholar 

  87. Paulsen J, Ali TML, Collas P. Computational 3D genome modeling using chrom3D. Nat Protoc. 2018;13(5):1137–52.

    Article  CAS  PubMed  Google Scholar 

  88. MacKay K, Kusalik A. Computational methods for predicting 3D genomic organization from high-resolution chromosome conformation capture data. Brief Funct Genomics. 2020;19(4):292–308.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  89. Liang J, Perez-Rathke A. Minimalistic 3D chromatin models: sparse interactions in single cells drive the chromatin fold and form many-body units. Curr Opin Struct Biol. 2021;71:200–14. https://doi.org/10.1016/j.sbi.2021.06.017.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  90. Di Stefano M, Nützmann H-W, Marti-Renom MA, Jost D. Polymer modelling unveils the roles of heterochromatin and nucleolar organizing regions in shaping 3D genome organization in Arabidopsis thaliana. Nucleic Acids Res. 2021;49(4):1840–58. https://doi.org/10.1093/nar/gkaa1275.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  91. Ulianov SV, Zakharova VV, Galitsyna AA, Kos PI, Polovnikov KE, Flyamer IM, Mikhaleva EA, Khrameeva EE, Germini D, Logacheva MD, et al. Order and stochasticity in the folding of individual Drosophila genomes. Nat Commun. 2021;12(1):41. https://doi.org/10.1038/s41467-020-20292-z.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  92. Esposito A, Bianco S, Chiariello AM, Abraham A, Fiorillo L, Conte M, Campanile R, Nicodemi M. Polymer physics reveals a combinatorial code linking 3D chromatin architecture to 1D chromatin states. Cell Rep. 2022;38(13):110601. https://doi.org/10.1016/j.celrep.2022.110601.

    Article  CAS  PubMed  Google Scholar 

  93. Di Stefano M, Paulsen J, Jost D, Marti-Renom MA. 4D Nucleome modeling. Curr Opin Genet Dev. 2021;67:25–32. https://doi.org/10.1016/j.gde.2020.10.004.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  94. Vivante A, Brozgol E, Bronshtein I, Levi V, Garini Y. Chromatin dynamics governed by a set of nuclear structural proteins. Genes Chromosomes Cancer. 2019;58(7):437–51. https://doi.org/10.1002/gcc.22719.

    Article  CAS  PubMed  Google Scholar 

  95. Nagano T, Lubling Y, Várnai C, Dudley C, Leung W, Baran Y, Cohen NM, Wingett S, Fraser P, Tanay A. Cell-cycle dynamics of chromosomal organization at single-cell resolution. Nature. 2017;547(7661):61–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  96. Maass PG, Barutcu AR, Weiner CL, Rinn JL. Inter-chromosomal contact properties in live-cell imaging and in Hi-C. Mol Cell. 2018;69(6):1039–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  97. Gibcus JH, Samejima K, Goloborodko A, Samejima I, Naumova N, Nuebler J, Kanemaki MT, Xie L, Paulson JR, Earnshaw WC, Mirny LA, Dekker J. A pathway for mitotic chromosome formation. Science. 2018;359(6376):6135. https://doi.org/10.1126/science.aao6135.

    Article  CAS  Google Scholar 

  98. Shinkai S, Nakagawa M, Sugawara T, Togashi Y, Ochiai H, Nakato R, Taniguchi Y, Onami S. Phi-C: deciphering Hi-C data into polymer dynamics. NAR Genom Bioinform. 2020;2(2):lqaa020.

    Article  PubMed  PubMed Central  Google Scholar 

  99. Zuo W, Chen G, Gao Z, Li S, Chen Y, Huang C, Chen J, Chen Z, Lei M, Bian Q. Stage-resolved Hi-C analyses reveal meiotic chromosome organizational features influencing homolog alignment. Nat Commun. 2021;12:5827. https://doi.org/10.1038/s41467-021-26033-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  100. Feig M, Sugita Y. Whole-cell models and simulations in molecular detail. Annu Rev Cell Dev Biol. 2019;35(1):191–211. https://doi.org/10.1146/annurev-cellbio-100617-062542. (PMID: 31299173).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  101. Brackley CA, Johnson J, Kelly S, Cook PR, Marenduzzo D. Simulated binding of transcription factors to active and inactive regions folds human chromosomes into loops, rosettes and topological domains. Nucleic Acids Res. 2016;44(8):3503–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  102. Fudenberg G, Imakaev M, Lu C, Goloborodko A, Abdennur N, Mirny LA. Formation of chromosomal domains by loop extrusion. Cell Rep. 2016;15(9):2038–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  103. Lifshitz IM, Grosberg AY, Khokhlov AR. Some problems of the statistical physics of polymer chains with volume interaction. Rev Mod Phys. 1978;50:683–713. https://doi.org/10.1103/RevModPhys.50.683.

    Article  CAS  Google Scholar 

  104. Kremer K, Grest G. Dynamics of entangled linear polymer melts: a molecular-dynamics simulation. J Chem Phys. 1990;92(8):5057–86. https://doi.org/10.1063/1.458541.

    Article  CAS  Google Scholar 

  105. Kim Y, Afanasyev AY, Tolokh IS, Sharakhov IV, Onufriev AV. The probability of chromatin to be at the nuclear lamina has no systematic effect on its transcription level. bioRxiv. 2022. https://doi.org/10.1101/2022.02.17.480932.

    Article  PubMed  PubMed Central  Google Scholar 

  106. Filion G, Van Bemmel J, Braunschweig U, Talhout W, Kind J, Ward L, Brugman W, Castro I, Kerkhoven R, Bussemaker H, van Steensel B. Systematic protein location mapping reveals five principal chromatin types in Drosophila cells. Cell. 2010;143:212–24. https://doi.org/10.1016/j.cell.2010.09.009.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  107. Rowley MJ, Nichols MH, Lyu X, Ando-Kuri M, Rivera ISM, Hermetz K, Wang P, Ruan Y, Corces VG. Evolutionarily conserved principles predict 3D chromatin organization. Mol Cell. 2017;67(5):837–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  108. Hochstrasser M, Mathog D, Gruenbaum Y, Saumweber H, Sedat JW. Spatial organization of chromosomes in the salivary gland nuclei of Drosophila melanogaster. J Cell Biol. 1986;102(1):112–23.

    Article  CAS  PubMed  Google Scholar 

  109. Fung JC, Marshall WF, Dernburg A, Agard DA, Sedat JW, John A, Sedat W, Biochemistry DO. Homologous chromosome pairing in Drosophila melanogaster proceeds through multiple independent initiations. J Cell Biol. 1998;141:5–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  110. Child M, Bateman JR, Jahangiri A, Reimer A, Lammers NC, Sabouni N, Villamarin D, Mckenzie-Smith GC, Johnson JE, Jost D, Garcia H. Live imaging and biophysical modeling support a button-based mechanism of somatic homolog pairing in Drosophila. eLife. 2021;10:64412. https://doi.org/10.7554/eLife.64412.

    Article  Google Scholar 

  111. Abed JA, Erceg J, Goloborodko A, Nguyen SC, Mccole RB, Saylor W, Fudenberg G, Lajoie BR, Dekker J, Mirny LA, Wu C-T. Highly structured homolog pairing reflects functional organization of the Drosophila genome. Nat Commun. 2019;10:4485. https://doi.org/10.1038/s41467-019-12208-3.

    Article  CAS  Google Scholar 

  112. Mirny L. The fractal globule as a model of chromatin architecture in the cell. Chromosome Res. 2011;19:37–51. https://doi.org/10.1007/s10577-010-9177-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  113. Thomas JO, Kornberg RD. An octamer of histones in chromatin and free in solution. Proc Natl Acad Sci. 1975;72(7):2626–30. https://doi.org/10.1073/pnas.72.7.2626.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  114. Campbell AM, Cotter RI. The molecular weight of nucleosome protein by laser light scattering. FEBS Lett. 1976;70(1–2):209–11.

    Article  CAS  PubMed  Google Scholar 

  115. Kornberg RD. Chromatin structure: a repeating unit of histones and DNA. Science. 1974;184(4139):868–71. https://doi.org/10.1126/science.184.4139.868.

    Article  CAS  PubMed  Google Scholar 

  116. Widom J. A relationship between the helical twist of DNA and the ordered positioning of nucleosomes in all eukaryotic cells. Proceed Natl Acad Sci. 1992;89(3):1095–9. https://doi.org/10.1073/pnas.89.3.1095.

    Article  CAS  Google Scholar 

  117. Routh A, Sandin S, Rhodes D. Nucleosome repeat length and linker histone stoichiometry determine chromatin fiber structure. Proc Natl Acad Sci. 2008;105(26):8872–7. https://doi.org/10.1073/pnas.0802336105.

    Article  PubMed  PubMed Central  Google Scholar 

  118. Zhao K, Harel A, Stuurman N, Guedalia D, Gruenbaum Y. Binding of matrix attachment regions to nuclear lamin is mediated by the rod domain and depends on the lamin polymerization state. FEBS Lett. 1996;380:161–4. https://doi.org/10.1016/0014-5793(96)00034-8.

    Article  CAS  PubMed  Google Scholar 

  119. Bondarenko SM, Sharakhov IV. Reorganization of the nuclear architecture in the Drosophila melanogaster Lamin B mutant lacking the CaaX box. Nucleus. 2020;11(1):283–98.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  120. Limbach H-J, Arnold A, Mann BA, Holm C. Espresso-an extensible simulation package for research on soft matter systems. Comput Phys Commun. 2006;174(9):704–27.

    Article  CAS  Google Scholar 

  121. Kremer K, Grest G. Molecular dynamics (MD) simulations for polymers. J Phys Condens Matter. 1990;2(S):295–8. https://doi.org/10.1088/0953-8984/2/s/045.

    Article  Google Scholar 

  122. Anandakrishnan R, Drozdetski A, Walker R, Onufriev A. Speed of conformational change: comparing explicit and implicit solvent molecular dynamics simulations. Biophys J. 2015;108(5):1153–64. https://doi.org/10.1016/j.bpj.2014.12.047.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  123. Weber S, Spakowitz A, Theriot J. Bacterial chromosomal loci move subdiffusively through a viscoelastic cytoplasm. Phys Rev Lett. 2010;104:238102. https://doi.org/10.1103/PhysRevLett.104.238102.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  124. Tamm MV, Nazarov LI, Gavrilov AA, Chertovich AV. Anomalous diffusion in fractal globules. Phys Rev Lett. 2015;114:178102. https://doi.org/10.1103/PhysRevLett.114.178102.

    Article  CAS  PubMed  Google Scholar 

  125. Brückner DB, Chen H, Barinov L, Zoller B, Gregor T. Stochastic motion and transcriptional dynamics of pairs of distal DNA loci on a compacted chromosome. bioRxiv. 2023. https://doi.org/10.1101/2023.01.18.524527.

    Article  PubMed  PubMed Central  Google Scholar 

  126. Di Pierro M, Zhang B, Aiden EL, Wolynes PG, Onuchic JN. Transferable model for chromosome architecture. Proc Natl Acad Sci. 2016;113(43):12168–73. https://doi.org/10.1073/pnas.1613607113.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  127. Collas P, Liyakat Ali TM, Brunet A, Germier T. Finding friends in the crowd: three-dimensional cliques of topological genomic domains. Front Genet. 2019;10:602. https://doi.org/10.3389/fgene.2019.00602.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  128. Paulsen J, Ali TML, Nekrasov M, Delbarre E, Baudement M-O, Kurscheid S, Tremethick D, Collas P. Long-range interactions between topologically associating domains shape the four-dimensional genome during differentiation. Nat Genet. 2019;51(5):835–43.

    Article  CAS  PubMed  Google Scholar 

  129. Flory PJ. Thermodynamics of heterogeneous polymers and their solutions. J Chem Phys. 1944;12(11):425–38.

    Article  CAS  Google Scholar 

  130. Schlosshauer M, Baker D. Realistic protein-protein association rates from a simple diffusional model neglecting long-range interactions, free energy barriers, and landscape ruggedness. Protein Sci. 2004;13(6):1660–9. https://doi.org/10.1110/ps.03517304.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  131. Yang T, Zhang F, Yardımcı GG, Song F, Hardison RC, Noble WS, Yue F, Li Q. HiCRep: assessing the reproducibility of Hi-C data using a stratum-adjusted correlation coefficient. Genome Res. 2017;27(11):1939–49. https://doi.org/10.1101/gr.220640.117.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  132. Dolfini S, Courgeon A, Tiepolo L. The cell cycle of an established line of Drosophila melanogaster cells in vitro. Experientia. 1970;26(9):1020–1.

    Article  CAS  PubMed  Google Scholar 

  133. Gatti M, Tanzarella C, Olivieri G. Analysis of the chromosome aberrations induced by X-rays in somatic cells of Drosophila melanogaster. Genetics. 1974;77(4):701–19.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  134. Ashburner M. Drosophila: a laboratory manual. Cold Spring Harbor, NY: Cold Spring Harbor Laboratory. 1989. https://books.google.com/books?id=FqfwAAAAMAAJ.

  135. Dernburg AF, Broman KW, Fung JC, Marshall WF, Philips J, Agard DA, Sedat JW. Perturbation of nuclear architecture by long-distance chromosome interactions. Cell. 1996;85(5):745–59.

    Article  CAS  PubMed  Google Scholar 

  136. Avşaroğlu B, Bronk G, Gordon-Messer S, Ham J, Bressan DA, Haber JE, Kondev J. Effect of chromosome tethering on nuclear organization in yeast. PLoS ONE. 2014;9(7):1–11. https://doi.org/10.1371/journal.pone.0102474.

    Article  Google Scholar 

  137. Camps J, Erdos MR, Ried T. The role of lamin B1 for the maintenance of nuclear structure and function. Nucleus. 2015;6(1):8–14. https://doi.org/10.1080/19491034.2014.1003510.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  138. Forsberg F, Brunet A, Ali TML, Collas P. Interplay of lamin a and lamin B lads on the radial positioning of chromatin. Nucleus. 2019;10(1):7–20. https://doi.org/10.1080/19491034.2019.1570810.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  139. Wang Q, Sun Q, Czajkowsky DM, Shao Z. Sub-kb Hi-C in D. melanogaster reveals conserved characteristics of TADs between insect and mammalian cells. Nat Commun. 2018;9(1):188. https://doi.org/10.1038/s41467-017-02526-9.

  140. Leemans C, Zwalm MVD, Brueckner L, Comoglio F, Schaik TV, Pagie L, Arensbergen JV, Steensel BV. Promoter-intrinsic and local chromatin features determine gene repression in lamina-associated domains. Cell. 2018;177:852–64. https://doi.org/10.1101/464081.

    Article  CAS  Google Scholar 

  141. Meuleman W, Peric-Hupkes D, Kind J, Beaudry J-B, Pagie L, Kellis M, Reinders M, Wessels L, van Steensel B. Constitutive nuclear lamina-genome interactions are highly conserved and associated with A/T-rich sequence. Genome Res. 2013;23(2):270–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  142. Rullens PMJ, Kind J. Attach and stretch: emerging roles for genome-lamina contacts in shaping the 3D genome. Curr Opin Cell Biol. 2021;70:51–7. https://doi.org/10.1016/j.ceb.2020.11.006.

    Article  CAS  PubMed  Google Scholar 

  143. Paulsen J, Sekelja M, Oldenburg AR, Barateau A, Briand N, Delbarre E, Shah A, Sørensen AL, Vigouroux C, Buendia B, et al. Chrom3D: three-dimensional genome modeling from Hi-C and nuclear lamin-genome contacts. Genome Biol. 2017;18(1):1–15.

    Article  CAS  Google Scholar 

  144. Lin W-Q, Ngian Z-K, Koh T-W, Ong C-T. Altered stability of nuclear lamin-B marks the onset of aging in male Drosophila. PLoS ONE. 2022;17(3):1–12. https://doi.org/10.1371/journal.pone.0265223.

    Article  CAS  Google Scholar 

  145. Egli D, Birkhoff G, Eggan K. Mediators of reprogramming: transcription factors and transitions through mitosis. Nat Rev Mol Cell biol. 2008;9:505–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  146. Güttinger S, Laurell E, Kutay U. Orchestrating nuclear envelope disassembly and reassembly during mitosis. Nat Rev Mol Cell Biol. 2009;10(3):178–91. https://doi.org/10.1038/nrm2641.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Raju Nadimpalli for help with ESPResSO input scripts, chromatin visualization and Hi-C map processing; Samira Mali for technical assistance and for producing the movie clip of the chromosomes dynamics; Yoonjin Kim for technical assistance; Frank Alber for a useful discussion and for providing a comprehensive set of configuration snapshots from Ref. [26].

Funding

This study was supported by the National Science Foundation [MCB-1715207], and, in part, by the National Institutes of Health [R01 GM144596] (to A.V.O).

Author information

Authors and Affiliations

Authors

Contributions

IT and AO designed the dynamical model with input from IS and NK. IT implemented the model with early input from NK. IT, AO and IS analyzed the results. IT, AO, IS and NK wrote the original manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Igor V. Sharakhov or Alexey V. Onufriev.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

All authors read and approved the manuscript.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1:

Figure S1. The number of TADs and LAD containing TADs (L-TADs) in each of the four epigenetic classes of TADs (Active, Null, PcG and HP1) are shown: the numbers of TADs – as solid color bars, the numbers of L-TADs – as bars with diagonal stripes. Figure S2. Determining the interaction strength between the special pairs of TADs. Shown is the logarithm of the ordered (index n) degree of enrichment for the contact probabilities of 268 specific pairs of remote chromatin loci [6] and approximating fit (Eq. S5) used to determine the attractive interaction LJ well depth for the interactions between the beads int hat specific pairs. Figure S3. Time dependence of mean squared displacement (MSD) of bead averaged over 9 beads and 18 trajectories of the model nuclei with different topologies and size, simulated with 3 different values of the friction parameter γ. The simulation data points (symbols with error bars) are fitted with Eq. S7 (dashed lines)– the equation corresponds to theexperimentally observed time dependence of MSD for a chromosomal loci [14]. The fit is used to estimate the simulation time scaling factor λ, which determines the mapping between the simulation time and biological time of a nucleus. Error bars in each curve are RMSD of 9-bead-average values in 18 trajectories. Figure S4. Cumulative normalized distributions of L-TADs (LADcontaining TADs) at different levels of attractive interaction with the nuclear envelope (NE) varying from 0.1 kT to 5 kT. The vertical dashed line at R=1.91 μm marks the center of the bead layer nearest to the NE. Experimental data [15] show that about 25% (horizontal dashed line at 0.75 mark) of LADs are on average located at the NE. Figure S5. Details of the temporal evolution of the WT model derived TAD-TAD contact probability (HiC) maps. The Hi-C maps are averaged over 5 min time intervals and over the ensemble of 18 system trajectories (see "Methods"). The top left panel is t=0 min Hi-C map, followed (from left to right) by the difference (from t=0 min) maps at 30 min, 3 h and 11 h. Figure S6. The effect of complete removal of the NE: chromatin de-compacts. Within 1 min, chromosome 4 dissociates and drifts away from the rest three de-condensed chromosomes. Figure S7. Top panel: Chromatin density distributions in the model nuclei at different levels of LAD-NE attractive interaction (from 0.1 to 5 kT). Bottom panel: Chromatin density distributions in the model nuclei at different levels of attractive interaction between Null TADs. Figure S8. Difference Hi-C map between Lamin mutant model Hi-C map (Fig.5, main text) and the corresponding WT model Hi-C map (Fig. 3, main text, bottom panel). Note that the intensity scale is 10 times smaller here than that of the original Hi-C map, pointing to relatively small differences in TAD-TAD contact probabilities. Most TAD-TAD contacts, including inter-chromosome contacts, are slightly enhanced (red areas) in the Lamin mutant compared to the WT nuclei. Small areas where the contact frequency decreases (blue spots) in the Lamin mutant are limited to close to diagonal intra-arm contacts. Figure S9. Increasing the cross-type TAD-TAD interactions from 0.5 kT (model selected) to 1.0 kT results in increased TAD-TAD contact probabilities, both in WT and Lamin mutant model nuclei. The Hi-C map differences between modified and selected model nuclei are shown. (Top panel:) for the WT model nuclei. (Bottom panel:) for the Lamin mutant model nuclei. Figure S10. The mobility of five select L-TADs over about two minutes time interval of the interphase, starting at t = 0. Shown are four snapshots from the supplementary movie that represents 20 minutes of time-evolution of model fruit fly nucleus. Each L-TAD is selected from a different chromosome. The L-TADs represented by orange, green, grey, red and blue are from Chr 2L, 2R, 3L, 3R and Chr X, respectively. All other beads, including the large heterochromatic beads, are shown as transparent grey spheres. The upper panel snapshots, from left to right, correspond t = 0.04(s) and t = 0.28(s), respectively. The lower panel snapshots, from left to right, correspond to t = 35.5(s) and t = 104(s), respectively. The coordinate frame vectors are shown in the left corner of each panel. Apparent variation in each L-TAD sphere size conveys the depth perception: spheres closer to the viewer are larger. Image Credit: Samira Mali. Rendering by VMD [17]. Visualization of motion of these L-TADs, from 3 different viewing angles, is available as supplementary movies (Additional files 2, 3 and 4). Figure S11. The sums of contact probabilities of Null TADs (B-type TADs) with Active TADs (A-type TADs) only, calculated from the model Hi-C maps for the WT nuclei (red circles) and Lamin mutant nuclei (green circles), for each Null TAD. The average relative increase of the contact probabilities for Lamin mutant Null TADs is 22%. Figure S12. Predicted cumulative frequencies of radial positions of three cytological regions for which the corresponding experimental data is available. Top panel: Cumulative frequencies of radial positions of the Null L-TAD #15 (cytological region 22A in Ref. [16]) in the WT and Lamin mutant nucleus models. Middle panel: Cumulative frequencies of radial positions of the Null L-TAD #120 (cytological region 36C in Ref. [16]) in the WT and Lamin mutant model nuclei. Bottom panel: Cumulative frequencies of radial positions of the PcG L-TAD #435 (cytological region60D in Ref. [16]) in the WT and Lamin mutant nucleus models.

Additional file 2: The mobility of five select L-TADs over about 20 minutes of the interphase. Each L-TAD is selected from a different chromosome. The L-TADs represented by orange, green, grey, red and blue are from Chr2L, 2R, 3L, 3R and Chr X, respectively. All other beads, including the large heterochromatic beads, are shown astransparent grey spheres. Apparent variation in each L-TAD sphere size conveys the depth perception: spherescloser to the viewer are larger. The coordinate frame vectors are shown in the left corner. Image Credit: Samira Mali. Rendering by VMD.

Additional file 3: The mobility of five select L-TADs over about 20 minutes of the interphase. Each L-TAD is selected from a different chromosome. The L-TADs represented by orange, green, grey, red and blue are from Chr 2L, 2R, 3L, 3R and Chr X, respectively. All other beads, including the large heterochromatic beads, are shown as transparent grey spheres. Apparent variation in each L-TAD sphere size conveys the depth perception: spheres closer to the viewer are larger. The coordinate frame vectors are shown in the left corner. Image Credit: Samira Mali. Rendering by VMD.

Additional file 4: The mobility of five select L-TADs over about 20 minutes of the interphase. Each L-TAD is selected from a different chromosome. The L-TADs represented by orange, green, grey, red and blue are from Chr2L, 2R, 3L, 3R and Chr X, respectively. All other beads, including the large heterochromatic beads, are shown astransparent grey spheres. Apparent variation in each L-TAD sphere size conveys the depth perception: spheres closer to the viewer are larger. The coordinate frame vectors are shown in the left corner. Image Credit: Samira Mali. Rendering by VMD.

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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tolokh, I.S., Kinney, N.A., Sharakhov, I.V. et al. Strong interactions between highly dynamic lamina-associated domains and the nuclear envelope stabilize the 3D architecture of Drosophila interphase chromatin. Epigenetics & Chromatin 16, 21 (2023). https://doi.org/10.1186/s13072-023-00492-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13072-023-00492-9

Keywords