The 3D genome is characterized by a complex organization made of genomic and epigenomic layers with profound implications on gene regulation and cell function. However, the understanding of the fundamental mechanisms driving the crosstalk between nuclear architecture and (epi)genomic information is still lacking. The plant Arabidopsis thaliana is a powerful model organism to address these questions owing to its compact genome for which we have a rich collection of microscopy, chromosome conformation capture (Hi-C) and ChIP-seq experiments. Using polymer modelling, we investigate the roles of nucleolus formation and epigenomics-driven interactions in shaping the 3D genome of A. thaliana. By validation of several predictions with published data, we demonstrate that self-attracting nucleolar organizing regions and repulsive constitutive heterochromatin are major mechanisms to regulate the organization of chromosomes. Simulations also suggest that interphase chromosomes maintain a partial structural memory of the V-shapes, typical of (sub)metacentric chromosomes in anaphase. Additionally, self-attraction between facultative heterochromatin regions facilitates the formation of Polycomb bodies hosting H3K27me3-enriched gene-clusters. Since nucleolus and heterochromatin are highly-conserved in eukaryotic cells, our findings pave the way for a comprehensive characterization of the generic principles that are likely to shape and regulate the 3D genome in many species.

Finding out how cells prepare for fate change during differentiation commitment was our task. To address whether the constitutive pericentromere-associated domains (PADs) may be involved, we used a model system with known transcriptome data, MCF-7 breast cancer cells treated with the ErbB3 ligand heregulin (HRG), which induces differentiation and is used in the therapy of cancer. PAD-repressive heterochromatin (H3K9me3), centromere-associated-protein-specific, and active euchromatin (H3K4me3) antibodies, real-time PCR, acridine orange DNA structural test (AOT), and microscopic image analysis were applied. We found a two-step DNA unfolding after 15-20 and 60 min of HRG treatment, respectively. This behavior was consistent with biphasic activation of the early response genes (c-fos – fosL1/myc) and the timing of two transcriptome avalanches reported in the literature. In control, the average number of PADs negatively correlated with their size by scale-free distribution, and centromere clustering in turn correlated with PAD size, both indicating that PADs may create and modulate a suprachromosomal network by fusing and splitting a constant proportion of the constitutive heterochromatin. By 15 min of HRG treatment, the bursting unraveling of PADs from the nucleolus boundary occurred, coinciding with the first step of H3K4me3 chromatin unfolding, confirmed by AOT. The second step after 60 min of HRG treatment was associated with transcription of long noncoding RNA from PADs and peaking of fosL1/c-myc response. We hypothesize that the bursting of PAD clusters under a critical silencing threshold pushes the first transcription avalanche, whereas the destruction of the PAD network enables genome rewiring needed for differentiation repatterning, mediated by early response bivalent genes.

The intrinsic dynamic nature of chromosomes is emerging as a fundamental component in regulating DNA transcription, replication, and damage-repair among other nuclear functions. With this increased awareness, reinforced over the last ten years, many new experimental techniques, mainly based on microscopy and chromosome conformation capture, have been introduced to study the genome in space and time. Owing to the increasing complexity of these cutting-edge techniques, computational approaches have become of paramount importance to interpret, contextualize, and complement such experiments with new insights. Hence, it is becoming crucial for experimental biologists to have a clear understanding of the diverse theoretical modeling approaches available and the biological information each of them can provide.

Bulk chromatin motion has not been analyzed at high resolution. We present Hi-D, a method to quantitatively map dynamics of chromatin and abundant nuclear proteins for every pixel simultaneously over the entire nucleus from fluorescence image series. Hi-D combines reconstruction of chromatin motion and classification of local diffusion processes by Bayesian inference. We show that DNA dynamics in the nuclear interior are spatially partitioned into 0.3–3-μm domains in a mosaic-like pattern, uncoupled from chromatin compaction. This pattern was remodeled in response to transcriptional activity. Hi-D can be applied to any dense and bulk structures opening new perspectives towards understanding motion of nuclear molecules.