Spatial genome organization: contrasting views from chromosome conformation capture and fluorescence in situ hybridization

1Fluorescence microscopy has been instrumental to our understanding of the spatial organization of the genome. Detection of chromosomes by fluorescence in situ hybridization (FISH) in either fixed cells or living cells reveals specific chromosome territories (CTs) that are nonrandomly organized in vertebrate nuclei. Gene-rich chromosomes generally cluster at the nuclear center, while gene-poor chromosomes and chromatin domains often locate near the nuclear periphery. This level of genome organization is stochastic: Loci at the nuclear periphery in a mother cell can localize to either the nuclear or nucleolar periphery in daughter nuclei.

2Cytologically, gene-dense chromatin containing active loci is observed to decondense and dynamically loop out from CT cores both at individual loci and chromosome-wide. A propensity for clustering of gene-rich chromatin domains from the same or different chromosomes is also seen.

3FISH also identifies differences in chromatin condensation at the submegabase level, between genomic regions in the same cell, for a given region during differentiation in vitro and in vivo, or between wild-type and mutant cells. FISH has also been used to examine tissue-specific colocalization of long-range enhancers and their target genes

4Although visually compelling, FISH and live-cell imaging are restricted to viewpoints corresponding to the regions detected by the probes used, are low-throughput assays, and have limited spatial resolution, although superresolution microscopy is improving the latter. In contrast, the chromosome conformation capture (3C) technique and its derivatives—including circular 3C (4C), 3C carbon copy (5C), and chromosome capture followed by high-throughput sequencing (Hi-C) (for review, see; )—offer a genome-wide perspective on genome organization. By digesting and ligating chromatin within formaldehyde cross-linked cells, these methods measure the frequency at which sequences are joined together by ligation. Unlike imaging in which distances are measured directly in individual cells, 3C methods infer physical proximity by considering that ligation frequency is more or less inversely proportional to the original spatial separation of the two sequences.

5The limited resolution of early Hi-C studies (∼1 Mb) gave coarse-grained views of nuclear organization that agree well with microscopy data, reproducing the existence of CTs, the spatial clustering of gene-rich chromosomes in the nuclear center, looping out of gene-rich domains from CTs, and the tendency for active genomic domains to associate. Whereas FISH is inherently a single-cell/allele technique, most Hi-C analyses give data sets of ligation frequencies averaged across all cells in a population. A single-cell Hi-C study recently overcame this limitation, albeit at the expense of resolution, and confirmed the stochastic nature of CT architecture and the general disposition of active chromatin to reside toward the outside of CTs.

6More recently, higher-resolution Hi-C and 5C studies have suggested that mammalian genomes are organized into globules or domains (topologically associating domains [TADs]) within which cross-linked associations are enriched. The largely invariant nature of TAD boundaries between different cell types suggests that they reflect a basic property of genome organization. Indeed, the average size of these domains (∼900 kb) is similar to the ∼1-Mb domains suggested from microscopy to form the basic building blocks (∼500-nm diameter) of CT architecture. FISH showed that sequences within a TAD intermingle to a much greater extent than those that span adjacent TADs.

7In addition to providing insight into global chromatin organization, 3C methods have been used to identify associations between specific sequences involved in gene regulation—including silent Hox gene clusters regulated by the polycomb-repressive complexes (PRCs), domains of clustered active genes (Hox and α-globin clusters), and “looping” between long-range enhancers and their target genes—and between sites bound by proteins involved in chromosome architecture.

8By supposing that high-frequency contacts are sequences close to each other in vivo, data from 3C experiments have been used to reconstruct the likely trajectory of chromatin fibers and relate chromatin folding to gene regulation. While the resulting models are informative, some are not consistent with the views assayed by FISH. This points to differences in what 3C methods and FISH can detect and demands a careful comparison of molecular and cytological data from the same experimental system. In most studies that have attempted to do this, interprobe FISH distances up 500 nm (sometimes even more) have been taken as evidence validating 3C interactions. However, such comparatively large physical separations do not seem compatible with interactions at the molecular scale. Here we investigate chromatin organization of the murine HoxD locus with 5C and three-dimensional (3D) FISH in different developmental and activity states and after perturbing epigenetic regulation. We identified some conditions under which the standard interpretation of 5C data as a measure of spatial proximity produces views of chromatin topography that are not compatible with those obtained by FISH. We conclude that results obtained by either high-resolution 3C-type experiments or FISH alone must be interpreted with caution when studying chromatin architecture.

