Loss of compartments and TADs in Prophase

Hi-C analysis showed that G₂-arrested cells display all features characteristic of vertebrate interphase cells (8). First, chromosomes form territories as indicated by relatively high levels of intra-chromosomal interactions (3). Second, chromosomes display the characteristic pattern of active A- and inactive B- compartments as revealed by the plaid pattern of Hi-C interactions (3) (Fig. 1B). The locations of these compartments in G₂ resembled those detected in exponentially growing cells, though the compartment signal strength was stronger and the pattern sharper in the synchronous cells likely as a result of uniformity in cell cycle stage. Third, TADs were readily visible in the Hi-C interaction maps as squares of relatively high interaction frequencies along the diagonal (Fig. 1C). TAD boundaries were similar in position and strength to those in non-synchronous cells as determined using an insulation score calculated from a 250 kb sliding window (32) (Fig. S5). Finally, we analyzed how contact frequency (P) of contacts between locus pairs depends on their genomic distance (s). P(s) decays with genomic distance and this relationship changes with cell cycle stages (8). For G₂ cells, we found P(s) to be highly similar to that observed previously in non-synchronous cells (Fig. S6, Fig. S7). Together, these analyses show that G₂ chromosomes, which are composed of two closely aligned and likely catenated sister chromatids, are organized similarly to G₁ chromosomes (8).

This interphase chromosomal organization was rapidly lost upon release of cells into prophase. As soon as 5 minutes after release we detected a marked reduction in the typical plaid pattern of long-range interactions, indicating a loss of compartments (Fig. 1B). By 10 minutes (late prophase) compartments were mostly gone. At the same time, TADs were also lost (Fig. 1C, Fig. S8).

We used eigenvector decomposition to quantify the disappearance of compartments (33). The first eigenvector readily captured compartments at t = 0 and 2.5 minutes, but starting at t = 5 minutes it explained progressively less of the variance in the Hi-C interaction maps, indicating weakening of the compartment structure. By t = 7.5 minutes the strength of the first eigenvector fell to 17% (from 80% at t = 0 minutes) and by t = 10 minutes, it no longer captured compartments. Loss of compartments was also quantified by calculating the ratio of A-to-A or B-to-B interactions over A-to-B interactions over the time course. From t = 0 to t = 2.5 minutes and onward this fraction decreased steadily, indicating that preferential interactions within compartments are lost (Fig. 1D, Fig. S9).

The strength of TADs can be quantified using the insulation score, which indicates the amount of contacts formed across a locus up to a certain distance (32). TAD boundaries have a low score (indicative of high insulation), whereas loci inside TADs show a high score value (no insulation). The genome-wide variance of insulation scores provides a quantitative measure of the presence of TADs (8). Starting at t = 2.5 minutes the variance of the insulation profiles progressively decreased, indicating loss of TADs (Fig. 1C, Fig. S5B). By t = 7.5 minutes the variance was reduced more than 2-fold and by t = 10 minutes no TADs were detected. This conclusion was confirmed by plotting the average Hi-C interaction pattern at and around TAD boundaries identified in G₂ at different time points during mitosis (Fig. 1C). Insulation is strongest in G₂ and by late prophase insulation values are near background levels (quantified in Fig. S10). We conclude that compartments and TADs disappear rapidly during early prophase.

By late prophase, when sister arms have resolved (11, 34), and around the time of nuclear envelope breakdown (t ~ 7.5–10 minutes), the Hi-C maps are characterized by a general decay of contact frequency P with genomic distance s (Fig. 2A). The shape of the P(s) curve changes as prophase progresses. In G₂ cells, it is shallow (P(s) ~ s^−0.5) up to a distance of several hundred kb, reflecting compaction within TADs (35, 36), but for larger distances the decay becomes steeper. During prophase, the initial shallow decay extends for longer-range interactions, with a steeper drop at 2 Mb at t = 10 minutes, suggesting a higher degree of compaction. As we demonstrate below, this decay and shape are consistent with formation of a linearly arranged, layered organization of the chromosome (8), where the size of each layer corresponds to the position of the steep drop in the P(s) curve.

Open in a separate window

Fig. 2

Prophase chromosomes fold as axially compressed loop arrays

