Gastrulation pharmacological screen and in silico modelingWe previously established that non-apical mitoses at the PS primarily give rise to basal delamination of nascent mesoderm cells13. Therefore, we set out to identify the factors controlling the occurrence of non-apical divisions and assess whether they are necessary and/or sufficient for cell delamination. Variations in the position within the epithelium at which a cell divides could come from differences in apicobasal polarity, actomyosin contractility or cell cycle parameters (such as total duration of the cycle or duration of the G2/M phase). Thus, we performed a pharmacological screen targeting these cellular characteristics. While genetic manipulation allows for spatial control, it lacks precision in intensity and duration, which can complicate the interpretation of a phenotype, particularly in a highly dynamic context. The pharmacological approach addresses these limitations as it allows establishing optimal posology and recording the consequences of short-term alterations. One caveat is that all cells (epiblast and PS) are simultaneously targeted. To circumvent this, we used a computational model of proliferating pseudostratified epithelium in which modifications of cell cycle parameters, attachment to the basement membrane, and fine regulation of apical and medial actomyosin can be adjusted at the single cell level.For biological experiments, we established a systematic protocol of measurements on embryos sectioned along the proximal-distal axis (Fig. 1a, b). Sections were positioned along the anterior-posterior axis based on anatomical landmarks such as the long axis of the oval embryo shape and the morphology of visceral endoderm cells, as well as molecular markers including collagen IV for the basement membrane. The epiblast was divided into three regions: anterior and posterior halves of the epiblast, and PS, defined as the subregion of the posterior epiblast where the basement membrane is discontinuous. Staining for mitosis (Phospho-Histone H3) and nuclei (DAPI) allowed calculating the mitotic and non-apical mitotic indexes (Fig. 1b). We quantified the percentage of basement membrane (collagen IV) degradation as a readout for EMT initiation and progression (Fig. 1a).Fig. 1: Calibration of pharmacological screening and in silico modeling.a Graphic representation of a transverse section from a E6.5 embryo. The epiblast is divided in three regions: anterior (pink), posterior (green) and the primitive streak (PS, blue). Anterior and posterior regions are defined by dividing the epiblast in halves. The PS region is the subsection of the posterior region where the basement membrane (BM) is discontinuous (Collagen IV, yellow). Mesoderm cells that exited the epiblast are in purple. The black double head arrow represents the distance in µm between the apical and the basal surface. b Graphic representation of a transverse section from a E6.5 embryo. Nuclei are in purple, apical mitosis in pink and non-apical mitosis in green. A mitosis is considered non-apical if it occurs at least 10 µm away from the apical surface (lumen), in a cell that retains apical and basal contact. c Diagram representing simulated cells controlled by various springs: apical cell adhesions (apical-to-apical spring), and the cytoplasmic cytoskeleton with basal cell-matrix adhesion (nucleus-to-basal spring and apical-to-nucleus spring). Cell-matrix adhesions are modeled by the attachment of the nucleus-to-basal spring to the basal line. Nuclei are modeled with two layers. The soft layer can adjust to the environment to model pseudo-stratified morphology, the hard core cannot be deformed and increases during mitosis. d Representation of a simulation setup (t = 0 h) for the in silico model adapted to the epiblast dynamics, with 19 PS (purple nuclei) cells and 40 epiblast cells (green nuclei). The red line represents the apical surface, and the black line the basal surface.For the modeling of pseudo-stratified epithelium, we relied on a model originally developed on the chick neuroepithelium20. We used an extended version that allows the control of each parameter in time and space at the single cell level, described in ref. 21. Briefly, each cell is approximated by a nucleus, an apical point, and a basal point. The two points are connected to the nucleus through dynamically adjustable springs representing the viscoelastic properties of the cytoplasm (Fig. 1c). The nucleus is composed of an outer soft core and an inner hard core. Hard cores of neighboring cells cannot overlap while soft cores can but are subjected to a repulsion force. This allows representing packing and compression occurring between cells in a dense tissue. Apical cell-cell adhesions are modeled by springs between the apical points of neighboring cells (aa). On the basal side, basal points are attached to a fixed basal line representing the basement membrane. Basal points cannot overlap or swap positions. When a cell detaches its basal point in the model, it is equivalent to a local fenestration of the matrix as it locally increases the mean distance between the basal points of the nearest neighbors. The cell cycle and INM movements are divided into three phases. The first phase corresponds to G1, S, and early G2, where movements are controlled by environmental constraints generated by neighboring cells. The second phase corresponds to late G2, where the pre-mitotic rapid apical movement (or PRAM) is simulated by contraction of the apical-nuclei string (aN) to bring nuclei to the apical pole, and extension of the basal-nuclei spring (bN) to accompany the movement. These settings are maintained in the third phase, the M phase, which results in the formation of two daughter cells. As the model represents a 2D view of a tissue in 3D, only a few daughter cells would actually divide exactly within the simulated 2D plane. As a dimensional correction, we neglect the second daughter cell in 95% of the cases. To adapt the model to the epiblast size and configuration, a maximum height was set to allow no more than three pseudo-layers of nuclei (Fig. 1d). A user-friendly online version of this model is available as a web-based stand-alone EMT simulator (sEMTor) at https://semtor.github.io/gastrulation/.Basement membrane degradation triggers non-apical mitosesIn the mouse embryo, perforations in the basement membrane separating epiblast and visceral endoderm are distributed uniformly prior to migration of the Anterior Visceral Endoderm (AVE), the organizer that determines the anterior-posterior axis and therefore the position of the streak. Upon axis specification, perforations become more abundant in the posterior half of the embryo. These holes in the basement membrane are generated by a Nodal-dependent increase in the expression of Matrix Metalloproteinases (MMP) 2 and 14, and are necessary for embryo growth, PS morphogenesis, and gastrulation5. In EMT, loss of basal attachment due to basement membrane degradation is frequently associated with a disruption in apical-basal polarity. However, during mouse gastrulation, polarity markers are present in all PS cells4,6,22. We wondered whether basement membrane degradation might nonetheless be required for the occurrence of non-apical mitosis, possibly through modification of mechanical constraints in the epiblast.To test this, we cultured pre-gastrulation E6 embryos for 6 hours with a combination of 20 µM Pan-MMPs inhibitor prinomastat hydrochloride and 100 µM (NSC405020) MMP14 inhibitor, as previously described5. After incubation, embryos were fixed, and transverse cryosections were immunostained for Phospho-Histone H3, Collagen IV, nuclei (DAPI), and F-actin (phalloidin). Additionally, adjacent sections were stained for Cerberus 1 (Cer1) to identify the location of the AVE. We confirmed the abrogation of basement membrane degradation in embryos treated with MMPs inhibitors (Fig. 2a–c). However, there was no significant difference in mitotic and non-apical mitotic indexes (Fig. 2d, e), nor in the non-apical mitosis ratio (Supplementary Fig. 1a) throughout the epiblast relative to untreated embryos.Fig. 2: Preventing basement membrane degradation does not affect mitosis position.Z-projection of transverse sections from representative E6 embryos cultured for 6 hours with vehicle or MMPs inhibitors, stained for (a) nuclei (DAPI, cyan), a mitosis (Phospho-histone H3 (Phh3), magenta), a, b basement membrane (collagen IV, yellow), and (a, b) F-actin (phalloidin, grey). Arrowheads indicate non-apical mitoses. Scale bars in main panels: 25 µm. Zooms in b show the basement membrane at the PS. Scale bars: 10 µm. c Violin plots representing the percentage of basement membrane degradation in embryos cultured with control vehicle or MMPs inhibitors. Inner doted lines indicate third quartile, median, and first quartile. Normality was assessed using a Shapiro–Wilk test followed by a two-tailed Mann Whitney test. Split violin plots representing the mitotic index (d) and non-apical mitotic index (e), in anterior, posterior and PS regions embryos cultured with control vehicle or MMPs inhibitors. Normality was assessed using a Shapiro–Wilk test followed by a two-tailed Mann-Whitney or unpaired t-test. Controls: n = 14 slides from 8 embryos; MMPs inhibitors: n = 21 slides from 7 embryos. f Graph from in silico simulations representing the percentage of non-apical mitoses depending on aN strength in epiblast cells (6 h cell cycle) and PS cells (3 h cell cycle), 100 simulations with 30 proliferating cells. g Table presenting the simulated parameters set for epiblast and PS cells in control conditions in the in silico epiblast model. Event B designates the detachment of cells’ basal point from the basal line. h Scatterplot representing the ratio of non-apical mitosis (%), in the region of the epiblast and PS in simulations ran with control settings or in situations where PS cells do not undergo event B. For each condition, 100 simulations with 19 PS cells surrounded by 40 control epiblast cells. Error bars: SEM. (i) Table recapitulating the settings of epiblast and PS cells in control and test simulations. For the simulation, only the modified parameters are indicated. A-P Anterior – Posterior, Ap-Bas Apical – Basal, PS Primitive Streak, CC cell cycle, aN apical-Nuclei spring, bN basal-Nuclei spring, aa apical adhesion spring.Before modelling the impact of detachment from the basement membrane in silico, we first had to calibrate the parameters of cell cycle and INM allowing to recapitulate the percentage of non-apical mitoses observed in vivo. As the velocity of PRAM is controlled by the contraction of the apical-to-Nucleus (aN) spring in the model, we tested a range of strength values for aN (Fig. 2f). We simulated 200 repetitions of 30 proliferating cells with a cell cycle of 6 ± 0.5 h or 3 ± 0.5 h, corresponding to measured values for epiblast and PS cells, respectively12,15. The percentages of non-apical mitoses observed at the PS (circa 40%) could only be reached at very low values of aN, suggesting that an almost complete loss of PRAM is necessary. In addition, we noted that a faster cell cycle favors non-apical mitoses even at higher PRAM intensity. Indeed, the frequency of non-apical mitoses for medium values of aN reached a plateau at 7% with a cell cycle length of 3 h, and only 1% with a length of 6 h. Based on these simulations, we used the following as default parameters: 6 ± 0.5 h of cell cycle and 3 of aN strength for epiblast cells, 3 ± 0.5 h of cell cycle and 0.25 of aN strength for PS cells (Fig. 2g, Supplementary Movie 1). We counted the number of PS cells on cross sections of stage E6.5 mouse embryos and found a mean population size of 19 PS cells per section13. Thus, we set the size of the PS population at 19 cells, surrounded by 20 epiblast cells on each side. Simulations start with a period of 6 hours during which all cells are attached basally. Then there is an 18 h long window of opportunity for PS cells to detach from the basal line, which is referred as event B21 (Fig. 2g) and represents basement membrane degradation21.We then ran simulations with or without detachment of PS cells from the basal line and plotted the rate of non-apical mitoses (Fig. 2h, i; Supplementary Movie 1). As observed in embryos treated for 6 h with MMPs inhibitors, impairing detachment from the basal line did not suppress the occurrence of non-apical mitoses at the PS (Fig. 2h, i).Collectively, these data indicate that preventing basement membrane degradation does not influence the rate of mitoses and non-apical mitoses at the PS.To determine if basement membrane degradation was sufficient to trigger the appearance of non-apical mitoses, we examined the position and frequency of mitoses in embryo deficient for FLRT3 (Fibronectin leucin-rich transmembrane protein 3). Knockout of FLRT3 in mouse leads to disorganization of the basement membrane between visceral endoderm and epiblast and the emergence of ectopic sites of epiblast delamination23. Since basement membrane degradation could no longer be used to define the region of the streak in mutant embryos, we approximated the position of the PS as the posterior 1/3rd of the epiblast, opposite the AVE. The rest of the epiblast, referred as anterior 2/3rd, approximates the anterior and posterior epiblast (Supplementary Fig. 2a, b). We detected no difference regarding the frequency of mitosis in WT vs FLRT3KO embryos (Supplementary Fig. 2c). However, there was a significant increase in the proportion of non-apical mitosis frequency in the anterior 2/3rd epiblast in FLRT3KO embryos, compared to wild-type embryos (Supplementary Fig. 2d). In ectopic sites of delamination in FLRT3KO embryos, 45% of mitoses happened apically, 30% non-apically, and 25% in delaminated cells (Supplementary Fig. 2e, f).These data indicate that ectopic basement membrane degradation is sufficient to trigger non-apical mitosis, associated with cell delamination, in regions of the epiblast that are presumably not exposed to the signaling pattern characteristic of the PS region.Actomyosin relaxation drives non-apical mitosesAs INM movements are usually driven by actomyosin in short pseudo-stratified epithelia16, we tested whether a regional modification in actomyosin contractility could be responsible for the increased rate of non-apical mitosis at the PS.The shape of nuclei in interphase can be used as a readout of the mechanical constraints exerted on cells24,25, which are believed to be particularly high in pseudo-stratified epithelia26,27. We marked nuclear membrane via staining for Lamin B1, the first Lamin subtype expressed during early development in mouse28,29, and calculated the nuclei aspect ratio by dividing the short axis length by the long axis length (Fig. 3a). A ratio close to 0 corresponds to an elongated nucleus, indicating that the cytoskeleton is contracted. A ratio close to 1 indicates a round nucleus, reflecting relaxation of the cytoskeleton network. Nuclei aspect ratios were recorded in the anterior, posterior, and PS regions of E6.25, E6.5, and E7 embryos (Fig. 3b-c”). Nuclei in the PS were rounder, compared to the rest of the epiblast, from E6.25 onwards, with a peak at E6.5 and maintenance until E7 (Fig. 3c-c”). This is compatible with a local relaxation of the cytoskeleton at the streak, leading to nuclei shape changes, upon the onset of gastrulation.Fig. 3: Actomyosin relaxation at the PS triggers the disharmonizing of cell cycle and pre-mitotic rapid apical movements.a Nucleus representation. b-b’ Z-projection of a transverse section from a E7 embryo, stained for Lamin B1 (green) (b), with zooms on nuclei from the anterior, posterior and PS regions (b’). Scale bars: 25 µm (b) and 10 µm. c-c” Violin plots representing nuclei aspect ratio at E6.25 (c), E6.5 (c’) and E7 (c”). Inner doted lines indicate third quartile, median, and first quartile. Normality was assessed using a Shapiro–Wilk test followed by a two-tailed Kruskall-Wallis test with Dunnet post hoc. E6.25: n = 4 embryos, E6.5: n = 4 embryos, E7: n = 4 embryos, with 20 nuclei per region. d Z-projection of transverse sections from E6.5 embryos cultured for 2 hours with vehicle or ROCK inhibitor at 40 µM, stained for nuclei (DAPI, cyan), basement membrane (collagen IV, yellow), mitosis (Phh3, magenta) and F-actin (phalloidin, grey). Arrowheads indicate non-apical mitoses. Scale bars: 25 µm. Split violin plots representing the mitotic index (e) and non-apical mitotic index (f), in embryos cultured with control vehicle or 40 µM ROCK inhibitor. Normality was assessed using a Shapiro–Wilk test followed by a two-tailed Mann-Whitney. Control: n = 26 slides from 12 embryos; 40 µM ROCK inhibitor: n = 14 slides from 7 embryos. g Violin plots representing nuclei aspect ratio in embryos cultured with vehicle or 40 µM ROCK inhibitor. Normality was assessed using a Shapiro–Wilk test followed by a two-tailed Mann-Whitney or unpaired t-test. Control: n = 11 embryos; 40 µM ROCK inhibitor: n = 5 embryos. h Scatterplot representing the percentage of non-apical mitosis in the region of the epiblast and PS in various simulations: control, relaxation of INM (aN/bN) springs, relaxation of all springs and relaxation of apical adhesion springs. For each condition, 100 simulations with 19 PS cells surrounded by 40 control cells. Error bars: SEM. i Table recapitulating the settings of epiblast and PS cells in control and test simulations. For the simulations, only the modified parameters are indicated. A-P Anterior – Posterior, Ap-Bas Apical – Basal, PS primitive streak, INM interkinetic nuclei migration, CC cell cycle, B detachment of basal adhesions from the basal line, aN apical-Nuclei spring, bN basal-Nuclei spring, aa apical adhesion spring.The small Rho-GTPase RhoA regulates both actomyosin contractility and cell cycle progression. In addition, during EMT, RhoA participates in the cellular changes required to acquire a mesenchymal-like phenotype and migrate30. Rho-associated kinase (ROCK), a downstream effector of RhoA, regulates actomyosin contractility via myosin light chain phosphorylation31,32. In the pseudostratified epithelium of the Drosophila imaginal discs, pharmacological treatment with ROCK inhibitor Y27632 leads to basal accumulation of nuclei and the appearance of non-apical mitoses33. To investigate the link between mitosis location and actomyosin dynamics in the mouse epiblast, we cultured E6.5 embryos for 2 h with 40 µM Y27632. Embryos were then fixed, and transverse cryosections were immunostained for Phospho-Histone H3, Collagen IV, nuclei (DAPI), and F-actin (phalloidin) (Fig. 3d). There was no change in the mitotic index (Fig. 3e). However, a 2-fold increase in the non-apical mitotic index (and a 10-fold increase of the non-apical mitosis ratio) were recorded in the anterior and posterior epiblast (Fig. 3f, Supplementary Fig. 1b). This is compatible with an actomyosin-driven INM in the epiblast: indeed, it is likely that disturbing actomyosin contractility affects the ability of cells to drive their nuclei to the apical surface in G2, which results in mitosis occurring in ectopic location. Interestingly, non-apical mitotic index remained unchanged at the PS (Fig. 3f), strongly suggesting that actomyosin contractility is intrinsically low in PS cells. Upon ROCK inhibition, nuclei shape was rounder in the anterior and posterior epiblast confirming it induced a relaxation of the actomyosin cytoskeleton, while the shape of nuclei in the PS did not change (Fig. 3g). The effect of the inhibitor was reversible, as embryos treated for 2 h with 40 µM of ROCK inhibitor then washed were able to adjust rapidly and survive overnight in culture (Supplementary Fig. 1e).To further explore the influence of cytoskeleton contractility on the occurrence of non-apical mitoses, we performed simulations with a reduced strength for all springs: aN and bN for the cell body/cytoplasm and aa for apical cell-cell adhesion (Fig. 3h, i; light purple). Interestingly, such global relaxation increased the frequency of non-apical mitoses in the epiblast and reduced it in the PS (Fig. 3h, Supplementary Movie 2). Then, we set out to assess the relative contributions of the cytoplasmic springs (aN, bN) that control INM/PRAM and of the apical adhesion spring (aa) by relaxing them independently (Fig. 3h, i; pink and orange; Supplementary Movie 2). Relaxing aN/bN strongly promoted non-apical mitoses in the epiblast but did not affect PS cells. Relaxing aa mildly promoted non-apical mitoses in the epiblast and lowered their frequency in the PS. It is interesting to note that in all simulations with aa relaxation the tissue experienced severe deformations (Fig. 3h, i; Supplementary Movie 2). The changes obtained in the scenario where all springs were relaxed are compatible with what we observed in embryos treated with a high dose of ROCK inhibitor (400 µM): those displayed severe morphological defects (Supplementary Fig. 3a, b); the mitotic index was decreased in all regions of the epiblast (Supplementary Fig. 3c) while the non-apical mitotic index was increased in the anterior and posterior epiblast but decreased in the PS (Supplementary Fig. 3d). We interpret that the effects observed at 40 µM of ROCK inhibitor correspond to a mild inhibition of the cell cytoskeleton contractility that can be modeled by the relaxation of only the INM springs, while the buckling effect seen at 400 µM requires stronger perturbation including that of the apical actomyosin.Collectively, these experimental and simulated data indicate that non-apical mitoses in the epiblast can be elicited by affecting cytoplasmic actomyosin activity rather than apical actomyosin. This suggests that INM is governed by actomyosin in the epiblast and that, comparatively to the rest of the epiblast, actomyosin contractility is weaker in PS cells, leading to the appearance of non-apical mitoses.Cell cycle regulation is less stringent in the PSThe complex formed by Cyclin-dependent Kinase 1 (CDK1) and CyclinD1 is critical for the control of cell entry into mitosis. In the avian neuroepithelium, overexpression of CyclinD1 leads to the appearance of non-apical mitoses34. In zebrafish retina, CDK1 is required for the rapid basal to apical movement of nuclei during G235. Therefore, we chose to pharmacologically target CDK1, negatively using a CDK1 inhibitor (R03366), and positively via an inhibitor (PD166285) of Wee1, a nuclear kinase that regulates mitosis by phosphorylating CDK1. Wee1 inhibition thereby enables to bypass the G2/M checkpoint36.A 2 h incubation with 175 µM of CDK1 inhibitor abrogated all mitosis, confirming that targeting CDK1 regulation can rapidly affect the cell cycle in the epiblast (Fig. 4a–c). Embryos were able to recover from the treatment and developed normally (Supplementary Fig. 1e). With the goal to increase the frequency of mitoses without major morphological anomaly, we performed a dose-response analysis for the Wee1 inhibitor (5, 10, 20, 40, and 80 µM for 2 hours). Embryos were fixed directly after incubation. We measured the percentage of basement membrane degradation as well as the mitotic and non-apical mitotic indexes on transverse cryosections immunostained for Phospho-Histone H3, Collagen IV, nuclei, and F-actin (Fig. 4d). In the anterior and posterior epiblast, we observed a dose-dependent increase of the mitotic index, indicating an acceleration of the cell cycle, as well as an increase of the non-apical mitotic index and non-apical mitosis ratio (Fig. 4e, f; Supplementary Fig. 1c). We previously showed that an acceleration of the cell cycle is sufficient to trigger non-apical mitoses (Fig. 2f), therefore some of the increase of non-apical mitoses in the epiblast upon treatment with Wee1 inhibitor might be linked to the global increase in mitotic index. However, as PRAM takes place in G2, Wee1 inhibition might also contribute to the higher rate of non-apical mitoses by shortening the window of opportunity for PRAM to occur. By contrast, the frequency and location of mitoses in the PS stayed constant except at high doses of drug (40–80 µM) at which most cells were entering mitosis (Fig. 4e, f). These data indicate that cells at the PS were less sensitive to Wee1 inhibition compared to the rest of the epiblast, strongly suggesting an intrinsically shorter cell cycle and G2 duration in PS cells.Fig. 4: Cell cycle regulation at the PS is less stringent than in the rest of the epiblast, which triggers the appearance of non-apical mitoses.a Z-projections of transverse sections from E6.5 embryos cultured for 2 hours with control vehicle or 175 µM of CDK1 inhibitor (R03366), stained for nuclei (DAPI, cyan), basement membrane (collagen IV, yellow), mitosis (phh3, magenta) and F-actin (phalloidin, grey). Scale bars: 25 µm. Violin plots representing the mitotic index (b) and non-apical mitotic index (c) in anterior, posterior, and PS regions in embryos cultured with control vehicle or 175 µM CDK1 inhibitor. Inner doted lines indicate third quartile, median, and first quartile. Normality was assessed using a Shapiro–Wilk test followed by a two-tailed Mann-Whitney. Control: n = 24 slides from 8 embryos; 175 µM CDK1 inhibitor: n = 11 slides from 7 embryos. d Z-projection of transverse sections from E6.5 embryos cultured for 2 hours in vehicle or with 20 µM, 40 µM or 80 µM of Wee1 inhibitor, stained for nuclei (DAPI, cyan), basement membrane (collagen IV, yellow), mitosis (phh3, magenta) and F-actin (phalloidin, grey). Arrow heads indicate non-apical mitoses. Scale bars: 25 µm. Violin plots representing the mitotic index (e) and non-apical mitotic index (f), in anterior, posterior, and PS regions of embryos cultured with control vehicle or 5 µM, 10 µM, 20 µM, 40 µM or 80 µM of Wee1 inhibitor for 2 hours. Ns nonsignificant, *P-value ≤ 0.05, **P-value ≤ 0.01 and ***P-value ≤ 0.001. Normality was assessed using a Shapiro–Wilk test followed by a two-tailed Mann-Whitney or unpaired t-test. Control: n = 27 slides from 13 embryos; 5 µM Wee1 inhibitor: n = 12 slides from 5 embryos; 10 µM Wee1 inhibitor: n = 16 slides from 10 embryos; 20 µM Wee1 inhibitor: n = 21 slides from 9 embryos; 40 µM Wee1 inhibitor: n = 9 slides from 4 embryos; 80 µM Wee1 inhibitor: n = 12 slides from 6 embryos. g Scatterplot representing the percentage of non-apical mitosis in the region of the epiblast and PS in various simulations: control, reduction of the total cell cycle duration (3 h), reduction of G2 and M phase duration (15 minutes) and reduction of both the total duration of the cell cycle (3 h) and of G2 and M duration (15 minutes). For each condition, 100 simulations with 19 PS cells surrounded by 40 control cells. Error bars: SEM. h Table recapitulating the settings of epiblast and PS cells in control and test simulations. For the simulation, only the modified parameters are indicated. A-P Anterior – Posterior, Ap-Bas Apical – Basal, PS primitive streak, CC cell cycle, B detachment of basal adhesions from the basal line, aN apical-Nuclei spring, bN basal-Nuclei spring, aa apical adhesion spring.To assess the relative impact of G2 and total cell cycle duration, we performed simulations in which we modulated their duration independently or simultaneously. The model has a simplified version of the cell cycle in which PRAM forces are applied through G2 and M. Thus, to mimic the effect of the Wee1 inhibitor in shortening biological G2 and the associated duration of PRAM forces, we reduced the total duration of simulated G2 + M by half. In the control situation, epiblast cells cycle every 6 ± 0.5 h and PS cells every 3 ± 0.5 h. All cells have a G2/M of 30 minutes each (Fig. 4g, h; dark purple). If all cells were forced to cycle every 3 ± 0.5 h at constant G2/M duration, it only generated a marginal increase of non-apical mitoses in the epiblast (Fig. 4g, h; light purple). If both cell types kept their respective cell cycle length but were forced to have a G2/M of 15 minutes each, there was a massive increase of non-apical mitoses in the epiblast (Fig. 4g, h; pink, Supplementary Movie 3). If both changes were combined (cell cycle 3 ± 0.5 h, G2/M 15 min in all cells), it did not generate a further increase of non-apical mitoses in the epiblast (Fig. 4g, h; orange; Supplementary Movie 3). By contrast, all these changes had negligible effects on the PS. These results indicate that specifically affecting G2/M duration has more impact on the rate of non-apical mitoses than shortening the whole cell cycle without affecting G2/M duration. Overall, these data indicate that an acceleration of the cell cycle and a shortened G2 duration contribute to the increased incidence of mitoses and non-apical mitoses observed at the PS.Non-apical mitoses occur independently of EMT progressionThere is evidence that mitosis can influence gene regulation and EMT progression, notably through the process of gene bookmarking37,38. Therefore, we asked whether mitosis frequency and location could favor EMT.In control conditions, non-apical mitoses are followed by basal delamination of mesodermal precursors13. Thus, we wondered if any of the pharmacological treatments that led to ectopic non-apical mitoses might be sufficient to promote ectopic EMT in the epiblast. We monitored markers of EMT progression, notably the expression of the EMT marker Snail, the pattern of E-cadherin and N-cadherin, and the percentage of basement membrane degradation, upon pharmacological inhibition of ROCK, Wee1, and CDK1.In embryos treated with 40 µM of ROCK inhibitor, immunostaining for Snail showed that cells at the PS were still able to transition, and no ectopic EMT site was detected (Fig. 5a). This was confirmed by staining for E-cadherin and N-cadherin, as there was no alteration in the cadherin switch (Fig. 5b). Similarly, we observed no significant change regarding the percentage, pattern, or localization of basement membrane degradation (Figs. 3d and 5c).Fig. 5: Dynamic changes observed at the PS are independent from epithelial-mesenchymal transition progression.Z-projections of transverse sections from E6.5 embryos cultured for 2 hours with control vehicle, 40 µM ROCK inhibitor, 20 µM Wee1 inhibitor or 175 µM CDK1 inhibitor, stained for (a) EMT marker Snail (yellow) and F-actin (phalloidin, grey), or (b) E-cadherin (cyan) and N-cadherin (magenta). Scale bars: 25 µm. Representative images from Control n = 9 embryos; 40 µM ROCK inhibitor n = 3 embryos; 20 µM Wee1 inhibitor: n = 5 embryos, 175 µM CDK1 inhibitor: n = 4 embryos. c-c” Violin plots representing the percentage of basement membrane degradation in embryos cultured with control vehicle, 40 µM ROCK inhibitor (c), 20 µM Wee1 inhibitor (c’), 175 µM CDK1 inhibitor (c”). Inner doted lines indicate third quartile, median, and first quartile. Normality was assessed using a Shapiro–Wilk test followed by a two-tailed unpaired t-test. Control: n = 27 slides from 13 embryos; 40 µM ROCK inhibitor: n = 21 slides from 9 embryos; 175 µM CDK1 inhibitor: n = 11 slides from 7 embryos. d In situ hybridization for Brachyury (T) (top), and Eomesodermin (bottom) on E7.5 WT, RhoAΔepiblast, Ect2Δepiblast embryos. Scale bars: 100 µm. Brachyury (T): RhoAΔepiblast n = 2; Ect2Δepiblast n = 4. Eomesodermin: RhoAΔepiblast n = 5; Ect2Δepiblast n = 5. e Z-selection from E7.5 WT (n = 7), RhoAΔepiblast(n = 4), Cdk1Δepiblast(n = 5), Ect2Δepiblast(n = 6) embryos stained for F-actin (Phalloidin, grey) and Vimentin (green). Scale bar: 50 µm.For Wee1, we focused on embryos treated with 20 µM of the inhibitor for 2 hours. Immunostaining for Snail, E-cadherin, and N-cadherin indicated that EMT was initiated normally at the PS (Fig. 5a, b). While these conditions led to a global increase of non-apical mitoses in all the epiblast (Fig. 4f; Supplementary Fig. 1c), there were no ectopic site of EMT nor significant change regarding the percentage of basement membrane degradation (Figs. 4d and 5a, b, c’). In addition, preventing entry in mitosis by culturing embryos for 2 hours with 175 µM of CDK1 inhibitor did not affect basement membrane degradation, the expression of Snail or the cadherin switch (Figs. 4a and 5a, b, c”), indicating that cell division per se is not necessary for EMT initiation.In addition, we examined the phenotypes of embryos carrying epiblast-specific deletion of the small GTPase RhoA, the Guanine exchange factor Ect2, and the kinase Cdk1, which are involved in cytoskeleton regulation, cytokinesis, and cell cycle regulation, respectively. Δepiblast embryos were obtained from the combination of RhoA, Ect2 and Cdk1 conditional alleles39,40,41 and a Cre-recombinase under the dependency of the Sox2 promoter, which is turned on in the prospective epiblast at E3.542 (Fig. 5d, e, Supplementary Fig. 4).Deletion of RhoA in the epiblast led to disruption of the pseudostratified arrangement, detectable from E6.5 onwards (Supplementary Fig. 4a–g). We found a slightly higher mitotic index at E6.5, possibly due to delays in cytokinesis (Supplementary Fig. 4c, d). Basement membrane degradation was detectable in the posterior region, and cells with a mesenchymal phenotype accumulated at the streak and into the lumen (Supplementary Fig. 4e). This is compatible with an adequate initiation of EMT but a defect in cell migration away from the PS, like what was observed in Rac1 epiblast-specific mutants22, or Rac1 and RhoA mesoderm-specific mutants43. By E7.5, RhoAΔepiblast mutants exhibited severe morphological defects. Lumens could be spotted within the accumulation of cells on the posterior side of the embryo, suggesting a disruption of polarity (Supplementary Fig. 4f). Similarly, epiblast-specific deletion of Cdk1 or the cytokinesis regulator Ect2 both resulted in a severe phenotype with developmental arrest (Fig. 6e, Supplementary Fig. 4h, i). In RhoAΔepiblast, as well as CDK1Δepiblast and ECT2Δepiblast mutants, mesoderm identity was specified, as shown by the expression of the PS and nascent mesoderm markers Brachyury (T), and Eomesodermin, (Fig. 5d), and the mesenchymal marker Vimentin (Fig. 5e).Fig. 6: Non-apical mitoses are positively correlated with basal extrusion at the primitive streak.a Scatterplot of the mean percentage of nuclei basal positioning in tested scenarios for epiblast and PS cells. Error bars: SD. b Legend for tested scenarios. Squares correspond to the region of the epiblast and circles to the PS. c Scatter plot of the percentage of non-apical mitoses and mean rate of nuclei basal positioning in tested scenarios for the epiblast and PS. d Scatterplot of the mean rate of basal extrusion and basal positioning in tested scenarios for the epiblast and the PS. e Scatter plot of the mean rate of basal extrusion and the percentage of non-apical mitosis in tested scenarios for the epiblast and PS. For each condition, 100 simulations with 19 PS cells surrounded by 40 control cells.While these data underscore the critical role for RhoA, CDK1 and ECT2 in maintaining epiblast integrity and promoting cell migration, they suggest that changes in actomyosin contractility and cell cycle parameters are not sufficient to trigger or block EMT progression.Non-apical mitosis frequency correlates with delaminationIn view of the results from pharmacological screening and in silico modeling, we conclude that relaxation of cytoplasmic actomyosin, a faster cell cycle, and a shortening of the G2 phase collectively contribute to the increased incidence of non-apical mitoses at the PS during gastrulation. Given that in vivo time-lapse imaging showed that cells undergoing non-apical mitosis primarily gave rise to basally delaminating mesodermal cells13, we wondered if those mechanisms might also be upstream of basal extrusion.Interestingly, computational modeling of EMT in a tall pseudostratified epithelium using the same model indicated that the position of nuclei correlates with the directionality of extrusion21. In particular, basal positioning (defined as more basal than the average position of nuclei) specifically favors basal extrusion. Therefore, some of the experimental and simulated experimental conditions we tested might have a broader impact on nuclei distribution in the tissue than the position of mitoses. However, in a comparatively flat epithelium like the gastrulating epiblast it is unclear whether basal positioning has any importance for delamination given that all cells have their nuclei relatively close to the basement membrane at any time point.To assess this, we plotted the rate of basal positioning of nuclei in epiblast and PS for all simulated conditions (Fig. 6a, b). There was a high rate of basal positioning in control PS cells. Preventing detachment from the basal line or enforcing normal actomyosin contractility (Fig. 6a, PS experimental conditions 2 to 5) prevented basal positioning of nuclei in PS cells. By contrast none of the treatments affecting contractility or cell cycle parameters were able to promote basal positioning in epiblast cells, where only a forced detachment from the basal line was sufficient to generate a high rate of basal positioning (Fig. 6a, epiblast experimental condition 2).Next, we wondered if overall, basal positioning of nuclei and the occurrence of non-apical mitoses were correlated (Fig. 6c). This revealed two different patterns: in epiblast cells, non-apical mitoses occurred in the absence of basal positioning unless loss of basal attachment was enforced (grey square with black outline). In PS cells both parameters were strongly correlated unless we imposed a relaxation of the springs throughout the tissue, which completely canceled basal positioning but only lowered the rate of non-apical mitoses (Fig. 6c, circles in hot colors). These data indicate that non-apical mitoses and overall basal positioning of nuclei can be uncoupled. This means that non-apical mitoses can occur even in absence of a global redistribution of nuclei in the tissue.Next, we looked at a putative correlation between basal positioning and basal extrusion (Fig. 6d) and found a similar relationship: there was a strong positive correlation in PS cells unless contractility was affected (circles in hot colors), and no correlation in the epiblast unless epiblast cells were forced to detach from the basal line (grey square with a black outline). Finally, we looked at the relationship between non-apical mitoses and basal extrusion. Here, the data were not split in two groups, and we found an overall positive correlation between the two factors (Pearson r = 0.59, p = 0.0157).These analyses revealed that non-apical mitoses can specifically occur in absence of overall tissue disorganization leading to massive basal positioning. In addition, non-apical mitoses are positively correlated with basal extrusion. However, it is important to note that extrusion itself only occurs significantly if cells detach from the basal line or if all springs or apical springs are relaxed. This means that actual delamination requires a condition in which epithelial integrity is impaired.Proliferation favors delamination in epiblast spheroidsTo further assess the putative relationship between the frequency/position of mitoses and exit from PS, we turned to a stem-cell-based epiblast model44. This allows for the study of specific events in a system where the chemical environment is controlled, and internal and external mechanical and biological constraints are reduced. Here, we used spheroids, three-dimensional cultures of self-renewing pluripotent epithelial cells (3D EpiSCs) that recapitulate the main transcriptional and architectural features of the gastrulating epiblast9. In this system, cells organize into a sphere lined by a short pseudostratified epithelium, from which cells spontaneously delaminate basally, then turn on the expression of Brachyury. As a large number of spheroids can easily be monitored overtime, they provide a unique opportunity to quantify delamination, which is very challenging in embryos.First, we checked if non-apical mitoses were occurring in mature epiblast spheroids. The spheroids were stained for Podocalyxin-1 (lumen), Phospho-Histone H3 (mitosis), Brachyury (PS), F-actin (phalloidin), and DAPI (nuclei) (Supplementary Fig. 5a). Well-organized spheroids presenting a pseudo-stratified epithelium around a Podocalyxin-1 positive lumen were selected for quantification. Non-apical mitoses accounted for approximately 20% of all mitoses (Supplementary Fig. 5b). Interestingly, out of all Brachyury-positive mitoses, the percentage of non-apical mitoses (64%) was higher than that of apical mitoses (34%) (Supplementary Fig. 5c), indicating that non-apical mitoses are associated with PS-like features in this model as well.3D EpiSCs were cultured with either control vehicle, 0.2 µM of Wee1 inhibitor or 30 µM of ROCK inhibitor, for 4 hours (Fig. 7a, a’). We analyzed the mitotic and non-apical mitotic indexes, as well as the proportion of epithelial and delaminated cells (Fig. 7b–e’). While Wee1 inhibition increased both the mitotic and non-apical mitotic indexes (Fig. 7b, c), ROCK inhibition did not affect the mitotic index, but raised the non-apical mitotic index (Fig. 7b’, c’). These data are consistent with the effects of the inhibitors on embryos. Furthermore, in 3D EpiSCs exposed to both drugs, the proportion of epithelial cells (Fig. 7d, d’) decreased to the benefit of cells that had exited the epithelium (“delaminated” cells, Fig. 7e, e’), particularly for spheroids treated with ROCK inhibitor (Fig. 7d’, e’). To explore this further, we performed live imaging of 3D EpiSCs spheroids expressing LifeAct-GFP, cultured in presence of a nuclear dye and either control vehicle or 0.2 µM of Wee1 inhibitor. This dynamic view illustrated the correlation between cell division and cell delamination (Supplementary Movie 4). In addition, while apical mitoses mostly gave rise to epithelial cells (Supplementary Fig. 5d; Supplementary Movie 5), non-apical mitosis could also lead to the direct exclusion of a daughter cell after cytokinesis (Supplementary Fig. 5e, f; Supplementary Movie 6). This is in line with what is observed in vivo at the PS13. These experiments show that increasing the mitotic and non-apical mitotic indexes favors the delamination of cells from the epithelium in 3D EpiSCs and suggests a similar functional role of proliferation in the mouse embryo at gastrulation.Fig. 7: Increased occurrence of non-apical mitoses in 3D EpiSCs promotes cell delamination.3D EpiSCs cultured 4 hours with vehicle, 0.2 µM of Wee1 inhibitor (a) or 30 µM of ROCK inhibitor (a’), stained for F-actin (phalloidin, grey) and mitosis (Phh3, magenta). White arrowheads indicate non-apical mitosis and red arrowheads indicate delaminated cells. Scale bars: 10 µm. b—c’ Violin plots representing the mitotic (b, b’) and non-apical mitotic index (c, c’) (%) in epithelial cells of 3D EpiSCs cultured 4 hours with control vehicle, 0.2 µM of Wee1 inhibitor (b, c) or 30 µM of ROCK inhibitor (b’, c’). Inner doted lines indicate third quartile, median, and first quartile. d-e’) Boxplots representing the ratio (%) of epithelial (d, d’) and delaminated (e, e’) cells in 3D EpiSCs cultured 4 hours with control vehicle, 0.2 µM of Wee1 inhibitor (d, e) or 30 µM of ROCK inhibitor (d, e’). Boxplots span from the lower quartile to the upper quartile, with the inner line representing the median; whiskers extend to the minimum and maximum values. Normality was assessed using a Shapiro–Wilk test followed by a two-tailed Mann-Whitney or unpaired t-test. Wee 1 inhibitor incubations: control n = 29, Wee1 inhibitor 0.2 µM n = 37. ROCK inhibitor incubations: control n = 36, ROCK inhibitor 30 µM n = 46.
