Multi-scale computational study of the mechanical regulation of cell mitotic rounding in epithelia

Reverse engineering morphogenesis through Bayesian optimization of physics-based models

Morphogenetic programs coordinate cell signaling and mechanical interactions to shape organs. In systems and synthetic biology, a key challenge is determining optimal cellular interactions for predicting organ shape, size, and function. Physics-based models defining the subcellular force distribution facilitate this, but it is challenging to calibrate parameters in these models from data. To solve this inverse problem, we created a Bayesian optimization framework to determine the optimal cellular force distribution such that the predicted organ shapes match the experimentally observed organ shapes. This integrative framework employs Gaussian Process Regression, a non-parametric kernel-based probabilistic machine learning modeling paradigm, to learn the mapping functions relating to the morphogenetic programs that maintain the final organ shape. We calibrated and tested the method on Drosophila wing imaginal discs to study mechanisms that regulate epithelial processes ranging from development to cancer. The parameter estimation framework successfully infers the underlying changes in core parameters needed to match simulation data with imaging data of wing discs perturbed with collagenase. The computational pipeline identifies distinct parameter sets mimicking wild-type shapes. It enables a global sensitivity analysis to support the regulation of actomyosin contractility and basal ECM stiffness to generate and maintain the curved shape of the wing imaginal disc. The optimization framework, combined with experimental imaging, identified that Piezo, a mechanosensitive ion channel, impacts fold formation by regulating the apical-basal balance of actomyosin contractility and elasticity of ECM. This workflow is extensible toward reverse-engineering morphogenesis across organ systems and for real-time control of complex multicellular systems.

Balancing competing effects of tissue growth and cytoskeletal regulation during Drosophila wing disc development

How a developing organ robustly coordinates the cellular mechanics and growth to reach a final size and shape remains poorly understood. Through iterations between experiments and model simulations that include a mechanistic description of interkinetic nuclear migration, we show that the local curvature, height, and nuclear positioning of cells in the Drosophila wing imaginal disc are defined by the concurrent patterning of actomyosin contractility, cell-ECM adhesion, ECM stiffness, and interfacial membrane tension. We show that increasing cell proliferation via different growth-promoting pathways results in two distinct phenotypes. Triggering proliferation through insulin signaling increases basal curvature, but an increase in growth through Dpp signaling and Myc causes tissue flattening. These distinct phenotypic outcomes arise from differences in how each growth pathway regulates the cellular cytoskeleton, including contractility and cell-ECM adhesion. The coupled regulation of proliferation and cytoskeletal regulators is a general strategy to meet the multiple context-dependent criteria defining tissue morphogenesis.

The Effect of Substrate Stiffness on Elastic Force Transmission in the Epithelial Monolayers over Short Timescales

Purpose

The importance of mechanical forces and microenvironment in guiding cellular behavior has been widely accepted. Together with the extracellular matrix (ECM), epithelial cells form a highly connected mechanical system subjected to various mechanical cues from their environment, such as ECM stiffness, and tensile and compressive forces. ECM stiffness has been linked to many pathologies, including tumor formation. However, our understanding of the effect of ECM stiffness and its heterogeneities on rapid force transduction in multicellular systems has not been fully addressed.

Methods

We used experimental and computational methods. Epithelial cells were cultured on elastic hydrogels with fluorescent nanoparticles. Single cells were moved by a micromanipulator, and epithelium and substrate deformation were recorded. We developed a computational model to replicate our experiments and quantify the force distribution in the epithelium. Our model further enabled simulations with local stiffness gradients.

Results

We found that substrate stiffness affects the force transduction and the cellular deformation following an external force. Also, our results indicate that the heterogeneities, e.g., gradients, in the stiffness can substantially influence the strain redistribution in the cell monolayers. Furthermore, we found that the cells' apico-basal elasticity provides a level of mechanical isolation between the apical cell-cell junctions and the basal focal adhesions.

Conclusions

Our simulation results show that increased ECM stiffness, e.g., due to a tumor, can mechanically isolate cells and modulate rapid mechanical signaling between cells over distances. Furthermore, the developed model has the potential to facilitate future studies on the interactions between epithelial monolayers and elastic substrates.

Supplementary Information

The online version of this article (10.1007/s12195-023-00772-0) contains supplementary material, which is available to authorized users.

The mechanics of cephalic furrow formation in the Drosophila embryo

