Metaphase kinetochore movements are regulated by kinesin-8 motors and microtubule dynamic instability.

During metaphase, sister chromatids are connected to microtubules extending from the opposite spindle poles via kinetochores to protein complexes on the chromosome. Kinetochores congress to the equatorial plane of the spindle and oscillate around it, with kinesin-8 motors restricting these movements. Yet, the physical mechanism underlying kinetochore movements is unclear. We show that kinetochore movements in the fission yeast Schizosaccharomyces pombe are regulated by kinesin-8-promoted microtubule catastrophe, force-induced rescue, and microtubule dynamic instability. A candidate screen showed that among the selected motors only kinesin-8 motors Klp5/Klp6 are required for kinetochore centering. Kinesin-8 accumulates at the end of microtubules, where it promotes catastrophe. Laser ablation of the spindle resulted in kinetochore movement toward the intact spindle pole in wild-type and klp5Δ cells, suggesting that kinetochore movement is driven by pulling forces. Our theoretical model with Langevin description of microtubule dynamic instability shows that kinesin-8 motors are required for kinetochore centering, whereas sensitivity of rescue to force is necessary for the generation of oscillations. We found that irregular kinetochore movements occur for a broader range of parameters than regular oscillations. Thus, our work provides an explanation for how regulation of microtubule dynamic instability contributes to kinetochore congression and the accompanying movements around the spindle center.

INTRODUCTION

During mitosis, proper segregation of genetic material requires that sister chromatids of each chromosome attach to microtubules (MTs) extending from the opposite spindle poles via kinetochores, protein complexes on the chromosome (Cheeseman and Desai, 2008). When a pair of sister kinetochores establishes this configuration, termed biorientation, they oscillate around the equatorial plane of the spindle (Skibbens ). During oscillations kinetochores are predominantly under tension (Nicklas, 1988; Waters ; Dumont ). Kinetochore oscillations have been observed in various organisms, including yeast, amphibians, and mammals (Wise ; Skibbens ; Nabeshima ).

In the fission yeast Schizosaccharomyces pombe, the spindle in prometaphase consists of roughly 20 MTs extending from each spindle pole body (SPB, centrosome equivalent in yeast), and each sister kinetochore interacts with two to four MTs (Ding ). In this phase, sister kinetochores move around the spindle center until they get separated and pulled to the opposite SPBs in anaphase (Funabiki ). Kinetochore movements in fission yeast consist of a series of stochastic back-and-forth movements (Nabeshima ), which are less regular than in vertebrate cells (Skibbens ).

Kinetochore oscillations are regulated by kinesin-8 motors in fission yeast (Garcia ; West ; Mary ; Gergely ) and in human cells (Mayr ; Gardner ; Stumpff ; Hafner ). Depletion of kinesin-8 leads to oscillations with a larger amplitude, suggesting that kinesin-8 improves chromosome congression to the metaphase plate. Fission yeast has two kinesin-8 motors, Klp5 and Klp6, which form a heterodimer (West ). Deleting either motor protein or both impairs MT dynamics in the same manner (West , 2002; Garcia ; Unsworth ; Tischer ), yet there are also particular mitotic defects specific to individual deletions (Gergely ). In addition to kinesin-8, other motor proteins influence kinetochore oscillations. Inactivation of the kinesin-5 motor protein Kif11/Eg5 leads to irregular oscillations with a small amplitude (Vladimirou ). Chromokinesins (Levesque and Compton, 2001; Stumpff ; Wandke ), the XMAP215 family member TOGp (Cassimeris ), the centromere protein H (CENP-H) (Amaro ), and the MT depolymerizing kinesin mitotic centromere-associated kinesin (MCAK) (Jaqaman ) are also required for regular chromosome oscillations.

Kinetochore movements on the spindle require forces. Kinetochores are attached to plus ends of MTs and thus depend on the MT plus-end dynamics (Gachet ), which includes MT growth and shrinkage, as well as the switch from growth to shrinkage, termed catastrophe, and the opposite switch, termed rescue. Shrinking MTs can exert pulling forces on kinetochores (Grishchuk ), which stabilize kinetochore-microtubule attachments (Li and Nicklas, 1995; Akiyoshi ). When pulling forces are exerted on shrinking MTs, their shrinkage velocity decreases, oppositely from growing MTs, whose growth velocity increases (Akiyoshi ; Volkov ). Growing MTs exert pushing forces when growing against an obstacle, which results in a decrease of growth velocity (Dogterom and Yurke, 1997).

MT dynamics is regulated by MT-associated proteins. In particular, kinesin-8 motors have been shown to accumulate at the MT plus end in a MT length-dependent manner, where they facilitate MT catastrophe and shrinkage in vitro and in vivo (Gupta ; Varga , 2009; Mayr ; Grissom ; Tischer ; Erent ; Leduc ). The influence of these motors on MT dynamics has been explored theoretically (Hough ; Varga ; Tischer ; Reese ; Johann ; Leduc ; Melbinger ). Length-dependent accumulation of kinesin-8 on MT ends and regulation of catastrophe is responsible for the centering of the nucleus by MT pushing against the cell edge in interphase cells (Tolic-Norrelykke, 2010; Gluncic ) and may also be important for kinetochore centering on the spindle.

The mechanism of kinetochore oscillations and centering has been explored in several theoretical studies. An early model of oscillations in higher eukaryotic cells is based on force-dependent kinetics of MT attachment to and detachment from the kinetochore (Joglekar and Hunt, 2002). Studies of oscillations in budding yeast introduced force-dependent MT dynamics as a key MT property needed to generate oscillations (Sprague ; Gardner ), which was also used in models for higher eukaryotic cells (Civelekoglu-Scholey , 2013; Liu ; Armond ; Banigan ). Contrary to the assumption that tension at the MT-kinetochore interface promotes detachment of the MT (Joglekar and Hunt, 2002), another study assumed that the absence of tension stimulates detachment of the MT from the kinetochore, describing the activity of Aurora B kinase in fission yeast (Gay ). This model was extended by including a MT length-dependent pulling force (Mary ). Moreover, a recent model included the effect of kinesin-8 motors on MT dynamics, kinetochore attachment, and sliding force in the spindle (Gergely ). More complex mechanisms have been proposed for higher eukaryotic cells (Civelekoglu-Scholey , 2013; Liu ; Armond ; Banigan ). However, the physical mechanism underlying metaphase chromosome movements is still under debate.

RESULTS

Kinesin-8 has a major role in restricting kinetochore movements around the spindle center

To study kinetochore movements in fission yeast, we used cells in which the SPBs were visualized by Sad1-dsRed and centromeres of chromosome 2 (Cen2) were marked by green fluorescent protein (GFP) (Figure 1A; see Materials and Methods). SPB labeling allowed us to follow the position of spindle ends and spindle length, whereas Cen2 labeling showed the position of a single pair of sister centromeres and therefore of the respective sister kinetochores assembled on these centromeres. We will refer to the Cen2 spots as kinetochores. By imaging mitosis from the beginning of prophase until the onset of anaphase, we observed that as the spindle elongates, which is visible as the separation of the two SPBs, the two sister kinetochores move along the spindle in a stochastic manner during metaphase and separate toward the opposite poles in anaphase (Figure 1A and Supplemental Movie S1). Quantification of kinetochore movements in metaphase (Figure 1B) showed that in 94% of all time frames, the midpoint between sister kinetochores was located in the central 50% of the spindle length (n = 2132 time frames in 21 cells; Figure 1, C and D). Thus, kinetochore movements are centered on the spindle in wild-type cells.

