Microtubule reorganization during female meiosis in C. elegans
Abstract
Most female meiotic spindles undergo striking morphological changes while transitioning from metaphase to anaphase. The ultrastructure of meiotic spindles, and how changes to this structure correlate with such dramatic spindle rearrangements remains largely unknown. To address this, we applied light microscopy, largescale electron tomography and mathematical modeling of female meiotic Caenorhabditis elegans spindles. Combining these approaches, we find that meiotic spindles are dynamic arrays of short microtubules that turn over within seconds. The results show that the metaphase to anaphase transition correlates with an increase in microtubule numbers and a decrease in their average length. Detailed analysis of the tomographic data revealed that the microtubule length changes significantly during the metaphasetoanaphase transition. This effect is most pronounced for microtubules located within 150 nm of the chromosome surface. To understand the mechanisms that drive this transition, we developed a mathematical model for the microtubule length distribution that considers microtubule growth, catastrophe, and severing. Using Bayesian inference to compare model predictions and data, we find that microtubule turnover is the major driver of the spindle reorganizations. Our data suggest that in metaphase only a minor fraction of microtubules, those closest to the chromosomes, are severed. The large majority of microtubules, which are not in close contact with chromosomes, do not undergo severing. Instead, their length distribution is fully explained by growth and catastrophe. This suggests that the most prominent drivers of spindle rearrangements are changes in nucleation and catastrophe rate. In addition, we provide evidence that microtubule severing is dependent on katanin.
Introduction
During meiosis, haploid gametes are produced from diploid progenitor cells in most sexually reproducing eukaryotes. The reduction of chromosome number in the gametes is essential for the development of the fertilized oocyte into an embryo. Therefore, the accuracy of meiosis is crucial as errors in chromosome segregation result in aneuploidies and lead to genetic disorders, which are often lethal (Hassold and Hunt, 2001).
In this paper, we study the structure and dynamics of female meiotic spindles in Caenorhabditis elegans. Female meiotic spindles, like mitotic and male meiotic spindles, are built from microtubules (MTs), motor proteins, and microtubuleassociated proteins. However, unlike other spindles in most animals, female meiotic spindles typically have no centrosomes due to centriole elimination during oogenesis. Thus, female meiotic spindles are assembled differently from those that rely on centrosomebased microtubuleorganizing centers. Specifically, this is true in humans (Holubcová et al., 2015), mice (Schuh and Ellenberg, 2007), and also in nematodes (Albertson and Thomson, 1993). It is an open question how their structure reflects the differences in assembly.
In C. elegans, female meiotic spindles undergo a drastic reorganization during the metaphasetoanaphase transition. During metaphase of meiosis I and II, microtubules form a pointed and elongated bipolar spindle, which then shortens during anaphase. In addition, the microtubule density shifts from the spindle poles to the midzone between the separating chromosomes. This midzone then extends until the end of anaphase (Dumont et al., 2010; Laband et al., 2017; Yu et al., 2019). The structural basis of these rearrangements and their role in meiotic chromosome segregation are not well understood. In particular, it remains unclear whether these structural rearrangements are driven by kataninmediated severing (Joly et al., 2016; Srayko et al., 2006), transport (Mullen and Wignall, 2017; Brugués et al., 2012), or changes in MT nucleation or polymerization dynamics (Brugués et al., 2012; Needleman et al., 2010).
Here, we study the mechanisms of structural microtubule arrangement by using largescale electron tomography (Redemann et al., 2017; Redemann et al., 2018; Fabig et al., 2020) to obtain the 3D ultrastructure of spindles in C. elegans oocyte meiosis. We find that female meiotic spindles are composed of arrays of short microtubules, 90% are shorter than half spindle length, which are highly dynamic and turnover within 10 s. During the transition from metaphase to anaphase, the number and lengths of microtubules change significantly. To understand the drivers of these changes, we developed a mathematical model for the length distributions of microtubules that includes lifetime (turnover), growth, and microtubule severing. We inferred the parameters of the model from data using Bayesian inference and Markov Chain Monte Carlo sampling. Our data suggest that the changes in microtubule number and lengths are mostly caused by changing microtubule turnover. While severing clearly occurs during metaphase, our data suggests that it affects only the small fraction of microtubules that contact the chromosomes. Here and in the following, we define that a microtubule contacts a chromosome if any point of its lattice is within 150 nm of a chromosomes surface. The rearrangements of spindle architecture towards anaphase are likely caused by changes in nucleation and not by any changes in cutting rates.
Results
Threedimensional reconstruction reveals structural changes of the spindle throughout female meiosis
Previous analyses of meiotic spindle microtubules using EM tomography were based on partial, rather than full, reconstructions of the spindle. This made it impossible to measure the total length of microtubules composing the spindle since most microtubules extend beyond the tomographic volume. We obtained the spindle ultrastructure of C. elegans oocytes by serial section electron tomography, which allows the reconstruction of whole spindles in 3D with singlemicrotubule resolution (Redemann et al., 2017). We reconstructed two meiotic spindles in metaphase I (T0243.5 and T0243.7), as well as four spindles at early (T0243.4), mid (T0208.1), and late stages of anaphase I (T0275.10 and T0369.8; Figure 1A, Figure 1—videos 1, 2; Figure 3—video 1) for meiosis I (Table 1). In addition, we reconstructed two metaphases (T0186.1 and T0209.7), one early anaphase (T0186.3), and one mid anaphase (T0234.6) stage of meiosis II oocytes (Figure 1B, Table 2). We had previously analyzed the modes of microtubule attachment to chromosomes (endon vs. lateral) and discussed the resulting role of microtubules in chromosome segregation during anaphase using the meiosis I metaphase (T0243.5), early anaphase (T0243.4), mid anaphase datasets (T0208.1), (Redemann et al., 2018) and late anaphase (T0275.10) (Yu et al., 2019) of the first meiotic division. Here, we used those datasets, together with newly generated additional datasets of meiosis I metaphase (T0243.7), and late anaphase (T0369.8), and two meiosis II metaphases (T0186.1, T0209.7), early anaphase (T0186.3), and mid anaphase (T0234.6), to quantitatively characterize the spindle structure and microtubule rearrangements during female meiosis.
For this, we compiled summary statistics that characterize spindle morphology. In Tables 1 and 2, we report the spindle length, the number of microtubules, the number of microtubules that are located within 150 nm from the chromosome surface, the average microtubule length, the maximum microtubule length, the total polymer length (added length of all microtubules), and the dimension of the tomogram and spindle dimensions for each stage of meiosis I (Table 1) and II (Table 2). We also indicated new and previously analyzed datasets in this table.
Female meiotic spindles are arrays of short microtubules
We first sought to characterize the overall spindle structure. Female meiotic spindles could be constructed mainly from long microtubules, reaching from near the spindle pole to chromosomes (as in C. elegans mitotic spindles; Redemann et al., 2017), or from short ones distributed throughout the spindle volume, as suggested by Srayko et al., 2006, and, for instance, meiotic Xenopus spindles (Brugués et al., 2012; Burbank et al., 2006; Yang et al., 2007).
Our data show that throughout meiosis I and meiosis II, roughly 40–50% of microtubules are shorter than 500 nm (Figure 1D). In contrast, the fraction of microtubules that are at least half spindle length or longer is only 10%. Throughout meiosis, the mean length of microtubules is approximately 20% of the spindle length (Tables 1 and 2).
Furthermore, we determined the position of the poleproximal ends of microtubules along the spindle axis. In metaphase, these putative microtubule minusends were found throughout the spindle (Figure 1E), with no marked preference for the poles. In anaphase, about 60–70% of microtubules are found within a distance of 150 nm from the chromosomes, and the others are distributed throughout the spindle (Figure 1F).
Together, these findings directly confirm that the acentrosomal female meiotic C. elegans spindle is an array of short microtubules, more similar to meiotic Xenopus spindles than to mitotic spindles in C. elegans, as previously inferred by others based on light microscopy (Brugués et al., 2012; Burbank et al., 2006; Yang et al., 2007). We conclude that the distance between poles in metaphase (5 µm) and chromosomes in anaphase (3–4 µm) is bridged by an array of short overlapping microtubules rather than long ones spanning the entire distance. This is an overall interesting finding as C. elegans meiotic spindles are rather short in comparison to Xenopus meiotic spindles and could easily generate microtubules long enough to bridge the entire spindle length, which is much harder to achieve in a 60 µm Xenopus spindle.
Meiotic spindles are composed of shortlived, fastmoving microtubules
We next sought to understand the role of microtubule movements and dynamics in the structural changes observed during the transition from metaphase to anaphase. Since electron tomography generates static snapshots, we investigated this question by light microscopy. Using fluorescence recovery after photobleaching (FRAP), we first measured microtubule turnover and motion in metaphase I spindles. By photobleaching a small stripe near the center as well as close to the poles of metaphase spindles (Figure 2A, arrows, Video 1), we observed a halftime of microtubule recovery in metaphase of 4.9 s ± 3.4 s (n = 6) in the spindle center and 4.1 s ± 2.7 s (n = 6) at the spindle poles (Figure 2B). In addition, we observed that the photobleached stripes close to the spindle poles showed a rapid poleward movement with a rate of 8.5 ± 2.2 µm/min (n = 6, Figure 2C).
We conclude that microtubules in C. elegans meiosis turn over rapidly and show substantial poleward motion. This is notably different from the first mitotic division in the early C. elegans embryos, where microtubule motion was not detected in similar experiments (Labbé et al., 2004; Redemann et al., 2017). Thus, in female meiosis, both nucleation/depolymerization and microtubule motion might be involved in shaping the spindle structure. Similar microtubule dynamics have been described in Xenopus meiotic spindles using fluorescent speckle microscopy (Burbank et al., 2006; Yang et al., 2007; Burbank et al., 2006; Yang et al., 2007).
Spindle rearrangements during meiosis correlate with substantial changes in microtubule number
To better understand the role of microtubule polymerization dynamics, we next asked whether the number of spindle microtubules changed during spindle rearrangements. From our reconstructions, we detected 3662 and 3812 microtubules in metaphase of meiosis I (3013 and 3808 in meiosis II), 7011 microtubules in early anaphase of meiosis I (5572 in meiosis II), 3317 in mid anaphase of meiosis I (2246 in meiosis II), and 1511 and 1306 microtubules in late anaphase of meiosis I (Tables 1 and 2). This shows that the number of microtubules almost doubles between metaphase and early anaphase. To get a more spatially resolved picture, we also measured the number of microtubules along the spindle axis in our datasets. This measurement showed that the number of microtubules increases by more than twofold in the center of the spindle from metaphase to early anaphase, while the number of microtubules at the spindle poles only slightly increases from metaphase to early anaphase and then decreases throughout anaphase (Figure 3A). During anaphase, the number of microtubules decreased on the poles and almost halved in the spindle center (Figure 3A, B). Thus, we conclude that changing the number of microtubules is correlated with the structural spindle rearrangements that occur during meiosis.
Microtubule rearrangements during meiosis are driven by changes in microtubule dynamics, not nucleation rates alone
We next sought to gain insight into the mechanisms that drive the observed changes in microtubule numbers. We first asked whether changes in nucleation rate alone could explain the data. If this were the case, microtubule growth, shrinkage, and severing rates would remain the same, and the average length of the microtubules would remain unchanged throughout meiosis. However, the total polymer amount, the length of all microtubules added together, would change in step with the microtubule number.
Interestingly, we found the opposite result: the total polymer length of microtubules remained almost constant from metaphase to anaphase even though the number of microtubules increases twofold (Tables 1 and 2). Furthermore, during the same time, the average length of microtubules went down by a factor of 2 (Figure 1D, Tables 1 and 2). Only in late anaphase, as the spindle continued to elongate, did the number of microtubules and the total polymer length start to decrease. Therefore, the observed changes in spindle morphology were not caused by changes in nucleation rates alone, but also by changes in microtubule disassembly. In principle, this could be caused by local or global changes in microtubule severing, changes in microtubule growth dynamics, or a combination thereof. In the following, we will try to disentangle these possibilities.
Microtubule dynamics during spindle rearrangement are globally regulated
We next investigated whether changes in microtubule length were uniform along the spindle axis by quantifying the average microtubule length along the spindle axis. This analysis revealed that, during the transition from metaphase to early anaphase, the average length of microtubules decreased everywhere in the spindle (Figure 3C, D). However, the magnitude of the decrease in length varies as a function of position. The average length of microtubules in the middle of the spindle decreased nearly twice as much as the average length of microtubules near the spindle poles (Figure 3C, D). The analysis of average microtubule length along the axis also showed that in all spindles, in metaphase and anaphase, the average length of microtubules is shorter near the spindle poles (Figure 3C, D), which in anaphase is located at the inner surface of the chromosomes.
We more closely mapped the position of microtubules according to their length (Figure 4A–G) by plotting the number of microtubules that are 500 nm and shorter as well as the number of microtubules between 0.5–1 µm and 1–1.5 µm along the spindle axis (Figure 4H–N). This analysis showed that short microtubules, 500 nm and below, are localized closer to the spindle poles, while longer microtubules, between 0.51 and 1.5 µm, are localized to the spindle center. We also found that the number of short microtubules increased from approximately 35% of all microtubules in metaphase to 50% of all microtubules in anaphase (Figure 4A–G).
Throughout anaphase, we noticed that there were more short microtubules in the spindle half closer to the cell cortex (Figures 3 and 4). This difference was most apparent in early and mid anaphase, where the average microtubule length in the cortical spindle half was 1.00 µm ± 0.6 (p=3.6 E6) and 1.17 µm ± 0.8 in the cytoplasmic half, and in early anaphase with an average length of 0.91 ± 0.53 (p=0.002) in the cortical half and 1.04 µm ± 0.6 in the cytoplasmic half. We do not know at this point what causes this asymmetry.
Together, our data suggest that microtubule dynamics change throughout the spindle during spindle rearrangements. While the magnitude of this effect is spatially varying, throughout the spindle nucleation and turnover rates increase.
Bayesian inference on a mathematical model for microtubule dynamics reveals the relative importance of microtubule severing and growth dynamics for the global spindle dynamics
We next sought to better understand the processes that drive the observed microtubule length changes during spindle rearrangements. For this, we constructed a mathematical model that solves for the expected microtubule length distribution given their growth velocity ${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$ their rate of undergoing catastrophe r, and the rate at which microtubules are cut per unit time and length ${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$. The stability of microtubule plus ends created by cutting is encoded in the parameter ${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$, which interpolates between the two extreme cases, where newly created plus ends immediately depolymerize (${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$ = 0) and where newly created plus ends are indistinguishable from preexisting ones (${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$ = 1; see Figure 5). The mathematical details of this model are given in Material and methods. In spirit, it is very similar to the model proposed in Kuo et al., 2020.
We used this mathematical model to infer the relative importance of microtubule cutting and changes in microtubule nucleation. Using Markov Chain Monte Carlo sampling (ForemanMackey et al., 2013; MacKay, 2003), we infer that the most likely values for the dimensionless ratios ${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$ and ${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$ predict the location of the average microtubule lengths ${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$, given the experimentally determined microtubule length distribution (see Materials and methods). Importantly this inference scheme only uses the experimentally measured microtubule length distribution as input and allows us to quantify the relative importance of cutting to turnover ${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$ independently of direct measurements of turnover rates and microtubule growth velocities, which are hard to do in these very small spindles. We note that additional experiments like FRAP allow us to give absolute time units since it directly measures the turnover rate. The results of the inference scheme are however independent of the FRAP measurements. For details, we refer the reader to Materials and methods.
Using this approach, we first sought to infer the relative importance of turnover to cutting throughout the spindle. In Figure 6, we show the inferred posterior probabilities of model parameters based on the data (Figure 6A–D, Figure 6—figure supplement 1A–D). These plots give the most likely values of the model parameters given the data, 95% confidence intervals for each of the parameters, and 2D error plots for each pair of parameters under the assumption that the third is at its most likely value. We also compared the predictions of the model parameterized by the expectation values of all parameters to our data (Figure 6E–H, Figure 6—figure supplement 1E–H).
We asked whether the transition of spindle structure from metaphase to early anaphase can be explained by microtubule turnover, which here means changes in growth rate and catastrophe, or microtubule cutting. We report our results in terms of the nondimensionalized turnover rate ${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$ and the nondimensionalized cutting rate ${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$, where ${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$ is the average microtubule length. Most importantly, the ratio of these two quantities ${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$ gives the length (in units of average MT length) a microtubule has to reach for it to be more likely to get severed than to undergo catastrophe. This quantifies the relative importance of turnover to cutting; for example, if ${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$ is much smaller than 1, then this indicates a length distribution that is dominated by cutting; if ${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$ is much larger than 1, it indicates a length distribution that is dominated by catastrophe and turnover. In metaphase (Figure 6A, B, E, F, Figure 6—figure supplement 1A, B, E, F) and early anaphase (Figure 6C, G, Figure 6—figure supplement 1C, G), we find that the nondimensionalized turnover rate ${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$ is equal to 1.0, while the nondimensionalized cutting rate ${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$ is about 0.1. Here, ${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$ is the average microtubule length. These numbers imply that the length at which a microtubule is more likely to be cut than to undergo catastrophe is ${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$, that is, cutting is a relatively rare event. In fact, there are no microtubules in metaphase that are longer than 10 times the average microtubule length (i.e., 10 µm) and only one microtubule was larger than 10 times the average length (i.e., 5 µm) in early anaphase. This also implies that the change in the microtubule length observed from metaphase (${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$ = 1 µm) to early anaphase (${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$ = 0.5 µm) is a consequence of either the microtubule turnover rate or changes in microtubule growth velocity. For error estimates and detailed likelihood plots of the inferred parameters, see Figure 6 and Figure 6—figure supplement 1, which give 95% confidence intervals and likelihood surfaces for all parameters. Note also that given that cutting occurs very rarely the inference is not informative on the parameter α, which describes the stability of newly cut microtubules.
We next use our estimate for the ratio ${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$ to infer the microtubule growth velocity in meiotic metaphase spindles. For this, we took the measured FRAP halflife in metaphase (see Figure 2B) of about 5 s as an estimate for microtubule turnover time 1/r and determine the average microtubule length of about 1 µm (Tables 1 and 2). Using these values, we estimated that the microtubule growth velocity (vg) during metaphase was close to 12 µm/min. Interestingly, this agrees with growth velocity estimates based on microtubule flux, and it is similar to values reported for microtubule growth rates in the dense regions of mitotic C. elegans spindles, reaching from 16 to 20 µm/min (Redemann et al., 2017). Notably, this value is slower than the growth rates reported for astral microtubules (44 µm/min) growing away from centrosomes reported by Redemann et al., 2017 and Srayko et al., 2005. This suggests that microtubule growth velocities are under spatial and temporal control within the cytoplasm, as was previously described by Geisterfer et al., 2019 and Walczak et al., 2016.
The drastic change in the average microtubule length from metaphase to anaphase could be due to a decrease in microtubule growth velocity to half its metaphase value (6 µm/min). Alternatively, it could also be caused by an increase in microtubule turnover rate. Given the large relative errors of the FRAP measurements in the very small meiotic spindles, we cannot definitely answer whether microtubule growth speed, turnover rate, or both are under biochemical control at this stage. In toto, our data argue that the transition from metaphase to early anaphase is best explained by modulating microtubule growth and polymerization dynamics, and not by increased microtubule cutting. Note that a role for katanin in amplification of microtubule mass and number has been suggested in vitro, where severing activity can result in the incorporation of GTP tubulin into the microtubule lattice, thus increasing the rescue frequency and stability of newly emerging microtubule ends (Vemu et al., 2018).
We finally asked the same question for the transition from early to mid anaphase. We find that later in anaphase the relative importance of cutting increases with ${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$ = 0.3 and ${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$ = 0.9 (see Figure 6D, H). This means that, at this later point, the length at which a microtubule is more likely to be cut than to undergo catastrophe decreases to 3 ${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$, which would be approximately 3% of the microtubule population.
Thus, our data suggest that the initial steps of the transition from metaphase to anaphase are due to changes in microtubule turnover rate and growth and not mediated by katanindependent severing. This is surprising since earlier work clearly demonstrated the importance of severing during spindle assembly (Srayko et al., 2006; McNally and McNally, 2011; McNally et al., 2014; Connolly et al., 2014; Joly et al., 2016). Estimates of the number of kataninmediated cutting events by counting lateral defects in partial EM reconstructions of meiotic spindles at earlier stages had found a large number of cutting events in wildtype spindles (Srayko et al., 2006). To test the predictions of our inference scheme, we decided to look for similar cutting sites in our datasets.
Analyzing the frequency of lateral defects in our tomographic data indicated a very low abundance (1.1% of microtubules in metaphase #2, 2.5% in early anaphase, and 1.8% in mid anaphase). We could also not detect an increased occurrence of lateral defects in longer microtubules or at distinct positions within the spindle. This is notably fewer events than what was reported for earlier spindles. We conclude that the initial steps of the transition from metaphase to anaphase are due to changes in microtubule turnover rate and growth and not mediated by katanindependent severing.
Cutting selectively occurs for microtubules in close chromosome contact
Our finding that global rearrangements of the spindle structure are mainly caused by changes in nucleation and turnover dynamics made us wonder about the role of kataninmediated cutting in these spindles. To investigate this, we separated out microtubules in close contact to chromosomes (defined as being closer than 150 nm) (Figure 7A, B) and analyzed them separately using our inference scheme. The results are shown in Figure 8 and Figure 8—figure supplement 1. In metaphase spindles, the fraction of microtubules in close contact to chromosomes is small. We found 13, 25, 16, and 23% for our four metaphase spindle datasets, respectively. However, in contrast to the global population, they showed strong evidence of cutting and ${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$ > 1 in all four examples. Furthermore, the length distribution of these microtubules, as depicted in Figure 7C, Figure 8E,,F, and Figure 8—figure supplement 1E, F, is clearly not exponential in metaphase, which is also indicative of cutting.
In contrast, in anaphase, the majority of microtubules (about 70%) are in close contact to chromosomes. Despite this, even microtubules in close contact with chromosomes showed a nearly exponential length distribution (Figure 7D, Figure 8G, H, and Figure 8—figure supplement 1G, H), and our inference scheme yields the same results for microtubules near chromosomes and the global microtubule population. The plots for microtubules that were located further than 150 nm are comparable to the plots for all microtubules (Figure 6) and are shown in Figure 8—figure supplement 2.
Together, these data suggest a role for cutting in maintaining the metaphase spindle. In particular, we hypothesize that katanin selectively cuts microtubules that are in close contact with chromosomes. This mechanism seems to be suppressed in late anaphase.
Katanin depletion reduces the fraction turnover but not the turnover rate of microtubules
We next sought to investigate the role of katanin on microtubule turnover in the metaphase spindle. As null mutants of the C. elegans katanin homologs MEI1 and MEI2 do not assemble bipolar spindles (Mains et al., 1990), we decided to work with the FM13 mei2(ct98) strain, which has an overall lower expression of MEI2 (Srayko et al., 2000; O'Rourke et al., 2011) and reduced microtubulesevering activity (McNally et al., 2006; McNally et al., 2014). We used FRAP experiments to measure the effect of katanin mutation on microtubule turnover using the FM13 mei2(ct98) strain (Figure 9A, B, Video 2). We found that the characteristic time scales for recovery were not significantly affected (4.5 s ± 1.8 s [n = 7] in control and 2.8 s ± 3.4 s [n = 11] in mei2(ct98), p=0.2434). However, strikingly we observed a significant decrease in the fraction of fluorescence recovery of microtubules in the mei2(ct98) embryos (Figure 9C). While in wildtype embryos 87% ± 10% of the fluorescent tubulin signal recovers, this is significantly reduced to 52% ± 26% (p=0.0039) in the katanin mutant. This suggests the presence of a population of stabilized microtubules in katanindepleted spindles. In line with other authors (Brinkley and Cartwright, 1975; Rieder, 1981; Bakhoum et al., 2009), we thus speculate that microtubules in close contact to chromosomes get stabilized by interacting with the chromosomes. Katanin keeps these microtubules from growing too long by cutting these microtubules, and thus, effectively enhances turnover in spindles. In line with this model, our inference data (Figure 8) suggests the new plus ends of newly cut microtubules are mostly unstable. In line with this model, our inference scheme predicts values for ${\partial}_{t}\mathrm{\psi}={\partial}_{l}\left({v}_{g}\mathrm{\psi}\right)\left(r+\mathrm{{\rm K}}l\right)\mathrm{\psi}+\mathrm{{\rm K}}\left(1+\alpha \right){\int}_{l}^{\infty}d{l}^{\text{'}}\psi \left({l}^{\text{'}}\right),$ that are smaller than 0.5 for chromosomeproximal microtubules in all four metaphase spindles. However, given the small number of microtubules in close contact to chromosomes, some significant uncertainty on this number remains.
Discussion
Acentrosomal meiotic spindles in C. elegans undergo a reorganization from a bipolar microtubule arrangement in metaphase to an interchromosomal microtubule array in anaphase. Along with this microtubule reorganization, chromosomes are segregated to extrude half of the genetic material as polar bodies. The underlying mechanism of the microtubule reorganization is not very well understood, and several mechanisms could be involved, for instance, kataninmediated severing as reported for C. elegans meiosis I (Joly et al., 2016; Srayko et al., 2006), transport of microtubules as reported for C. elegans and Xenopus meiotic spindles (Mullen and Wignall, 2017; Brugués et al., 2012), or changes in microtubule polymerization dynamics a shown for Xenopus meiotic spindles (Brugués et al., 2012; Needleman et al., 2010). Here we have generated complete 3D reconstructions of meiotic spindles in C. elegans at different stages and combined this ultrastructural analysis with light microscopy and simulations to investigate the rearrangement of microtubules during meiosis.
Based on light microscopy, previous publications suggested that the microtubule rearrangement from metaphase to anaphase in meiosis I and II in C. elegans could be driven by an initial inward sliding of antiparallel microtubules by kinesin 14 family members (McNally et al., 2016) and a subsequent depolymerization of microtubules at the spindle pole due to severing by katanin (McNally et al., 2006). This was based on the observation that the initial phase of spindle shortening is accompanied by an increase of microtubule density, followed by a further shortening and depolymerization of microtubules at the spindle poles, resulting in a decrease of microtubule density. In agreement with this, our tomographic reconstructions show an initial increase in microtubule number and density during early anaphase, which is followed by a decrease in microtubule number and density in mid anaphase. However, while inward sliding of microtubules would result in an increased density, it does not explain the observed twofold increase of microtubule number. In addition, our data showed a polewarddirected movement of microtubules, contradicting an inward sliding. This suggests that inward sliding is unlikely to contribute to spindle shortening and cannot account for the appearance of microtubules between the chromosomes that characterizes anaphase.
The reorganization of microtubules could alternatively be driven by a selective depolymerization of microtubules at the spindle poles. However, our data suggests that the spindle rearrangement is mainly driven by changes in microtubule nucleation and turnover. A more local analysis of the changes in microtubule length and number, which revealed differences between microtubules near the chromosomes (within 150 nm) and those further away (Figure 7), as well as investigations of local changes using the mathematical model (Figure 8) showed spatial differences in microtubule properties and dynamics.
Our data showed that spindles are made from arrays of dynamic short microtubules, which turnover within 5 s. We further showed that the dramatic structural rearrangements observed from metaphase to anaphase are correlated with drastic changes in the microtubule number and length distribution. Meiotic metaphase spindles are composed of fewer but longer microtubules, while spindles in early anaphase are made of more but shorter microtubules.
These observations led us to ask whether severing of microtubules by katanin, or an increase in microtubule dynamics, allowing more nucleation and/or higher rates of catastrophe, would better explain our observations. For this, we developed a mathematical model that predicts the microtubule length distribution from cutting rates and turnover dynamics. We inferred parameters for the model using the microtubule length distribution and numbers found in electron tomography. This analysis severely constrains the possible mechanisms for spindle restructuring from metaphase to anaphase. Our data suggest that cutting of microtubules in the vicinity of chromosomes is important for maintaining metaphase spindles. It is tempting to speculate that microtubules near the chromosomes are more stable, similar to kinetochore microtubules in mitosis, and that katanin supports the turnover of those more stable microtubules by severing. For the transition from metaphase to anaphase, however, cutting seems to play less of a role. Globally the fraction of microtubules that show signs of cutting is small. Moreover, cutting seems to be much less prominent, even in chromosomeproximal microtubules during anaphase.
Interestingly, a role for microtubulesevering proteins, including katanin, in microtubule amplification has recently been suggested based on in vitro studies. Here, katanin induces nanoscale damage to the microtubule lattice, resulting in the incorporation of GTP tubulin (Vemu et al., 2018; Schaedel et al., 2015; Schaedel et al., 2019) and stabilization of microtubules. The incorporation of GTP tubulin into the lattice is thought to have two effects: promoting rescue as well as the stability of a new plus end formed upon severing. This would result in an amplification of microtubules. Consistent with this, our data shows a twofold increase of microtubule number during the transition from metaphase to early anaphase. Similarly, Srayko et al., 2006 reported a decrease in microtubule number in C. elegans embryos depleted of katanin.
In summary, by combining light microscopy with electron tomography and mathematical modeling we analyzed the reorganization of microtubules during the transition from metaphase to anaphase in C. elegans meiotic embryos. Our data suggests that this reorganization is driven by changes in microtubule growth and/or turnover and that katanin promotes microtubule turnover by severing microtubules near the chromosomes.
Materials and methods
Worm strains and gene silencing by RNAi
Worm strains
Request a detailed protocolThe following C. elegans strains were used in this study: strain MAS91 (unc119(ed3) III; ItIs37[pAA64; pie1::mCherry::HIS58]; ruIs57[pie1::GFP::tubulin+unc119(+)]) for livecell imaging and correlative light microscopy/electron tomography, strain SA250 (tjIs54 [pie1p::GFP::tbb2+pie1p::2xmCherry::tbg1+unc119(+)]; tjIs57 [pie1p::mCherry::his48+unc119(+)]), MAS91 and FM13 (mei2(ct98) I; ruls57 [pAZ147: pie1^{promoter}::tubulin::GFP; unc119(+)]; itls37 [unc119(+) pie1^{promoter}::mCherry::H2B]; him8(e1489) IV), for FRAP experiments. Strains were cultured and maintained at 16°C as described (Brenner, 1974).
Light microscopy
Sample preparation for light microscopy
Request a detailed protocolOocytes for livecell imaging were prepared as described previously (Woog et al., 2012). For FRAP measurements, meiotic spindles in oocytes were observed in the uterus of adult hermaphrodites (strain SA250) mounted on a thin 4% Agarose pad between a slide and a coverslip. Polystyrene microspheres (Microspheres 0.10 μm, Polysciences, Inc) were added to the agar solution before specimen mounting to immobilize the living worms.
Spinning disk confocal fluorescence imaging
Request a detailed protocolLive imaging was performed using a spinning disk confocal microscope (Nikon Ti2000, Yokugawa CSUX1), equipped with 488 nm and 561 nm diode lasers, an EMCCD camera (Hamamatsu), and a 60× waterimmersion objective (CFI Plan Apo VC 60X WI, NA 1.2, Nikon). Acquisition parameters were controlled using a homedeveloped LabVIEW program (LabVIEW, National Instruments). Images were acquired every 1 s with a single zplane.
FRAP experiments
Request a detailed protocolThe photobleaching system was constructed on the abovementioned spinning disk confocal microscope. 80 MHz femtosecond pulsed laser with 0.3 nJ pulse energy and 800 nm center wavelength was used for performing photobleaching and generated from a Ti:sapphire pulsed laser (MaiTai, SpectraPhysics, Mountain View, CA). The photobleaching laser was focused through the same objective for imaging, and photobleaching was performed by moving the sample on a piezostage (P545 PIano XYZ, Physik Instrumente) in three dimensions controlled by a homedeveloped LabVIEW program (LabVIEW, National Instruments). Scanning linebleaching with zsteps was created by moving the stage perpendicular to the poletopole axis back and forth on the focal plane while lowering the stage in the z direction. The parameter for bleaching in length by depth was 6 × 2 µm. The moving speed of the stage was 50 µm/s.
mei2(ct98) FRAP experiment
Request a detailed protocolWorms from the mei2(ct98ts) strain as well as MAS91 worms were kept at 16× until the time of imaging. The FRAP experiment was conducted using a Yokogawa CSUW1 SoRa dual cam spinning disk confocal. This microscope is equipped with seven lasers for imaging (405, 445, 488, 514, 561, 594, 640 nm) that allow nearoptimal excitation of common fluorescent proteins and minimize crosstalk. This spinning disk confocal is mounted on a Nikon Ti2 inverted microscope with a motorized stage for timelapse imaging of multiple stage positions, piezo zdrive for rapid zstack acquisition, and a Perfect Focus module to compensate for stage drift during timelapse imaging. For image acquisition, the microscope is equipped with two types of cameras optimal for different types of experiments, which are mounted on the two camera ports of the CSUW1 dual cam head: (1) a Hamamatsu ORCA Fusion camera for highresolution imaging and (2) dual Photometrics Prime 95B cameras on a Cairn TwinCam for higher speed imaging and/or simultaneous acquisition of two channels. Hardware triggering allows rapid acquisition of zstacks. The microscope is equipped with an Acal BFi UV OptiMicroscan point scanner that we used for the FRAP experiment. This system is integrated with Nikon NIS Elements software for seamless experimental setup and data acquisition. The movies were acquired with a 60 × 1.27 NA water objective and 2.8× SoRa magnifier with 100 ms exposure times and 250 ms intervals (4 frames/s).
Analysis of fluorescence recovery after photobleaching and microtubule poleward flux
Request a detailed protocolFRAP and rate of microtubule poleward flux were calculated with a combination of Fiji (Schindelin et al., 2012) and MATLAB (MATLAB and Statistics Toolbox Release 2012, The MathWorks, Nitick, USA). Timelapse images of spindles expressing GFP::tubulin and mCherry::histone (corresponding to chromosomes) were realigned in a routine for matching, rotation, and translation using Rigid Body of Fiji’s StackReg plugin, so that the random displacement of the spindle due to the spontaneous motion of the worm was corrected.
Poleward flux and recovery of photobleached makers were tracked using a program written in MATLAB (MATLAB and Statistics Toolbox Release 2012, The MathWorks). Line scans of GFPlabeled tubulin along the metaphase spindle were extracted over the course of anaphase. To track the microtubule poleward flux, each line scan at each time point can be divided into two halves by the middle plane of the spindle. The half of the line profile with the bleached mark was subtracted from the other half of the nonbleached profile by a reflection of symmetry around the middle plane of the spindle. This profile subtraction was used to remove spatial variations in the background fluorescence, a valid procedure assuming mirror symmetry of the spindle around its middle plane. A Gaussian function was used to fit the subtracted profile to locate the center of the bleached mark, and thus the position of the bleached mark versus time was extracted. A straight line was fitted to the position of the bleached mark versus time to retrieve the rate of the bleached mark. The recovery time of GFP::tubulin after photobleaching was determined by using an exponential fit. The fluorescence intensities of photobleached marks were calculated by summing intensities over three pixels (~0.5 µm) around the center of the photobleach mark.
Electron microscopy
Sample preparation
Request a detailed protocolSamples for electron tomography were prepared as described (Woog et al., 2012). Briefly, hermaphrodites were dissected in Minimal Edgar’s Growth Medium (Edgar, 1995) and embryos in early meiosis were selected and transferred to cellulose capillary tubes (Leica Microsystems, Vienna, Austria) with an inner diameter of 200 μm. The embryos were observed with a stereomicroscope, transferred to membrane carriers at appropriate stages, and immediately cryoimmobilized using an EMPACT2 highpressure freezer (Leica Microsystems) equipped with a rapid transfer system (Pelletier et al., 2006). Freeze substitution was performed over 3 days at −90°C in anhydrous acetone containing 1% OsO_{4} and 0.1% uranyl acetate using an automatic freeze substitution machine (EM AFS, Leica Microsystems). Epon/Aralditeinfiltrated samples were then embedded in a thin layer of resin and polymerized for 3 days at 60°C. Embedded embryos were remounted on dummy blocks and serial semithick (300 nm) sections were cut using an Ultracut UCT Microtome (Leica Microsystems). Sections were collected on Formvarcoated copper slot grids and poststained with 2% uranyl acetate in 70% methanol followed by Reynold’s lead citrate.
Electron tomography
Request a detailed protocolFor dualaxis electron tomography (Mastronarde, 1997), 15 nm colloidal gold particles (SigmaAldrich) were attached to both sides of semithick sections to serve as fiducial markers for subsequent image alignment. Series of tilted views were recorded using a TECNAI F30 transmission electron microscope (FEI Company, Eindhoven, The Netherlands) operated at 300 kV. Images were captured every 1.0° over a ± 60° range at a pixel size of 2.3 nm using a Gatan US1000 2K × 2K CCD camera. Using the IMOD software package, a montage of 2 × 1 (meiosis I: metaphase #1, metaphase #2, anaphase [late] #1, anaphase [late] #2; meiosis II: metaphase #1, metaphase #2) or 2 × 2 (meiosis I: anaphase [early]; meiosis II: anaphase [late]) frames was collected and combined for each serial section to cover the lengths of the meiotic spindles (Kremer et al., 1996; Mastronarde, 1997).
For image processing, the tilted views were aligned using the positions of the fiducials. Tomograms were computed for each tilt axis using the Rweighted backprojection algorithm (Gilbert, 1972). In order to cover the entire volume of each spindle, we acquired tomograms of about 8–12 consecutive sections per sample. In total, we recorded 10 wildtype spindles in meiosis I and II.
Threedimensional reconstruction and automatic segmentation of microtubules
Request a detailed protocolWe used the IMOD software package (http://bio3d.colorado.edu/imod) for the calculation of electron tomograms (Kremer et al., 1996). We applied the Amira software package for the segmentation and automatic tracing of microtubules (Stalling et al., 2005). For this, we used an extension to the filament editor of the Amira visualization and data analysis software (Redemann et al., 2017; Redemann et al., 2014; Weber et al., 2012). We also used the Amira software to stitch the obtained 3D models in z to create full volumes of the recorded spindles (Redemann et al., 2017; Weber et al., 2014). The automatic segmentation of the spindle microtubules was followed by a visual inspection of the traced microtubules within the tomograms. Correction of the individual microtubule tracings included manual tracing of undetected microtubules, connection of microtubules from section to section, and deletions of tracing artifacts (e.g., membranes of vesicles). Approximately 5% of microtubules needed to be corrected (Redemann et al., 2017).
Data and error analysis
Request a detailed protocolData analysis was performed using either the Amira software package or by exporting the geometric data of the traced microtubules followed by an analysis using MATLAB (MATLAB and Statistics Toolbox Release 2012, The MathWorks). In our analysis of spindle structure, the following errors were considered (Redemann et al., 2017). Briefly, during the data preparation and the imaging process, the tomograms are locally distorted. Furthermore, the exposure of the electron beam causes a shrinking of the sample. During the reconstruction of the microtubules, however, the most important errors occur in the tracing and matching process. In addition, the data is again distorted in all directions to align the tomograms. We assumed that this distortion primarily compensates the distortion of the imaging process. For the tracing, the error was previously analyzed for reconstructions of C. elegans centrosomes (Weber et al., 2012). We assumed that the error lies in the same range of 5–10%. In addition, the traced microtubules were manually verified. It is more difficult to estimate the error of the matching algorithm (Weber et al., 2014) since it depends on the local density and properties of the microtubules. The quality of our analysis should be influenced only by minor 3D distortions.
Quantification of tomographic data
Microtubule length and positioning
Request a detailed protocolThe reconstruction algorithm for microtubules in serial section electron microscopy represents each microtubule as a 3D piecewise linear curve of its centerline. Thus, the length of the microtubule is given by the sum of the lengths of the line segments.
Mathematical model and inference
Mathematical model for microtubule length distributions
Request a detailed protocolThe distribution ψ(l) of microtubule lengths l is set by the microtubule growth velocity v_{g}, their turnover rate r, and the rate at which microtubules are severed by proteins like katanin κ. It obeys
Note that the dimensions of κ are per time per length, and thus the longer microtubules are more likely to get cut. The dimensionless parameter α, which can take values from 0 to 1, determines the likelihood that a new microtubule plus end created by severing is stable. After some basic algebra, we find that the steady states of Equation (1) obey
with the boundary conditions
which state that the creation of new microtubules by nucleation and severing balances the loss of microtubules by turnover such that the total number of microtubules N stays a constant.
Here, (l) = $\left(fl\right)=\frac{1}{N}{\int}_{0}^{\infty}dl\text{'}{\psi}^{\text{'}}\left({l}^{\text{'}}\right){l}^{I}$ denotes the average microtubule length.
For given parameters r, κ, α, we can solve for the microtubule length ψ(l) distribution numerically. Our custom written code uses secondorder finite differences and is available from the authors upon reasonable request.
Bayesian inference method for the relative importance of cutting and catastrophe
Request a detailed protocolWe seek to infer the probability distribution P (model data) of the model parameters (r, α, κ), from the electron tomography data. This can be done using Bayes’ formula
where P (datamodel) is the likelihood of the observed data given the model, P (model) is the prior, and P (data) is the finally marginalized likelihood. Since P (data) is in general not easily calculated, we use Markov Chain Monte Carlo sampling – specifically the toolbox (ForemanMackey, 2020; https://github.com/dfm/emcee) – to approach this problem. This has the advantage that we need not provide an expression for P (data). We next discuss the expressions that we use for P (datamodel) and P (model), respectively.
Calculating P (datamodel)
Request a detailed protocolThis term calculates the likelihood that the measured data (the observed microtubule lengths in our case) are a draw from the probability distribution function predicted by the model. To calculate this term, we first solve Equation 2 with the model parameters and obtain the probability distribution $p\left(l\right)=\psi \left(l\right)/{\int}_{0}^{\infty}dl\psi \left(l\right)$, which characterizes the length distribution in the model.
To obtain the likelihood of the data given the model, we sort the data in to 256 bins of equal size distributed from length 0 to 1.1 times the length of the longest microtubule seen in experiment. We then use p(l) to determine the expectation value λ_{i} for the number of microtubules found in bin i. (Note that we constrained the overall nucleation rate in the model to the value which would make the experimentally observed number of microtubules equal to the expected total number of microtubules.) With this, the probability of finding m_{i} microtubules in bin i is Poisson distributed, and thus
Note that during this procedure we ignore data for microtubules that are shorter than 150 nm long since we suspect that many of these are artifacts from the tomography reconstruction method rather than real data points.
Choosing a prior P (model)
Request a detailed protocolWe seek to not overconstrain our search space and thus use a socalled uninformative prior. This means that we assign all models within the parameter ranges $\overline{r}>0,\text{}\overline{\kappa}0$ and $0<\alpha <1$ the same a priori likelihood.
Test on generated test data
View detailed protocolTo test our inference scheme, we generated datasets from known distributions. We generated datasets that have comparable numbers and length of distributions as our experimental data and ran our inference scheme. In Figures 10 and 11, we show the results of this procedure and compare against ground truth for one case without and one case with cutting.
Data availability
Electron microscopy models of microtubules and chromosome surfaces will be made available on Dryad under https://doi.org/10.5061/dryad.x3ffbg7k5. Example data and analysis code is available at https://github.com/SebastianFuerthauer/SpindleRerrangement (copy archived at https://archive.softwareheritage.org/swh:1:rev:d558f2000a3186de466f6fa3491570298d4a3950).
References

Coldlabile and coldstable microtubules in the mitotic spindle of mammalian cellsAnnals of the New York Academy of Sciences 253:428–439.https://doi.org/10.1111/j.17496632.1975.tb19218.x

A new method reveals microtubule minus ends throughout the meiotic spindleJournal of Cell Biology 175:369–375.https://doi.org/10.1083/jcb.200511112

Blastomere culture and analysisMethods in Cell Biology 48:303–321.https://doi.org/10.1016/s0091679x(08)61393x

emcee: the MCMC hammerPublications of the Astronomical Society of the Pacific 125:306.

The reconstruction of a threedimensional structure from projections and its application to electron microscopy. II. direct methodsProceedings of the Royal Society of London. Series B, Biological Sciences 182:89–102.https://doi.org/10.1098/rspb.1972.0068

To err (meiotically) is human: the genesis of human aneuploidyNature Reviews Genetics 2:280–291.https://doi.org/10.1038/35066065

Computer visualization of threedimensional image data using IMODJournal of Structural Biology 116:71–76.https://doi.org/10.1006/jsbi.1996.0013

Chromosome segregation occurs by microtubule pushing in oocytesNature Communications 8:1–11.https://doi.org/10.1038/s41467017015398

BookInformation Theory, Inference and Learning AlgorithmsCambridge university press.

Dualaxis tomography: an approach with alignment methods that preserve resolutionJournal of Structural Biology 120:343–352.https://doi.org/10.1006/jsbi.1997.3919

Katanin controls mitotic and meiotic spindle lengthJournal of Cell Biology 175:881–891.https://doi.org/10.1083/jcb.200608117

A novel chromosome segregation mechanism during female meiosisMolecular Biology of the Cell 27:2576–2589.https://doi.org/10.1091/mbc.e16050331

The spindle assembly function of Caenorhabditis elegans katanin does not require microtubulesevering activityMolecular Biology of the Cell 22:1550–1560.https://doi.org/10.1091/mbc.e10120951

The segmentation of microtubules in electron tomograms using AmiraMethods in Molecular Biology 1136:261–278.https://doi.org/10.1007/9781493903290_12

A Switch in Microtubule Orientation during C. elegans MeiosisCurrent Biology 28:2991–2997.https://doi.org/10.1016/j.cub.2018.07.012

Microtubules selfrepair in response to mechanical stressNature Materials 14:1156–1163.https://doi.org/10.1038/nmat4396

Lattice defects induce microtubule selfrenewalNature Physics 15:830–838.https://doi.org/10.1038/s4156701905424

Fiji: an opensource platform for biologicalimage analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019

MEI1/MEI2 kataninlike microtubule severing activity is required for Caenorhabditis elegans meiosisGenes & Development 14:1072–1084.

BookAmira: A highly interactive system for visual data analysisIn: Hansen C. D, Johnson C. R, editors. The Visualization Handbook. Burlington: Elsevier. pp. 749–767.https://doi.org/10.1016/B9780123875822/50040X

Spatial regulation of astral microtubule dynamics by Kif18B in PtK cellsMolecular Biology of the Cell 27:3021–3030.https://doi.org/10.1091/mbc.e16040254

Automated tracing of microtubules in electron tomograms of plastic embedded samples of Caenorhabditis elegans embryosJournal of Structural Biology 178:129–138.https://doi.org/10.1016/j.jsb.2011.12.004

Architectural dynamics of the meiotic spindle revealed by singlefluorophore imagingNature Cell Biology 9:1233–1242.https://doi.org/10.1038/ncb1643

Centralspindle microtubules are strongly coupled to chromosomes during both anaphase A and anaphase BMolecular Biology of the Cell 30:2503–2514.https://doi.org/10.1091/mbc.E19010074
Decision letter

Thomas SurreyReviewing Editor; Centre for Genomic Regulation (CRG), Spain

Anna AkhmanovaSenior Editor; Utrecht University, Netherlands

Erin TranfieldReviewer; Gulbenkian Institute, Portugal
In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.
Acceptance summary:
Using a combination of electron tomography, fluorescence microscopy and theoretical modelling, this study provides detailed information about microtubule numbers, lengths and turnover in the spindle of the nematode C. elegans at different stages of meiosis. Microtubules are short compared to spindle length, and their number peaks in early anaphase, although the total mass of polymerised tubulin decreases only later. This work provides new insight into the spindle structure at 'single microtubule resolution' and into microtubule dynamics during meiosis.
Decision letter after peer review:
Thank you for submitting your article "Microtubule reorganization during female meiosis in C. elegans" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, one of whom is a member of our Board of Reviewing Editors, and the evaluation has been overseen by Anna Akhmanova as the Senior Editor. The following individual involved in review of your submission has agreed to reveal their identity: Erin Tranfield (Reviewer #2).
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
As the editors have judged that your manuscript is of interest, but as described below that additional experiments are required before it is published, we would like to draw your attention to changes in our revision policy that we have made in response to COVID19 (https://elifesciences.org/articles/57162). First, because many researchers have temporarily lost access to the labs, we will give authors as much time as they need to submit revised manuscripts. We are also offering, if you choose, to post the manuscript to bioRxiv (if it is not already there) along with this decision letter and a formal designation that the manuscript is "in revision at eLife". Please let us know if you would like to pursue this option. (If your work is more suitable for medRxiv, you will need to post the preprint yourself, as the mechanisms for us to do so are still in development.)
Summary:
In this manuscript, findings from tomographic datasets of 10 C. elegans meiotic spindles from metaphase and anaphase (early, mid, and late) spindles (6 MI and 4 MII) are presented. The focus of the manuscript is on the observation that the transition from metaphase to anaphase involves a significant reorganization of the structure in which the number of MTs increases and the mean length decreases 2fold. The authors develop a mathematical model to assess the relative contributions of 1) changes in MT dynamics, and 2) increased MT severing activity to the reorganization phenomenon. The model explains the data by a global change in MT dynamics and, in fact, indicates that MT severing makes hardly any contribution to the MT shortening observed in anaphase. The work is timely and the topic is of great interest; the quality of the EM data is excellent and these data can be expected to become a valuable resources in the field.
Essential revisions:
1. To compare the model with the data, the authors "average away" a large amount of detailed information present in the EM data and make additional simplifying assumptions that may be questioned. For example, it may be an oversimplification to assume monomodal length distributions in the model that can be described by averages. In figure 1B the metaphase spindles look like there are two populations. The situation in anaphase looks even more complicated particularly if there is a surge of nucleation at the start of anaphase generating new short MT. The detailed 3D data sets are simplified down to a single spatial dimension (the spindle axis) and single length estimator (the average). The authors should provide some evidence/do some tests to validate their approach. How sensitive are the predictions of the model to the simplifying assumptions made and to the averaging out of detail?
2. A major weakness of the manuscript is considered to be the lack of experimental test of the prediction of the model which the authors present as their main conclusion. It should be possible to perform FRAP experiments to test the effect of katanin mutants on microtubule turnover to confirm or contradict the main conclusion that the authors derive from their model and that in part argues against pervious work. There is a wellcharacterized (and fast acting) ts allele of MEI1 called mei1(or642) (O'Rourke et al., PLoS One 2011 and McNally et al., MBoC, 2014) that could be used to test the effect of katanin on microtubule turnover by FRAP.
3. The authors should please be a bit clearer which EM data sets are new and which ones were reused from previous work (for example including this information in Table 1). The expectation would be that the information from the new datasets is also used in the theoretical analysis presented in this manuscript. The context to previous work by others could be explained more clearly by being more specific when presenting background in the introduction so that it will be easier to understand what's new and different here compared to previous work (particularly compared to Yu et al. 2019 and Srayko et al. 2006).
4. Technical concerns:
4.1. FRAP analysis: To which extent does flux versus microtubule polymerization/depolymerization contribute to recovery. Is using a monoexponential function to fit the recovery curves justified given that flux may contribute to recovery? How are the FRAP data used in the model? Is the contribution from flux to recovery considered separately from the contribution of polymerization/ depolymerization?
4.2. p.10, 2nd paragraph: Is the observed decrease in average microtubule length really independent of position? What is the factor of decrease as a function of position? Is the notion of global vs local change really fully supported by the data?
4.3. Model: Are microtubule minus ends considered stable after severing? Alpha is introduced, but the authors do not seem to come back to it later. What is it?
4.4. Model: Throughout, it would be useful to provide confidence intervals for the values that the authors extract from their model or provide some other statistical measure for the reliability of the prediction.
4.5. Does the model make the same predictions for meiosis II spindles or is turnover regulation different there?
4.6. Model: On page 9, lines 1517 – the authors claim that if all the dynamic parameters except nucleation do not change then the length distribution should not change. However, if there is a change in nucleation, there will be a shortterm increase in short MT, thereby shifting the length distribution.
https://doi.org/10.7554/eLife.58903.sa1Author response
Essential revisions:
1. To compare the model with the data, the authors "average away" a large amount of detailed information present in the EM data and make additional simplifying assumptions that may be questioned. For example, it may be an oversimplification to assume monomodal length distributions in the model that can be described by averages. In figure 1B the metaphase spindles look like there are two populations. The situation in anaphase looks even more complicated particularly if there is a surge of nucleation at the start of anaphase generating new short MT. The detailed 3D data sets are simplified down to a single spatial dimension (the spindle axis) and single length estimator (the average). The authors should provide some evidence/do some tests to validate their approach. How sensitive are the predictions of the model to the simplifying assumptions made and to the averaging out of detail?
We apologize to the referees who have been misled (by us) into thinking that we only use averaged data to infer rates of cutting and turnover from the data. In this paper we used the measured lengths of all microtubules to infer the relative importance of cutting and microtubule catastrophe for setting the microtubule length distribution. The theory that we use to calculate the length distribution for given rates of cutting and catastrophe is generic. The Bayesian inference method (Markov Chain Monte Carlo Sampling) that we adapted to infer parameters of the theory uses the lengths of all microtubules, and not just the average length. We now better explain this procedure in an expanded methods section.
The referees are however correct in stating that we assumed that all microtubules obeyed the same growth rules (same cutting and catastrophe rates). In the original manuscript we had tested this assumption by including inferred rates separately for filaments whose center of masses where closer to the poles and closer to the chromosomes, and showing that these where not statistically different. We now further subdivided the data, in particular we separated out all microtubules that approach the chromosomes closer than 150nm. This turns out to be a very small subpopulation of microtubules in metaphase, less than 25% of microtubules. Interestingly this small subpopulation has a markedly different length distribution than the other MTs, which shows clear signatures of cutting. We thank the referees for motivating us to push our analysis further in this regard.
2. A major weakness of the manuscript is considered to be the lack of experimental test of the prediction of the model which the authors present as their main conclusion. It should be possible to perform FRAP experiments to test the effect of katanin mutants on microtubule turnover to confirm or contradict the main conclusion that the authors derive from their model and that in part argues against pervious work. There is a wellcharacterized (and fast acting) ts allele of MEI1 called mei1(or642) (O'Rourke et al., PLoS One 2011 and McNally et al., MBoC, 2014) that could be used to test the effect of katanin on microtubule turnover by FRAP.
We thank the referees for this suggestion. We have conducted FRAP experiments in the FM13 mei2(ct98) strain, which has a reduced microtubule severing rate, and included those results in the paper (Figure 9). This experiment revealed that katanin does promote the turnover of microtubules in meiotic C. elegans spindles, possibly by acting on microtubules near the chromosomes (150nm). Together with our new data on the different length distribution of microtubules within 150nm from the chromosomes this additional experiment has helped to significantly improve this manuscript.
3. The authors should please be a bit clearer which EM data sets are new and which ones were reused from previous work (for example including this information in Table 1). The expectation would be that the information from the new datasets is also used in the theoretical analysis presented in this manuscript. The context to previous work by others could be explained more clearly by being more specific when presenting background in the introduction so that it will be easier to understand what's new and different here compared to previous work (particularly compared to Yu et al. 2019 and Srayko et al. 2006).
We thank the referees for this suggestion and have included this information in Table 1 and Table 2. We have also combined the data from meiosis I and II in the main figures to clarify that the information of all available datasets was used for this manuscript. In addition to this we updated the relevant text passages to better explain previous analysis that were done on some of the datasets.
4. Technical concerns:
4.1. FRAP analysis: To which extent does flux versus microtubule polymerization/depolymerization contribute to recovery. Is using a monoexponential function to fit the recovery curves justified given that flux may contribute to recovery? How are the FRAP data used in the model? Is the contribution from flux to recovery considered separately from the contribution of polymerization/ depolymerization?
We thank the reviewers for this comment. In order to distinguish possible effects of microtubule flux from recovery due to microtubule polymerization/depolymerization we have reanalyzed our FRAP data using larger regions for the analysis. For this we analyzed the recovery rate in a 2µm wide box. Within the recovery time of approximately 10s a microtubule could slide 1µm (based on the velocity of the measured poleward flux). This distance of 1µm is much smaller than the box size, arguing that the recovery is mainly caused by microtubule polymerization/depolymerization.
The FRAP results complement the inference scheme, but are not used in it. Our Bayesian inference scheme infers the relative importance of microtubule cutting and catastrophe for the overall length distribution, i.e. the ratio of two time scales. It does so using only the measured lengths of microtubules. In contrast FRAP measures the absolute time scale of microtubule turnover. This additional information can be used to estimate absolute turnover time scales. We clarify this in the model description.
4.2. p.10, 2nd paragraph: Is the observed decrease in average microtubule length really independent of position? What is the factor of decrease as a function of position? Is the notion of global vs local change really fully supported by the data?
We thank the reviewers for raising this point. We have reanalyzed our data and updated the figure 3 accordingly. As the spindles change significantly in length from metaphase to late anaphase we have normalized the spindle length across meiosis in these new plots (Figure 3) to better determine the decrease as a function of position. We still find that the length of microtubules changes everywhere along the spindle axis, however microtubule in the center of the spindles show a bigger change in length than microtubules at the poles.
4.3. Model: Are microtubule minus ends considered stable after severing? Alpha is introduced, but the authors do not seem to come back to it later. What is it?
We included alpha in the theory and sought to infer it from length data. For the global population, given that cutting turns out to be rare, the value of alpha cannot be substantially constrained from the data. This is shown in Figure 6, Figure 6 Supplementary Figure 1 In essence, all values of alpha are equally consistent with the global data.
However, in the revised manuscript we also separately analyzed the MTs contacting the chromosomes (closer than 150nm) and find that these show clear signs of cutting. For these microtubules we infer that alpha is likely closer to 0 than to 1 – meaning that newly created minus ends are most likely unstable, see Figure 8 and Figure 8 – Supplementary Figure 1 and 2 We now discuss this in the appropriate section of the paper.
4.4. Model: Throughout, it would be useful to provide confidence intervals for the values that the authors extract from their model or provide some other statistical measure for the reliability of the prediction.
We give the errors as standard error to the mean in the text. The figures also include values for the 5% and 95% quantile of the distribution of all model parameters. Finally, the density plots in Figure 6 and 8 and Figure 6 – Supplementary Figure 1, Figure 8 – Supplementary Figure 1 and 2 show the confidence surfaces of parameters in their three dimensional phase space spanned by rbar, kappabar, alpha. We now take more time to explain these plots. We also emphasize that this way of reporting error goes beyond the usual standards of our field.
4.5. Does the model make the same predictions for meiosis II spindles or is turnover regulation different there?
We now include the analysis and data for meiosis two. In short, all our results hold for meiosis I and II alike.
4.6. Model: On page 9, lines 1517 – the authors claim that if all the dynamic parameters except nucleation do not change then the length distribution should not change. However, if there is a change in nucleation, there will be a shortterm increase in short MT, thereby shifting the length distribution.
The referee is right that a sudden change in the nucleation rate will in the short term change the length distribution. This change would subsist for time scales that are comparable to MT turnover times – i.e. seconds. In contrast the reconstruction of the filament network progresses slowly – over several minutes. We thus operate under the assumption that whatever we measure is the reflection of a quasisteady state. In other words, we assume that the rate of change of parameters is slow compared to a microtubule lifetime.
https://doi.org/10.7554/eLife.58903.sa2Article and author information
Author details
Funding
Deutsche Forschungsgemeinschaft (MU 1423/31)
 Ina Lantzsch
Deutsche Forschungsgemeinschaft (MU 1423/32)
 Ina Lantzsch
Deutsche Forschungsgemeinschaft (MU 1423/81)
 Erik Szentgyoergyi
Technische Universität Darmstadt (Frauenhabilitation)
 Stefanie Redemann
Natural Sciences and Engineering Research Council of Canada
 Martin Srayko
National Science Foundation (DMR0820484)
 CheHang Yu
National Science Foundation (NeuroNex #1934288)
 CheHang Yu
National Institutes of Health (1R01GM10497601)
 CheHang Yu
Human Frontier Science Program (RGP 0034/201)
 YuZen Chen
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors are thankful to Th. MüllerReichert for continuous support and to S Tulok and Dr A Walther (Core Facility Cellular Imaging, Faculty of Medicine Carl Gustav Carus, TU Dresden) for help with light microscopy. The authors are also grateful to the members of the electron and light microscopy facility at the Max Planck Institute of Molecular Cell Biology and Genetics (MPICBG, Dresden) for technical assistance, J Baumgart (MPIPKS, Dresden) for initial help in data analysis, and F McNally (UC Davis) for the mei2(ct98) C. elegans strain. The authors would like to thank Drs E O’Toole and R McIntosh (Boulder) for a critical reading of the manuscript. IL and ES were supported by funds from the Deutsche Forschungsgemeinschaft (MU 1423/31, 32, and 81 to T MR) and the Human Frontier Science Programme (RGP 0034/2010 to DN and TMR). SR received funding from the Faculty of Medicine Carl Gustav Carus of the TU Dresden (Frauenhabilitationsstipendium). The work of CHY was supported by the National Science Foundation grant (DMR0820484 and NeuroNex #1934288) and National Institutes of Health Grant (1R01GM10497601).
Senior Editor
 Anna Akhmanova, Utrecht University, Netherlands
Reviewing Editor
 Thomas Surrey, Centre for Genomic Regulation (CRG), Spain
Reviewer
 Erin Tranfield, Gulbenkian Institute, Portugal
Version history
 Received: May 14, 2020
 Accepted: May 24, 2021
 Accepted Manuscript published: June 11, 2021 (version 1)
 Version of Record published: June 24, 2021 (version 2)
 Version of Record updated: February 23, 2023 (version 3)
Copyright
© 2021, Lantzsch et al.
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

 2,187
 Page views

 295
 Downloads

 6
 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

 Cell Biology
 Structural Biology and Molecular Biophysics
Previously we showed that 2D template matching (2DTM) can be used to localize macromolecular complexes in images recorded by cryogenic electron microscopy (cryoEM) with high precision, even in the presence of noise and cellular background (Lucas et al., 2021; Lucas et al., 2022). Here, we show that once localized, these particles may be averaged together to generate highresolution 3D reconstructions. However, regions included in the template may suffer from template bias, leading to inflated resolution estimates and making the interpretation of highresolution features unreliable. We evaluate conditions that minimize template bias while retaining the benefits of highprecision localization, and we show that molecular features not present in the template can be reconstructed at high resolution from targets found by 2DTM, extending prior work at lowresolution. Moreover, we present a quantitative metric for template bias to aid the interpretation of 3D reconstructions calculated with particles localized using highresolution templates and fine angular sampling.