Skip to main content

Advertisement

Initial high-resolution microscopic mapping of active and inactive regulatory sequences proves non-random 3D arrangements in chromatin domain clusters

Article metrics

Abstract

Background

The association of active transcription regulatory elements (TREs) with DNAse I hypersensitivity (DHS[+]) and an ‘open’ local chromatin configuration has long been known. However, the 3D topography of TREs within the nuclear landscape of individual cells in relation to their active or inactive status has remained elusive. Here, we explored the 3D nuclear topography of active and inactive TREs in the context of a recently proposed model for a functionally defined nuclear architecture, where an active and an inactive nuclear compartment (ANC–INC) form two spatially co-aligned and functionally interacting networks.

Results

Using 3D structured illumination microscopy, we performed 3D FISH with differently labeled DNA probe sets targeting either sites with DHS[+], apparently active TREs, or DHS[−] sites harboring inactive TREs. Using an in-house image analysis tool, DNA targets were quantitatively mapped on chromatin compaction shaped 3D nuclear landscapes. Our analyses present evidence for a radial 3D organization of chromatin domain clusters (CDCs) with layers of increasing chromatin compaction from the periphery to the CDC core. Segments harboring active TREs are significantly enriched at the decondensed periphery of CDCs with loops penetrating into interchromatin compartment channels, constituting the ANC. In contrast, segments lacking active TREs (DHS[−]) are enriched toward the compacted interior of CDCs (INC).

Conclusions

Our results add further evidence in support of the ANC–INC network model. The different 3D topographies of DHS[+] and DHS[−] sites suggest positional changes of TREs between the ANC and INC depending on their functional state, which might provide additional protection against an inappropriate activation. Our finding of a structural organization of CDCs based on radially arranged layers of different chromatin compaction levels indicates a complex higher-order chromatin organization beyond a dichotomic classification of chromatin into an ‘open,’ active and ‘closed,’ inactive state.

Background