9In terms of chromatin organization, the Hox gene clusters are among the most extensively studied regions in mammalian cells. By FISH, they have a visibly compact conformation in mouse ESCs (mESCs), where they are maintained in an inactive state by the PRCs. This is also seen in mouse embryonic development. These observations are compatible with the condensed appearance on polytene chromosomes of the bithorax complex, one of the two Hox clusters in Drosophila, mitigating concerns about artifacts due to the FISH procedure. A relatively compact chromatin domain with multiple interactions within it is also predicted from 3C analysis of silent HOX loci in human embryonal carcinoma (EC) cells and 4C analysis of silent Hox loci in mESCs and mouse embryos.

10In contrast, the lack of long-range interactions captured by 5C at silent parts of HOXA in human fibroblasts was suggested to reflect an extended linear chromatin domain. Barring any source of technical bias, this result either reflects tissue-specific organization of silent HOX genes or might point to differences between what is sometimes captured with 3C methods and by FISH. We first investigated whether FISH and 5C invariably yield similar views of chromatin organization by examining a ∼670-kb region containing the inactive HoxD cluster (Fig. 1A) in undifferentiated mESCs. A prominent feature is the high level of 5C interactions detected within the HoxD–Evx2 region (Fig. 1B; Supplemental Figs. S1A, S6A). The pattern of these 5C contacts (Fig. 1C; Supplemental Fig. S1B) is similar to the one that we reported for the silent HOXD cluster in human EC cells, supporting a conserved chromatin conformation at HoxD in pluripotent stem cells. If interaction frequency is viewed as inversely proportional to spatial distance, our 5C analysis indicates that the region encompassing the Hoxd and Evx2 genes is folded into a structure that is the most compact part of the region under study (Fig. 1D). This is consistent with our previous observation that the silent HoxD cluster in mESCs adopts a visibly compact structure by FISH.

11The extent of the strong 5C long-range contacts across HoxD in mESCs corresponds well to the domain of histone H3 Lys27 trimethylation (H3K27me3) deposited by the PRC2 complex and to the sub-TAD identified at Evx2–HoxD by Hi-C in mESCs (Fig. 2A; ). At our sequencing depth, we did not see evidence of the two TADs that extend 3′ and 5′ of HoxD in the published Hi-C analysis and saw little evidence of other significant long-range 5C interactions in the regions flanking HoxD.

12Retinoic acid (RA)-directed differentiation of mESCs results in activation of Hoxd genes and visible chromatin decompaction of the locus reminiscent of the puffed appearance of the active bithorax locus in Drosophila polytene chromosomes. However, this contrasts with views of activated Hox loci extrapolated from 5C analysis of human fibroblasts and 4C in mouse embryos. From those data, it was suggested that active Hox genes form compact domains.

13To determine whether these contrasting models stem from inherent differences in what FISH and 3C techniques can detect, we examined chromatin organization across HoxD with 5C and FISH in mESCs differentiated for 3 d (D3) with RA (Fig. 2A; Supplemental Figs. S1A, S6A). Under these conditions, gene activation and loss of H3K27me3 mainly occur in the 3′ part of HoxD (Fig. 2A;; ). The high 5C interactions across HoxD that we detected in undifferentiated mESCs are strongly reduced in D3 cells but, like H3K27me3, are not completely lost (Fig. 2A,B; Supplemental Fig. S1C). This is consistent with the visible unfolding of HoxD chromatin seen upon mESC differentiation (Fig. 2C,D) and the loss of long-range 3C and 5C contacts reported upon HOXA activation.