(A) Genome-wide curves of contact frequency P(s) vs genomic distance s, normalized to unity at s=100 kb. The curves are derived from prophase Hi-C data at the indicated time points after release from G₂ arrest. The dotted line indicates P(s) = s^−0.5 observed for mitotic chromosomes (8). (B) Overview of the coarse grained model of prophase chromosomes. The chromosome is compacted into a series of consecutive loops and compressed into a cylindrical shape. The loop bases form a scaffold at the chromosomal axis, each loop occupies a cylindrical sector of height h and angular size , oriented at angle ϴ_i. The coarse-grained model predicts the P(s) curve to have three distinct regions: an intra loop (I), intra layer (II) and inter layer (III) regions. (C) The best fitting P(s) predictions by the coarse grained model for late prophase (t = 7.5 minutes) under two different assumptions on loop orientations: (top panels) uncorrelated and (bottom panels) correlated orientations of consecutive loops. Uncorrelated angular loop orientations lead to a plateau in P(s) in the intra-layer, whereas correlated angles lead to the experimentally observed P(s) = s^−0.5 (right panels). (D) Polymer models of prophase chromosomes. Chromatin fibers are modeled as chains of particles (dark grey circles), compacted into arrays of consecutive loops (loop bases indicated in orange). Chromosomes are compacted into a cylinder with a density of one nucleosome per 11×11×11nm cube (lower right). (E) Goodness of fit for simulated vs experimental P(s). Polymer simulations were performed for a range of loop densities and loop lengths, and for each simulation P(s) was calculated. The heatmap shows the quality of a match between the predicted and experimental P(s) curves at late prophase (t = 7.5 minutes). (F) P(s) derived from late prophase Hi-C experiments (green line) and the best fitting polymer models (grayscale lines). Average loop size and linear density of loops along the chromosome axis are listed.. (G) Top and side view of the best fitting polymer model of late prophase chromosomes. Loops bases are shown in red and several loops rendered in different colors. (H) The average loop size and linear density of the 3 best-fitting models of prophase chromosomes at different time points.

Appearance of a second diagonal band in Hi-C maps from prometaphase cells

At t = 15 minutes, when cells have entered prometaphase, the Hi-C maps produce a P(s) curve with a drop at 2 Mb. Strikingly, a distinct second diagonal band appears, running in parallel with the primary diagonal for all loci and chromosomes (Fig. 1B, Fig. S11). This second diagonal represents increased interaction frequencies between any pair of loci separated by around 3 Mb. At 15 min, this feature is clearly observed in P(s) plots as a local peak at ~3 Mb (Fig. 3A, Fig. S6, Fig. S7). As cells progress through prometaphase, the position of the drop in P(s) and the position of the second diagonal migrate to larger genomic distances (Fig. 3A; Fig. S6, S7). By t = 60 minutes, when compact metaphase chromosomes have formed, the second diagonal is positioned at ~12 Mb and appears more diffuse. The second diagonal appears in all chromosomal maps and its position is independent of chromosome size over two orders of magnitude (Fig. S11). The appearance and movement of the second diagonal is not dependent or affected by nocodazole: no nocodazole was added to the t = 15 minute sample, and a replicate Hi-C dataset obtained from a culture collected at t = 30 minutes in the absence of nocodazole was nearly identical to the data obtained in the presence of nocodazole (Table S2, Fig. S7, Fig. S8). Together, these P(s) curves reveal a periodicity of interactions that reflects chromosome structure at the scale of megabases.

Open in a separate window

Fig. 3

Helical organization of prometaphase chromosomes

(A) Genome-wide curves of contact frequency P(s) vs genomic distance (separation, s), normalized to unity at s=100 kb. The curves are derived from Hi-C data obtained from prometaphase cells (t = 10–60 minutes after release from G₂ arrest). The dashed line indicates P(s) = s^−0.5. Arrows indicate positions of a local peak in P(s) representing the second diagonal band observed in Hi-C interaction maps. (B) The coarse grained model of prometaphase chromosomes with staircase loop arrangement. Left, top: the staircase loop arrangement implies that loops rotate in genomic order around a central scaffold (see Supplemental Materials). Left, bottom: angles of adjacent loops are correlated and steadily increasing, reflecting helical arrangement of loops. Right: this helical arrangement can be observed as gyres by DNA staining and a helical scaffold can be observed in cells expressing GFP-tagged condensins. (C) The best fitting P(s) predictions by the staircase coarse grained model for late prometaphase t = 30 minutes (30min; left panel) and t = 60 minutes (60min, right panel) after release from G₂ arrest (Hi-C data: colored lines; model; gray lines). (D) Polymer model of prometaphase chromosomes. Chromosomes are modeled as arrays of consecutive nested loops with a helical scaffold (outer loops in red, inner loops in blue, also indicated diagrammatically bottom right). (E) Goodness of fit for simulated vs experimental P(s). Polymer simulations were performed varying the helix height (nm), the size of a helical turn (Mb), and the sizes of inner and outer loops. For each simulation P(s) was calculated. The heatmaps show the quality of the best match between the predicted and experimental P(s) at prometaphase (t = 30 minutes), when two out of four parameters were fixed to the specified values. (F) P(s) derived from prometaphase Hi-C experiments (colored lines) and the best fitting polymer models (gray lines). Left panel: t = 30 minutes, right panel t = 60 minutes after release from G₂ arrest. Average size of outer and inner loops, the length of a helix turn and the helical pitch are indicated. (G) Parameters of the helical scaffolds from the best fitting polymer models. X-axis: ratio of the radius of the helical scaffold to that of the whole chromatid; Y-axis: ratio of the pitch to the helix radius. The dashed lines show the corresponding values (0.46 and 2.5122) for the optimal space-filling helix (84). Classical solenoid configurations are predicted to be in sector III, while the “spiraling staircase” configurations are in I and II. On the right, three examples of models of type I, II and III are shown with loops bases in red and several individual loops rendered in different colors. Also shown is a schematic of a prometaphase chromosome with the helical winding of loops indicated by arrow around the loop array. (H) Parameters of the best 3 models of prometaphase chromosomes at different time points.

