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

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.

The adhesion GPCR and PCP component flamingo (FMI-1) alters body size and regulates the composition of the extracellular matrix

The extracellular matrix (ECM) is a network of macromolecules that presents a vital scaffold for cells and enables multiple ways of cellular communication. Thus, it is essential for many physiological processes such as development, tissue morphogenesis, homeostasis, the shape and partially the size of the body and its organs. To ensure these, the composition of the ECM is tissue-specific and highly dynamic. ECM homeostasis is therefore tightly controlled by several mechanisms. Here, we show that FMI-1, the homolog of the Adhesion GPCR Flamingo/CELSR/ADGRC in the nematode Caenorhabditis elegans, modulates the composition of the ECM by controlling the production both of ECM molecules such as collagens and also of ECM modifying enzymes. Thereby, FMI-1 affects the morphology and functionality of the nematode´s cuticle, which is mainly composed of ECM, and also modulates the body size. Mechanistic analyses highlight the fact that FMI-1 exerts its function from neurons non-cell autonomously (trans) solely via its extracellular N terminus. Our data support a model, by which the activity of the receptor, which has a well-described role in the planar cell polarity (PCP) pathway, involves the PCP molecule VANG-1, but seems to be independent of the DBL-1/BMP pathway.

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.

Modeling the extracellular matrix in cell migration and morphogenesis: a guide for the curious biologist

The extracellular matrix (ECM) is a highly complex structure through which biochemical and mechanical signals are transmitted. In processes of cell migration, the ECM also acts as a scaffold, providing structural support to cells as well as points of potential attachment. Although the ECM is a well-studied structure, its role in many biological processes remains difficult to investigate comprehensively due to its complexity and structural variation within an organism. In tandem with experiments, mathematical models are helpful in refining and testing hypotheses, generating predictions, and exploring conditions outside the scope of experiments. Such models can be combined and calibrated with in vivo and in vitro data to identify critical cell-ECM interactions that drive developmental and homeostatic processes, or the progression of diseases. In this review, we focus on mathematical and computational models of the ECM in processes such as cell migration including cancer metastasis, and in tissue structure and morphogenesis. By highlighting the predictive power of these models, we aim to help bridge the gap between experimental and computational approaches to studying the ECM and to provide guidance on selecting an appropriate model framework to complement corresponding experimental studies.

Combined computational modeling and experimental study of the biomechanical mechanisms of platelet-driven contraction of fibrin clots

While blood clot formation has been relatively well studied, little is known about the mechanisms underlying the subsequent structural and mechanical clot remodeling called contraction or retraction. Impairment of the clot contraction process is associated with both life-threatening bleeding and thrombotic conditions, such as ischemic stroke, venous thromboembolism, and others. Recently, blood clot contraction was observed to be hindered in patients with COVID-19. A three-dimensional multiscale computational model is developed and used to quantify biomechanical mechanisms of the kinetics of clot contraction driven by platelet-fibrin pulling interactions. These results provide important biological insights into contraction of platelet filopodia, the mechanically active thin protrusions of the plasma membrane, described previously as performing mostly a sensory function. The biomechanical mechanisms and modeling approach described can potentially apply to studying other systems in which cells are embedded in a filamentous network and exert forces on the extracellular matrix modulated by the substrate stiffness.

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.

Mechanobiological implications of age-related remodelling in the outer retina

The outer retina consists of the light-sensitive photoreceptors, the pigmented epithelium, and the choroid, which interact in a complex manner to sustain homeostasis. The organisation and function of these cellular layers are mediated by the extracellular matrix compartment named Bruch's membrane, situated between the retinal epithelium and the choroid. Like many tissues, the retina experiences age-related structural and metabolic changes, which are relevant for understanding major blinding diseases of the elderly, such as age-related macular degeneration. Compared with other tissues, the retina mainly comprises postmitotic cells, making it less able to maintain its mechanical homeostasis over the years functionally. Aspects of retinal ageing, like the structural and morphometric changes of the pigment epithelium and the heterogenous remodelling of the Bruch's membrane, imply changes in tissue mechanics and may affect functional integrity. In recent years, findings in the field of mechanobiology and bioengineering highlighted the importance of mechanical changes in tissues for understanding physiological and pathological processes. Here, we review the current knowledge of age-related changes in the outer retina from a mechanobiological perspective, aiming to generate food for thought for future mechanobiology studies in the outer retina.