FIGURE 1:

Kinesin-8 has a major role in kinetochore centering in S. pombe. (A–D) Kinetochore movements and positioning in wild-type cells. (A) Selected images (top) and kymograph (bottom) from a movie of mitosis in a wild-type cell (strain PD01) expressing Cen2-GFP (green) and Sad1-dsRed (magenta). Time is given in min:s; scale bar, 1 μm. See Supplemental Movie S1. (B) Distance from the spindle center of Cen2 (green) and SPBs (magenta) as a function of time for 21 wild-type cells (strain PD01). (C) Histogram of kinetochore position along the spindle (0 denotes spindle center, 1 denotes the SPB). The positions of the center between sister kinetochores in all time frames for spindles longer than 0.8 µm until the onset of anaphase were used; n is the number of cells. (D) Cumulative histogram of the data in C. (E–H) Kinetochore movements and positioning in klp5Δ cells (strain AK06, n = 20), legend as in A–D, respectively. See Supplemental Movie S2. (I) Histograms of kinetochore position along the spindle in the indicated genetic backgrounds (strains AK30-AK39), legend as in C. (J) SD of kinetochore position around the spindle center for the strains from panel I, calculated when the spindle length was normalized to the interval [–1,1].

Movie S1

Kinetochore oscillations in a wild-type cell (strain PD01) expressing Cen2-GFP (green) and Sad1-dsRed (magenta). Time is given in min:s; scale bar, 1 μm.

To explore the effect of kinesin-8 on kinetochore movements and positioning, we analyzed cells lacking Klp5 and found that kinetochores move in a stochastic manner as in wild type. However, contrary to wild type, kinetochores in klp5Δ cells are found along the whole spindle (Figure 1E and Supplemental Movie S2). Once moving, the kinetochores move to a larger extent than in wild type (Figure 1F). In only 67% of all time frames, the midpoint between sister kinetochores is found in the central 50% of the spindle length (n = 2080 time frames in 20 cells; Figure 1, G and H). To test whether the less efficient centering in klp5Δ cells is due to longer spindles in comparison with wild type, we analyzed in both strains spindles smaller than 1.47 µm (mean plus SD of spindle lengths in wild type). The SD of kinetochore position around the spindle center was 0.50 and 0.35 for klp5Δ and wild-type cells, respectively, calculated when the spindle length was normalized to the interval [-1,1]. These SDs are similar to those measured on spindles of all lengths, 0.53 and 0.34 for klp5Δ and wild-type cells, respectively, suggesting that the difference in centering efficiency is not due to differences in spindle length. Thus, Klp5 restricts kinetochore movements around the spindle center. This finding is consistent with previous studies in fission yeast (Mary ; Gergely ).

Movie S2

Kinetochore oscillations in a klp5Δ cell (strain AK06) expressing Cen2-GFP (green) and Sad1-dsRed (magenta). Time is given in min:s; scale bar, 1 μm.

Next, we explored the role of other MT motor proteins in kinetochore movements and positioning during metaphase. We used a candidate approach in which we deleted individually all genes in the S. pombe genome that had been identified to have homology to the kinesin motor domain (Schoch ) and dynein heavy chain dhc1 (Courtheoux ). We also deleted the nonmotor cross-linker ase1 (Loiodice ; Yamashita ). Note that the deletion of cut7, a kinesin-5 family member, was not studied because this deletion is lethal and inactivation of Cut7p in a temperature-sensitive mutant leads to spindle collapse (Hagan and Yanagida, 1992; Syrovatkina ). In each mutant, as well as in the respective wild-type strain, Cen2 was labeled with tdTomato and tubulin with GFP (Supplemental Table S1). We found that in all the mutants except the kinesin-8 klp6Δ the distribution of kinetochore position (Figure 1, I and J) and the general appearance of kinetochore movements (Supplemental Figure S1) were as in wild type. In the klp6Δ mutant, the distribution of kinetochore position was almost uniform along the spindle (Figure 1I), as in the klp5Δ mutant described above (Figure 1G). This result suggests that deletions of these two motor proteins have a similar effect on kinetochore movements and positioning, in agreement with a previous study (Mary ), and thus we focused on Klp5 in other experiments. We conclude that in S. pombe kinesin-8 motors Klp5/Klp6 have a major role in limiting kinetochore movements to the central region of the spindle, whereas other kinesins, dynein and Ase1, do not influence these movements significantly.

Polar microtubules undergo catastrophe less often and shrink faster in klp5Δ cells than in wild type

We next aimed to determine the effect of Klp5 on MT dynamics during mitosis. We focused on polar MTs, i.e., intranuclear MTs that grow from the SPB at an oblique angle with respect to the spindle during prometaphase (Zimmerman ; Vogel ; Kalinina ). The advantage of polar MTs for the measurement of MT dynamics is that their plus end is located away from the spindle and thus can be followed without interference of the signal from interpolar and kinetochore MTs. Moreover, polar MTs are free of forces generated at kinetochores or in overlap zones, which may affect MT dynamics, contrary to the MTs extending between the SPBs. Note that the structure that appears as a polar MT may be a bundle of a few MTs; thus our measurements correspond to the longest MT in the bundle.

We imaged polar MTs in wild-type and klp5Δ cells expressing mCherry-atb2 to visualize MTs (Figure 2, A and B) and analyzed their length in time (Figure 2C). We found that the catastrophe frequency is lower in klp5Δ cells than in wild type (Figure 2D). In klp5Δ cells the MTs grew to a longer length before undergoing catastrophe than in wild type (Figure 2D). Sixty-seven percent of wild-type polar MTs underwent catastrophe at a length that corresponds to the central third of the spindle (average spindle length 2.9 ± 0.9 µm, central third: 1.0−1.9 µm, n = 36 cells). In klp5Δ cells, only 45% of polar MTs underwent catastrophe at a length corresponding to the central third of klp5Δ spindles (average spindle length 4.1 ± 1.1 µm, central third: 1.3−2.7 µm, n = 33 cells). Moreover, the shrinkage speed of polar MTs in klp5Δ cells was higher than in wild type, whereas the growth speeds were similar (Figure 2D). In summary, these data indicate that Klp5 promotes catastrophe and slows shortening of polar MTs.

FIGURE 2:

Klp5 promotes catastrophe of polar MTs and accumulates on the plus ends in a MT length-dependent manner. (A, B) Time-lapse images of the mitotic spindle with polar MTs in a wild-type (A) and klp5Δ (B) cell. Cells express mCherry-atb2 (strains AK19 and AK20). Scale bars, 1 μm. (C) Length of polar MTs (distance between the plus end and SPB) in wild-type (black) and klp5Δ (magenta) cells as a function of time. The data are aligned so that the distance and time at MT catastrophe are 0. (D) Summary of polar MT dynamics in wild-type and klp5Δ cells; p = 4 × 10–5 (****), p = 7 × 10–4 (***). (E) Time-lapse images (sum projections) of a mitotic spindle in a cell expressing Klp5-GFP (white) and Sad1-dsRed (magenta, strain AK40). The image of SPBs (magenta) was recorded shortly before recording the time series of Klp5-GFP. The images of Klp5-GFP were bleach-corrected (see Materials and Methods). Scale bar, 1 µm. See Supplemental Movie S3. (F) Deconvolved time-lapse wide-field fluorescence images (maximum projections) of a mitotic spindle in a cell expressing Klp5-GFP (green) and Atb2-mCherry (magenta, strain AK15). Merged images (top), Klp5-GFP (middle) and Atb2-mCherry (bottom) are shown. Scale bar, 1µm, time interval 2.8 s. (G) Klp5-GFP intensity (bleach corrected) at the tip of a polar MT and (H) length of polar MTs as a function of time, for 21 polar MTs. Each MT is shown in a color corresponding to its average length, from blue for short MTs to red for long MTs. Time 0 is defined as the time of the maximal Klp5-GFP intensity. The black curve shows the average values and the gray area denotes SEM. (I) Klp5-GFP intensity (bleach corrected) at the tip of a polar MT as a function of polar MT length, for the same polar MTs shown in the same colors as in panels G and H. Data for individual MTs are plotted from start of the measurements until the peak of Klp5-GFP intensity. The black curve shows the average values of the binned data; error bars, SEM. Note that the outlier (the longest MT) was not included in the average. (J) Wide-field fluorescence image (sum projection) of S. pombe cells expressing Klp5-GFP (strain AK18, red arrows show Klp5-GFP on polar MTs) together with S. cerevisiae cells expressing Cse4-GFP (strain KBY7006, yellow arrows show Cse4-GFP at spindle poles in anaphase). Scale bar, 1 µm. (K) Histograms of Cse4-GFP signal intensities at anaphase poles (yellow) and Klp5-GFP at polar MT tips (red). Mean ± SEM for Cse4-GFP was 3.4 ± 1.0 a.u. (n = 66), and 1.4 ± 0.8 a.u. (n = 73) for Klp5-GFP.

Movie S3

Movie of a cell expressing Klp5-GFP (strain AK40). Note Klp5-GFP at the tip of polar microtubules. Time is given in min:s; scale bar, 1 μm.

Klp5 accumulates on the plus end of polar microtubules in a microtubule length-dependent manner, reaching a number of roughly 30 molecules

To understand how kinesin-8 motors are distributed along MTs in mitosis, we used cells expressing Klp5-GFP. We found Klp5 in a spotted pattern along the spindle (Figure 2, E and F). Because it was not possible in our assay to relate these Klp5 spots with the length of the MT to which Klp5 is bound, we again turned to polar MTs. We observed Klp5-GFP accumulating on the growing plus end of polar MTs (Figure 2E). Klp5-GFP signal intensity reached a maximum before the polar MT started to shrink (examples are shown in Figure 2, E and F, and quantification in Figure 2, G and H). Klp5-GFP signal intensity increased with increasing length of the polar MTs (Figure 2I), in agreement with previous measurements (Mary ). The increase in intensity was much faster for individual MTs than the increase averaged over several MTs (compare the slopes of colored data sets representing individual polar MTs and the slope of the black line representing the average in Figure 2I).

To estimate the average number of Klp5-GFP molecules on MT tips, we compared the signal intensity of Klp5-GFP spots with a signal intensity for which the corresponding number of molecules is known, following the approach from previous work (Ananthanarayanan ). We imaged exponentially growing fission yeast cells expressing Klp5-GFP and budding yeast cells expressing GFP-tagged histone H3-like centromere protein Cse4 in the same field of view (Figure 2J). In each budding yeast cell there are four to five Cse4 molecules per kinetochore (Coffman ; Lawrimore ). During late anaphase, kinetochores are found in two equal clusters, each close to one spindle pole. Each cluster contains 16 kinetochores and thus 64–80 Cse4 molecules. We found that the average signal intensity of Klp5-GFP on the tip of a polar MT was 2.43 lower than the signal intensity of Cse4-GFP at the spindle pole in anaphase (average Cse4-GFP intensity: 3.4 ± 1.0 a.u., n = 66; average Klp5-GFP intensity: 1.4 ± 0.8 a.u., n = 73; Figure 2K). Hence, we estimate that, on average, 26–33 Klp5-GFP molecules are located on the tip of a polar MT. Thus, our analysis of Klp5-GFP shows that Klp5 accumulates at the plus end of polar MTs in a MT length-dependent manner, typically reaching a number of 30 molecules.

Chromosome movements along the spindle are caused by pulling forces

To examine the forces acting on the kinetochores during metaphase in wild-type and klp5Δ cells, we performed laser ablation experiments of the spindle. We reasoned that if we destroy the MTs on one side of the sister kinetochores, their movement after ablation would provide information about the forces that act on the kinetochore. The metaphase spindle in wild-type cells is only 2–3 µm long, leaving not much space for ablation and for the subsequent movement of kinetochores. To increase the spindle size, we treated the cells with a low dose of hydroxyurea (Kim and Huberman, 2001). This treatment increased the spindle from 2.9 ± 0.9 µm to 3.6 ± 0.7 µm (n = 36 and 13 untreated and hydroxyurea-treated cells, respectively). This treatment did not affect the growth and shrinkage speed, catastrophe rate, and length at the time of catastrophe of polar MTs (Supplemental Figure S2), in agreement with previous results on interphase MTs (Tischer ).

We severed all spindle MTs on one side of the kinetochore pair, including interpolar and kinetochore MTs (Figure 3A). Spindle breakage and inward movement of the spindle poles were evidence for successful severing (Figure 3A), as previously described (Khodjakov ; Tolic-Norrelykke ; Raabe ; Maghelli and Tolic-Norrelykke, 2010, 2011). We found that in wild-type cells the kinetochores moved toward the SPB on the intact side of the spindle after ablation, irrespective of the kinetochore movement before ablation (n = 14 of 14 cells; Figure 3A, and Supplemental Figure S3A). This result is similar to previous observations that merotelic kinetochores move toward the intact side of the spindle after spindle ablation (Courtheoux ; Rumpf ; Cojoc ). We found that, after spindle ablation in klp5Δ cells, the kinetochores also moved toward the SPB on the intact side of the spindle (n = 10 of 10 cells; Figure 3B and Supplemental Figure S3A). These results suggest that kinetochores are moved by pulling forces in both wild-type and klp5Δ cells.

FIGURE 3:

Spindle cutting shows that chromosome movements are driven by pulling forces. (A) Laser cutting of the spindle in a wild-type cell. Top, time-lapse images of a hydroxyurea-treated wild-type cell expressing Cen2-GFP (green) and mCherry-atb2 (magenta, strain AK20). Laser ablation is marked by the yellow lightning sign; time is given in seconds; scale bar, 1 μm. See Supplemental Movie S4. Below, kymographs of the area around the spindle (shown to the left of the kymograph) of the same cell are displayed (merged channels, mCherry-atb2, and Cen2-GFP from top to bottom). Yellow lines mark the time of laser ablation. Vertical and horizontal scale bars are 1 μm and 10 s, respectively. (B) Laser cutting of the spindle in a klp5Δ cell (strain AK19), legend as in A. See Supplemental Movie S5. (C) Kinetochore position along the spindle, with respect to its position at the time of ablation, before and after ablation in wild-type (n = 14, left) and klp5Δ (n = 10, right) cells. The data are aligned to the time of ablation (0 on abscissa) and position of the kinetochore one frame before ablation (0 on ordinate). Laser ablation is marked by the yellow lightning sign. The plot shows the closer sister kinetochore in cases where two kinetochores were detectable. Thick lines represent the mean. (D) Poleward movement of kinetochores toward the intact side of the spindle in wild-type (n = 14, left) and klp5Δ (n = 10, right) cells. The data are aligned to the time and kinetochore position just before the poleward movement starts (0 on abscissa and ordinate, respectively). For wild-type cells, a linear fit yielded a poleward speed of 2.2 ± 0.5 µm/min, and for klp5Δ cells 2.5 ± 0.7 µm/min (mean ± SD, p = 0.14). (E) Top, time-lapse images of a hydroxyurea-treated wild-type cell expressing Cen2-GFP (green) and mCherry-atb2 (magenta, strain AK20). Laser ablation is marked by the yellow lightning sign; time is given in seconds; scale bar, 1 μm, 0 is time of ablation. Kymographs of the same cell are displayed below (merged channels, mCherry-atb2, and Cen2-GFP from top to bottom); vertical and horizontal scale bars are 1 μm and 10 s, respectively.