14FISH using probe pairs across the whole region interrogated by 5C (Supplemental Table S1) confirmed a visible chromatin decompaction during differentiation and that this is restricted to the region overlapping HoxD and Evx2 (Hoxd13–Hoxd3; Prox–Hoxd10) (Fig. 2D; Supplemental Tables S2, S3). However, although FISH indicates that the HoxD cluster itself is by far the most unfolded part of this entire 700-kb region in D3 cells, the capture frequencies of many 5C fragment pairs within HoxD are still comparable with, or greater than, those in the flanking regions. For example, the global control region (GCR)–Lnp region 5′ of HoxD is more visibly compact than HoxD in D3 cells (Fig. 2D) but has lower 5C reads (Fig. 2A, arrow; Supplemental Fig. S1A).

15Apart from the loss of interactions within HoxD, the other major change in the 5C profile is the increased interactions captured during differentiation between 3′ Hoxd genes and the region extending through Mtx2 and into the gene desert beyond (Fig. 2B; Supplemental Fig. S1C). However, the distribution of FISH interprobe distances across these regions (Hoxd3–Mtx2P and Mtx2P–3′Mtx2) remains largely unchanged upon differentiation (Fig. 2D; Supplemental Fig. S2; Supplemental Tables S2–4), suggesting no major change in chromatin compaction.

16These data caution against the simple assumption that spatial proximity between genomic regions can always be inferred at high resolution from either contact frequencies in 3C-based assays or FISH alone.

17To gain more insight into views of long-range chromatin conformation from 3C and FISH, we examined HoxD in another developmental context. Digit growth and patterning depend on the correct regulation of expression from 5′ Hoxd genes, initiated in the posterior part of the distal limb bud during the late phase of limb development. This is controlled by multiple long-range regulatory elements/enhancers spread throughout a large gene desert 5′ of HoxD (Fig. 3A;; ). In 4C assays from embryonic day 12.5 (E12.5) distal limb buds, these enhancers could be captured cross-linked to both the promoter of Hoxd13 and each other. It was therefore suggested that these elements are brought close to each other and Hoxd13 by multiple chromatin loops.

18Consistent with this model, FISH in the E11 distal forelimb reveals a posterior-specific high (30%) colocalization (≤200 nm) frequency of Hoxd13 alleles, with the best-characterized regulatory element 180 kb upstream of Hoxd13—the GCR (Fig. 3B,C; Supplemental Table S5; ). This is recapitulated in mesenchymal cell lines derived from the anterior (A2) or posterior (P2) of dissected E10.5 forelimb buds (Fig. 3B,C), which also preserve the anterior–posterior differences in 5′ Hoxd expression (Supplemental Fig. S3; ). The Hoxd13–GCR colocalization frequency (35%) in P2 cells, which strongly express Hoxd13, is significantly higher (P = 0.005) than that in A2 cells (Fig. 3C; Supplemental Table S5). However, as we reported previously in limb tissue, there is little colocalization between Hoxd13 and another of the regulatory elements, island III (IslIII), or between GCR and IslIII in A2 and P2 cells, and there is no elevated colocalization frequency in P2 cells compared with A2 (Fig. 3C; Supplemental Table S5). This argues against a model in which the different enhancers are brought together by discrete contacts into a regulatory hub. Rather, the low median Hoxd13–IslIII and GCR–IslIII interprobe FISH distances (221–250 nm) (Fig. 3D; Supplemental Table S6) suggest that the whole region from IslIII to Hoxd13 is configured into a generally compact structure in both cell lines and tissue from the distal limb.

19A key attribute of 4C is that ligation products are detected looking out from single viewpoints, and so the conformation of whole chromosomal domains cannot be extrapolated from these restricted data. We therefore used 5C to compare all of the cross-linked associations captured between sequences in the 1.2-Mb region containing HoxD and its 5′ regulatory domain (Fig. 3A) in A2 and P2 cells. Two interaction domains are clearly identified in both anterior and posterior limb cells (Fig. 4A; Supplemental Figs. S4A, S6B), demarcated by a boundary at Hoxd9-10. This corresponds well to the major TAD boundary reported in mESC Hi-C data (Fig. 4B, bottom heat map; ). The 5′ Hoxd region, particularly around Hoxd13, has a high cross-linking efficiency with Evx2 and Lnp and throughout the 5′ gene desert. In contrast, the 3′ end of HoxD predominantly interacts with the 3′ gene desert.