Cephalic furrow formation (CFF) is a major morphogenetic movement during gastrulation in Drosophila melanogaster embryos that gives rise to a deep, transitory epithelial invagination. Recent studies have identified the individual cell shape changes that drive the initiation and progression phases of CFF; however, the underlying mechanics are not yet well understood. During the progression phase, the furrow deepens as columnar cells from both the anterior and posterior directions fold inwards rotating by 90°. To analyze the mechanics of this process, we have developed an advanced two-dimensional lateral vertex model that includes multinode representation of cellular membranes and allows us to capture the membrane curvature associated with pressure variation. Our investigations reveal some key potential mechanical features of CFF, as follows. When cells begin to roll over the cephalic furrow cleft, they become wedge shaped as their apical cortices and overlying membranes expand, lateral cortices and overlying membranes release tension, internal pressures drop, and basal cortices and membranes contract. Then, cells reverse this process by shortening apical cortices and membranes, increasing lateral tension, and causing internal pressures to rise. Since the basal membranes expand, the cells recover their rotated columnar shape once in the furrow. Interestingly, our findings indicate that the basal membranes may be passively reactive throughout the progression phase. We also find that the smooth rolling of cells over the cephalic furrow cleft necessitates that internalized cells provide a solid base through high levels of membrane tension and internal pressure, which allows the transmission of tensile force that pulls new cells into the furrow. These results lead us to suggest that CFF helps to establish a baseline tension across the apical surface of the embryo to facilitate cellular coordination of other morphogenetic movements via mechanical stress feedback mechanisms.

Reverse engineering morphogenesis through Bayesian optimization of physics-based models

Morphogenetic programs direct the cell signaling and nonlinear mechanical interactions between multiple cell types and tissue layers to define organ shape and size. A key challenge for systems and synthetic biology is determining optimal combinations of intra- and inter-cellular interactions to predict an organ's shape, size, and function. Physics-based mechanistic models that define the subcellular force distribution facilitate this, but it is extremely challenging to calibrate parameters in these models from data. To solve this inverse problem, we created a Bayesian optimization framework to determine the optimal cellular force distribution such that the predicted organ shapes match the desired organ shapes observed within the experimental imaging data. This integrative framework employs Gaussian Process Regression (GPR), a non-parametric kernel-based probabilistic machine learning modeling paradigm, to learn the mapping functions relating to the morphogenetic programs that generate and maintain the final organ shape. We calibrated and tested the method on cross-sections of Drosophila wing imaginal discs, a highly informative model organ system, to study mechanisms that regulate epithelial processes that range from development to cancer. As a specific test case, the parameter estimation framework successfully infers the underlying changes in core parameters needed to match simulation data with time series imaging data of wing discs perturbed with collagenase. Unexpectedly, the framework also identifies multiple distinct parameter sets that generate shapes similar to wild-type organ shapes. This platform enables an efficient, global sensitivity analysis to support the necessity of both actomyosin contractility and basal ECM stiffness to generate and maintain the curved shape of the wing imaginal disc. The optimization framework, combined with fixed tissue imaging, identified that Piezo, a mechanosensitive ion channel, impacts fold formation by regulating the apical-basal balance of actomyosin contractility and elasticity of ECM. This framework is extensible toward reverse-engineering the morphogenesis of any organ system and can be utilized in real-time control of complex multicellular systems.

Constriction imposed by basement membrane regulates developmental cell migration

The basement membrane (BM) is a specialized extracellular matrix (ECM), which underlies or encases developing tissues. Mechanical properties of encasing BMs have been shown to profoundly influence the shaping of associated tissues. Here, we use the migration of the border cells (BCs) of the Drosophila egg chamber to unravel a new role of encasing BMs in cell migration. BCs move between a group of cells, the nurse cells (NCs), that are enclosed by a monolayer of follicle cells (FCs), which is, in turn, surrounded by a BM, the follicle BM. We show that increasing or reducing the stiffness of the follicle BM, by altering laminins or type IV collagen levels, conversely affects BC migration speed and alters migration mode and dynamics. Follicle BM stiffness also controls pairwise NC and FC cortical tension. We propose that constraints imposed by the follicle BM influence NC and FC cortical tension, which, in turn, regulate BC migration. Encasing BMs emerge as key players in the regulation of collective cell migration during morphogenesis.

A multiscale chemical-mechanical model predicts impact of morphogen spreading on tissue growth