The only known regular periodic structural feature of chromosomes is helical coiling, which was first described in 1880 (37) and can be observed in certain chromosome preparations (10, 38–40). Experimentally induced banding of the chromosome arms is much more irregular (41). In helical chromosomes with a pitch (the length of a complete turn) of ~3 Mb, each locus is in relatively close proximity to loci located one turn up or down the chromosome, i.e. 3 Mb up or downstream along the DNA. The progressive movement of the second diagonal band to larger distances during prometaphase would reflect an increased “winding up” and shortening of the helix.

Although this is the first Hi-C data revealing the helical coiling of chromosomes, the chicken DT40 late prometaphase (t = 60 minutes) Hi-C contact maps strongly resemble those for mitotic human HeLa S3 chromosomes that we had reported earlier (8). In fact, re-analysis of mitotic HeLa S3 Hi-C data in more detail by deeper sequencing revealed a weak second diagonal band at ~10 Mb distance (Fig. S12), suggesting this periodic folding is a conserved feature of vertebrate mitotic chromosomes.

Testing models of chromosomes

Previous studies suggested that mitotic chromosomes are organized as arrays of consecutive loops emanating from a condensin-rich scaffold, forming a polymer bottlebrush (42, 43), with a layered organization of loops (6, 19, 44). To understand chromosome organization at different stages of compaction we built coarse-grained models of chromosomes as arrays of loops aiming to reproduce P(s) curves of Hi-C data, separately for each time point. In these models, a chromosome is represented by a cylinder with an axial scaffold; loop bases are arranged consecutively along the scaffold, and each chromatin loop, emanating from the scaffold in a particular direction, is represented by a blob of loci (Fig. 2B, see Supplemental Material). Loops are regularly placed along the axis, with angular positions determined by a stochastic model; loop sizes are exponentially distributed and bases of loops are not positioned at defined genomic sequences or loci (8). Analysis of condensin ChIP data for DT40 cells (45) supports sequence-independent positioning for >95% of loops (Supplemental Material, Fig. S13; see discussion). For specific models of loop arrangements, presented below, the P(s) curve can be found analytically as the return probability of a stochastic process describing angular positions of loops (see Supplemental Material, section “Coarse-grained model of contact probability decay in mitotic chromosomes”). The resulting P(s) always has three regions (Fig. 2B): (i) the intra-loop region at short separations, where two loci are likely to be within the same loop and P(s) reflects the internal organization of loops; (ii) the “intra-layer” region at larger genomic separations, where loci are located on different loops within the same axial layer of the cylindrical chromosome and P(s) reflects the specific arrangement of loops relative to each other; (iii) the “inter-layer” that appears as a steep drop in contact frequency at large genomic distances, where loci are located in loops that are so distant along the scaffold that their blobs can no longer overlap. In the P(s) plot of experimental Hi-C data throughout mitosis, the intra-layer region and the drop-off can be readily discerned (Fig. 2A, 2C).

Prophase chromosomes

The coarse-grained models show that the relative orientation of consecutive loops strongly affects the shape of the P(s) curve in the intra-layer region. If the orientations of consecutive loops are independent of each other, the contact frequency P(s) does not decay with genomic distance in the intra-layer region, as any pair of loops within a layer are equally likely to interact (Fig. 2C). In contrast, introducing correlations between orientations of consecutive loops, i.e. forcing neighboring loops to project in similar directions, makes them follow an angular random walk. The angular random walk is a 1D random walk on a circle and has a return probability of P(s)~s^−0.5 until the full turn is made by the walk. The P(s)~s^−0.5 decay followed by a drop is in good agreement with the late prophase Hi-C data (t = 7–10 minutes – Fig. 2C). Taken together, these results suggest that by late prophase chromosomes are already organized into arrays of consecutive loops with correlated angular orientations.

