Embryonic development depends on the precise coordination of cell fate specification, patterning and morphogenesis. Although great strides have been made in the molecular understanding of each of these processes, how their interplay governs the formation of complex tissues remains poorly understood. New techniques for experimental manipulation and image quantification enable the study of development in unprecedented detail, resulting in new hypotheses on the interactions between known components. By expressing these hypotheses in terms of rules and equations, computational modeling and simulation allows one to test their consistency against experimental data. However, new computational methods are required to represent and integrate the network of interactions between gene regulation, signaling and biomechanics that extend over the molecular, cellular and tissue scales. In this thesis, I present a framework that facilitates computational modeling of multiscale multicellular systems and apply it to investigate pancreatic development and the formation of vascular networks. This framework is based on the integration of discrete cell-based models with continuous models for intracellular regulation and intercellular signaling. Specifically, gene regulatory networks are represented by differential equations to analyze cell fate regulation; interactions and distributions of signaling molecules are modeled by reaction-diffusion systems to study pattern formation; and cell-cell interactions are represented in cell-based models to investigate morphogenetic processes. A cell-centered approach is adopted that facilitates the integration of processes across the scales and simultaneously constrains model complexity. The computational methods that are required for this modeling framework have been implemented in the software platform Morpheus. This modeling and simulation environment enables the development, execution and analysis of multi-scale models of multicellular systems. These models are represented in a new domain-specific markup language that separates the biological model from the computational methods and facilitates model storage and exchange. Together with a user-friendly graphical interface, Morpheus enables computational modeling of complex developmental processes without programming and thereby widens its accessibility for biologists. To demonstrate the applicability of the framework to problems in developmental biology, two case studies are presented that address different aspects of the interplay between cell fate specification, patterning and morphogenesis. In the first, I focus on the interplay between cell fate stability and intercellular signaling. Specifically, two studies are presented that investigate how mechanisms of cell-cell communication affect cell fate regulation and spatial patterning in the pancreatic epithelium. Using bifurcation analysis and simulations of spatially coupled differential equations, it is shown that intercellular communication results in a multistability of gene expression states that can explain the scattered spatial distribution and low cell type ratio of nascent islet cells. Moreover, model analysis shows that disruption of intercellular communication induces a transition between gene expression states that can explain observations of in vitro transdifferentiation from adult acinar cells into new islet cells. These results emphasize the role of the multicellular context in cell fate regulation during development and may be used to optimize protocols for cellular reprogramming. The second case study focuses on the feedback between patterning and morphogenesis in the context of the formation of vascular networks. Integrating a cell-based model of endothelial chemotaxis with a reaction-diffusion model representing signaling molecules and extracellular matrix, it is shown that vascular network patterns with realistic morphometry can arise when signaling factors are retained by cell-modified matrix molecules. Through the validation of this model using in vitro assays, quantitative estimates are obtained for kinetic parameters that, when used in quantitative model simulations, confirm the formation of vascular networks under measured biophysical conditions. These results demonstrate the key role of the extracellular matrix in providing spatial guidance cues, a fact that may be exploited to enhance vascularization of engineered tissues. Together, the modeling framework, software platform and case studies presented in this thesis demonstrate how cell-centered computational modeling of multi-scale and multicellular systems provide powerful tools to help disentangle the complex interplay between cell fate specification, patterning and morphogenesis during embryonic development.

The cell is the central functional unit of life. Cell behaviors, such as cell division, movements, differentiation, cell death as well as cell shape and size changes, determine how tissues change shape and grow during regeneration and development. However, a generally applicable framework to measure and describe the behavior of the multitude of cells in a developing tissue is still lacking. Furthermore, the specific contribution of individual cell behaviors, and how exactly these cell behaviors collectively lead to the morphogenesis and growth of tissues are not clear for many developmental and regenerative processes. A promising strategy to fill these gaps is the continuing effort of making developmental biology a quantitative science. Recent advances in methods, especially in imaging, enable measurements of cell behaviors and tissue shapes in unprecedented detail and accuracy. Consequently, formalizing hypotheses in terms of mathematical models to obtain testable quantitative predictions is emerging as a powerful tool. Tests of the hypotheses involve the comparison of model predictions to experimentally observed data. The available data is often noisy and based on only few samples. Hence, this comparison of data and model predictions often requires very careful use of statistical inference methods. If one chooses this quantitative approach, the challenges are the choice of observables, i.e. what to measure, and the design of appropriate data-driven models to answer relevant questions. In this thesis, I applied this data-driven modeling approach to vertebrate morphogenesis, growth and regeneration. In particular, I study spinal cord and muscle regeneration in axolotl, muscle development in zebrafish, and neuron development and maintenance in the adult human brain. To do so, I analyzed images to quantify cell behaviors and tissue shapes. Especially for cell behaviors in post-embryonic tissues, measurements of some cell behavior parameters, such as the proliferation rate, could not be made directly. Hence, I developed mathematical models that are specifically designed to infer these parameters from indirect experimental data. To understand how cell behaviors shape tissues, I developed mechanistic models that causally connect the cell and tissue scales. Specifically, I first investigated the behaviors of neural stem cells that underlie the regenerative outgrowth of the spinal cord after tail amputation in the axolotl. To do so, I quantified all relevant cell behaviors. A detailed analysis of the proliferation pattern in space and time revealed that the cell cycle is accelerated between 3-4 days after amputation in a high-proliferation zone, initially spanning from 800 µm anterior to the amputation plane. The activation of quiescent stem cells and cell movements into the high-proliferation zone also contribute to spinal cord growth but I did not find contributions by cellular rearrangements or cell shape changes. I developed a mathematical model of spinal cord outgrowth involving all contributing cell behaviors which revealed that the acceleration of the cell cycle is the major driver of spinal cord outgrowth. To compare the behavior of neural stem cells with cell behaviors in the regenerating muscle tissue that surrounds the spinal cord, I also quantified proliferation of mesenchymal progenitor cells and found similar proliferation parameters. I showed that the zone of mesenchymal progenitors that gives rise to the regenerating muscle segments is at least 350 µm long, which is consistent with the length of the high-proliferation zone in the spinal cord. Second, I investigated shape changes in developing zebrafish muscle segments by quantifying time-lapse movies of developing zebrafish embryos. These data challenged or ruled out a number of previously proposed mechanisms. Motivated by reported cellular behaviors happening simultaneously in the anterior segments, I had previously proposed the existence of a simple tension-and-resistance mechanism that shapes the muscle segments. Here, I could verify the predictions of this mechanism for the final segment shape pattern. My results support the notion that a simple physical mechanism suffices to self-organize the observed spatiotemporal pattern in the muscle segments. Third, I corroborated and refined previous estimates of neuronal cell turnover rates in the adult human hippocampus. Previous work approached this question by combining quantitative data and mathematical modeling of the incorporation of the carbon isotope C-14. I reanalyzed published data using the published deterministic neuron turnover model but I extended the model by a better justified measurement error model. Most importantly, I found that human adult neurogenesis might occur at an even higher rate than currently believed. The tools I used throughout were (1) the careful quantification of the involved processes, mainly by image analysis, and (2) the derivation and application of mathematical models designed to integrate the data through (3) statistical inference. Mathematical models were used for different purposes such as estimating unknown parameters from indirect experiments, summarizing datasets with a few meaningful parameters, formalizing mechanistic hypotheses, as well as for model-guided experimental planning. I venture an outlook on how additional open questions regarding cell turnover measurements could be answered using my approach. Finally, I conclude that the mechanistic understanding of development and regeneration can be advanced by comparing quantitative data to the predictions of specifically designed mathematical models by means of statistical inference methods.

The cell is the central functional unit of life. Cell behaviors, such as cell division, movements, differentiation, cell death as well as cell shape and size changes, determine how tissues change shape and grow during regeneration and development. However, a generally applicable framework to measure and describe the behavior of the multitude of cells in a developing tissue is still lacking. Furthermore, the specific contribution of individual cell behaviors, and how exactly these cell behaviors collectively lead to the morphogenesis and growth of tissues are not clear for many developmental and regenerative processes. A promising strategy to fill these gaps is the continuing effort of making developmental biology a quantitative science. Recent advances in methods, especially in imaging, enable measurements of cell behaviors and tissue shapes in unprecedented detail and accuracy. Consequently, formalizing hypotheses in terms of mathematical models to obtain testable quantitative predictions is emerging as a powerful tool. Tests of the hypotheses involve the comparison of model predictions to experimentally observed data. The available data is often noisy and based on only few samples. Hence, this comparison of data and model predictions often requires very careful use of statistical inference methods. If one chooses this quantitative approach, the challenges are the choice of observables, i.e. what to measure, and the design of appropriate data-driven models to answer relevant questions. In this thesis, I applied this data-driven modeling approach to vertebrate morphogenesis, growth and regeneration. In particular, I study spinal cord and muscle regeneration in axolotl, muscle development in zebrafish, and neuron development and maintenance in the adult human brain. To do so, I analyzed images to quantify cell behaviors and tissue shapes. Especially for cell behaviors in post-embryonic tissues, measurements of some cell behavior parameters, such as the proliferation rate, could not be made directly. Hence, I developed mathematical models that are specifically designed to infer these parameters from indirect experimental data. To understand how cell behaviors shape tissues, I developed mechanistic models that causally connect the cell and tissue scales. Specifically, I first investigated the behaviors of neural stem cells that underlie the regenerative outgrowth of the spinal cord after tail amputation in the axolotl. To do so, I quantified all relevant cell behaviors. A detailed analysis of the proliferation pattern in space and time revealed that the cell cycle is accelerated between 3-4 days after amputation in a high-proliferation zone, initially spanning from 800 µm anterior to the amputation plane. The activation of quiescent stem cells and cell movements into the high-proliferation zone also contribute to spinal cord growth but I did not find contributions by cellular rearrangements or cell shape changes. I developed a mathematical model of spinal cord outgrowth involving all contributing cell behaviors which revealed that the acceleration of the cell cycle is the major driver of spinal cord outgrowth. To compare the behavior of neural stem cells with cell behaviors in the regenerating muscle tissue that surrounds the spinal cord, I also quantified proliferation of mesenchymal progenitor cells and found similar proliferation parameters. I showed that the zone of mesenchymal progenitors that gives rise to the regenerating muscle segments is at least 350 µm long, which is consistent with the length of the high-proliferation zone in the spinal cord. Second, I investigated shape changes in developing zebrafish muscle segments by quantifying time-lapse movies of developing zebrafish embryos. These data challenged or ruled out a number of previously proposed mechanisms. Motivated by reported cellular behaviors happening simultaneously in the anterior segments, I had previously proposed the existence of a simple tension-and-resistance mechanism that shapes the muscle segments. Here, I could verify the predictions of this mechanism for the final segment shape pattern. My results support the notion that a simple physical mechanism suffices to self-organize the observed spatiotemporal pattern in the muscle segments. Third, I corroborated and refined previous estimates of neuronal cell turnover rates in the adult human hippocampus. Previous work approached this question by combining quantitative data and mathematical modeling of the incorporation of the carbon isotope C-14. I reanalyzed published data using the published deterministic neuron turnover model but I extended the model by a better justified measurement error model. Most importantly, I found that human adult neurogenesis might occur at an even higher rate than currently believed. The tools I used throughout were (1) the careful quantification of the involved processes, mainly by image analysis, and (2) the derivation and application of mathematical models designed to integrate the data through (3) statistical inference. Mathematical models were used for different purposes such as estimating unknown parameters from indirect experiments, summarizing datasets with a few meaningful parameters, formalizing mechanistic hypotheses, as well as for model-guided experimental planning. I venture an outlook on how additional open questions regarding cell turnover measurements could be answered using my approach. Finally, I conclude that the mechanistic understanding of development and regeneration can be advanced by comparing quantitative data to the predictions of specifically designed mathematical models by means of statistical inference methods.

Embryonic development depends on the precise coordination of cell fate specification, patterning and morphogenesis. Although great strides have been made in the molecular understanding of each of these processes, how their interplay governs the formation of complex tissues remains poorly understood. New techniques for experimental manipulation and image quantification enable the study of development in unprecedented detail, resulting in new hypotheses on the interactions between known components. By expressing these hypotheses in terms of rules and equations, computational modeling and simulation allows one to test their consistency against experimental data. However, new computational methods are required to represent and integrate the network of interactions between gene regulation, signaling and biomechanics that extend over the molecular, cellular and tissue scales. In this thesis, I present a framework that facilitates computational modeling of multiscale multicellular systems and apply it to investigate pancreatic development and the formation of vascular networks. This framework is based on the integration of discrete cell-based models with continuous models for intracellular regulation and intercellular signaling. Specifically, gene regulatory networks are represented by differential equations to analyze cell fate regulation; interactions and distributions of signaling molecules are modeled by reaction-diffusion systems to study pattern formation; and cell-cell interactions are represented in cell-based models to investigate morphogenetic processes. A cell-centered approach is adopted that facilitates the integration of processes across the scales and simultaneously constrains model complexity. The computational methods that are required for this modeling framework have been implemented in the software platform Morpheus. This modeling and simulation environment enables the development, execution and analysis of multi-scale models of multicellular systems. These models are represented in a new domain-specific markup language that separates the biological model from the computational methods and facilitates model storage and exchange. Together with a user-friendly graphical interface, Morpheus enables computational modeling of complex developmental processes without programming and thereby widens its accessibility for biologists. To demonstrate the applicability of the framework to problems in developmental biology, two case studies are presented that address different aspects of the interplay between cell fate specification, patterning and morphogenesis. In the first, I focus on the interplay between cell fate stability and intercellular signaling. Specifically, two studies are presented that investigate how mechanisms of cell-cell communication affect cell fate regulation and spatial patterning in the pancreatic epithelium. Using bifurcation analysis and simulations of spatially coupled differential equations, it is shown that intercellular communication results in a multistability of gene expression states that can explain the scattered spatial distribution and low cell type ratio of nascent islet cells. Moreover, model analysis shows that disruption of intercellular communication induces a transition between gene expression states that can explain observations of in vitro transdifferentiation from adult acinar cells into new islet cells. These results emphasize the role of the multicellular context in cell fate regulation during development and may be used to optimize protocols for cellular reprogramming. The second case study focuses on the feedback between patterning and morphogenesis in the context of the formation of vascular networks. Integrating a cell-based model of endothelial chemotaxis with a reaction-diffusion model representing signaling molecules and extracellular matrix, it is shown that vascular network patterns with realistic morphometry can arise when signaling factors are retained by cell-modified matrix molecules. Through the validation of this model using in vitro assays, quantitative estimates are obtained for kinetic parameters that, when used in quantitative model simulations, confirm the formation of vascular networks under measured biophysical conditions. These results demonstrate the key role of the extracellular matrix in providing spatial guidance cues, a fact that may be exploited to enhance vascularization of engineered tissues. Together, the modeling framework, software platform and case studies presented in this thesis demonstrate how cell-centered computational modeling of multi-scale and multicellular systems provide powerful tools to help disentangle the complex interplay between cell fate specification, patterning and morphogenesis during embryonic development.