The exact mechanism controlling cell growth remains a grand challenge in developmental biology and regenerative medicine. The Drosophila wing disc tissue serves as an ideal biological model to study mechanisms involved in growth regulation. Most existing computational models for studying tissue growth focus specifically on either chemical signals or mechanical forces. Here we developed a multiscale chemical-mechanical model to investigate the growth regulation mechanism based on the dynamics of a morphogen gradient. By comparing the spatial distribution of dividing cells and the overall tissue shape obtained in model simulations with experimental data of the wing disc, it is shown that the size of the domain of the Dpp morphogen is critical in determining tissue size and shape. A larger tissue size with a faster growth rate and more symmetric shape can be achieved if the Dpp gradient spreads in a larger domain. Together with Dpp absorbance at the peripheral zone, the feedback regulation that downregulates Dpp receptors on the cell membrane allows for further spreading of the morphogen away from its source region, resulting in prolonged tissue growth at a more spatially homogeneous growth rate.

The Mechanical Influence of Densification on Initial Epithelial Architecture

Epithelial tissues are the most abundant tissue type in animals, lining body cavities and generating compartment barriers. The function of a monolayer epithelium - whether protective, secretory, absorptive, or filtrative -relies on regular tissue architecture with respect to the apical-basal axis. Using an unbiased 3D analysis pipeline developed in our lab, we previously showed that epithelial tissue architectures in culture can be divided into distinct developmental categories, and that these are intimately connected to cell density: at sparse densities, cultured epithelial cell layers have a squamous morphology (Immature); at intermediate densities, these layers develop lateral cell-cell borders and rounded cell apices (Intermediate); cells at the highest densities reach their full height and demonstrate flattened apices (Mature). These observations prompted us to ask whether epithelial architecture emerges from the mechanical constraints of densification, and to what extent a hallmark feature of epithelial cells, namely cell-cell adhesion, contributes. In other words, to what extent is the shape of cells in an epithelial layer a simple matter of sticky, deformable objects squeezing together? We addressed this problem using a combination of computational modeling and experimental manipulations. Our results show that the first morphological transition, from Immature to Intermediate, can be explained simply by cell crowding. Additionally, we identify a new division (and thus transition) within the Intermediate category, and find that this second morphology relies on cell-cell adhesion.

Cell cycle dynamics control fluidity of the developing mouse neuroepithelium

As developing tissues grow in size and undergo morphogenetic changes, their material properties may be altered. Such changes result from tension dynamics at cell contacts or cellular jamming. Yet, in many cases, the cellular mechanisms controlling the physical state of growing tissues are unclear. We found that at early developmental stages, the epithelium in the developing mouse spinal cord maintains both high junctional tension and high fluidity. This is achieved via a mechanism in which interkinetic nuclear movements generate cell area dynamics that drive extensive cell rearrangements. Over time, the cell proliferation rate declines, effectively solidifying the tissue. Thus, unlike well-studied jamming transitions, the solidification uncovered here resembles a glass transition that depends on the dynamical stresses generated by proliferation and differentiation. Our finding that the fluidity of developing epithelia is linked to interkinetic nuclear movements and the dynamics of growth is likely to be relevant to multiple developing tissues.

Combined computational modeling and experimental analysis integrating chemical and mechanical signals suggests possible mechanism of shoot meristem maintenance

Stem cell maintenance in multilayered shoot apical meristems (SAMs) of plants requires strict regulation of cell growth and division. Exactly how the complex milieu of chemical and mechanical signals interact in the central region of the SAM to regulate cell division plane orientation is not well understood. In this paper, simulations using a newly developed multiscale computational model are combined with experimental studies to suggest and test three hypothesized mechanisms for the regulation of cell division plane orientation and the direction of anisotropic cell expansion in the corpus. Simulations predict that in the Apical corpus, WUSCHEL and cytokinin regulate the direction of anisotropic cell expansion, and cells divide according to tensile stress on the cell wall. In the Basal corpus, model simulations suggest dual roles for WUSCHEL and cytokinin in regulating both the direction of anisotropic cell expansion and cell division plane orientation. Simulation results are followed by a detailed analysis of changes in cell characteristics upon manipulation of WUSCHEL and cytokinin in experiments that support model predictions. Moreover, simulations predict that this layer-specific mechanism maintains both the experimentally observed shape and structure of the SAM as well as the distribution of WUSCHEL in the tissue. This provides an additional link between the roles of WUSCHEL, cytokinin, and mechanical stress in regulating SAM growth and proper stem cell maintenance in the SAM.

