Notch controls the cell cycle to define leader versus follower identities during collective cell migration
Abstract
Coordination of cell proliferation and migration is fundamental for life, and its dysregulation has catastrophic consequences, such as cancer. How cell cycle progression affects migration, and vice versa, remains largely unknown. We address these questions by combining in silico modelling and in vivo experimentation in the zebrafish trunk neural crest (TNC). TNC migrate collectively, forming chains with a leader cell directing the movement of trailing followers. We show that the acquisition of migratory identity is autonomously controlled by Notch signalling in TNC. High Notch activity defines leaders, while low Notch determines followers. Moreover, cell cycle progression is required for TNC migration and is regulated by Notch. Cells with low Notch activity stay longer in G1 and become followers, while leaders with high Notch activity quickly undergo G1/S transition and remain in S-phase longer. In conclusion, TNC migratory identities are defined through the interaction of Notch signalling and cell cycle progression.
Editor's evaluation
Using a combination of in vivo and in silico approaches, the authors have demonstrated how cell-fate decisions are orchestrated at the level of leader vs. follower cells in collective cell migration of trunk neural crest cells. They highlight the role of Notch signaling and cell cycle progression, showing how these traits differ between the leader and follower cells. The findings are of wide interest, as collective cell migration is a fundamental process critical for embryonic development as well as invasion of various cancers.
https://doi.org/10.7554/eLife.73550.sa0Introduction
The harmonious coupling of cell proliferation with migration is fundamental for the normal growth and homeostasis of multicellular organisms. A prominent consequence of the dysregulation of these processes is cancer. Uncontrolled cell proliferation leads to primary tumours, and the acquisition of migratory capacities leads to the formation of secondary tumours, the most common cause of cancer deaths. Metastatic cells can migrate collectively, which endows them with more aggressive behaviours (Nagai et al., 2020). Collective cell migration refers to the movement of a group of cells that maintain contact and read guidance cues cooperatively (Rorth, 2009). This mechanism has been studied in several contexts, such as wound healing, angiogenesis, and neural crest (NC) migration. However, how cell proliferation impacts collective cell migration, and vice versa, remains largely unknown. The molecular signals that may couple these two fundamental processes remain equally unclear.
The NC is a mesenchymal cell population that arises early in development and migrates throughout the body, giving rise to a variety of cell types (neurons, glia, pigment cells, etc.). The NC’s stereotypical migratory behaviour (Gammill and Roffers-Agarwal, 2010) and similarity to metastatic cells (Maguire et al., 2015) make this cell type an ideal model to study the mechanisms of collective cell migration in vivo. Our previous work has shown that zebrafish trunk neural crest (TNC) migrate collectively forming single-file chains (Richardson et al., 2016). One cell at the front of the chain, the leader, is the only cell capable of instructing directionality to the group, while follower cells trail the leader. This division of roles into leaders and followers has been observed in other collectively migrating systems (Theveneau and Linker, 2017). Moreover, histopathological studies from cancer samples and cell lines show clear morphological and molecular differences between the invasive front, leaders, and the lagging cells, followers (Pandya et al., 2017). One outstanding question from these studies is what are the signals that determine leader versus follower migratory identities?
Notch signalling is a cell-cell communication pathway that directly translates receptor activation at the membrane into gene expression changes. Notch receptors are activated by membrane-bound ligands of the Delta/Serrate/Lag2 family. Upon ligand binding, Notch receptors are cleaved by γ-secretases releasing their intracellular domain (NICD). Subsequently, NICD translocates to the nucleus, binds the CBF1/Su(H)/Lag-1 complex, and initiates transcription (Bray, 2016). Among the direct Notch targets are members of the Hes gene family, which encode transcriptional repressors able to antagonise the expression of specific cell fate determinants and Notch ligands, generating a negative feedback loop in which cells with high Notch receptor activity downregulate the expression of Notch ligands, and cannot activate the pathway in their neighbours. Hence, adjacent cells interacting through the Notch pathway typically end up with either low or high levels of Notch activity and adopt distinct fates, a mechanism known as lateral inhibition (Lewis, 1998). Interestingly, Notch signalling has also been implicated in cell migration (Giniger, 1998; Leslie et al., 2007; Timmerman et al., 2004) and promotes invasive behaviours during cancer progression (Reichrath and Reichrath, 2012). Furthermore, lateral inhibition is implicated in the allocation of migratory identities during angiogenesis (Phng and Gerhardt, 2009), trachea formation in Drosophila (Caussinus et al., 2008), and in cell culture (Riahi et al., 2015). Whether Notch signalling plays a similar role in the context of mesenchymal cell migration is unknown. Notch signalling is required for NC induction (Cornell and Eisen, 2005), and its components and activity remain present in migrating NC (Liu et al., 2015; Rios et al., 2011). Nevertheless, the role of Notch during NC migration remains unclear. Cardiac NC are reported to develop normally under lack of Notch signalling (High et al., 2007). However, using different genetic tools, it has been shown that both gain and loss of Notch function led to the lack of NC derivatives (Mead and Yutzey, 2012). Moreover, in Xenopus the loss of Notch effectors leads to aberrant NC migration (Vega‐López et al., 2015).
The Notch pathway has not only been implicated in cell fate allocation, but it is also important for cell proliferation. Depending on the context, Notch can inhibit or promote cell cycle progression (Campos et al., 2002; Carlson et al., 2008; Devgan et al., 2005; Fang et al., 2017; Georgia et al., 2006; Mammucari et al., 2005; Nguyen et al., 2006; Nicoli et al., 2012; Noseda et al., 2004; Ohnuma et al., 1999; Park et al., 2005; Patel et al., 2016; Rangarajan et al., 2001; Riccio et al., 2008; Zalc et al., 2014). Indeed, Notch target genes include important cell cycle regulators such as CyclinD1, p21 and MYC (Campa et al., 2008; Guo et al., 2009; Joshi et al., 2009; Palomero et al., 2006; Ronchini and Capobianco, 2001).
Using a combination of in vivo and in silico approaches, we have established that differences in Notch activity between premigratory TNC select the leader cell. Cells with high levels of Notch signalling adopt a leader identity, while cells that lack Notch activity become followers. Our data show that a single progenitor cell in the premigratory area divides asymmetrically, giving rise to a large prospective leader and smaller follower cell. We propose that this original small asymmetry generates differences in Notch activity between TNC that are thereafter enhanced by cell-cell communication through Notch lateral inhibition. Differences in Notch activity in turn drive distinct cell cycle progression patterns and regulate the expression of phox2bb. Leader cells undergo the G1/S transition faster and remain in S-phase for longer than follower cells. Moreover, continuous progression through the cell cycle is required for TNC migration. Taken together, our results support a model in which the interaction between Notch and the cell cycle defines leader and follower migratory behaviours.
Results
Notch signalling is required for TNC migration
NC cells are induced at the border of the neural plate early during development. The prospective NC expresses Notch components, and Notch activity is required for NC induction (Cornell and Eisen, 2005). Our analysis reveals that Notch components remain expressed in NC after induction, suggesting that Notch signalling may also be involved in later aspects of NC development (Figure 1—figure supplement 1). Moreover, analysis of the Notch activity reporter line 12xNRE:egpf (Moro et al., 2013) shows that Notch signalling levels vary widely between premigratory TNC (Figure 1), suggesting that Notch may play a role after TNC induction. To explore the role of Notch in TNC development, we first aimed to define the stage at which NC induction becomes independent of Notch signalling. To this end, we treated embryos with the γ-secretase inhibitor DAPT (Richter et al., 2017) and assessed expression of NC marker. Our results showed that Notch inhibition impairs TNC induction up to 11 hours post-fertilization (hpf; Figure 2) and confirmed previous reports that induction of the cranial and vagal NC populations is independent of Notch signalling (Cornell and Eisen, 2000). Next, we analysed the effect of Notch inhibition at 12 hpf on the development of TNC derivatives. We found a reduction in all TNC derivatives (neurons, glia, and pigment cells; Figure 3A–F) upon Notch inhibition, suggesting that Notch activity is important in a process subsequent to induction, yet prior to differentiation. We next explored whether TNC migration is affected by Notch inhibition. Analysis of crestin expression showed a reduction in the number of TNC cell chains formed and in their ventral advance upon DAPT treatment (Figure 3G–J), which likely explains the lack of TNC derivatives at later stages. We then asked whether these results are due to a delay or a halt of migration. To this end, embryos were treated with DAPT from 12 hpf for 6–12 hr and processed for crestin expression. Decreased numbers of migratory chains were observed at all timepoints, but as embryos developed new chains were formed, indicating that the blockade of Notch signalling delays TNC migration (Figure 3K). Comparable results were obtained by inhibiting Notch genetically in embryos where the dominant-negative form of Suppressor of Hairless is under the control of a heat shock element (Latimer et al., 2005; hs:dnSu(H); Figure 3L). We reasoned that if Notch inhibition delays the onset of TNC migration, its overactivation might lead to TNC migrating earlier, leading to an increased number of chains. To test this, we induced NICD expression in all tissues by heat shock of hs:Gal4;UAS:NICD embryos (Scheer and Campos-Ortega, 1999). To our surprise, Notch gain of function (GOF) and loss of function (LOF) resulted in almost identical phenotypes, both showing a similar reduction of TNC chain numbers (Figure 3L). Taken together, these results show that precise regulation of Notch signalling levels is required for TNC migration.
In vivo Notch activity allocates TNC migratory identity
Interestingly, Notch signalling is required during collective migration to define distinct identities (Phng and Gerhardt, 2009; Caussinus et al., 2008; Riahi et al., 2015). To test whether Notch plays a similar role in TNC migration, we performed live-imaging analysis of TNC migration under lack (inhibition and LOF) or overactivation (GOF) of Notch signalling (Figure 4, Figure 4—videos 1 and 2). Our previous work defined a leader as the cell that retains the front position of the chain throughout migration, advancing faster and in a more directional manner than followers (Richardson et al., 2016). Under Notch inhibition (treatment with γ-secretase inhibitor Compound E; Richter et al., 2017), TNC remain motile with a single cell initiating the movement of the chain, but in contrast to control treatment (DMSO) the leader cell is unable to retain the front position and is overtaken by one or several followers (Figure 4A and C, Figure 5A and B, Figure 4—video 1). The overtaking follower cell, in turn, is not always able to retain the front position and can be overtaken by cells further behind in the chain. This loss of group coherence corresponds with a reduction in ventral advance, with most leader cells unable to move beyond the neural tube/notochord boundary (NT/not; Figure 4C, Figure 5A and C). This behaviour leads to an accumulation of cells at the NT/not, where some cells repolarise moving anterior or posteriorly and crossing the somite boundary and, in some cases, joining adjacent chains. Analysis of single-cell tracking showed that under Notch inhibition leader cells also have decreased speed and directionality (Figure 5D and E). Similar results were observed when Notch inhibition was achieved genetically by driving overexpression of dnSu(H) through heat shock in the entire embryo (not shown; hs:dnSu(H) line). Together, these results strongly suggest that upon lack of Notch signalling the TNC population is formed solely by follower cells that are unable to coordinate the movement of the group. Nevertheless, Notch signalling is important for the development of tissues surrounding TNC that act as a substrate for migration, raising the possibility that Notch signalling does not act cell-autonomously in TNC and instead the phenotypes observed are simply the consequence of somite and/or neural tube malformations. However, this appears unlikely as somite development (formation, patterning, and differentiation) and neuron formation are not affected by Notch inhibition at the axial level analysed (Figure 4—figure supplement 1). Next, we directly tested whether Notch signalling is autonomously required in TNC by inhibiting Notch activity exclusively in NC at the time of migration. To this end, we generated a new UAS:dnSu(H) line and crossed it with Sox10:Kalt4 fish (Alhashem et al., 2021). In the resultant embryos, all NC express Gal4 fused to the oestrogen receptor binding region (Gal4-ER) and are fluorescently labelled by nuclear-RFP. Under normal conditions, Gal4-ER is maintained inactive in the cytoplasm, whilst upon addition of tamoxifen, Gal4-ER is translocated to the nucleus activating transcription from the UAS:dnSu(H) transgene (Figure 4—figure supplement 2). We found that autonomous inhibition of Notch signalling in NC phenocopies the chemical inhibition. Leader cells are unable to retain the front position, being overtaken by followers, and ventral advance is reduced with cells accumulating at the NT/not boundary (Figure 4D, Figure 5A–C, Figure 4—video 2). Moreover, leader cells adopt followers’ migratory parameters, showing decreased speed and directionality (Figure 5D and E), confirming that Notch activity is autonomously required in TNC for identity allocation, and suggest that in the absence of Notch signalling a homogenous group of followers is established. In view of these results, we hypothesised that a homogeneous group of leaders would be formed upon Notch overactivation. Using a similar strategy, Notch overactivation was induced in the whole embryo (not shown, hs:Gal4;UAS:NICD; Scheer and Campos-Ortega, 1999), or exclusively in NC (Sox10:Kalt4;UAS:NICD), and migration was analysed by live imaging. Similar results were obtained in both experimental conditions: group coherence is lost, leader cells are overtaken by followers, and ventral advance is impaired (Figure 4F, Figure 5A–C, Figure 4—video 2). Interestingly, in Notch GOF conditions follower cells adopt leaders’ characteristics, moving with increased speed, but all cells in the chain follow less directional trajectories, which hinders the ventral advance of the group (Figure 5D and E), indicating that all cells in the chain migrate as leaders. Next, we tested whether the behavioural changes observed upon Notch alterations were mirrored by molecular changes by using the leader marker phox2bb. In control conditions, phox2bb transcripts are highly enriched in the leader cells from early stages of migration (Figure 6A, B, and G; Alhashem et al., 2022a). Consistent with expectations, upon Notch overactivation phox2bb is expressed by all the cells in the chain (Figure 6C, D, and G), while its expression is absent when Notch is inhibited (Figure 6E–G). These data show that Notch activity controls phox2bb expression and allocates TNC migratory identity.
In summary, our in vivo and molecular data show that Notch signalling is required autonomously in TNC for migratory identity allocation. TNC with high levels of Notch express phox2bb and become leaders, while cells with low Notch activity migrate as followers. Alterations of Notch signalling lead to a homogeneous TNC group with a single migratory identity that is unable to undergo collective migration. Taken together, these data suggest Notch lateral inhibition as the mechanism responsible for TNC migratory identity acquisition.
In silico modelling predicts that more than one leader is required for TNC migration
Our in vivo analysis shows that upon both Notch inhibition and overactivation TNC are unable to undergo collective migration due to lack of group coherence. On the other hand, our molecular analysis shows that upon Notch inhibition an all-followers group is established, while Notch overactivation leads to the formation of an all-leaders group. To gain a better understanding of these paradoxical results, we took an in silico approach, developing a discrete element model of TNC migration. Cells were simulated as 2D particles moving into a constrained space and endowed with intrinsic motility. Four variables control cell movement in the model: contact inhibition of locomotion (CIL) and co-attraction (co-A) define movement directionality and group cohesion, while volume exclusion regulates cell overlap, intuitively understood as cell size, while a noise element (zeta) was added to the cell’s trajectory (Figure 7A). A multi-objective scoring system, based on in vivo measurements, was developed to evaluate how close simulations with different underlying mechanisms matched chain behaviours. The scores were (1) chain cohesion, a maximum distance of 57 μm is allowed between adjacent cells; (2) single file migration for at least 80% of the simulation; (3) followers undergo rearrangements, while (4) leaders retain the front position, and (5) the chain should advance to the end of the migratory path (Figure 7B). Using this analysis and a parsimonious modelling approach, we attempted to match in vivo TNC migration with the simplest form of the model, only adding complexity incrementally in an effort to find the minimal set of predicted mechanisms required. We first simulated chains composed of homogeneous cells and systematically covaried all parameters. We found no parameter combination able to match all scores, confirming our previous findings that cell heterogeneity is required for TNC migration (Figure 7C; Richardson et al., 2016). Evidence from other systems (Astin et al., 2010; Bentley et al., 2014; Parkinson and Edwards, 1978; Theveneau and Mayor, 2013) led us to hypothesise that differences in the CIL response between cells may be at play. Thus, we simulated chains in which only cells of different identities present CIL (Diff CIL; Figure 7A). These simulations match several scores, but chains are unable to reach the end of the migratory path (Figure 7C, Figure 4—video 3). Next, we varied Diff CIL intensity, co-A, and cell size (volume exclusion) for leader cells. Interestingly, the model is only able to recapitulate control conditions when the difference between leaders and follower is maximal for all variables. Nevertheless, it is unable to recapitulate Notch GOF and LOF phenotypes (Figure 7C). Our previous results show that differences in Notch signalling establish migratory identities, suggesting that lateral inhibition may be the mechanism at play. To explore whether different outcomes of lateral inhibition may allow the model to simulate Notch altered conditions (GOF and LOF), different ratios of leader/follower cells were simulated. We first tested a 1:1 ratio, surprisingly this chain architecture over-migrates, moving beyond the end of the pathway (Figure 7C, Figure 4—video 3). Interestingly, we found that several parameter combinations from the 1:2 and 1:3 leader/follower ratios were able to recapitulate in vivo control condition, as well as the loss of group coherence and ventral advance observed in Notch GOF (all leader simulation) and LOF (all follower simulation; Figure 4B, E, and G, Figure 5, Figure 4—video 3). In these simulations, the six parameter combinations that match all in vivo scores had followers at the low setting, while leaders’ CIL intensity took medium or high values, cell size took medium or low values, and co-attraction took all levels. Nevertheless, all these parameter combinations endow the leader with enhanced migratory behaviour.
Next, we used a linear discriminant analysis (LDA) to study which of the model parameters bear most weight in the definition of leader and follower identity. LDA is a dimensionality reduction method that projects the data onto a lower dimensional space minimizing the variation within classes (e.g. between leaders) and maximizing the variation between classes (leaders versus followers), allowing the hierarchical ordering of the factors that best explain the class separation. First, we used the in vivo data to determine whether leaders and followers were properly separated by LDA. A visual inspection of the data makes clear that LDA works well to classify migratory identities (Figure 7D). Moreover, the LDA shows that ventral distance is the most important variable separating leaders from followers, with speed and directionality playing a less dominant role (Figure 7E). Next, we used this method to assess the importance of each of the model parameters. CIL intensity appears to be the parameter that most differ between leader and follower cells, while heterogeneity in the other parameters is not essential (Figure 7F). Taken together, the in silico data confirms our previous conclusion that TNC chains are a heterogeneous group. Remarkably, it also predicts CIL intensity to be the most important distinction between leaders and followers. Finally, the model anticipates that TNC chains are formed of leaders and followers in a 1:2 or 1:3 ratio.
Leader cells arise from the asymmetric division of a progenitor cell
Cell size is a prominent characteristic distinguishing leader from follower cells. Leaders are almost twice as big as followers during migration and this difference is evident before migration initiation (Richardson et al., 2016), suggesting that size disparity arises at birth or shortly thereafter. Interestingly, differential cell size emerged as an important parameter in our in silico analysis, contributing to more realistic leader/follower coordination behaviours. To understand the origin of these size differences, we investigated whether leader and follower cells share a common progenitor, and at which point differences in size become apparent. To this end, we imaged FoxD3:mCherry;H2aFVA:H2a-GFP embryos. The FoxD3:mCherry reporter (Hochgreb-Hägele and Bronner, 2013; Lukoseviciute et al., 2018) labels NC from early stages and allows us to define TNC identity at later stages by their migratory position. Moreover, the nuclear marker H2aFVA:H2a-GFP (Pauls et al., 2001) was used to track single cells and their divisions. Tracking analysis shows that the asymmetric division of a single progenitor cell in each body segment gives rise to a larger cell that becomes a leader (102 ± 20 µm2), and a smaller sibling that migrates as follower (72 ± 9 µm2; Figure 8A and B, Figure 8—video 1). In contrast, all other progenitors divide symmetrically, giving rise to two follower cells (87 ± 27 µm2; Figure 8C and D). We also noticed that the leader progenitors’ divisions are spatially restricted to the anterior quarter of the premigratory area in each segment, while the followers’ progenitor divisions take place across the premigratory area (Figure 8E).
We then reasoned that leader cells, being bigger, may undergo the next division in a shorter time span than follower cells, and in consequence, mitotic figures would be observed at different, but consistent, positions in their trajectory. Indeed, we found two different patterns of divisions in respect to migration: (1) cells that first Divide and then Migrate (D→M) or (2) cells that first Migrate and then Divide (M→D; Figure 8F and G, Figure 8—video 2). Interestingly, we found that the patterns of cell division correlate with cell identity. Most leader cells divide during migration (M→D: 86%), while the bulk of follower cells divide before migration initiation (D→M: 90%, Figure 8H). These patterns result in leader and follower cells dividing at distinct positions, 74% of leaders divide at the NT/not boundary (65.3 ± 9.6 µm), while 85% of followers divide mostly within the premigratory area or in the dorsal-most region of the somite (42 ± 12.4 µm; Figure 8I). Together, these results show that leader cells arise from the asymmetric division of a progenitor. Thereafter, leader and follower cells show distinct locations and patterns of division, suggesting that leaders and followers progress asynchronously through the cell cycle, which may influence their migratory behaviour.
Cell cycle progression is required for TNC migration
To test the role of cell cycle progression in TNC migration directly, we used inhibitory drugs. The S-phase inhibitor aphidicolin blocks over 94.7% ± 4.5% of mitotic figures after 3 hr of treatment, while the G2/M inhibitor genistein prevents 90% ± 10% of divisions within 6 hr, but neither treatment affects NC induction (Figure 9—figure supplement 1). Inhibition of cell cycle progression by either of the treatments resulted in reduced numbers of migratory chains and decreased ventral advance (control 19 ± 2, genistein 10 ± 3, aphidicolin 6 ± 2 chains; Figure 9A–H). This result was not due to the loss of cell motility as premigratory TNC cells actively extend protrusions and move along the anteroposterior axis but are unable to migrate ventrally (Figure 9—video 1). Importantly, these effects were not a consequence of cell death or the permanent impairment of motility as TNC reinitiate migration and form new chains upon drug withdrawal (Figure 9G and H), showing that active cell cycle progression is required for migration. Next, we directly analysed TNC cell cycle progression in vivo. To this end, we imaged Sox10:FUCCI embryos (Rajan and Gallik, 2018), in which TNC nuclei are RFP-labelled during G1 and GFP-labelled during S and G2. Tracking analysis shows differential cell cycle progression, with most leader cells initiating migration in S-phase (79%), while followers start movement during G1 (77%; Figure 9I and J, Figure 9—video 2). These results show that cell cycle progression is required for migration and that leader and follower cells initiate movement at different points of the cell cycle, suggesting an intimate connection between cell growth and movement.
Leader and follower cells progress through the cell cycle at different rates
Next, we studied TNC cell cycle progression in detail. First, we asked whether leaders and followers differ in the total length of their cell cycle. Measurements of the time span between two consecutive mitoses showed no significant differences in the total length of the cell cycle between leaders and followers (13.6 ± 1.2 and 13.3 ± 1.4 hr, respectively; Figure 10B). Next, we examined the length of each phase of the cell cycle by imaging the characteristic nuclear labelling pattern of the PCNA-GFP fusion protein (Leung et al., 2011). Sox10:Kalt4 embryos, in which all NC can be recognised by nuclear RFP expression, were injected with PCNA-GFP mRNA and live imaging was performed. PCNA-GFP shows uniform nuclear GFP labelling during G1, intense fluorescent nuclear puncta characterise the S-phase, these puncta dissipate during G2 restoring homogeneous nuclear fluorescence, at the onset of mitosis PCNA is degraded and TNC are recognised solely by nuclear RFP (Figure 10A, Figure 10—video 1). In these embryos, leader cells initiate migration during S-phase and followers in G1, confirming our FUCCI results and establishing that PCNA overexpression does not introduce artefacts to cell cycle progression (Figure 10—figure supplement 1). Using this tool, we measured the length of the cell cycle phases in TNC. We found striking differences in the time spent in G1- and S-phase between leader and follower cells. Leaders present a short G1 (3.2 ± 0.6 hr) but remain for twice as long in S-phase (8.7 ± 1.3 hr). Followers, on the other hand, present the opposite distribution, remaining for twice as long in G1 (7.4 ± 2.7 hr) than in S-phase (4.6 ± 2.8 hr; Figure 10C and D). No significant differences were observed in the length of G2 (leaders 1.6 ± 0.4 hr; followers 1.5 ± 0.3 hr) or M (leaders 0.6 ± 0.1 hr; followers 0.5 ± 0.1 hr). These data show that leader and follower cells present marked differences in the length of G1- and S-phase, suggesting that cell cycle progression may regulate their migratory behaviour.
Notch signalling regulates TNC cell cycle progression
Our data show that Notch signalling allocates leader and follower identities, cell cycle progression is necessary for TNC migration, and leader and follower cells progress through the cell cycle at different rates. Does Notch signalling regulate cell cycle progression, thus differentiating leader from follower cells? To investigate this question, we measured the total length of the cell cycle and the length of each phase under control and Notch-inhibited conditions. Neither the total cell cycle length (Figure 11A) nor the number of TNC (Figure 9—figure supplement 1) were affected by alterations of Notch signalling. Remarkably, we found significant differences in the length of G1- and S-phase upon Notch inhibition. Leader cells lose their characteristic cell cycle progression pattern and behave as followers, with a long G1 and a short S-phase (Figure 11B). Furthermore, Notch inhibition abolishes the size difference between migratory leader and follower cells, with all cells presenting the average follower’s area (Figure 11C and D). These data show that Notch activity defines TNC migratory identity by regulating cell cycle progression, cells with low Notch activity remain for longer in G1 behaving as followers. Interestingly, we noticed that Notch inhibition also changes the cell cycle behaviour of the followers’ population. While the followers’ average length of cell cycle phases is not altered, the dispersion of this population is significantly reduced, with standard deviations cut almost by half (from 2.7 hr to 1.42 hr for G1 and from 2.8 hr to 1.38 hr for S; Figure 11B). This prompted us to analyse the frequency distribution of cell cycle phases length. In control conditions, leader cells show a normal distribution with a single peak for G1- and S-phase, as expected for a homogeneous population. Followers, on the other hand, present a bimodal distribution, with the smaller peak coinciding with that of leader cells, and accounting for 26% of followers in G1- and 31% in S-phase (Figure 11E and F). Strikingly, these results fulfil the predictions of our in silico model that best recapitulates TNC migration when chains are composed of leaders and followers in a 1:2 or 1:3 ratio. Furthermore, upon Notch inhibition the bimodal distribution of the follower population is lost, with all cells grouped at the major mean (Figure 11G and H). Consistent with these data, closer analysis (at higher magnification) of normal phox2bb expression shows increased expression in followers at position 3 in addition to that in leaders (Figure 11I–N). Taken all together, our data demonstrate that the levels of Notch activity in TNC allocate migratory identity by controlling cell cycle progression and that migratory chains are formed of one leader cell for every three followers.
Discussion
Collective migration plays an important role in embryogenesis, wound healing, and cancer. The acquisition of specific migratory identities has proven fundamental to angiogenesis, trachea development in Drosophila, and cancer metastasis. TNC migrate collectively, forming chains with a leader cell at the front of the group that direct the migration, while follower cells form the body of the chain that trails the leaders. TNC leader and follower identities are established before migration initiation and remain fixed thereafter (Richardson et al., 2016). Herein, we have addressed the mechanism that establishes leader and follower identities and can propose the following model (Figure 12): (A) premigratory TNC progenitors arise at the dorsal part of the neural tube. The leader’s progenitor divides asymmetrically, giving rise to a large prospective leader cell and a small sibling that migrates as a follower. Other progenitors divide symmetrically, giving rise to follower cells. (B) Interactions via Notch signalling results in the prospective leader cell accumulating higher levels of Notch activity, which induces phox2bb expression. (C) The combination of high Notch activity and a larger cell size prompts the prospective leader cell to rapidly undergo the G1/S transition, entering S-phase and initiating migration earlier than its follower siblings, which are smaller and initiate migration whilst in G1. (D) Premigratory cells that have not been in contact with the prospective leader cell or that have lost contact with it due to its ventral advance maintain communication with surrounding premigratory TNC through Notch and undergo a new round of leader cell selection. This working model of TNC migration is supported by both our experimental data and our in silico modelling, and provides a useful conceptual framework for future studies to build upon.
Notch signalling is a seemingly simple pathway that directly transduces receptor activation into changes in gene expression. Nevertheless, its outcomes in terms of cellular patterning are very diverse, from the generation of gene expression boundaries to temporal oscillations, or from the induction of similar fates in neighbouring cells to forcing adjacent cells into alternative fates. The latter function, known as lateral inhibition, is characterised by an intercellular negative feedback loop regulating the expression of Notch ligands. The activation of the Notch receptor in a ‘signal-receiving’ cell leads to the downregulation of Notch ligands expression, making it less able to act as a ‘signal-sending’ cell. The signature 2D patterning outcome of lateral inhibition is a mosaic of signal-sending cells with low Notch activity, surrounded by signal-receiving cells with high Notch levels. This is the case during the selection of sensory organ precursor cells in the epidermis of Drosophila (Lewis, 1998) or the formation of the mosaic of hair cells and supporting cells in the sensory organs of the inner ear (Daudet and Żak, 2020). In general, however, lateral inhibition operates among cells subjected to extensive rearrangements and its patterning outcome is not a salt-and-pepper mosaic of cells (Bocci et al., 2020). For example, during angiogenesis, cells with low Notch signalling become tip or leaders, while cells with high Notch activity differentiate as stalk or followers (Phng and Gerhardt, 2009). In this context, leaders are interspersed with various numbers of followers. Several models have been proposed to explain how signal-sending (leader/tip) cells can exert a long-lasting or long-range inhibition on signal-receiving (follower/stalk) cells. These take into account the modulation of Notch signalling that arise from heterogeneity in Notch receptor levels, tension, Notch-regulators, and interaction with other pathways (Bentley and Chakravartula, 2017; Hadjivasiliou et al., 2019; Koon et al., 2018; Kur et al., 2016; Venkatraman et al., 2016). Our data show that TNC deviate from the classical mosaic pattern, forming chains with one leader every two or three followers. Further studies will be required to define whether the aforementioned mechanisms are responsible for this architecture.
In the case of the TNC, however, the most striking divergence from the classic lateral inhibition model (or indeed angiogenesis) is the fact that the leader cell identity is associated with higher intrinsic Notch activity. In other words, there are more signal-sending cells than signal-receiving cells. This apparent inversion in the ratio of the cell types produced is surprising. Explanation of this conundrum may arise from the fact that Notch lateral inhibition, dynamics and outcomes, can be modulated by ‘cis-inhibition’, a process whereby Notch ligands cell-autonomously interfere with the activation of Notch receptors (Bray, 2016; del Álamo et al., 2011). Computational models show that an increase in the strength of cis-inhibition can result in the inversion of the salt-and-pepper pattern (signal-sending to signal-receiving cells ratio), with the production of one cell with high Notch activity for every three cells with low Notch levels (Formosa-Jordan and Ibañes, 2014), a scenario that is congruent with the leader/follower ratio we observe in TNC. The detailed dynamics of lateral inhibition and whether cis-inhibition is at work in TNC remain to be investigated and will require direct visualisation at the single-cell level of Notch activity in live embryos.
Our data show that active progression through the cell cycle is required for TNC migration. This is consistent with studies in chicken embryos, showing that progression through G1/S is required for TNC delamination, and that NC continue cycling as they migrate (Burstyn-Cohen and Kalcheim, 2002; Théveneau et al., 2007). Our data extend these findings by showing that leader and follower cells progress through the cell cycle at different rates. Leader cells, which are larger and more motile, initiate migration in S-phase and spend twice as long in this phase as followers. It is possible that these differences arise from the fact that leaders are larger than followers. It has been shown that the timing of G1/S transition depends on cell size and the dilution of the nuclear retinoblastoma protein (Zatulovskiy and Skotheim, 2020). Due to the larger volume of their cytoplasm, leader cells could be primed for a rapid G1/S phase transition. The initiation of S-phase may in turn enhance leaders’ migratory characteristics through the interaction of cyclins and cyclin/CDK inhibitors (CDKI) with small GTPases. Cyclins B and D have been shown to phosphorylate cytoskeleton regulators, resulting in increased cell migration and tumour invasion (Blethrow et al., 2008; Chen et al., 2020; Chi et al., 2008; Hirota et al., 2000; Li et al., 2006; Manes et al., 2003; Song et al., 2008; Zhong et al., 2010). Furthermore, Rac1 activity, which is required for migration, oscillates during the cell cycle being highest at S-phase when cells are most invasive (Kagawa et al., 2013; Walmod et al., 2004). CDKIs, on the other hand, interact with RhoA- and ROCK-enhancing motility (Bendris et al., 2015; Creff and Besson, 2020; Yoon et al., 2012). Interestingly, enhanced motility increases actin branching, which in turn can accelerate the G1/S transition (Molinie et al., 2019). These factors could therefore generate a positive feedback loop in which slightly larger leader cells are prone to undergo the G1/S transition, in turn the activation of S-phase cyclins and CDKIs may enhance motility, reinforcing S-phase initiation.
Our data also show that TNC cell cycle progression is under the control of Notch signalling. Upon Notch inhibition, all TNC present cell cycle phase lengths typical of follower cells. Notch has been shown to regulate cell cycle in a context-dependent manner. Depending on the cell type, Notch can regulate cell cycle through the transcriptional induction of cyclins A and D, and the inhibition of CDKIs (Campa et al., 2008; Dabral et al., 2016; Ridgway et al., 2006; Rizzo et al., 2008; Rowan et al., 2008). Conversely, cell cycle progression can impact on Notch signalling. Notch activity is enhanced at the G1/S transition, while cells become refractory to Notch during G2/M (Ambros, 1999; Carrieri et al., 2019; Hunter et al., 2016; Nusser-Stein et al., 2012). Hence, the combination of large volumes and higher Notch activity levels could act synergistically to promote leaders’ G1/S transition.
In this study, we have uncovered new functional interactions between Notch signalling, cell cycle dynamics, and the migratory behaviour of leader and follower cells in the TNC. These complex and intricate interactions, which remain to be fully characterised at a molecular level, could apply to other cell types exhibiting collective migration. For example, studies in cancer cell lines have shown that activation or inhibition of Notch signalling hinders migration, similar to what we observe in TNC (Konen et al., 2017), while the maintenance of collective migration depends on the regulation of cell proliferation during angiogenesis (Costa et al., 2016). In view of our work, it is important to revisit the assumption that migratory phenotypes are in conflict with cell cycle progression (Kohrman and Matus, 2017) and consider the possible implication for cancer therapies.
Materials and methods
Resource availability
Request a detailed protocolFurther information and requests for resources and reagents should be directed to and will be fulfilled by the lead contact, Claudia Linker (claudia.linker@kcl.ac.uk).
Materials availability
Request a detailed protocolNewly generated materials from this study are available by request from the lead contact, Claudia Linker (claudia.linker@kcl.ac.uk), except for computational tools to be requested from Katie Bentley (katie.bentley@crick.ac.uk).
Data and code availability
Request a detailed protocolThe model code is accessible at https://github.com/Bentley-Cellular-Adaptive-Behaviour-Lab/NeuralCrestCpp, (copy archived at swh:1:rev:cdf63f3786390fb7905092717456cc69f5657ddc; Alhashem, 2022b). The code used to perform the LDA is accessible in the supplementary files. All numerical data used in the figures are accessible in the supplementary data source file.
Zebrafish lines and injections
Request a detailed protocolZebrafish were maintained in accordance with UK Home Office regulations UK Animals (Scientific Procedures) Act 1986, amended in 2013 under project licence P70880F4C. Embryos were obtained from the following strains: wild type, AB strain; Sox10:mG, Tg(–4.9sox10: Hsa.HIST1H2BJ-mCherry-2A-GLYPI-EGFP); Sox10:Fucci, Tg(–4.9sox10:mAGFP-gmnn-2A-mCherry-cdt1); hs:dnSu(H), vu21Tg (hsp70l:XdnSu(H)-myc); hs:Gal4, kca4Tg Tg(hsp70l:Gal4)1.5kca4 (1); UAS:NICD, Tg(UAS:myc-Notch1a-intra)kca3; Sox10:Kalt4, Tg(–4.9sox10: Hsa.HIST1H2BJ-mCherry-2A-Kalt4ER); UAS:dnSu(H), Tg(UAS:dnSu(H)-myc); Tg(h2afva:GFP)kca13; 12XNRE:egfp. Embryos were selected based on anatomical/developmental good health and the expression of fluorescent reporters when appropriate, split randomly between experimental groups and maintained at 28.5°C in E3 medium. Genotyping was performed by PCR of single embryos after imaging when required (UAS:NICD; UAS:dnSu(H); hs:dnSu(H)). Injections were carried at 1–4-cell stage with 30 pg of PCNA-GFP mRNA in a volume of 1 nl. mRNA was synthesised from pCS2 + PCNA GFP plasmid, kindly provided by C. Norden (IGC, Portugal), linearised with NotI and transcribed with the SP6 mMessage Machine Kit (Thermo Fisher Scientific, Cat# AM1340).
Live imaging and tracking
Request a detailed protocolImaging and analysis were carried out as in Alhashem et al., 2021. In short, embryos were mounted in 1% agarose/E3 medium plus 40 µM Tricaine. Segments 6–12 were imaged in lateral views every 5′ from 16 hpf for 16–18 hr in an upright PerkinElmer Ultraview Vox system using a ×40 water immersion objective. 70 μm z-stacks with 2 μm z-steps were obtained. Image stacks were corrected using Correct 3D Drift Fiji and single-cell tracking performed with View5D Fiji plugin. Tracks were displayed using the MTrackJ and Manual Tracking Fiji plugins. Cell area measurements were done in Fiji using the freehand selection tool to draw around cell membranes in 3D stacks using the orientation that best recapitulated the cell morphology (as in Richardson et al., 2016). Cell speed measurements were calculated from 3D tracks using the following formula: ((SQRT((X1-X2)^2+(Y1-Y2)^2+(Z1-Z2)^2))/T)*60, where X, Y, and Z are the physical coordinates and T is the time step between time-lapse frames. Ventral distances were measured in a straight line from dorsal edge of the embryo to the cell position at the end of the movie. Cell directionality measurements were calculated using a previously published Excel macro (Gorelik and Gautreau, 2014). Total duration of the cell cycle was measured between two mitotic events. Cell cycle phase duration were measured using the characteristic nuclear pattern of PCNA-GFP, in movies where only TNC (expressing RFP and GFP) are shown using this custom Fiji macro:
macro "Segment Nuclei [s]" {
title = getTitle();
run("Split Channels");
selectWindow("C1-" + title); //select window with C1 in its name, nuclei should be C1
getDimensions(width, height, channelCount, slices, frames);
run("Subtract Background...", "rolling = 200 sliding stack");
setAutoThreshold("Default dark");
run("Threshold...");
setThreshold(5, 255); //change as appropriate for your cells
setOption("BlackBackground", false);
run("Convert to Mask", "method = Default background = Dark");
run("Close");
run("Fill Holes", "stack");
run("Despeckle", "stack");
run("Dilate", "stack");
run("Dilate", "stack");
//now go over every frame and slice
for(frame = 1; frame ≤ frames; frame++){
for(slice = 1; slice ≤ slices; slice++){
selectWindow("C1-" + title);
setSlice(slice);
Stack.setFrame(frame);
run("Create Selection");
selectWindow("C2-" + title);
setSlice(slice);
Stack.setFrame(frame);
run("Restore Selection");
setBackgroundColor(0, 0, 0);
run("Clear Outside", "slice");
In situ hybridisation, immunostaining, and sectioning
Request a detailed protocolThe whole-mount in situ hybridisation protocol was adapted from https://wiki.zfin.org/display/prot/Whole-Mount+In+Situ+Hybridization. In short, embryos were fixed overnight (O/N) in 4% paraformaldehyde (PFA), dehydrated in 100% methanol, then rehydrated, digested with proteinase K for different times depending on the stage and pre-hybridised for 2 hr at 65°C. Riboprobes were added, and embryos incubated at 65°C O/N. Probes were removed and embryos washed and equilibrated to PBS. Embryos were incubated in blocking solution for 2 hr and in anti-dig antibody O/N (Sigma-Aldrich, Cat# 11093274910), washed 5 × 30′ and NBT/BCIP colour reaction performed. Riboprobes for notch1a, dlb (deltaB), dld (deltaD), her4, cb1045 were kindly provided by J. Lewis (CRUK); crestin, mbp, bdh, myoD by S. Wilson (UCL, UK). After the in situ colour development, embryos were processed for sections, washed 5 × 10′ with PBS, embedded in OCT, frozen by dipping the blocks in dry ice-cold 70% ETOH, and sectioned to 12–15 μm using a cryostat. Sections were thawed at room temperature (RT), incubated with blocking solution for 30′ (10% goat serum, 2% BSA, 0.5% Triton, 10 mM sodium azide in PBS) and in anti-GFP antibody ON at 4°C (Merck Millipore, Cat# 06-896). Sections were washed with PBST 5 × 5′ (0.5% Triton-PBS) and incubated with secondary antibody for 2 hr at RT, mounted in ProLong Gold Antifade Mountant (Molecular Probes, Cat# P10144) and imaged. Whole-mount antibody staining was performed in embryos fixed for 2 hr in 4% PFA, washed 4 × 10′, incubated in blocking solution for 2 hr and in primary antibodies O/N at 4°C (anti-myc, Cell Signaling, Cat# 2276S; F59 and Znp1, Developmental Studies Hybridoma Bank; acetylated tubulin, Sigma-Aldrich, Cat# MABT868). Embryos were washed 5 × 30′, incubated in secondary antibodies O/N at 4°C, washed 6 × 30′, and mounted in 1% agarose for imaging. Imaging of sectioned and whole-mount antibody-stained samples was performed in PerkinElmer Ultraview Vox system.
RNAScope (RNAscope Fluorescent Multiplex Reagent Kit, Cat# 320850) experiments were performed as in Alhashem et al., 2022a. In short, embryos were fixed with 4% PFA O/N at 4°C and dehydrated in 100% methanol and stored at –20°C until processing. All methanol was removed, and embryos were air dried at RT for 30′, permeabilised with Proteinase Plus for 10′ at RT (provided in kit), washed with PBS-Tween 0.01%, and incubated with probes for egfp and sox10 or phox2bb at 1:100 dilution at 60°C O/N. Probes were recovered, embryos washed three times with SSCT 0.2× for 15′. We followed manufacturer’s instructions for amplification steps AMP 1–3 and HRP C1–C4. Opal dyes 520, 570, and 650 (Akoya Biosciences, Cat# FP1487001KT, Cat# FP1488001KT, and Cat# FP1496001KT) were added at 1:3000 dilution followed by HRP blocker. Washes in between steps were performed with SSCT 0.2× for 10′ twice. Primary a-GFP-chicken (1:750) and a-RFP-rabbit (1:750; TFS, Cat# A10262, and MBL, Cat# PM005) antibodies diluted in blocking solution (PBS-Tween 0.1%, goat serum 5%, DMSO 1%) were added and incubated O/N at 4°C. Samples were washed three times in PBS-Tween 0.1% for 1 hr and then incubated in secondary antibodies, a-chicken-Alexa Fluor488 and a-rabbit-Alexa Fluor546 (TFS, Cat# A11039 and Cat# A11010) both in a 1:1000 dilution in blocking solution, for 3 hr at RT. Samples were washed six times with PBS-Tween 0.1% for 30′. For counterstaining, DAPI was added (1:1000) in the third wash (Roche, Cat# 10236276001, 2 mg/ml). Embryos were cleared in 50% glycerol/PBS an mounted in glass-bottom Petri dishes and imaged using Zeiss Laser Scanner Confocal Microscope 880 (405, 488, 514, 561, and 633 lasers).
Drug treatments and gene expression induction
Request a detailed protocolEmbryos were treated by adding cell cycle inhibitors to the media from 11 hpf and incubated for 3–12 hr at 28.5°C. 20 μM hydroxyurea (Sigma-Aldrich, Cat# H8627), 300 μM aphidicolin (Sigma-Aldrich, Cat# A0781), 100 μM genistein (Calbiochem, Cat# 345834), teniposide (Sigma-Aldrich, Cat# SML0609), or 1% DMSO as control (Sigma-Aldrich, Cat# D8418). Notch signalling was inhibited at 11 hpf by adding 100 μM DAPT (Sigma-Aldrich, Cat# D5942-25MG) or 50 μM of Compound E (Abcam, Cat# ab142164). The latter reagent was used to perform live imaging, which is difficult to do with DAPT as it generates an interfering precipitate. 1% DMSO was added as control. Gene expression was induced by addition of 2.5 μM of tamoxifen (Sigma-Aldrich, Cat#H7904) to the media at 11 hpf of Sox10:Kalt4 embryos, or by heat shock at 11 hpf in hs:Gal4 and hs:dnSu(H) embryos by changing the media to 39°C E3, followed by 1 hr incubation at this temperature, thereafter embryos were grown at 28.5°C to the desired stage.
Generation of UAS:dnSu(H) transgenic line
Request a detailed protocolUsing the Infusion cloning system (Takara), the following construct was inserted into the Ac/Ds vector (Chong-Morrison et al., 2018): 5xUAS sequence (Tol2Kit, http://tol2kit.genetics.utah.edu/index.php/Main_Page) flanked at the 3′ and 5′ ends by rabbit β-globin intron sequence. At the 3′ end, GFP followed by SV40polyA sequence was cloned to generate the Ac/Ds dUAS:GFP vector. The cmlc2:egfp transgenesis marker (Tol2Kit) was cloned after GFP in the contralateral strand to prevent interaction between the UAS and the cmnl sequences. The Xenopus dnSu(H)-myc sequence (Latimer et al., 2005) was cloned into the Ac/Ds dUAS:GFP vector at the 5′ end of the 5xUAS sequence, followed by the SV40polyA sequence (Figure 4—figure supplement 2). Transgenesis was obtained by injecting Sox10:Kalt4 embryos with 1 nl containing 50 pg of DNA plus 30 pg of Ac transposase mRNA at 1-cell stage. Embryos carrying the transgene were selected by heart GFP expression at 24 hpf. Upon Gal4ER activation by tamoxifen, dnSu(H)-myc protein was readily detected with anti-Myc antibody (Figure 4—figure supplement 2). GFP fluorescence driven by UAS was never observed.
Statistical analysis
Request a detailed protocolAll graphs and statistical analysis were carried out in GraphPad Prism 9. All numbers in the texts are mean ± standard deviation. Every sample was tested for normality using the d’Agostino–Pearson, followed by Shapiro–Wilk tests. Samples that passed both tests were compared using either unpaired two-tailed t-test or one-way ANOVA. Those without a normal distribution were compared through a Mann–Whitney U-test, Kruskal–Wallis test, or Brown–Forsythe and Welch’s ANOVA tests. For all analyses, p-values < 0.05 were deemed statistically significant, with ****p<0.0001, ***p<0.001, **p<0.01, and *p<0.05. Full statistical analysis of data in Figure 5 is presented in Supplementary file 1.
Computational model
Request a detailed protocolThe computational model used in this study is described in Appendix 1.
Standard LDA was carried out using the sklearn package in Python (see supplementary code files).
Appendix 1
Computer modelling methods
A minimal discrete element model of TNC migration was developed in which each cell is modelled as an infinitesimal particle moving in 2D space. A network of neighbours within the particle system is identified by a Delaunay triangulation (Appendix 1—figure 1).
Thus defined, the system exhibits Brownian dynamics as described by the over-damped Langevin equation (Equation 2), such that the velocity of each particle is proportional to the resultant force applied to it (), plus a stochastic component (ζ).
Components of the resultant force on each cell arise from cell-cell interactions, cell-boundary interactions, and cell-autonomous motion.
Tissue environment (boundary)
Cells move into permissive space between the neural tube/notochord and the somites. Boundary locations are specified before any simulation. The boundary is implemented as a region of space that applies strong repulsion to nearby cells (Appendix 1—figure 1). Any cell that moves within a cell radius of the boundary experiences a force given by the gradient of the same Morse potential used in cell-cell interactions, such that the repulsion of any cell from the boundary depends upon the cell volume exclusion and increases exponentially as the cell approaches the boundary (Appendix 1—figure 1).
The size and shape of the boundary represent a space for the pre-migratory cells at the top, a space in the middle where the notochord and neural tube meet (midline) and a vertical space where the chain can proceed downwards. The dimensions of the environment boundary were calibrated to in vivo measurements (Appendix 1—figure 1, showing micron-scale dimensions on the boundary).
The system is setup in a ‘T’ shape, which is interrupted in the middle by a space of horizontal mobility, because in vivo cells regularly move into this space. The wider region at the top represents the premigratory zone (PMZ) at the top of each migratory chain. Cells are able to filter in from the sides to mimic the continuous clustering of cells above migration chains.
Cell properties/behaviours
Contact inhibition and autonomous motion
Cells exhibit autonomous motion in a direction determined by their internal polarisation. This polarisation is influenced by interaction with the cell’s neighbours, such that the cell will try to move into empty space. We introduce contact inhibition into the model as a term in the Langevin equation (Equation 2), with magnitude determined by a user-defined parameter. The direction of autonomous magnitude for a given cell is found by identifying all adjacent nearest neighbours surrounding the cell, calculating the angle subtended by each adjacent pair, and bisecting the largest such angle (Appendix 1—figure 1). The magnitude of this autonomous velocity component is proportional to the user-defined parameter (aMag) and the square of the maximum subtended angle, representing the combined effect of greater polarisation and more free space to move into. Any cell that moves beyond a threshold distance from its nearest neighbour will stop autonomous motion, modelling the loss of polarisation when losing contact with neighbouring cells.
Cell volume exclusion
Cells exhibit volume exclusion (two cells repel from one another if they get closer than an equilibrium distance). This simply models how two cells cannot occupy the same space at the same time. The extent to which volume exclusion is exhibited can be thought of as the level of cell stiffness. Low k means cells are squishier. This is modelled using the k term in the Morse potential calculation (Equation 1).
Co-attraction (co-A)
When cells drift more than the equilibrium distance apart, they are drawn back towards their neighbours with a force calculated by the Morse potential curve (Equation 1).
Migratory identity
Leader and follower migratory identities were allocated to cells according to the order in which these enter the chain. That is, the first cell becomes leader then the next X many cells become follower cells before the cell after that becomes leader. A sensitivity analysis on leader cells frequency was performed by spacing parameter S.
Simulation procedure
The simulation follows the process steps listed below Appendix 1—figure 2 and was simulated on CAMP – the Francis Crick Institute’s Linux-based high-performance computing system. Parameter combination/experimental condition pairs were run 100 times in parallel across 10 nodes.
Initialise cells within PMZ
Loop until maximum time has been reached:
2.1 Perform Delaunay triangulation on cells
2.2 Identify neighbours
2.3 Calculate forces between nearest neighbours
2.4 Apply boundary forces to all cells
2.5 Calculate autonomous motion velocties as determined by contact inhibition
2.6 Calculate noise
2.7 Update system
2.8 Apply experimental conditions for the next time step
Model pseudocode overview.
A predefined number of cells is initialised in the premigratory zone (PMZ). Thereafter, the system enters a loop for every time step up to tmax. In this loop, forces are linearly summed to obtain each cell’s velocity vector for that time step: (1) a Delaunay triangulation is performed on cells. (2) Each cell’s nearest neighbours are identified. (3) Local forces between cells are calculated according to the Morse potential (Appendix 1—figure 1). (4) Boundary forces are applied to each cell. (5) Autonomous motion and contact inhibition are calculated for each cell. (6) Gaussian noise is added to each cell’s velocity vector. (7) The system’s clock is updated, as well as each cell’s position. (8) Experimental conditions are applied for the next time step (e.g. giving certain cells leader qualities).
Parameterisation
Time was calibrated as follows: in vivo control cells tend to migrate to approximately 120 µm from dorsal midline on average (Figure 5C). The total time of migration is on average 11.64 hr long (~700 min). In 2000 time steps, the control case (with differential CIL, heterogenous migratory identities and S = ‘1:3’) also migrates to approximately 120 µm. We gathered data every 20 time steps, which means in our simulation movies there are 100 frames (i.e. 1 frame = 7 min).
Where possible, parameters were calibrated to values measured in vivo (Appendix 1—table 1). Model-specific parameters unable to be linked directly to in vivo values were set to values that produced realistic bounds of behaviour.
Sensitivity analysis
In the grid search calibration approach, we fixed follower cells’ properties to be at their low levels. Next, we looked at how changes to leader cell physical properties affected ventral distance (Appendix 1—figure 2). This shows a strong effect in leader aMag, whereby low leader aMag resulted in cells not migrating much beyond 100 µm no matter the level of co-attraction or cell volume exclusion. aMag had to be varied across a wider range to see a clear effect. Through this, aMag has a dominating effect on ventral distance: higher aMag is associated with higher ventral distance.
Data availability
The model code is accessible at https://github.com/Bentley-Cellular-Adaptive-Behaviour-Lab/NeuralCrestCpp, (copy archived at swh:1:rev:cdf63f3786390fb7905092717456cc69f5657ddc). The code used to perform the LDA analysis is accessible in the supplementary files. All numerical data used in the figures is accessible in the supplementary data source file.
References
-
Zebrafish Neural Crest: Lessons and Tools to Study In Vivo Cell MigrationMethods in Molecular Biology (Clifton, N.J.) 2179:79–106.https://doi.org/10.1007/978-1-0716-0779-4_9
-
Cell cycle and cell fate in C. elegansDevelopment (Cambridge, England) 126:1947–1958.
-
Cell cycle, cytoskeleton dynamics and beyond: the many functions of cyclins and CDK inhibitorsCell Cycle (Georgetown, Tex.) 14:1786–1798.https://doi.org/10.1080/15384101.2014.998085
-
The role of differential VE-cadherin dynamics in cell rearrangement during angiogenesisNature Cell Biology 16:309–321.https://doi.org/10.1038/ncb2926
-
The temporal basis of angiogenesisPhilosophical Transactions of the Royal Society of London. Series B, Biological Sciences 372:20150522.https://doi.org/10.1098/rstb.2015.0522
-
Notch signalling in contextNature Reviews. Molecular Cell Biology 17:722–735.https://doi.org/10.1038/nrm.2016.94
-
Notch activates cell cycle reentry and progression in quiescent cardiomyocytesThe Journal of Cell Biology 183:129–141.https://doi.org/10.1083/jcb.200806104
-
Delta signaling mediates segregation of neural crest and spinal sensory neurons from zebrafish lateral neural plateDevelopment (Cambridge, England) 127:2873–2882.https://doi.org/10.1242/dev.127.13.2873
-
Notch in the pathway: the roles of Notch signaling in neural crest developmentSeminars in Cell & Developmental Biology 16:663–672.https://doi.org/10.1016/j.semcdb.2005.06.009
-
Asymmetric division coordinates collective cell migration in angiogenesisNature Cell Biology 18:1292–1301.https://doi.org/10.1038/ncb3443
-
Functional Versatility of the CDK Inhibitor p57Kip2Frontiers in Cell and Developmental Biology 8:584590.https://doi.org/10.3389/fcell.2020.584590
-
Notch1 signalling regulates endothelial proliferation and apoptosis in pulmonary arterial hypertensionThe European Respiratory Journal 48:1137–1149.https://doi.org/10.1183/13993003.00773-2015
-
Notch Signalling: The Multitask Manager of Inner Ear Development and RegenerationAdvances in Experimental Medicine and Biology 1218:129–157.https://doi.org/10.1007/978-3-030-34436-8_8
-
Mechanism and Significance of cis-Inhibition in Notch SignallingCurrent Biology: CB 21:R40–R47.https://doi.org/10.1016/j.cub.2010.10.034
-
Division of labor during trunk neural crest developmentDevelopmental Biology 344:555–565.https://doi.org/10.1016/j.ydbio.2010.04.009
-
An essential role for Notch in neural crest during cardiovascular development and smooth muscle differentiationThe Journal of Clinical Investigation 117:353–363.https://doi.org/10.1172/JCI30070
-
Coordinated control of Notch/Delta signalling and cell cycle progression drives lateral inhibition-mediated tissue patterningDevelopment (Cambridge, England) 143:2305–2310.https://doi.org/10.1242/dev.134213
-
Divide or Conquer: Cell Cycle Regulation of Invasive BehaviorTrends in Cell Biology 27:12–25.https://doi.org/10.1016/j.tcb.2016.08.003
-
her9 promotes floor plate development in zebrafishDevelopmental Dynamics 232:1098–1104.https://doi.org/10.1002/dvdy.20264
-
Endothelial signalling by the Notch ligand Delta-like 4 restricts angiogenesisDevelopment (Cambridge, England) 134:839–844.https://doi.org/10.1242/dev.003244
-
Apical migration of nuclei during G2 is a prerequisite for all nuclear motion in zebrafish neuroepitheliaDevelopment (Cambridge, England) 138:5003–5013.https://doi.org/10.1242/dev.071522
-
Notch signalling and the control of cell fate choices in vertebratesSeminars in Cell & Developmental Biology 9:583–589.https://doi.org/10.1006/scdb.1998.0266
-
Cyclin D1 functions in cell migrationCell Cycle (Georgetown, Tex.) 5:2440–2442.https://doi.org/10.4161/cc.5.21.3428
-
Second-generation Notch1 activity-trap mouse line (N1IP::CreHI) provides a more comprehensive map of cells experiencing Notch1 activityDevelopment (Cambridge, England) 142:1193–1202.https://doi.org/10.1242/dev.119529
-
Tumors of the neural crest: Common themes in development and cancerDevelopmental Dynamics 244:311–322.https://doi.org/10.1002/dvdy.24226
-
Alpha(v)beta3 integrin expression up-regulates cdc2, which modulates cell migrationThe Journal of Cell Biology 161:817–826.https://doi.org/10.1083/jcb.200212172
-
Notch pathway regulation of neural crest cell development in vivoDevelopmental Dynamics 241:376–389.https://doi.org/10.1002/dvdy.23717
-
Generation and application of signaling pathway reporter lines in zebrafishMolecular Genetics and Genomics 288:231–242.https://doi.org/10.1007/s00438-013-0750-z
-
Tactics of cancer invasion: solitary and collective invasionJournal of Biochemistry 167:347–355.https://doi.org/10.1093/jb/mvaa003
-
Cross-regulation between Notch and p63 in keratinocyte commitment to differentiationGenes & Development 20:1028–1042.https://doi.org/10.1101/gad.1406006
-
Notch activation induces endothelial cell cycle arrest and participates in contact inhibition: role of p21Cip1 repressionMolecular and Cellular Biology 24:8813–8822.https://doi.org/10.1128/MCB.24.20.8813-8822.2004
-
Cell-cycle regulation of NOTCH signaling during C. elegans vulval developmentMolecular Systems Biology 8:618.https://doi.org/10.1038/msb.2012.51
-
Actomyosin contractility and collective migration: may the force be with youCurrent Opinion in Cell Biology 48:87–96.https://doi.org/10.1016/j.ceb.2017.06.006
-
Oligodendrocyte specification in zebrafish requires notch-regulated cyclin-dependent kinase inhibitor functionThe Journal of Neuroscience 25:6836–6844.https://doi.org/10.1523/JNEUROSCI.0981-05.2005
-
A zebrafish histone variant H2A.F/Z and A transgenic H2A.F/Z:GFP fusion protein for in vivo studies of embryonic developmentDevelopment Genes and Evolution 211:603–610.https://doi.org/10.1007/s00427-001-0196-x
-
Angiogenesis: A team effort coordinated by notchDevelopmental Cell 16:196–208.https://doi.org/10.1016/j.devcel.2009.01.015
-
Tracking neural crest cell cycle progression in vivoGenesis (New York, N.Y.: 2000) 56:e23214.https://doi.org/10.1002/dvg.23214
-
Notch Signaling in Embryology and Cancer, Advances in Experimental Medicine and BiologyNotch Signaling in Embryology and Cancer, Notch Signaling in Embryology and Cancer, Advances in Experimental Medicine and Biology, New York, NY, Springer Science+Business Media, 10.1007/978-1-4614-0899-4.
-
Collective cell migrationAnnual Review of Cell and Developmental Biology 25:407–429.https://doi.org/10.1146/annurev.cellbio.042308.113231
-
Notch signaling regulates growth and differentiation in the mammalian lensDevelopmental Biology 321:111–122.https://doi.org/10.1016/j.ydbio.2008.06.002
-
Use of the Gal4-UAS technique for targeted gene expression in the zebrafishMechanisms of Development 80:153–158.https://doi.org/10.1016/s0925-4773(98)00209-3
-
Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
-
Collective cell migration of epithelial and mesenchymal cellsCellular and Molecular Life Sciences 70:3481–3492.https://doi.org/10.1007/s00018-012-1251-7
-
Functional analysis of Hairy genes in Xenopus neural crest initial specification and cell migrationDevelopmental Dynamics 244:988–1013.https://doi.org/10.1002/dvdy.24295
-
Cell-cycle-dependent regulation of cell motility and determination of the role of Rac1Experimental Cell Research 295:407–420.https://doi.org/10.1016/j.yexcr.2004.01.011
-
Cell cycle regulation by the intrinsically disordered proteins p21 and p27Biochemical Society Transactions 40:981–988.https://doi.org/10.1042/BST20120092
-
Antagonistic regulation of p57kip2 by Hes/Hey downstream of Notch signaling and muscle regulatory factors regulates skeletal muscle growth arrestDevelopment (Cambridge, England) 141:2780–2790.https://doi.org/10.1242/dev.110155
-
On the Molecular Mechanisms Regulating Animal Cell Size HomeostasisTrends in Genetics 36:360–372.https://doi.org/10.1016/j.tig.2020.01.011
Article and author information
Author details
Funding
Medical Research Council (G1000080/1)
- Claudia Linker
Royal Society (2010/R1)
- Claudia Linker
Wellcome Trust (207630/Z/17/Z)
- Claudia Linker
Eucine Kennedy Shiver National Institute of Child Health & Human Development of the National Institutes of Health (T32HD055164)
- Manuel Rocha
Eucine Kennedy Shiver National Institute of Child Health & Human Development of the National Institutes of Health (F31HD097957)
- Manuel Rocha
Cancer Research UK (FC001751)
- Dylan Feldner-Busztin
- Katie Bentley
Medical Research Council (FC001751)
- Dylan Feldner-Busztin
- Katie Bentley
Wellcome Trust (FC001751)
- Dylan Feldner-Busztin
- Katie Bentley
Biotechnology and Biological Sciences Research Council (BB/S015906/1)
- Robert N Kelsh
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication. For the purpose of Open Access, the authors have applied a CC BY public copyright license to any Author Accepted Manuscript version arising from this submission.
Acknowledgements
In memory of Julian Lewis. We are especially grateful to N Daudet for his scientific and personal support. To the KCL fish facility staff, particularly to J Glover. We are grateful to Y Hinits and S Wilson for sharing reagents. This project was funded by MRC G1000080/1, Royal Society 2010/R1 and Wellcome Trust 207630/Z/17/Z to CL; MR was supported by the Eunice Kennedy Shriver National Institute of Child Health & Human Development of the National Institutes of Health under awards T32HD055164 and F31HD097957. KB and DF-B were supported by the Francis Crick Institute core funding from CRUK (FC001751), MRC (FC001751), and Wellcome Trust (FC001751); BBSRC BB/S015906/1 to RNK.
Ethics
Zebrafish were maintained in accordance with UK Home Office regulations UK Animals (Scientific Procedures) Act 1986, amended in 2013 under project license P70880F4C.
Copyright
© 2022, Alhashem 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,781
- views
-
- 521
- downloads
-
- 15
- 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
-
- Developmental Biology
Dysfunction of primary cilia leads to genetic disorder, ciliopathies, which shows various malformations in many vital organs such as brain. Multiple tongue deformities including cleft, hamartoma, and ankyloglossia are also seen in ciliopathies, which yield difficulties in fundamental functions such as mastication and vocalization. Here, we found these tongue anomalies in mice with mutation of ciliary protein. Abnormal cranial neural crest-derived cells (CNCC) failed to evoke Hh signal for differentiation of mesoderm-derived cells into myoblasts, which resulted in abnormal differentiation of mesoderm-derived cells into adipocytes. The ectopic adipose subsequently arrested tongue swelling formation. Ankyloglossia was caused by aberrant cell migration due to lack of non-canonical Wnt signaling. In addition to ciliopathies, these tongue anomalies are often observed as non-familial condition in human. We found that these tongue deformities could be reproduced in wild-type mice by simple mechanical manipulations to disturb cellular processes which were disrupted in mutant mice. Our results provide hints for possible future treatment in ciliopathies.
-
- Developmental Biology
- Genetics and Genomics
Cilia defects lead to scoliosis in zebrafish, but the underlying pathogenic mechanisms are poorly understood and may diverge depending on the mutated gene. Here, we dissected the mechanisms of scoliosis onset in a zebrafish mutant for the rpgrip1l gene encoding a ciliary transition zone protein. rpgrip1l mutant fish developed scoliosis with near-total penetrance but asynchronous onset in juveniles. Taking advantage of this asynchrony, we found that curvature onset was preceded by ventricle dilations and was concomitant to the perturbation of Reissner fiber polymerization and to the loss of multiciliated tufts around the subcommissural organ. Rescue experiments showed that Rpgrip1l was exclusively required in foxj1a-expressing cells to prevent axis curvature. Genetic interactions investigations ruled out Urp1/2 levels as a main driver of scoliosis in rpgrip1 mutants. Transcriptomic and proteomic studies identified neuroinflammation associated with increased Annexin levels as a potential mechanism of scoliosis development in rpgrip1l juveniles. Investigating the cell types associated with annexin2 over-expression, we uncovered astrogliosis, arising in glial cells surrounding the diencephalic and rhombencephalic ventricles just before scoliosis onset and increasing with time in severity. Anti-inflammatory drug treatment reduced scoliosis penetrance and severity and this correlated with reduced astrogliosis and macrophage/microglia enrichment around the diencephalic ventricle. Mutation of the cep290 gene encoding another transition zone protein also associated astrogliosis with scoliosis. Thus, we propose astrogliosis induced by perturbed ventricular homeostasis and associated with immune cell activation as a novel pathogenic mechanism of zebrafish scoliosis caused by cilia dysfunction.