Growth anisotropy of the extracellular matrix shapes a developing organ

Final organ size and shape result from volume expansion by growth and shape changes by contractility. Complex morphologies can also arise from differences in growth rate between tissues. We address here how differential growth guides the morphogenesis of the growing Drosophila wing imaginal disc. We report that 3D morphology results from elastic deformation due to differential growth anisotropy between the epithelial cell layer and its enveloping extracellular matrix (ECM). While the tissue layer grows in plane, growth of the bottom ECM occurs in 3D and is reduced in magnitude, thereby causing geometric frustration and tissue bending. The elasticity, growth anisotropy and morphogenesis of the organ are fully captured by a mechanical bilayer model. Moreover, differential expression of the Matrix metalloproteinase MMP2 controls growth anisotropy of the ECM envelope. This study shows that the ECM is a controllable mechanical constraint whose intrinsic growth anisotropy directs tissue morphogenesis in a developing organ.

Bi-directional crosstalk between cells and extracellular matrix leads to network morphogenesis in multi-layered tissues

Cell-generated mechanical forces drive many cellular and tissue-level movements and rearrangements required for the tissue or organ to develop its shape1, 2, 3, 4, 5. The prevalent view of tissue morphogenesis relies on epithelial folding resulting in compressed epithelial monolayers, overlooking the involvement of stroma in morphogenesis1, 4, 6, 7. Here, we report a giant web-like network formation of stromal cells in the epithelium-stroma interface, resulting from a multi-scale mechano-reciprocity between migrating cells and their extracellular environment. In multi-layered tissues, surface wrinkles form by a stromal cell-mediated tensional force exerted at the basement membrane. The topographical cue is transmitted to the stromal cell, directing its protrusion and migration along the wrinkles. This inductive movement of the cells conveys traction forces to its surrounding extracellular matrix, remodeling the local architectures of the stroma. In this manner, stromal cells and wrinkles communicate recursively to generate the cellular network. Our observation provides a rational mechanism for network formation in living tissues and a new understanding of the role of cellular-level tensional force in morphogenesis.

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.

The wing imaginal disc

The Drosophila wing imaginal disc is a tissue of undifferentiated cells that are precursors of the wing and most of the notum of the adult fly. The wing disc first forms during embryogenesis from a cluster of ∼30 cells located in the second thoracic segment, which invaginate to form a sac-like structure. They undergo extensive proliferation during larval stages to form a mature larval wing disc of ∼35,000 cells. During this time, distinct cell fates are assigned to different regions, and the wing disc develops a complex morphology. Finally, during pupal stages the wing disc undergoes morphogenetic processes and then differentiates to form the adult wing and notum. While the bulk of the wing disc comprises epithelial cells, it also includes neurons and glia, and is associated with tracheal cells and muscle precursor cells. The relative simplicity and accessibility of the wing disc, combined with the wealth of genetic tools available in Drosophila, have combined to make it a premier system for identifying genes and deciphering systems that play crucial roles in animal development. Studies in wing imaginal discs have made key contributions to many areas of biology, including tissue patterning, signal transduction, growth control, regeneration, planar cell polarity, morphogenesis, and tissue mechanics.

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.

Dynamic changes in epithelial cell packing during tissue morphogenesis

