Skip to main content

Thank you for visiting You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Human neural tube morphogenesis in vitro by geometric constraints


Understanding human organ formation is a scientific challenge with far-reaching medical implications1,2. Three-dimensional stem-cell cultures have provided insights into human cell differentiation3,4. However, current approaches use scaffold-free stem-cell aggregates, which develop non-reproducible tissue shapes and variable cell-fate patterns. This limits their capacity to recapitulate organ formation. Here we present a chip-based culture system that enables self-organization of micropatterned stem cells into precise three-dimensional cell-fate patterns and organ shapes. We use this system to recreate neural tube folding from human stem cells in a dish. Upon neural induction5,6, neural ectoderm folds into a millimetre-long neural tube covered with non-neural ectoderm. Folding occurs at 90% fidelity, and anatomically resembles the developing human neural tube. We find that neural and non-neural ectoderm are necessary and sufficient for folding morphogenesis. We identify two mechanisms drive folding: (1) apical contraction of neural ectoderm, and (2) basal adhesion mediated via extracellular matrix synthesis by non-neural ectoderm. Targeting these two mechanisms using drugs leads to morphological defects similar to neural tube defects. Finally, we show that neural tissue width determines neural tube shape, suggesting that morphology along the anterior–posterior axis depends on neural ectoderm geometry in addition to molecular gradients7. Our approach provides a new route to the study of human organ morphogenesis in health and disease.

This is a preview of subscription content

Access options

Buy article

Get time limited or full article access on ReadCube.


All prices are NET prices.

Fig. 1: A reproducible human stem-cell model of neural-tube morphogenesis.
Fig. 2: Formation of a neural and non-neural ectoderm bilayer.
Fig. 3: Towards modelling neural tube defects.
Fig. 4: Neural plate size determines neural tube shape.

Data availability