20This aspect of HoxD compartmentalization is consistent with 4C results from E9.5 and E12.5 limbs. To examine this more closely, we extracted the interactome of IslIII, GCR, Hoxd13, Hoxd9–11, and Hoxd1 from our 5C data sets and generated “virtual 4C” profiles by considering each of the fixed viewpoints as pseudobaits (Fig. 4C; Supplemental Fig. S4C). This is reminiscent of the approach used to compare interaction profiles derived from 4C and Hi-C in mESCs. Although virtual 4C confirms the generally higher Hoxd9–10 and Hoxd1 interaction frequencies with the region 3′ of HoxD compared with the 5′ flanking region, for the fixed viewpoint of Hoxd13 interaction, frequencies up to 400 kb into the 5′ regulatory domain are not substantially higher than those extending 400 kb in the other direction. This argues against strictly exclusive conformational compartments for 5′ Hoxd and more 3′ regions in the limb and is not consistent with the 4C signals enriched throughout the region 5′ of Hoxd13, compared with the 3′ direction, during limb development.

215C heat maps (Fig. 4AB; Supplemental Figs. S4AB) and virtual 4C (Fig. 4C; Supplemental Fig. S4C) also do not identify specific interactions between the regulatory elements of the 5′ gene desert and Hoxd13 that would be indicative of discrete enhancer–promoter loops. Instead, the data point to chromatin conformations with enriched interaction frequencies between the enhancers (e.g., IslI, IslIII, and GCR) and the whole 5′ regulatory domain, including Hoxd13 (Fig. 4A,B, larger area with dotted outline). This is consistent with the view obtained by FISH of enhancers and target genes brought into proximity by a general compaction of the whole region rather than by specific loops (Fig. 3D; ). For the region up to and just beyond the GCR, these interactions are generally more frequent in Hoxd13-expressing P2 cells than in A2 cells (Fig. 4B,C, smaller dotted area; Supplemental Fig. S4B,C). FISH also captures this difference (Fig. 3). In contrast, 5C interaction frequencies between IslIII and the rest of the regulatory domain extending toward Hoxd13 are no different between the two cell lines (Fig. 4C), and this is also reflected in the FISH data.

22Limb cells do not lend themselves readily to further dissection of the mechanisms that control chromatin structure. In contrast, in mESCs, the visibly compact chromatin structures of the silent HoxB and HoxD loci have been shown to be dependent on the PRC complexes PRC1 and PRC2.

23H3K27me3 is catalyzed by the PRC2 component Ezh2, whose chromatin immunoprecipitation (ChIP) profile across HoxD in mESCs also parallels that of H3K27me3 (Fig. 5A; ). The distribution of H3K27me3 and Ezh2 also corresponds to the extent of the strong 5C interaction domain, suggesting that PRC2 might be important for chromatin conformation at HoxD (Fig. 5A; Supplemental Figs. S5A, S7). To directly demonstrate a role of PRC2 in chromatin compaction, we carried out both 5C and 3D-FISH on ESCs mutant for PRC2 (Eed−/−). Ezh2 is degraded, and H3K27me3 is lost globally from these cells. Ring1B (PRC1) binding is also lost from Hox loci in PRC2-null cells, likely due to the dependence on H3K27me3 recognition by the Cbx subunits of canonical PRC1. The 5C interaction domain over HoxD was virtually eliminated in Eed−/− cells (Fig. 5B; Supplemental Fig. S5B), coincident with a visible decompaction that is detected by FISH and restricted to HoxD (Fig. 5D,E, Hoxd13–Hoxd3; Supplemental Table S7, wt/EED−/−, P < 0.0001; ). Therefore, both FISH and 5C appear to show that PRC2 is required to maintain Hox loci in a compact silent chromatin conformation in mESCs.