Cell packing - the spatial arrangement of cells - determines the shapes of organs. Recently, investigations of organ development in a variety of model organisms have uncovered cellular mechanisms that are used by epithelial tissues to change cell packing, and thereby their shapes, to generate functional architectures. Here, we review these cellular mechanisms across a wide variety of developmental processes in vertebrates and invertebrates and identify a set of common motifs in the morphogenesis toolbox that, in combination, appear to allow any change in tissue shape. We focus on tissue elongation, folding and invagination, and branching. We also highlight how these morphogenetic processes are achieved by cell-shape changes, cell rearrangements, and oriented cell division. Finally, we describe approaches that have the potential to engineer three-dimensional tissues for both basic science and translational purposes. This review provides a framework for future analyses of how tissues are shaped by the dynamics of epithelial cell packing.

The origin and the mechanism of mechanical polarity during epithelial folding

Epithelial tissues are sheet-like tissue structures that line the inner and outer surfaces of animal bodies and organs. Their remarkable ability to actively produce, or passively adapt to, complex surface geometries has fascinated physicists and biologists alike for centuries. The most simple and yet versatile process of epithelial deformation is epithelial folding, through which curved shapes, tissue convolutions and internal structures are produced. The advent of quantitative live imaging, combined with experimental manipulation and computational modeling, has rapidly advanced our understanding of epithelial folding. In particular, a set of mechanical principles has emerged to illustrate how forces are generated and dissipated to instigate curvature transitions in a variety of developmental contexts. Folding a tissue requires that mechanical loads or geometric changes be non-uniform. Given that polarity is the most distinct and fundamental feature of epithelia, understanding epithelial folding mechanics hinges crucially on how forces become polarized and how polarized differential deformation arises, for which I coin the term 'mechanical polarity'. In this review, five typical modules of mechanical processes are distilled from a diverse array of epithelial folding events. Their mechanical underpinnings with regard to how forces and polarity intersect are analyzed to accentuate the importance of mechanical polarity in the understanding of epithelial folding.

A multiscale model via single-cell transcriptomics reveals robust patterning mechanisms during early mammalian embryo development

During early mammalian embryo development, a small number of cells make robust fate decisions at particular spatial locations in a tight time window to form inner cell mass (ICM), and later epiblast (Epi) and primitive endoderm (PE). While recent single-cell transcriptomics data allows scrutinization of heterogeneity of individual cells, consistent spatial and temporal mechanisms the early embryo utilize to robustly form the Epi/PE layers from ICM remain elusive. Here we build a multiscale three-dimensional model for mammalian embryo to recapitulate the observed patterning process from zygote to late blastocyst. By integrating the spatiotemporal information reconstructed from multiple single-cell transcriptomic datasets, the data-informed modeling analysis suggests two major processes critical to the formation of Epi/PE layers: a selective cell-cell adhesion mechanism (via EphA4/EphrinB2) for fate-location coordination and a temporal attenuation mechanism of cell signaling (via Fgf). Spatial imaging data and distinct subsets of single-cell gene expression data are then used to validate the predictions. Together, our study provides a multiscale framework that incorporates single-cell gene expression datasets to analyze gene regulations, cell-cell communications, and physical interactions among cells in complex geometries at single-cell resolution, with direct application to late-stage development of embryogenesis.

To form and function: on the role of basement membrane mechanics in tissue development, homeostasis and disease

The basement membrane (BM) is a special type of extracellular matrix that lines the basal side of epithelial and endothelial tissues. Functionally, the BM is important for providing physical and biochemical cues to the overlying cells, sculpting the tissue into its correct size and shape. In this review, we focus on recent studies that have unveiled the complex mechanical properties of the BM. We discuss how these properties can change during development, homeostasis and disease via different molecular mechanisms, and the subsequent impact on tissue form and function in a variety of organisms. We also explore how better characterization of BM mechanics can contribute to disease diagnosis and treatment, as well as development of better in silico and in vitro models that not only impact the fields of tissue engineering and regenerative medicine, but can also reduce the use of animals in research.