Impact of contact inhibition on collective cell migration and proliferation

Contact inhibition limits migration and proliferation of cells in cell colonies. We consider a multiphase field model to investigate the growth dynamics of a cell colony, composed of proliferating cells. The model takes into account the mechanism of contact inhibition of proliferation by local mechanical interactions. We compare nonmigrating and migrating cells, in order to provide a quantitative characterization of the dynamics and analyze the velocity of the colony boundary for both cases. Additionally, we measure single cell velocities, number of neighbor distributions, as well as the influence of stress and age on positions of the cells and with respect to each other.

Adhesion-regulated junction slippage controls cell intercalation dynamics in an Apposed-Cortex Adhesion Model

Cell intercalation is a key cell behaviour of morphogenesis and wound healing, where local cell neighbour exchanges can cause dramatic tissue deformations such as body axis extension. Substantial experimental work has identified the key molecular players facilitating intercalation, but there remains a lack of consensus and understanding of their physical roles. Existing biophysical models that represent cell-cell contacts with single edges cannot study cell neighbour exchange as a continuous process, where neighbouring cell cortices must uncouple. Here, we develop an Apposed-Cortex Adhesion Model (ACAM) to understand active cell intercalation behaviours in the context of a 2D epithelial tissue. The junctional actomyosin cortex of every cell is modelled as a continuous viscoelastic rope-loop, explicitly representing cortices facing each other at bicellular junctions and the adhesion molecules that couple them. The model parameters relate directly to the properties of the key subcellular players that drive dynamics, providing a multi-scale understanding of cell behaviours. We show that active cell neighbour exchanges can be driven by purely junctional mechanisms. Active contractility and cortical turnover in a single bicellular junction are sufficient to shrink and remove a junction. Next, a new, orthogonal junction extends passively. The ACAM reveals how the turnover of adhesion molecules regulates tension transmission and junction deformation rates by controlling slippage between apposed cell cortices. The model additionally predicts that rosettes, which form when a vertex becomes common to many cells, are more likely to occur in actively intercalating tissues with strong friction from adhesion molecules.

Seven challenges in the multiscale modeling of multicellular tissues

The growth and dynamics of multicellular tissues involve tightly regulated and coordinated morphogenetic cell behaviors, such as shape changes, movement, and division, which are governed by subcellular machinery and involve coupling through short- and long-range signals. A key challenge in the fields of developmental biology, tissue engineering and regenerative medicine is to understand how relationships between scales produce emergent tissue-scale behaviors. Recent advances in molecular biology, live-imaging and ex vivo techniques have revolutionized our ability to study these processes experimentally. To fully leverage these techniques and obtain a more comprehensive understanding of the causal relationships underlying tissue dynamics, computational modeling approaches are increasingly spanning multiple spatial and temporal scales, and are coupling cell shape, growth, mechanics, and signaling. Yet such models remain challenging: modeling at each scale requires different areas of technical skills, while integration across scales necessitates the solution to novel mathematical and computational problems. This review aims to summarize recent progress in multiscale modeling of multicellular tissues and to highlight ongoing challenges associated with the construction, implementation, interrogation, and validation of such models. This article is categorized under: Reproductive System Diseases > Computational Models Metabolic Diseases > Computational Models Cancer > Computational Models.

Pinching and pushing: fold formation in the Drosophila dorsal epidermis

Epithelial folding is a fundamental morphogenetic process that shapes planar epithelial sheets into complex three-dimensional structures. Multiple mechanisms can generate epithelial folds, including apical constriction, which acts locally at the cellular level, differential growth on the tissue scale, or buckling because of compression from neighboring tissues. Here, we investigate the formation of dorsally located epithelial folds at segment boundaries during the late stages of Drosophila embryogenesis. We found that the fold formation at the segment boundaries occurs through the juxtaposition of two key morphogenetic processes: local apical constriction and tissue-level compressive forces from posterior segments. Further, we found that epidermal spreading and fold formation are accompanied by spatiotemporal pulses of Hedgehog (Hh) signaling. A computational model that incorporates the local forces generated from the differential tensions of the apical, basal, and lateral sides of the cell and active forces generated within the whole tissue recapitulates the overall fold formation process in wild-type and Hh overexpression conditions. In sum, this work demonstrates how epithelial folding depends on multiple, separable physical mechanisms to generate the final morphology of the dorsal epidermis. This work illustrates the modularity of morphogenetic unit operations that occur during epithelial morphogenesis.