24Binding of the Ring1B PRC1 subunit across HoxD parallels that of both H3K27me3 and Ezh2 in wild-type mESCs (Fig. 5A; ). PRC2 is necessary but not sufficient for chromatin compaction at HoxD. In Ring1B−/− mESCs, HoxD visibly decompacts relative to the wild type (P < 0.0001) to the same extent as that seen in the absence of PRC2; the interquartile range of interprobe distances in Ring1B−/− cells (300–569 nm) is similar to that in Eed−/− cells (300–539 nm), and both are different from wild-type mESCs (192–300 nm) (Fig. 5D,E; Supplemental Table S7). Hox genes are also activated despite the persistence of H3K27me3 and PRC2 at these loci (Fig. 5C;; ). Thus, we expected that 5C data from Ring1B−/− mESCs would resemble those from PRC2 mutant cells. Instead, high 5C contacts frequencies remain across HoxD in Ring1B−/− mESCs in a pattern similar to that seen in wild-type cells (Fig. 5C, arrow; Supplemental Figs. S5C, S7). Thus, if interaction frequency is considered as inversely proportional to spatial distance, the 5C data would suggest that HoxD remains at least partially folded in the absence of Ring1B, which is in conflict with the decompact HoxD domain seen by FISH (Figs. 5D,E, 6). Moreover, in Ring1B−/− cells, 5C data appear to indicate that the HoxD locus itself is the most highly folded part of the entire region under investigation (Fig. 5C), whereas, to the contrary, FISH suggests that HoxD is the most unfolded part of the whole region in these mutant cells (Fig. 5E).

25We previously showed that Ring1B reintroduced into Ring1B−/− cells recapitulates the binding profile of the endogenous protein across HoxD and restores chromatin compaction as measured by FISH. Strong and reproducible 5C contacts are also restored to a pattern comparable with wild-type cells at HoxD in the cells rescued with wild-type Ring1B (Fig. 5F, arrow; Supplemental Fig. S5D,S7). Chromatin folding and compaction at HoxD does not require the enzymatic activity of Ring1B, since complementation of Ring1B−/− cells with the ubiquitination-deficient Ring1B I53A mutant also restores both the levels and the pattern of 5C contacts (Fig. 5G, arrow). These results are consistent with our previous study demonstrating that Ring1B visibly compacts chromatin structure at Hox loci in ESCs independent of its histone ubiquitination activity.

26With rather few exceptions, there have been limited systematic and detailed comparisons between 3C-type data sets and cytological assays of 3D chromatin organization. Here we directly compared 5C and FISH data across a region of the mouse genome encompassing the HoxD locus in different cell types and activity states and in the absence of major epigenetic regulators of the locus. We identified some situations in which the 5C and FISH data are in agreement, as, for example, in undifferentiated mESCs (Figs. 1, 2; Supplemental Fig. S1) and in the absence of the PRC2 complex (Fig. 5; Supplemental Figs. S5, S7). However, in other cases, such as upon gene activation during ESC differentiation (Fig. 2; Supplemental Fig. S1) or mutation of the PRC1 complex (Fig. 5; Supplemental Fig. S5), domains of elevated 5C cross-linking frequencies conflict with the visible decompaction of these same regions observed by FISH (Fig. 6). This suggests that, at high resolution, 3C interaction frequency might not always simply reflect physical distances. This is important to consider if 3C-type information is ever to be used in modeling genome organization at the level of sub-TADs and beyond.

27In undifferentiated mESCs, 5C data identify a domain of robust interactions across HoxD that is consistent with its visibly compact structure relative to the surrounding gene deserts (Fig. 1B; Supplemental Fig. S1A; ). Upon RA-induced differentiation and activation of 3′ Hoxd genes, 5C interactions within HoxD are reduced compared with undifferentiated mESCs (Fig. 2; Supplemental Figs. S1, S6). This is consistent with FISH analysis showing chromatin unfolding strictly at HoxD and not at the flanking regions (Fig. 2). However, the 3′ Hoxd region simultaneously gains 5C contacts with the upstream gene desert region, while there are no significant reductions in distances detectable by FISH, which might have been expected if enhanced cross-linking simply arose from spatial proximity (Fig. 2D).