We developed detailed polymer models to test whether specific classes of conformations can reproduce experimental Hi-C data, though they do not prove mechanisms by which these structures form. Further, all our simulations produce equilibrium models and do not reflect kinetics of chromosome folding. In these models, chromatin is represented as a 10 nm fiber (46, 47), where one monomer corresponds to one nucleosome (Fig. 2D), allowing us to simulate up to 40 Mb of chromatin. Prophase chromosomes are modeled as arrays of consecutive loops of exponentially distributed length and random genomic locations, emanating from a flexible scaffold, as would result from a loop extrusion process (48). The loop array is further condensed by imposing poor solvent conditions to the density observed in electron microscopy (one nucleosome per 11×11×11nm cube, i.e. ~40% volume fraction) (49), while preserving the overall cylindrical shape of the chromosome (Fig. 2D). We systematically varied two parameters: the average loop size and the linear loop density along the chromosomal scaffold (Fig. 2E). For all combinations, we generated equilibrium conformations, simulated a Hi-C experiment, and evaluated its ability to reproduce P(s) curves from Hi-C data for different time points during prophase (Fig. 2E–H).

These polymer models can accurately reproduce P(s) (20 kb<s<4 Mb) for all prophase time points, in agreement with the prediction of the coarse-grained model (Fig. 2C). The best matching models for later prophase time points, when sister chromatids are separate and lie side-by-side, have gradually increasing average loop size: from 40–50 kb at t = 5 minutes to ~60–70 kb at t = 10 minutes (Fig. 2H), reproducing the gradually shifting position of the drop-off from 2 to 3.5Mb (i.e. increase of the layer size), while maintaining about the same ~50 loops per layer and ~250 loops per µm. These results are consistent with a model where loop arrays are formed early in prophase, and loop sizes grow gradually, e.g. by merging smaller adjacent loops (25). Thus, both coarse grained modes and polymer simulations indicate that by late prophase, chromosomes fold as dense arrays of loops, with consecutive loops positioned with correlated radial orientations.

Prometaphase spirals

A striking feature of prometaphase Hi-C data is the appearance of the second diagonal band, which appears as a distinct peak on the P(s) curves (Fig. 3A). This feature cannot be explained by interactions between sister chromatids, as these become minimal in prometaphase, and simulations show that no amount of overlap between sisters gives rise to such periodic pattern in interactions (Fig. S14). As argued above, periodic interactions seen by Hi-C are most readily explained by a helical organization of mitotic loop arrays, which has been observed microscopically (6, 10, 40, 50). Two classes of chromosome architecture can give rise to periodicity in contact frequencies: an “external” helix when the whole chromosome is folded into a solenoid (50) (the solenoid model), and a “staircase” model in which consecutive loops wind in a helical order around a centrally located scaffold (“internal” helix). We note that by “scaffold” we do not necessarily imply a solid integrated structure stretching from one end of the chromosome to the other. The “scaffold” could equally be a dynamic association of smaller complexes linked by DNA that pack with helical symmetry. By modeling we examined these classes of architectures and the continuum of models in between them.

To explore whether an internal helix can arise through reorganization of loop orientations, while preserving the cylindrical morphology of the whole chromosome, we extended our coarse-grained prophase model (Fig. 2B) by adding a preferred angular orientation for each loop: (1) as in prophase, the orientation of each loop is correlated with its neighbors; (2) these loops have preferred, but not fixed, orientations that follow a helical path, thus winding around the chromosomal scaffold (Fig. 3B). Loops in this spiral staircase model follow an angular Ornstein-Uhlenbeck random walk with bias toward preferred positions, and P(s) can be found analytically (51) (Supplemental Material, section “Coarse-grained model of contact probability decay in mitotic chromosomes – loops with spiral staircase orientation”). This coarse-grained model yields a P(s) curve that closely follows the experimental prometaphase P(s) and displays both the P(s)~s^−0.5 decay and the narrow peak corresponding to the second diagonal band (Fig. 3C). These results indicate that, (i) the emergent second diagonal band in Hi-C data can result from a spiral organization, and (ii) such organization can arise from preferred orientations of loops around the central scaffold. Further support for a centrally located spiraling scaffold is provided by analysis of chromosome shape and SMC2-mAID-GFP, CAP-H-mAID-GFP or CAP-H2-mAID-GFP localization along mitotic chromosomes from colchicine arrested DT40 cells (Fig. 3B, Fig. S15, Fig. S16). We observe a pattern of condensin localization that is consistent with a helical path of the scaffold.