Cell-substrate mechanics guide collective cell migration through intercellular adhesion: a dynamic finite element cellular model

During the process of tissue formation and regeneration, cells migrate collectively while remaining connected through intercellular adhesions. However, the roles of cell-substrate and cell-cell mechanical interactions in regulating collective cell migration are still unclear. In this study, we employ a newly developed finite element cellular model to study collective cell migration by exploring the effects of mechanical feedback between cell and substrate and mechanical signal transmission between adjacent cells. Our viscoelastic model of cells consists many triangular elements and is of high resolution. Cadherin adhesion between cells is modeled explicitly as linear springs at subcellular level. In addition, we incorporate a mechano-chemical feedback loop between cell-substrate mechanics and Rac-mediated cell protrusion. Our model can reproduce a number of experimentally observed patterns of collective cell migration during wound healing, including cell migration persistence, separation distance between cell pairs and migration direction. Moreover, we demonstrate that cell protrusion determined by the cell-substrate mechanics plays an important role in guiding persistent and oriented collective cell migration. Furthermore, this guidance cue can be maintained and transmitted to submarginal cells of long distance through intercellular adhesions. Our study illustrates that our finite element cellular model can be employed to study broad problems of complex tissue in dynamic changes at subcellular level.

Epithelial organ shape is generated by patterned actomyosin contractility and maintained by the extracellular matrix

Epithelial sheets define organ architecture during development. Here, we employed an iterative multiscale computational modeling and quantitative experimental approach to decouple direct and indirect effects of actomyosin-generated forces, nuclear positioning, extracellular matrix, and cell-cell adhesion in shaping Drosophila wing imaginal discs. Basally generated actomyosin forces generate epithelial bending of the wing disc pouch. Surprisingly, acute pharmacological inhibition of ROCK-driven actomyosin contractility does not impact the maintenance of tissue height or curved shape. Computational simulations show that ECM tautness provides only a minor contribution to modulating tissue shape. Instead, passive ECM pre-strain serves to maintain the shape independent from actomyosin contractility. These results provide general insight into how the subcellular forces are generated and maintained within individual cells to induce tissue curvature. Thus, the results suggest an important design principle of separable contributions from ECM prestrain and actomyosin tension during epithelial organogenesis and homeostasis.

The regulation and maintenance of an organ’s shape is a major outstanding problem in developmental biology. An iterative approach combining multiscale computational modelling and quantitative experimental approach was used to decouple direct and indirect roles of subcellular mechanical forces, nuclear positioning, and extracellular matrix in shaping the major axis of the wing pouch during the larval stage in fruit flies, which serves as a prototypical system for investigating epithelial morphogenesis. The research findings in this paper demonstrate that subcellular mechanical forces can effectively generate the curved tissue profile, while extracellular matrix is necessary for preserving the bent shape even in the absence of subcellular mechanical forces once the shape is generated. The developed integrated multiscale modelling environment can be readily extended to generate and test hypothesized novel mechanisms of developmental dynamics of other systems, including organoids that consist of several cellular and extracellular matrix layers.

The Mechanics of Mitotic Cell Rounding

When animal cells enter mitosis, they round up to become spherical. This shape change is accompanied by changes in mechanical properties. Multiple studies using different measurement methods have revealed that cell surface tension, intracellular pressure and cortical stiffness increase upon entry into mitosis. These cell-scale, biophysical changes are driven by alterations in the composition and architecture of the contractile acto-myosin cortex together with osmotic swelling and enable a mitotic cell to exert force against the environment. When the ability of cells to round is limited, for example by physical confinement, cells suffer severe defects in spindle assembly and cell division. The requirement to push against the environment to create space for spindle formation is especially important for cells dividing in tissues. Here we summarize the evidence and the tools used to show that cells exert rounding forces in mitosis in vitro and in vivo, review the molecular basis for this force generation and discuss its function for ensuring successful cell division in single cells and for cells dividing in normal or diseased tissues.

Mesoscopic dynamic model of epithelial cell division with cell-cell junction effects