28What does visibly change upon mESC differentiation is the nuclear localization of both the 3′ Hoxd genes and the upstream gene desert as they gain the ability to loop out of the territory occupied by the bulk of chromosome 2. Looping out from CTs is dynamic, and, genome-wide, the ability of sequences to become cross-linked and captured to distant sequences—even those located on other chromosomes—by Hi-C corresponds remarkably well to the looping-out frequency. We suggest that access to this nuclear compartment outside of the CT core rather than actual physical distance per se may be responsible for the high 5C contact frequencies between 3′ Hoxd genes and the upstream genomic region in differentiated mESCs. The fact that 3C ligation products originate mainly from insoluble aggregates within unlysed swollen nuclei rather than on soluble chromatin raises the possibility that nuclear structures enhance or promote 3C associations, as was suggested for the active β-globin gene and its enhancer. This may also explain why higher 3C contact frequencies and inferred structurally compact domains were previously found associated with activated Hox loci in human fibroblasts and mouse embryos.

29The compact structure of Hox loci in mESCs depends on the PRCs, and the HoxD 5C interaction domain in mESCs (Fig. 1) almost exactly corresponds to that of Ezh2 (PRC2) binding and H3K27me3 enrichment (Fig. 5A; ). We demonstrated the role of PRC2 in maintaining this compact domain by showing that, in PRC2 Eed−/− mutant mESCs, 5C domain interactions are completely lost and that chromatin unfolding visualized by FISH is restricted to the HoxD cluster itself (Fig. 5B; Supplemental Fig. S5B, S7). To our knowledge, this is the first instance in which the genetic deletion of a major epigenetic regulator has been shown to have such a profound effect of almost eliminating a 3C interaction domain, and, in this instance, both 5C and FISH results are concordant (Fig. 6). However, this is not the case for the data from PRC1 Ring1B−/− mutant mESCs, where FISH indicates an unfolding of HoxD to an extent similar to that seen in PRC2-null cells, yet a 5C interaction domain—albeit weakened compared with the wild type—remains over HoxD (Fig. 5C; Supplemental Fig. S5C). These data suggest that 5C can capture substantial cross-linked ligation products from within a visibly unfolded chromatin domain.

30Chromatin composition might provide one explanation for this. While HoxD remains coated with PRC2 proteins in Ring1B−/− cells, it is completely devoid of them in Eed−/− cells. Formaldehyde preferentially forms cross-links with lysine, tryptophan, and cysteine side chains in proteins. Just considering the core components of PRC2 (Eed, EzH2, and Suz12), each molecule of the complex provides >200 extra of these reactive residues, including >150 lysine residues alone. While the stoichiometry of PRC2 per nucleosome is not known, given the blanket ChIP signal for PRC2 across HoxD, it is plausible that the thousands of extra formaldehyde cross-linkable amino acid side chains present in PRC1 mutant cells due to the persistent binding of PRC2 are responsible for the elevated 5C signals in a structure that is much more open than the data would suggest. Similarly, the residual levels of H3K27me3 at the 5′ end of HoxD could also explain why, despite FISH indicating that HoxD is the most decondensed part of the region under study in differentiated ESCs, 5C signals across HoxD generally exceed those outside of this region (Fig. 2). Another confounding factor could be the low frequency of ligation products produced in 3C reactions, even for enhancer–promoter interactions, relative to the total number of fragments present. If this represents the bona fide frequency of interaction between two genomic loci, then it may be below the level at which FISH is able to distinguish this from background. While we raised these issues from our study of HoxD using 5C and FISH, there may be different factors that could influence the interpretation of FISH with other 3C data (4C or Hi-C) and the outcome of these studies at other genomic loci with different compaction states and chromatin flavors.

31We conclude that, at high resolution, products captured by 3C do not always simply represent spatial proximity or molecular interaction between two DNA sequences but may arise from indirect cross-linking of chromatin fragments that may be hundreds of nanometers apart in nuclear space in the original cells, perhaps exacerbated by the partial chromatin decondensation that accompanies the initial digestion and SDS steps of 3C prior to ligation. This can be compounded by differences in the nuclear environment of the regions and/or the protein composition of the chromatin fiber. This is unlikely to be an issue when only considering low-resolution views of spatial nuclear organization. Of course, we cannot exclude the possibility that the FISH technique is also contributing to the apparent conflict in the data at high resolution. This could be due to an inability of FISH to capture weak or transient interactions or to perturbations to chromatin ultrastructure as a result of the heat denaturation step. Nonetheless, we think that FISH quite faithfully reflects the chromatin compaction state because, in all cases in which a direct comparison has been made, the data inferred from FISH are completely consistent with analysis of the same regions in vivo in living cells and with the biophysical properties of chromatin from those same regions assayed in vitro. Moreover, a careful quantitative analysis by superresolution microscopy before and after 3D-FISH shows evidence that FISH preserves many aspects of chromatin structure and organization.