Movie S4

Spindle ablation in a hydroxyurea-treated wild-type cell (strain AK20) expressing Cen2-GFP (green) and mCherry-atb2 (magenta). Time is given in min:s; scale bar, 1 μm.

Movie S5

Spindle ablation in a hydroxyurea-treated klp5Δ cell (strain AK19) expressing Cen2-GFP (green) and mCherry-atb2 (magenta). Time is given in min:s; scale bar, 1 μm.

We then asked whether Klp5 motors affect the dynamics by which the kinetochores are pulled back toward the SPB after ablation. Poleward movement of the kinetochores started either immediately after ablation or after a short pausing time (Figure 3C). We did not detect a significant difference in the pausing time between wild-type and klp5Δ cells (13 ± 8 s and 16 ± 12 s, respectively; p = 0.49). Furthermore, poleward velocity of the kinetochores after ablation was similar in both cell types (–2.2 ± 0.5 µm/min for wild-type, –2.5 ± 0.7 µm/min for klp5Δ, p = 0.14; Figure 3D), suggesting a similar shrinkage speed of kinetochore MTs in the presence or absence of Klp5.

Finally, we tested whether pulling forces that move the kinetochores poleward after ablation are generated at the depolymerizing MT tip or along the MT. After ablation, we typically observed that as the kinetochore moved toward the SPB on the intact side of the spindle, the MT signal intensity between the kinetochore and that SPB was high, whereas the signal between the kinetochore and the ablation site was low (Figure 3E; additional examples are shown in Supplemental Figure S3B). We interpret the high MT signal as that of kinetochore MTs attached to the Cen2-labeled kinetochores together with other spindle MTs, whereas the low MT signal corresponds to the latter group only. Thus, the transition from high to low MT signal intensity visualizes kinetochore MT depolymerization in time, indicating that depolymerizing MT tips pull kinetochores toward the SPB.

Theoretical model for kinetochore congression and movements

To explore the physical basis of chromosome congression and movements around the metaphase plate, we develop a one-dimensional theoretical model, which is depicted in Figure 4A. The main idea of our approach is to use Langevin description of the fluctuations in the number of growing and shrinking MTs attached to kinetochores. We choose Langevin approach because it is suitable for systems with fluctuations in the number of particles (Van Kampen, 1992) and allows us to explore the ability of the system to generate movements around the center. In our model, which is formulated in Figure 4B and Materials and Methods, movements of kinetochores are calculated by the force-balance equations. The numbers of growing and shrinking MTs dynamically change, where tension increases MT rescue and decreases MT catastrophe (Akiyoshi ). Catastrophe is increased by kinesin-8 motors at the plus end of the MT (Tischer ), and the dynamics of motors along the MT length is described by the mean-field approximation, as for the motors on interphase MTs (Tischer ; Gluncic ).

FIGURE 4:

Theoretical model for metaphase kinetochore movements. (A) Scheme of the model. Kinetochores (green circles) connected by elastic spring (spring motive) are positioned at xR,L. MTs (magenta lines) of lengths LR,L, extend from the spindle poles (gray circles) to the kinetochores. Molecular kinesin-8 motors (black pictograms) attach at the MTs from nucleoplasm by constant rate kon and walk toward MTs plus ends with velocity υm. (B) Model equations. (C) Parameters of the model, see Materials and Methods, Choice of parameter values.

Our model has 15 parameters, 10 of which are taken from the literature or measured here (Figure 4C). To estimate the unknown values of the parameters that describe MT rescue, kr0 and fr, and to explore the behavior of the model, we use the model without kinesin-8 motors, c = 0. Numerical simulations show that a kinetochore pair moves along the spindle while the direction of the movement changes stochastically (Figure 5A). The values of the unknown parameters were obtained by comparison with experimentally measured distribution of kinetochore position (Figures 5B and 1F and Materials and Methods). Our model predicts that without kinesin-8 motors, kinetochores spend a similar amount of time under compression (~60%) and under tension (~40%) (Figure 5C).

FIGURE 5:

Results of the theoretical model for metaphase kinetochore movements. (A) Positions of the right (green) and left (blue) kinetochore shown as a function of time for the parameters given in Figure 4C for the model without the kinesin-8 motors, c = 0. (B) Distribution of the relative distance from spindle center of the central point between kinetochores. The SD of kinetochore position around the spindle center was 0.60, calculated when the spindle length was normalized to the interval [–1,1]. (C) Force at the kinetochore calculated for the model without kinesin-8 motors, c = 0. Negative force values correspond to the compression (light blue), while positive forces correspond to tension (red). (D) Position of the right (green) and left (blue) kinetochore shown as a function of time for the parameters given in Figure 4C. (E) Distribution of the relative distance from spindle center of the central point between kinetochores. The SD of kinetochore position around the spindle center was 0.37, calculated when the spindle length was normalized to the interval [–1,1]. (F) Force at the kinetochore calculated for the model with kinesin-8 motors. Color code as in C.

To determine the parameters that describe kinesin-8 motors, konckc1 and kd1, we used a model with kinesin-8 motors. The parameter values were obtained by reproducing experimentally measured lengths of polar MTs and the number of kinesin-8 motors at their plus ends (Supplemental Figure S4, B–E, and Materials and Methods). Here, to describe polar MTs, we modified the model by setting the force at the MT plus end to zero. For the parameters in Figure 4C, our model provides a prediction for the movement and positioning of the kinetochores in cells with kinesin-8 motors. Numerical solutions of the model show that a kinetochore pair moves along the spindle with stochastic changes in the direction of movement (Figure 5D). However, in comparison with the model without motors (Figure 5A), the movement is more restricted to the central region of the spindle. This centering effect of motors is reflected in the reduced SDs of kinetochore position as compared with the model without motors (compare Figure 5, B and E). The centering effect is due to accumulation of motors at the MT plus end in a MT length-dependent manner (Supplemental Figure S4D), which results in more catastrophes of longer MTs and thus less frequent excursions away from the spindle center. Kinetochores are predominantly under tension (98% of time, Figure 5F), in contrast to the model without kinesin-8 motors (see Figure 5C).

We compare the prediction from the model with experiments. The movements of kinetochores in experiments on wild-type cells are similar to those from the model (compare kymograph in Figure 1A with Figure 5D). The model also reproduces the histogram of kinetochore positions from experiments on wild-type cells (compare Figures 1C and 5E), providing a quantitative confirmation of the model prediction.

