The interplay between asymmetric and symmetric DNA loop extrusion
Abstract
Chromosome compaction is essential for reliable transmission of genetic information. Experiments suggest that ∼1000fold compaction is driven by condensin complexes that extrude chromatin loops, by progressively collecting chromatin fiber from one or both sides of the complex to form a growing loop. Theory indicates that symmetric twosided loop extrusion can achieve such compaction, but recent singlemolecule studies (Golfier et al., 2020) observed diverse dynamics of condensins that perform onesided, symmetric twosided, and asymmetric twosided extrusion. We use simulations and theory to determine how these molecular properties lead to chromosome compaction. High compaction can be achieved if even a small fraction of condensins have two essential properties: a long residence time and the ability to perform twosided (not necessarily symmetric) extrusion. In mixtures of condensins I and II, coupling twosided extrusion and stable chromatin binding by condensin II promotes compaction. These results provide missing connections between singlemolecule observations and chromosomescale organization.
Introduction
During mitosis in metazoan cells, each chromosome is linearly compacted ~1000fold from a ~1mmlong chromatin polymer globule into a ~1 µm array of chromatin loops (Paulson and Laemmli, 1977; Marsden and Laemmli, 1979; Earnshaw and Laemmli, 1983; Maeshima et al., 2005). This remarkable reorganization is primarily driven by the protein complex condensin (Hirano and Mitchison, 1994; Saitoh et al., 1994; Saka et al., 1994; Strunnikov et al., 1995; Hirano et al., 1997; Shintomi et al., 2015; Gibcus et al., 2018), which is one of a class of DNAbinding complexes known as structural maintenance of chromosomes (SMC) complexes. Condensin is believed to compact the chromatin fiber by an active process known as ‘loop extrusion’ (Yatskevich et al., 2019; Banigan and Mirny, 2020). In the loop extrusion model, a loopextruding factor (LEF), such as a condensin motor, binds the chromosome and progressively grows a DNA/chromatin loop by translocating along and processively collecting the nearby chromatin fiber (Riggs, 1990; Alipour and Marko, 2012; Goloborodko et al., 2016b). DNA loop extrusion by condensins (Ganji et al., 2018; Golfier et al., 2020; Kong et al., 2020; Kim et al., 2020) and other SMC complexes (Kim et al., 2019; Davidson et al., 2019; Golfier et al., 2020) has recently been observed in singlemolecule experiments in vitro. However, it has not been established how condensins with the properties observed in vitro can attain the high degree of linear compaction required for mitotic chromosome compaction in vivo.
In the first singlemolecule experiments that directly imaged loop extrusion, yeast condensins were observed to extrude DNA loops in an asymmetric, ‘onesided’ manner (Ganji et al., 2018). In this mode of loop extrusion, part of the condensin complex remains anchored to DNA (Kschonsak et al., 2017), while condensin extrudes DNA from one side of the complex into a loop (i.e., collecting DNA either upstream or downstream of the bound site, but not both) (Ganji et al., 2018). Importantly, this contrasts with most models for loop extrusion by condensin, in which each loopextruding factor performs symmetric ‘twosided’ extrusion, growing loops by gathering DNA/chromatin from both sides of the protein complex (Alipour and Marko, 2012; Goloborodko et al., 2016b; Fudenberg et al., 2017; Banigan and Mirny, 2020). Furthermore, theoretical arguments and computational modeling predict that the observed onesided activity is insufficient to generate the 1000fold linear compaction expected for metazoan mitotic chromosomes (Banigan and Mirny, 2019; Banigan et al., 2020; Banigan and Mirny, 2020).
Subsequent singlemolecule experiments with condensins revealed different and more diverse properties for loop extrusion, particularly in higher eukaryotes. Recent experiments show that condensins from human and Xenopus cells can perform both onesided and twosided loop extrusion (Golfier et al., 2020; Kong et al., 2020). However, twosided extrusion by Xenopus condensins proceeds asymmetrically rather than symmetrically (Golfier et al., 2020). Chromosome compaction by condensins performing asymmetric but twosided loop extrusion has not yet been systematically and quantitatively investigated.
An additional complication is that higher eukaryotes have two types of condensin, condensins I and II (Ono et al., 2003), each of which has different properties, including residence times and extrusion speeds (Gerlich et al., 2006a; Walther et al., 2018; Kong et al., 2020). Each of these condensins plays a role in mitotic chromosome compaction (Ono et al., 2003; Hirota et al., 2004; Ono et al., 2004; Green et al., 2012; Ono et al., 2017; Gibcus et al., 2018; Hirano, 2016; Kalitsis et al., 2017; Takahashi and Hirota, 2019), but the linear compaction abilities of mixtures of loopextruding condensins with different dynamic properties has not been systematically explored.
Using simulations and theory, we investigated whether asymmetric twosided extrusion or a mixture of one and twosided loopextruding factors (LEFs) with different dynamics can generate the high degree of linear compaction observed for mitotic chromosomes in vivo. We find that asymmetric twosided extrusion can eliminate unlooped gaps between neighboring LEFs and compact chromosomes $>1000$fold. Importantly, compaction can be achieved even with the relatively large asymmetries that are observed in vitro, provided that LEFs are twosided. We also perform simulations and develop a theory that show that mixtures of one and twosided LEFs can achieve high levels of compaction, provided that the twosided LEFs have sufficiently long residence times. Furthermore, the simulations suggest that mitotic chromosome compaction may require a tight coupling between stable chromatin binding and twosided extrusion by condensin II complexes, while condensins I and II that dynamically exchange may perform onesided extrusion. This result suggests that condensin II complexes may dimerize in vivo to promote chromosome compaction. These models provide the first demonstration of how loopextruding condensins with the properties observed in singlemolecule experiments could generate the linear compaction required to form metazoan mitotic chromosomes.
Model
In the model, LEFs representing SMC complexes perform loop extrusion on a polymer fiber representing the chromosome (Alipour and Marko, 2012; Goloborodko et al., 2016b; Fudenberg et al., 2016; Banigan et al., 2020). Each LEF is composed of two subunits or ‘sides’, which may have different translocation abilities, so the entire LEF may be either symmetric or asymmetric. A subunit may be either active or inactive. An active subunit processively translocates at speed $v$ along the polymer fiber, thus creating and enlarging the polymer loop between the subunits. In this work, inactive subunits are immobile. We refer to LEFs with two active subunits as ‘twosided’. LEFs with one active subunit and one inactive subunit are referred to as ‘onesided’. Each LEF subunit is assumed to act as a barrier to the translocation of other LEFs so that an active subunit cannot pass through another LEF subunit. Thus, pseudoknots or ‘Zloops’ (Kim et al., 2020) are prohibited; the scenario in which onesided LEFs may traverse each other has been considered previously (Banigan et al., 2020) (moreover, compaction by mixtures of one and twosided LEFs in that model would simply rescale the mean processivity compared to the pure onesided LEF model). LEFs bind with equal probability to any location on the polymer fiber, following a previous analysis of condensin localization that found condensin loading to be largely sequenceindependent (Gibcus et al., 2018). Furthermore, each onesided LEF has two possible binding orientations ($\leftarrow$ or $\to$), which determines the direction in which extrusion proceeds along the polymer. Following in vitro experiments (Ganji et al., 2018; Golfier et al., 2020) and lacking molecular evidence that binding to chromatin in a particular orientation could be favored, the extrusion orientation for each LEF is selected randomly. Each LEF stochastically unbinds at rate $k$, which releases the corresponding polymer loop. Altogether, these LEF dynamics lead to a dynamic statistical steady state in which loops formed by LEFs stochastically appear, grow, and vanish. The steadystate fold linear compaction, $\mathcal{F}\mathcal{C}$, is calculated from the fraction, $f$, of the fiber that is extruded into loops as $\mathcal{F}\mathcal{C}={(1f)}^{1}$ (Banigan and Mirny, 2019). This quantity primarily depends on the ratio, $\lambda /d$, of the processivity ($\lambda =pv/k$, where $p=1$ or 2 for one or twosided LEFs, respectively, and in mixtures, $\lambda $ denotes the populationaveraged processivity) to the mean separation ($d$) between LEFs (Goloborodko et al., 2016b; Banigan and Mirny, 2019). Further details and a public link to the simulation code are provided in the Materials and methods section.
Previous computational models of loop extrusion generally assumed that all LEFs have the same average unbinding rate and that all active LEF subunits translocate at the same average speed (Alipour and Marko, 2012; Sanborn et al., 2015; Goloborodko et al., 2016b; Fudenberg et al., 2016; Miermans and Broedersz, 2018; Banigan and Mirny, 2019; Banigan et al., 2020). Since experiments observed SMC complex dynamics that are contrary to these assumptions, we consider models in which these assumptions are relaxed. In the Results section, we first present results for systems with LEFs that perform asymmetric twosided loop extrusion. Second, we investigate mixtures of one and twosided LEFs with different extrusion velocities, mixtures of one and twosided LEFs with different residence times, and mixtures of only onesided LEFs with different residence times. Third, we specifically consider mixtures of LEFs with the properties measured for condensins I and II in experiments. Additionally, the simulation results for mixtures of one and twosided LEFs are explained by theoretical arguments, which are presented in detail in Appendix 1.
Results
Asymmetric twosided extrusion can linearly compact mitotic chromosomes
To determine the ability of asymmetric twosided loopextruding condensins (Golfier et al., 2020) to compact mitotic chromosomes, we simulated an asymmetric variant of the twosided loop extrusion model. In these simulations, each LEF has two active subunits. One of these is a fast subunit that extrudes at speed $v}_{\mathrm{f}\mathrm{a}\mathrm{s}\mathrm{t}$, while the other is a slow subunit that extrudes at speed $v}_{\mathrm{s}\mathrm{l}\mathrm{o}\mathrm{w}$. LEFs unbind from the chromatin polymer fiber at rate $k=1/\tau $, where $\tau $ is the mean residence time. Thus, the mean processivity (i.e. how large a loop a LEF can extrude before unbinding) is $\lambda =({v}_{\mathrm{fast}}+{v}_{\mathrm{slow}})\tau $. A prototypical trajectory, a schematic illustration, and an arch diagram are shown in Figure 1a.
Eliminating or avoiding unlooped gaps between LEFs is critical to achieving a high degree of linear compaction (Banigan and Mirny, 2019; Banigan et al., 2020). Onesided LEFs cannot close all gaps because one quarter of all pairs of neighboring LEFs are in a divergent orientation ($\leftarrow \to $); thus, they extrude loops by collecting chromatin from opposite directions, while leaving the chromatin between the LEFs unlooped (Banigan and Mirny, 2019). In contrast, at sufficiently high processivities ($\lambda $) and linear densities ($1/d$), symmetric twosided LEFs eliminate unlooped gaps (Goloborodko et al., 2016b; Banigan et al., 2020). Based on the idea of closing gaps between LEFs, we previously argued that asymmetric twosided extrusion could fully linearly compact mitotic chromosomes provided that the residence time is sufficiently long; in particular, we require ${\lambda}_{\mathrm{slow}}/d\gg 1$(Banigan and Mirny, 2019; Figure 1b).
In simulations, we controlled asymmetry by varying the relative extrusion speeds of the active subunits, quantified by the dimensionless ratio ${v}_{\mathrm{slow}}/{v}_{\mathrm{fast}}$. For all cases of asymmetric twosided extrusion (i.e. ${v}_{\mathrm{slow}}>0$) over the simulated range (${10}^{4}\le {v}_{\mathrm{slow}}/{v}_{\mathrm{fast}}<1$), chromosomes can be linearly compacted 1000fold, provided that ${\lambda}_{\mathrm{fast}}/d$ (and thus, ${\lambda}_{\mathrm{slow}}/d$) is sufficiently large (Figure 1c). Specifically, fold compaction, $\mathcal{F}\mathcal{C}$, grows rapidly for ${\lambda}_{\mathrm{slow}}/d>1$ (Figure 1—figure supplement 1a), as predicted. This occurs because gaps between LEFs are eliminated even if the gap is bordered by two slow LEF subunits (Figure 1b and Figure 1—figure supplement 1b). As shown in Figure 1c, 1000fold compaction can be achieved for all asymmetries of twosided LEFs, notably including the asymmetries and $\lambda /d$ in the range of expected physiological values (see below and Materials and methods).
We next determined whether 1000fold linear compaction can be achieved with the asymmetries (${v}_{\mathrm{slow}}/{v}_{\mathrm{fast}}$) observed in singlemolecule experiments (Golfier et al., 2020) and previously estimated values of $\lambda /d$ for condensin (Goloborodko et al., 2016b; Banigan et al., 2020). We performed three sets of simulations with asymmetric twosided LEFs (see Materials and methods): (1) simulations with a 50/50 mixture of asymmetric twosided LEFs and onesided LEFs, which has the same distribution of ‘symmetry scores’ ($S=({v}_{\mathrm{fast}}{v}_{\mathrm{slow}})/({v}_{\mathrm{fast}}+{v}_{\mathrm{slow}})$) as measured for Xenopus condensins by Golfier et al., 2020, (2) simulations with only asymmetric twosided LEFs that also reproduce the symmetry score distributions measured by Golfier et al., 2020, and (3) simulations with asymmetric twosided LEFs that can both grow and shrink loops according to our reanalysis of the experiments of Golfier et al., 2020.
These three scenarios lead to two qualitatively different outcomes. The mixture of asymmetric twosided LEFs and onesided LEFs (scenario 1) can achieve only ~60fold linear compaction (blue circles in Figure 1d). Because half of the LEFs are onesided in this scenario, a relatively large number of unlooped gaps remain in steady state, which limits linear compaction. In contrast, simulations with different distributions of asymmetric twosided extrusion (scenarios 2 and 3) do not have this limitation, which results in $>1000$fold linear compaction for plausible values of $\lambda /d$ ($<1000$) (red squares and purple triangles in Figure 1d). Although a significant fraction of LEFs are highly asymmetric ($\ge 20\%$ of LEFs with ${v}_{\mathrm{slow}}<0.1{v}_{\mathrm{fast}}$), they are typically able to close gaps within their residence times. Thus, we conclude that even highly asymmetric twosided LEFs can close gaps and compact chromosomes, while a modest amount of strictly onesided LEFs significantly inhibits compaction.
Compaction by mixtures of one and twosided LEFs depends on their relative dynamic properties
Model for mixtures of one and twosided LEFs with different dynamics
Previous modeling predicts that a large majority ($>84\%$) of LEFs must perform twosided extrusion in order to sufficiently compact a mitotic chromosome (Banigan and Mirny, 2019), but experiments only observe symmetric twosided extrusion by 20–80% of condensins (Kong et al., 2020; Golfier et al., 2020). However, previous analyses of mixtures of one and twosided LEFs made a potentially important simplifying assumption; they considered only mixtures in which every active subunit translocates along the chromatin fiber at the same speed and every LEF has the same mean residence time (Banigan and Mirny, 2019; Banigan et al., 2020). In contrast, experimental measurements indicate that the condensins I and II, both of which compact mitotic chromosomes (Ono et al., 2003; Hirota et al., 2004; Ono et al., 2004; Green et al., 2012; Ono et al., 2017; Gibcus et al., 2018; Hirano, 2016; Kalitsis et al., 2017; Takahashi and Hirota, 2019), have different mean residence times (Gerlich et al., 2006a; Walther et al., 2018) and extrusion speeds (Kong et al., 2020).
In FRAP experiments, condensin I and condensin II have markedly different residence times on mammalian mitotic chromosomes. Condensin I typically remains bound to chromosomes for 2–3 min. Condensin II, in contrast, exhibits two types of turnover kinetics; 15–40% of condensins have a mean residence time of 5–8 min, while the remaining complexes are stably bound for longer durations (Gerlich et al., 2006a; Walther et al., 2018). Furthermore, a recent estimate based on HiC analysis and computational modeling suggests a 2hr residence time for condensin II (Gibcus et al., 2018).
Condensins I and II also have different extrusion velocities in vitro. Recent singlemolecule experiments (Kong et al., 2020), observed that loops extruded by condensin I grow at approximately twice the speed of those extruded by condensin II. Intriguingly, different extrusion speeds are also observed for loopextruding cohesins, depending on whether they perform onesided or twosided extrusion (Golfier et al., 2020). Together, these results demonstrate that the dynamics of SMC complexes may depend on their associated proteins (as with condensin I versus condensin II) or their mode of extrusion (as with onesided versus twosided extrusion).
To determine whether experimentally observed fractions of twosided condensins can achieve 1000fold linear compaction, we developed simulation and theory models for mixtures of LEFs with different mean velocities and/or residence times. We primarily consider mixtures of onesided and twosided LEFs, where each population has a distinct residence time and extrusion velocity. Further details are provided in the Materials and methods section. Figure 2a shows schematic drawings of the LEFs and an arch diagram for an example system.
Theoretical analysis of mixtures of one and twosided LEFs with different dynamics
We developed a theoretical model to quantitatively predict the degree of compaction expected with mixtures of one and twosided LEFs with different residence times and extrusion speeds (denoted by subscripts 1 and 2, respectively; see Appendix 1 for the full theory). The theoretical analysis predicts that the ratio of the extrusion speeds, ${v}_{2}/{v}_{1}$, does not affect the maximum fold linear compaction, $\mathcal{F}{\mathcal{C}}_{\mathrm{max}}$, because the speeds do not affect gap formation (for simulations, see Appendix 2). However, the theory predicts differences between mixtures with very longlived twosided LEFs (${\tau}_{2}\gg {\tau}_{1}$) and mixtures with very shortlived twosided LEFs (${\tau}_{2}\ll {\tau}_{1}$).
In the scenario with longlived twosided LEFs (${\tau}_{2}>{\tau}_{1}$), the shortlived onesided LEFs act as transient barriers to extrusion by the twosided LEFs (Figure 2b, left panel). After a barrier unbinds, the twosided LEF can extrude beyond that barrier, potentially closing an unlooped gap between onesided LEFs. The presence of transient barriers reduces the effective speed of the longlived twosided LEFs, in turn reducing the effective processivity of the those LEFs to ${\lambda}_{2}^{\mathrm{eff}}$. We then numerically compute the fraction of the fiber that we expect to be compacted by a system with only twosided LEFs at processivitytodensity ratio ${\lambda}_{2}^{\mathrm{eff}}/{d}_{2}$ (where $d}_{2$ is the mean separation between twosided LEFs). The remaining fraction that is not compacted by the twosided LEFs is assumed to be ~90% compacted by the onesided LEFs (i.e. the shortlived onesided LEFs compact the remaining fiber ~10fold [Banigan and Mirny, 2019]). The theory predicts that mixtures with longlived twosided LEFs (${\tau}_{2}>{\tau}_{1}$) compact more effectively than mixtures with a single mean residence time (${\tau}_{2}={\tau}_{1}$).
In the scenario with shortlived twosided LEFs (${\tau}_{2}<{\tau}_{1}$), longlived onesided LEFs act as permanent barriers to extrusion by twosided LEFs (Figure 2b, right panel). In this case, the processivity of the twosided LEFs is effectively limited to the mean separation, $d}_{1$, between onesided LEFs. Following the previous calculation (see Appendix 1), we find that compaction by mixtures with shortlived twosided LEFs is lower than in mixtures with ${\tau}_{2}={\tau}_{1}$.
These theoretical limits, along with theory developed previously for populations of LEFs with a single residence time (${\tau}_{2}={\tau}_{1}$) (Banigan and Mirny, 2019), establish predictions for the simulations described below; theoretical results are shown in Figure 2d.
Longlived twosided LEFs enhance compaction by mixtures of LEFs
We next used simulations to test the prediction that increasing the residence time of the twosided LEFs relative to that of the onesided LEFs could increase compaction in mixtures of LEFs. We hypothesized that longlived twosided LEFs might be able to further compact chromatin by two mechanisms. First, increasing the residence time increases the effective processivity, and thus potentially, the loop size. Second, increasing the residence time additionally allows twosided LEFs to remain on the chromatin fiber longer than the onesided LEFs; consequently, each transient gap between onesided LEFs may eventually be extruded into a loop by a twosided LEF.
In both qualitative and quantitative agreement with the theory, the maximum fold linear compaction, $\mathcal{F}{\mathcal{C}}_{\mathrm{max}}$, in simulations increases with increasing twosided LEF residence times (Figure 2c–e, Figure 2—figure supplement 2a, and Figure 2—figure supplement 3). As predicted, unlooped gaps that are formed by pairs of neighboring onesided LEFs are shortlived, and thus, they are less common for larger ${\tau}_{2}/{\tau}_{1}$ (Figure 2—figure supplement 2b). Therefore, mixtures of LEFs can achieve $>1000$fold compaction with longlived twosided LEFs because almost all of the chromatin fiber can be extruded into loops by the twosided LEFs.
The compaction abilities of mixtures of one and twosided LEFs with different residence times is summarized by the phase diagram in Figure 2e. The maximum fold compaction, $\mathcal{F}{\mathcal{C}}_{\mathrm{max}}$, depends on both the composition, ${\varphi}_{1}$, and the relative residence times, ${\tau}_{2}/{\tau}_{1}$. With very longlived twosided LEFs (${\tau}_{2}/{\tau}_{1}>10$), twosided LEFs can extrude most of the fiber because gaps formed by the onesided LEFs are relatively transient and infrequent. In this case, 1000fold compaction can be achieved even with fairly large fractions, ${\varphi}_{1}$, of onesided LEFs; for example, with ${\tau}_{2}/{\tau}_{1}\approx 40$, up to ~70% of LEFs may be onesided. In contrast, with shortlived twosided LEFs (${\tau}_{2}/{\tau}_{1}<1$), a large fraction, ${\varphi}_{2}$, of twosided LEFs is required to achieve 1000fold compaction because the twosided LEFs are frequently impeded by the longlived onesided LEFs; therefore, many twosided LEFs are needed to fully extrude the gaps between onesided LEFs (e.g. ${\varphi}_{1}=0.7$ now results in $\mathcal{F}\mathcal{C}<25$, and $\mathcal{F}\mathcal{C}=1000$ requires ${\varphi}_{1}<0.16$). Between these limits ($1<{\tau}_{2}/{\tau}_{1}<10$), systems with moderate fractions of onesided LEFs (e.g. $0.16<{\varphi}_{1}<0.5$) can achieve 1000fold linear compaction. These results establish that mixtures of LEFs can fully compact chromosomes provided that high fractions of onesided LEFs are adequately offset by long residence times for twosided LEFs.
LEFs with the dynamics of condensins I and II can compact chromosomes
We next considered our results for relative residence times (${\tau}_{2}/{\tau}_{1}$) and fractions of onesided extruders (${\varphi}_{1}$) measured by and estimated from experiments. FRAP experiments and HiC experiments and modeling suggest a range of $2\le {\tau}_{2}/{\tau}_{1}\le 60$ for the ratio of condensin II to condensin I residence times (Gerlich et al., 2006a; Gibcus et al., 2018; Walther et al., 2018) (yellow region in Figure 2d and e); singlemolecule experiments with metazoan condensins suggest a fraction of onesided condensins in the range $0.2\le {\varphi}_{1}\le 0.8$ green box in Figure 2d and yellow region in Figure 2e (Golfier et al., 2020; Kong et al., 2020). These ranges suggest that 1000fold compaction can be achieved within plausible physiological ranges of ${\tau}_{2}/{\tau}_{1}$ and ${\varphi}_{1}$.
However, the actual experimental situation is considerably more complicated. In vitro singlemolecule experiments with human condensins I and II show that condensins of both types may be either onesided or twosided (Kong et al., 2020); thus, some onesided condensins may be longlived and some twosided condensins may be shortlived. To address this scenario, we simulated mixtures of LEFs with the properties of condensins I and II (Figure 3a, blue box). LEFs representing condensin I were shortlived, with residence time ${\tau}_{\mathrm{I}}$, and LEFs representing condensin II were longerlived, with residence time ${\tau}_{\mathrm{II}}=3{\tau}_{\mathrm{I}}$. Furthermore, each population of condensin is itself a mixture of onesided and twosided LEFs; respectively, 80% and 50% of condensin I and condensin II LEFs were twosided, as in singlemolecule experiments (Kong et al., 2020). To match experimental measurements in Xenopus and HeLa cells (Ono et al., 2003; Shintomi and Hirano, 2011; Walther et al., 2018), we assume 80% of LEFs are condensin I (as a result, 74% of all LEFs are twosided).
With this base model for mixtures of condensins I and II, we simulated chromosome compaction. We first found that the mixtures of condensins described above can generate only ~260fold linear compaction, less than the 1000fold compaction required for human mitotic chromosomes (blue circles in Figure 3b).
We then noted that in vivo, condensin II has two subpopulations with different residence times; 15–40% of all condensin II dynamically exchange with a mean residence time of 5–8 min, while the remaining 60–85% of condensin II complexes are stably bound for a much longer residence time (Gerlich et al., 2006a; Walther et al., 2018). Therefore, we simulated a modified condensin model in which 50% of the condensin II LEFs are more stably bound (‘extrastable’) with a mean residence time of ${\tau}_{\mathrm{II},\ast}=20{\tau}_{\mathrm{II}}=60{\tau}_{\mathrm{I}}$ (the underestimate of 50% stably bound is for simplicity; see scenario below) (Figure 3a, red box). In these simulations, condensin I and II mixtures can generate up to ~540fold linear compaction, still short of our expectation for mitotic chromosomes (red squares in Figure 3b).
To further increase the degree of compaction, we assumed that stably bound condensin II performs twosided extrusion, while dynamic condensin II performed onesided extrusion (Figure 3a, purple box). In this scenario, mixtures of condensins I and II can generate $>1000$fold linear compaction (purple triangles in Figure 3b). Extrastable condensin II LEFs form large loops, while shortlived condensins I LEFs form smaller, nested loops (Figure 3—figure supplement 1). These results hold for simulations modeling other cell types with higher or lower levels of condensin I relative to condensin II (Figure 3—figure supplement 2), which model mitotic chromosomes in other types of cells (Ono et al., 2003; Ohta et al., 2010). Altogether, the simulations demonstrate that a coupling between long residence times and twosided extrusion (and between shorter residence times onesided extrusion) can enhance the attainable degree of compaction, including in experimentally relevant scenarios.
Discussion
A key outstanding question for loopextruding SMC complexes is how predominantly asymmetric extrusion, observed in vitro, can generate the high degree of linear compaction observed for mitotic chromosomes in vivo. We previously argued that effectively twosided extrusion or a strong targeted loading bias is needed to compact and organize chromosomes (Banigan and Mirny, 2019; Banigan et al., 2020; Banigan and Mirny, 2020). Recent experiments provide evidence that condensins might perform twosided extrusion (Golfier et al., 2020; Kong et al., 2020), albeit not precisely in the manner envisioned in previous theoretical arguments (Banigan and Mirny, 2020). Our present work establishes how metazoan mitotic chromosomes can be linearly compacted 1000fold by condensins performing asymmetric twosided extrusion or by condensins in a predominantly onesided mixture (Figure 4).
First, LEFs performing asymmetric twosided extrusion as observed in Xenopus extracts (Golfier et al., 2020) could compact mitotic chromosomes 1000fold if their ‘slow sides’ extrude quickly enough to eliminate unlooped gaps (Figures 1 and 4b). Second, mixtures of one and twosided LEFs in which the twosided LEFs have relatively long residence times can linearly compact chromosomes 1000fold, even with large fractions of onesided LEFs (Figures 2 and 4c). Third, in order to achieve 1000fold compaction with mixtures of condensins I and II, we predict that stable chromatin binding by condensin II complexes (Gerlich et al., 2006a; Walther et al., 2018) must be coupled to twosided extrusion (Figures 3 and 4d).
Our results for asymmetric LEFs show that LEFs performing asymmetric twosided extrusion, as in Xenopus (Golfier et al., 2020), could fully compact mitotic chromosomes (Figure 1d). The magnitudes and distribution of asymmetries observed in vitro (${v}_{\mathrm{slow}}/{v}_{\mathrm{fast}}>{10}^{2}$) are quantitatively consistent with the condition that gaps between LEFs must be closed within a single residence time (${\lambda}_{\mathrm{slow}}/d>1$; Figure 1b–d and Figure 1—figure supplement 1). Thus, asymmetric twosided extrusion with one rapidly extruding side (${v}_{\mathrm{fast}}\sim 1\mathrm{kb}/\mathrm{s}$) and one slowly extruding side ($10\mathrm{bp}/\mathrm{s}<{v}_{\mathrm{slow}}<{v}_{\mathrm{fast}}$) could compact mitotic chromosomes (Figure 4b, right).
In mixtures of one and twosided LEFs, a longer residence time for twosided LEFs allows those LEFs to extrude the gaps between onesided LEFs into loops after the onesided LEFs unbind (Figure 2a and b and right panel of Figure 4c). Thus, 1000fold compaction can be achieved even if the fraction of onesided LEFs exceeds the ${\varphi}_{1}\approx 0.16$ threshold fraction previously predicted for systems with a single mean LEF residence time (Banigan and Mirny, 2019; Figure 2d and e). In systems with only onesided LEFs, differences in the mean residence time can enhance linear compaction, but only to a $\sim 90$fold limit, which is still too small (Appendix 1, Appendix 2, and Figure 2—figure supplement 4b and c). Importantly, for mixtures of one and twosided LEFs, 1000fold compaction may be achieved within ranges of residence times expected from experimental observations and computational modeling (Gerlich et al., 2006a; Walther et al., 2018; Gibcus et al., 2018; Figure 2d and e).
These results are particularly relevant for condensins I and II, which have different residence times. According to a recent singlemolecule experiment, both condensins I and II may perform onesided and twosided extrusion (Kong et al., 2020). As a result, a difference between the residence times of condensins I and II is insufficient to explain chromosome compaction by itself because some onesided condensins could be longlived (Figure 3a and left panel of Figure 4d). However, if twosided extrusion by condensin II is coupled to the extremely long residence times observed for a subpopulation of condensin II in vivo (Gerlich et al., 2006a; Walther et al., 2018), complete mitotic chromosome compaction may be achieved (Figure 3b, Figure 3—figure supplement 2, and right panel of Figure 4d).
These results demonstrate the importance of long residence times together with twosided extrusion – whether it be symmetric or asymmetric – to robust chromosome compaction. They have several possible implications for the molecular mechanisms of loop extrusion and chromosome organization by SMC complexes.
Mechanisms for asymmetric twosided extrusion
It is unclear how a condensin complex could perform asymmetric twosided extrusion, while still compacting chromosomes in vivo. For instance, diffusive sliding of one side combined with active translocation by the other is unable to form large DNA loops in vitro (Ganji et al., 2018), and it is insufficient to consistently eliminate chromatin gaps and achieve 1000fold compaction in simulations (Banigan et al., 2020). Instead, complete compaction requires directed translocation to collect chromatin from both sides of the condensin complex. Directed translocation could be governed by one or more of several mechanisms.
One possibility is that asymmetric twosided extrusion could occur due to asymmetry that is intrinsic to the complex. The degree of asymmetry of a condensin could be controlled by the species of the kleisin or Hawks (HEAT proteins associated with kleisins) associated with the complex. Consistent with this hypothesis, it has been observed that the yeast kleisin Brn1 and Hawk Ycg1 can act as a ‘safety belt’ that anchors condensin to DNA (Kschonsak et al., 2017). However, anchoring by the safety belt combined with condensin translocation (Terakawa et al., 2017) generates pure onesided loop extrusion (Ganji et al., 2018), which is incompatible with the high degree of mitotic chromosome compaction observed in metazoan cells (Banigan and Mirny, 2019). Conceivably, a looser safety belt combined with the ability of the complex to perform twosided extrusion could generate asymmetric twosided extrusion in higher eukaryotes.
Alternatively, complexes performing asymmetric twosided extrusion could be dimers of condensins, for which asymmetries could arise by several mechanisms. The extruding complex could contain condensins with two different protein compositions. Each side of the resulting complex might have different extrusion dynamics. However, this possibility is diminished by the strong tendency for particular kleisins and HAWKs to associate together (i.e., form either condensin I or II) (Ono et al., 2003) and the different patterns of spatial localization of these proteins (Ono et al., 2003; Shintomi et al., 2017; Walther et al., 2018). A more likely possibility is that the compositions of the condensins within a dimer are identical, but asymmetry arises if the components are differentially regulated by posttranslational modifications. This would be consistent with the observation of asymmetric twosided extrusion by Xenopus condensins in their native context (Golfier et al., 2020). Yet another possibility is that dimerization of condensins generates an asymmetric complex. This could arise due to chirality within the joined molecules or through differential conformational changes that are required to form the dimer. These scenarios are not mutually exclusive; several of these mechanisms could act together to generate asymmetric twosided extrusion.
Mechanisms for stability of twosided condensins
What could give rise to the predicted stability of twosided SMC complexes? One hypothesis is that condensin I performs onesided extrusion, while condensin II performs twosided extrusion (Banigan and Mirny, 2020). The mean residence time of dynamic condensin II is about three times longer than that of condensin I, and there is an immobile subpopulation of condensin II with a much longer residence time (Gerlich et al., 2006a; Walther et al., 2018). Thus, twosided complexes (condensin II in this scenario) would be longlived compared to onesided complexes (condensin I). Such a mixture of condensins I and II could generate 1000fold chromosome compaction (Figures 2 and 4c). Further supporting this hypothesis, yeast condensin, which is evolutionarily conserved as condensin I (Hirano, 2012), performs onesided loop extrusion in vitro (Ganji et al., 2018). Thus, it is appealing to think that condensin II is a twosided counterpart to onesided condensin I in metazoan cells.
However, singlemolecule experiments with condensins I and II suggest otherwise, at least for human cells (Kong et al., 2020). It has been observed that human condensins I and II can both perform either onesided or twosided extrusion in vitro. Moreover, twosided extrusion events are observed more frequently for condensin I as compared to condensin II. Interestingly, photobleaching experiments show that twosided loopextruding condensin I complexes are dimers of condensin I molecules (Kong et al., 2020). This observation suggests that perhaps individual condensins are onesided extruders, but they can dimerize to form twosided loopextruding complexes.
Taking the in vitro experiments (Kong et al., 2020) together with our simulation results (Figure 2e and Figure 3b), we propose that condensin II complexes dimerize to form twosided loopextruding complexes with very long residence times. Dimerization could facilitate twosided extrusion by combining two onesided complexes, while also stabilizing binding of condensin II to DNA by increasing (or otherwise modifying) the condensinDNA binding surface. This combination of effects would naturally generate a tight coupling between twosided extrusion and long residence times. Thus, dimerization could simultaneously allow twosided extrusion and increase the residence time. This combination would generate the necessary conditions for chromosome compaction by mixtures of one and twosided condensins.
Spatial organization by condensins I and II
Metazoan chromosomes are composed of large (~400 kb) chromatin loops with condensin II at their bases, with smaller (~80 kb) loops mediated by condensin I nested within (Gibcus et al., 2018; Walther et al., 2018). Consequently, condensin II is tightly localized to the central axis of the rodlike chromosome, while condensin I is localized along the axis in a broader pattern (Ono et al., 2003; Shintomi et al., 2017; Walther et al., 2018).
This hierarchy of loop nesting naturally emerges from the loop extrusion model for mixtures of LEFs with two different residence times. LEFs with longer residence times, $\tau $, have a larger processivities, $\lambda =v\tau $, and thus tend to form larger loops. In 3D, the bases of these loops localize along a central axis (Goloborodko et al., 2016a; Gibcus et al., 2018). This hierarchy emerges in mixtures of LEFs with the properties of condensins I and II (Figure 3—figure supplement 1), mixtures of one and twosided LEFs (Figure 2b, bottom left), mixtures of two populations of onesided LEFs (Figure 2—figure supplement 4e), and mixtures of two populations of twosided LEFs (Gibcus et al., 2018). While not all combinations of LEF symmetries and dynamics can fully compact mitotic chromosomes (twosided extrusion is required; Figure 4), differences in residence times can generate the nesting structure required for the patterns of spatial localization of condensins I and II observed in vivo. Consistent with this idea, the residence time of condensin II in vivo is longer than that of condensin I (Gerlich et al., 2006a; Walther et al., 2018).
Furthermore, our results suggest specific properties of condensins that could regulate mitotic chromosome morphology. Previously, it was shown that mitotic chromosome morphology can be controlled by the relative ratio of condensin I to condensin II (Shintomi and Hirano, 2011). Generally, condensin II is responsible for lengthwise ‘axial’ compaction along the central axis of a chromatid, while condensin I drives ‘lateral’ compaction, reducing the width of a chromatid (Ono et al., 2003; Shintomi and Hirano, 2011; Green et al., 2012; Bakhrebah et al., 2015; Shintomi et al., 2017; Hirano, 2016; Kalitsis et al., 2017). We find that morphology and the degree of linear compaction could alternatively be controlled by the relative residence times of condensins I and II and the fraction of condensins that perform onesided (instead of twosided extrusion). For example, a perturbation (such as a posttranslational modification) that increases the residence time of condensin II could lead to greater lengthwise compaction due to the larger loops that would be formed by condensin II. Similarly, perturbations that promote twosided extrusion, especially by condensin II, could also increase the degree of lengthwise compaction. Moreover, molecular perturbations affecting the coupling between symmetry and residence time could have a significant effect on compaction and morphology. These types of biomolecular perturbations would enable cells to regulate chromosome compaction without requiring global changes to levels of condensin expression.
Kinetics of compaction by condensins I and II
In the loop extrusion model, LEFs can linearly compact chromosomes by extruding chromatin into loops within approximately one LEF residence time. Chromosome morphology, as quantified by loop sizes, equilibrates over longer timescales of approximately 5–10 residence times (Goloborodko et al., 2016a; Goloborodko et al., 2016b; Banigan et al., 2020). In mixtures of LEFs with different dynamics, there are multiple timescales that could govern compaction kinetics. Linear compaction into a series of loops may still occur within ~1–2 residence times of the shortlived LEFs, but equilibration of large loops formed by the longlived LEFs will occur over several residence times of the longlived LEFs.
With the experimentally measured residence times of condensins I and II (Gerlich et al., 2006a; Walther et al., 2018), loop extrusion could compact chromosomes into loop arrays within a few minutes during prophase (Gibcus et al., 2018). However, loops should continue to merge and grow as mitosis progresses since the residence time for stably bound condensin II is comparable to the duration of mitosis (Gerlich et al., 2006a; Walther et al., 2018). Consequently, chromosomes should axially shorten and undergo 3D compaction as mitosis progresses, as observed in vivo (Nagasaka et al., 2016; Gibcus et al., 2018). Simultaneously, chromosomes should be further compacted by progressive loading of condensins, which increases throughout mitosis (Walther et al., 2018) and hypercompacts chromosomes when mitosis is stalled (Sun et al., 2018). Thus, even with mixtures of condensins with long residence and equilibration times, 1000fold linear compaction is achievable during mitosis.
Expectations for extrusion in other scenarios
The principles for mitotic chromosome compaction developed here are generalizable to other chromosome organization scenarios. In particular, we previously showed that the physical principles underlying mitotic chromosome formation by loop extrusion are relevant for topologically associated domain (TAD) formation in interphase (Banigan et al., 2020), where extrusion by the SMC complex cohesin can be paused by barriers such as CTCF (Sanborn et al., 2015; de Wit et al., 2015; Fudenberg et al., 2016; Busslinger et al., 2017; Nora et al., 2017; Wutz et al., 2017). Thus, we expect that asymmetric twosided LEFs and mixtures of LEFs might be able to form TADs. There are similar requirements for juxtaposition of bacterial chromosome arms by bacterial SMC complexes, but there also are additional constraints due to the specific loading site near the origin of replication (Banigan et al., 2020). Accordingly, asymmetric twosided LEFs and mixtures of LEFs with different residence times might be able to form TADs but not be able to juxtapose chromosome arms.
Topologically associated domains
We previously showed that much like mitotic chromosome compaction, the formation of major features of interphase chromosomes, such as TADs, ‘dots’, and ‘stripes’ requires avoiding unlooped gaps, either between LEFs or between LEFs and TAD boundaries. Onesided extrusion can form TADs and stripes by enhancing local chromatin contacts (Banigan et al., 2020), as observed in HiC experiments (Dixon et al., 2012; Nora et al., 2012; Sexton et al., 2012; Rao et al., 2014; Fudenberg et al., 2016; Vian et al., 2018; Barrington et al., 2019). However, dots (Rao et al., 2014; Krietenstein et al., 2020) can only be generated by ‘effectively twosided’ loop extrusion because such extrusion can reliably bring together TAD boundaries (e.g. convergently oriented CTCF binding sites [Rao et al., 2014; Guo et al., 2015; Sanborn et al., 2015; de Wit et al., 2015; Vietri Rudan et al., 2015]). In TADs, asymmetric twosided LEFs should be able to eliminate unlooped gaps if the slow side of each LEF is fast enough. Specifically, TAD boundaries could be brought together if the processivity, ${\lambda}_{\mathrm{slow}}$, of the slow side is larger than either the mean distance between LEFs ($d$) or the TAD size (${L}_{\mathrm{TAD}}$). We expect $\lambda ={\lambda}_{\mathrm{fast}}+{\lambda}_{\mathrm{slow}}\sim 1001000\mathrm{kb}$, $d\sim 100200\mathrm{kb}$, and ${L}_{\mathrm{TAD}}\sim 1001000\mathrm{kb}$, based on previous simulations (Fudenberg et al., 2016; Banigan et al., 2020), measurements of cohesin’s properties, (Davidson et al., 2019; Kim et al., 2019; Golfier et al., 2020; Gerlich et al., 2006b; Kueng et al., 2006; Tedeschi et al., 2013; Hansen et al., 2017; Wutz et al., 2017; Cattoglio et al., 2019; Holzmann et al., 2019), and HiC maps (Dixon et al., 2012; Nora et al., 2012; Sexton et al., 2012; Rao et al., 2014). These values suggest that asymmetric twosided loop extrusion by cohesin could generate TADs, dots, and stripes for moderate asymmetries (${v}_{\mathrm{slow}}/{v}_{\mathrm{fast}}>0.1$). Consistently, DNA loop extrusion by cohesin in vitro is largely symmetric (Davidson et al., 2019; Kim et al., 2019; Golfier et al., 2020).
Mixtures of one and twosided cohesins with different residence times should also be able to form TADs. We previously showed that mixtures with relatively high fractions of onesided LEFs (${\varphi}_{1}\sim 0.5$) could form TADs. Long residence times for the twosided extruders could enhance the ability of mixtures to bring together TAD boundaries; as in simulations of mitotic chromosomes, shortlived onesided LEFs would merely form transient barriers to twosided extrusion. Such a scenario, however, remains largely hypothetical since extrusion by cohesin is mainly twosided (Davidson et al., 2019; Kim et al., 2019; Golfier et al., 2020).
Juxtaposition of bacterial chromosome arms
Juxtaposition of bacterial chromosome arms can be achieved in a more limited set of loop extrusion scenarios. Loading of SMC complexes near the origin of replication breaks the translational symmetry of the system; thus, LEFs must extrude loops symmetrically (or nearly so) (Banigan et al., 2020).
Consequently, we do not expect that asymmetric twosided LEFs could produce the patterns observed in HiC maps of Bacillus subtilis and Caulobacter crescentus chromosomes (Umbarger et al., 2011; Le et al., 2013; Marbouty et al., 2015). There, symmetric twosided extrusion produces a secondary diagonal that is perpendicular to the main (selfcontact) diagonal (Umbarger et al., 2011; Le et al., 2013; Marbouty et al., 2015; Wang et al., 2017; Miermans and Broedersz, 2018; Banigan et al., 2020). For bacterial chromosomes, each asymmetric LEF would juxtapose sites separated by different genomic distances, ${s}_{1}\ne {s}_{2}$, from the loading site. Asymmetric juxtaposition by many LEFs would thus generate two secondary diagonals in HiC maps, where both diagonals would not be perpendicular to the main diagonal.
For similar reasons, mixtures of one and twosided LEFs with different residence times generally will not juxtapose bacterial chromosome arms. Each onesided LEF brings one chromosomal arm into contact with the loading site and interferes with juxtaposition by LEFs that bind subsequently (Banigan et al., 2020). Therefore, any substantial level of onesided extrusion would disrupt chromosomal arm juxtaposition. Such interference in mixtures of one and twosided LEFs might be partially mitigated in scenarios in which LEFs may traverse each other (e.g. form Zloops). However, this possibility requires further detailed investigation since loop extrusion with LEF traversal may lead to a variety of complicated bacterial chromosome structures in simulations and in vivo (Brandão et al., 2020) and is subject to ongoing investigation (Anchimiuk et al., 2020).
Conclusion
Two ingredients are essential for mitotic chromosome compaction by condensins: sufficiently long residence times and some amount of either symmetric or asymmetric (effectively) twosided extrusion. Strikingly, the presence of even a small (~20%) fraction of such condensins in a mixture with other, purely onesided condensins could be sufficient to achieve compaction under physiological conditions. It remains to be determined what mechanisms are responsible for different types of extrusion dynamics and what factors might facilitate a transition from onesided to twosided loop extrusion in some molecules in vivo. However, our analysis suggests that twosided extrusion should be tightly coupled to stable chromatin binding. In particular, twosided extrusion by stably bound condensin II is sufficient to linearly compact mitotic chromosomes in simulations. We thus hypothesize that condensins, particularly condensin II, may bind chromosomes as a dimer of condensin complexes in vivo. Such a dimer might have a longer residence time due to a larger proteinDNA interface, while also performing twosided loop extrusion via its two protein motors.
Further singlemolecule, biochemical, and structural studies could also help to understand how kleisins, Hawks, and posttranslational modifications might generate diverse SMC complex symmetries and dynamics, and thus, functions. Such experiments, together with the principles established by our models, could clarify how the molecular properties of loopextruding SMC complexes compact and organize chromosomes throughout the cell cycle.
Materials and methods
Simulations
Request a detailed protocolStochastic simulations of LEFs on a chromatin fiber are performed as previously described (Goloborodko et al., 2016b; Banigan et al., 2020) with adaptations as described in the Model section. The chromatin fiber is a onedimensional lattice of length $L=60000$ sites, each of which is taken to be $a=0.5$ kb. LEFs bind to chromatin at rate $k}_{b$, and the two LEF subunits initially occupy two adjacent lattice sites upon binding the fiber. Each active subunit of the LEF may translocate. Translocation occurs in a directed manner away from the sites originally occupied by the LEF unless otherwise noted. Simulations typically consist of $N=10002000$ LEFs (i.e. $d=L/N=3060$ lattice sites or $d=1530\mathrm{kb}$). Each simulation is run for a duration of ${t}_{\mathrm{total}}=400{\tau}_{\mathrm{longest}}$, where ${\tau}_{\mathrm{longest}}$ is the longest mean residence time. For each simulation, 100 data points are collected long after achieving steady state (which occurs after $t\approx 10{\tau}_{\mathrm{longest}}$) from the time interval $300{\tau}_{\mathrm{longest}}\le t\le 400{\tau}_{\mathrm{longest}}$. Each simulation is run at least twice. The fraction, $f$, of the fiber compacted into loops is thus measured for each set of parameters with a standard error that is $<5\%$ of the mean. The simulation code is publicly available at (https://github.com/mirnylab/one_sided_extrusion/tree/master/mitotic/; Banigan, 2020; copy archived at swh:1:rev:b27012e95d354e8deaac5bcfdcb3c36b375626ce).
Estimation of physiological values of $\lambda /d$
Request a detailed protocolThe physiological range of the ratio of the processivity to the mean separation is estimated as $10<\lambda /d<1000$, as calculated previously (Banigan et al., 2020). Processivity, $\lambda $, was estimated from experimental measurements of condensin’s extrusion speed of $\sim 1\mathrm{kb}/\mathrm{s}$ in vitro (Ganji et al., 2018; Golfier et al., 2020; Kong et al., 2020) and measured and estimated residence times of order 1–100 min in vivo and in vitro minutes (Gerlich et al., 2006a; Terakawa et al., 2017; Walther et al., 2018; Gibcus et al., 2018). Mean separation, $d$, was determined by measured linear densities of 1 per 10–100 kb in vivo (Takemoto et al., 2004; Fukui and Uchiyama, 2007; Walther et al., 2018).
Analysis of asymmetric extrusion
Request a detailed protocolFor the general estimate of the physiological range of asymmetries, the expected range of asymmetries is calculated from the symmetry scores measured by Golfier et al., 2020. The symmetry score is given by $S=({v}_{\mathrm{fast}}{v}_{\mathrm{slow}})/({v}_{\mathrm{fast}}+{v}_{\mathrm{slow}})$, where $S=0$ indicates perfectly symmetric twosided extrusion and $S=1$ indicates (completely asymmetric) onesided extrusion. Using experimentally measured symmetry scores, we calculate that ${v}_{\mathrm{slow}}>0.025{v}_{\mathrm{fast}}$ for two thirds of condensin loop extrusion events (green box in Figure 1c).
We also considered three scenarios of asymmetric extrusion to model experiments by Golfier et al., 2020. In the first scenario, all LEFs perform asymmetric twosided extrusion, but there is a distribution of asymmetries; half of the LEFs have symmetry scores, $S$, uniformly randomly selected from the interval $[0,0.8]$ and half of LEFs have $S$ randomly selected from $[0.8,1]$. In the second scenario, only half of the LEFs perform asymmetric twosided extrusion with asymmetries in the interval $[0,0.8]$; the other half of the population performs onesided extrusion (i.e., $S=1$).
In the third scenario, we calculated loopgrowing and loopshrinking velocities for condensins from singlemolecule experiments by Golfier et al., 2020. Ten trajectories from the experiments were smoothed with a SavitskyGolay filter with a second order polynomial and a window of 63 frames as in the previous analysis by Golfier et al., 2020. For each condensin trajectory, we then calculated the mean size of loopgrowing and loopshrinking steps and computed mean loopgrowing and loopshrinking speeds for each of the two sides. We simulated twosided asymmetric extrusion with these ten sets of velocities, with each set of velocities assigned to one tenth of the LEFs.
Appendix 1
Theoretical analysis of mixtures of LEFs with different dynamics
Mixtures of one and twosided LEFs
We theoretically analyzed the ability of mixtures of one and twosided LEFs with different mean speeds and residence times to linearly compact the chromatin fiber. To compute the maximum achievable compaction, we considered two scenarios beyond the case of uniform speeds and residence times considered previously (Banigan and Mirny, 2019): (1) the case in which twosided LEFs are longlived compared to onesided LEFs (i.e. ${\tau}_{2}\gg {\tau}_{1}$) and (2) the case in which twosided LEFs are shortlived compared to onesided LEFs (i.e.${\tau}_{2}\ll {\tau}_{1}$). We consider speeds $v}_{1$ and $v}_{2$ that are sufficiently fast for a LEF to extrude all nearby chromatin into a loop before the LEF unbinds (i.e. $\lambda /d\gg 1$).
In both scenarios, it is sufficient to calculate the fractions, $f}_{1$ and $f}_{2$, of the fiber that would be extruded, respectively, by the onesided LEFs and the twosided LEFs for their respective linear densities and tobedetermined effective processivities. This calculation allows us to compute the fraction of the fiber that remains in unlooped gaps:
The total fraction of the fiber that is compacted into loops is then given by:
while the fold linear compaction is:
From previous meanfield theoretical calculations, we know that onesided LEFs linearly compact a fraction ${f}_{1}\approx 0.895$ of the fiber that is accessible to them (leaving $1{f}_{1}\approx 0.105$ in unlooped gaps) (Banigan and Mirny, 2019). The fraction, $f}_{2$, compacted by the twosided LEFs is determined below.
For ${\tau}_{2}\gg {\tau}_{1}$, twosided LEFs extrude on a fiber loaded with transient barriers (the onesided LEFs), while onesided LEFs act to compact a fiber that has already been partially compacted (by the twosided LEFs) (Figure 2b, left). To compute the fraction, $f}_{2$, compacted by the twosided LEFs, we consider that each twosided LEF extrudes processively at speed $v}_{2$, but periodically encounters onesided LEFs that act as transient barriers. The barriers limit the mean speed of the twosided LEFs. Since the barriers are separated by a mean distance of $d}_{1$, disappear in mean time ${\tau}_{1}$, the twosided LEFs have an effective velocity of:
which is the average distance that a twosided LEF travels between barriers over the mean lifetime that each barrier is present. The resulting ratio of effective processivity (${\lambda}_{2}^{\mathrm{eff}}={v}_{2}^{\mathrm{eff}}{\tau}_{2}$) to mean separation between twosided LEFs is:
This ratio is larger for a longer relative twosided LEF lifetime (larger ${\tau}_{2}/{\tau}_{1}$) because the twosided LEFs outlast and extrude past a greater number of transient barriers. The ratio also grows for increasing distances between barriers (larger $d}_{1$) because twosided LEFs can freely translocate for greater distances along the fiber.
The fraction, $f}_{2$, of the fiber that is compacted by the twosided LEFs is equivalent to the fraction, ${f}_{2}({\lambda}_{2}^{\mathrm{eff}}/{d}_{2})$, of the fiber compacted in a system of twosided LEFs with a single residence time at the processivitytoseparation ratio ${\lambda}_{2}^{\mathrm{eff}}/{d}_{2}$. The fraction ${f}_{2}({\lambda}_{2}^{\mathrm{eff}}/{d}_{2})$ is obtained by simulations of twosided LEFs (Goloborodko et al., 2016b; Banigan et al., 2020).
In contrast to the scenario above, with ${\tau}_{2}\ll {\tau}_{1}$, twosided extrusion occurs on a chromatin fiber that has been partially compacted by onesided extrusion. We approximate this condition as twosided LEFs performing loop extrusion on a background of effectively permanent barriers (longlived onesided LEFs) (Figure 2b, right).
Again, to determine the total fraction, $f$, of the fiber that is linearly compacted into loops, we first compute the effective velocity of the twosided LEFs. The twosided LEFs translocate until encountering barriers, spaced $d}_{1$ apart, after which they are permanently stalled (because the barriers are longlived compared to the residence time of a twosided LEF). Therefore, the effective velocity of a twosided LEF is:
Now the ratio of effective processivity to separation is:
This ratio is independent of the speed and residence times of the twosided LEFs, because the onesided LEFs appear to be permanent barriers to the twosided LEFs, irrespective of $v}_{2$ and ${\tau}_{2}$ (provided that ${\tau}_{2}\ll {\tau}_{1}$). Once again, the fraction of the fiber compacted by the onesided LEFs is given by the fraction ${f}_{2}({\lambda}_{2}^{\mathrm{eff}}/{d}_{2})$, obtained from simulations of twosided LEFs.
Combining the above results, we find that the total fraction of the fiber that is compacted into loops is:
where ${\lambda}_{2}^{\mathrm{eff}}={v}_{2}^{\mathrm{eff}}{\tau}_{2}$ is set by either Equation 5 or 7, depending on ${\tau}_{2}/{\tau}_{1}$. The fraction, $f$, of the fiber compacted (Equation 8) can be expressed in terms of the fraction, ${\varphi}_{1}$, of onesided LEFs in the mixture by noting that ${d}_{1}/{d}_{2}=(1{\varphi}_{1})/{\varphi}_{1}$. Then, the ratio of the effective processivity to the mean separation is given by:
We interpolate between these two scenarios using the result from the previously developed meanfield theory (for ${\tau}_{2}={\tau}_{1}$, $f=(1({\varphi}_{1}/2{)}^{2}+\mathrm{ln}(4/({\varphi}_{1}{)}^{2}))/(1+\mathrm{ln}(4/({\varphi}_{1}{)}^{2}))$) (Banigan and Mirny, 2019).
In summary, we expect the maximum achievable linear compaction to grow rapidly with increasing ${\tau}_{2}/{\tau}_{1}$ for mixtures with relatively longlived twosided LEFs. Compaction can be depressed by large fractions of onesided LEFs. For mixtures with shortlived twosided LEFs, compaction is insensitive to ${\tau}_{2}/{\tau}_{1}$, but it decreases as the fraction, ${\varphi}_{1}$, of onesided LEFs is increased. Interestingly, the theory predicts that the ratio of relative velocities, ${v}_{2}/{v}_{1}$, does not alter the maximum achievable linear compaction. The theoretical results are shown in Figure 2d, Figure 2—figure supplement 1, and Figure 2—figure supplement 2, where they are compared to the simulation results.
Systems with two populations of onesided LEFs with different dynamics
Computing compaction from the effective processivities
We adapted the theory for mixtures of one and twosided LEFs to describe systems with onesided LEFs with two different residence times, ${\tau}_{\mathrm{S}}$ and ${\tau}_{\mathrm{L}}$, which are short and long, respectively. Similar to Equation 1, the total fraction of the fiber that remains in unlooped gaps is:
where ${f}_{\mathrm{L}}$ and ${f}_{\mathrm{S}}$ are the fractions of the fiber that would be compacted by, respectively, the longlived and shortlived LEFs at their respective linear densities and tobedetermined effective processivities.
As before, shortlived LEFs act as barriers to longlived LEFs, and thus limit the processivity of the longlived LEFs:
Therefore, the fraction of fiber that the longlived LEFs may compact is the fraction, $f}_{1$, compacted by onesided LEFs at ${\lambda}_{\mathrm{L}}^{\mathrm{eff}}/{d}_{\mathrm{L}}$. That is given by ${f}_{\mathrm{L}}={f}_{1}(\frac{{\tau}_{\mathrm{L}}{d}_{\mathrm{S}}}{{\tau}_{\mathrm{S}}{d}_{\mathrm{L}}})$, which is obtained from simulations of onesided LEFs (Banigan and Mirny, 2019; Banigan et al., 2020).
The fraction, ${f}_{\mathrm{S}}$, compacted by the shortlived onesided LEFs can be computed analogously to the calculation of compaction by shortlived twosided LEFs in the previous section. Longlived onesided LEFs are effectively permanent barriers to shortlived onesided LEFs, which have an effective processivity:
Then, the amount compacted by the shortlived LEFs is ${f}_{\mathrm{S}}={f}_{1}(\frac{{d}_{\mathrm{L}}}{{d}_{\mathrm{S}}})$.
Combining the expressions for ${f}_{\mathrm{S}}$ and ${f}_{\mathrm{L}}$ and noting that ${d}_{\mathrm{S}}/{d}_{\mathrm{L}}={\varphi}_{\mathrm{L}}/(1{\varphi}_{\mathrm{L}})$ (where ${\varphi}_{\mathrm{L}}$ is the fraction of longlived LEFs), we have:
As usual, the resulting fold compaction is given by $\mathcal{F}\mathcal{C}=1/(1f)=1/{f}_{\mathrm{gap}}$; thus, minimizing $f}_{\mathrm{g}\mathrm{a}\mathrm{p}$ in Equation 13 maximizes fold compaction, $\mathcal{F}\mathcal{C}$. Equation 13 shows the dependence of the maximum fold compaction, $\mathcal{F}{\mathcal{C}}_{\mathrm{max}}$, on several physical variables. The predictions of Equation 13 are shown as a function of ${\tau}_{\mathrm{L}}/{\tau}_{\mathrm{S}}$ for several values of ${\varphi}_{\mathrm{L}}$ in Figure 2—figure supplement 4b.
As in the case of one and twosided mixtures, the theory predicts that the maximum linear compaction by two populations of onesided LEFs is independent of their relative velocities. The velocities considered are large enough to close gaps between neighboring LEFs (if properly oriented), and therefore, altering the velocities does not change the rates of gap formation and closure (Figure 2—figure supplement 4d).
As the ratio of residence times, ${\tau}_{\mathrm{L}}/{\tau}_{\mathrm{S}}$, increases, the fraction, $f}_{\mathrm{g}\mathrm{a}\mathrm{p}$, of the fiber in gaps decreases. Increasing ${\tau}_{\mathrm{L}}/{\tau}_{\mathrm{S}}$ increases the effective processivity, ${\lambda}_{\mathrm{L}}^{\mathrm{eff}}$ (Equation 11). This effect increases the fraction, ${f}_{\mathrm{L}}$, of the fiber compacted by the longlived LEFs, and consequently, the total fraction, $f$, of the fiber that is compacted.
Interestingly, $\mathcal{F}{\mathcal{C}}_{\mathrm{max}}$ increases as the fraction, ${\varphi}_{\mathrm{L}}$, of longlived onesided LEFs decreases. This contrasts with mixtures of one and twosided LEFs, for which compaction increases as the fraction of longlived twosided LEFs increases (Equations 8 and 9). Decreasing ${\varphi}_{\mathrm{L}}$ increases the distance, ${\lambda}_{\mathrm{S}}^{\mathrm{eff}}$, that the shortlived LEFs may travel (i.e. the processivity, Equation 12), which increases their compaction ability. Moreover, while decreasing ${\varphi}_{\mathrm{L}}$ may decrease ${\lambda}_{\mathrm{L}}^{\mathrm{eff}}$, this decrease can be offset by increasing the ratio of residence times, ${\tau}_{\mathrm{L}}/{\tau}_{\mathrm{S}}$ (which increases ${\lambda}_{\mathrm{L}}^{\mathrm{eff}}$; Equation 13). Altogether, in the limit of very small but nonzero ${\varphi}_{\mathrm{L}}$ and very large ${\tau}_{\mathrm{L}}/{\tau}_{\mathrm{S}}$, the maximum fold compaction is $\mathcal{F}{\mathcal{C}}_{\mathrm{max}}\approx 90$ (each population of LEFs compacts $\sim 9.5$fold).
Computing compaction by counting unlooped gaps
To better understand the system with two types of onesided LEFs, we developed an alternative theory based on counting the number of unlooped gaps between loops in the limit of a very large disparity between LEF residence times (very large ${\tau}_{\mathrm{L}}/{\tau}_{\mathrm{S}}$). By counting the number of gaps and estimating their sizes by the mean separation between LEFs, we determine the fraction of the fiber that is not compacted into loops; using Equations 2 and 3, we compute the fold compaction.
First, because longlived LEFs remain on the chromatin fiber longer than the transient barriers (the shortlived LEFs), they can achieve their maximal $\sim 10$fold compaction when their processivitytodensity ratio, ${\lambda}_{\mathrm{L}}/{d}_{\mathrm{L}}$, is sufficiently large. We can thus use the previously developed mean field theory of onesided LEFs (Banigan and Mirny, 2019) to count the number of gaps, ${N}_{g,\mathrm{L}}$, between longlived LEFs:
where ${N}_{p,\mathrm{L}}$ is the number of longlived LEFs that form ‘‘parent’’ loops (as opposed to being nested as ‘‘child’’ loops). Furthermore, the previous mean field theory shows that:
where ${N}_{\mathrm{L}}$ is the total number of longlived LEFs.
The longlived gaps are partially extruded into loops formed by shortlived LEFs. Similar to above, we expect the number of shortlived gaps, ${N}_{g,\mathrm{S}}$, to be related to the number of shortlived parent loops, ${N}_{p,\mathrm{S}}$, by:
However, we expect that only a small fraction of the shortlived LEFs will reside within the longlived gaps and form parent loops, depending on the fraction, ${f}_{\mathrm{L}}$, compacted by the longlived gaps. Therefore, the number of shortlived parent loops is:
As previously calculated (Banigan and Mirny, 2019), ${f}_{\mathrm{L}}=(3+4\mathrm{ln}4)/(4+4\mathrm{ln}4)=0.895$. The gaps between shortlived LEFs are not extruded into loops, and they have size ${g}_{\mathrm{S}}={d}_{\mathrm{S}}$.
Finally, we account for the unlooped gaps between longlived LEFs and shortlived LEFs. We make the meanfield assumption that these gaps occur at the each of the boundaries of the longlived gaps with 50% probability and have size ${g}_{b}=d$, where $d=L/({N}_{\mathrm{L}}+{N}_{\mathrm{S}})$, where $L$ is the length of the chromatin fiber. These boundary gaps occur, on average, once per longlived gap, or equivalently, once per four longlived loops:
The total fraction of the chromatin fiber that is not compacted into loops is then given by:
Combining Equations 14–19, the total fraction of the fiber compacted into loops is:
where we have used the fact that the numbers of long and shortlived LEFs can be written in terms of the fraction, ${\varphi}_{\mathrm{L}}$, of longlived LEFs as ${N}_{\mathrm{L}}={\varphi}_{\mathrm{L}}({N}_{\mathrm{L}}+{N}_{\mathrm{S}})$ and ${N}_{\mathrm{S}}=(1{\varphi}_{\mathrm{L}})({N}_{\mathrm{L}}+{N}_{\mathrm{S}})$, respectively.
We may now determine the degree of compaction that can be achieved as a function of the fraction of longlived LEFs, ${\varphi}_{\mathrm{L}}$; the prediction is shown by the black curve in Figure 2—figure supplement 4c. In the limit of ${\varphi}_{\mathrm{L}}\to 0$, we find:
or 91fold compaction, close to the theoretical prediction in the previous section. In the limit of ${\varphi}_{\mathrm{L}}\to 1$, we find:
or 8.6fold compaction. This is close, but not precisely equal, to the 9.5fold compaction ($f=0.895$) for a single population of LEFs calculated previously (Banigan and Mirny, 2019). This discrepancy can be remedied by more carefully counting the boundary gaps, as described in the following subsection.
Refined counting of boundary gaps
The meanfield counting of the number of unlooped gaps (Equations 16 and 18) fails in two ways if the number of shortlived parent loops is comparable to or smaller than the number of longlived gaps.
First, in the above argument, we made a simplifying assumption in counting the number of boundary gaps between shortlived and longlived parent loops (Equation 18). We assumed that the number of boundary gaps equals the number of gaps between longlived parent loops, ${N}_{g,\mathrm{b}}={N}_{g,\mathrm{L}}$; this assumption is clearly violated if the number of shortlived parent loops falls below the number of longlived gaps (${N}_{p,\mathrm{S}}<{N}_{g,\mathrm{L}}$), which occurs for ${\varphi}_{\mathrm{L}}>{(2+\mathrm{ln}4)}^{1}=0.295$ in the above theory. This issue is corrected by:
Second, we assumed that ${N}_{g,\mathrm{S}}={N}_{p,\mathrm{S}}/4$ holds exactly, even though it clearly must be violated when there are fewer than two shortlived parent loops in a longlived gap (${N}_{p,\mathrm{S}}<2{N}_{g,\mathrm{L}}$), because a shortlived gap cannot be formed. This assumption is violated for ${\varphi}_{\mathrm{L}}>{(3+2\mathrm{ln}4)}^{1}=0.173$. We correct this issue by:
Altogether, we now find that the fraction of the fiber that remains uncompacted in gaps is:
After simplification, the fraction compacted is:
As shown by the gray curve in Figure 2—figure supplement 4c, this modified theory has similar behavior to the meanfield theory in the previous section. This refined counting argument better estimates the ratio of the total number of gaps, ${N}_{g,\mathrm{L}}+{N}_{g,\mathrm{S}}$, to total number of parent loops, ${N}_{\mathrm{\ell}}={N}_{p,\mathrm{L}}+{N}_{p,\mathrm{S}}$ (inset to Figure 2—figure supplement 4c). In addition, at ${\varphi}_{\mathrm{L}}=1$, this theory precisely reproduces the meanfield theory for a single population of onesided LEFs (Banigan and Mirny, 2019).
Appendix 2
Simulations of other types of mixtures of LEFs
Compaction by a mixture of LEFs does not depend on relative velocities
We speculated that rapid extrusion by twosided LEFs (${v}_{2}>{v}_{1}$) might generate large twosided loops, leading to nesting by the onesided LEFs that form smaller loops; consequently, there might be fewer gaps formed by divergently oriented ($\leftarrow \to $) LEFs. Consistent with this idea, previous modeling with two varieties of twosided LEF showed that LEFs with large processivities ($\lambda =v/k$) formed loops that were split by the extruded loops of LEFs with small processivities (Gibcus et al., 2018). On the other hand, theory (see Appendix 1) predicts that the relative extrusion velocities, set by the ratio ${v}_{2}/{v}_{1}$, do not set the maximum compaction. To test the theoretical predictions for linear compaction against the possibility of forming a hierarchy of loops, we simulated mixtures of one and twosided LEFs with different velocities, but identical mean residence times.
We find that the maximum achievable linear compaction is insensitive to the relative velocities, ${v}_{2}/{v}_{1}$, as predicted by the theory (Figure 2—figure supplement 1a). This is consistent with the previous theoretical model (Banigan and Mirny, 2019), which assumed rapid extrusion and closure of all possible gaps between LEFs in steady state. Altering the velocities of the LEFs does not alter the steadystate rate at which LEFs unbind from and rebind to the chromatin fiber; consequently, unlooped gaps are not eliminated any more effectively than in the scenario with a single mean velocity (Figure 2—figure supplement 1b). As a result, the maximum achievable fold compaction, $\mathcal{F}{\mathcal{C}}_{\mathrm{max}}$, for each fraction, ${\varphi}_{1}$, of onesided LEFs is given by meanfield theoretical limit calculated for mixtures of LEFs with uniform velocities (Banigan and Mirny, 2019; Figure 2—figure supplement 1a and c). Thus, irrespective of the relative velocities of one and twosided LEFs, a minimum fraction of twosided LEFs of ${\varphi}_{2}\approx 0.84$ LEFs is required to achieve 1000fold linear compaction.
Compaction by onesided LEFs is ineffective even with a longlived subpopulation
We explored whether having a two populations with different mean residence times would enhance linear compaction for a chromosome with only pure onesided LEFs. As described in Appendix 1, we adapted the theory for mixtures of one and twosided LEFs to systems with two populations of onesided LEFs, each with different residence times, ${\tau}_{\mathrm{S}}$ and ${\tau}_{\mathrm{L}}$ (denoting short and long). Once again, compaction is independent of the relative velocities of the two populations (Figure 2—figure supplement 4a). The theory predicts and simulations show that increasing the ratio of residence times, ${\tau}_{\mathrm{L}}/{\tau}_{\mathrm{S}}$, increases the maximum fold compaction, $\mathcal{F}{\mathcal{C}}_{\mathrm{max}}$ beyond the $\sim 10$fold limit for onesided LEFs with single mean residence time (${\tau}_{\mathrm{L}}/{\tau}_{\mathrm{S}}=1$) (Figure 2—figure supplement 4b and c). Interestingly, the compaction increases as the fraction of longlived LEFs, ${\varphi}_{\mathrm{L}}$, is decreased. This effect occurs because decreasing the number of longlived large loops increases the effective processivity of the shortlived LEFs (alternatively, the relationship can be understood by counting gaps and loops; see Appendix 1). Nevertheless, we find that linear fold compaction is still limited to a maximum of $\mathcal{F}{\mathcal{C}}_{\mathrm{max}}\approx 90$ (i.e. 9.5fold compaction from each population of onesided LEFs).
Data availability
Software used to perform simulations is publicly and freely available at https://github.com/mirnylab/one_sided_extrusion/tree/master/mitotic (copy archived at https://archive.softwareheritage.org/swh:1:rev:b27012e95d354e8deaac5bcfdcb3c36b375626ce/). Data analyzed from singlemolecule experiments was previously published as part of Golfier et al. eLife 9:e53885 (2020).
References

Selforganization of domain structures by DNAloopextruding enzymesNucleic Acids Research 40:11202–11212.https://doi.org/10.1093/nar/gks925

Disruption of a conserved CAPD3 threonine alters condensin loading on mitotic chromosomes leading to chromosome hypercondensationJournal of Biological Chemistry 290:6156–6167.https://doi.org/10.1074/jbc.M114.627109

Softwareone_sided_extrusion, version swh:1:rev:b27012e95d354e8deaac5bcfdcb3c36b375626ceSoftware Heritage.

Limits of chromosome compaction by loopextruding motorsPhysical Review X 9:031007.https://doi.org/10.1103/PhysRevX.9.031007

Loop extrusion: theory meets singlemolecule experimentsCurrent Opinion in Cell Biology 64:124–138.https://doi.org/10.1016/j.ceb.2020.04.011

CTCF binding polarity determines chromatin loopingMolecular Cell 60:676–684.https://doi.org/10.1016/j.molcel.2015.09.023

Architecture of metaphase chromosomes and chromosome scaffoldsJournal of Cell Biology 96:84–93.https://doi.org/10.1083/jcb.96.1.84

Formation of chromosomal domains by loop extrusionCell Reports 15:2038–2049.https://doi.org/10.1016/j.celrep.2016.04.085

ConferenceEmerging evidence of chromosome folding by loop extrusionCold Spring Harbor Symposia on Quantitative Biology. pp. 45–55.

Chromosome compaction by active loop extrusionBiophysical Journal 110:2162–2168.https://doi.org/10.1016/j.bpj.2016.02.041

Contrasting roles of condensin I and condensin II in mitotic chromosome formationJournal of Cell Science 125:1591–1604.https://doi.org/10.1242/jcs.097790

Condensins: universal organizers of chromosomes with diverse functionsGenes & Development 26:1659–1678.https://doi.org/10.1101/gad.194746.112

Distinct functions of condensin I and II in mitotic chromosome assemblyJournal of Cell Science 117:6435–6445.https://doi.org/10.1242/jcs.01604

Condensin, master organizer of the genomeChromosome Research 25:61–76.https://doi.org/10.1007/s1057701795530

Human cohesin compacts DNA by loop extrusionScience 366:1345–1349.https://doi.org/10.1126/science.aaz4475

Bacterial chromosome organization by collective dynamics of SMC condensinsJournal of the Royal Society Interface 15:20180495.https://doi.org/10.1098/rsif.2018.0495

Sister chromatid resolution is an intrinsic part of chromosome organization in prophaseNature Cell Biology 18:692–699.https://doi.org/10.1038/ncb3353

Spatial and temporal regulation of condensins I and II in mitotic chromosome assembly in human cellsMolecular Biology of the Cell 15:3296–3308.https://doi.org/10.1091/mbc.e04030242

Condensin II plays an essential role in reversible assembly of mitotic chromosomes in situMolecular Biology of the Cell 28:2875–2886.https://doi.org/10.1091/mbc.e17040252

DNA methylation and late replication probably aid cell memory, and type 1 DNA reeling could aid chromosome folding and enhancer functionPhilosophical Transactions of the Royal Society of London. Series B, Biological Sciences 326:285–297.https://doi.org/10.1098/rstb.1990.0012

Reconstitution of mitotic chromatids with a minimum set of purified factorsNature Cell Biology 17:1014–1023.https://doi.org/10.1038/ncb3187

The relative ratio of condensin I to II determines chromosome shapesGenes & Development 25:1464–1469.https://doi.org/10.1101/gad.2060311

Folding the genome into mitotic chromosomesCurrent Opinion in Cell Biology 60:19–26.https://doi.org/10.1016/j.ceb.2019.03.005

Cell cycledependent phosphorylation, nuclear localization, and activation of human condensinJournal of Biological Chemistry 279:4551–4559.https://doi.org/10.1074/jbc.M310925200

A quantitative map of human condensins provides new insights into mitotic chromosome architectureJournal of Cell Biology 217:2309–2328.https://doi.org/10.1083/jcb.201801048

Organization of chromosomal DNA by SMC complexesAnnual Review of Genetics 53:445–482.https://doi.org/10.1146/annurevgenet112618043633
Decision letter

Adèle L MarstonReviewing Editor; University of Edinburgh, United Kingdom

Kevin StruhlSenior Editor; Harvard Medical School, United States

Sarah SwygertReviewer; Fred Hutchinson Cancer Research Center, United States

Jan BruguésReviewer; Max Planck Institute of Molecular Cell Biology and Genetics, Germany
In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.
Congratulations, we are pleased to inform you that your article, "The interplay between asymmetric and symmetric DNA loop extrusion", has been accepted for publication in eLife.
Reviewer #2:
Banigan and Mirny revisit their modeling approach previously published in eLife (Banigan et al., 2020) to incorporate new experimental findings that in vivo condensins likely utilize asymmetric twosided loop extrusion in addition to the onesided and symmetrical twosided extrusion examined in the previous work. By adding variable residence time and asymmetric extrusion, the authors are able to model physiologically relevant compaction consistent with experimentally determined levels and properties of Condensins I and II. These findings represent a significant advancement over the previous work, the conclusions are welljustified, and a wellwritten and conceived Discussion gives a valuable perspective of the current state of the field. As a result, I support the publication of this manuscript in eLife without additional data.
Reviewer #3:
In this paper, Banigan and Mirny study the effects of different loopextrusion symmetries to the overall compaction of chromatin. This study is motivated by recent work showing that condensin loop extrusion is not onesided as shown purely in vitro, but rather displays a mixed population of symmetries. From a theory point of view, onesided loop extrusion cannot condense DNA sufficiently, so since the first in vitro demonstration of onesided loop extrusion it has not been clear how to marry theory and experiments with the requirements of DNA condensation in cells. Based on the results from Xenopus laevis egg extracts, the authors extend previous models of loop extrusion by allowing populations of loop extruders with different residence times, and active loop extruding subunits that translocate at the different speed. They show that the distribution of symmetries observed in these experiments (even though most of loop extrusion is onesided) can account for over 1000 fold condensation (as in cells). They also considered the effects of two condensins (I and II) corresponding to purely one and two sided extrusion. Interestingly they show that a longer residence time of a twosided condensin would also account for over 1000 fold condensation. As a consequence, their model suggests a testable prediction for the activities of condensin I and II. Overall I found the paper very clear and systematic, and their findings are an important contribution to understanding how loop extruding factors may contribute to the overall organization of chromatin, even though this work is an extension of their previous models. As such I think it fits as a research advances and support publication as is.
https://doi.org/10.7554/eLife.63528.sa1Article and author information
Author details
Funding
National Institutes of Health (U54DK107980)
 Edward J Banigan
 Leonid A Mirny
National Institutes of Health (U54CA193419)
 Edward J Banigan
 Leonid A Mirny
National Institutes of Health (GM114190)
 Edward J Banigan
 Leonid A Mirny
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Hugo Brandão, Jan Brugués, and Maxim Imakaev for helpful discussions and Job Dekker and John Marko for critically reading the manuscript. This work was supported by the NIH Center for 3D Structure and Physics of the Genome of the 4DN Consortium (U54DK107980), the NIH Physical SciencesOncology Center (U54CA193419), and NIH grant GM114190.
Senior Editor
 Kevin Struhl, Harvard Medical School, United States
Reviewing Editor
 Adèle L Marston, University of Edinburgh, United Kingdom
Reviewers
 Sarah Swygert, Fred Hutchinson Cancer Research Center, United States
 Jan Brugués, Max Planck Institute of Molecular Cell Biology and Genetics, Germany
Version history
 Received: October 1, 2020
 Accepted: November 30, 2020
 Accepted Manuscript published: December 9, 2020 (version 1)
 Version of Record published: January 8, 2021 (version 2)
Copyright
© 2020, Banigan and Mirny
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics

 1,325
 Page views

 235
 Downloads

 9
 Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading

 Chromosomes and Gene Expression
 Neuroscience
Mathys et al. conducted the first singlenucleus RNAseq (snRNAseq) study of Alzheimer’s disease (AD) (Mathys et al., 2019). With bulk RNAseq, changes in gene expression across cell types can be lost, potentially masking the differentially expressed genes (DEGs) across different cell types. Through the use of singlecell techniques, the authors benefitted from increased resolution with the potential to uncover cell typespecific DEGs in AD for the first time. However, there were limitations in both their data processing and quality control and their differential expression analysis. Here, we correct these issues and use bestpractice approaches to snRNAseq differential expression, resulting in 549 times fewer DEGs at a false discovery rate of 0.05. Thus, this study highlights the impact of quality control and differential analysis methods on the discovery of diseaseassociated genes and aims to refocus the AD research field away from spuriously identified genes.

 Chromosomes and Gene Expression
 Genetics and Genomics
Spermatogenesis in the Drosophila male germline proceeds through a unique transcriptional program controlled both by germlinespecific transcription factors and by testisspecific versions of core transcriptional machinery. This program includes the activation of genes on the heterochromatic Y chromosome, and reduced transcription from the X chromosome, but how expression from these sex chromosomes is regulated has not been defined. To resolve this, we profiled active chromatin features in the testes from wildtype and meiotic arrest mutants and integrate this with singlecell gene expression data from the Fly Cell Atlas. These data assign the timing of promoter activation for genes with germlineenriched expression throughout spermatogenesis, and general alterations of promoter regulation in germline cells. By profiling both active RNA polymerase II and histone modifications in isolated spermatocytes, we detail widespread patterns associated with regulation of the sex chromosomes. Our results demonstrate that the X chromosome is not enriched for silencing histone modifications, implying that sex chromosome inactivation does not occur in the Drosophila male germline. Instead, a lack of dosage compensation in spermatocytes accounts for the reduced expression from this chromosome. Finally, profiling uncovers dramatic ubiquitinylation of histone H2A and lysine16 acetylation of histone H4 across the Y chromosome in spermatocytes that may contribute to the activation of this heterochromatic chromosome.