Source data are provided with this paper. scRNA-seq data have been deposited and are publicly available in the NCBI Gene Expression Omnibus (GEO; under accession GSE173492Source data are provided with this paper.

Code availability

The Matlab code for the simulation used in this manuscript is deposited in GitHub (


  1. 1.

    Wallingford, J. B., Niswander, L. A., Shaw, G. M. & Finnell, R. H. The continuing challenge of understanding, preventing, and treating neural tube defects. Science 339, 1222002 (2013).

    Article  Google Scholar 

  2. 2.

    Lee, S. & Gleeson, J. G. Closing in on mechanisms of open neural tube defects. Trends Neurosci. 43, 519–532 (2020).

    CAS  Article  Google Scholar 

  3. 3.

    Kim, J., Koo, B.-K. & Knoblich, J. A. Human organoids: model systems for human biology and medicine. Nat. Rev. Mol. Cell Biol. 21, 571–584 (2020).

    CAS  Article  Google Scholar 

  4. 4.

    Quadrato, G. & Arlotta, P. Present and future of modeling human brain development in 3D organoids. Curr. Opin. Cell Biol. 49, 47–52 (2017).

    CAS  Article  Google Scholar 

  5. 5.

    Haremaki, T. et al. Self-organizing neuruloids model developmental aspects of Huntington’s disease in the ectodermal compartment. Nat. Biotechnol. 37, 1198–1208 (2019).

    CAS  Article  Google Scholar 

  6. 6.

    Britton, G., Heemskerk, I., Hodge, R., Qutub, A. A. & Warmflash, A. A novel self-organizing embryonic stem cell system reveals signaling logic underlying the patterning of human ectoderm. Development 146, dev179093 (2019).

    CAS  Article  Google Scholar 

  7. 7.

    Ybot-Gonzalez, P., Cogram, P., Gerrelli, D. & Copp, A. J. Sonic hedgehog and the molecular regulation of mouse neural tube closure. Development 129, 2507–2517 (2002).

    CAS  Article  Google Scholar 

  8. 8.

    Nikolopoulou, E., Galea, G. L., Rolo, A., Greene, N. D. E. & Copp, A. J. Neural tube closure: cellular, molecular and biomechanical mechanisms. Development 144, 552–566 (2017).

    CAS  Article  Google Scholar 

  9. 9.

    Nikolaev, M. et al. Homeostatic mini-intestines through scaffold-guided organoid morphogenesis. Nature 585, 574–578 (2020).

    Article  ADS  Google Scholar 

  10. 10.

    Xue, X. et al. Mechanics-guided embryonic patterning of neuroectoderm tissue from human pluripotent stem cells. Nat. Mater. 17, 633–641 (2018).

    CAS  Article  ADS  Google Scholar 

  11. 11.

    Zheng, Y. et al. Dorsal-ventral patterned neural cyst from human pluripotent stem cells in a neurogenic niche. Sci. Adv. 5, eaax5933 (2019).

    CAS  Article  ADS  Google Scholar 

  12. 12.

    Sahni, G. et al. A micropatterned human‐specific neuroepithelial tissue for modeling gene and drug‐induced neurodevelopmental defects. Adv. Sci. 8, 2001100 (2021).

    CAS  Article  Google Scholar 

  13. 13.

    Developmental Stages in Human Embryos (eds O'Rahilly, R. & Müller, F.) (Carnegie Institute of Washington, 1987).

  14. 14.

    Ray, H. J. & Niswander, L. A. Dynamic behaviors of the non-neural ectoderm during mammalian cranial neural tube closure. Dev. Biol. 416, 279–285 (2016).

    CAS  Article  Google Scholar 

  15. 15.

    Rolo, A. et al. Regulation of cell protrusions by small GTPases during fusion of the neural folds. eLife 5, e13273 (2016).

    Article  Google Scholar 

  16. 16.

    Martins-Green, M. Origin of the dorsal surface of the neural tube by progressive delamination of epidermal ectoderm and neuroepithelium: implications for neurulation and neural tube defects. Development 103, 687–706 (1988).

    CAS  Article  Google Scholar 

  17. 17.

    Molè, M. A. et al. Integrin-mediated focal anchorage drives epithelial zippering during mouse neural tube closure. Dev. Cell 52, 321–334.e6 (2020).

    Article  Google Scholar 

  18. 18.

    Ji, Y. et al. Single cell transcriptomics and developmental trajectories of murine cranial neural crest cell fate determination and cell cycle progression. Preprint at (2021).

  19. 19.

    Vega, F. M., Fruhwirth, G., Ng, T. & Ridley, A. J. RhoA and RhoC have distinct roles in migration and invasion by acting through different targets. J. Cell Biol. 193, 655–665 (2011).

    CAS  Article  Google Scholar 

  20. 20.

    Baldassarre, M. et al. Filamins regulate cell spreading and initiation of cell migration. PLoS ONE 4, e7830 (2009).

    Article  ADS  Google Scholar 

  21. 21.

    Grande-García, A. et al. Caveolin-1 regulates cell polarization and directional migration through Src kinase and Rho GTPases. J. Cell Biol. 177, 683–694 (2007).

    Article  Google Scholar 

  22. 22.

    de Almeida, P. G., Pinheiro, G. G., Nunes, A. M., Gonçalves, A. B. & Thorsteinsdóttir, S. Fibronectin assembly during early embryo development: A versatile communication system between cells and tissues. Dev. Dyn. 245, 520–535 (2016).

    Article  Google Scholar 

  23. 23.

    Haigo, S. L., Hildebrand, J. D., Harland, R. M. & Wallingford, J. B. Shroom induces apical constriction and is required for hingepoint formation during neural tube closure. Curr. Biol. 13, 2125–2137 (2003).

    CAS  Article  Google Scholar 

  24. 24.

    Chen, Z., Kuang, L., Finnell, R. H. & Wang, H. Genetic and functional analysis of SHROOM1–4 in a Chinese neural tube defect cohort. Hum. Genet. 137, 195–202 (2018).

    CAS  Article  Google Scholar 

  25. 25.

    Butler, M. B. et al. Rho kinase-dependent apical constriction counteracts M-phase apical expansion to enable mouse neural tube closure. J. Cell Sci. 132, jcs230300 (2019).

    CAS  Article  Google Scholar 

  26. 26.

    Schoenwolf, G. C. & Smith, J. L. Mechanisms of neurulation: traditional viewpoint and recent advances. Development 109, 243–270 (1990).

    CAS  Article  Google Scholar 

  27. 27.

    Hernandez, I. et al. A farnesyltransferase inhibitor activates lysosomes and reduces tau pathology in mice with tauopathy. Sci. Transl. Med. 11, eaat3005 (2019).

  28. 28.

    Karzbrun, E., Khankhel, A. & Streichan, S. J. Recapitulating neural tube morphogenesis with human pluripotent stem cells. Protoc. Exch. (2021).

  29. 29.

    Li, L. et al. Ectodermal progenitors derived from epiblast stem cells by inhibition of Nodal signaling. J. Mol. Cell. Biol. 7, 455–465 (2015).

    Article  Google Scholar 

  30. 30.

    Wilson, P. A., Lagna, G., Suzuki, A. & Hemmati-Brivanlou, A. Concentration-dependent patterning of the Xenopus ectoderm by BMP4 and its signal transducer Smad1. Development 124, 3177–3184 (1997).

    CAS  Article  Google Scholar 

  31. 31.

    Boel, N. M.-E., Hunter, M. C. & Edkins, A. L. LRP1 is required for novobiocin-mediated fibronectin turnover. Sci. Rep. 8, 11438 (2018).

    Article  ADS  Google Scholar 

  32. 32.

    Robert, E. & Guibaud, P. Maternal valproic acid and congenital neural tube defects. Lancet 320, 937 (1982).

    Article  Google Scholar 

  33. 33.

    Hughes, A., Greene, N. D. E., Copp, A. J. & Galea, G. L. Valproic acid disrupts the biomechanics of late spinal neural tube closure in mouse embryos. Mech. Dev. 149, 20–26 (2018).

    CAS  Article  Google Scholar 

  34. 34.

    Hamburger, V. & Hamilton, H. L. A series of normal stages in the development of the chick embryo. Dev. Dyn. 195, 231–272 (1992).

    CAS  Article  Google Scholar 

  35. 35.

    Theiler, K. The House Mouse: Atlas of Embryonic Development (Springer Science & Business Media, 1989).

  36. 36.

    Butler, A., Hoffman, P., Smibert, P., Papalexi, E. & Satija, R. Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat. Biotechnol. 36, 411–420 (2018).

    CAS  Article  Google Scholar 

  37. 37.

    Stuart, T. et al. Comprehensive integration of single-cell data. Cell 177, 1888–1902.e21 (2019).

    CAS  Article  Google Scholar 

  38. 38.

    McGinnis, C. S., Murrow, L. M. & Gartner, Z. J. DoubletFinder: doublet detection in single-cell RNA sequencing data using artificial nearest neighbors. Cell Syst. 8, 329–337.e4 (2019).

    CAS  Article  Google Scholar 

  39. 39.

    Hafemeister, C. & Satija, R. Normalization and variance stabilization of single-cell RNA-seq data using regularized negative binomial regression. Genome Biol. 20, 296 (2019).

    CAS  Article  Google Scholar 

  40. 40.

    Zhou, Y. et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat. Commun. 10, 1523 (2019).

    Article  ADS  Google Scholar 

  41. 41.

    Tirosh, I. et al. Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq. Science 352, 189–196 (2016)

    CAS  Article  ADS  Google Scholar 

Download references


S.J.S. acknowledges NIH grant R21 HD099598-0. E.K. acknowledges the Human Frontier Science Program for the long-term postdoctoral fellowship (LT000629/2018-L). A.H.K. acknowledges support of the National Science Foundation Graduate Research Fellowship Program (grant no. 1650114). E.D.S. acknowledges NSF grant 2013131. A.W. acknowledges Welch Foundation grant no C-2021. S.M.K.G. was recipient of an Early Postdoc Mobility grant and a Postdoc Mobility grant from the Swiss National Science Foundation (P2ZHP3_174753 and P400PB_186800). B.I.S. acknowledges NSF Physics grant 1707973. K.S.K. acknowledges support from Dr Miriam and Sheldon G. Adelson Medical Research Foundation. We thank Y. Huang and the Gladstone institute for sharing CTR2 hiPS cell line, and Brivanlou laboratory and the Rockefeller University for sharing RUES2 embryonic and reporter lines. We thank UCLA Technology Center for Genomics and Bioinformatics for sequencing of 10X libraries. We thank the Streichan and Kosik laboratory members, as well as the UCSB Stem Cell Center and associate director C. Arnold for fruitful discussions and advice. Finally, we thank A. M. Tayar for valuable comments on the manuscript.

Author information




E.K. and S.J.S. designed research. E.K., A.H.K., S.M.K.G., Y.W. and G.B. performed experiments. E.K., A.H.K. and H.C.M., analysed the data. A.W. and E.D.S. provided valuable information on data interpretation. E.K., H.C.M. and B.I.S. developed and analysed the physical model. S.M.K.G. and K.S.K. performed single-cell sequencing experiment and analysis. E.K., A.H.K. and S.J.S. wrote the manuscript. S.J.S. supervised the project.

Corresponding authors

Correspondence to Eyal Karzbrun or Sebastian J. Streichan.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review information Nature thanks Gioele La Manno, Matthias Lutolf and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data figures and tables

Extended Data Fig. 1 hPS cell-derived neural tube protocol.

(a) Scheme showing stem-cell derived neural tube protocol. (b) Scheme showing shape-controlled ECM pattern deposited on glass surface. (c) Seeding of hPS cells onto micropatterns results in two-dimensional cultures which are restricted to the micropattern geometry. ZO1-GFP indicates tight epithelia. (d) Adding 4% Matrigel to the media results in three-dimensional stem-cell cultures containing a single large lumen. (e) The hPS cell lumen forms a tight biochemical niche. hPS cell reporter line endogenously expressed tight junction protein ZO1 tagged with GFP (AICS-0023) is used. The 3D culture is exposed to 10kda dextran tagged with Texas Red fluorophore. Horizontal sections show that tight junctions are localized to the inner surface facing the lumen (green). Dextran is visible outside the tissue, but the lumen is devoid of dextran (red). (f) The fluorescence of dextran is quantified inside and outside the lumen, showing a 10-60 fold decrease in the lumen dextran concentration. Data are presented as mean values +/- SD. (N = 9). (g) Plot of the relative fraction of closed lumens as a function of Matrigel concentration, after 24hrs incubation. Data are presented as mean values +/- SD. Total N = 88 samples. (h) Scheme showing micropatterned hPS cell before and after exposure to Matrigel. (i) Vertical and (j) horizontal sections of pluripotent markers NANOG, OCT4 and SOX2, as well as epithelial marker ZO1. Images taken before and after lumen formation. (k) Horizontal section of an array of 2D (left) and 3D (right) stem-cell cultures immunostained with OCT4. (l) Lumen formation success rate as a function of micropattern size. Total N  = 510 samples. (m) Segmentation of single nuclei is used to count the number of cells in each sample. Each nucleus is labeled with a different color. (n) Total cell number in each sample (green curve) scales linearly with pattern area while cell density (black curve) remains constant. Data are presented as mean values +/- SD. Total N = 194 samples. Scale bars are 50μm (d,e,i,j), and 500μm (k).

Source data

Extended Data Fig. 2 In vitro morphogenesis Road Map.

(a) Schematic representation of in vitro morphogenesis in four differentiation protocols. Observed cell markers are indicated. (i) Exposure to BMP, without neural induction, results in an amnion-like tissue containing cells from three germ layers without formation of a neural fold. (ii) Neural tube morphogenesis is observed when neural induction is followed by exposure to BMP4. (iii) Homogenous expression of forebrain markers is observed under exposure to neural induction without BMP. (iv) Dorsoventral axis formation is observed in response to SHH activation (SAG) and WNT downregulation (IWP2). (b) Experimental timeline for the four protocols. Neural induction media includes N2 supplement and TGF-β inhibitor SB-431542 (SB). (c-f) Vertical sections of immunostained samples from the four differentiation protocols. (g) Horizontal section close to the glass surface and (h) an array of 20 samples. FOXA2+ cells always appear at the midline. Scale bar is 50μm (c,f-h), 50μm (d,e).

Extended Data Fig. 3 Folding morphogenesis is reproducible.

(a) Horizontal sections of 16 circular cultures from a single array exhibit stereotypic fate-patterning and morphology. (b) Bar plot showing success rate of neural tissue folding morphogenesis, in 13 experiments, with total N = 100 samples. Data are presented as mean values +/- SD. (c) Neural opening area as a function of time. Data are presented as mean values +/- SD. N = 3. (d) Neural closure period in the stem-cell system compared to chick, mouse, and human. (e) Images of six samples at 96hrs showing closed neural tissue (red, N-Cadherin) covered by non-neural ectoderm (cyan, E-cadherin). Samples were generated on rectangle  shaped micropatterns. Scale bar is 500μm (a), 50μm in (c,e).

Source data

Extended Data Fig. 4 Folding morphogenesis is not observed in 2D micropatterned cultures.

(a) Induction of neural pattern formation in 3D stem-cell cultures with a lumen triggers folding morphogenesis. In contrast triggering pattern formation in 2D micropatterned cultures does not trigger folding morphogenesis. Vertical sections of a 3D immunostained sample (top) and 2D immunostained sample (bottom). Imaged 72hrs after BMP. (b) Vertical section at 36hrs post BMP, and before folding. The neural tissue thickness and width are indicated. (c) Quantification of neural tissue thickness and width as a function of micropattern size analyzed 36hrs after BMP exposure. Data are presented as mean values +/- SD, Total N = 23. (d) Vertical section of immunostained samples during neural fold morphogenesis. Arrows indicate hinge points. (e) Distance between hinge points as a function of time after BMP exposure. Data are presented as mean values +/- SD, N = 3. Scale bars are 50μm.

Source data

Extended Data Fig. 5 Neural hinge points and apical surface characterization.

(a) Vertical sections of a hinge points showing changes in nuclei orientation across the hinge. (b) Neural crest membrane marker P75 appears in a layer of cells in between the neural tissue and surface ectoderm (ECAD). (c) Immunostaining of shroom3 and (d) F-actin at the neural apical surface. (e) Image projection of apical marker PKCζ and primary cilia marker ARL13B. (f) Sparse labeling of cells with CAAX-GFP. Labeled cells arise from a single clone indicating motility of cells following proliferation. (g) Live-cell imaging showing active filopodia (arrows) of neural cells as they change their position and shape over time. The continuous spectrum of colors indicates distance from glass (red 0μm -green - blue 35μm). Scale bars are 25μm (a,f) 10μm (b-e,g).

Extended Data Fig. 6 Live imaging of neural crest cells.

Time-lapse images showing motion of neural crest cells (SOX10-GFP). (a) Scheme indicating region in which live imaging was performed. (b) Low magnification image showing neural crest cells are organized in a stream. (c) Neural crest cell motion along the stream. A single cell is tracked and artificially colored in red. Scale bar 25μm.

Extended Data Fig. 7 Extracellular localization and expression during neural folding.

(a) Scheme showing ECM composition in three interfaces: surface-ectoderm/glass, surface-ectoderm/ neural-ectoderm and neural-ectoderm/glass. (b) Collagen and fibronectin fluorescence intensity compared at the three interfaces. Data are presented as mean values +/- SD. N = 3 independent samples (c) Vertical sections show fibronectin. (d) Horizontal sections near the glass interface show collagen V and fibronectin. Dashed line indicates border between neural and non-neural ectoderm. (e) UMAP projection and (f) violin plots of ECM genes and ECM modulators are localized to the surface ectoderm clusters. Average expression was 1-2 orders of magnitude higher in surface ectoderm vs neural ectoderm. (g) Vertical sections showing immunostaining of fibronectin matrix in control sample, ROCK inhibition (Y-27632, 10μM), and Novobiocin (100 μM). (h) Fluorescence intensities of immunostained fibronectin were measured at the neural/surface-ectoderm interface and at the neural/glass interface in all three experimental conditions. Immunostaining and imaging conditions were identical for all conditions. Fluorescence values were measured by averaging the raw fluorescence intensities in a region of interest in the image. The number of independent samples N is indicated within the plot. (i) Cell velocity of non-neural ectoderm cells measured by live-imaging of H2B-RFP cell line. No significant difference is observed in ROCK inhibited samples. N = 6. Data are presented as mean values +/- SD (b,h,i). Scale bars are 25 microns. Two-sided t-test was applied (b,h).

Source data

Extended Data Fig. 8 Single cell expression of key cell type markers.

(a) UMAP plots of single cell gene expression color coded by cluster identity (top) and experimental time point (bottom). Four main cell clusters are observed. (b) Temporal evolution of gene expression in the neural crest population indicate a transition from neural plate border/neural crest specification to migration and divarication. (c) Inferred cell-cycle populations for three clusters: neural ectoderm (NE), surface ectoderm (SE), and neural crest (NC). SE exhibits increase in G1 over time, and NC exhibits a decrease in G1. (d) Expression of key cell type markers from the four main clusters. (e) Sub-clustering of NE and SE. Neural cells sub-cluster mainly by experimental time and cell-cycle phase, which correspond to UMAP 1 and 2. In contrast, surface ectoderm cells sub-clusters do not cluster by cell-cycle or experimental time, suggesting additional differences in cell identity. (f) Scheme showing cellular process in each cluster. Titles are enriched GO pathways, and example genes appearing within each pathway.

Extended Data Fig. 9 Folding morphogenesis occurs in the absence of mesendoderm tissue.

(a) Experimental timeline to examine effect of neural induction on cell fate and folding morphogenesis. (b-d) Vertical (left) and horizontal (middle, right) sections showing that neural fates (NCAD) are upregulated with longer neural induction, whereas mesendoderm fates (Brachyury) are downregulated. (e,f) A small number of Brachyury+ cells (<10) is present in the protocol used to for neural tube morphogenesis. Total number of cells in the tissue is ~5000 cells. Scale bar is 50μm. Dome present in all cases.

Extended Data Fig. 10 Computational model of neural folding recapitulates in vitro morphology.

(a) Vertex model scheme. Vertices are connected by Hookean springs of length \({\triangle r}_{i}\). \({\vec{r}}_{i}^{a}\) and \({\vec{r}}_{i}^{b}\) are the apical and basal vertices (x and z) of cell \(i\). Neural cells have additional myosin dependent positive apical line tension, and there is a negative line tension between neural and non-neural basal surfaces (see Supplementary Notes details). (b) Time sequence of simulation showing initial cell pattern (neural tissue in red, non-neural in blue), followed by neural bending, folding and hinge formation, and closure. Each time unit corresponds to 1000 steps in the simulation. Apical lines in the neural domain strongly contract due to myosin, but never reach zero values. (c) Energy of the tissue reduces over time and reaches a minimum. (d) Interface radius and (e) Number of cells on the interface as a function of time. Both reach a steady-state value. (f) Phase diagram of the dimensionless parameters interface energy (\(\epsilon \)), neural contractility (\(\alpha \)). Each point represents a simulation that reached a steady state. Color code represents the final configuration of the simulations: neural fold closed (red), formation of neural/non-neural interface without closure (blue), and no interface formation (black). Shaded colors were drawn by hand to highlight regime boundaries. (g) Final configuration images (steady state shapes) from the simulation in different regions i) high neural apical contractility and high interface energy leads to a closed fold; ii) Low neural apical contractility and high interface energy leads to the formation of the interface but an open fold; iii) High neural apical contractility and low interface energy leads to a flat tissue. (h) The interface radius of curvature as a function of neural apical contractility. (i) The number of cells in an interface as a function of interface energy. (j) Bar plots showing the neural apical curvature, neural apical area, and interface size for the three conditions whose final configuration are shown in g. (k) Three simulations with a 40 non-neural cells, and a varying number of neural cells N = 10,30,40. Left: images of the simulation at time zero. Right: Images of the simulation at the hinge formation stage. Images were taken at this stage for consistency with experimental data (Fig. 4d, 72hrs). We observe the formation of a single hinge point for N = 10, and two lateral hinges for N = 30,40.