Next we explore how the model generates stochastic movements of kinetochores around the spindle center. To study the importance of the stochastic terms describing fluctuations in the number of growing and shrinking MTs, we investigate the movement of kinetochores in the case with and without stochastic terms (Figure 6). In the case without stochastic terms, our system attains stable oscillations around the spindle center or a fixed position at the center for different values of the motor concentration parameter (Figure 6A).

FIGURE 6:

Influence of the Langevin terms on dynamic behavior of the model and parameter space. (A) Solutions of the model without the stochastic terms for different values of konckc1 and initial conditions. Top, positions of the right (green) and left (blue) kinetochore shown as s function of time. Bottom, phase plot of the number of kinesin-8 motors at the left MTs plus ends and the position of the left kinetochore (gray line), while stable limit cycle is marked by blue line and stable fixed point by blue circle. Parameters and initial conditions are (I) konckc1 = 50 µm-1min-2, Q = 0; (II, left) konckc1 = 250 µm-1min-2, Q = 0; (II, right) konckc1 = 250 µm-1min-2, Q = 400; (III) konckc1 = 400 µm-1min-2, Q = 0. (I)–(III) Other parameters are as in Figure 4C. Initial conditions are xR, L = ± 0.1µm, , , , and . (B) Bifurcation diagram for the model without the stochastic terms. The amplitude of right kinetochore movement is plotted for different values of konckc1. Hysteresis occurs in interval konckc1 ≈ 60 - 275 µm-1min-2. Stable pole-to-pole kinetochore oscillations (solid line) and the unstable fixed position (dashed line) occur in region I. Stable oscillations and the stable fixed position coexist in region II. Amplitude of oscillations decreases from pole-to-pole (white) to a smaller values (light green) in region II. Stable fixed position occurs in region III. Other parameters are as in Figure 4C. (C, D) Parameter space for the model without and with the stochastic terms, respectively. Colors denote values of the SD of the central point between kinetochores calculated when the spindle length was normalized to the interval [–1,1], as a function of konckc1 and fr. For each point, equations are integrated for t = 1000 min. Grid size is 100 × 100. The same color bar applies to C and D. The black circle marks the position in the parameter space corresponding to Figure 4C. In C, arrowhead corresponds to the values of parameters investigated in B. In the dark blue area for fr < 2 pN kinetochores oscillate in antiphase. Stochastic terms are described under Materials and Methods, section Theoretical model. Other parameters are as in Figure 4C.

For small motor concentrations, kinetochores oscillate from one pole to the other (Figure 6A (I)). At intermediate concentrations, oscillations have a smaller amplitude and coexist together with a stable fixed point at the spindle center (Figure 6A (II)). For large concentrations, only a stable fixed point exists (Figure 6A (III)). Thus, for increasing concentrations an unstable fixed point becomes stable and our system undergoes a subcritical Hopf bifurcation (Figure 6B and Materials and Methods). The concentration range in which kinetochores oscillate without covering the whole spindle, similarly to kinetochore movement in wild-type cells, is narrow (Figure 6B, green). We conclude that our model with vanishing noise and biologically relevant parameters operates close to a subcritical Hopf bifurcation.

To explore the parameter space, we choose to vary the parameter related to motor concentration and the sensitivity of rescue to force. We find that the largest regions in the parameter space correspond to pole-to-pole oscillations (Figure 6C, red) and stable position in the spindle center (Figure 6C, blue). The range of motor concentrations for which the oscillations occur is broader for smaller values of the rescue critical force. Thus, the sensitivity of rescue to force promotes oscillations.

To test the influence of the stochastic terms on kinetochore movement, we explore the parameter space of the same parameters, without and with noise (Figure 6, C and D). Similarly to Figure 6B, in the case without noise the region in the parameter space corresponding to biologically relevant movements is narrow (Figure 6C, between red and blue regions). Interestingly, addition of the stochastic terms expands this region significantly (Figure 6D, orange to light blue). Thus, stochastic terms in our model ensure biologically relevant movements of kinetochores around the spindle center even for a large parameter variability.

DISCUSSION

Our model together with experiments suggests that kinesin-8 motors are required for kinetochore centering at the metaphase plate, whereas sensitivity of rescue to force is necessary for the generation of oscillations. Kinetochore centering is a consequence of accumulation of motors at the MT plus end in a MT length-dependent manner (Varga ), resulting in more catastrophes of longer MTs (Tischer , 2010) and thus less frequent movements of kinetochore away from the spindle center.

In our model we use Langevin approach to describe stochastic nature of MT dynamic instability, which is the main difference in comparison to previous simulations of metaphase kinetochore movements in fission yeast (Gay ; Mary ; Gergely ) and in human cells (Civelekoglu-Scholey , 2013; Liu ; Armond ; Banigan ). To describe stochasticity in MT dynamics, these previous models mainly used Monte Carlo simulations. By using Langevin approach, we found that kinetochore movement in the vicinity of metaphase plate occurs for a broad range of parameters only if MT dynamics is described in a stochastic manner. Similarly, the stochastic nature of the system ensures robustness of the dynamical behavior in rock-paper-scissors games (Reichenbach ). In our system the observed robustness of the kinetochore movements suggests that our model provides a plausible explanation for how a biological system with large variability generates typical kinetochore movements around the metaphase plate.

Our model can be extended to describe kinetochore congression to the metaphase plate and their oscillations in mammalian cells. In that case kinetochore fibers consist of roughly 20 MTs (McDonald ; Wendell ), for which we expect smaller influence of noise and therefore more regular oscillations. In addition to kinesin-8, which has been shown to enhance kinetochore centering, it will be necessary to include the bridging fiber in the model, because this fiber is laterally connected to the kinetochore fiber by passive cross-linkers (Vukusic ) and balances the tension between kinetochores (Kajtez ). A generalization of our model will reveal how these mechanisms contribute to kinetochore congression and the accompanying oscillations.

MATERIALS AND METHODS

Theoretical model

In our one-dimensional model we describe two oppositely oriented bundles of MTs connected to a pair of sister kinetochores (Figure 4A). The positions of the right and left kinetochore with respect to the cell center are denoted xR and xL, respectively. Subscripts R and L denote the right and left sides, respectively.

The movements of kinetochores at time t are calculated from balances of forces for each of them,

Elastic force Fel = k∆x, is generated by the cohesin protein that connects sister kinetochores and is described as Hookean spring of stiffness k and extension ∆x = (xR - xL) - xrest, where xrest denotes kinetochores rest length (Chacon ). Total forces are generated by MT bundles in which contributions of growing and shrinking MTs are taken into account, . Here and denote number of growing and shrinking MTs, respectively. Total number of MTs attached to each kinetochores is same, . Forces by which MTs act on the kinetochores are generated by MTs polymerization, , and MTs depolymerization, . These forces occur at MTs plus ends and are described by a linear force–velocity relationship,