Cell division is central for embryonic development, tissue morphogenesis, and tumor growth. Experiments have evidenced that mitotic cell division is manipulated by the intercellular cues such as cell-cell junctions. However, it still remains unclear how these cortical-associated cues mechanically affect the mitotic spindle machinery, which determines the position and orientation of the cell division. In this paper, a mesoscopic dynamic cell division model is established to explore the integrated regulations of cortical polarity, microtubule pulling forces, cell deformability, and internal osmotic pressure. We show that the distributed pulling forces of astral microtubules play a key role in encoding the instructive cortical cues to orient and position the spindle of a dividing cell. The present model can not only predict the spindle orientation and position, but also capture the morphological evolution of cell rounding. The theoretical results agree well with relevant experiments both qualitatively and quantitatively. This work sheds light on the mechanical linkage between cell cortex and mitotic spindle, and holds potential in regulating cell division and sculpting tissue morphology.

Decoding Calcium Signaling Dynamics during Drosophila Wing Disc Development

The robust specification of organ development depends on coordinated cell-cell communication. This process requires signal integration among multiple pathways, relying on second messengers such as calcium ions. Calcium signaling encodes a significant portion of the cellular state by regulating transcription factors, enzymes, and cytoskeletal proteins. However, the relationships between the inputs specifying cell and organ development, calcium signaling dynamics, and final organ morphology are poorly understood. Here, we have designed a quantitative image-analysis pipeline for decoding organ-level calcium signaling. With this pipeline, we extracted spatiotemporal features of calcium signaling dynamics during the development of the Drosophila larval wing disc, a genetic model for organogenesis. We identified specific classes of wing phenotypes that resulted from calcium signaling pathway perturbations, including defects in gross morphology, vein differentiation, and overall size. We found four qualitative classes of calcium signaling activity. These classes can be ordered based on agonist stimulation strength Gαq-mediated signaling. In vivo calcium signaling dynamics depend on both receptor tyrosine kinase/phospholipase C γ and G protein-coupled receptor/phospholipase C β activities. We found that spatially patterned calcium dynamics correlate with known differential growth rates between anterior and posterior compartments. Integrated calcium signaling activity decreases with increasing tissue size, and it responds to morphogenetic perturbations that impact organ growth. Together, these findings define how calcium signaling dynamics integrate upstream inputs to mediate multiple response outputs in developing epithelial organs.

Cell-Based Model of the Generation and Maintenance of the Shape and Structure of the Multilayered Shoot Apical Meristem of Arabidopsis thaliana

One of the central problems in animal and plant developmental biology is deciphering how chemical and mechanical signals interact within a tissue to produce organs of defined size, shape, and function. Cell walls in plants impose a unique constraint on cell expansion since cells are under turgor pressure and do not move relative to one another. Cell wall extensibility and constantly changing distribution of stress on the wall are mechanical properties that vary between individual cells and contribute to rates of expansion and orientation of cell division. How exactly cell wall mechanical properties influence cell behavior is still largely unknown. To address this problem, a novel, subcellular element computational model of growth of stem cells within the multilayered shoot apical meristem (SAM) of Arabidopsis thaliana is developed and calibrated using experimental data. Novel features of the model include separate, detailed descriptions of cell wall extensibility and mechanical stiffness, deformation of the middle lamella, and increase in cytoplasmic pressure generating internal turgor pressure. The model is used to test novel hypothesized mechanisms of formation of the shape and structure of the growing, multilayered SAM based on WUS concentration of individual cells controlling cell growth rates and layer-dependent anisotropic mechanical properties of subcellular components of individual cells determining anisotropic cell expansion directions. Model simulations also provide a detailed prediction of distribution of stresses in the growing tissue which can be tested in future experiments.

Reverse-engineering organogenesis through feedback loops between model systems

Biological complexity and ethical limitations necessitate models of human development. Traditionally, genetic model systems have provided inexpensive routes to define mechanisms governing organ development. Recent progress has led to 3D human organoid models of development and disease. However, robust methods to control the size and morphology of organoids for high throughput studies need to be developed. Additionally, insights from multiple developmental contexts are required to reveal conserved genes and processes regulating organ growth and development. Positive feedback between quantitative studies using mammalian organoids and insect micro-organs enable identification of underlying principles for organ size and shape control. Advances in the field of multicellular systems engineering are enabling unprecedented high-content studies in developmental biology and disease modeling. These will lead to fundamental advances in regenerative medicine and tissue-engineered soft robotics.