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

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.

Chromosome-nuclear envelope attachments affect interphase chromosome territories and entanglement.

BACKGROUND: It is well recognized that the interphase chromatin of higher eukaryotes folds into non-random configurations forming territories within the nucleus. Chromosome territories have biologically significant properties, and understanding how these properties change with time during lifetime of the cell is important. Chromosome-nuclear envelope (Chr-NE) interactions play a role in epigenetic regulation of DNA replication, repair, and transcription. However, their role in maintaining chromosome territories remains unclear. RESULTS: We use coarse-grained molecular dynamics simulations to study the effects of Chr-NE interactions on the dynamics of chromosomes within a model of the Drosophila melanogaster regular (non-polytene) interphase nucleus, on timescales comparable to the duration of interphase. The model simulates the dynamics of chromosomes bounded by the NE. Initially, the chromosomes in the model are prearranged in fractal-like configurations with physical parameters such as nucleus size and chromosome persistence length taken directly from experiment. Time evolution of several key observables that characterize the chromosomes is quantified during each simulation: chromosome territories, chromosome entanglement, compactness, and presence of the Rabl (polarized) chromosome arrangement. We find that Chr-NE interactions help maintain chromosome territories by slowing down and limiting, but not eliminating, chromosome entanglement on biologically relevant timescales. At the same time, Chr-NE interactions have little effect on the Rabl chromosome arrangement as well as on how chromosome compactness changes with time. These results are rationalized by simple dimensionality arguments, robust to model details. All results are robust to the simulated activity of topoisomerase, which may be present in the interphase cell nucleus. CONCLUSIONS: Our study demonstrates that Chr-NE attachments may help maintain chromosome territories, while slowing down and limiting chromosome entanglement on biologically relevant timescales. However, Chr-NE attachments have little effect on chromosome compactness or the Rabl chromosome arrangement.

Quantified effects of chromosome-nuclear envelope attachments on 3D organization of chromosomes.

We use a combined experimental and computational approach to study the effects of chromosome-nuclear envelope (Chr-NE) attachments on the 3D genome organization of Drosophila melanogaster (fruit fly) salivary gland nuclei. We consider 3 distinct models: a Null model - without specific Chr-NE attachments, a 15-attachment model - with 15 previously known Chr-NE attachments, and a 48-attachment model - with 15 original and 33 recently identified Chr-NE attachments. The radial densities of chromosomes in the models are compared to the densities observed in 100 experimental images of optically sectioned salivary gland nuclei forming "z-stacks." Most of the experimental z-stacks support the Chr-NE 48-attachment model suggesting that as many as 48 chromosome loci with appreciable affinity for the NE are necessary to reproduce the experimentally observed distribution of chromosome density in fruit fly nuclei. Next, we investigate if and how the presence and the number of Chr-NE attachments affect several key characteristics of 3D genome organization: chromosome territories and gene-gene contacts. This analysis leads to novel insight about the possible role of Chr-NE attachments in regulating the genome architecture. Specifically, we find that model nuclei with more numerous Chr-NE attachments form more distinct chromosome territories and their chromosomes intertwine less frequently. Intra-chromosome and intra-arm contacts are more common in model nuclei with Chr-NE attachments compared to the Null model (no specific attachments), while inter-chromosome and inter-arm contacts are less common in nuclei with Chr-NE attachments. We demonstrate that Chr-NE attachments increase the specificity of long-range inter-chromosome and inter-arm contacts. The predicted effects of Chr-NE attachments are rationalized by intuitive volume vs. surface accessibility arguments.

Investigation of the chromosome regions with significant affinity for the nuclear envelope in fruit fly--a model based approach.

Three dimensional nuclear architecture is important for genome function, but is still poorly understood. In particular, little is known about the role of the "boundary conditions"--points of attachment between chromosomes and the nuclear envelope. We describe a method for modeling the 3D organization of the interphase nucleus, and its application to analysis of chromosome-nuclear envelope (Chr-NE) attachments of polytene (giant) chromosomes in Drosophila melanogaster salivary glands. The model represents chromosomes as self-avoiding polymer chains confined within the nucleus; parameters of the model are taken directly from experiment, no fitting parameters are introduced. Methods are developed to objectively quantify chromosome territories and intertwining, which are discussed in the context of corresponding experimental observations. In particular, a mathematically rigorous definition of a territory based on convex hull is proposed. The self-avoiding polymer model is used to re-analyze previous experimental data; the analysis suggests 33 additional Chr-NE attachments in addition to the 15 already explored Chr-NE attachments. Most of these new Chr-NE attachments correspond to intercalary heterochromatin--gene poor, dark staining, late replicating regions of the genome; however, three correspond to euchromatin--gene rich, light staining, early replicating regions of the genome. The analysis also suggests 5 regions of anti-contact, characterized by aversion for the NE, only two of these correspond to euchromatin. This composition of chromatin suggests that heterochromatin may not be necessary or sufficient for the formation of a Chr-NE attachment. To the extent that the proposed model represents reality, the confinement of the polytene chromosomes in a spherical nucleus alone does not favor the positioning of specific chromosome regions at the NE as seen in experiment; consequently, the 15 experimentally known Chr-NE attachment positions do not appear to arise due to non-specific (entropic) forces. Robustness of the key conclusions to model assumptions is thoroughly checked.