Source data

Supplementary information

Supplementary Information

This file contains a Supplementary Note on the physical model of neural folding.

Reporting Summary

Supplementary Table 1

This table contains top differentially expressed genes for cell clusters identified by scRNA-seq.

Supplementary Table 2

This table contains GO Analysis for cell clusters identified by scRNA-seq.

Supplementary Table 3

This table contains a list of media and reagents used in this work.

Supplementary Table 4

This table contains a list of primary and secondary antibodies used in this work.

Supplementary Video 1

Interkinetic nuclear migration. Live-cell imaging of interkinetic nuclear migrations. Time-lapse microscopy video using a 40× objective in the region of the neural fold from 72–94 h post BMP using a live reporter cell line (RUES2, H2B-RFP).

Supplementary Video 2

Neural progenitors. Live-cell imaging of CAAX-RFP reporter line mixed into a background non-fluorescent cell line. Color indicates distance from glass (red 0 μm–blue 35 μm). 72 h post-BMP. Elongated neural progenitor cells are observed exhibiting motion and filopodia. Frame rate: one frame every 30 min.

Supplementary Video 3

Primary cilia. 63× image slices of a fixed stem cell-derived neural tube at 120 h post-BMP. Apical marker PKCζ (red), and primary cilia marker ARL13B (yellow) are shown at the apical surface of neural tissue region.

Supplementary Video 4

Surface ectoderm cell migration. Live-cell imaging of H2B-RFP reporter line showing directional migration of non-neural cells in parallel to the culture long axis. Frame rate: one frame every 25 min.

Supplementary Video 5

Neural crest cells 3D. Three-dimensional view of a fixed stem cell-derived neural tube at 96 h post-BMP. Neural crest (SOX10) and non-neural ectoderm (ECAD) cells shown.

Supplementary Video 6

Neural crest migration. Live-cell imaging of SOX10-GFP reporter cell line showing a neural crest cell migrating and proliferating. Frame rate: one frame every 10 min.

Supplementary Video 7

Live-cell imaging of H2B-RFP reporter line demonstrating cell motility in control sample 72 h post-BMP. Frame rate: one frame every 5 min.

Supplementary Video 8

Live-cell imaging of H2B-RFP reporter line demonstrating cell motility in ROCK-inhibited sample. 72 h post-BMP. Frame rate: one frame every 5 min.

Source data

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Karzbrun, E., Khankhel, A.H., Megale, H.C. et al. Human neural tube morphogenesis in vitro by geometric constraints. Nature 599, 268–272 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI:


By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.


Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing