Almost all animal cells maintain a thin layer of actin filaments and associated proteins underneath the cell membrane. The actomyosin cortex is subject to internal stress patterns which result from the spatiotemporally regulated activity of non-muscle myosin II motors in the actin network. We study how these active stresses drive changes in cell shape and flows within the cortical layer, and how these cytoskeletal deformations and flows govern processes such as cell migration, cell division and organelle transport. Following a continuum mechanics approach, we develop theoretical descriptions for three different cellular processes, to obtain - in collaboration with experimental groups - a detailed and quantitative understanding of the underlying cytoskeletal mechanics. We investigate the forces and cortex flows involved in adhesion-independent cell migration in confinement. Many types of cell migration rely on the extension of protrusions at the leading edge, where the cells attach to the substrate with specific focal adhesions, and pull themselves forward, exerting stresses in the kPa range. In confined environments however, cells exhibit migration modes which are independent of specific adhesions. Combining hydrodynamic theory, microfluidics and quantitative imaging of motile, non-adherent carcinosarcoma cells, we analyze the mechanical behavior of cells during adhesion-independent migration. We find that the accumulation of active myosin motors in the rear part of these cells results in a retrograde cortical flow as well as the contraction of the cell body in the rear and expansion in the front, and we describe how both processes contribute to the translocation of the cells, depending on the geometric and mechanical parameters of the system. Importantly, we find that the involved propulsive forces are several orders of magnitude lower than during adhesive motility while the achieved migration velocities are similar. Moreover, the distribution of forces on the substrate during non-adhesive migration is fundamentally different, giving rise to a positive force dipole. In contrast to adhesive migration modes, non-adhesive cells move by exerting pushing forces at the rear, acting to expand rather than contract their substrate as they move. These differences may strongly affect hydrodynamic and/or deformational interactions between collectively migrating cells. In addition to the work outlined above, we study contractile ring formation in the actin cytoskeleton before and during cell division. While in disordered actin networks, myosin motor activity gives rise to isotropic stresses, the alignment of actin filaments in the cortex during cell division introduces a preferred direction for motor-filament interactions, resulting in anisotropies in the cortical stress. Actin filaments align in myosin-dependent shear flows, resulting in possible feedback between motor activity, cortical flows and actin organization. We investigate how the mechanical interplay of these different cortical properties gives rise to the formation of a cleavage furrow during cell division, describing the level of actin filament alignment at different points on the cortex with a nematic order parameter, in analogy to liquid crystal physics. We show that cortical anisotropies arising from shear-flow induced alignment patterns are sufficient to drive the ingression of cellular furrows, even in the absence of localized biochemical myosin up-regulation. This mechanism explains the characteristic appearance of pseudocleavage furrows in polarizing cells. Finally, we study the characteristic nuclear movements in pseudostratified epithelia during development. These tissues consist of highly proliferative, tightly packed and elongated cells, with nuclei actively travelling to the apical side of the epithelium before each cell division. We explore how cytoskeletal properties act together with the mechanics of the surrounding tissue to control the shape of single cells embedded in the epithelium, and investigate potential mechanisms underlying the observed nuclear movements. These findings form a theoretical basis for a more detailed characterization of processes in pseudostratified epithelia. Taken together, we present a continuum mechanics description of the actomyosin cell cortex, and successfully apply it to several different cell biological processes. Combining our theory with experimental work from collaborating groups, we provide new insights into different aspects of cell mechanics.

Comment les cellules établissent leurs formes et organisations internes est un problème biologique fondamental. Au cours de cette thèse, j’ai étudié le développement de la forme cellulaire et de la polarité chez la cellule de levure fissipare. Ces études sont fondées sur l’exploration de la façon dont les petites spores symétriques de levures se développent et s’organisent pour briser la symétrie pour la définition de leur tout premier axe de polarité. Dans une première partie, j’ai étudié les couplages entre la mécanique de surface de la paroi cellulaire des spores et la stabilité de domaines de polarité de Cdc42 qui contrôlent les aspects spatio-temporelles de la brisure de symétrie de ces spores. Dans une seconde partie, j’ai étudié les mécanismes par lesquels ces domaines de polarité contrôlent leur taille et l'adapte à la géométrie de la cellule, un processus vraisemblablement pertinents pour comprendre comment des domaines fonctionnels corticaux s’adaptent à la taille des cellules. Globalement, ces nouvelles recherches focalisant sur la façon dont les cellules développent dynamiquement leur forme et polarité de novo, permettent de mettre en évidence des couplages complexes dans la morphogenèse qui ne peuvent pas être testés en regardant les cellules à « l’état stationnaire» ou avec des outils génétiques. / How cells establish their proper shapes and organization is a fundamental biological problem. In this thesis, I investigated the dynamic development of cellular form and polarity in the rod-shape fission yeast cell. These studies are based on monitoring how small symmetric fission yeast spores grow and self-organize to break symmetry for the definition of their very first polarity axis. In a first part, I studied interplays between surface mechanics of the spore cell wall and the stability of Cdc42-based polarity domains which control spatio-temporal aspects of spore symmetry breaking. In a second part, I studied mechanisms by which these polarity domains control their width and adapt it to cell surface geometry, a process likely relevant to understand how functional cortical domains scale to cell size. Overall these novel investigations focusing on how cells dynamically develop their form and polarity de novo highlight complex feedbacks in morphogenesis that cannot be evidenced by looking at cells at “steady state” or with genetics.

Comment les cellules établissent leurs formes et organisations internes est un problème biologique fondamental. Au cours de cette thèse, j’ai étudié le développement de la forme cellulaire et de la polarité chez la cellule de levure fissipare. Ces études sont fondées sur l’exploration de la façon dont les petites spores symétriques de levures se développent et s’organisent pour briser la symétrie pour la définition de leur tout premier axe de polarité. Dans une première partie, j’ai étudié les couplages entre la mécanique de surface de la paroi cellulaire des spores et la stabilité de domaines de polarité de Cdc42 qui contrôlent les aspects spatio-temporelles de la brisure de symétrie de ces spores. Dans une seconde partie, j’ai étudié les mécanismes par lesquels ces domaines de polarité contrôlent leur taille et l'adapte à la géométrie de la cellule, un processus vraisemblablement pertinents pour comprendre comment des domaines fonctionnels corticaux s’adaptent à la taille des cellules. Globalement, ces nouvelles recherches focalisant sur la façon dont les cellules développent dynamiquement leur forme et polarité de novo, permettent de mettre en évidence des couplages complexes dans la morphogenèse qui ne peuvent pas être testés en regardant les cellules à « l’état stationnaire» ou avec des outils génétiques. / How cells establish their proper shapes and organization is a fundamental biological problem. In this thesis, I investigated the dynamic development of cellular form and polarity in the rod-shape fission yeast cell. These studies are based on monitoring how small symmetric fission yeast spores grow and self-organize to break symmetry for the definition of their very first polarity axis. In a first part, I studied interplays between surface mechanics of the spore cell wall and the stability of Cdc42-based polarity domains which control spatio-temporal aspects of spore symmetry breaking. In a second part, I studied mechanisms by which these polarity domains control their width and adapt it to cell surface geometry, a process likely relevant to understand how functional cortical domains scale to cell size. Overall these novel investigations focusing on how cells dynamically develop their form and polarity de novo highlight complex feedbacks in morphogenesis that cannot be evidenced by looking at cells at “steady state” or with genetics.

Das Kriechen von Zellen über Oberflächen spielt eine entscheidende Rolle bei lebenswichtigen Prozessen wie der Embryonalentwicklung, der Immunantwort und der Wundheilung, aber auch bei der Metastasenbildung. Die Zellbewegung erfolgt über die Bildung einer flachen Ausstülpung der Zellmembran, des Lamellipodiums. In dieser Arbeit wird ein mathematisches Modell entwickelt, das die Bildung, Stabilität und Stärke des Lamellipodiums, sowie die Dynamik der Zellvorderkante beschreibt. Dabei werden zwei Bereiche innerhalb des Lamellipodiums unterschieden. Im Hauptteil besteht es aus einem dichten Netzwerk von Aktinfilamenten, dem sogenannten Aktingel. An der Vorderkante wachsen die Enden der Aktinfilamente durch Polymerisation und bilden einen dynamischen Grenzbereich, die semiflexible Region. Das Gleichgewicht zwischen den Filamentkräften in der semiflexiblen Region und den viskosen sowie den äußeren Kräften bestimmt die Geschwindigkeit der Zellvorderkante. Eine Stabilitätsanalyse liefert Bedingungen für die Existenz stabiler Lamellipodien. Im Parameterbereich mit Filamentdichte Null können keine stabilen Lamellipodien existieren, aber aufgrund von Anregbarkeit trotzdem vorrübergehend gebildet werden. Hier beschreibt das Modell sehr gut das in Epithelzellen gemessene aufeinanderfolgende Vorschieben und Zurückziehen von Lamellipodien. Es zeigt, dass für die Zyklen prinzipiell keine Änderung in der Konzentration von Signalmolekülen innerhalb der Zelle notwendig ist. Das Modell wird auch auf die gemessene Kraft-Geschwindigkeits-Beziehung von Fischkeratozyten angewandt. Aufgrund der guten Übereinstimmung zwischen Experiment und Simulationen wird ein Mechanismus vorgeschlagen, der die charakteristischen Merkmale der Beziehung erklärt. Es wird gezeigt, dass die gemessene Kraft-Geschwindigkeits-Beziehung ein dynamisches Phänomen ist. Eine stationäre Beziehung, die unter der Bedingung gilt, dass die Zellen einer konstanten Kraft ausgesetzt sind, wird vorhergesagt. / Many cells move over surfaces during embryonic development, immune response, wound healing or cancer metastasis by protruding flat lamellipodia into the direction of migration. In this thesis, a mathematical model is developed that describes the formation of lamellipodia, their stability, strength and leading edge dynamics. Two regions inside the lamellipodium are distinguished in the model. The bulk contains a dense cross-linked actin filament network called actin gel. The newly polymerized tips of the actin filaments form a highly dynamic boundary layer at the leading edge called semiflexible region. The balance of filament forces on the membrane in the semiflexible region with viscous and external forces determines the velocity of the leading edge. A stability analysis defines criteria for the existence of stable lamellipodia. No stable lamellipodium can exist in the parameter regime with a filament density of zero. However, due to excitability, lamellipodia can still form transiently. The measured subsequent protrusions and retractions of lamellipodia in epithelial cells are very well reproduced by the excitable behavior. The modeling results show that in principle no signaling is necessary for cycles of protrusion and retraction. Furthermore, they are fitted to the force-velocity relation of keratocytes, which has been measured by placing a flexible cantilever into the cell''s path of migration. Due to the good agreement between experiment and simulations, a mechanism leading to the characteristic features of the force-velocity relation is suggested. Moreover, properties of the structure of the stable keratocyte lamellipodium, like the length of actin filaments and the branch point density, can be concluded. It is shown that the force-velocity relation measured with the cantilever is a dynamic phenomenon. A stationary force-velocity relation is predicted that should apply if cells experience a constant force, e.g. exerted by surrounding tissue.

Le fuseau mitotique assure la ségrégation des chromatides sœurs et le maintien de la poïdie des cellules filles. Le fuseau est composé de microtubules dynamiques (qui polymérisent et dépolymérisent continuellement), de nombreux moteurs moléculaires, d'agents de réticulations et de régulateurs. Bien que la structure du fuseau au niveau moléculaire soit connue, son fonctionnement reste délicat à comprendre, et nécessite la prise en compte de la dynamique de ses composants et leurs interactions. Les approches utilisées pour répondre à ces problématiques sont jusqu'à maintenant plutôt des approches in silico et in vitro. Il manque aujourd'hui une caractérisation de la mécanique du fuseau dans son contexte physiologique. Nous proposons une méthode non invasive basée sur de l'analyse d'image, combiné à une modélisation heuristique pour mesurer les paramètres mécaniques durant toute la division. Nous suivons les pôles du fuseau marqués par protéine fluorescente avec un taux acquisition rapide et une bonne résolution spatiale ce qui nous permet d'accéder aux fluctuations de longueur du fuseau in vivo. Avec la transformée de Fourier aux temps courts, nous calculons leurs densités spectrales de puissances — leurs signatures mécaniques. Ces spectres sont alors ajustés avec un modèle Kelvin — Voigt avec inertie (un ressort, un amortisseur et un terme inertiel en parallèle). Nous avons validé la méthode par des expériences numériques où nous retrouvons les évolutions des paramètres sur des données simulées et la calibration a été réalisée par l'utilisation de la rupture du fuseau induite par micro chirurgie laser ou par la génétique. Nous avons caractérisé le fuseau de l'embryon unicellulaire du nématode C. elegans. La méthaphase apparaît dominée par l'amortisseur, ce qui est cohérent avec la lente élongation du fuseau que nous observons. Mais contraste l'idée répandue de l'existence d'un mécanisme de maintien de la longueur du fuseau durant la métaphase. Au passage en anaphase, les trois paramètres mécaniques chutent, avant de réaugmenter environ 50 secondes après la transition pour réatindre un régime dominé de nouveau par l'amortisseur, ce qui suggère que les microtubules interpolaires jouent un rôle mineur durant l'élongation du fuseau en début d'anaphase. Dans la perspective de comprendre le lien entre la mécanique du fuseau et les interactions des acteurs moléculaires, nous avons partiellement supprimé un gène par sous-structure du fuseau. Nous avons alors retrouvé des comportements connus avec une perspective augmentée offerte par notre méthode. Cette méthode, ne va pas seulement permettre la compréhension fondamentale de la mécanique du fuseau, en remplaçant la modélisation du fuseau basé uniquement sur la longueur, mais aussi d'aller vers la prise en compte de la robustesse de fonctionnement du fuseau mitotique face aux défauts tel que la polyou l'aneuploïdie. / The mitotic spindle ensures the correct segregation of the sister chromatids to maintain ploidy in daughter cells. The spindle comprises dynamical microtubules (alternating polymerizing and depolymerizing), a variety of molecular motors, crosslinker and the regulators. Although the molecular grounds of spindle structure is well known, the link to its functions remain elusive, calling for including the dynamics of its components and their interactions. These questions were mostly investigated by in silico or in vitro approaches. But a detailed characterizing of spindle mechanics, in physiological conditions, is missing. We propose an image processing based, non invasive, method combined to an heuristic model to measure mechanical parameters of the mitotic spindle along time. We tracked fluorescently labeled spindle pole at high temporal and spatial resolution and measured the variations of spindle length, in vivo. We computed their power density spectrum using short time Fourier transform (sliding window) — a blueprint of spindle mechanics. Such a spectrum is then fitted with a Kelvin —Voigt model with inertia (a spring, a damper, an inertial element in parallel). We validated this method by recovering the mechanical parameters over time from simulated data and calibrated it uses laser and genetically induced spinlde cut. We characterized the mitotic spindle of the one-cell embryo of nematode C. elegans. Metaphase appeared dominated by damping element, consistent with the slow spindle elongation observed. But in contrast with the common thought that a mechanism maintains the spindle length during metaphase. At anaphase onset, all three parameters collapsed, before increasing about 50s later to reach a regime where damping dominated again, suggesting the overlapping spinlde microtubules may play a minor role in early anaphase spinlde elongation. In perspective of understanding how spindle mechanics emerge of molecular players interactions, we depleted one gene per splindle sub-structure — overlapped microtubules, kinetochore microtubules, central spindle and astral microtubules. We succefully recovered some known behavior but with the augmented insight offered by our method. This method paves the way not only towards understanding the fundamentals of spindle mechanics, superseding the degenerated modeling based on the sole spindle length but also towards acounting for spindle functional robustness towards defect as polyor aneuploidy.

La motilité des cellules est un phénomène fondamental en biologie souvent étudié sur des surfaces planes, conditions peu physiologiques. Nous avons analysé la migration cellulaire dans une matrice cellulaire 3D contenant de la fibronectine fluorescente. Nous démontrons que les cellules y sont confinées, et déforment leur environnement de manière cyclique avec une période de ~14 min avec deux centres de contractions à l’avant et à l’arrière de la cellule qui contractent avec un déphasage de ~3.5 min. Une perturbation de ces cycles entraîne une réduction de la motilité. Par l’utilisation d’inhibiteurs spécifiques, nous avons identifié l’acto-myosine comme étant l’acteur principal de ce phénomène. En imposant des contractions-relaxations locales par ablations laser, nous avons déclenché la motilité cellulaire ce qui confirme notre hypothèse. L’ensemble de cette étude met en évidence un nouveau mécanisme fondamental de dynamique cellulaire impliqué dans le mouvement des cellules. / Cell motility is an important process in Biology. It is mainly studied on 2D planar surfaces, whereas cells experience a confining 3D environment in vivo. We prepared a 3D Cell Derived Matrix (CDM) labeled with fluorescently labeled fibronectin, and strikingly cells managed to deform the matrix with specific patterns : contractions occur cyclically with two contraction centers at the front and at the back of the cell, with a period of ~14 min and a phase shift of ~3.5 min. These cycles enable cells to optimally migrate through the CDM, as perturbation of cycles led to reduced motility. Acto-myosin was established to be the driving actor of these cycles, by using specific inhibitors. We were able to trigger cell motility externally with local laser ablations, which supports this framework of two alternating contractions involved in motion. Altogether, this study reveals a new mechanism of dynamic cellular behaviour linked to cell motility.

L'objet de cette thèse est l'étude comparative de la réponse de cellules saines et malignes à la micro- et la nano-topographie de surface. L'interaction avec des stries de profondeur nanométrique est étudiée grâce à une méthode statistique. Nous démontrons que les cellules saines s'alignent plutôt sur des stries profondes, et que les cellules cancéreuses sont plus sensibles aux stries peu profondes. L'analyse des noyaux révèle qu’ils suivent l'alignement des corps cellulaires plus fidèlement dans le cas des cellules cancéreuses et que les noyaux de ces dernières sont plus sensibles aux stries de faible profondeur. Sur des micro-piliers nous démontrons que les cellules d’ostéosarcomes sont capables de se déformer et de faire adopter à leurs noyaux la forme de l'espace entre les piliers. Ceci ne se produit que durant la phase initiale d'adhésion pour les cellules saines. Les cellules immortalisées présentent un niveau intermédiaire de déformation. Quand l'espacement entre piliers est réduit, des différences de déformation sont révélées entre les lignées cancéreuses testées. La déformation est aussi liée au caractère cancéreux de kératinocytes et à l'expression de Cdx2 dans des lignées d'adénocarcinomes. Nous avons tenté d'expliquer ce mécanisme de déformation en l'attribuant au cytosquelette grâce à des analyses en microscopie confocale et avec des inhibiteurs du cytosquelette. L'imagerie de cellules vivantes a permis d'observer que les cellules sont très mobiles même quand elles sont déformées, que la mitose nécessite la perte de la déformation et que la déformation après mitose est plus rapide que la déformation pendant l'adhésion initiale des cellules. / This thesis deals with the differential response of healthy and cancerous cells to surface topography at the nanoscale and the microscale. Using a statistical method we developed we studied the interactions of cells with grooves of nanoscale depth. We demonstrate that healthy cells have a greater ability to align with deeper grooves, whereas cancerous cells are more sensitive to shallow grooves. Analysis reveals that the nucleus follows the alignment of the cell body more closely in cancerous cells, and that the nucleus of cancerous cells is more sensitive to shallow grooves.On microscale pillars we demonstrate for the first time that osteosarcoma cells deform to adopt the surface topography and that the deformation extends to the interior of the cell and in particular to the nucleus. We show that healthy cells only deform during the initial stages of adhesion and that immortalized cells show intermediate deformation between the healthy and cancerous cells. When the spacing between the pillars is reduced, differences in the deformation of different cancerous cell lines are detected. Deformation was also found to be related to the malignancy in keratinocytes, and related to the expression of Cdx2 in adenocarcinoma. The mechanism of deformation is tentatively attributed to the cytoskeleton and attempts to identify the main actors of deformation were performed using confocal microscopy and cytoskeleton inhibitors. Live cell imaging experiments reveal that the deformed cells are very mobile on the surfaces, loss of deformation is necessary for mitosis to occur and deformation after mitosis is more rapid than initial deformation upon adhesion to surfaces.

Much experimental evidence shows that the cytoskeleton is a downstream target and effector during cell death in numerous neurodegenerative diseases, including Parkinson's, Huntington's, and Alzheimer's diseases. However, recent evidence indicates that cytoskeletal dysfunction can also trigger neuronal death, by mechanisms as yet poorly understood. We studied a mathematical model of cytoskeleton-induced neuron death in which assembly control of the neuronal cytoskeleton interacts with both cellular stress levels and cytosolic free radical concentrations to trigger neurodegeneration. This trigger mechanism is further modulated by the presence of cell interactions in the form of a diffusible toxic factor released by dying neurons. We found that, consistent with empirical observations, the model produces one-hit exponential and sigmoid patterns of cell dropout. In all cases, cell dropout is exponential-tailed and described accurately by a gamma distribution. The transition between exponential and sigmoidal is gradual, and determined by a synergetic interaction between the magnitude of fluctuations in cytoskeleton assembly control and by the degree of cell coupling. We concluded that a single mechanism involving neuron interactions and fluctuations in cytoskeleton assembly control is compatible with the experimentally observed range of neuronal attrition kinetics. We also studied the transit of neurons through states intermediate between initial viability and cell death. We found that the stochastic flow of neuron fate, from viability to cell death, self-organizes into two distinct temporal phases. There is a rapid relaxation of the initial neuron population to a more disordered phase that is long-lived, or metastable, with respect to the time scales of change in single cells. Strikingly, cellular egress from this metastable phase follows the one-hit kinetic pattern of exponential decline now established as a principal hallmark of cell death in neurodegenerative disorders. Intermediate state metastability may therefore be an important element in the systems biology of one-hit neurodegeneration. Further, we studied the full spatiotemporal dynamics of death factor pulses released from dying neurons to emphasize the effects of the cell-to-cell coupling strength on neuron death rates. The rate of neuron cell loss monotonically increased with increased diffusion-dependent intercellular communication. Death factor diffusion effects may therefore be important moderators of one-hit neurodegeneration.

Much experimental evidence shows that the cytoskeleton is a downstream target and effector during cell death in numerous neurodegenerative diseases, including Parkinson's, Huntington's, and Alzheimer's diseases. However, recent evidence indicates that cytoskeletal dysfunction can also trigger neuronal death, by mechanisms as yet poorly understood. We studied a mathematical model of cytoskeleton-induced neuron death in which assembly control of the neuronal cytoskeleton interacts with both cellular stress levels and cytosolic free radical concentrations to trigger neurodegeneration. This trigger mechanism is further modulated by the presence of cell interactions in the form of a diffusible toxic factor released by dying neurons. We found that, consistent with empirical observations, the model produces one-hit exponential and sigmoid patterns of cell dropout. In all cases, cell dropout is exponential-tailed and described accurately by a gamma distribution. The transition between exponential and sigmoidal is gradual, and determined by a synergetic interaction between the magnitude of fluctuations in cytoskeleton assembly control and by the degree of cell coupling. We concluded that a single mechanism involving neuron interactions and fluctuations in cytoskeleton assembly control is compatible with the experimentally observed range of neuronal attrition kinetics. We also studied the transit of neurons through states intermediate between initial viability and cell death. We found that the stochastic flow of neuron fate, from viability to cell death, self-organizes into two distinct temporal phases. There is a rapid relaxation of the initial neuron population to a more disordered phase that is long-lived, or metastable, with respect to the time scales of change in single cells. Strikingly, cellular egress from this metastable phase follows the one-hit kinetic pattern of exponential decline now established as a principal hallmark of cell death in neurodegenerative disorders. Intermediate state metastability may therefore be an important element in the systems biology of one-hit neurodegeneration. Further, we studied the full spatiotemporal dynamics of death factor pulses released from dying neurons to emphasize the effects of the cell-to-cell coupling strength on neuron death rates. The rate of neuron cell loss monotonically increased with increased diffusion-dependent intercellular communication. Death factor diffusion effects may therefore be important moderators of one-hit neurodegeneration.

This thesis presents a whole-embryo finite element model of neurulation -- the first of its kind. An advanced, multiscale finite element approach is used to capture the mechanical interactions that occur across cellular, tissue and whole-embryo scales. Cell-based simulations are used to construct a system of constitutive equations for embryonic tissue fabric evolution under different scenarios including bulk deformation, cell annealing, mitosis, and Lamellipodia effect. Experimental data are used to determine the parameters in these equations. Techniques for obtaining images of live embryos, serial sections of fixed embryo fabric parameters, and material properties of embryonic tissues are used. Also a spatial-temporal correlation system is introduced to organize and correlate the data and to construct the finite element model. Biological experiments have been conducted to verify the validity of this constitutive model. A full functional finite element analysis package has been written and is used to conduct computational simulations. A simplified contact algorithm is introduced to address the element permeability issue. Computational simulations of different cases have been conducted to investigate possible causes of neural tube defects. Defect cases including neural plate defect, non-neural epidermis defect, apical constriction defect, and convergent extension defect are compared with the case of normal embryonic development. Corresponding biological experiments are included to support these defect cases. A case with biomechanical feedbacks on non-neural epidermis is also discussed in detail with biological experiments and computational simulations. Its comparison with the normal case indicates that the introduction of biomechanical feedbacks can yield more realistic simulation results.