Here fg denotes the force that stalls the MT growth and fs denotes the force that stalls the MT shrinkage. Also, υg denotes MT growth velocity, while υg denotes MT shrinkage velocity in the absence of force. This linear force–velocity relationship is an approximation of the exponential force–velocity relationship measured in Akiyoshi for both growth and shrinkage. Our linear force–velocity relationship approximates quantitatively the regime with forces smaller than the stall force, which is typical for our system: for example, when the left kinetochore is moving to the right, dx/dt > 0, the number of growing MTs is much larger than the number of shrinking MTs, , and the forces exerted by the growing MTs are much smaller than the stall force, (Supplemental Figure S4A). An analogous argumentation holds for the case when the left kinetochore is moving to the left. MTs in our model are rigid and straight, and thus their lengths are calculated as L = D/2 ∓ x, i = R, L, where D is spindle length. Growing MTs that are longer than spindle length, L > D, undergo catastrophes, whereas shrinking MTs of negative lengths L < 0 undergo rescue.

Numbers of growing and shrinking MTs in bundles change dynamically,

The rate of switching from shrinkage to growth, called rescue rate, increases when the force at the kinetochore increases (Akiyoshi ),

Here kr0 denotes dependence of the rescue rate on the tension and fr is the characteristic force of rescue. The rate of switching from growth to shrinkage, called catastrophe rate, is given as a function of the force and the number of kinesin-8 motors, . Using a Taylor expansion of the function g(Q) to the first order, we approximate g(Q) = 1 + kc1Q. Including the dependence of the catastrophe rate on force at the kinetochore (Akiyoshi ), we calculate the catastrophe rate as

Parameter kc0 denotes the catastrophe rate without force and fc is characteristic force of the catastrophe. The second term describes increase of the catastrophe rate with the number of the kinesin-8 motors at the plus ends of the microtubules, Q and kc1 denotes to which extent 1 motor increases catastrophe rate. The Langevin terms in Eqs. 3a and 3b describe fluctuations in rescue and catastrophe events and are calculated as . Function ξ(t) denotes Gaussian white noise, obeying ξ(t) = 0 and ξ(t)ξ(t′) = δ(t - t′), where δ is the Dirac delta function. Note that in the case without the stochastic terms, , our model becomes deterministic. We used reflective boundary condition in case if number of MTs decreases below zero or increases above total number of MTs at each kinetochore, Ntot.

To describe the distributions of kinesin-8 motors, we use a mean-field description by considering the motor density along the MT, ρ(u, t). For a coordinate system with the origin at the kinetochore and the positive direction of the u coordinate in the direction of the plus end, in the low-density limit the motor density obeys

where the motors bind to and unbind from the MTs at the rates kon and koff, respectively (Tischer ). The concentration of the motors in the nucleoplasm is denoted c and the velocity of the motors is denoted υm. The Heaviside step function θ ensures that the motors bind only along the MT. In the case when the motor velocity is larger than the MT growth velocity (υm > υg), the solution of the Eq. 6 reads

We use a limit of small motor unbinding rate in which this solution simplifies to

The number of motors at the plus end of the MTs is described as

based on the experimentally observed accumulation of kinesin-8 motors at the plus end of growing MTs (Figure 2, E–G). As the intensity of the kinesin-8 motors decreases roughly three times during 10 s after reaching the maximum values, while MT length decreases 20% (Figure 2, G and H), we use an exponential dependence of the motor detachment rate on the shrinkage velocity, kd, = kd0 exp (-kd1dL/dt), as a simplest function that changes faster than its input. Here kd0 denotes rate of motor detachment from the plus end of MT of a constant length, and kd1 denotes dependence of the rate of the detachment on the change of the MT length. This term ensures fast detachment of the motors from the MT plus end once the MT starts to shrink. We introduce a relative number of motors at the plus end as q = Q/konc, for which Eq. 7 becomes independent of parameters kon and c. In this case, the second term in Eq. 5 reads kc1Q = konckc1q. In our model, the parameters appear as one product, konckc1, which we use as one parameter. Differential equations 1–7 were integrated numerically in C++ programming language with time step ∆t = 0.1 s for initial conditions given in Figure 4A.

Derivation of the Langevin terms.

The Langevin terms are calculated for N MTs that are independent of each other and that stochastically switch between the growing and the shrinking state. In this case, we can represent the transition between the two states as random events that happen at certain rates. We denote by kres the rate at which one MT switches from shrinking to growing and by kcat the reverse transition. We assume that those transitions occur at constant rates, corresponding to an adiabatic approximation. This means that there is a probability kresn∆t that one of the n MTs in the shrinking state switches from shrinking to growing in the time interval (t, t + ∆t). If there are N MTs in a bundle, then one of (N - n) MTs in the growing state can switch from growing to shrinking with a probability of kcat(N - n)∆t. By taking infinitesimally small time intervals in which only one event is possible (i.e., ∆t → dt, (∆t)2 → 0), we obtain the N + 1 differential equations describing the probability that there are n MTs in the shrinking state at time t, valid for n = 0, 1, ..., N,

Here we introduce global rates of switching between global states with n MTs in a shrinking state

Equation 8 is a master equation for the probability that the system is in the state with n shrinking MTs, p, where the first term is the loss of probability from state n to either n - 1 or n + 1, whereas the second and third terms represent the gain of probability from states with n - 1 or n + 1 MTs, respectively. We define the linear “step” operator E for any function f(n) that depends on an integer variable n, as and its powers as , for . Using operator E we rewrite master equation or define Langevin equation for the change of the number of MTs in the shrinking state as Eq. 10, where for , and are expected value and variance, respectively. The term corresponds to the Langevin terms in Eqs. 3a and 3b (Van Kampen, 1992)

Please note that the Eqs. 8 and 10 are obtained from fundamental principles and do not include additional parameters.

Choice of parameter values

The parameters of the model were taken from the literature, measured here, or determined as follows (Figure 4C). Values for k and xrest are taken from Chacon and are within the range of values used in Gay . Growth velocity in the absence of force, vg, is measured here (Figure 2D). Note that this value is smaller than the one measured in Kalinina and Blackwell , most likely because of a different number of MTs in the bundle and the different strains and temperatures used in those studies. Shrinkage velocity in the absence of force, vs, is taken from ablation experiments, because the measured velocity of kinetochores in those experiments corresponds to the shrinkage velocity of kinetochore MTs in the absence of force from the opposite side (Figure 3D). Note that this value is smaller than the shrinkage velocity of polar MTs measured here (Figure 2D) and in Sagolla , Kalinina , and Blackwell . We choose MT shrinkage stall force to be the same as MT growth stall force. This choice does not influence the model behavior because the system is in the regime of forces smaller than the stall force. The catastrophe rate in the absence of force and kinesin-8 motors, kc0, is taken from kinesin-8 deletion experiments (Figure 2D). Note that this value is smaller than the value from Blackwell . Other parameters that should be determined are kc0, fr, kd1 and the product konckc1. By comparing our experimental and theoretical results for the distributions of kinetochore positions, in the case of kinesin-8 deletion (Figures 1 and 5), c = 0, and by comparing the lengths of the polar MTs in the case with kinesin-8 motors (see Theory of polar MTs), we obtain the values given in Figure 4C. Our values of the parameters kr0 and fr are within the range of in vitro measurements, kr0 = 1.43 ± 0.47 min-1 and fr = 6.4 ± 4.2 pN in Akiyoshi . Note that our value of kr0 is smaller than the value from Mary and larger than the one from Blackwell .

Theory of polar MTs.