32A decade ago, it was realized that the formation of cross-linked nuclear aggregates or “nuclear crumbs” presented problems in extrapolating data from formaldehyde cross-linking to nuclear organization, and our data support that view. We suggest that visual and molecular approaches are complementary to each other and that models of 3D genome organization should be extrapolated from data validated by independent methods.

33Wild-type (WTF18), Ring1B−/−, and Eed mutant (17Rn5-3354SB; Eed−/− B1.3 and Eed−/− G8.1) mESCs were grown on mitomycin C-treated primary embryonic fibroblasts (PEFs) derived from E12.5 mouse embryos as described. OS25 ESCs were cultured and differentiated for 3 d (2 d with RA) as described. For FISH and 5C analysis, ESCs were trypsinized, and the PEFs were removed by allowing them to reattach to the tissue culture plastic twice for 30 min in LIF-containing medium.

34Mesenchymal cell lines derived from Immortomouse E10.5 anterior (A2) and posterior (P2) distal forelimb buds were cultured as described previously.

35Paraformaldehyde (pFA)-fixed cells were permeabilized in 0.5% Triton X-100, washed in PBS, and stored at −80°C. 3D-FISH was carried out as previously described. Fosmid clones used as FISH probes are listed in Supplemental Table S1.

36Slides were imaged and analyzed as described previously. The statistical significance of differences in interprobe distances was assessed using the nonparametric Mann-Whitney U-test. Each data set consisted of between ∼80 and 120 loci for each cell line and each probe combination.

37Cells were fixed with 1% formaldehyde for 10 min at room temperature. Cross-linking was stopped with 125 mM glycine for 5 min at room temperature followed by 15 min on ice. Cells were centrifuged at 400g for 10 min at 4°C, supernatants were removed, and cell pellets were flash-frozen on dry ice.

38Cell pellets were treated as previously described. Briefly, 10 million to 20 million fixed cells were incubated for 15 min on ice in 200 μL of lysis buffer (10 mM Tris at pH 8.0, 10 mM NaCl, 0.2% NP40, supplemented with fresh protease inhibitor cocktail). Cells were then disrupted on ice with a dounce homogenizer (pestle B; 2 × 20 strokes); cell suspensions were transferred to Eppendorf tubes and centrifuged at 2000g for 5 min. Supernatants were removed, the cell pellets were washed twice with 100 μL of 1× EcoRI buffer (New England Biolabs), and the cell pellet was resuspended in 100 μL of 1× EcoRI buffer and divided into two Eppendorf tubes. We added 1× EcoRI buffer (337 μL) to each tube, and the mixture was incubated for 10 min at 65°C with 0.1% SDS. Forty-four microliters of 10% Triton X-100 was added before overnight digestion with 400 U of EcoRI. The restriction enzyme was then inactivated by adding 86 μL of 10% SDS and incubation for 30 min at 65°C. Samples were then individually diluted into 7.62 mL of ligation mix (750 μL of 10% Triton X-100, 750 μL of 10× ligation buffer, 80 μL of 10 mg/mL of BSA, 80 μL of 100 mM ATP, 3000 cohesive end units of T4 DNA ligase) and incubated at for 2 h 16°C.

393C libraries were incubated overnight at 65°C with 50 μL of Proteinase K (10 mg/mL) and an additional 50 μL of Proteinase K the following day for 2 h. The DNA was purified by one phenol and one phenol–chloroform extraction and precipitated with 0.1 vol (800 μL) of 3 M NaOAc (pH 5.2) and 2.5 vol of cold EtOH (20 mL). After at least 1 h at −80°C, the DNA was centrifuged at 20,000g for 25 min at 4°C, and the pellets were washed with cold 70% EtOH. DNA was resuspended in 400 μL of TE (pH 8.0) and transferred to Eppendorf tubes for another phenol–chloroform extraction and precipitation with 40 μL of 3 M NaOAc (pH 5.2) and 1.1 mL of cold EtOH. DNA was recovered by centrifugation and washed eight times with cold 70% EtOH. Pellets were then dissolved in 100 μL of TE (pH 8.0) and incubated with1 μL of 10 mg/mL RNase A for 15 min at 37°C.