Detailed polymer modeling allowed us to explore a broader range of architectures, with both external and internal helices, and to obtain quantitative estimates of loop sizes and other aspects of organization. Two aspects of the prometaphase organization must be captured by any model: (i) a higher linear density of chromatin of up to 50–70 Mb/µm, necessitating an evolution of the loop architecture; and (ii) spiraling of the scaffold. The higher density of loops can be achieved by a nested loop organization where several smaller (inner) loops are organized consecutively within each larger (outer) loop whose bases form the central axis (Fig. 3D). The presence of nested loops is an essential feature for prometaphase models, as models with a single layer of loops could not reproduce Hi-C P(s) curves even when other parameters were varied (Fig. S14B). To model helical architecture we made the scaffold follow a helical path in 3D, while allowing loops to adopt their equilibrium conformations within an otherwise cylindrical chromosome (Fig. 3D).

We systematically varied the model parameters, such as geometry of the spiral scaffold and loop sizes (Fig. 3E). This also probed different lengths and widths of chromosomes as the volume density was kept constant. For t = 30 minutes the best agreement was achieved for a relatively narrow internal spiral staircase-like scaffold (R=30–60nm) (Fig. 3F–G). This spiral is much more narrow than the ~300nm diameter of the chromatid, and has a small pitch (the height of one turn, 100–200nm) (Fig. 3H). Interestingly, this spiral arrangement of loop bases can achieve helical winding of loops that reproduces the second diagonal in the interaction maps and the peak on the P(s) curves for t = 15, 30 and 60 minutes (Fig. 3F–H). Wider spiraling of the scaffold (Fig. 3G–III) approximating external helix architectures (50) failed to accurately reproduce P(s) (Fig. S14C). Taken together, coarse-grained and polymer models that agree with Hi-C data overwhelmingly support the spiral staircase (“internal” helix) architecture of the scaffold (I and II in Fig. 3G), and subsequent helical winding of loops.

Fitting consecutive time points probed by Hi-C, we found that the linear chromatin density of the best-matching models continued to grow throughout prometaphase, in agreement with the observed steadily shortening of mitotic chromosomes (Fig. 1A, ​,3H,3H, Fig. S18). Simulations show that shifting the peak in P(s) to larger genomic distances representing the second diagonal in consecutive time points (Fig. 3A) can be achieved by increasing the radius of the helical scaffold from 30 to 100 nm, the radius of the chromatid from 300 to 360 nm, and increasing the pitch from 100 to 250 nm, while maintaining a constant outer and inner loop size (~400 kb and ~80 kb respectively) (Fig. 3H). These changes lead to the increase in amount of DNA (Mb) per turn of the spiral from ~3 Mb in early prometaphase up to ~12 Mb by late prometaphase.

Comparison of the dimensions of our best models with direct microscopic measurements of prometaphase chromosomes prepared according to the Hi-C protocol reveals good agreement between experiment and predictions (Fig. S17).

Condensins are critical for prophase chromosome morphogenesis

To determine the role of condensin complexes in chromosome morphogenesis, we fused a minimal auxin inducible degron domain (mAID) to SMC2 (Supplemental Materials, Table S3). In the presence of the plant F-box protein osTIR1, addition of auxin induces rapid proteasome-dependent degradation of the SMC2-mAID protein, thus disrupting both condensin I and II complexes (18, 26, 27). Incubation of cells for 3 hours in the presence of auxin during the 1NM-PP1-induced G₂–arrest (Supplementary Materials) reduced SMC2 levels to <5% (Fig. S19). Remarkably, this did not affect global chromosome organization as compartments and TADs were comparable to those in WT G₂ arrested cells (Fig. 4A, Fig. S20). Cells entered prophase rapidly after washout of 1NM-PP1, and the onset of NEBD, as indicated by DAPI staining, occurred as in wild type at ~7–10 minutes (Fig. S21).

Open in a separate window

Fig. 4

Defects in chromosome morphogenesis in condensin depleted cells

(A–C). Hi-C interaction frequency maps (binned at 100 kb) for chromosome 7 at indicated time points (top right in each heatmap) after release from G₂ arrest. The first plot below each Hi-C interaction map displays the compartment signal (Eigenvector 1). The bottom graph shows the insulation score (TADs; binned at 50 kb). (A) SMC2-mAID cells were treated with auxin for three hours prior to release from G₂ arrest to deplete SMC2. SMC2+: Hi-C interaction map for G₂-arrested cells prior to auxin treatment. SMC2-: Hi-C interaction map for G₂-arrested cells after three hours of auxin treatment. (B) Hi-C data for CAP-H2-mAID cells treated for three hours with auxin prior to release from G₂ arrest to deplete CAP-H2. (C). Hi-C data for CAP-H-mAID cells treated for three hours with auxin prior to release from G₂ arrest to deplete CAP-H.