The parameter that describes kinesin-8–dependent catastrophe, konckc1, was set to reproduce the experimentally measured lengths of the polar MTs (Supplemental Figure S4C). In our theory, polar MTs were modeled by setting parameters Fel = 0, Ntot = 1, and initial conditions , LL = 0, and qL = 0. Note that we use Ntot = 1 because in our experiments we track most likely a single MT, namely the longest one in the bundle. We do not include forces generated by the interaction between the nuclear envelope and MTs to avoid extension of the model for the particular case of polar MTs. Equations were integrated without the stochastic terms. Using the experimental result that there are ∼30 kinesin-8 motors at the MT tip before entering catastrophe, we were able to fix the value of parameters kc1 = 2 min-1 and the product konc = 50 µm-1 min-1, giving the model parameter konckc1. The parameter that describes kinesin-8 motors detachment from the MTs, kd1, was fixed to reproduce the experimental dependence of the number of kinesin-8 motors on MT length (Supplemental Figure S4, D and E). We find that for kd1 = 2 min/µm (Supplemental Figure S4D) the agreement with experiments (Figure 2I) was better than for kd1 = 1 min/µm and kd1 = 5 min/µm (Supplemental Figure S4E).

Strains and media

Fission yeast strains used in this study are listed in Supplemental Table S1. Cells were grown on yeast extract (YE) or Edinburgh minimal medium (EMM) with appropriate supplements at 25 ± 0.5°C in a Heraeus incubator (Thermo Scientific) (Forsburg and Rhind, 2006).

To construct strain AK06, klp5 was deleted using the PCR gene-targeting method (Bahler ). The klp5 open reading frame (ORF) was hereby exchanged by homologous recombination with a NatMX-cassette (nourseothricin resistance). The primers were designed using the web tool www.bahlerlab.info/resources/ (Penkett ). Klp5-deletion-forward-primer: 5′GGTTATTGGA­ATTGTTGTTTGCGAGTGCCCTTCTTTTCATCTGTTTTGGATGATTGCTTGCATGGAAATTTTGCTTTTCACGGATCCCCGGGTTAATTAA3′

Klp5-deletion-reverse-primer: 5′ATCATCAAGCTTA­TCCGTTTTTTTTTTTTAAATATACCCAACAGGATATTTAGAGGATTCGTATTTGAATATACGGACCTGAATTCGAGCTCGTTTAAAC3′

The primers contain 80 base pairs homologous to the flanking sequences directly before and after the ORF of Klp5 and 20 base pairs homologous to template plasmid pFa6a-NatMX. DNA fragments including the homology domains of Klp5 upstream and downstream region and the NatMX cassette were PCR amplified. The product of the PCR was purified using the QIAquick PCR Purification Kit for microcentrifuges (QIAGEN, Hilden, Germany) and then concentrated by isopropanol precipitation. The strain PD01 was transformed with the 15 µg Klp5-deletion-cassette, diluted in 15 µl 10 mM Tris-Cl, pH 8.5, by electroporation at 1.5 kV using an Eppendorf Electroporator 2510 (Eppendorf) (Forsburg and Rhind, 2006). Klp5-deleted strains were selected by growth on YE + supplements + nourseothricin medium. Deletion was checked by PCR after DNA isolation (Forsburg and Rhind, 2006) using the following primers: 5′-CCATTCAAAACAAAGTGCGTTA-3′ (forward); 5′- GGA­CCCATCCAGTGCCTCGAT-3′ (reverse inside NatMX-­casette), as also DNA sequencing (sequencing primer: 5′-CCATTCAAAACAAAGTGCGTTA-3′).

Strains AK31-AK39 were generated by crossing of strains of the S. pombe Deletion Mutant Library (Bioneer) with strain PP682 (see Strains and Plasmid, Supplemental Table S1) following the Random Spore Analysis protocol (Forsburg and Rhind, 2006). For ectopic expression of GFP-atb2, cells were transformed with pREP81-GFP-atb2 using a chemical transformation method based on lithium acetate (Forsburg and Rhind, 2006). pREP81-GFP-atb2-transformed strains were selected with and kept on EMM lacking l-leucine, supplemented with 20 µM thiamine to suppress nmt81-promotor regulated GFP-atb2 expression.

Sample preparation

Cells were grown in YE liquid medium including the appropriate supplements overnight and the next day diluted 3–5 h before the measurements to achieve an exponentially growing culture. The cells were incubated on lectin-coated (L2380; Sigma-Aldrich) 35-mm glass-bottom culture dishes (MatTek Corporation) for 10 min for attachment and then washed with and kept in EMM supplemented with adenine, l-histidine, l-leucine, uracil, and l-lysine for live-cell imaging. All strains carrying pREP81-GFP-atb2 were grown in EMM lacking l-leucine supplemented with 1 µM thiamine during preparation and imaging. For hydroxyurea-treatment, hydroxyurea (H8627; Sigma-Aldrich) was added at a final concentration of 3 mM during preparation and imaging, according to the conditions used in a previous study (Tischer ).

Microscopy and laser ablation

Spinning disk confocal microscopy setup 1.

Live-cell images were taken using an Andor Revolution Spinning Disk System (Andor Technology), consisting of a Yokogawa CSU-X1 spinning-disk scan head (Yokogawa Electric Corporation), connected to an Olympus IX81 inverted microscope (Olympus). The microscope was equipped with a Prior ProScanIII xy- scanning stage (Prior Scientific), a Prior NanoScanZ (Prior Scientific) and an Olympus UPlanSApo x100/1.35 NA oil objective (Olympus). The microscope was equipped with an iXon EM+ DU-897 BV back-illuminated electron-multiplying charge-coupled device (CCD) camera, cooled to –80°C (Andor Technology). The resulting xy-pixel size in the images was 81 nm. The system was controlled by Andor iQ2 software version 2.9 (Andor Technology). Excitation for imaging was done using a sapphire 488-nm solid-state laser (50 mW; Coherent) and a 561-nm solid-state laser (50 mW; Cobolt Jive). The laser power was controlled using the acousto-optic tunable filter in the Andor Revolution Laser Combiner (Andor Technology).

For recording chromosome oscillations in strains lacking kinesins, dhc1, or ase1, image stacks of 11 z-planes of 500-nm distance were taken at room temperature (22–24°C) with a time interval of 10 s, exposing sequentially for 100 ms with 15% of the 488-nm laser and 250 ms with 15% of the 560-nm laser within the time interval.

For recording Klp5-GFP on polar MTs, image stacks of 7 z-planes of 500-nm distance were taken at room temperature (22–24°C) with a time-interval of 1 s, exposing each plane for 100 ms with 50% of the 488-nm laser within the time interval. The Sad1-dsRed signal was recorded before and after recording the time-lapse series of Klp5-GFP, exposing for 200 ms with 27% of the 560-nm laser. A razor edge 488LP emission filter (Semrock) was used for the detection of the Klp5-GFP signal and a razor edge 568LP emission filter (Semrock) for the Sad1-dsRed signal.

Spinning disk confocal microscopy setup 2 with a laser ablation module.

