Nuclei determine the spatial origin of mitotic waves
Abstract
Traveling waves play an essential role in coordinating mitosis over large distances, but what determines the spatial origin of mitotic waves remains unclear. Here, we show that such waves initiate at pacemakers, regions that oscillate faster than their surroundings. In cell-free extracts of Xenopus laevis eggs, we find that nuclei define such pacemakers by concentrating cell cycle regulators. In computational models of diffusively coupled oscillators that account for nuclear import, nuclear positioning determines the pacemaker location. Furthermore, we find that the spatial dimensions of the oscillatory medium change the nuclear positioning and strongly influence whether a pacemaker is more likely to be at a boundary or an internal region. Finally, we confirm experimentally that increasing the system width increases the proportion of pacemakers at the boundary. Our work provides insight into how nuclei and spatial system dimensions can control local concentrations of regulators and influence the emergent behavior of mitotic waves.
Introduction
Traveling waves are often used in nature to transmit information quickly and reliably over large distances (Cross and Hohenberg, 1993; Tyson and Keener, 1988; Gelens et al., 2014; Beta and Kruse, 2017; Deneke and Di Talia, 2018). For example, action potentials are well known to propagate along the axon of a neuron (Hodgkin and Huxley, 1952), but a wealth of other biological processes have been shown to be coordinated via traveling waves (Winfree, 1987; Dawson et al., 1999; Loose et al., 2008; Chang and Ferrell, 2013; Deneke et al., 2016; Prindle et al., 2015; Bement et al., 2015; Fukujin et al., 2016). In particular, cell cycle oscillations also self-organize via mitotic waves in a spatially extended system (Chang and Ferrell, 2013; Deneke et al., 2016). Such waves that coordinate cell division in space are especially relevant in the large developing eggs (ranging from ≈100 µm to ≈1 mm in diameter) that are laid externally by insects, amphibians, and fish, because they are too large to be synchronized by diffusion alone (see Box 1). While several studies have addressed the potential biochemical mechanisms of mitotic waves (Chang and Ferrell, 2013; Deneke et al., 2016; Vergassola et al., 2018), what determines the spatial origin of mitotic waves remains unclear.
Box 1.
Spatial cell cycle coordination in early frog and fly embryos.
Insects, amphibians, and fish lay their eggs externally. After fertilization, all of these organisms need to go from a single, large cell to several thousands of somatic-sized cells that then further develop into an adult animal. They do so by carrying out multiple rounds of rapid cleavages following fertilization (Foe and Alberts, 1983; Olivier et al., 2010; Farrell and O'Farrell, 2014; Anderson et al., 2017). While in the Xenopus laevis frog embryo a new cell membrane is formed around each nucleus (green dots in figure), this is not the case in the Drosophila fly embryo, leading to a multinucleated syncytium. Due to the large size of these embryos, diffusion is not fast enough to spatially coordinate the cell cycle (Gelens et al., 2014; Deneke and Di Talia, 2018). There is thus a need for an alternative mechanism to coordinate cell cycle processes over large distances. Mitotic waves have been observed in the shared cytoplasm of large Xenopus cells (Chang and Ferrell, 2013) and the syncytium of the Drosophila embryo (Deneke et al., 2016; Vergassola et al., 2018). Such wave-like propagation of the mitotic state is believed to help spatially coordinate cell cycle progression, yet the spatial origin of this mitotic wave remains unclear.
Here, we address this open question using cell-free extracts made from eggs of the frog Xenopus laevis, which exhibit biochemical cell cycle oscillations in vitro that are similar to those found in vivo (Murray, 1991). We find that mitotic waves originate at nuclei, which act as so-called pacemakers, regions that oscillate faster than their surroundings (Kuramoto, 1984). While previous studies have suggested centrosomes or nuclei to serve as pacemakers (Chang and Ferrell, 2013; Ishihara et al., 2014), their role in organizing mitotic waves has not been empirically demonstrated. We provide evidence that nuclei serve as pacemakers, both in the absence and presence of centrosomes. Having the nucleus setting the pace of the cell cycle may help ensure proper DNA replication prior to initiation of mitosis. If the pacemaker were elsewhere, the decision to divide might be decoupled from DNA replication, leading to division occurring before DNA replication completes. We postulate that nuclei can concentrate cell cycle regulators, leading to faster cell cycle oscillations at those nuclear locations. Nuclei and their spatial positioning, which is affected by the spatial dimensions of the system, determine how the cell cycle is coordinated in space and time.
By monitoring mitotic waves in Teflon tubes using time-lapse microscopy (see Box 2), we find that pacemakers are often located near nuclei that are brighter due to increased import of exogeneously added GFP-NLS. We show that the generation of such pacemakers does not require centrosomes and explore the influence of nuclear density and nuclear import strength on cell cycle period and pacemaker wave formation. Based on these observations, we then develop a theoretical model where nuclei play an active role in concentrating cell cycle regulators. This concentration decreases the period of oscillation around the nuclei. Our modeling shows that the distribution of regulators depends on the nuclear positioning and spatial dimensions of the system, with thicker tubes having a larger tendency to concentrate cell cycle regulators at the boundaries (i.e. outer edges of the tube). Using both numerical simulations and experiments, we go on to show that mitotic waves can originate from the system interior or from the system boundary, depending on the spatial dimensions of the system. These observed dynamics are the result of competition between waves originating from different pacemaker regions, where the relative strength of the pacemakers in the interior and at the boundary is determined by the system dimensions.
Box 2.
Reconstituting cell cycle oscillations using cell-free extracts.
Cell-free cycling extracts can be made from thousands of unfertilized Xenopus laevis frog eggs, following the protocol by Murray, 1991. Cycling extracts can be supplemented with green fluorescent protein with a nuclear localization signal (GFP-NLS) and demembranated sperm nuclei. We load them in Teflon tubes of varying diameters, and image them with a confocal microscope. Under these conditions biochemical oscillations persist and drive the spontaneous formation of nuclei in the extract. Such cell cycle oscillations can be observed by the fluorescent nuclei (importing GFP-NLS) that periodically appear (interphase) and disappear (mitotic phase). Similar oscillations can be observed in bright-field and/or by using fluorescently labeled microtubules (HiLyte Fluor 488).
Results
Nuclei serve as pacemakers to organize mitotic waves
We reconstituted mitotic waves in vitro according to Chang and Ferrell (Chang and Ferrell, 2013; Chang and Ferrell, 2018). We loaded cycling extracts in a 100 µm wide Teflon tube and used green fluorescent protein with a nuclear localization signal (GFP-NLS) to image mitotic waves (see Box 2). This approach allows visualization of regular oscillations between interphase and mitotic phase. In interphase, nuclei form spontaneously in the extract supplemented with sperm chromatin. These nuclei then import GFP-NLS. In mitosis, the nuclear envelope breaks down and GFP is no longer localized to nuclei. Mitotic waves can be observed by the disappearance of nuclei in a wave-like fashion. Waves become apparent after a couple of cell cycles and they self-organize so that they emerge from more clearly defined foci (see Figure 1A, Figure 1—video 1). The origin of the wave (point P) was determined as the intersection of straight lines drawn through the points where the nuclei disappear (see orange curve and Figure 1—figure supplement 1). The wave at cell cycle 5–6 was found to propagate with a speed of ∼ 20 µm/min.
We noticed that the mitotic wave originated close to a nucleus that is considerably brighter than the surrounding nuclei (Figure 1A). We hypothesized that a region with higher GFP-NLS intensity correlates with a higher local oscillation frequency, serving as a pacemaker that organizes the mitotic wave. We therefore analyzed the spatial GFP-NLS intensity profile, the spatial profile of cell cycle periods, and the internuclear distance (Figure 1B). As a brighter nucleus has taken up more GFP-NLS, we reasoned that it similarly concentrates cell cycle regulators that lead to a local increase in the cell cycle frequency. We directly correlated this with the local period, which indeed showed that this region oscillated faster (Figure 1B). To further understand why certain nuclei were brighter, we explored whether their environment had any particular characteristics. We characterized the distance between the different nuclei and found that they were typically separated by 150–200 µm (Figure 1—figure supplement 2). However, we found that the brightest nucleus is also most separated from its neighboring nuclei (Figure 1B). This finding is consistent with the idea that nuclei increase their oscillation frequency by concentrating cell cycle regulators, as they have a larger pool of regulators in their surroundings to import. We analyzed the spatial GFP-NLS intensity profile and the internuclear distance for nine other experiments where we could clearly identify nuclei and mitotic waves. Overall, in 90% of the analyzed experiments the pacemaker location was well predicted by the region with the highest GFP-NLS intensity and/or the region where nuclei were most separated from their neighboring nuclei (Figure 1A,B, Figure 1—figure supplement 3, Figure 1—figure supplement 4). The total nuclear GFP-NLS intensity was also found to be a better indicator of the pacemaker location than the nuclear size as indicated by Hoechst staining, or than the GFP-NLS intensity normalized to the Hoechst signal (Figure 1—figure supplement 4).
In order to further test the role of nuclei as pacemakers, we explored alternative markers of mitotic entry that do not rely on the nuclei themselves. We repeated the experiment with a microtubule reporter, using fluorescently labeled tubulin (HiLyte Fluor 488). Figure 1C and Figure 1—video 2 show that mitotic waves are also observed using such a microtubule reporter, as well as in bright-field. With these tools in hand, we set out to test how critical system parameters such as nuclear density and nuclear import strength influence the mitotic wave dynamics.
Nuclear density and nuclear import strength control cell cycle period and mitotic wave speed
We repeated the experiment in tubes of 100 and 200 µm width for two different concentrations of added demembranated sperm nuclei (approx. 60 and 250 nuclei/µl) (Figure 2A,B). We found that extracts with less added sperm nuclei had a faster cell cycle (Figure 2B). Mitotic waves were similarly observed, but the wave speeds were initially faster than in tubes with a higher nuclear density (Figure 2A). The waves then slowed down to similar speeds as in the case with the higher concentration of sperm nuclei. For both nuclear densities we also found that the average cell cycle period increases over time (Figure 2B). Such a correlation of mitotic wave speed with cell cycle duration is consistent with a transition from sweep waves to trigger waves as the cell cycle slows down (Vergassola et al., 2018). An increase in cell cycle period has been linked to a decrease in ATP supply over time (Guan et al., 2018). An additional explanation could be that an increase in cell cycle period is related to increasing levels of DNA as it is replicated (Dasso and Newport, 1990). This would also explain the decreasing period when reducing the concentration of added sperm nuclei.
Interestingly, a decrease in nuclear density did not lead to a big change in the internuclear distance (Figure 1—figure supplement 2I). Instead, it created more and larger regions where nuclei were absent (Figure 1—figure supplement 3), and pacemakers were predominantly found close to these regions (Figure 1—figure supplement 3). Cheng and Ferrell observed a similar transition from a regular pattern of equidistantly spaced nuclei to a system with holes in Xenopus interphase egg extracts when decreasing the concentration of added sperm nuclei (Cheng and Ferrell, 2019). Next, we further decreased the nuclear density (approx. 30 nuclei/µl), such that only few nuclei remained in an entire tube. Here, we used the fluorescent microtubule reporter to visualize the spatial coordination of mitotic entry, while bright-field images were used to track the location of nuclei (Video 1). Mitotic waves were found to originate at the few nuclei present in the tube, and they traveled through the whole tube (several mm) at a speed of approx. 60 µm/min (Video 1, Figure 2C). In the absence of any nuclei in the tube (no added demembranated sperm nuclei), we still observed cell cycle oscillations with periods similar to extracts with low concentrations of demembranated sperm nuclei (Figure 2—figure supplement 1). However, no mitotic waves were observed (Video 1). These experiments underscore the critical role that nuclei play in changing the cell cycle period and organizing mitotic waves.
Centrosomes have also been suggested to serve as pacemakers (Chang and Ferrell, 2013; Ishihara et al., 2014), potentially by concentrating pro-mitotic factors such as Cdc25 and cyclin B (Bonnet et al., 2008; Jackman et al., 2003). Demembranated sperm nuclei are known to have associated centrioles, which give rise to centrosomes that can generate microtubule asters. In order to test whether such centrosomes are critical to generate pacemakers, we added purified DNA to the extracts, which assembled into nuclei (Newmeyer et al., 1986). Mitotic waves were still observed indicating that DNA alone is sufficient to create pacemaker-generated mitotic waves without a need for centrosomes (Figure 2D, Figure 2—video 1).
As we hypothesize that the import of cell cycle regulators into the nucleus locally changes the cell cycle period, we decided to manipulate the nuclear import strength. We used the nuclear import inhibitor importazole, which is an inhibitor of importin- transport receptors. Increasing levels of importazole were found to increase the cell cycle period and slowed down the formation of nuclei (Figure 2E,F). Mitotic waves were still observed with similar speeds for lower concentrations of importazole, while concentrations higher than 60 µM abolished the formation of nuclei and mitotic waves. Increasing inhibition of nuclear import was also found to lead to smaller nuclei with dimmer levels of GFP-NLS (Figure 2G). When nuclei became very small (i.e. for 40 µM importazole), it took long for the extract to start cycling and mitotic waves were lost (Figure 2F, Figure 2—video 2). We also indirectly manipulated nuclear formation by inhibiting the kinesin Eg5 using S-Trityl-L-cysteine (STLC), which interferes with the proper formation of microtubule structures. We found that increasing concentrations of STLC gradually increased the average cell cycle period (Figure 2—figure supplement 2). Here too, nuclei no longer formed and mitotic waves were no longer observed when STLC was present in too high concentrations (approx. 40 µM STLC). Overall, these findings confirm that nuclear import processes are important in organizing mitotic waves. They ensure that nuclei are able to introduce sufficient spatial heterogeneity in cell cycle period to generate clear mitotic waves.
A computational model where nuclei spatially redistribute cell cycle regulators predicts the location of pacemaker regions
Based on our experimental observation showing that brighter nuclei serve as pacemakers, we set out to develop a theoretical model that describes how GFP-NLS and other proteins can be spatially redistributed by nuclei. A sketch illustrating such a model is shown in Figure 3A. The system toggles between interphase and mitosis with a fixed period. During interphase, nuclei form and nuclear proteins (such as GFP-NLS) are actively imported into the nucleus. During mitosis, the nuclear envelope breaks down and proteins are free to diffuse away. We implemented the competing import and diffusion processes using a generic partial differential equation (PDE) model that describes the evolution of the concentration of nuclear protein, such as GFP-NLS (for details on this model, see Appendix 1). These competing processes are relevant for all proteins that localize to the nucleus. For example, it is known that APC/C is mostly localized in the nucleus, and Wee1 and Cdc25 are actively transported between cytoplasm and nucleus during the cell cycle (Baldin and Ducommun, 1995; Arnold et al., 2015). Such relocalization of cell cycle regulators can locally change the cell cycle oscillation frequency. Note that different different proteins can have opposing effects. For example, while increasing activity of Wee1 and APC/C tend to increase the cell cycle oscillation period, increasing Cdc25 activity leads to faster oscillations (Novak and Tyson, 1993; Tsai et al., 2014). Our experiments thus suggest that the overall effect of increasing nuclear import is to decrease the cell cycle period.
We start by studying the simplest case of a single nucleus in the center of a one-dimensional domain (see Figure 3B). We defined the spatial range of attraction around the nucleus to be approx. 100 µm, such that it is consistent with the so-called nuclear domain, a subdomain of the cytoplasm in which spatial constraints show an effect on nuclear growth (Hara and Merten, 2015). In Xenopus cell-free extract, this nuclear domain has a diameter of approx. 170 µm (Hara and Merten, 2015). The term nuclear domain was originally introduced to describe the surroundings of evenly spaced nuclei in syncytial muscle fibers and Drosophila embryos (Landing et al., 1974; Telley et al., 2012). In Drosophila embryos, the nuclear domain (also called energid) is approx. 30 µm for nuclei which are approx. 5–10 µm in diameter (Chen et al., 2012; Telley et al., 2012). In our experiments we find an internuclear distance of approx. 150 µm for nuclei of approx. 40 µm in diameter (Figure 1—figure supplement 2G-I). Figure 3B shows that proteins quickly build up in the nuclear region in the early phase of the import period and then the proteins quickly disperse after nuclear envelope breakdown. As expected, when averaging the concentration profile over one cell cycle, we find that the time-averaged concentration peaks around the nucleus (see red area in Figure 3B), defining a pacemaker at the nucleus.
While a typical cell contains a single nucleus, the cell-free extract experiment shown in Figure 1 consists of many distributed nuclei. From the experimental data, we calculated how far different nuclei are separated from each other, finding that the distance between neighboring nuclei is typically around 150 µm (Figure 1—figure supplement 2G–I). Note that this distance is consistent with the typical size of a nuclear domain in Xenopus cell-free extract as mentioned before. Moreover, the internuclear distance is also consistent with the size of the recently characterized cell-like compartments that self-organize from homogenized interphase egg cytoplasm (Cheng and Ferrell, 2019). Using this information, we carried out simulations where many nuclei are equidistantly distributed over the whole domain. Such a simulation with 15 nuclei in a domain of 2.4 mm is shown in Figure 3C. Similarly as in the case of a single nucleus, the concentration increases during interphase at each nuclear location, while it quickly decreases during mitosis. However, nuclei close to the boundary are found to have a higher average concentration (see blue shaded area in Figure 3C), which corresponds to a stronger pacemaker region at the boundary.
The build-up of regulators at the boundary is mainly attributed to the fact that nuclei in the interior of the domain compete with neighboring nuclei to attract the available proteins, while nuclei close to the boundary only have one such ‘competitor’. In Figure 3D we verify how the number of nuclei in the system affects the average distribution of regulators, keeping the distance between the outer nuclei and the system boundary constant. Starting from the situation with 15 nuclei in Figure 3C (blue), we gradually decreased the number of nuclei in the system. Figure 3D shows that for decreasing numbers of nuclei (nine in green, five in orange, and three in red), the build-up of regulators at the boundary gradually decreases. When only having three nuclei in the system (red), the central nucleus is found to be dominant and the boundary effect is completely lost. Apart from this competition for regulators between neighboring nuclei, the location of the boundary itself could play an important role. We quantified this boundary effect by changing the distance from the outer nuclei to the system boundary (), while keeping the number of nuclei in the system fixed (15 nuclei). Figure 3E shows that initially an increase in the distance to the boundary leads to a larger build-up of regulators at the boundary, but this increase saturates as becomes larger (Figure 3—figure supplement 1). Although the extent to which regulators build up close to the boundary also depends on the model parameters and on the exact nuclear distribution (see Figure 3—figure supplement 2, Figure 3—figure supplement 3, Figure 3—figure supplement 4), it was found to be a robust phenomenon. Interestingly, however, randomly removing a few nuclei within the domain could abolish the build-up of regulators at the boundary. Instead, proteins build up close to the nuclei adjacent to the gaps (Figure 3—figure supplement 3).
Finally, we expanded our model to two spatial dimensions. We considered rectangular domains of varying aspect ratios, keeping one side fixed in length, while varying the other side in width. The long side was chosen the same as in Figure 3C in which we again define 15 nuclei. We then explored the effect of different widths with increasing rows of nuclei, see Figure 3F. The number of rows of nuclei was based on the experimental observation that wider systems support more nuclei and that those nuclei are separated by the same internuclear distance as in the thin tubes (Figure 3—figure supplement 5). Similarly as in the one-dimensional case, we observe that nuclear cell cycle regulators build up at the edges of the domain. This effect was particularly strong along the longest side of the rectangle, and strikingly, it became more pronounced as the width of the domain increased (see Figure 3F, Figure 3—figure supplement 6).
Multiple pacemakers compete to define the direction of mitotic waves
Based on the model in the previous section, we were able to make predictions of how different nuclear patterns can lead to well-defined spatial distributions of cell cycle regulators. However, transitions between interphase (nuclear import) and mitotic phase (nuclear envelope breakdown and diffusion) occurred with a fixed period. Here, we expand the model by introducing a dependence of the cell cycle period on the local concentration of cell cycle regulators (see details in Appendix 1). In this way a spatial heterogeneity in the concentration of cell cycle regulators leads to a corresponding spatial frequency profile. In general, one expects that such spatial heterogeneities in the cell cycle period create multiple waves. These waves typically propagate into the surrounding medium and compete with each other until the pacemaker with the highest frequency ultimately entrains the whole system (Kuramoto, 1984).
We used this model to explore the dynamics of a pattern of 20 equidistantly distributed nuclei in a domain of 4.2 mm. Figure 4A,D shows that on average cell cycle regulators build up close to the boundary, similarly as in Figure 3C. In the current model, however, this build-up of regulators also leads to a decreased cell cycle period at the boundary. Such a pacemaker region close to the boundary then sends out waves that gradually control the whole domain and they travel more quickly for larger diffusion strengths . We then gradually increased the strength of nuclear import of the three most central nuclei, which on average led to an increased concentration of cell cycle regulators here. For moderate increases in nuclear import strength, two waves compete with one another. A boundary-driven wave and a wave coming from the interior of the domain coexist (Figure 4B). Further increasing the nuclear import strength, waves no longer emerged from the boundary and were entirely controlled by the central region of ”bright" nuclei (Figure 4C).
Next, we removed a nucleus from the center of the domain. Previously, for fixed cell cycle periods, we found that removing nuclei abolished the build-up of regulators at the boundary and proteins localized close to the nuclei adjacent to the gaps (Figure 3—figure supplement 3). Figure 4E indeed illustrates that there is an increased concentration of regulators close to the central gap, but a build-up of regulators close to the boundary also persisted, such that two competing waves were found. We then removed two more nuclei from the center (Figure 4F), which caused the central pacemaker region to send out a wave that controlled the whole domain. The fact that increasing nuclear import strengths and the absence of nuclei within a nuclear pattern both lead to the creation of waves from a nearby location is consistent with the experimental observations reported in Figure 1, Figure 1—figure supplement 3).
We wondered whether these dynamics of competing pacemakers are specific to this particular computational model that includes nuclear import and diffusion processes. Therefore, we also implemented known PDE models of cell cycle oscillations (Appendix 2), where we define two pacemaker regions (see Figure 4—figure supplement 1G–I): an internal pacemaker and a boundary pacemaker region. We carried out simulations continuously changing the relative strength of both pacemaker regions by increasing the difference in cell cycle period. We found a gradual transition from boundary-driven dynamics to internal pacemaker-driven dynamics (Figure 4—figure supplement 1). Similar results were found by using the FitzHugh-Nagumo oscillator model, a general model for relaxation-type oscillatory systems (Figure 4—figure supplement 1K, Figure 4—figure supplement 2). Moreover, we found that even more sinusoidal oscillations preserved boundary-driven waves (Figure 4—figure supplement 2). This suggests that the generation of boundary-driven waves is largely independent of the type of oscillations, as long as the oscillation period is decreased close to the boundary.
Our findings underscore the generic character of the dynamics of multiple competing pacemakers. Pacemaker-driven traveling waves, also often referred to as target patterns, have been widely studied and they form thanks to spatial heterogeneities that locally increase the oscillation frequency. The majority of such pacemaker waves were initially observed in chemical reaction-diffusion systems where heterogeneities were introduced as dust particles that locally modified the properties of the medium (Zaikin and Zhabotinsky, 1970; Zhabotinsky and Zaikin, 1973; Tyson and Fife, 1980). These experimental observations triggered many other studies on both traveling waves (Tyson and Fife, 1980; Kopell, 1981; Hagan, 1981; Kuramoto, 1984; Jakubith et al., 1990; Bugrim et al., 1996; Bub et al., 2005; Stich and Mikhailov, 2006) and spiral waves (Jakubith et al., 1990; Bub et al., 2002; Bub et al., 2005) triggered by a pacemaker. The interaction of multiple pacemaker waves has also been analyzed (Kuramoto, 1984; Walgraef et al., 1983; Mikhailov and Engel, 1986; Lee et al., 1996; Kheowan et al., 2007). In general, they propagate into the surrounding medium and compete with each other until the pacemaker with the highest frequency ultimately entrains the whole system (Kuramoto, 1984). The existence of the transition region is therefore somewhat surprising. However, simulating the system for increasingly longer transient times, we find that the transition region where boundary-driven waves and internal pacemaker-driven waves coexist shrinks, suggesting that after infinitely long transients one pacemaker indeed controls the whole domain. Such infinite transient times are, however, less biologically relevant as the early embryonic cell cycle oscillations only persist for about 13 cycles (Box 2). Therefore, one would expect to observe the full range of transient pacemaker dynamics in actual biological systems.
Wider systems lead to boundary-driven mitotic waves
Our modeling leads to several predictions. First, wider systems lead to higher concentrations of cell cycle regulators at the boundary. Such a local decrease of the cell cycle period leads to boundary-driven mitotic waves. Second, systems with intermediate width allow both internally- and boundary-driven pacemakers. Third, sparsely distributed nuclei favor internal pacemakers. Based on these three predictions, we set out to verify them experimentally.
We repeated the experiment in Figure 1 for varying diameters of the Teflon tubes (approximately 100, 200, 300 and 560 µm) for a nuclear concentration of ≈ 250 nuclei/µl. A representative selection of videos corresponding to this set of experiments is shown in Video 2 (for corresponding kymographs, see Figure 5—figure supplement 1). While the thinnest tube shows mitotic waves coordinated by internal pacemakers, mitotic waves are boundary-driven over the whole domain in the thickest tube. This is consistent with the first theoretical prediction that wider systems lead to boundary-driven mitotic waves. Furthermore, Video 2 illustrates that in tubes of intermediate width (200 and 300 µm), boundary-driven waves coexist with mitotic waves that are driven by internal pacemakers. This is consistent with the second theoretical prediction. By analyzing experiments of 49 tubes of varying widths, we found these findings to be consistent (see Figure 5A). While the thinnest tubes have the lowest probability of finding boundary-driven mitotic waves, all of the experiments with the thickest 560 µm tubes showed boundary-driven waves (Figure 5—figure supplement 2). The fraction of experiments with boundary-driven wave dynamics increased smoothly with the tube width. For a more detailed analysis, see Figure 5—figure supplement 3.
Next, we repeated the experiments using a lower concentration of added sperm nuclei (≈ 60 nuclei/µl). This strongly decreased the probability for mitotic waves to originate from the boundary (see Figure 5B). We noticed that the regularity of the nuclear pattern was disrupted due to the decreased amount of nuclei. Consistent with our third theoretical prediction, the absence of neighboring nuclei was found to strengthen nearby pacemaker regions and decreased the likelihood of having pacemaker regions at the boundary (see Figure 1—figure supplement 3).
As boundary-driven waves were especially clear in the thickest Teflon tubes (560 µm), we wondered whether it was important for the system to be wide enough in all three spatial dimensions. In principle, the theory we developed predicts boundary-driven waves to be present in one-dimensional (Figure 3, Figure 4) and two-dimensional (Figure 3F) spatial systems. Therefore, we carried out experiments with droplets of cycling cell-free extracts on Teflon-coated glass slides, providing a thin structure, yet wide in diameter (≈ 1 mm). All such experiments showed that mitosis was coordinated via mitotic waves that originate at the boundary, consistent with the theoretical predictions (see Video 2).
Finally, we analyzed each individual experiment in more detail with the goal to directly link the presence of pacemaker regions (be it at the boundary or internally) to a local increase in GFP-NLS intensity. This analysis confirmed that there is a higher build-up of GFP-NLS intensity towards the boundaries in wider tubes (see Figure 5—figure supplement 4).
Our findings illustrate that the spatial environment has a strong influence on how biological processes self-organize. In particular, increasing the spatial dimensions of the system leads to a higher probability of observing mitotic waves that originate at the boundary of the system. Other studies have also stressed the importance of system size, boundaries, and geometry on self-organization processes. For example, using cell-free frog extracts, cytoplasmic volume was demonstrated to determine the spindle size (Good et al., 2013; Hazel et al., 2013) and the size of the nucleus (Hara and Merten, 2015).
System boundaries (Kopell et al., 1991; Haim et al., 1996; Rabinovitch et al., 2001; McNamara et al., 2016; Bernitt et al., 2017) and system geometry (Wettmann et al., 2018) have been shown to affect the dynamics of traveling waves. In the widely studied amoeba Dictyostelium discoideum, the origin of cAMP waves have been studied in inhomogeneous systems. Waves appear spontaneously in areas of higher cell density with the oscillation frequency of these centers depending on their density (Vidal-Henriquez and Gholami, 2019). In the presence of advection, a boundary-induced instability was found to periodically excite a cAMP wave near the boundary (Vidal-Henriquez et al., 2017). Another well-characterized model organism is the bacterium Escherichia coli, where Min-protein wave patterns help select the site of cell division (Hu and Lutkenhaus, 1999; Raskin and de Boer, 1999). Wave patterns and the location of cell division have been shown to strongly depend on the system size and geometry, both in vivo by deforming cell shape (Männik et al., 2012; Wu et al., 2015; Wettmann et al., 2018) and in vitro by reconstituting Min oscillations in open and enclosed compartments (Zieske and Schwille, 2014; Zieske et al., 2016; Caspi and Dekker, 2016; Wettmann et al., 2018). As thin compartments were gradually increased in length, multiple regions of oscillations were observed (Zieske and Schwille, 2014; Zieske et al., 2016; Caspi and Dekker, 2016; Wettmann et al., 2018). For more complex geometries, many more wave patterns have been observed, such as standing waves, traveling planar and spiral waves, and coexisting stable stationary distributions (Zieske and Schwille, 2014; Zieske et al., 2016; Caspi and Dekker, 2016; Wettmann et al., 2018). While there are similarities with our findings in the Xenopus cell-free extracts, one important difference is that the wave patterns in the Min system are mainly controlled by the spatial dimensions and geometry. In contrast, in our findings the influence of the spatial dimensions are, at least partially, mediated by the nuclei within the oscillatory medium that serve as pacemakers.
Discussion
A crucial task that a developing cell needs to accomplish is the replication of its DNA and, subsequently, cell division. In large cells, which demand spatial coordination in order to accomplish this task, mitotic waves can organize the process. We have demonstrated that nuclei act as pacemakers generating the mitotic waves in Xenopus cell-free extracts. Pacemakers are regions that oscillate faster than their environment, and, as such, initiate traveling waves (Kuramoto, 1984). A nucleus becomes a pacemaker by its ability to import factors into the nucleus and, presumably, concentrate cell cycle regulators. Indeed, we found that pacemakers are often located near nuclei that are brighter due to increased import of exogeneously added GFP-NLS. We built a generic computational model, which showed that the distribution of cell cycle regulators also depends on the nuclear positioning and spatial dimensions of the system. We tested this idea by experimentally exploring the mitotic wave dynamics in cell-free extracts in which we changed the nuclear density and nuclear import strength. In cell-free extracts with only few nuclei, we found that mitotic waves originated at those nuclei and spread through the parts of the extract devoid of nuclei. In the absence of any nuclei in the system, no mitotic waves were observed. Decreasing the nuclear import strength similarly avoided the formation of mitotic waves. Finally, we changed the spatial dimensions of the system, and found that thicker tubes have a larger tendency to concentrate cell cycle regulators at the boundaries, leading to mitotic waves originating at the outer edges of the tubes. Thus, nuclei are central hubs that organize this complex cellular process.
One advantage to having the nucleus control the timing of mitosis is that it allows the cell to ensure that DNA replication has completed before initiating mitosis. While DNA checkpoints are largely silenced in the early Xenopus embryo (Newport and Dasso, 1989), in Drosophila DNA content is known to activate the DNA-replication checkpoint and alter the cell cycle period (Farrell and O'Farrell, 2014; Deneke et al., 2016). A failure in the correct regulation of mitosis is associated with polyploidy, which plays a key role in nonmalignant physiological and pathological processes (Fox and Duronio, 2013). In the absence of a proper pacemaker, or if the pacemaker were to be located elsewhere, linking DNA replication to mitosis would be more complicated and, perhaps, more prone to error.
Previous studies have pointed to the critical role of the nucleus in spatial redistributing cell cycle regulators (Gavet and Pines, 2010; Santos et al., 2012). In particular, the nuclear import of Cyclin B has been shown to lead to spatial positive feedback, ensuring a robust and irreversible mitotic entry (Santos et al., 2012). Nuclei have also been found to be crucial in ensuring cell cycle oscillations in the Drosophila embryo (Huang and Raff, 1999; Deneke et al., 2019). Interestingly, although previous reports have suggested that centrosomes serve as pacemakers (Chang and Ferrell, 2013; Ishihara et al., 2014), we found that they are dispensable. After treating extracts with purified DNA, which lacks centrosomes, we still observed mitotic waves.
We also found that the interaction of multiple nuclei in a shared cytoplasm can lead to unexpected behavior. Nuclei self-organize in regular spatial patterns within a tube of Xenopus cell-free extract. The measured regular spacing between neighboring nuclei was found to be approximately 150 µm, which coincides with the nuclear subdomain of the cytoplasm in which spatial constraints show an effect on nuclear growth as studied in syncytial muscle fibers (Landing et al., 1974), Drosophila embryos (Telley et al., 2012), and cell-free frog extracts (Hara and Merten, 2015). It is also consistent with the size of cell-like compartments that spontaneously form in homogenized interphase cell-free frog extracts (Cheng and Ferrell, 2019). We found that such regularity in the nuclear distribution led to a build-up of cell cycle regulators towards the boundary of the system, such that the collective behavior of many nuclei creates a pacemaker region at the boundary of the oscillatory medium. This boundary effect was stronger with increasing widths of the tubes, in the presence of more extended regular nuclear patterns. We consistently observed more boundary-driven waves in such wider tubes.
Mitotic waves in the early Drosophila embryo also often originate at the boundary (Foe and Alberts, 1983). During nuclear cycles 10–13 in the syncytial blastoderm of these early embryos, nuclei enter (and exit) mitosis in waves that originate from the opposite anterior and posterior poles of the embryo and terminate in its mid-region. While mitotic waves are associated to so-called trigger waves in the Xenopus embryo (Chang and Ferrell, 2013; Gelens et al., 2014), they have been shown to be so-called sweep waves in the Drosophila embryo (Vergassola et al., 2018). We find, by computational modeling, that sweep waves are also able to generate boundary-driven waves in a syncytium, and that they propagate faster than trigger waves as predicted by Vergassola et al. (2018); Figure 4—figure supplement 2. However, the internuclear distance of our simulations is significantly larger than the one observed in the more crowded Drosophila embryo, so it remains unclear whether our results can directly extend to that system. Despite the limitations of the model, our work is expected to be relevant for all coenocytes (Ondracka et al., 2018), where waves of mitosis have also been observed (Sears, 1967; Brown et al., 2003).
Nuclei are a natural choice of pacemaker for mitotic waves because they allow for a natural way to link one biological process, DNA replication, with another, mitosis. We hope that our work will further trigger new studies into the origin of pacemakers as the initiation of biological decisions mediated by traveling waves seem to be key in the proper coordination of a biological process. Traveling waves have, for example, also been found to propagate apoptosis (Cheng and Ferrell, 2018), action potentials (Hodgkin and Huxley, 1952), and calcium signals (Stricker, 1999) over large distances. In these systems, defective mitochondria, signals from neighboring neurons, or fertilization serve as the initial trigger to locally activate a wave.
Materials and methods
Numerical integration
Request a detailed protocolAll PDE models are solved by numerical integration using custom-made Fortran scripts. Discretization in time is done with a forward Euler method, while discretization in space is carried out with a central difference method. Data is written to .txt files which are then analyzed in Mathematica. The ODE models (CCO and FHN for Figure 4—figure supplement 1) are directly solved in Mathematica, since computational time is limited to seconds. The numerical codes that were used are available through GitHub (Nolet, 2020).
Experimental setup
Request a detailed protocolWe reconstitute cell cycle oscillations in vitro in cell-free cycling extracts made from unfertilized Xenopus laevis frog eggs, following the protocol by Murray, 1991; Box 2). Female Xenopus laevis frogs are injected subcutaneously with 500 injection units (IU) human chorionic gonadotropin (MSD Animal Health) to induce ovulation, after prior priming with 100 IU pregnant mare’s serum gonadotropin (MSD Animal Health). The obtained eggs are rinsed with deionized water and subsequently their jelly coat is removed by incubation in a 2% w/v cysteine in 1 × XB salts solution. Dejellied eggs are now susceptible to activation with the calcium ionophore A23187 (0.5 µg/mL in 0.2 × Marc’s Modified Ringer’s buffer, Sigma-Aldrich) for 2 min to start the biochemical processes of the cell cycle. After a packing step, the activated eggs are crushed in an ultracentrifuge (XPN90, Optima) at 16,000 × g at 2°C for 10 min. This allows the collection of the cytoplasmic fraction to which the protease inhibitors leupeptin, pepstatin and chymostatin (Sigma-Aldrich) are added to a final concentration of 10 µg/mL. Cytochalasin B (10 µg/mL, Sigma-Aldrich) is also added to inhibit actin assembly and thus gelation-contraction, keeping the extract fluid at room temperature (Field et al., 2011).
Finally, the extract is supplemented with GFP-NLS (∼ 25 µM), green fluorescent protein with a nuclear localization signal, and sperm chromatin (using two different concentrations: ∼ 63 or 250 nuclei/µL extract). The construct for GFP-NLS was kindly provided by James Ferrell (Stanford Univ., USA). Sperm chromatin was prepared according the protocol by Murray, 1991. The supplemented extracts are then loaded in Teflon tubes (Cole-Parmer PTFE, 06417–11), through aspiration, and imaged at 24°C on a Leica TCS SPE confocal fluorescence microscope. This approach allows to visualize regular oscillations between interphase and mitotic phase. In interphase, nuclei form spontaneously in the extract supplemented with sperm chromatin. These nuclei then import GFP-NLS (see Box 2). In mitosis, the nuclear envelope breaks down and GFP is no longer localized to nuclei. Here, we use this experimental system to explore the influence of system size by varying the width of the Teflon tubes. The tubes were approximately 100, 200, 300, and 560 µm in width (the actual inner diameters are 102, 203, 305, and 559 µm). Furthermore, we change the amount of nuclear material and its distribution by considering two different concentrations of added sperm chromatin.
In addition, DNA was purified from the sperm chromatin. This was done using a GenElute Mammalian Genomic DNA kit (Sigma-Aldrich), with the use of proteinase K (Sigma-Aldrich) to release the DNA from the histones and give a higher yield. After purification, the concentration of DNA was determined using a NanoDrop spectrophotometer. Purified DNA was added to the extract at final concentrations of 5, 10, 15, 20, 25, 45 and 60 ng/µL.
Nuclear import was inhibited by adding importazole (Sigma-Aldrich), an inhibitor of importin- transport receptors. Final concentrations of 5, 10, 20, 40, and 60 µM were tested.
Microtubule dynamics was disrupted by adding S-Trityl-L-cysteine (STLC, Acros Organics), a kinesin Eg5 inhibitor. Final concentrations of 10, 20, 30, 40, and 50 µM were tested.
In some of the experiments fluorescent reporters other than GFP-NLS were used. These included a green microtubule reporter (Tubulin porcine HiLyte 488; Cytoskeleton, Inc) at 1 µM final concentration and DNA staining (Hoechst 33342) at 5 µg/mL final concentration.
Image analysis
Microscope data
Request a detailed protocolWe used a Leica TCS SPE confocal fluorescence microscope (5x objective) in confocal mode to excite the GFP-NLS with a 488 nm solid state laser, and capture the emission from 493 to 600 nm. In the non-confocal experiments we used the Leica EL6000 metal halide external fluorescence light source for excitation of the fluorophores. The different filter cubes used were the L5 (excitation 480/40 nm bandpass, emission 527/30 nm bandpass) for GFP-NLS and HiLyte Fluor 488; and the A4 (excitation 360/40 nm bandpass, emission 470/40 nm bandpass) for the Hoechst 33342 staining. First, we fixed imaging positions at different (x,y) locations of the Teflon tubes, ensuring overlap between subsequent positions to capture the whole tubes. Within a tube, the z-position was fixed, but could differ between tubes to be able to image the central plane of the tubes. We then captured time-lapse images of these different positions during 18 hr, creating image stacks for each position in a .lif (Leica Image File) format. The .lif files belonging to one tube were then imported in Fiji (Schindelin et al., 2012). The maximum intensity of the different image stacks was put at the same level. Then, using the overlap between subsequent image positions, the image stacks were stitched pairwise (Preibisch et al., 2009). Subsequently, the images were cropped and saved as separate .tiff files per timepoint, an .avi file and a kymograph were made.
Data analysis from images
Request a detailed protocolThe .tiff files are imported in Mathematica and for all the maximum intensity over the width is calculated. This allows us to have a one-dimensional intensity profile for each time, see Figure 1—figure supplement 1C. Kymographs as in Figure 1A and Figure 1—figure supplement 3 were made from these profiles over time. Lines are drawn through the points of mitotic entry (disappearance of nuclei), for every visible cycle. This is done by manually detecting the start- and endpoints of the wave, as depicted in the sketch of Figure 1—figure supplement 1D. The lines are drawn through those points automatically and periods and wave speeds are then calculated based on these lines. The period is calculated by taking 20 points on these lines and determining the time to the next line. This gives an average period (and standard deviation) for each cycle. The wave speed is calculated by taking the derivative of the lines. For the full cycle, the wave speed is only reported if the wave travels a large enough (> 600 µm) distance (to only include well-formed waves and to reduce noise), and if multiple waves are present, the minimum speed is reported. The locations of the nuclei (one-dimensional) are extracted from the kymographs at the last one or two lines (if nuclei are well-separated). For each nucleus the average distance to their neighbors (left and right) is calculated which is also plotted in Figure 1C and Figure 1—figure supplement 3. For the last two cycles, the maximum intensity over the cycle is calculated at every , yielding an intensity profile at each cycle.
Processing for specific analyses
Request a detailed protocolWhen calculating properties of individual nuclei (e.g. size, location, intensity), the Ilastik software was used to automatically recognize nuclei in a series of .tiff files. This program relies on machine learning software which makes recognition a lot faster than manual tracking. The files are imported in Ilastik, where we provided three labels (’nucleus’, ’background’ or ’outside of the tube’) to train the implemented random forest classifier to recognize the labels in the images (Sommer et al., 2011). After the training phase, we exported the results as a .hdf5 file, which contains the probability of each pixel to be ’nucleus’, ’background’ and ’outside of the tube’ for each timepoint. The .hdf5 files were imported in Mathematica for further analysis. The data of these files was binarized by defining all pixels with a high probability (≥ 75%) as nuclei (1) and others as background (0). Adjacent pixels were grouped together and the separate groups were recognized as the nuclei. Noise was reduced by ignoring nuclei consisting only of a few pixels. This resulted in a binarized picture, such as in Figure 1—figure supplement 1A. Of all recognized nuclei (orange), information as location (center) and size is extracted with Mathematica. In order to obtain continuous-time kymographs (such as in Figure 1A and Figure 5—figure supplement 1, we overlayed the binarized matrix with the original .tiff and integrated over the width. In this way intensity differences were still visible.
Analysis of the pacemaker strength of internal regions and the boundary regions
Request a detailed protocolThe GFP-NLS intensity profile of the experiments is analyzed in order to calculate the strength of the boundary and of internal pacemakers (Figure 5—figure supplement 1 and Figure 5—figure supplement 4). An example of such an intensity profile is shown in Figure 1—figure supplement 1B. The averaged intensity profile is filtered using a low-pass filter, to obtain a ‘background’ signal . This is the red line in Figure 1—figure supplement 1B. All frequencies higher than a threshold are filtered out. The obtained background profile does of course depend on the parameter . The position of the minimum of is denoted by , that is
From the background profile, we calculate two measures for the GFP build-up at the boundary, by
and
These correspond to the GFP build-up in the blue areas in Figure 5—figure supplement 1.
A second parameter, , is introduced and defines the boundary width. In other words, the intervals and are the boundary domains and is the internal domain. The background profile might over- or underestimate GFP build-up in the boundary domains. This is compensated by calculating the second type of measures, and . These are defined by
and
The GFP build-up at the boundary, denoted by , of this intensity profile is now defined as
The internal GFP build-up (i.e. by nuclei located internally) is defined by those areas where the intensity is higher than the background profile . This internal GFP build-up is calculated by
which correspond to the orange areas in Figure 5—figure supplement 1.
Figure 5—figure supplement 4 shows the GFP build-up at the boundary and internally, vs. , for 20 experiments. This is done for various values of and . Since and depend on these parameters, the figure will change with those parameters. However, we see that qualitatively differences are small.
Data availability
Request a detailed protocolAll the data generated during the study are summarized and provided in the manuscript and supporting files. Source files have been provided for Figure 1, Figure 1—figure supplement 4, Figure 2, Figure 5—figure supplement 1, Box 2, Video 1 and Video 2 in the format of microscopy videos. Additionally, representative microscopy videos of all different conditions are provided as a Zenodo dataset (http://doi.org/10.5281/zenodo.3736728). The numerical codes that were used, together with an overview table of the performed experiments, are available through GitHub (Nolet, 2020; copy archived at https://github.com/elifesciences-publications/eLife_paper).
Appendix 1
Modeling import and diffusion processes
A model of import and diffusion of cell cycle regulators with a fixed cell cycle period
We implemented the competing import and diffusion processes through a simple PDE equation that describes the evolution of the concentration of a nuclear protein:
with a diffusion constant, a periodic function that periodically switches between 0 (in mitosis) and 1 (in interphase),
with the cell cycle period, the ratio of interphase to the cell cycle period and a potential function to define nuclear attraction at different locations. For simplicity, we use a sum of Gaussian potential wells throughout this work,
with attracting positions at (the ‘nuclei’), and parameters that define the amplitude and the width of the potential wells, respectively. The analysis is also repeated with a different potential function, giving qualitatively similar results.
Rewriting the one-dimensional model equation as
we decrease numerical errors by explicitly implementing the first and second derivative of . With basic calculus we find for a Gaussian potential function,
and
For the appropriate boundary conditions we need to calculate the flux. Let be the flux of , then we have
and when rewriting the model equation as
we find the flux to be
and the zero-flux boundary condition is then , with the normal vector at the boundary. In one dimension we will have a boundary condition at given by
for all and for the numerical derivative of this becomes
where denotes the numerical grid distance. This yields
for which we use the explicit formula for . Similarly, one can find at ,
A model of import and diffusion of cell cycle regulators with a concentration-dependent cell cycle period
The previous model is uniform in time, that is the cell cycle is regulated by the function with a given length of S phase and M phase. Here we let the period of each nucleus depend on the local concentration of cell cycle regulators. This is done by including a separate function in the potential for every nucleus and letting the period depend on . This means the period is not a fixed constant anymore, every nucleus has its own that depends on (and therefore on time).
The model equation becomes
and recall that every nucleus defines a potential and together they form the full potential function
via superposition. The single-nucleus potential can now be written as
i.e. it can be separated in a time-dependent and a space-dependent function. The function is the same Gaussian well as before, that is
where is a parameter that defines the width of the potential function (i.e. the attraction range of the nucleus) and is the amplitude of the potential (i.e. the attraction strength). The parameter is nucleus-dependent and can be used to give different strengths to different nuclei. Normally, when all nuclei have the same attraction strength, for all . The function regulates the potential over time, essentially by turning it on and off. The definition is
where and are time-dependent via
for , that is denotes the last for which and the last for which . Therefore the function depends on its history and we define . The parameter determines the fraction of the cycle where nuclei are in attracting phase. The length of the cycle, denoted by the period function , is not constant but changes in time. This is because it is linked to the concentration at the position of the nucleus via
where is a reference period (chosen the same as in the previous model), a reference concentration and is a parameter that determines how strong the period is coupled to the concentration at . This time-dependent period is concentration-dependent (which on itself is time-dependent) and has as a consequence that nuclei with a higher concentration have a shorter period (for ). The underlying biological assumption is that nuclei that import (attract) more cell cycle regulators have a shorter cell cycle length than other nuclei.
Boundary conditions are the same as before. The standard values of the model parameters are given in the table below.
Par. | Value | Unit |
---|---|---|
D | 600 | µm2/min |
T | 40 | min |
0.7 | ||
1 | ||
1 | ||
1.1 | ||
4 · 104 | µm3/min | |
60 | µm |
Appendix 2
Cell cycle models
Biochemical cell cycle oscillator
Mathematically, the cell cycle oscillator (CCO) based on the regulatory network in Figure 4—figure supplement 1A–C can be described using the following ODEs for cyclin B concentrations ([cyc]) and concentrations of the cyclin B-Cdk1 complex in its active state ([cdk1]),
with the synthesis rate of cyclin B, and parameters corresponding to the influence of Cdc25 (), Wee1 () and degradation via APC/C (). For a more detailed discussion of this model, including an experimental motivation and/or measurement of the different parameter values, we refer to Chang and Ferrell, 2013 and references therein.
Generic FitzHugh-Nagumo oscillator
Although the CCO model (Equation (28)) includes known biochemical interactions, much of the mitotic wave dynamics solely results from the fact that the cell cycle is driven by a relaxation-type oscillator, such as the generic FHN oscillator (Fitzhugh, 1961). The FHN system was originally constructed as a simplified version of the Hodgkin-Huxley model, describing the activation and deactivation dynamics of a spiking neuron (Hodgkin and Huxley, 1952). Similarly, the FHN model has later been used to describe pulse dynamics in the heart (Aliev and Panfilov, 1996) and cell cycle oscillations (Gelens et al., 2014; Gelens et al., 2015). Using such a simple model to describe complex biological processes has the advantage of being numerically more effective, allowing for analytical estimations, and having fewer parameters.
We use a generalized form of the classical FHN oscillator, including two ODEs for variables and ,
with parameters and . When taking , the original FHN equations are recovered. Figure 4—figure supplement 1E,F show that the relaxation oscillations in the FHN are indeed very similar to the cell cycle oscillations, after applying a linear mapping (chosen by trial-and-error to obtain visual correspondence).
Spatially extended models
In order to capture spatial dynamics, we extended the CCO model (Equation (28)) and FHN model (Equation (29)) by including free diffusion of proteins (with diffusion constants µm2/min), leading to the following coupled partial differential equations (PDEs) for the CCO and FHN system, respectively:
and
Parameter values
Standard parameters of both the CCO (top) and FHN (bottom) models are given in the table below. When other values are used for certain figures, this is stated in the corresponding caption. Note that since the FHN is a dimensionless model, no units are reported.
Parameter | Value | Unit | Parameter | Value | Unit |
---|---|---|---|---|---|
0.8 | min−1 | 4 | min−1 | ||
0.4 | min−1 | 2 | min−1 | ||
0.01 | min−1 | 0.06 | min−1 | ||
35 | nM | 11 | |||
30 | nM | 3.5 | |||
32 | nM | 17 | |||
1.5 | nM/min | 600 | µm2/min | ||
−0.85 | 0.05 | ||||
1.2 | 0.5 | ||||
0.01 | 600 |
Modeling sweep waves in the FHN model
Sweep waves are caused by moving the nullcline in time, which increases the frequency and can also increase the wave speed in an already oscillating model. In the FHN model, this is implemented by adding an extra term, , as follows:
This term in the -equation moves the nullcline (Figure 4—figure supplement 1F) to the left for positive . The term depends on time: it is zero in M phase and increases in S phase. We distinguish two types of sweep waves, local and global. In the local case, is altered at every , that is . In the global case, only depends on time but is uniform in space. The increase in S phase is with a certain speed, and then takes the form . The fraction 1/30 is for convenience (since the duration of S phase is around 30 min one can easily check how far the nullcline is sweeped).
Data availability
All the data generated during the study are summarized and provided in the manuscript and supporting files. Source files have been provided for Figure 1, Figure 1-figure supplement 3, Figure 2, Figure 5-figure supplement 1, Box 2, Video 1 and Video 2 in the format of microscopy videos. Additionally, representative microscopy videos of all different conditions are provided as a Zenodo dataset (http://doi.org/10.5281/zenodo.3736728). The numerical codes that were used, together with an overview table of the performed experiments, are available through GitHub (https://github.com/felixnolet/eLife_paper; copy archived at https://github.com/elifesciences-publications/eLife_paper).
References
-
A simple two-variable model of cardiac excitationChaos, Solitons & Fractals 7:293–301.https://doi.org/10.1016/0960-0779(95)00089-5
-
Insights into the cellular mechanism of the yeast ubiquitin ligase APC/C-Cdh1 from the analysis of in vivo degronsMolecular Biology of the Cell 26:843–858.https://doi.org/10.1091/mbc.E14-09-1342
-
Subcellular localisation of human wee1 kinase is regulated during the cell cycleJournal of Cell Science 108 ( Pt 6:2425–2432.
-
Intracellular oscillations and wavesAnnual Review of Condensed Matter Physics 8:239–264.https://doi.org/10.1146/annurev-conmatphys-031016-025210
-
Spiral wave generation in heterogeneous excitable mediaPhysical Review Letters 88:058101.https://doi.org/10.1103/PhysRevLett.88.058101
-
Global organization of dynamics in oscillatory heterogeneous excitable mediaPhysical Review Letters 94:028105.https://doi.org/10.1103/PhysRevLett.94.028105
-
Heterogeneous sources of target patterns in reaction−diffusion systemsThe Journal of Physical Chemistry 100:19017–19022.https://doi.org/10.1021/jp961603t
-
Robustly cycling Xenopus laevis Cell-Free Extracts in Teflon ChambersCold Spring Harbor Protocols 2018:pdb.prot097212.https://doi.org/10.1101/pdb.prot097212
-
Pattern formation outside of equilibriumReviews of Modern Physics 65:851–1112.https://doi.org/10.1103/RevModPhys.65.851
-
Chemical waves in cell and developmental biologyJournal of Cell Biology 217:1193–1204.https://doi.org/10.1083/jcb.201701158
-
Actin behavior in bulk cytoplasm is cell cycle regulated in early vertebrate embryosJournal of Cell Science 124:2086–2095.https://doi.org/10.1242/jcs.082263
-
Studies of nuclear and cytoplasmic behaviour during the five mitotic cycles that precede gastrulation in Drosophila embryogenesisJournal of Cell Science 61:31–70.
-
Self-organization of chemoattractant waves in Dictyostelium depends on F-actin and cell–substrate adhesionJournal of the Royal Society Interface 13:20160233.https://doi.org/10.1098/rsif.2016.0233
-
Activation of cyclin B1-Cdk1 synchronizes events in the nucleus and the cytoplasm at mitosisThe Journal of Cell Biology 189:247–259.https://doi.org/10.1083/jcb.200909144
-
Spatial trigger waves: positive feedback gets you a long wayMolecular Biology of the Cell 25:3486–3493.https://doi.org/10.1091/mbc.e14-08-1306
-
Target patterns in reaction-diffusion systemsAdvances in Applied Mathematics 2:400–416.https://doi.org/10.1016/0196-8858(81)90042-7
-
Breathing spots in a Reaction-Diffusion systemPhysical Review Letters 77:190–193.https://doi.org/10.1103/PhysRevLett.77.190
-
Organization of early frog embryos by chemical waves emanating from centrosomesPhilosophical Transactions of the Royal Society B: Biological Sciences 369:20130454.https://doi.org/10.1098/rstb.2013.0454
-
Active cyclin B1-Cdk1 first appears on centrosomes in prophaseNature Cell Biology 5:143–148.https://doi.org/10.1038/ncb918
-
Wave mediated synchronization of nonuniform oscillatory mediaPhysical Review Letters 98:074101.https://doi.org/10.1103/PhysRevLett.98.074101
-
Target pattern solutions to reaction-diffusion equations in the presence of impuritiesAdvances in Applied Mathematics 2:389–399.https://doi.org/10.1016/0196-8858(81)90041-5
-
On chains of oscillators forced at one endSIAM Journal on Applied Mathematics 51:1397–1417.https://doi.org/10.1137/0151070
-
BookChemical Oscillations, Waves, and TurbulenceBerlin: Springer-Verlag.https://doi.org/10.1007/978-3-642-69689-3
-
Competing patterns of signaling activity in Dictyostelium discoideumPhysical Review Letters 76:1174–1177.https://doi.org/10.1103/PhysRevLett.76.1174
-
Multiple target pattern creation and synchronization phenomenaPhysics Letters A 117:257–260.https://doi.org/10.1016/0375-9601(86)90088-5
-
BookCell Cycle ExtractsIn: Kay B. K, Peng H. B, editors. Methods in Cell Biology, 36. Academic Press. pp. 581–605.https://doi.org/10.1016/S0091-679X(08)60298-8
-
On the coupling between DNA replication and mitosisJournal of Cell Science 1989:149–160.https://doi.org/10.1242/jcs.1989.Supplement_12.13
-
Modeling the cell division cycle: m-phase trigger, oscillations, and size controlJournal of Theoretical Biology 165:101–134.https://doi.org/10.1006/jtbi.1993.1179
-
Inwards propagating waves in a limit cycle mediumPhysical Review Letters 87:084101.https://doi.org/10.1103/PhysRevLett.87.084101
-
Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
-
ConferenceIlastik: interactive learning and segmentation toolkit. in:biomedical imaging: from nano to macro, 2011IEEE International Symposium on IEEE.https://doi.org/10.1109/ISBI.2011.5872394
-
Target patterns in two-dimensional heterogeneous oscillatory reaction–diffusion systemsPhysica D: Nonlinear Phenomena 215:38–45.https://doi.org/10.1016/j.physd.2006.01.011
-
Comparative biology of calcium signaling during fertilization and egg activation in animalsDevelopmental Biology 211:157–176.https://doi.org/10.1006/dbio.1999.9340
-
Aster migration determines the length scale of nuclear separation in the Drosophila syncytial embryoThe Journal of Cell Biology 197:887–895.https://doi.org/10.1083/jcb.201204019
-
Target patterns in a realistic model of the Belousov–Zhabotinskii reactionThe Journal of Chemical Physics 73:2224–2237.https://doi.org/10.1063/1.440418
-
Singular perturbation theory of traveling waves in excitable media (a review)Physica D: Nonlinear Phenomena 32:327–361.https://doi.org/10.1016/0167-2789(88)90062-0
-
Convective instability and boundary driven oscillations in a reaction-diffusion-advection modelChaos: An Interdisciplinary Journal of Nonlinear Science 27:103110.https://doi.org/10.1063/1.4986153
-
Spontaneous center formation in Dictyostelium discoideumScientific Reports 9:3935.https://doi.org/10.1038/s41598-019-40373-4
-
Chemical waves in a two‐dimensional oscillating systemThe Journal of Chemical Physics 78:3043–3051.https://doi.org/10.1063/1.445266
-
BookWhen Time Breaks Down: The Three-Dimensional Dynamics of Chemical Waves and Cardiac ArrhythmiasPrinceton, NJ: Princeton University Press.https://doi.org/10.1063/1.2811674
-
Symmetry and scale orient min protein patterns in shaped bacterial sculpturesNature Nanotechnology 10:719–726.https://doi.org/10.1038/nnano.2015.126
-
Autowave processes in a distributed chemical systemJournal of Theoretical Biology 40:45–61.https://doi.org/10.1016/0022-5193(73)90164-1
-
Protein Patterns and Oscillations on Lipid Monolayers and in MicrodropletsAngewandte Chemie International Edition 55:13455–13459.https://doi.org/10.1002/anie.201606069
Article and author information
Author details
Funding
Research Foundation - Flanders (GOA5317N)
- Lendert Gelens
KU Leuven Research Fund (C14/18/084)
- Lendert Gelens
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Jim Ferrell, Sophie De Buyl, and Jan Rombouts for valuable feedback on the manuscript. We also thank Jonás Noguera López, Ine Vlaeminck and Virginia Tsiouri for their help in the lab and stitching movies. This work was supported by the Research Foundation - Flanders (grant GOA5317N) and the KU Leuven Research Fund (C14/18/084).
Ethics
Animal experimentation: This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the KU Leuven. All of the animals were handled according to approved institutional animal care and use committee (IACUC) protocols of the KU Leuven. The protocol was approved by the Committee on the Ethics of Animal Experiments of the KU Leuven (ECD permit Number: P165/2016/).
Copyright
© 2020, Nolet 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,949
- views
-
- 365
- downloads
-
- 30
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
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
-
- Computational and Systems Biology
- Neuroscience
The basolateral amygdala (BLA) is a key site where fear learning takes place through synaptic plasticity. Rodent research shows prominent low theta (~3–6 Hz), high theta (~6–12 Hz), and gamma (>30 Hz) rhythms in the BLA local field potential recordings. However, it is not understood what role these rhythms play in supporting the plasticity. Here, we create a biophysically detailed model of the BLA circuit to show that several classes of interneurons (PV, SOM, and VIP) in the BLA can be critically involved in producing the rhythms; these rhythms promote the formation of a dedicated fear circuit shaped through spike-timing-dependent plasticity. Each class of interneurons is necessary for the plasticity. We find that the low theta rhythm is a biomarker of successful fear conditioning. The model makes use of interneurons commonly found in the cortex and, hence, may apply to a wide variety of associative learning situations.
-
- Cancer Biology
- Computational and Systems Biology
Effects from aging in single cells are heterogenous, whereas at the organ- and tissue-levels aging phenotypes tend to appear as stereotypical changes. The mammary epithelium is a bilayer of two major phenotypically and functionally distinct cell lineages: luminal epithelial and myoepithelial cells. Mammary luminal epithelia exhibit substantial stereotypical changes with age that merit attention because these cells are the putative cells-of-origin for breast cancers. We hypothesize that effects from aging that impinge upon maintenance of lineage fidelity increase susceptibility to cancer initiation. We generated and analyzed transcriptomes from primary luminal epithelial and myoepithelial cells from younger <30 (y)ears old and older >55y women. In addition to age-dependent directional changes in gene expression, we observed increased transcriptional variance with age that contributed to genome-wide loss of lineage fidelity. Age-dependent variant responses were common to both lineages, whereas directional changes were almost exclusively detected in luminal epithelia and involved altered regulation of chromatin and genome organizers such as SATB1. Epithelial expression of gap junction protein GJB6 increased with age, and modulation of GJB6 expression in heterochronous co-cultures revealed that it provided a communication conduit from myoepithelial cells that drove directional change in luminal cells. Age-dependent luminal transcriptomes comprised a prominent signal that could be detected in bulk tissue during aging and transition into cancers. A machine learning classifier based on luminal-specific aging distinguished normal from cancer tissue and was highly predictive of breast cancer subtype. We speculate that luminal epithelia are the ultimate site of integration of the variant responses to aging in their surrounding tissue, and that their emergent phenotype both endows cells with the ability to become cancer-cells-of-origin and represents a biosensor that presages cancer susceptibility.