Chromosomes in SMC2-depleted cells did not form well-resolved chromatids as cells progressed to prometaphase, confirming previous observations (Fig. S21) (18, 52–54). Chromatin in such cells lacks functional condensin (55, 56), but nonetheless achieves a normal degree of chromatin compaction despite the absence of individualized chromosomes (57). FACS analysis confirmed that these cells are incapable of normal mitotic exit. They ultimately undergo mitotic slippage, forming tetraploid interphase cells (Fig. S2).

Hi-C analysis revealed that in the absence of SMC2, interphase compartments and TADs were still present and largely unaffected by late prophase, at a time when they were completely disassembled in WT (t = 10 minutes, Fig. 4A, Fig. S20A). NEBD and spindle assembly did occur; indicating cells progressed to physiological prometaphase. In prometaphase (t = 45 minutes and t = 75 minutes), compartments and TADs progressively weakened, but remained detectable (Fig. 4A, Fig. S22–23). No second diagonal, characteristic for WT prometaphase ever appeared in Hi-C maps (Fig. 4A), instead P(s) curves show little change from G₂ (Figs. S20A, S24–26). Preferential A-to-A interactions and B-to-B interactions became progressively weaker (Fig. S20A, S22A). Analysis of the variation of the insulation score along chromosomes indicated that TAD boundaries were reduced in strength but not eliminated (Fig. 4A, Fig. S23). Further, removal of cohesin (SMC1/3) and CTCF from chromatin, as assessed by chromatin enrichment for proteomics (ChEP) (58) was delayed and reduced compared to WT (Fig. S4D). This may explain the incomplete loss of TAD boundaries. Combined, these data reveal that condensin is not required for TAD and compartment architecture during interphase. In its absence, mitotic chromatin is compacted but chromosomes do not become individualized or acquire the normal mitotic morphology, while partially preserving elements of interphase architecture. This indicates (i) that compaction and formation of rod shaped mitotic chromosomes are two separate processes, as assumed by our model; and (ii) a critical role for condensin is in the formation of proper morphology and internal organization of mitotic chromosomes, and in disassembly of the interphase architecture (59).

Condensin I and II play distinct roles in chromosome morphogenesis

Next, we determined the roles of condensin I and II separately. Condensin I and II bind chromatin independently (52, 56, 60, 61), and recent in vitro mitotic chromosome assembly experiments show that they can act independently (62). Therefore, depletion of one condensin complex is unlikely to affect the other, though we cannot rule out more subtle interplay between the complexes. We fused auxin inducible degron domains to the condensin II-specific kleisin CAP-H2 (CAP-H2-mAID) or the condensin I-specific kleisin CAP-H (CAP-H-mAID) in CDK1as DT40 cells (Supplemental Methods). Addition of auxin lead to >95% protein depletion in G₂–arrested CAP-H-mAID or CAP-H2-mAID cells (Fig. S19). Cells were then released from the G₂ block and chromosome conformation was determined by microscopy and Hi-C as cells progressed through mitosis.

Depleting either condensin I or II alone led to less severe phenotypes than depleting both together (Fig. S21). In contrast to cells lacking both condensin I and II (SMC2-mAID), these cells exited mitosis within 3 hours after entry into prophase (Fig. S2).

Comparison of Hi-C interaction matrices (Fig. 4BC, Fig. S20BC) and P(s) curves (Fig. 5A, 5B, Fig. S24BC, S25–26) for CAP-H and CAP-H2 depleted cells in late prometaphase (t = 30 minutes and 60 minutes) shows that they capture different aspects of the WT architecture. The P(s) curve for CAP-H2-depleted cells, where only condensin I remains active, matches that of the intra-layer organization of WT up to ~6 Mb, and lacks the second diagonal band (Fig. 5B). The P(s) curve for CAP-H depleted cells (active condensin II), matches that of WT only for the long-range organization (6–20 Mb), including the second diagonal band (Fig. 5B). CAP-H depleted cells have a much lower contact frequency between loci separated less than 6 Mb than do WT and CAP-H2 depleted cells. Thus, condensin I and II play distinct roles at different structural levels, in mitotic chromosome morphogenesis, providing a mechanistic explanation for earlier microscopic studies (60, 61, 63, 64).

Open in a separate window

Fig. 5

Distinct roles for condensin I and II in mitotic chromosome formation