Laser ablation experiments and imaging of basic chromosome oscillations of cen2-GFP/Sid4-dsRed-labeled strains were performed on a Yokogawa CSU10 spinning disk system (Yokogawa Electric Corporation) connected to an Olympus IX71 inverted microscope (Olympus) with an iXon EM+ DU-897 BV back-illuminated electron-multiplying CCD (Andor Technology), cooled to –80°C, electron multiplication gain 300. For excitation, a sapphire 488-nm solid-state laser (75 mW; Coherent) and a Jive 561-nm solid-state laser (75 mW; Cobolt) were used for GFP and tdTomato/dsRed, respectively. The laser intensity was controlled using the acousto-optic tunable filter inside the Andor Revolution Laser Combiner (Andor Technology). For sequential imaging, the emission wavelength was selected using respective emission filters BL 525/30 (Semrock) and ET 605/70 (Chroma) mounted in a fast, motorized filter wheel (Lambda-10B; Sutter Instrument Company). All images were recorded with an Olympus UPlanSApo × 100/1.4 NA oil objective (Olympus). The resulting xy-pixel size in the images was 174 nm; the z-distance between optical sections was 500 nm. The system was controlled by Andor iQ software version 2.9 (Andor Technology). The temperature was kept at 25°C using a Warner Heating Chamber (Warner Instruments).

For recording the oscillations of cen2-lacOp/lacI-GFP and Sid4-dsRed signals in wild-type and klp5Δ cells, image stacks of 15 z-planes of 500-nm distance were taken with a time interval of 10 s, exposing sequentially for 140 ms with 10% of the 488-nm laser and 210 ms with 10% of the 561-nm laser within the time interval.

For recording the oscillations of cen2-lacOp/lacI-GFP and pRep81-GFP-atb2 in limited candidate screen, image stacks of 11 z-planes of 500 nm distance were taken with a time interval of 10 s, exposing sequentially for 100 ms with 15% of the 488-nm laser and 200 ms with 18% of the 561-nm laser within the time interval.

For recording dynamics of atb2-mCherry-labeled polar MT of strains AK19 and AK20, image stacks of 7 z-planes of 500-nm distance were taken with a time interval of 2 s, exposing for 250 ms with 14.5% of the 561-nm laser. Cells were kept fixed to the glass slide using a 4% agarose in EMM + all amino acid agarose pad. AK19 and AK20 express the cen2-lacOp/lacI-GFP signal, which was checked in each cell before measurement to confirm metaphase.

For recording chromosome movements before and after ablation, image stacks of 5 z-planes of 500-nm distance were taken with a time interval of 2 s, exposing sequentially for 100 ms with 22% of the 488-nm laser and 100 ms with 20% of the 561-nm laser within the time-interval. Spindle severing was performed using a MicroPoint (Andor Technology) with 408-nm dye resonator cell.

Imaging of cells expressing Klp5-GFP together with S. cerevisiae cells expressing Cse4-GFP.

Budding yeast cells (strain KBY7006) were grown on yeast peptone dextrose agar with kanamycin at 25 ± 0.5°C. Exponentially growing cells of strains AK18 and KBY7006 were mixed and imaged in a poly-l-lysine–coated glass bottom dish using a DeltaVision Microscope, equipped with an Olympus 100 × /1.40, UPLS Apo oil objective, using a GFP-filter cube, a 50% neutral density filter and an exposure time of 300 ms and 4x camera gain. Pixel size was 64 nm. Z-stacks of seven planes with 500-nm distance were recorded.

Widefield imaging of Klp5-GFP and Atb2-mCherry.

Exponentially growing cells of strain AK15 were imaged in poly-l-lysine–coated glass bottom dishes using a DeltaVision Microscope, equipped with a UPLS Apo × 100/1.40 NA oil objective. Klp5-GFP was imaged using a GFP-filter cube, Atb2-mCherry a mCherry filter cube. A 32% neutral density filter, an exposure time of 300 ms, and 4x camera gain was used for both channels.

Image processing and data analysis

Analysis of Klp5-GFP on polar MTs.

For evaluation of the Klp5-GFP signal, sum projections of seven z-planes were bleach corrected using the Bleach Correction tool (Kota Miura and Jens Rietdorf) in Fiji (Schindelin ) using the exponential fitting method on the entire field of view. The first nine time frames were excluded from the data analysis due to unproportionally strong bleaching, presumably due to autofluorescent organelles. The tip of the polar MT in the bleach-corrected sum projection was tracked using the Low Light Tracking tool in Fiji (Krull ), giving the center and the intensity of the distal Klp5-GFP spot. For measuring the length of the polar MT, the position of the SPB (Sad1-dsRed) originating the polar MT was determined within the sum projection using the Low Light Tracking tool (Krull ), and the two-dimensional distance between SPB and polar MT tip was calculated. Data were further analyzed and plotted using scripts written in MATLAB (The Mathworks).

SPB and kinetochore tracking.

Tracking of the centromeres/kinetochores of chromosome 2 (cen2-lacOp/lacI-GFP) and the spindle pole bodies (Sid4-GFP)/the ends of the spindle (pRep81-GFP-atb2) was performed on maximum projections of 15 z-planes (imaging conditions see above) using the Low Light Tracking tool in Fiji (Krull ). Sigma was set to 0.8 constant for centromere signals to 1.0 constant for SPB/spindle signals. Only cells that underwent complete metaphase and entered anaphase were analyzed. Data was further analyzed and plotted using scripts written in MATLAB (The Mathworks). The position of the kinetochore of chromosome 2 were determined relative to the center of the mitotic spindle, as a projection of the Cen2 position onto the spindle axis.

Histograms of kinetochore position.

We calculated the midpoint between the sister kinetochores and the relative position of this midpoint along the spindle, defining the spindle center as 0 and the position of the closer SPB as +1. Evaluation was performed for all time frames from the time point when the spindle had reached a length of 0.8 µm until the onset of anaphase. The onset of anaphase was defined as the time point when the two sister kinetochores segregate completely, moving toward the SPBs without approaching each other again.

Analysis of polar MT dynamics.

Z-stacks were maximum projected and the tip of the polar MT and the extreme ends of the spindle were manually tracked in Fiji using the Manual Tracking tool in Fiji (Fabrice Cordelières, Institut Curie, Orsay, France). Data were further analyzed and plotted using scripts written in MATLAB (The Mathworks). In a few cases (strain AK19: 5 of 52 events, strain AK20: 4 of 51 events), we measured a growth speed lower than 0.5 µm/min, which we considered as pausing, and we excluded those data from the analysis of growth and shrinkage dynamics. Catastrophe frequency was calculated considering the time spent from the start of growth until catastrophe.

Analysis of laser ablation experiments.

The position of the kinetochores of chromosome 2 and the intact spindle end, that is, the end of the atb2-mCherry signal that is connected to the kinetochore, was tracked using the Low Light Tracking tool in Fiji (Krull ). Data were further analyzed and plotted using scripts written in MATLAB (The Mathworks). The distance between intact spindle pole body and the closer sister kinetochore was calculated. The speed of the poleward movement was calculated from the first time point after pausing (if present) until the kinetochore had reached a distance of 420 nm to the spindle end. The pausing time was manually determined. Only cells in which the spindle was clearly cut into two fragments were evaluated.

Quantification of the Klp5-GFP signal intensity by comparing it with the Cse4-GFP signal intensity in S. cerevisiae.

Average pixel intensities of the sum projections were measured in a circular selection of 11 pixels diameter (708 nm), centered on the tip of polar MTs (S. pombe) or at the anaphase pole (S. cerevisiae).

Statistical analysis

Data are presented as mean ± SEM. For comparison of data sets, t test was used, with a significance level of 0.05. Data were analyzed in MATLAB (The Mathworks). Figures were assembled in Adobe Illustrator (Adobe Systems).

Click here for additional data file.