405C primers covering the USP22 (mm9, chr11: 60,917,307–61,017,307) and HoxD (mm9, chr2: 73,750,000–74,910,000) regions were designed using my5C.primer and the following parameters: optimal primer length of 30 nucleotides (nt), optimal TM of 65°C, and default primer quality parameters (mer: 800, U-blast: 3, S-blasr: 50). Primers were not designed for large (>20-kb) and small (<100-bp) restriction fragments, low-complexity and repetitive sequences, or when there were sequence matches to more than one genomic target. The USP22 regions was used to assess the success of each 5C experiment but was not used for further data normalization or quantification.

41The universal A-key (CCATCTCATCCCTGCGTGTCTCCGACTCAG-[5C-specific]) and the P1-key tails ([5C-specific]-ATCACCGACTGCCCATAGAGAGG) were added to the forward and reverse 5C primers, respectively. Reverse 5C primers were phosphorylated at their 5′ ends. An alternating design consisting of 133 primers in the HoxD region (66 forward and 67 reverse primers) was used for analysis of ESCs. An extended alternating design of 203 primers (99 forward and 104 reverse primers) was used for the limb cells. Primer sequences are listed in Supplemental Table S8.

425C libraries were prepared and amplified with the A-key and P1-key primers as described previously. Briefly, 3C libraries were first titrated by PCR for quality control (single band, absence of primer dimers, etc.) and to verify that contacts were amplified at frequencies similar to that usually obtained from comparable libraries (same DNA amount from the same species and karyotype). In general, we used 1–11 μg of 3C library per 5C ligation reaction.

435C primer stocks (20 μM) were diluted individually in water on ice and mixed to a final concentration of 0.002 μM. Mixed diluted primers (1.7 μL) were combined with 1 μL of annealing buffer (10× NEBuffer 4 [New England Biolabs]) on ice in reaction tubes. Salmon testis DNA (1.5 μg) was added to each tube, followed by the 3C libraries and water to a final volume of 10 μL. Samples were denatured for 5 min at 95°C and annealed for 16 h at 48°C. Ligation with 10 U of Taq DNA ligase was performed for 1 h at 48°C. One-tenth (3 μL) of each ligation was then PCR-amplified individually with primers against the A-key and P1-key primer tails. We used 26 or 28 cycles based on dilution series showing linear PCR amplification within that cycle range. The products from two to four PCR reactions were pooled before purifying the DNA on MinElute columns (Qiagen).

445C libraries were quantified on agarose gels and diluted to 0.0534 ng/μL (for Xpress template kit version 2.0) or 0.0216 ng/μL (for ion PGM template OT2 200 kit). One microliter of diluted 5C library was used for sequencing with an ion PGM sequencer. Samples were sequenced onto ion 316 chips following either the ion Xpress template kit version 2.0 and ion sequencing kit version 2.0 protocols or the ion PGM template OT2 200 kit and ion PGM sequencing 200 kit version 2.0 protocols as recommended by the manufacturer (Life Technologies).

45Analysis of the 5C sequencing data was performed as described earlier. The sequencing data were processed through a Torrent 5C data transformation pipeline on Galaxy (https://main.g2.bx.psu.edu). Posterior-enriched and D3-enriched 5C interactions were obtained by subtracting anterior and undifferentiated 5C sequencing data, respectively. Data were normalized by dividing the number of reads of each 5C contact by the total number of reads from the corresponding sequence run. All scales correspond to this ratio multiplied by 103. The numbers of total reads and used reads are provided for each experiment in Supplemental Table S9. The unprocessed heat maps of the normalized 5C data sets can be found in Supplemental Figures S6 and S7. 5C data sets (accession no. GSE61814) can be downloaded from the Gene Expression Omnibus (GEO) (http://www.ncbi.nlm.nih.gov/geo).