Transcription regulatory elements (TREs) such as promoters, enhancers or insulators comprise non-coding sequences located within or in close vicinity to a gene but are also found up to ~1 Mb distally from their genes [1, 2]. Active TREs are characterized by a local ‘open’ chromatin conformation [3,4,5] associated with nucleosome displacement and specific histone signatures [6, 7] and were shown to have an increased sensitivity to DNAse I digestion, constituting DNAse I hypersensitive sites (DHS[+] of ~100–1000 bp [8,9,10,11]. Genome-wide profiling in various human cell types averaged over large cell populations identified >10,000 DHS[+] clusters. Approximately 10% were found to be cell type specific [2, 10, 12]. These sites encompass all experimentally validated cis-regulatory sequences; thus, DHS[+] clusters typically signify the location of active TREs in a genome.

The spatial organization of TREs has been addressed in terms of their contact frequencies with defined chromatin segments by chromosome conformation capture (Hi-C) or ChIP-Seq analyses [13,14,15,16,17], where genome-wide detection of DHS[+] sites in single cells was recently achieved by an ultrasensitive DNase sequencing strategy [18]. Still, little is known about the global 3D and 4D organization of regulatory sequences within the nuclear landscape of individual cells. Accordingly, it is unknown to which extent ‘open’ and ‘closed’ chromatin configurations may guide and constrain the accessibility of transcription factors (TFs) or chromatin modifiers to regulatory sequences. This issue has gained interest in the context of strong evidence for a distinct, but also dynamic nuclear architecture [19,20,21,22,23,24].

Microscopic investigations have demonstrated a structural organization of chromosome territories (CTs) built up from ~1-Mb chromatin domains (CDs) [25,26,27]. In particular, fluorescence labeling of specific genomic regions combined with 3D super-resolved microscopy has provided unprecedented opportunities to study nuclear arrangements of specific chromatin structures at the single-cell level and their cell-to-cell variability [28,29,30,31]. These studies have indicated that ~1-Mb CDs are composed of smaller subdomains (subCDs) and also form larger chromatin domain clusters (CDCs) [21, 26]. Genome-wide chromosome conformation capture methods have confirmed the territorial organization of chromosomes in mammalian cell nuclei and led to the discovery of ~1-Mb-sized topologically associating domains (TADs) [32,33,34,35]. TADs are built from smaller subdomains [32, 36] but also form larger units, called ‘metaTADs’ [37]. CDs and TADs reflect higher-order chromatin entities [38, 39], which provide the structural backbone for tissue-specific regulatory interactions [23, 40, 41].

In our present study, we explored the feasibility of 3D-FISH, 3D structured illumination microscopy (3D-SIM) and new 3D image analysis tools to determine the 3D nuclear topography of active and inactive TREs defined by their DHS[+] or DHS[−] status in human fibroblasts (BJ1 cells) and an adenocarcinoma cell line (A549). Below we use the terms DHS[+] interchangeably for sites with active TREs and DHS[−] for sites lacking (active) TREs. 3D FISH experiments were performed with an appropriately adapted protocol, which was previously shown to preserve key characteristics of the nuclear ultrastructure discernible at the resolution level of 3D-SIM [30, 42, 43], which is set at ~120 nm lateral and 250–300 nm axial [44, 45]. Highly resolved, quantitative measurements of DAPI intensities after DNA staining of nuclei were used as a proxy for local differences in chromatin compaction [46]. We describe the 3D topography of regulatory sequences in the context of a recently proposed model for a functionally defined nuclear landscape, where an active and an inactive nuclear compartment (termed ANC and INC) form two spatially co-aligned and functionally interacting 3D networks (for review, see [21] and Additional file 1 for illustration). The INC is formed by compact chromatin domains with low transcriptional activity for coding genes, which forms the interior core of a chromatin domain cluster (CDC). This compact core of CDCs is lined by a peripheral layer of decondensed, transcriptionally active chromatin domains, enriched in marks for transcriptionally competent chromatin termed perichromatin region (PR) [43, 47, 48]. The PR lines a contiguous channel system, the interchromatin compartment (IC), which starts at nuclear pores and permeates between CDs/CDCs [43, 48]. In addition to its potential role in nuclear import/export functions, the IC harbors nuclear bodies required for functions occurring within the PR. Accordingly, IC and PR together are considered as the active nuclear compartment (ANC).

According to the ANC–INC network model, we expected an enrichment of active TREs within the ANC, whereas a location of inactive TREs seemed possible within either the ANC or the INC. In case of inactive TREs embedded within the INC, their activation would correlate with a relocation toward the ANC. Our study argues for a non-random distinct distribution of targeted sites: Segments harboring active TREs are typically exposed at the outer periphery of CDCs constituting the active nuclear compartment (ANC), whereas inactive TREs are enriched toward the more compacted interior of CDCs, constituting the INC.

Results

Semi-automated quantitative 3D mapping of DNA sequences on chromatin compaction-defined nuclear landscapes

For our investigation, we developed a semi-automated approach for 3D mapping of FISH signals on chromatin compaction-defined nuclear landscapes [46] based on optical serial sections of DAPI-stained nuclei recorded with 3D-SIM. 3D mapping of specific hybridization signals is exemplified in Fig. 1a with two contiguous, differentially labeled 6-kb probes (see below for details and application of these probes). For an unbiased 3D assessment of probe signals, defined preset values such as minimal target size, relative signal intensity and a maximal distance between the centroid position of a given signal to its nearest differently labeled signal were applied for signal segmentation (Fig. 1b). An algorithm for chromatin compaction classification of DAPI-stained nuclei was employed to generate seven DAPI intensity classes with equal intensity variance as a measure for chromatin compaction [46]. Classes of chromatin compaction can be visualized as a color-coded heat map (Fig. 1c). Voxels assigned to class 1 (blue) depict regions with or close to background DAPI intensities, representing the largely DNA-free interchromatin compartment (IC), classes 2–3 (purple and deep red) comprise chromatin with low DAPI intensity, representing decondensed (‘open’) chromatin. Class 4 (dark orange) is considered as an intermediate zone, classes 5–7 (orange, yellow, white) comprise ‘closed’ chromatin. 3D coordinates of segmented FISH signals are spatially mapped on chromatin compaction classes (Fig. 1d) and plotted with their relative distributions on the respective classes. Figure 1e exemplifies a highly non-random distribution of 3D FISH targets within the seven chromatin compaction classes. Classes 1–3 represent the active nuclear compartment, class 4 an intermediate zone, classes 5–7 the INC (Fig. 1f, for review, see [21]).

Fig. 1
figure1

Workflow for quantitative mapping of specific FISH signals on 3D chromatin compaction maps. a Representative part of a section from an original 3D-SIM image stack of a whole nucleus acquisition shown by the example of a BJ1 nucleus. Chromatin counterstained with DAPI (gray), two adjacent DNA targets visualized by green and red fluorescent signals (arrow). The arrowhead points to an additional small green fluorescent (background) signal. Scale bar 2 µm. b Segmented fluorescent voxels of FISH signals within DAPI mask after defined parameter settings. The small isolated green signal seen in (arrowhead in a) is discarded due to a distance >0.5 µm from the nearest red signal centroid, set as limit between two differently labeled targets. c Same section after classification of DAPI signals into seven intensity classes as proxy for chromatin compaction visualized as color heat map. d Inset magnification from framed area in b and c with outlined green and red segmented signals. e Relative distribution of green and red fluorescent voxels in this nucleus mapped on DAPI intensity-defined chromatin compaction classes (gray). f Assignment of the active and inactive nuclear compartment (ANC/INC) linked to chromatin compaction classes, for explanation see text

Multilayered shell-like organization of chromatin domain clusters

Color heat maps of DAPI-stained nuclear SIM sections (as exemplified in Fig. 1c, d) suggest a multilayered shell-like chromatin organization of CDCs with compact chromatin (classes 5–7) typically located in their interior, surrounded by a decondensed peripheral layer (classes 2–3) lined by the IC (class 1). The visual impression of a radially arranged compartmentalization of chromatin layers is supported by quantitation of nearest-neighbor voxels of a given chromatin compaction class in 3D SIM serial sections of BJ1 (N = 45) and A549 (N = 30) cell nuclei (Fig. 2a). Most nearest neighbors belong to the same intensity class, a smaller fraction to the next higher or lower class and only rare voxels to remote classes. For a rough estimate of minimal distances required for potential movements of any target DNA from the most compact (interior) to the most decondensed (peripheral) part of CDCs, the distance for each voxel assigned to class 2 to the nearest voxel assigned to all other classes was measured. Minimal average distances indicate that movements of ~100 nm may suffice for a relocation of a sequence between the most inner and the most outer layer of CDCs (Fig. 2b).

Fig. 2
figure2

Quantitative assessment for a compartmentalized organization of distinct chromatin compaction classes. a Frequency of chromatin compaction class for each nearest-neighbor voxel of a given compaction class in BJ1 (N = 45) and A549 (N = 30). b Minimal average distances of each voxel of class 2 to the nearest voxel assigned to all other classes

3D topography of active and inactive TREs in chromatin domain clusters

For the coverage of DHS[+] and DHS[−] sites located on different chromosomes in BJ1 cells, we used two differently labeled DNA probe sets of pooled fosmid clones with each fosmid carrying a human sequence of ~40 kb (Fig. 3). DHS profiles of the respective sites were identified on availability of NIH Roadmap Epigenomics Mapping Consortium data (www.roadmapepigenomics.org hg19). Fosmid pool 1 comprises six genomic targets located on chromosome bands 1p33.1, 2p13.3, 2q37.3, 3p13, 5q35.3 and 12q24.21 peppered with numerous DHS[+] clusters expanding over several kb (Fig. 3a). These targets contain different types of active TREs (for type and ‘open’ chromatin marks, see Additional file 2). Fosmid pool 2 (Fig. 3b) comprises two DHS[−] segments on chromosome bands 3p22 and 13q21.31. Several TREs were identified in these segments providing evidence for regulatory potential in these regions [49, 50]. However, in line with their DHS[−] status none of these regions shows ‘open’ chromatin marks in fibroblasts (Additional file 2).

Fig. 3
figure3

DHS profiles of targeted genomic regions representing DHS[+] or DHS[−] sites on different chromosomes in BJ1 nuclei and fosmid clones used for their delineation by 3D-FISH. a Selected regions with interspersed clusters of DHS[+] sites, DHS profiles shown in black (browser shots adopted from http://encodeproject.org/). Genomic position and assignment to DHS profiles of fosmid clones used in fosmid pool 1 are indicated by green lines. b Selected regions representing DHS[−] sites, DHS profiles shown in black (browser shots adopted from http://encodeproject.org/). Genomic position and assignment to DHS profiles of fosmid clones used in fosmid pool 2 are indicated by red lines. Using pairs of (partially) overlapping clones for both fosmid pools ensures optimal hybridization efficiency. Asterisks mark the approximate location of probes described in Fig. 5

BJ1 is a diploid cell strain (46, XY, Additional file 3). Accordingly, up to twelve distinct segments carrying DHS[+] sites can be targeted in a BJ1 nucleus by fosmid pool 1 and up to four distinct DHS[−] segments by fosmid pool 2 (Additional file 4: movie_1). Representative SIM sections of a BJ1 nucleus demonstrate the preferential localization of DHS[+] targets in low chromatin compaction classes at the periphery of CDCs, and of DHS[−] targets in more interior and compacted regions (Fig. 4a). Quantitative image analysis confirms the highly significant enrichment of DHS[+] segments indicating active TREs in low chromatin compaction classes 2 and 3 and of DHS[−] segments (inactive TREs) in high compaction classes 5 and 6 (p < 0.001 for classes 2/3 and 5/6) (Fig. 4b, c). This difference is consistently seen in all single-cell profiles considered as a series of ‘snapshots’ (Additional file 5).

Fig. 4
figure4

3D nuclear topography and quantitative mapping of ~40-kb targets of DHS[+] and DHS[−] regions in BJ1 nuclei. a Part of a SIM light-optical section from a whole nucleus acquisition with framed areas indicating representative inset magnifications 1–3. DAPI-stained DNA after intensity classification shown as gray gradations and color heat maps, respectively. Segmented signals delineating targets of fosmid pool 1 (DHS[+], green) and fosmid pool 2 (DHS[−], red) show a preferential localization of pool 1 signals DHS[+] in zones of low DAPI intensity and of pool 2 (DHS[−]) within the more compacted core of CDCs as shown by outlined signals in color heat maps: pool 1 = green, pool 2 = black). Scale bar 2 µm, insets 0.5 µm. b Quantified distributions (N = 25 nuclei) of fosmid pool 1 (DHS[+], green) and pool 2 (DHS[−], red) within respective chromatin compaction classes (all classes shown in gray) confirm the significantly distinct topography for DHS[+] and DHS[−] associated signals with a shift of DHS[−] sites toward higher compaction classes. c Quantified levels of relative enrichment (positive values) or depletion (negative values) of fosmid pool 1 and pool 2 signals within chromatin compaction classes. Error bars = standard deviation of the mean.*p ≤ 0.05; **p ≤ 0.01; ***p ≤ 0.001

Both fosmid pools were also hybridized on the chromosomally rearranged adenocarcinoma cell line A549. This cell line was used for comparison with BJ1 cells for two reasons: First, its DHS profile is clearly distinct from BJ1 cells in that all genomic targets delineated by fosmid pool 1 are DHS[−] in A549 (Additional file 6) indicating the lack of active TREs as supported by a lack of ‘open’ chromatin marks (Additional file 2). Second, their flat nuclear shape with z-diameters of 4–5 µm (data not shown) facilitates 3D-SIM acquisition [51]. Numerical and structural rearrangements in A549 cells (see Additional file 7: for karyotype) permit up to 19 distinct hybridization sites for fosmid pool 1 and up to 5 sites for fosmid pool 2 (Additional file 8: movie_2). In contrast to BJ1, A549 nuclei show a fairly similar 3D nuclear topography of targets delineated by fosmid pools 1 and 2 that are both DHS[−]. A representative A549 nucleus and quantitative evaluation of 10 nuclei (as provided in Additional file 9) show the highest enrichment of signals for both fosmid pools in chromatin compaction classes 4 and 5 (for single-cell profiles, see Additional file 10; for detailed information on statistical values for all measurements, see Additional file 11).

To dissect the topography of a specific smaller segment in BJ1 cells harboring a cluster of DHS[+] peaks over its entire length, we performed 3D-FISH with a probe delineating a ~6-kb segment on chromosome 2q37 located within the 17-kb-long first intron of the COL6A3 gene, termed ‘6-kb probe 1’ (Fig. 5A). This segment contains several TREs such as an annotated transcription start site, enhancers and CTCF-binding sites (Additional file 2) with apparently high activity in BJ1 cells. The topography of this segment was compared in A549 cells where the respective sites are DHS[−] over the entire probe length (Fig. 5A). In all experiments, probe 1 was co-hybridized with a differently labeled probe which delineates an adjacent 6-kb segment with no known content of TREs [52] and accordingly DHS[−] status in both BJ1 and A549 cells, termed ‘6-kb probe 2’ (Fig. 5A). Co-hybridization of the two contiguous probes was performed for the following reasons: First, we wanted to test whether we could detect a distinct localization in different chromatin compaction classes of DHS[+] and DHS[−] sites at this length scale. Second, using a probe set composed of adjacent differently labeled probes for quantitative assessment helps to exclude in an unbiased way any dotted unspecific background signals, which can become a challenge for very small single-copy 3D FISH signals, reflecting with a signal volume of ~0.005 µm3 the volumetric resolution limit of 3D-SIM. Signals with a similar volume can arise from background fluorophores below this resolution limit. Such background cannot be entirely avoided even under most meticulous experimental conditions. Signals were accepted as true hybridization events only in case of centroid distances ≤500 nm between green and red signals, taking into account that the length of 10 kb as a fully extended 10-nm fiber is approximately 500 nm [53].

Fig. 5
figure5

Scheme and topography of two contiguous 6-kb segments targeting a DHS[+] and adjacent DHS[−] site. A Scheme of probe 1 (green; DHS[+] in BJ1, DHS[−] in A549) and probe 2 (red; DHS[−] both in BJ1 and in A549 cells) in relation to their DHS profiles (black; browser shots adopted from http://encodeproject.org/). BD SIM light-optical sections from 3D acquisitions of different BJ1 nuclei and representative inset magnifications delineating DAPI-stained DNA after intensity classification (shown as gray gradations and as color heat maps) and segmented probe signals (probe 1 green; probe 2 red). Signal positions are shown as outlines in the color heat maps (red signals outlined in black for better visibility). Magnifications reveal the localization of both red and green signals in compartments of low chromatin compaction. Note variability of signal conformation: B elongated signal pair; C signal pair with extended red signal lined by two separate green signals suggestive for ongoing replication; D two compact separate signal pairs <0.5 µm apart suggest two chromatids after replication. A second signal pair (see Additional files 12, 13, 14) is seen in a different section of each nucleus. E, F Same probe setup in A549 cells shows the preferential topology of targeted sites (both DHS[−]) toward the compacted core of a chromatin domain cluster. Scale bar 2 µm, inset 0.5 µm

Representative SIM sections of hybridized BJ1 (Fig. 5B–D) and A549 nuclei (Fig. 5E, F) reveal signals with a lateral diameter of ~120–150 nm. This length is at the diffraction limit of 3D-SIM, so their true diameters may be smaller. In BJ1 nuclei, both signal pairs of 6-kb probes 1 and 2 are consistently noted in low chromatin compaction zones at the periphery of chromatin domain clusters (CDCs). In A549 nuclei where both probes delineate DHS[−] sites, signals are closer toward the compact core of CDCs, yet excluded from the most compact zones in the interior of CDCs. Examples shown in Fig. 5C, D are suggestive of replicated DNA in close proximity (C) or separated by ~400 nm (D), since a second site of (replicated) signal pairs is seen elsewhere in these nuclei (Additional files 12, 13, 14: movies_3, 5; for a presumable G1 nucleus, see Additional file 15: movie_6). Replicated signal pairs hint to a consistent orientation since 3D distances between centroids of 6-kb probe 1 signals are found significantly smaller compared to 6-kb probe 2 centroid distances in both cell types (Additional file 16). Quantitation in BJ1 cells confirms enrichment of the 6-kb DHS[+] segment delineated by probe 1 in the low chromatin compaction class 2 and a depletion in classes 4–6 (Fig. 6a, b left panels), while the corresponding (DHS[−]) segment in A549 cells shows a significant enrichment in higher compaction classes (Fig. 6a, b right panels, Fig. 6c shows the direct comparison between BJ1 and A549 cells). Notably—except for few nuclei (see Additional file 17: for single-cell profiles)—probe 2 largely mirrors the 3D topography of probe 1 in both cell types. This was expected for A549 cells where both probes delineate DHS[−] segments. In BJ1 cells, the similar topography of both segments within the ANC may be explained as a passive consequence of the functionally important looping out of a DHS[+] segment with an active TRE enforcing the concomitant movement of an adjacent DHS[−] segment at this length scale. Minimal distance measurements between centroids of 6-kb probe 1 and 2 signals show a slightly higher, although not significant extension of the entire ~12-kb segment in BJ1 compared to A549 nuclei (Fig. 6d). This finding may suggest a difference in chromatin compaction of the respective target sites in the two cell lines. A detailed information on statistical values for all measurements is given in Additional file 18.

Fig. 6
figure6

Quantitative 3D mapping of two contiguous 6-kb DHS[+] and DHS[−] targets on DAPI intensity classes. a Relative signal distribution of probe 1 (green) and probe 2 (red) on respective DAPI intensity classes (all classes shown in gray) in BJ1 (left; N = 20) and A549 nuclei (right, N = 10). Note similar distribution patterns for both probes (n.s. at p < 0.5) within a cell type but distinct distribution between cell types. b Quantified levels of relative enrichment/overrepresentation (positive values) or depletion/underrepresentation (negative values) of probe signals relative to the classified DAPI signals show a highly significant overrepresentation in low-density class 2 and an underrepresentation in high-density classes 4–6. In A549 cells, both probes are significantly underrepresented in low-density classes 1 and 2 and overrepresented in classes 3–4. Error bars standard deviation of the mean. c Comparison between BJ1 (light green/light red) and A549 (dark green/dark red) confirms the distinct topography of both probes in distinct chromatin density compartments. Error bars standard deviation of the mean. d Minimal distances (nearest-neighbor analysis) between all differently labeled fluorescent signals for probe 1 and probe 2 in BJ1 and A549 cells, n.s. at p < 0.5 level. *p ≤ 0.05; **p ≤ 0.01; ***p ≤ 0.001

Discussion

In this study, we further explored the structural and functional organization of chromatin domain clusters (CDCs) at the level of super-resolution microscopy and compared for the first time the 3D nuclear topography of selected DHS[+] and DHS[−] sites which typically reflect active or inactive TREs, respectively. In two human types of cultured cells (BJ1 and A549), we found an enrichment of active TREs (DHS[+]) in the ANC, i.e., at the periphery of CDCs extending into the IC, whereas inactive TREs (DHS[−]) were enriched toward the more compacted interior of CDCs, constituting the ‘inactive’ nuclear compartment (INC). Further studies with other cell types and species are indicated to test whether our current results present a general feature of the ANC–INC model, which was based on studies of a variety of normal and cancer cell types from different mammalian species [43, 47, 48]. In case that the ANC–INC model stands further experimental scrutiny, the spatial separation of inactive TREs in the INC could be explained as additional protection against their inappropriate activation, but experimental evidence for this is lacking. The distinctly different topography of active and inactive TREs in the ANC and INC, respectively, suggests a dynamic organization of CDCs, which allows positional changes of TREs between the two compartments depending on their functional state.

Structural and functional organization of CDCs

Our study indicates a radially arranged structural organization of CDCs based on layers of different chromatin compaction with most decondensed, transcriptionally competent chromatin at the CDC periphery to the most compact chromatin within the CDC interior. This layered organization suggests transition zones of chromatin compaction between CDs with a fully ‘closed’ and a fully ‘open’ configuration [6, 54, 55]. Hi-C experiments revealed two higher-order compartments A and B, respectively, of ~1-Mb CDs corresponding to open (transcriptionally competent) and closed (transcriptionally silent) chromatin [33], but their relationship to CDCs has remained elusive. Current microscopic evidence demonstrates a significant enrichment of transcriptionally competent chromatin at the periphery and of repressed chromatin in interior of CDCs [43, 47, 48]. The term INC suggests a heterochromatic nature of compact CDs. According to the classical view of heterochromatin, facultative heterochromatic contains repressed genes, while constitutive heterochromatin is built up from repetitive blocks without interspersed coding genes. In both cases, the paucity of transcription was formerly considered as a hallmark of heterochromatin. Recent evidence, however, has shown that non-coding RNAs, called hetRNA, are transcribed from heterochromatin at a previously unexpected level, including pericentric and intergenic major satellite repeats [56].

Our initial assessment of the extent of minimal movements required for shifts of regulatory sequences between the decondensed periphery of CDCs and the compact interior shows that positional changes <100 nm may suffice. However, covered distances may be substantially larger since movements of genomic elements within the nuclear landscape typically occur along individual trajectories resembling an anomalous diffusion rather than the shortest possible path [19]. A considerable local dynamic of nucleosomes in living cells that drive chromatin accessibility was reported [57]. An in-depth exploration of the dynamic nature of CDCs, including movements of CDs and chromatin loops, requires live cell imaging with resolution <100 nm. The necessity to study large numbers of sites of specific types of TREs with special scrutiny of housekeeping, developmental and cell-type-specific genes adds another methodological challenge for future studies.

Perspectives of a functional interplay between CDCs and IC channels

3D/4D super-resolution microscopic studies are limited by the fact that only a few targets can be studied in each experiment. For our initial study, we selected eight targets on seven chromosomes (compare Fig. 3) and altogether mapped the 3D positions of about 300 DHS[+] and 400 DHS[−] sites within a total of 45 BJ1 and 30 A549 cell nuclei. Hi-C of cell populations with millions of fixed cells allows the genome-wide detection of billions (109) of pairwise 3D DNA–DNA contacts yielding a 1-kb resolution of topologically associating domains (TADs) [32, 33, 36]. In line with microscopic evidence, Hi-C led to the discovery of ~1-Mb-sized TADs carrying smaller sub-TADs and larger meta-TADs [32, 36, 37]. In contrast to microscopic studies, Hi-C provides the major advantage of direct comparisons with other genome-wide data sets, such as gene expression profiles, histone signatures or DHS sites (see, e.g., Roadmap Epigenomics Mapping Consortium data). Recent advancements of Hi-C have made it possible to explore the 3D organization of the entire genome in individual cell nuclei as well, although at the cost of a strongly reduced number of 3D contact sites detected in individual cell nuclei compared with Hi-C of cell populations [24, 58].

Considering these powerful possibilities, one may argue that Hi-C may replace microscopic studies, which is, however, not the case as outlined below. Hi-C experiments support a CT organization based on a 3D multiloop aggregate/rosette chromatin architecture [59]. At face value such an organization may leave sufficient inter-chromatin space to allow diffusion of macromolecules directly through open loops. Maeshima and colleagues have challenged this view with their proposal of a liquid drop model of chromatin domain organization [60, 61]. According to this model, chromatin domains are composed of irregularly folded, highly compacted 10-nm nucleosome fibers. The diffusion of non-coding RNAs and of single transcription factors into the interior of these compact CDs is possible, yet highly constrained [62]. Based on Monte Carlo simulations, the authors have proposed that small gene-specific transcription factors with a size of ~50 kD can penetrate into compact chromatin domains and search their target sequences, whereas large transcription complexes are excluded.

Taking into account a high level of compaction, the ANC–INC model [21] predicts that transcription factors and other functional proteins, which enter the nucleus via nuclear pores, reach their sites of action preferentially by constrained diffusion along routes provided by IC channels, which start at nuclear pores and pervade the nuclear interior between CDCs with finest branches extending into their interior. Similarly, the IC-channel network may serve for the rapid intranuclear distribution of mobile, non-coding RNAs involved in gene expression and for the export of mRNAs. Our observation of a shell-like organization of CDCs argues for the possibility that CDs with different chromatin compaction levels coexist in individual CDCs. Further studies are necessary to explore the diffusional constraints inflicted by the true physical compaction of individual CDs and TADs, respectively, as well as by structural entities, which may provide ‘stumbling blocks’ within open loops [21]. We expect that electron and super-resolved fluorescence microscopy will remain the methods of choice to analyze the true 3D geometry of CDCs, including the size, shape and extent of compaction or decondensation of individual CDs. Ongoing efforts to achieve multicolor visualization of specific DNA targets in live cells will help to address dynamic aspects of CDC organization [28, 63].

Conclusions

This study demonstrates the non-random distribution of active and inactive transcription regulatory elements (TREs) within the higher-order chromatin landscape of cell nuclei studied in a diploid human fibroblast line (BJ1) and an aneuploid human lung cancer cell line (A549). Data were obtained by 3D FISH with differently labeled DNA probe sets targeting sites with apparently active or inactive TREs and 3D quantitative image analyses of DAPI-stained nuclei recorded with 3D structured illumination microscopy (3D-SIM). Our results indicate a 3D organization of chromatin domain clusters (CDCs) with radially arranged layers of increasing chromatin compaction from the periphery toward the CDC core. Segments harboring active TREs are significantly enriched at the decondensed periphery of CDCs, while segments with inactive TREs are embedded within the more compacted interior layers. This difference suggests positional changes of TREs within CDCs depending on their functional state. Live cell studies with high resolution are required to directly observe a relocation of TREs within CDCs in line with their state of activity. A further improvement of resolution beyond the reach of SIM, achieved, for example, by single-molecule localization microscopy [64] is essential both with respect to the precision of target localization and with respect to the quantitative measurements of 3D chromatin domain compaction within CDCs, which depends critically on the resolution limit. Comparisons of individual cells in their living state and after 3D-FISH are further necessary to quantify at a given resolution the extent of potential changes of the 3D chromatin landscape due to fixation and DNA denaturation required for 3D-FISH [30, 65].

Materials and methods

Cells and culture conditions

BJ1-tert skin fibroblasts (ATCC # CRL-2522) and A549 cells Human Lung Carcinoma Epithelial Cells (ATCC #: CCL-185) were grown under conditions used for DNase-seq analysis. BJ1 cells: MEM, supplemented with 1.5 g/L sodium bicarbonate, 1 mM sodium pyruvate, 2 mM l-glutamine, 1× non-essential amino acids, 10% FBS, Pen-Strep (1×). A549 cells: F-12 K Medium, supplemented with 10% FBS, Pen-Strep (1×).

Generation of complex DNA probe sets

DHS profiles in BJ1 and A549 cells were identified on availability of NIH Roadmap Epigenomics Mapping Consortium data (www.roadmapepigenomics.org hg19). Genomic segments with or void of DHS[+] sites in BJ1 and A549 cells, respectively, were selected as DNA targets for the generation of two probe sets. Fosmid clones for fosmid pools 1 and 2 are described in Fig. 3. For optimal target representation, a pair of overlapping fosmid clones was selected for each region from the WIBR-2 human fosmid library (http://www.ncbi.nlm.nih.gov/clone/library/genomic/) and purchased from the BACPAC Resources Center. For clone ID and genomic positions, see Additional file 19. In total, 6-kb probes 1 and 2, located within fosmid clone G248P85778F6 (as detailed description in Fig. 5A), were generated by PCR of 1-kb subfragments and subsequent pooling of amplification products (for position and primer_seq, see Additional file 20). PCR was performed in 10 mM Tris–HCl with pH 8.3, 50 mM KCl, 2 mM MgCl2, 200 µM each dNTP, 1 µM each primer, 10 ng template DNA using 25 cycles of 94 °C for 30 s/56 °C for 30 s/72 °C for 1 min. Equal amounts of fosmid DNA- or PCR-amplified DNA assigned for labeling with either biotin or digoxigenin in the respective probe set were pooled and the pooled samples labeled by either biotin or digoxigenin by standard nick translation. Forty nanograms of labeled probe together with 20-fold excess of human COT-1 DNA was dissolved per 1 µl hybridization solution (2 × SSC/10% dextran sulfate/50% formamide). Hybridization efficiency and specificity of probes were verified on human metaphase chromosomes (Additional file 21).

Regulatory element annotation

The following data sets were converted from hg19 into hg38 using the liftover chain downloaded from the UCSC web-page and the liftover-tool implemented in the R package rtracklayer.

To identify regulatory elements, we accessed the Ensembl Regulation Database (v 88) [52] via the R package biomaRt. We used the data sets of Human Other Regulatory Regions, which only identified FANTOM annotations [49] as well as Human Regulatory Features to find annotated regulatory elements and Human Regulatory Evidence to get a better idea, what kind of chromatin marks or binding sites were found [66]. Type and sequence location of TREs targeted by targeted in this study are summarized in Additional file 2.

Pretreatment of cells for 3D-FISH and hybridization/detection setup

Unsynchronized cells grown up to ~60–70% confluency on high-precision borosilicate glass coverslips (170 ± 5 µm thickness) were fixed in 4% para-formaldehyde/PBS (10 min) followed by a stepwise replacement with PBS/0.05% Tween 20 and subsequent quenching of free aldehydes by 20 mM glycine (10 min). Cell and nuclear membranes were permeabilized by 0.5% Triton X-100/PBS (10 min), repeated freezing/thawing of cells in liquid N2 and subsequent gradual incubation in 0.1 N HCl (5 min). Cells were equilibrated in 2× SSC; at this step, RNA was removed by RNAse I treatment (100 μg/ml, 1 h at 37°). Cells were incubated in 50% formamide/2× SSC (pH = 7.0) at 4 °C until hybridization, at least overnight.

After simultaneous denaturation of cells and probes (76°/2 min), hybridization was performed at 37 °C for at least 48 h. Stringent washings in 0.1× SSC at 60 °C were followed by extensive blocking in 2% bovine serum albumin/0.5% Fish gelatin/4× SSCT for at least 2 h at RT. Probe detection was performed with avidin-Alexa488 (Molecular probes) and mouse-antidigoxigenin (Sigma) followed by an Alexa594-conjugated anti-mouse IgG (Molecular Probes). Cells were postfixed for 10 min in 4% formaldehyde/PBS, and DNA was counterstained with 1 µg/ml 4′,6-diamidino-2-phenylindole (DAPI) in 2× SSC. Samples were mounted in Vectashield antifade mounting medium (Vector Lab). A detailed protocol for 3D-FISH meeting the requirements for 3D-SIM is provided in [42].

3D-SIM

Super-resolution structured illumination imaging was performed on a DeltaVision OMX V3 system (Applied Precision Imaging/GE Healthcare) equipped with a 100×/1.4 UPlan S Apo oil immersion objective (Olympus), Cascade II:512 EMCCD cameras (Photometrics) and 405, 488 and 593 nm lasers (for detailed description, see [67]). Raw data image stacks were acquired with 15 raw images per plane (5 phases, 3 angles) and an axial distance of 125 nm and then computationally reconstructed with a Wiener filter setting of 0.002 and channel-specific optical transfer functions (OTFs) using SoftWoRx (Applied Precision). The reconstruction process generates 32-bit data sets with the pixel number doubled in the lateral axes, resulting in the pixel size being halved from 79–39.5 nm in order to meet the Nyquist sampling criterion. The level of spherical aberration was minimized and matched to the respective OTFs using immersion oil of different refractive indices (RI). Best results were typically obtained with OTFs measured on red, green (both 110 nm diameter) and blue (170 nm diameter) fluorescent FluoSpheres (Invitrogen) using RI 1.512, and sample acquisition with RI 1.512 for depth adjustment in the region of optimal reconstruction a few µm into the sample. Images from the different color channels were corrected for chromatic aberration in SoftWoRx with alignment parameters obtained from calibration measurements with 0.2-µm-diameter TetraSpeck beads (Invitrogen). To normalize all image stacks for subsequent image processing and data analysis, the original 32-bit images were shifted to positive values and transformed to 16-bit. All further image processing was performed in ImageJ (http://rsb.info.nih.gov/ij/). For a detailed description of methodological image quality assessment survey, see [68].

Chromatin compaction classification by 3D assessment of DAPI intensity classes

Nuclei voxels were identified automatically from the DAPI channel intensities using Gaussian filtering and automatic threshold determination. For chromatin quantification, a 3D mask was generated in ImageJ to define the nuclear space considered for the segmentation of DAPI signals into seven classes with equal intensity variance by an in-house algorithm described previously [43, 48], available on request. Briefly, a hidden Markov random field model classification was used, combining a finite Gaussian mixture model with a spatial model (Potts model), implemented in the statistics software R [69, 70]. This approach allows threshold-independent signal intensity classification at the voxel level, based on the intensity of an individual voxel. Color or gray value heat maps of the seven intensity classes in individual nuclei were performed in ImageJ. For a detailed description, see [46].

Semi-automatic segmentation of hybridization signals and their allocation on 3D chromatin compaction classes

Individual voxels of FISH signals of the respective marker channels were segmented using a semi-automatic thresholding algorithm (using custom-built scripts for the open-source statistical software R http://www.r-project.org, available on request). Xyz-coordinates of segmented voxels were mapped to the seven DNA intensity classes. The relative frequency of intensity weighted signals mapped on each DAPI intensity class was used to calculate the relative distribution of signals over chromatin classes. For 3D mapping of two contiguous differentially labeled 6-kb DNA probes (6-kb probes 1 and 2) relative to chromatin compaction classes, any fluorescent dot with a distance >0.5 µm from the nearest signal centroid of a differentially labeled target was attributed to background and eliminated from further consideration after signal segmentation with appropriate parameter settings. For each studied nucleus, the total number of voxels counted for each intensity class and the total number of voxels identified for the respective FISH signals were set to 1.

As an estimate of over/under representations (relative depletion/enrichment) of marker signals in the respective intensity classes, we calculated the difference between the percentage points obtained for the fraction of voxels for a given DAPI intensity class and the corresponding fraction of voxels calculated for the FISH signals. Calculations were performed on single-cell level and average values over all nuclei used for evaluation and plotting. For a detailed description, see [46].

Nearest-neighbor/minimal distance measurements

Nearest-neighbor/minimal distance measurements between centroid xyz coordinates of differently labeled segmented FISH ‘objects’ were taken using the TANGO Plugin for ImageJ/Fiji [43, 71, 72]. Mode-subtracted, 16-bit transformed 3D-SIM image stacks were imported into TANGO. Xyz centroid coordinates from segmented objects were extracted based on the geometrical gravity center of the segmented 3D foci and subsequently used for centroid mapping and nearest-neighbor analysis. Nearest-neighbor distances of different experiments were analyzed by pairwise t test comparison with Bonferroni correction of level of significance. For the measurements of minimal absolute distances between DAPI intensity classes, distances between voxels were calculated from their centroid. For each class, the distance from a voxel of this class to the nearest voxel of each other class was calculated.

Statistical evaluation

GraphPad Prism 6 was used for plots and statistical evaluations. Statistical differences were tested using the Wilcoxon rank sum test with continuity correlation as well as Student’s t test (two-tailed, p < 0.05). For interpolation models of DHS[+] and DHS[−] distributions, a second-order polynomial fit was used. SD and SEM were used for error bars, as indicated.

Data access

DHS profiling: www.roadmapepigenomics.org hg19

Fosmid library: http://www.ncbi.nlm.nih.gov/clone/library/genomic/

Abbreviations

3D FISH:

3D fluorescence in situ hybridization

3D SIM:

3D structured illumination microscopy

ANC/INC:

active/inactive nuclear compartment

CT:

chromosome territory

CD(C):

chromatin domain (cluster)

CTCF:

CCCTC-binding factor

DAPI:

4′,6-diamidino-2-phenylindole

DHS[+]:

DNAse I hypersensitive site indicates the presence of active regulatory sequences

DHS[−]:

DNAse I non-sensitive site in a given cell, although the site is DHS[+] in at least one other cell type or cell line, indicates the presence of inactive regulatory sequences (exception: probe 2 delineating a DHS[−] segment lacks the content of TREs)

DNase-seq:

DNAse I digestion and sequencing

Hi-C:

chromosome conformation capturing combined with deep sequencing

IC:

interchromatin compartment

OTF:

optical transfer function

PR:

perichromatin region

TAD:

topologically associating domain

TF:

transcription factor

TRE:

transcription regulatory element

References

  1. 1.

    Maston GA, Evans SK, Green MR. Transcriptional regulatory elements in the human genome. Annu Rev Genomics Hum Genet. 2006;7:29–59.

  2. 2.

    Sakabe NJ, Nobrega MA. Genome-wide maps of transcription regulatory elements. Wiley Interdiscip Rev Syst Biol Med. 2010;2:422–37.

  3. 3.

    Bell O, Tiwari VK, Thoma NH, Schubeler D. Determinants and dynamics of genome accessibility. Nat Rev Genet. 2011;12:554–64.

  4. 4.

    Boyle AP, Davis S, Shulha HP, Meltzer P, Margulies EH, Weng Z, Furey TS, Crawford GE. High-resolution mapping and characterization of open chromatin across the genome. Cell. 2008;132:311–22.

  5. 5.

    Cockerill PN. Structure and function of active chromatin and DNase I hypersensitive sites. FEBS J. 2011;278:2182–210.

  6. 6.

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

  7. 7.

    Pundhir S, Bagger FO, Lauridsen FB, Rapin N, Porse BT. Peak-valley-peak pattern of histone modifications delineates active regulatory elements and their directionality. Nucleic Acids Res. 2016;44:4037.

  8. 8.

    Crawford GE, Davis S, Scacheri PC, Renaud G, Halawi MJ, Erdos MR, Green R, Meltzer PS, Wolfsberg TG, Collins FS. DNase-chip: a high-resolution method to identify DNase I hypersensitive sites using tiled microarrays. Nat Methods. 2006;3:503–9.

  9. 9.

    Crawford GE, Holt IE, Whittle J, Webb BD, Tai D, Davis S, Margulies EH, Chen Y, Bernat JA, Ginsburg D, et al. Genome-wide mapping of DNase hypersensitive sites using massively parallel signature sequencing (MPSS). Genome Res. 2006;16:123–31.

  10. 10.

    Mercer TR, Edwards SL, Clark MB, Neph SJ, Wang H, Stergachis AB, John S, Sandstrom R, Li G, Sandhu KS, et al. DNase I-hypersensitive exons colocalize with promoters and distal regulatory elements. Nat Genet. 2013;45:852–9.

  11. 11.

    Sabo PJ, Humbert R, Hawrylycz M, Wallace JC, Dorschner MO, McArthur M, Stamatoyannopoulos JA. Genome-wide identification of DNaseI hypersensitive sites using active chromatin sequence libraries. Proc Natl Acad Sci USA. 2004;101:4537–42.

  12. 12.

    Thurman RE, Rynes E, Humbert R, Vierstra J, Maurano MT, Haugen E, Sheffield NC, Stergachis AB, Wang H, Vernot B, et al. The accessible chromatin landscape of the human genome. Nature. 2012;489:75–82.

  13. 13.

    Ay F, Bailey TL, Noble WS. Statistical confidence estimation for Hi-C data reveals regulatory chromatin contacts. Genome Res. 2014;24:999–1011.

  14. 14.

    Kieffer-Kwon KR, Tang Z, Mathe E, Qian J, Sung MH, Li G, Resch W, Baek S, Pruett N, Grontved L, et al. Interactome maps of mouse gene regulatory domains reveal basic principles of transcriptional regulation. Cell. 2013;155:1507–20.

  15. 15.

    Ma W, Ay F, Lee C, Gulsoy G, Deng X, Cook S, Hesson J, Cavanaugh C, Ware CB, Krumm A, et al. Fine-scale chromatin interaction maps reveal the cis-regulatory landscape of human lincRNA genes. Nat Methods. 2015;12:71–8.

  16. 16.

    Ramani V, Cusanovich DA, Hause RJ, Ma W, Qiu R, Deng X, Blau CA, Disteche CM, Noble WS, Shendure J, Duan Z. Mapping 3D genome architecture through in situ DNase Hi-C. Nat Protoc. 2016;11:2104–21.

  17. 17.

    van de Werken HJ, Landan G, Holwerda SJ, Hoichman M, Klous P, Chachik R, Splinter E, Valdes-Quezada C, Oz Y, Bouwman BA, et al. Robust 4C-seq data analysis to screen for regulatory DNA interactions. Nat Methods. 2012;9:969–72.

  18. 18.

    Jin W, Tang Q, Wan M, Cui K, Zhang Y, Ren G, Ni B, Sklar J, Przytycka TM, Childs R, et al. Genome-wide detection of DNase I hypersensitive sites in single cells and FFPE tissue samples. Nature. 2015;528:142–6.

  19. 19.

    Bronshtein I, Kanter I, Kepten E, Lindner M, Berezin S, Shav-Tal Y, Garini Y. Exploring chromatin organization mechanisms through its dynamic properties. Nucleus. 2016;7:27–33.

  20. 20.

    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:8044.

  21. 21.

    Cremer T, Cremer M, Hubner B, Strickfaden H, Smeets D, Popken J, Sterr M, Markaki Y, Rippe K, Cremer C. The 4D nucleome: evidence for a dynamic nuclear landscape based on co-aligned active and inactive nuclear compartments. FEBS Lett. 2015;589:2931–43.

  22. 22.

    Gaspar-Maia A, Alajem A, Meshorer E, Ramalho-Santos M. Open chromatin in pluripotency and reprogramming. Nat Rev Mol Cell Biol. 2011;12:36–47.

  23. 23.

    Sexton T, Cavalli G. The role of chromosome domains in shaping the functional genome. Cell. 2015;160:1049–59.

  24. 24.

    Stevens TJ, Lando D, Basu S, Atkinson LP, Cao Y, Lee SF, Leeb M, Wohlfahrt KJ, Boucher W, O’Shaughnessy-Kirwan A, et al. 3D structures of individual mammalian genomes studied by single-cell Hi-C. Nature. 2017;544:59–64

  25. 25.

    Cremer T, Cremer C. Chromosome territories, nuclear architecture and gene regulation in mammalian cells. Nat Rev Genet. 2001;2:292–301.

  26. 26.

    Cremer T, Cremer M. Chromosome territories. Cold Spring Harb Perspect Biol. 2010;2:a003889.

  27. 27.

    Lanctot 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:104–15.

  28. 28.

    Anton T, Leonhardt H, Markaki Y. Visualization of genomic loci in living cells with a fluorescent CRISPR/Cas9 system. Methods Mol Biol. 2016;1411:407–17.

  29. 29.

    Ma H, Naseri A, Reyes-Gutierrez P, Wolfe SA, Zhang S, Pederson T. Multicolor CRISPR labeling of chromosomal loci in human cells. Proc Natl Acad Sci USA. 2015;112:3002–7.

  30. 30.

    Markaki Y, Smeets D, Fiedler S, Schmid VJ, Schermelleh L, Cremer T, Cremer M. The potential of 3D-FISH and super-resolution structured illumination microscopy for studies of 3D nuclear architecture: 3D structured illumination microscopy of defined chromosomal structures visualized by 3D (immuno)-FISH opens new perspectives for studies of nuclear architecture. BioEssays. 2012;34:412–26.

  31. 31.

    Beliveau BJ, Boettiger AN, Avendano MS, Jungmann R, McCole RB, Joyce EF, Kim-Kiselak C, Bantignies F, Fonseka CY, Erceg J, et al. Single-molecule super-resolution imaging of chromosomes and in situ haplotype visualization using Oligopaint FISH probes. Nat Commun. 2015;6:7147.

  32. 32.

    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:376–80.

  33. 33.

    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:289–93.

  34. 34.

    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:381–5.

  35. 35.

    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:458–72.

  36. 36.

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

  37. 37.

    Fraser J, Ferrai C, Chiariello AM, Schueler M, Rito T, Laudanno G, Barbieri M, Moore BL, Kraemer DC, Aitken S, et al. Hierarchical folding and reorganization of chromosomes are linked to transcriptional changes in cellular differentiation. Mol Syst Biol. 2015;11:852.

  38. 38.

    Boulos RE, Drillon G, Argoul F, Arneodo A, Audit B. Structural organization of human replication timing domains. FEBS Lett. 2015;589:2944–57.

  39. 39.

    Pope BD, Ryba T, Dileep V, Yue F, Wu W, Denas O, Vera DL, Wang Y, Hansen RS, Canfield TK, et al. Topologically associating domains are stable units of replication-timing regulation. Nature. 2014;515:402–5.

  40. 40.

    Neems DS, Garza-Gongora AG, Smith ED, Kosak ST. Topologically associated domains enriched for lineage-specific genes reveal expression-dependent nuclear topologies during myogenesis. Proc Natl Acad Sci USA. 2016;113:E1691–700.

  41. 41.

    Yang R, Kerschner JL, Gosalia N, Neems D, Gorsic LK, Safi A, Crawford GE, Kosak ST, Leir SH, Harris A. Differential contribution of cis-regulatory elements to higher order chromatin structure and expression of the CFTR locus. Nucleic Acids Res. 2016;44:3082–94.

  42. 42.

    Markaki Y, Smeets D, Cremer M, Schermelleh L. Fluorescence in situ hybridization applications for super-resolution 3D structured illumination microscopy. Methods Mol Biol. 2013;950:43–64.

  43. 43.

    Smeets D, Markaki Y, Schmid VJ, Kraus F, Tattermusch A, Cerase A, Sterr M, Fiedler S, Demmerle J, Popken J, et al. Three-dimensional super-resolution microscopy of the inactive X chromosome territory reveals a collapse of its active nuclear compartment harboring distinct Xist RNA foci. Epigenetics Chromatin. 2014;7:8.

  44. 44.

    Schermelleh L, Heintzmann R, Leonhardt H. A guide to super-resolution fluorescence microscopy. J Cell Biol. 2010;190:165–75.

  45. 45.

    Demmerle J, Wegel E, Schermelleh L, Dobbie IM. Assessing resolution in super-resolution imaging. Methods. 2015;88:3–10.

  46. 46.

    Schmid VJ, Cremer M, Cremer T. Quantitative analyses of the 3D nuclear landscape recorded with super-resolved fluorescence microscopy. Methods. 2017;123:33–46

  47. 47.

    Hubner B, Lomiento M, Mammoli F, Illner D, Markaki Y, Ferrari S, Cremer M, Cremer T. Remodeling of nuclear landscapes during human myelopoietic cell differentiation maintains co-aligned active and inactive nuclear compartments. Epigenetics Chromatin. 2015;8:47.

  48. 48.

    Popken J, Brero A, Koehler D, Schmid VJ, Strauss A, Wuensch A, Guengoer T, Graf A, Krebs S, Blum H, et al. Reprogramming of fibroblast nuclei in cloned bovine embryos involves major structural remodeling with both striking similarities and differences to nuclear phenotypes of in vitro fertilized embryos. Nucleus. 2014;5:555–89.

  49. 49.

    Arner E, Daub CO, Vitting-Seerup K, Andersson R, Lilje B, Drablos F, Lennartsson A, Ronnerblad M, Hrydziuszko O, Vitezic M, et al. Transcribed enhancers lead waves of coordinated transcription in transitioning mammalian cells. Science. 2015;347:1010–4.

  50. 50.

    Lawrence M, Huber W, Pages H, Aboyoun P, Carlson M, Gentleman R, Morgan MT, Carey VJ. Software for computing and annotating genomic ranges. PLoS Comput Biol. 2013;9:e1003118.

  51. 51.

    Popken J, Dahlhoff M, Guengoer T, Wolf E, Zakhartchenko V. 3D structured illumination microscopy of mammalian embryos and spermatozoa. BMC Dev Biol. 2015;15:46.

  52. 52.

    Zerbino DR, Wilder SP, Johnson N, Juettemann T, Flicek PR. The ensembl regulatory build. Genome Biol. 2015;16:56.

  53. 53.

    Belmont AS, Bruce K. Visualization of G1 chromosomes: a folded, twisted, supercoiled chromonema model of interphase chromatid structure. J Cell Biol. 1994;127:287–302.

  54. 54.

    Boeger H, Shelansky R, Patel H, Brown CR. From structural variation of gene molecules to chromatin dynamics and transcriptional bursting. Genes (Basel). 2015;6:469–83.

  55. 55.

    Over RS, Michaels SD. Open and closed: the roles of linker histones in plants and animals. Mol Plant. 2014;7:481–91.

  56. 56.

    Saksouk N, Simboeck E, Dejardin J. Constitutive heterochromatin formation and transcription in mammals. Epigenetics Chromatin. 2015;8:3.

  57. 57.

    Hihara S, Pack CG, Kaizu K, Tani T, Hanafusa T, Nozaki T, Takemoto S, Yoshimi T, Yokota H, Imamoto N, et al. Local nucleosome dynamics facilitate chromatin accessibility in living mammalian cells. Cell Rep. 2012;2:1645–56.

  58. 58.

    Nagano T, Lubling Y, Stevens TJ, Schoenfelder S, Yaffe E, Dean W, Laue ED, Tanay A, Fraser P. Single-cell Hi-C reveals cell-to-cell variability in chromosome structure. Nature. 2013;502:59–64.

  59. 59.

    Knoch TA, Wachsmuth M, Kepper N, Lesnussa M, Abuseiris A, Ali Imam AM, Kolovos P, Zuin J, Kockx CE, Brouwer RW, et al. The detailed 3D multi-loop aggregate/rosette chromatin architecture and functional dynamic organization of the human and mouse genomes. Epigenetics Chromatin. 2016;9:58.

  60. 60.

    Maeshima K, Ide S, Hibino K, Sasai M. Liquid-like behavior of chromatin. Curr Opin Genet Dev. 2016;37:36–45.

  61. 61.

    Maeshima K, Hihara S, Eltsov M. Chromatin structure: does the 30-nm fibre exist in vivo? Curr Opin Cell Biol. 2010;22:291–7.

  62. 62.

    Maeshima K, Kaizu K, Tamura S, Nozaki T, Kokubo T, Takahashi K. The physical size of transcription factors is key to transcriptional regulation in chromatin domains. J Phys: Condens Matter. 2015;27:064116.

  63. 63.

    Shao S, Zhang W, Hu H, Xue B, Qin J, Sun C, Sun Y, Wei W. Long-term dual-color tracking of genomic loci by modified sgRNAs of the CRISPR/Cas9 system. Nucleic Acids Res. 2016;44:e86.

  64. 64.

    Kirmes I, Szczurek A, Prakash K, Charapitsa I, Heiser C, Musheev M, Schock F, Fornalczyk K, Ma D, Birk U, et al. A transient ischemic environment induces reversible compaction of chromatin. Genome Biol. 2015;16:246.

  65. 65.

    Hubner B, Cremer T, Neumann J. Correlative microscopy of individual cells: sequential application of microscopic systems with increasing resolution to study the nuclear landscape. Methods Mol Biol. 2013;1042:299–336.

  66. 66.

    Andersson R, Gebhard C, Miguel-Escalada I, Hoof I, Bornholdt J, Boyd M, Chen Y, Zhao X, Schmidl C, Suzuki T, et al. An atlas of active enhancers across human cell types and tissues. Nature. 2014;507:455–61.

  67. 67.

    Dobbie IM, King E, Parton RM, Carlton PM, Sedat JW, Swedlow JR, Davis I. OMX: a new platform for multimodal, multichannel wide-field imaging. Cold Spring Harb Protoc. 2011;2011:899–909.

  68. 68.

    Ball G, Demmerle J, Kaufmann R, Davis I, Dobbie IM, Schermelleh L. SIMcheck: a toolbox for successful super-resolution structured illumination microscopy. Sci Rep. 2015;5:15915.

  69. 69.

    Pau G, Fuchs F, Sklyar O, Boutros M, Huber W. EBImage: an R package for image processing with applications to cellular phenotypes. Bioinformatics. 2010;26:979–81.

  70. 70.

    R Core Team.  A language and environment for statistical computing.  R Foundation for Statistical Computing, Vienna, Austria; 2013. http://www.R-project.org/.

  71. 71.

    Ollion J, Cochennec J, Loll F, Escude C, Boudier T. TANGO: a generic tool for high-throughput 3D image analysis for studying nuclear organization. Bioinformatics. 2013;29:1840–1.

  72. 72.

    Cerase A, Smeets D, Tang YA, Gdula M, Kraus F, Spivakov M, Moindrot B, Leleu M, Tattermusch A, Demmerle J, et al. Spatial separation of Xist RNA and polycomb proteins revealed by superresolution microscopy. Proc Natl Acad Sci USA. 2014;111:2235–40.

Download references

Authors’ contributions

TC, JS and MC initiated and conceived the study; SJ provided sequence data for PCR-generated probes and fosmid probes; MC planned and performed all 3D-FISH experiments shown in this study; YM initiated and performed initial experiments; VJS developed and performed quantitative 3D image analyses; FK performed quantitative 3D image analyses; IH analyzed regulatory elements in relevant databases; HL provided input for the 3D imaging part; AM carried out 3D SIM microscopy and image processing; MC and TC wrote the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We thank Anna Jauch (University of Heidelberg, Germany) for kindly performing M-FISH analysis on A549 cells. We are grateful to Shreeram Akilesh (Dept of Pathology, University of Washington Seattle, USA) and to Tobias Ragoczy (Altius Institute for Biomedical Sciences, Seattle, USA) for critical reading the manuscript. The Center of Advanced Light Microcopy (CALM) of the LMU Biocenter (headed by H. Leonhardt/H. Hartz) was essential for the 3D-SIM recording.

Competing interests

The authors declare to have no competing interests.

Availability of supporting data

All data generated or analyzed during this study are included in this published article and its Additional files.

Consent for publication

All authors read and approved the manuscript.

Ethical approval and consent to participate

Not applicable.

Funding

This work was supported by funding of the Ludwig Maximilians-Universität (LMU Munich) to TC.

Publisher’s Note

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

Author information

Correspondence to Marion Cremer or Thomas Cremer.

Additional files

13072_2017_146_MOESM1_ESM.jpg

Additional file 1. ANC–INC network model of nuclear organization based on spatially co-aligned active and inactive nuclear compartments (for detailed information, see [21]). Nuclear organization according to co-aligned 3D networks of an active (ANC) and an inactive nuclear compartment (INC). The ANC is a composite structural and functional entity of a 3D-channel network, the ‘Inter-chromatin-compartment’ (IC) together with the decondensed periphery of a higher-order chromatin network, which is built up from ~1-Mb chromatin domains (CDs), representing basic units of larger chromatin domain clusters (CDCs). The decondensed periphery of CDCs is known as the perichromatin region (PR). According to this model, the PR harbors regulatory and coding sequences of active genes and represents the preferential nuclear subcompartment for transcription, RNA-splicing, and possibly also for DNA replication and repair. Small chromatin loops expand from the perichromatin region into the interior of IC channels which start/end at nuclear pore complexes. Nuclear bodies are located within the IC, which serves as a transport system for macromolecule complexes. The INC is represented by the compacted core of CDCs enriched in markers for silent chromatin (Fig. modified from [21]).

13072_2017_146_MOESM2_ESM.xlsx

Additional file 2. Overview of type, sequence location and open chromatin marks of TREs targeted by probe sets used in this study (fosmid pools 1 and 2; 6-kb probes 1 and 2). Data are based on [52]. For the segment covered by 6-kb probe 2, no TREs were identified in the used databases. Sequence coordinates represent the coordinates for hg38 after conversion of hg19 into hg38 by liftover. The used databases do not provide data on ‘open chromatin marks’ H3K9ac and H3K4me3 as additional indirect information for their state of activity at the respective loci in BJ1 cells. In the sheet ‘open chromatin marks,’ available data for IMR90 fibroblasts are shown instead. IMR90 cells show an almost identical DHS[+] profile to BJ1 cells (www.roadmapepigenomics.org); accordingly, similar epigenetic signatures between both cell lines can be assumed.

13072_2017_146_MOESM3_ESM.jpg

Additional file 3. Metaphase spreads of BJ1 cells: Five out of six randomly selected Giemsa stained metaphase spreads reveal an inconspicuous diploid chromosome set of n = 46, XY (n = 45 in metaphase 2 is likely due to loss of one chromosome during preparation).

13072_2017_146_MOESM4_ESM.mov

Additional file 4. Movie_1 entire image stack of a BJ1 nucleus hybridized with fosmid pools 1 and 2. DAPI-stained DNA after intensity classification shown in gray. Fosmid pool 1 (green), fosmid pool 2 (red).

13072_2017_146_MOESM5_ESM.jpg

Additional file 5. Single-cell profiles for target sites of fosmid pools 1 and 2 mapped to chromatin compaction classes in BJ1 cells for illustration of intercellular variability. (A) Mapping profiles from ten randomly chosen individual nuclei illustrate consistent distinct distribution profiles of fosmid pool 1 (green) toward low chromatin compaction classes and of pool 2 (red) toward higher chromatin compaction classes. (B) Standard deviations of relative probe signal distributions of all evaluated nuclei (compare Fig. 4 for standard errors of the mean (SEM). (C) Standard deviations of DAPI signal distribution on classes (compare Fig. 4 for SEM).

13072_2017_146_MOESM6_ESM.jpg

Additional file 6. DHS profiles and fosmid clones used for target regions delineating DHS[−] sites on different chromosomes in A549 nuclei. (A) Selected regions with clones of fosmid pool 1 (green) and (B) with clones of fosmid pool 2 (red). DHS profile in black (browser shots adopted from http://encodeproject.org/). Note: probe sets are identical to probe sets shown in Fig. 3.

13072_2017_146_MOESM7_ESM.pdf

Additional file 7. M-FISH karyotype analysis of A549 cells. (A) Representative karyotype obtained by M-FISH after combinatorial labeling of chromosome-specific paint probes with seven fluorochromes. (B) Quantitation of 20 metaphases reveals a karyotype with 62–66 chromosomes with consistent structural rearrangements involving chromosomes 1,2,3,6,8,11,15,19,20. This constellation allows for fosmid pool 1 up to 19, for fosmid pool 2 up to five distinct hybridization sites in a nucleus (compare Additional file 6).

13072_2017_146_MOESM8_ESM.mov

Additional file 8. Entire image stack of A549 nucleus after hybridization with fosmid pools 1 and 2. DAPI-stained DNA after intensity classification shown in gray. Fosmid pool 1 (green), fosmid pool 2 (red).

13072_2017_146_MOESM9_ESM.jpg

Additional file 9. 3D nuclear topography and quantitative mapping of ~40 kb targets of DHS[−] regions in A549 nuclei. (A) Part of a SIM light-optical section from a whole nucleus acquisition with representative inset magnifications. DAPI-stained DNA after intensity classification shown as gray gradations and color heat map, respectively. Segmented signals delineating targets both of fosmid pool 1 (green) and fosmid pool 2 (red) show a similar location with regard to chromatin compaction classes (asterisks in color heat maps, pool 1 (green), pool 2 (black). Scale bar 2 µm, insets 0.5 µm. (B) Quantified distributions (N = 10 nuclei) of fosmid pools 1 (green) and 2 (red) within respective chromatin compaction classes (all classes shown in gray). (C) Quantified levels of relative enrichment (positive values) or depletion (negative values) of fosmid pool 1 and pool 2 signals show an enrichment of signals in higher compaction classes. Error bars = standard deviation of the mean *p ≤ 0.05, **p ≤ 0.01, ***p ≤ 0.001.

13072_2017_146_MOESM10_ESM.jpg

Additional file 10. Single-cell profiles for target sites of fosmid pools 1 and 2 mapped to chromatin compaction classes in A549 cells for illustration of intercellular variability. (A) Mapping profiles of A549 nuclei (N = 10) illustrate for most nuclei fairly similar distribution profiles of fosmid pool 1 (green) and fosmid pool 2 (red). (B) Standard deviations of relative probe signal distributions of all evaluated nuclei (compare Additional file 8 for standard errors of the mean (SEM). (C) Standard deviations of DAPI signal distribution on classes (compare Additional file 8 for SEM).

13072_2017_146_MOESM11_ESM.pdf

Additional file 11. Significance values for relative signal distributions of fosmid pools 1 and 2 in BJ1 and A549 cells.

13072_2017_146_MOESM12_ESM.mov

Additional file 12. Entire image stack of the nucleus shown in Fig. 5B. DAPI-stained DNA after intensity classification (blue). Segmented probe signals (green/red, seen in section z19 and 23) reveal two sites of double signal pairs, suggestive of ongoing/post-replication. Note that the respective section shown in Fig. 5B is horizontally flipped for arrangement.

13072_2017_146_MOESM13_ESM.mov

Additional file 13. Entire image stack of the nucleus shown in Fig. 5C. DAPI-stained DNA after intensity classification (blue). Segmented probe signals (green/red, seen in sections z11 and 15) reveal two sites of double signal pairs, suggestive of ongoing/post-replication.

13072_2017_146_MOESM14_ESM.mov

Additional file 14. Entire image stack of the nucleus shown in Fig. 5D. DAPI-stained DNA after intensity classification (blue). Segmented probe signals (green/red seen in sections z27 and 34) reveal two sites of double signal pairs, suggestive of ongoing/post-replication.

13072_2017_146_MOESM15_ESM.mov

Additional file 15. Entire image stack of a presumable G1 nucleus. DAPI-stained DNA after intensity classification (blue). Segmented probe signals (green/red seen in sections z24 and 35) reveal two sites of a single signal pair, suggestive of G1.

13072_2017_146_MOESM16_ESM.jpg

Additional file 16. 3D distances between centroids of 6-kb probe 1 (green) and centroids of 6-kb probe 2 in BJ1 (left) and in A549 cells (right). Distance measurements are restricted to distances <500 nm presumably comprising only sister chromatids of S/G2 nuclei. The smaller distances between centroids of green signals (probe 1; DHS[+] in BJ1 cells, DHS[−] in A549 cells) compared to distances between centroids of red signals (probe 2; DHS[−] both in BJ1 and A549 cells) hint to a consistent orientation of these segments irrespective of DNAse I sensitivity.

13072_2017_146_MOESM17_ESM.jpg

Additional file 17. Single-cell profiles of 6-kb probes 1 and 2 targets mapped to chromatin compaction classes in BJ1 (AC) and A549 cells (DF) for illustration of intercellular variability. (A) Mapping profiles from ten randomly chosen individual BJ1 nuclei for illustration of intercellular variabilities and similarities of relative signal distribution of probe 1 (green) and probe 2 (red) within DAPI intensity classes. Note an only marginal signal representation in classes 5–7. (B) Standard deviations of relative probe signal distributions of all evaluated nuclei (compare Fig. 6 for standard errors of the mean (SEM). (C) Standard deviations of DAPI signal distribution on classes (compare Fig. 6 for standard errors of the mean (SEM). (D) Respective mapping profiles from 10 individual A549 nuclei of relative signal distribution of probe 1 (green) and probe 2 (red) within DAPI intensity classes. Profiles show an overall broader distribution range compared to BJ1 nuclei. (E) Standard deviations of relative probe signal distributions of all evaluated nuclei (compare Fig. 6 for standard errors of the mean (SEM). (F) Standard deviations of DAPI signal distribution on classes (compare Fig. 6 for standard errors of the mean (SEM).

13072_2017_146_MOESM18_ESM.pdf

Additional file 18. Significance values for relative signal distributions of 6-kb probes 1 and 2 in BJ1 and A549 nuclei.

13072_2017_146_MOESM19_ESM.pdf

Additional file 19. fosmid ID (G248 library #) and sequence alignment of fosmids used. Data are based on hg19.

13072_2017_146_MOESM20_ESM.xls

Additional file 20. Position and primer sequences used for amplification of ~1-kb subfragments for assembling of 6-kb probes 1 and 2. Data are based on hg19.

13072_2017_146_MOESM21_ESM.jpg

Additional file 21. FISH of fosmid pairs on normal human metaphases for verification of specificity. (A) Human metaphase ideogram with marked positions of tested fosmids. (BE) Metaphase spreads after FISH with (B) fosmid pairs G248P8092D1/G248P89035F6 mapped on 1p and G248P80020B1/G248P8977D10 mapped on 12q, (C) G248P83004C6/G248P82547F4 mapped on 2p and G248P87313E8/G248P85778F6 mapped on 2q, (D) G248P8631F6/G248P88483C3 mapped on 3p and G248P87150D8/G248P89650D7 mapped on 5q, (E) G248P83624H8/G248P83627E4 mapped on 3p and G248P80223H2/G248P84663H7 mapped on 13q. All tested probes show a specific hybridization signal at the expected chromosomal position.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Cremer, M., Schmid, V.J., Kraus, F. et al. Initial high-resolution microscopic mapping of active and inactive regulatory sequences proves non-random 3D arrangements in chromatin domain clusters. Epigenetics & Chromatin 10, 39 (2017) doi:10.1186/s13072-017-0146-0

Download citation

Keywords

  • Transcription regulatory sequences
  • DNAse I hypersensitive sites
  • Super-resolution microscopy
  • Chromatin domain
  • Nuclear architecture
  • Active and inactive nuclear compartment
  • Chromatin compaction