(A) Genome-wide curves of contact frequency P(s) vs genomic distance s, normalized to unity at s=100 kb. The curves are P(s) derived from Hi-C data obtained from CAP-H2-depleted (left panel) and CAP-H-depleted cells (right panel), at t = 7–60 minutes after release from G₂ arrest. Dashed line indicates P(s) = s^−0.5. (B) Overlayed P(s) curves of WT, CAP-H-and CAP-H2-depleted chromosomes show independent contributions of two condensin complexes to short- and long-distance contacts. (C) Polymer models of CAP-H2 (top) and CAP-H (bottom) depleted chromosomes. Top: depletion of CAP-H2 is modeled via removal of outer loops and relaxation of the helix.. Bottom: depletion of CAP-H is modeled via removal of the inner loops, while preserving the helical arrangement of the scaffold. Condensin II loop anchors are shown in red, condensin I loop anchors are shown in blue. (D) P(s) derived from late prometaphase CAP-H2 depletion Hi-C experiments (red line) and the three best fitting polymer models (grayscale lines). The average loop size and linear density of loops along the chromosome axis are indicated. (E) The average loop size and linear DNA density of the 3 best-fitting models of CAP-H2-depleted chromosomes at different time points. (F) P(s) derived from late prometaphase CAP-H depletion Hi-C experiments (blue line) and the best fitting polymer models with and without nested inner loops (grayscale lines). The average size of outer and inner loops, the length of a helix turn in Mb and the helical pitch are indicated. (G) The best fitting P(s) predictions by the staircase coarse grained model for late prometaphase CAP-H depletion Hi-C experiments at t = 30 minutes after release of G₂ arrest (gray lines; experimental P(s): red lines. Top: loop size is 200 kb, bottom: loop size is 400 kb.

Helical winding during prometaphase requires condensin II

In condensin II-depleted cells, both A- and B-compartments and TADs were lost starting around the prophase-prometaphase transition (t = 10 – 15 minutes; Fig. 4B, Fig. S20, Fig. S22). In late prometaphase (t = 30 – 60 minutes), chromosomes in these cells were longer and narrower than WT chromosomes, as previously observed (60, 63, 64) (Fig. S21). P(s) curves for t = 10 and 15 minutes (early prometaphase) resembled those in WT for late prophase (t = 10 minutes; compare Fig. 5A with Fig. 2A), displaying a mild decay followed by a steep drop that is characteristic for a densely packed loop array (Fig. 2B). Most strikingly, CAP-H2 depletion prevented emergence of the second diagonal band in prometaphase in Hi-C contact frequency maps and P(s) plots (Fig. 4B, ​,5B;5B; Fig. S20B, S24B).

The close similarity between CAP-H2 prometaphase and WT prophase Hi-C, and lack of the second diagonal, allowed us to model CAP-H2 chromosomes as a prophase-like array of a single layer of loops emanating from a flexible, non-helical scaffold. By systematically varying the loop size and the degree of linear compaction, we obtained excellent agreement with experimental P(s) curves for ~40–60 kb loops, and a linear density of 15 Mb/µm for all prometaphase time points (Fig. 5D, 5E). This linear density is 3–4 times smaller than that of WT prometaphase chromosomes (50–70 Mb/µm). These simulations indicate that in the absence of condensin II, prometaphase chromosomes form extended prophase-like loop arrays and do not progress to further longitudinal shortening and helical winding.

A mitotic chromosome morphogenesis pathway

The data and modeling presented here suggest a chromosome morphogenesis pathway by which cells convert interphase chromosome organization into compacted mitotic chromosomes (Fig. 6). Together, our imaging and Hi-C data, coarse-grained models and polymer simulations, and previous observations (11) reveal that upon entry into prophase, interphase features such as compartments and TADs are lost within minutes in a condensin-dependent process and by late prophase, chromosomes are organized as radial loop arrays. The mechanism by which TADs and compartments are lost is not known. Our data show that condensin is required. Additional contributing factors could include loss of CTCF and cohesin binding (Fig. S4B, S4C) and increased levels of loop extrusion that can erase boundaries even when CTCF is still bound (Supplemental Materials, Fig. S27). Importantly, activation of the mitotic kinase cascade is not sufficient to disassemble interphase chromatin organization without the action of condensin.

Open in a separate window

Fig. 6

A mitotic chromosome morphogenesis pathway

In prophase, condensin II compacts chromosomes into arrays of consecutive loops and sister chromatids split along their length. The scaffold of condensin II-mediated loop bases is indicated in red. Upon nuclear envelope breakdown and entry into prometaphase, condensin II-mediated loops become increasingly large as they split into smaller ~80 kb loops by condensin I. Chromosomes are shown as arrays of loops (only inner loops can be observed microscopically; top: cross-section, bottom: side view. For clarity, loops are indicated as separate entities pointing in one direction, while in reality loops are unstructured and can mix). The nested arrangements of centrally located condensin II-mediated loop bases and more peripherally located condensin I-mediated loop bases are indicated in red and blue respectively. During prometaphase central scaffold acquires a helical arrangement with loops rotating around the scaffold as steps in a “spiral staircase” (helical path of loops is indicated by arrows). As prometaphase progresses outer loops grow and the number of loops per turn increases and chromosomes shorten to form the mature mitotic chromosome.

Our models that achieve best agreement with Hi-C data show that during prophase, condensin II-dependent loops grow from 30–40 kb to 60 kb in size, leading to a ~2-fold increase in linear chromatin density from ~7 Mb/µm to 15 Mb/µm. Condensins at loop bases form a chromosomal scaffold (19, 62), which may be a dynamic, rather than static structure, and loops are arranged consecutively along it (one loop every ~5 nm of the axis). Interestingly, the radial arrangement of loops around the central flexible scaffold is not random, with consecutive loops projecting in similar directions i.e. with an angularly correlated arrangement.

Chromosomes shorten along their longitudinal axis and become wider during prometaphase. Our simulations show that condensin II loops continue to grow to 200–400 kb by 30 min and 400–700 kb by 60 min, accompanied by an increase in the linear chromatin density, which reaches 60 Mb/µm. However, two important reorganizations take place during prometaphase. First, large condensin II-mediated loops are subdivided into smaller 80 kb loops in a condensin I-dependent process, thus producing a nested loop arrangement with ~400 kb outer loops and ~80 kb inner loops. Second, the loop array acquires a helical arrangement as evidenced by the appearance of a second diagonal band in Hi-C maps for all loci and chromosomes. Models show that this helical arrangement of loops can be achieved if the scaffold forms a narrow helical “spiral staircase” inside an otherwise homogeneous cylindrical chromosome. Interestingly, the period, radius and pitch of this helix continue to grow through prometaphase, and this growth is to some extent restrained by condensin I. An emerging model of the prometaphase chromosome thus has a central helical scaffold formed by condensin II (62) that organizes 200–400 kb outer loops, that are further subdivided into 80 kb condensin Imediated inner loops in order to achieve a high volume density.

Comparison to previous and classical studies

While specific details of this model emerge from an unbiased fitting of models to the data, the emerging organization and its quantitative characteristics agree with earlier studies. First, the 60–80 kb sizes of the inner loops are remarkably similar to values suggested by an extensive survey of the literature (44), measurements from electron microscopy (6, 19) and Hi-C analysis of mitotic HeLa cells (8). Similarly, changes of linear density from prophase to prometaphase in the best models (from 15 Mb/µm to 50 Mb/µm) are consistent with prophase chromosomes being at least two-fold longer than metaphase chromosomes (11, 61).

Second, helical prometaphase chromosomes have long been observed in certain chromosome preparations (10, 37, 38, 40), and this has led to diverse models for how mitotic chromosomes are folded. Our analysis of Hi-C data indicates that the prometaphase chromosome is organized around a helical central region or scaffold: loops emanate with helical packing from a centrally located “spiral staircase” scaffold. Modeling shows that other helical arrangements of loop arrays, e.g. coiling of the entire loop array itself (40, 50, 65, 66), are not consistent with our Hi-C data.

Our helical scaffold/loop model unifies a range of models and observations made over the years. It explains how a helical chromatin packing arrangement can be achieved while scaffold proteins such as condensins and topoisomerase II are localized centrally (15–17), within a cylindrical chromatid that is not obviously helical when visualized with a DNA dye such as DAPI (67). Interestingly, by late prometaphase we estimate the height of one helical turn to be around 200 nm, which is also the size of the layer (12 Mb layer at linear density 60 Mb/um), and is consistent with microscope observations suggesting that consecutive genomic loci follow a helical gyre with a pitch of ~250 nm within the cylindrical shape of chromatids (68).

Supported by grants from the National Human Genome Research Institute (HG003143) and the National Institutes of Health Common Fund (DK107980) to J. Dekker and L. Mirny. J. Dekker is an investigator of the Howard Hughes Medical Institute. Work in the Mirny lab was also supported by the National Science Foundation (Physics of Living Systems, 1504942), and the National Institute of General Medical Sciences (GM114190). Work in the Earnshaw lab was funded by Wellcome, of which W.C.E. is a Principal Research Fellow (Grant 107022). The Wellcome Centre for Cell Biology is supported by core funding from the Wellcome Trust [203149]. M.T.K. was supported by JSPS KAKENHI Grant (16K5095), and research grants from Mochida Memorial Foundation for Medical and Pharmaceutical Research, SGH Foundation, The Sumitomo Foundation, and The Canon Foundation. Sorting of GFP positive and negative cells was performed in the Flow Cytometry facility, Institute of Immunology & Infection Research, Edinburgh with assistance of Dr. Martin Waterfall. We thank Dr. Hakan Ozadam for assistance with Hi-C data processing. All Hi-C data has been submitted to GEO and will be publicly available upon publication (accession number GSE102740). This paper is dedicated to the memory of Donald S. Coffey, who pioneered the study of chromatin loops in the nucleus.