Note: This response was posted by the corresponding author to Review Commons. The content has not been altered except for formatting.
Learn more at Review Commons
Reply to the reviewers
We appreciate the time and effort the reviewers have invested in providing constructive feedback on our manuscript. Below, we’ve detailed additional work, corrections, and improvements that we will complete during the revision process.
Reviewer #1 (Evidence, reproducibility and clarity (Required)):
Summary
Folding is a major morphogenetic process that shapes tissues and organs in three dimensions. The mechanisms underlying tissue folding have been extensively explored and are often driven by actomyosin-based apical constriction. Here, the authors describe changes in cell geometry and mechanics during mouse neural tube formation. They build on quantitative fixed imaging and live junction ablation to extract cell geometry and junctional tension. These analyses are performed at different developmental stages and in both male and female embryos to propose a mechanical mechanism for neural tube elevation in the brain.
Major comments
The authors report quantitative data on cell geometry and junctional tension inferred from laser ablation. Overall, there are numerous statements that require stronger support from the experimental data. To substantiate several of their claims, the authors need to provide a larger number of data points-or at least comparable numbers across experimental conditions-for the tension measurements. Additional statistical analyses are required throughout to support the conclusions.
Figure 1
- Does the projection algorithm account for tissue curvature when computing cell geometrical parameters such as area and anisotropy?
At present, our projection algorithm does not correct for tissue curvature. Curvature in the tissue can make larger cells appear smaller in projections, skew the angle of cell orientations, and change aspect ratios. The largest curvature in the midbrain neural tube samples that we analyze is found in the transition region from the midline and lateral regions (~10-30% of tissue width) of 5 ss and 8ss embryos. The regions at the midline and more laterally are relatively flat. Therefore, distortion from curvature will not dramatically alter our key conclusions. We will apply a curvature correction using existing tools (Herbert S., et al (2021) BMC Biology) to sample images and determine if there are substantial differences in curvature-sensitive cells shape metrics. These will be included in a supplement to Figure 1. If there is a significant difference, we will expand the correction to all images that we analyze and update our analysis.
The authors should provide information on the accuracy and reliability of the cell segmentation.
We can provide a supplement to Figure 1 to demonstrate the accuracy of the segmentation. We have used F-Actin to segment cells in our images, which is enriched along the cell junctions but can also form medial cables that cross the cell surface. Junctional actomyosin is notably brighter than medial cables, and segmentation with our trained CellPose model is consistently able to distinguish the junctions. We also checked segmentation and performed manual corrections to ensure accuracy. To demonstrate this for our readers, we will prepare samples stained with both F-actin and ZO-1, a tight junction component that is localized to cell junctions. We will then segment the image twice in CellPose, once using the F-actin signal and once using the ZO-1 signal. The resulting cell outlines will then be digitally superimposed to show how much the signals overlap, and we will plot out the cell frequency as a function of area to determine if F-actin segmentations can segment with the same fidelity as ZO-1. Recent work by a co-author has shown excellent corroboration of neuroepithelial apical cell areas segmented using F-Actin and ZO-1 (Ampartzidis I., et al. eLife 2026). We are confident that our data will show a similar result.
The authors indicate that the rate of apical constriction differs between male and female embryos. However, apical sizes differ only at specific positions along the ML axis (Fig. 2H, I).
In Figure 2H, we show that at 5 ss males have larger apical areas than females at the midline, adjacent lateral cells, and at the surface ectoderm-neural epithelium border. By 8ss (Figure 2I), cells at the midline are smaller in males than females, while cells in more lateral regions are now equivalent between sexes. This change in apical area over time suggests that males have faster rates of constriction than females at the midline and adjacent lateral region where male cells become smaller or equivalent in size to female cells, respectively. We will perform statistical analysis (see comment #4) to determine if there are regions with significant differences in rate and amend our language to clarify that these differences are region specific as appropriate.
The authors should provide statistical analyses for the rates shown in Fig. 2J. Are these rates significantly different between males and females, and between medial and lateral regions?
Currently we calculate our rates using the difference in population averages of apical area at each stage shown in Figure 2H and 2I for each sex, and dividing by the number of somite stages, 3. As a result, there is only one rate value at each midline-lateral bin for each sex which is not amenable to statistical analysis. To correct this, we will calculate rates by subtracting the average apical area of each embryo at 8 ss from the population average of embryos at 5 ss. This will create 5 rates for both females and males at each 10% midline-lateral bin. We plan to perform a two-way ANOVA to determine if there are statistical differences in rates between males and females at each bin position and between medial and lateral regions. We will also add a section describing these calculations to the “Statistical Analysis” portion of the methods.
Please clearly state the main novelty of this study relative to the work published by Brooks et al.
Our study builds on the work of Brooks ER, et al. (2020) eLife. Brooks demonstrates that cells in a region of the lateral neural folds undergo apical constriction (Figure 1) and that cells at the midline do not (Figure 2). We expand and improve upon this work in the following ways:
- A) As required by our funding sources at the NIH (NOT-OD-15-102) we have collected, analyzed, and reported on sex as a biological variable of interest. In doing so, we have shown that there are clear sex differences in apical area in the neural tube that were not previously shown. We also show that there is apical constriction within the neural tube midline in a sex dependent manner. Brooks et al do not address sex in their work.
- B) We have provided more complete and spatially precise information on midline-lateral patterns of apical area and apical constriction. To show changes in apical area of lateral cells, Brooks selects a 100 x 100 µm region of interest in the midbrain (Figure 1E-F, Figure 2A) but does not specify the midline-lateral or rostral-caudal location of this region of interest or standardize it between embryos of different ages and dimensions. In our study, we’ve standardized our measurements to a 100 µm wide band across the midbrain adjacent to the midbrain/hindbrain boundary (Figure 2A-C). We also standardize positions as a percent distance from midline to account for differences in width between embryos and ages. This allows us to consistently compare similar populations of cells along the midline-lateral axis and determine changes in apical area over time.
- C) We connect patterns of apical area and constriction to F-actin and Myosin-IIB density. Though Brooks et al report some analysis of F-actin in lateral cells (Figure 6), they do not analyze the midline cells or explore the relationship between cell shape and actomyosin.
- D) Finally, we tested the mechanical properties of the tissue through laser ablation in living mouse embryos. From these ablations we’ve found that tension at the midline is less than in more lateral regions. Work in the neural tubes of frog (Haigo S., et al. (2003) Current Biology, Baldwin AT., et al. (2022) eLife, Matsuda M., et al. (2023) Nature Communications) and chicken (Kinoshita N., (2008) * Cell, Nishimura T., et al. (2012) Cell) embryos has conclusively shown that enriched midline actomyosin promotes apical contractility and drives hinge formation. It was therefore largely believed that a similar contractile hinge was employed in mammals (Copp AJ. and Green NDE. (2010) J. Pathol, Nikolopoulo E., et al. (2017) Development). Collectively, our work is the first to demonstrate that such a contractile hinge is not present in the mammalian brain neural tube.
Figure 3*
The authors need to provide statistical support for the claim that large midline cells exhibit reduced F-actin and Myosin IIB levels.
We will conduct a two-way ANOVA to determine if there are statistical differences in F-actin and Myosin IIB density at the midline and more lateral regions in both males and females. We will update our language in the text and plots as appropriate from these results.
F-actin and Myosin IIB intensities should be plotted as a function of cell area to support the proposed anticorrelation between apical area and actomyosin levels.
We will make plots of cell areas vs. F-actin or Myosin IIB density for cells in each embryo. We will then fit a line to determine the R value for each embryo to determine if there is a negative correlation between cell area and actomyosin intensity. We will also adjust our language in the text as appropriate based on the results of these tests.
Statistical analyses are missing to substantiate the increase in F-actin levels between stages ss5 and ss8.
We will perform an F-test to determine homogeneity of variance between F-actin at 5 ss and 8 ss followed by the appropriate t-test to determine if there is a statistical increase in F-actin over time. We will also amend our language in the text to reflect the results of this test.
Figure S3 should be supported by plots showing Myosin II and F-actin intensity as a function of position along the ML axis, together with appropriate statistics.
In Figure 3A-D, we show representative images of F-Actin and Myosin IIB density in female embryos. These are plotted as the purple lines in Figure 3 E-H. Figure 3 Supplement 1 shows representative images of F-actin and Myosin IIB density in male embryos. These are plotted as the green lines in Figure 3 E-H. We will add a line in the caption of Figure 3 Supplement 1 indicating that these samples are represented and plotted in Figure 3. We also noted a typo in the respective captions, incorrectly indicating male or females were shown in the figure. We will correct these typos as well. Additionally, we will perform the statistical tests indicated under comment #6.
Figure 4
The authors state that lateral tension in male embryos is not different from midline tension, yet the number of data points is much lower than in females. To support this claim, the number of ablations should be comparable across sexes.
As part of this study we performed 270 ablations in the neural tubes of 83 mouse embryos: an exceptional scale of ablations that is the first of its kind in early embryos. We conducted our initial recoil velocity analysis blinded to information on sex. Male embryos were statistically underrepresented in our data set because male embryos develop faster than their female littermates (Seller MJ. and Perkins-Cole KJ. (1987) J. Reprod. Fert.). As such, the neural folds of male embryos were too elevated to ablate. At present we do not have the resources or justification to perform laser ablations on additional animals to obtain the number of male embryos needed to supplement the already exceptionally large data set. We will instead perform a power analysis to determine if: 1) we have a sample size large enough to detect a biologically-meaningful difference with suitable power, 2) the sample size required to detect the observed difference is so large that the difference would not be biologically meaningful, or 3) we do not have a sample size large enough to detect a difference confidently. With the results of this analysis, we will amend our language in the text to reflect the most accurate claims that can be made.
Is lateral tension different between males and females?
In Figure 4G we show that females have statistically different tension between the lateral and midline regions, while males do not. However, we do not test if the lateral or midline tension is different between females and males. We will perform an F-test and t-test to determine if there are statistical differences between males and females in this region.
Similarly, the data in Fig. S4 used to claim no change in tension over time are not supported by sufficient data points.
As discussed in comment #10, the scale of ablations is already substantial, and the initial recoil velocities were analyzed blinded to information on embryo age. We will calculate a best fit line for these plots to demonstrate if there is a trend in recoil velocity over time. We will then adjust our language in the text as appropriate with this added information.
Would the medial and lateral tensions reported in Fig. 4G remain unchanged if the authors perform statistical analyses on 10-15 ablations per condition?
We do not have a justification for removal or exclusion of any of the laser ablations analyzed in this study. We will instead perform a power analysis, as indicated in comment # 10, and adjust the language in the text as appropriate given the results of that analysis.
Figure 5
The number of data points in Fig. 5J and L is insufficient to support claims of no difference. The only detectable difference arises in the comparison with much higher sample size (Fig. 5L, ML vs RC).
In Figure 5J we disaggregate ablations performed at the midline by directionality (midline-lateral or rostral-caudal). We were unable to detect a statistically significant difference based on the direction of initial recoil velocity in either sex, though N’s for all categories are comparable. As discussed in comments #10 and #12, the scale of ablations conducted in this study is uniquely substantial. We will perform a power analysis for our anisotropy measurements in the lateral region of the tissue to determine if we have a sample size large enough to have detect a biologically-relevant difference with high confidence or if the sample size required to detect the observed difference is so large that the difference would not be biologically meaningful. Given the results of this analysis, we will amend our language in the text to reflect the most accurate claims that can be made.
The authors conclude that males have higher ML tension than RC tension, but given the limited data this conclusion should be amended to "no detectable difference."
In Figure 5L, we disaggregate ablations performed in the lateral regions, by directionality (midline-lateral or rostral-caudal). We find a statistical difference in the directionality of initial recoil velocity in females. In males, though we can observe a difference in the initial recoil velocity means, we are unable to detect a statistical difference, likely due to the smaller male sample size. As discussed in comments #10 and #12, the scale of ablations conducted in this study is uniquely substantial and was conducted blinded to embryo sex. Given that males develop faster than their female littermates (Seller MJ. and Perkins-Cole KJ. (1987) J. Reprod. Fert.) we were unable to obtain more males in our data set. We will perform a power analysis for our anisotropy measurements in the lateral region of the tissue to determine if: 1) we have a sample size large enough to detect a biologically-meaningful difference with suitable power, 2) the sample size required to detect the observed difference is so large that the difference would not be biologically meaningful, or 3) we do not have a sample size large enough to detect a difference confidently. With the results of this analysis, we will amend our language in the text to reflect the most accurate claims that can be made.
Code availability
The authors should provide access to the code used to generate the projections.
We are committed to ensuring open access to all code used as part of this study, including components of the projection workflow, data analysis, and figure creation. We are in the process of assembling a GitHub repository containing these files as well as documentation to allow for use by other members of the research community and public. We will publicly publish this documentary upon completion of the repository or at time of publication, whichever comes first.
Reviewer #1 (Significance (Required)):
The authors propose a mechanical model for neural tube elevation based on analyses of cell geometry and tension at two developmental stages. The reported differences in cell geometry or actomyosin levels do not appear to explain the differences in geometry or tension suggested between male and female embryos. This raises questions about the relationship between these measurements and their relevance for understanding the mechanisms of neural tube elevation.
If the major concerns outlined above are rigorously addressed, the manuscript will offer a valuable descriptive characterization of neural tube cell geometry and mechanical stress during morphogenesis. Such datasets could form a foundation for future studies investigating the mechanisms driving neural tube elevation.
Reviewer #2 (Evidence, reproducibility and clarity (Required)):
The manuscript investigates the role of apical constriction and actomyosin organization in shaping the mouse brain neural epithelium during neural tube elevation, with particular emphasis on sex-specific differences. The authors develop an imaging and analysis pipeline to reconstruct the apical surface of the neural plate in three dimensions and perform quantitative measurements of apical cell area, actin, and myosin IIB distributions. Targeted laser ablation experiments are used to infer regional tissue tension.
The main findings can be summarized as follows. First, the authors identify a mediolateral gradient in apical cell area, with larger cells at the midline and smaller cells on the lateral neural folds, which inversely correlates with actomyosin density. Laser ablation experiments suggest that apical tension is lower and isotropic at the midline, whereas it is higher and anisotropic on the lateral folds, particularly in females. Second, sex-dependent differences in apical cell area, constriction rates, and actomyosin levels are reported at early somite stages, preceding previously described sex biases in neural tube defects.
The experimental work is technically solid, and the imaging and quantification pipeline represents a useful advance for analyzing large, curved epithelial surfaces. However, the study feels incomplete in its current form. Despite addressing neural tube elevation, the manuscript does not provide a comprehensive analysis of the folding process itself. Key aspects such as three-dimensional tissue morphology, curvature evolution, or global shape changes of the neural folds are not quantified. In addition, other potentially relevant cellular behaviors, such as proliferation, cell rearrangements, or contributions from neighboring tissues, are not examined, nor are they compared systematically between sexes.
Conceptually, the study focuses narrowly on correlations between apical cell area, actomyosin density, and inferred tension. While these measurements are carefully performed, the relationship between differential actomyosin contractility and three-dimensional tissue folding remains largely descriptive. No mechanical model or simulation framework is provided to link changes in actomyosin organization and cell shape to the emergence of neural folds and hinge formation. As a result, it is difficult to assess whether the measured differences in tension (on the order of ~40%) are sufficient to account for the proposed mechanical behavior of the tissue.
The central hypothesis advanced by the authors is that a relatively "soft" midline, flanked by stiffer, tension-bearing lateral folds, facilitates hinge formation during brain neurulation. However, this hypothesis is not directly tested by perturbation. For example, experimentally increasing contractility or stiffness at the midline (e.g., via optogenetic activation of apical constriction machinery) would provide a more direct test of causality. As it stands, the data demonstrate correlation rather than necessity or sufficiency.
Relatedly, alternative interpretations are not fully addressed. Large apical cell areas and low actomyosin levels at the midline could arise as a consequence of tissue geometry, contact with underlying structures such as the notochord, or extrinsic mechanical constraints, rather than being the primary cause of hinge formation. Similarly, anisotropic stresses generated at the tissue or embryo scale could align cells and actomyosin cables, producing the observed patterns without requiring locally specified apical tension differences as the initiating mechanism. The manuscript does not clearly distinguish whether apical tension asymmetries are a driver of folding or an emergent outcome of folding dynamics.
Finally, while the identification of sex differences is intriguing, it remains unclear what mechanistic insight is gained beyond establishing that such differences exist. The functional consequences of these differences for neural tube closure, robustness, or failure are not explored, nor is it clear how they integrate into the proposed lateral tension model.
In summary, this study provides high-quality measurements of apical cell geometry, actomyosin organization, and inferred tension in the mouse neural epithelium. However, the lack of direct perturbations, mechanical modeling, and quantitative analysis of three-dimensional tissue deformation limits the strength of the mechanistic conclusions. Addressing these gaps would substantially strengthen the manuscript and clarify the causal role of apical tension patterns in neural fold formation.
__ __The reviewer makes an excellent point, that direct perturbation of the system would enable us to test our hypothesis and inform whether the reduced contractility at the midline is essential for neural tube elevation. However, at present the technology needed to conduct an optogenetic experiment like that described by the reviewer does not exist. As with the laser ablations, an optogenetic experiment requires access to live and healthy embryos. Currently, mouse embryos can be cultured for several days in roller culture, where they are continuously rotated, or for several hours in static culture (Aguilera-Castrejon A. and Hanna JH. (2021) J. Vis. Exp.). Both techniques require that the yolk and amniotic sacs remain intact around the embryo. To access the apical surface of the brain neural tube for imaging, both sacs must be breached, after which the embryo has about 30 minutes before it begins to exhibit altered cellular morphology and tissue integrity and ultimate embryo death.
The neural tube elevates over several hours and closes fully after more than a day (Jacobson AG. and Tam PPL. (1982) The Anatomical Record). Even if we did acquire mice expressing photoactivatable constructs, the support membranes of the embryos would need to be breached to activate protein interactions. The embryos would die before any meaningful progress in neural tube elevation could be evaluated. Conducting an experiment like this would greatly advance our understanding of the system, and we hope that the needed technologies are developed to enable future work of this nature. The Galea lab previously purchased a photo-activatable Cre line, but was unable to induce deletion of a protein of interest using this allele before closure of the neural tube was completed (and the blue light needed to activate the cre was photo-toxic).
At present, there is some experimental evidence to suggest that lack of apical constriction at the midline if important for proper neural tube closure. Brooks ER, et al. (2020) eLife shows that a truncated Ift122 mutant, leads to abnormal constriction of the midline cells but does not disrupt lateral cell apical constriction, leading to a failure in brain neural tube closure in these embryos. Ift122 regulates trafficking and signaling proteins in cilia, which in turn regulates Sonic hedgehog signaling which Brooks ER, et al. also demonstrates regulates apical constriction. While this disruption is clearly multifaceted and nuanced, it provides some genetic support for the lack of apical constriction at the midline being important for neural tube closure.
Major Comments
Figure quality. Figure 1 contains very low-resolution images, which makes it difficult to evaluate the segmentation quality and tissue morphology. Higher-resolution versions should be provided.
In Figure 1, we outline the conceptual strategy and approach used to create and analyze shell projections of the curved neural tube. As much of our analysis builds from segmentation of cells in the projections, being able to assess segmentation quality from high resolution images is critical to evaluating the quality of the data shown. As discussed in comment #2, we will create a supplement to Figure 1 to demonstrate the accuracy of the segmentation. This will include high resolution images of both the label used to segment and the resulting segmentation, with corresponding overlays.
Cell segmentation strategy and validation. The authors segment cell areas using Myosin II and F-actin signals. This approach may introduce inaccuracies, as actomyosin cables can traverse the apical surface of individual cells and do not always coincide with cell boundaries. Segmentation based on junctional markers such as ZO-1 may be more appropriate. At minimum, the authors should provide a quantitative validation of segmentation accuracy, for example by overlaying segmentation results on raw images together with a nuclear marker (e.g., DAPI or H2B-GFP), to demonstrate that the number of segmented cells corresponds to the number of nuclei.
We will provide a supplement to Figure 1 to demonstrate the accuracy of the segmentation. We have used F-Actin to segment cells in our images. F-actin is enriched along junctions but cells can also have medial pools and F-actin cables, which might lead to errors. Though we understand the reviewer’s logic in asking to align segmentations with marked nuclei, the morphology of the neural epithelium makes this approach infeasible. The neural epithelium is pseudostratified, and nuclear position varies along the apical-basal axis depending on the cell cycle phase of each cell. As a result, an apical shell projection of nuclei would not capture all nuclei and a maximum intensity projection in Z of all nuclei would be uninterpretable as there would be substantial XY overlap between nuclei. Instead, we will create a supplement to Figure 1 to demonstrate the accuracy of the segmentation as discussed in comment #2. We will segment samples stained for both F-Actin and junctional markers like ZO-1. We will then create overlays of the resulting cell outlines and a cell area frequency plot for both segmentations to evaluate if F-actin based segmentation deviates from tight junction-based segmentation.
Lack of cross-sectional views of neural tube morphology. The manuscript would benefit from the inclusion of cross-sectional images of the neural tissue at different developmental stages. This would serve two purposes: (i) to demonstrate that the authors have a comprehensive understanding of the full three-dimensional folding process during neural tube closure, including medial and lateral hinge formation, and (ii) to allow readers to visualize the tissue geometry corresponding to the analyzed projection datasets (e.g., at 5 ss and 8 ss).
A key component of our model states that the changes in cell-level morphology and features correspond to changes in tissue level morphology (Figure 6). Specifically, that lateral apical constriction coincides with the flattening and elevation of the dorsal bulges on the lateral neural folds. We agree that it is beneficial to include additional visuals of tissue morphology. We plan to add an additional figure at the start of manuscript that details both the dorsal and relevant cross-sectional views of the somite stages analyzed. These visuals will take the form of graphical illustrations along with 3D confocal microscopy images and optical reconstructions of samples.
Sex-specific differences in overall neural plate morphology. The authors report that at 5 ss, males consistently have larger apical cell areas than females. It is unclear whether this difference reflects a global difference in neural plate morphology. Showing representative images of female and male neural plates would help readers directly assess whether there are overt morphological differences beyond those revealed by quantitative analysis.
If one sex has larger cells than the other, it would be reasonable to expect that the neural folds may be wider as well. In Figure 2B-C, we show representative images of male embryos at 5 and 8 ss. As part of the additions we indicated in comment #19, we will also include dorsal and cross-sectional views of both male and female embryos at the stages analyzed. If there is a difference in tissue morphology between sexes, we will also quantify these differences in tissue size, curvature, etc.
Cell number analysis. The authors state, based on prior literature, that cell numbers do not change between 5 and 8 ss. Given that the tissue is already segmented in the current study, this claim should be directly verified using the authors' own data. This analysis should be straightforward and would strengthen the conclusions.
We agree and will determine the number of cells analyzed for each embryo to test if there are changes in cell numbers at different stages and between sexes, along with appropriate statistical tests.
Relation between tissue curvature and cellular properties. It would be highly informative to extract the three-dimensional morphology of the neural plate, in particular its curvature, and examine how curvature correlates with two-dimensional cell anisotropy, apical area, and F-actin/Myosin intensity. For example, at 8 ss the authors report a U-shaped dependence of cell area along the mediolateral axis. How does this pattern relate to local tissue curvature?
We agree with this assessment and will create optical reslices in the midbrain adjacent to but excluding the midbrain hindbrain boundary. We will then divide the apical surface into 10% bins and fit a circle to the apical surface of the neural epithelium in each to calculate the local radius of curvature, which is the reciprocal of curvature for the surface. We can then correlate these values with two-dimensional cell shape and actomyosin density metrics.
Visualization of sex differences in medial actin levels (Figure 3). In Figure 3, the reported female-male difference in medial actin levels would benefit from visualization of the raw data. A zoomed-in inset of the midline region, shown separately for females and males, would help substantiate this claim.
In Figure 3, we demonstrate patterns of the whole-cell apical F-actin (Fig. 3A, B) and Myosin IIB (Fig. 3C, D) density. We find that there is no difference in F-Actin density between males and females (Fig. 3E, F), but a significant difference in midline Myosin IIB density at 5 ss that is mostly absent by 8 ss (Fig. 3G, H). We currently provide representative images for female and male myosin IIB expression across the midline-lateral axis in Figure 3C, D, and Figure 3-Supplement 1C and D. We can provide a close-up image of Myosin IIB in the midline region for both sexes as part of Figure 3, with additional annotations on existing representative image to indicate their origin.
Typographical error. Line 143: please correct "cell are" to "cell area".
We thank the reviewer for pointing out this error and will correct this typo and perform additional editing to correct any other typos present in the manuscript.
Quantitative correlation analysis between cell area and actomyosin. The authors qualitatively discuss the relationship between cell area dynamics and actomyosin levels. It would strengthen the analysis to directly compute and report correlations between these variables, and to explicitly test whether actin and myosin levels are anti-correlated with apical cell area.
As discussed in comment #6, we will plot cell area vs. F-actin or Myosin IIB density for each embryo and fit a line to calculate their correlation coefficient. From there, we will determine if there is a negative correlation between cell area and actomyosin intensity.
Interpretation of anti-correlation and contractile hinge mechanism. In lines 143-157, the authors state that the observed anti-correlation between actomyosin and cell area argues against a contractile hinge mechanism. However, this anti-correlation could also suggest that apical cell area is determined by local mechanical or geometric constraints rather than by local actomyosin contractility. The authors should clarify and discuss this alternative interpretation.
Within the neural epithelium of mice and other vertebrates, F-actin and myosin-IIB are enriched on the apical surface relative to other regions of the cell (Sadler TW, et al. (1982) Science, Matsuda M., et al. (2023), Nat. Communication, Röper, K. (2013) BioArchitecture). This poises the actomyosin network to be able to selectively constrict the apical surface relative to the basal side of the cells. Apical constriction is observed to actively facilitate the formation of hinges in folding tissues (Chanet S, et al. (2017), Nature Communications, Nishimura T., et al. (2012) Cell, Chistrodoulou N., et al. (2015) Cell Reports) in what we term the contractile hinge model of tissue folding. Tissues that employ this model of folding are expected to have small apical areas and apical enrichment of contractile actomyosin at the hinge point during folding. We observe large apical areas, low apical actomyosin density, and low apical tension at the midline hinge of the mouse midbrain neural tube, which are all inconsistent with a contractile hinge mechanism being employed in this tissue folding process. We agree with the reviewer that “cell shape does not always match [acto]myosin contractility levels, because cell shape depends on extrinsic, as well as intrinsic forces” (Line 147-149). We also agree that anticorrelation of actomyosin density and apical cell area does not per se argue against the contractile hinge model and will amend our language to be clearer. We will also further elaborate on potential extrinsic factors that may lead to the observed cell behaviors at the midline in the discussion.
Statistical robustness of laser ablation results (Figure 4). The differences in recoil velocity between regions appear small, with substantial overlap between the distributions. In addition, the sample sizes for lateral versus midline ablations appear unequal (with visibly more data points in the lateral condition). These factors raise concerns about the robustness and statistical significance of the reported differences, which should be addressed more carefully.
In Figure 4E, we show initial recoil velocities binned only by region: lateral vs. midline and report a 3.03 μm/s vs. 2.40 μm/s, or 26% difference between the two regions. We then show in Figure 4G that by considering another relevant variable, sex, we find initial recoil to be 3.15 μm/s vs. 2.30 μm/s, or 37% difference in females and 2.68 μm/s vs. 2.57 μm/s, or 4% difference in males. We go on to show In Figure 5L that within the lateral region that recoils also vary by direction, with a 38% difference. Ultimately the final conclusions that we draw regarding tissue tension that we present in our model are derived from the most finely disaggregated data in Figure 5. Our goal in presenting a stepwise disaggregation of the data was to demonstrate which variables had the greatest impact on the variance within our data set. We agree with the reviewer that a more precise statistical analysis of this data set is warranted that accounts for the complexity and multitude of variables that can influence our conclusions. In addition to the power analysis described in comment #10, we plan to conduct a mixed-effect model analysis of our data that considers factors including sex, age, cut direction, cut region, cut number, and embryos to determine which factors explain the most variance in the population. We will add this analysis as a supplement to Figure 4 alongside a description of the tests performed in the Statistical Analysis section of the methods. We will also adjust our language in the text to clearly state the limitations of the data as presented and qualify conclusions as appropriate.
Speculative statement regarding anisotropic tension in males. Line 278: "We believe that both sexes demonstrate anisotropic tension, given that males have cell aspect ratios and orientations in the lateral neural folds similar to females." This statement is speculative. Either anisotropic tension in males should be directly measured and reported, or this statement should be removed.
As discussed in comment # 15, in Figure 5L, though we can observe a difference in the initial recoil velocity means, we are unable to detect a statistical difference. Ablations were conducted blinded to embryo sex, but fewer male embryos were suitable for ablation because males develop faster than their female littermates (Seller MJ. and Perkins-Cole KJ. (1987) J. Reprod. Fert.). We were therefore unable to obtain more males in our data set. At present we do not have the resources to perform additional laser ablations to supplement the existing data set. We will instead perform a power analysis for our anisotropy measurements in the lateral region of the tissue to determine if: 1) we have a sample size large enough to detect a biologically-meaningful difference with suitable power, 2) the sample size required to detect the observed difference is so large that the difference would not be biologically meaningful, or 3) we do not have a sample size large enough to detect a difference confidently. With the results of this analysis, we will amend our language in the text to reflect the most accurate claims that can be made.
Reviewer #2 (Significance (Required)):
This study provides high-quality measurements of apical cell geometry, actomyosin organization, and inferred tension in the mouse neural epithelium. However, the lack of direct perturbations, mechanical modeling, and quantitative analysis of three-dimensional tissue deformation limits the strength of the mechanistic conclusions. Addressing these gaps would substantially strengthen the manuscript and clarify the causal role of apical tension patterns in neural fold formation. At the end of the day, the authors suggest an hypothesis that is not well support by their data, which is of high quality.
Reviewer #3 (Evidence, reproducibility and clarity (Required)):
Summary
This manuscript by De La O et al addresses a long-standing question of how actomyosin contributes mechanically to cranial neural tube elevation in the mouse, a system in which classical midline contractile hinge models appear insufficient. The authors develop an image-processing and analysis pipeline that enables reconstruction and quantitative analysis of the apical actomyosin network across the large, curved dorsal surface of the mouse brain neuroepithelium. Using this approach, combined with laser ablation-based tension measurements in live embryos, they report a medio-lateral gradient of apical cell area and an inverse gradient of actomyosin density. Contrary to contractile hinge models described in frog, chick and invertebrate systems, they find that the midline exhibits low, isotropic tension, while the lateral neural folds show higher anisotropic apical tension, consistent with their proposal of a "lateral tension" mechanism for neural tube elevation.
The work provides an important reframing of actomyosin function in mammalian cranial neurulation, supported by extensive quantitative imaging and mechanical measurements. The finding that lateral, rather than midline, actomyosin networks dominate tissue tension is compelling and helps reconcile previous observations that midline hinge formation in mouse can proceed despite actomyosin perturbation. The study is technically sophisticated and addresses a biologically important process with clear relevance to neural tube defect etiology. However, several aspects of the statistical treatment, interpretation of laser ablation data, and mechanistic framing require clarification or tempering to fully support the authors' conclusions.
Major comments
Statistical unit and pseudo-replication in cell-based analyses (Figures 2-3)
In Figures 2 and 3, it is unclear whether statistical comparisons were performed at the level of individual cells or embryos. Because cells are nested within embryos, treating cells as independent observations raises concerns about pseudo-replication and inflated statistical significance, particularly for sex-dependent effects. While the color-coded maps are visually compelling, they may overstate confidence in differences between conditions if embryo-to-embryo variability is not explicitly accounted for.
Clarification is needed as to whether statistical testing was performed on embryo-level summary values (e.g., one value per embryo per positional bin), or whether hierarchical or mixed-effects models were used with embryo treated as a random effect. Providing embryo-level summary plots would also help readers assess inter-embryo variability. Addressing this point is important for confidence in both the reported medio-lateral gradients and the sex differences.
We agree with the reviewer that it is inappropriate to calculate statistics based on measurements of individual cells. As indicated under the ‘Statistical Analysis and Figure Assembly’ section of our methods “For fixed images, cell shape and protein intensity analysis (Figure 2H-J, Figure 3E-H, Figure 5E-H), N = 5 embryos for all conditions and n, or the number of cells in each 10% bin, is ≥ 150 cells for each embryo” (Line 556 – 558). The average and SD between embryos are shown in these plots and is calculated at the embryo level, not the cell level. We chose to consolidate this information in the methods section as the same data set is used across the three figures. We will add a line to the figure captions that N values for all experiments can be found in this section of the methods. We will also provide supplementary plots showing the bin averages for each individual embryo, color coded by embryo to show the distribution of the data set.
Interpretation of actomyosin density as a proxy for contractility (Figure 3)
The descriptive correlation between apical cell area and actomyosin density is clear and consistent. However, actomyosin abundance alone does not necessarily equate to force generation, particularly in the absence of measurements of myosin activation state (e.g., pMLC), actomyosin dynamics, or direct perturbations linking actomyosin levels to mechanical output. Although the authors appropriately note that cell shape does not always reflect intrinsic contractility, actomyosin density is nevertheless used to argue against a contractile hinge mechanism.
While the subsequent laser ablation experiments address tissue tension more directly, the mechanistic conclusions drawn from actomyosin density measurements alone would benefit from more careful qualification. Tempering language that equates actomyosin enrichment with contractile output, or explicitly acknowledging these limitations, would strengthen the interpretation.
It is largely believed that apical pools of actomyosin are active and that apical localization of actomyosin is dependent on activation. Shroom3, an actin-binding protein, is localized to the apical adherens junctions in the neural tube (Haigo SL., et al. (2003) Curr. Biol., Hildebrand JD. and Soriano P. (1999) Cell), where it can recruit Rho kinases (ROCKs) that in turn phosphorylate and activate Myosin IIB (Nishimura T. and Takeichi M. (2008) Dev.). Mutations in Shroom3 lead to neural tube close defects and its overexpression in the neural tube can induce apical constriction and increased apical accumulation of Myosin II tube (Haigo SL., et al. (2003) Curr. Biol., Hildebrand JD. (2005) J. Cell Sci.). In the mouse neural tube, Myosin IIB intensity is greater in cells that can apically constrict than in those that cannot constrict (Galea GL., et al. (2021) Nat Commun). Additionally, inhibition of ROCK reduces apical tension, presumably by reduction of activated Myosin II (Butler MB., et al. (2019) J. Cell Sci.). We agree with the reviewer’s assessment that to definitively state that the apical pools of Myosin IIB and F-actin are promoting apical contractility, a demonstration of the phosphorylation state of the Myosin II regulatory light chain (pMLC) or observations/perturbations in live embryos is necessary. We will adjust our language to reflect this limitation. We will also provide information on the relationship between apically localized actomyosin and contractility.
Statistical and biological independence of laser ablation measurements (Figures 4-5)
The Methods indicate that 155 laser ablations were analyzed across 71 embryos, implying that multiple ablations were performed per embryo. It would be helpful to clarify how this hierarchical data structure was handled statistically. Specifically, were recoil velocities averaged per embryo, paired with embryos for ML vs. RC comparisons, or analyzed using hierarchical/mixed-effects models?
Our laser ablation data set captures variables including embryo sex, age, cut location, cut direction, and cut number. Therefore, we did not feel it appropriate to average recoils within the same embryo as these cuts were intentionally in different regions (lateral vs midline) or in different orientations (i.e. a rostral-caudal cut and midline-lateral cut on opposite lateral folds), which our analysis has shown would lead to averaging out potential differences. Ablations were far apart from each other, and we had checked that ablation order did not predict changes in recoil. However, we agree with the reviewer that a more precise statistical analysis of this data set is warranted that accounts for the complexity of variables potentially influencing initial recoil velocities. As discussed in comment #27, we plan to conduct a mixed-effect model analysis of our data that considers the above and add this analysis as a supplement to figure 4. We will include a description of this in the methods and our language in the text to clearly state the limitations of the data as presented and qualify conclusions as appropriate.
In addition, embryos were subjected up to 5 ablations within a short time window. Because laser ablation disrupts tissue integrity and can induce rapid cytoskeletal remodeling, it is unclear whether later ablations represent independent measurements of the native tension state. Clarification is needed regarding whether the authors tested for effects of ablation order (e.g., first vs. later cuts), ensured sufficient spatial separation between ablation sites, or verified that repeated ablations did not systematically alter recoil measurements. Demonstrating that initial recoil velocity is independent of cut number would substantially strengthen confidence in the mechanical conclusions.
We agree with the reviewer that laser ablations cause disruptions to tissue, and these disruptions can impact the results of additional ablations performed near the site of prior ablations. The average embryo in our data set has three ablations: one on either neural fold and one at the midline, with hundreds of µm distances from each other. In embryos that had more than 3 ablations made far away from each other (additional ablations were performed in the hindbrain rhombomeres, rhombomere boundaries, at the neuroepithelium and surface ectoderm boundary, or at the zipper point, but n numbers of these are insufficient for analysis). We will supplement the methods text describing the laser ablations to clarify this for readers. Additionally, after an ablation, displacement is not detectable further than 3-5 cell lengths away from the cut even after several seconds post ablation. We will provide visual examples of these cuts after ablation to demonstrate this phenomenon. As discussed in comment #27 and #32, we will also perform mixed-effect modeling to determine if cut number impacts observed initial recoil velocities. We will also provide plots demonstrating relevant examples of these comparisons (e.g. sequential lateral cuts made in the same direction).
Interpretation of sex-dependent tension differences (Figures 4-5)
Figure 4 shows a clear lateral-greater-than-midline tension difference in females, whereas this pattern is not detected in males under initial analysis. Later, Figure 5 reveals directional anisotropy in the lateral neural folds of both sexes. As currently framed, this creates some ambiguity regarding whether the proposed lateral tension mechanism is sex-specific, sex-biased in magnitude, or sex-general but masked by directional averaging in males.
Clarifying this distinction, both in figure presentation and in the text, would strengthen the mechanistic interpretation and prevent confusion. In particular, it would be helpful to more clearly explain how directional anisotropy reconciles the apparent absence of regional tension differences in males in Figure 4.
We appreciate the reviewer taking the time to indicate this point of confusion. We ultimately conclude that the lateral tension model of neural tube elevation is agnostic of sex. Though there are nuanced differences in some of the details regarding Myosin IIB density, midline apical constriction and tension anisotropy, we do not believe these differences would fundamentally change the mechanical model used between sexes. With specific regards to masking of the lateral neural fold tension in males, we briefly address this in the discussion: “The averaging of [Rostra-Caudal] and [Midline-Lateral] [Initial Recoil Velocities] likely masked tension differences between the midline hinge and lateral neural folds, creating the false impression that males did not have high tension on the lateral neural folds” (Line 280-282). We will adjust the text in the results and discussion section to clearly indicate that are lateral tension model applies to both sexes, though some differences in specific details exist, and that averaging may have led to the result in Figure 4G.
Causal overreach in mechanistic interpretation of anisotropic tension
While the laser ablation data convincingly demonstrates spatial and directional differences in recoil consistent with patterned mechanical anisotropy, the manuscript frequently treats anisotropic apical tension as a mechanistic explanation for neural tube elevation. The presented experiments do not directly test whether anisotropic apical tension is necessary or sufficient for tissue bending, nor whether isotropic tension at the midline plays a causal role. Initial recoil velocity reflects not only pre-existing tension but also tissue geometry and viscoelastic properties, which may differ between midline and lateral regions.
As such, statements suggesting that anisotropic lateral tension "explains" neural fold elevation should be tempered or reframed. The data strongly support spatial patterning of mechanical properties but do not yet establish causal primacy. Recasting the model as a mechanically consistent framework rather than a definitive mechanism would better align conclusions with the data.
Our lateral tension model proposes that a regionalized difference in tension, with high tension in the lateral neural folds and low tension at the midline, is needed to enable neural tube elevation and ultimate closure. We agree with the reviewer, our work demonstrates that the results of our laser ablation experiments, along with measurements of cell shapes and protein density, are consistent with the lateral tension model that we propose. Our model is also supported by past work that shows that perturbations that disrupt actomyosin contractility leads to defects in brain neural tube elevation and closure but not midline hinge formation. For example, chemical perturbation of actin polymerization with Cytochalasin D (Ybot-Gonzalez P. and Copp AJ. (1999) Dev. Dyn), and genetic perturbations of Shroom 3, which apically localizes actomyosin (Hildebrand JD. And Soriano P. (1999) Cell) or Fhod3, which promotes actin polymerization (Sulistomo HW, et al. (2019) J. Biol. Chem.) all have brain neural tube closure defects but form a midline hinge. However, since we do not directly perturb tension, we have only demonstrated consistency rather than causality or sufficiency. We will adjust and temper our language accordingly in the relevant sections of the results and discussion.
Minor comments
Manuscript length and clarity
The manuscript is longer and more complex than necessary for its central message. Several sections of the Results, particularly methodological validation and somite-stage stratification, could be streamlined.
We agree with the reviewer and will continue editing the manuscript, prioritizing clarity, brevity, and precision of language so that readers are able to quickly understand the key points of the manuscript.
Sex differences section
The section on sex differences is interesting but somewhat tangential. Clarifying whether these findings are intended as mechanistic insight or observational motivation for future work could improve focus.
We intended this section to offer up perspectives that inform and motivate future work to continue to track, analyze, and report on sex difference during development. We will edit this section of the discussion to improve clarity and brevity so that the reader can easily acquire this takeaway. Sex differences in the penetrance of exencephaly is an active area of research and our manuscript provides the first cell-level measurements which will guide the field in disaggregating future analyses by embryo sex.