Sizing it up: The mechanical feedback hypothesis of organ growth regulation

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 Interaction of Mechanics and the Hippo Pathway in Drosophila melanogaster

Drosophila melanogaster has emerged as an ideal system for studying the networks that control tissue development and homeostasis and, given the similarity of the pathways involved, controlled and uncontrolled growth in mammalian systems. The signaling pathways used in patterning the Drosophila wing disc are well known and result in the emergence of interaction of these pathways with the Hippo signaling pathway, which plays a central role in controlling cell proliferation and apoptosis. Mechanical effects are another major factor in the control of growth, but far less is known about how they exert their control. Herein, we develop a mathematical model that integrates the mechanical interactions between cells, which occur via adherens and tight junctions, with the intracellular actin network and the Hippo pathway so as to better understand cell-autonomous and non-autonomous control of growth in response to mechanical forces.

Transducing compressive forces into cellular outputs in cancer and beyond

In living organisms, cells sense mechanical forces (shearing, tensile, and compressive) and respond to those physical cues through a process called mechanotransduction. This process includes the simultaneous activation of biochemical signaling pathways. Recent studies mostly on human cells revealed that compressive forces selectively modulate a wide range of cell behavior, both in compressed and in neighboring less compressed cells. Besides participating in tissue homeostasis such as bone healing, compression is also involved in pathologies, including intervertebral disc degeneration or solid cancers. In this review, we will summarize the current scattered knowledge of compression-induced cell signaling pathways and their subsequent cellular outputs, both in physiological and pathological conditions, such as solid cancers.

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.

The multimodal action of G alpha q in coordinating growth and homeostasis in the Drosophila wing imaginal disc

Background

G proteins mediate cell responses to various ligands and play key roles in organ development. Dysregulation of G-proteins or Ca ²⁺ signaling impacts many human diseases and results in birth defects. However, the downstream effectors of specific G proteins in developmental regulatory networks are still poorly understood.

Methods

We employed the Gal4/UAS binary system to inhibit or overexpress Gαq in the wing disc, followed by phenotypic analysis. Immunohistochemistry and next-gen RNA sequencing identified the downstream effectors and the signaling cascades affected by the disruption of Gαq homeostasis.

Results

Here, we characterized how the G protein subunit Gαq tunes the size and shape of the wing in the larval and adult stages of development. Downregulation of Gαq in the wing disc reduced wing growth and delayed larval development. Gαq overexpression is sufficient to promote global Ca ²⁺ waves in the wing disc with a concomitant reduction in the Drosophila final wing size and a delay in pupariation. The reduced wing size phenotype is further enhanced when downregulating downstream components of the core Ca ²⁺ signaling toolkit, suggesting that downstream Ca ²⁺ signaling partially ameliorates the reduction in wing size. In contrast, Gαq -mediated pupariation delay is rescued by inhibition of IP ₃ R, a key regulator of Ca ²⁺ signaling. This suggests that Gαq regulates developmental phenotypes through both Ca ²⁺ -dependent and Ca ²⁺ -independent mechanisms. RNA seq analysis shows that disruption of Gαq homeostasis affects nuclear hormone receptors, JAK/STAT pathway, and immune response genes. Notably, disruption of Gαq homeostasis increases expression levels of Dilp8, a key regulator of growth and pupariation timing.

Conclusion

Gαq activity contributes to cell size regulation and wing metamorphosis. Disruption to Gαq homeostasis in the peripheral wing disc organ delays larval development through ecdysone signaling inhibition. Overall, Gαq signaling mediates key modules of organ size regulation and epithelial homeostasis through the dual action of Ca ²⁺ -dependent and independent mechanisms.

MAPPER: An Open-Source, High-Dimensional Image Analysis Pipeline Unmasks Differential Regulation of Drosophila Wing Features

Phenomics requires quantification of large volumes of image data, necessitating high throughput image processing approaches. Existing image processing pipelines for Drosophila wings, a powerful genetic model for studying the underlying genetics for a broad range of cellular and developmental processes, are limited in speed, precision, and functional versatility. To expand on the utility of the wing as a phenotypic screening system, we developed MAPPER, an automated machine learning-based pipeline that quantifies high-dimensional phenotypic signatures, with each dimension quantifying a unique morphological feature of the Drosophila wing. MAPPER magnifies the power of Drosophila phenomics by rapidly quantifying subtle phenotypic differences in sample populations. We benchmarked MAPPER’s accuracy and precision in replicating manual measurements to demonstrate its widespread utility. The morphological features extracted using MAPPER reveal variable sexual dimorphism across Drosophila species and unique underlying sex-specific differences in morphogen signaling in male and female wings. Moreover, the length of the proximal-distal axis across the species and sexes shows a conserved scaling relationship with respect to the wing size. In sum, MAPPER is an open-source tool for rapid, high-dimensional analysis of large imaging datasets. These high-content phenomic capabilities enable rigorous and systematic identification of genotype-to-phenotype relationships in a broad range of screening and drug testing applications and amplify the potential power of multimodal genomic approaches.

From spikes to intercellular waves: Tuning intercellular calcium signaling dynamics modulates organ size control

Information flow within and between cells depends significantly on calcium (Ca²⁺) signaling dynamics. However, the biophysical mechanisms that govern emergent patterns of Ca²⁺ signaling dynamics at the organ level remain elusive. Recent experimental studies in developing Drosophila wing imaginal discs demonstrate the emergence of four distinct patterns of Ca²⁺ activity: Ca²⁺ spikes, intercellular Ca²⁺ transients, tissue-level Ca²⁺ waves, and a global “fluttering” state. Here, we used a combination of computational modeling and experimental approaches to identify two different populations of cells within tissues that are connected by gap junction proteins. We term these two subpopulations “initiator cells,” defined by elevated levels of Phospholipase C (PLC) activity, and “standby cells,” which exhibit baseline activity. We found that the type and strength of hormonal stimulation and extent of gap junctional communication jointly determine the predominate class of Ca²⁺ signaling activity. Further, single-cell Ca²⁺ spikes are stimulated by insulin, while intercellular Ca²⁺ waves depend on Gαq activity. Our computational model successfully reproduces how the dynamics of Ca²⁺ transients varies during organ growth. Phenotypic analysis of perturbations to Gαq and insulin signaling support an integrated model of cytoplasmic Ca²⁺ as a dynamic reporter of overall tissue growth. Further, we show that perturbations to Ca²⁺ signaling tune the final size of organs. This work provides a platform to further study how organ size regulation emerges from the crosstalk between biochemical growth signals and heterogeneous cell signaling states.

Calcium (Ca²⁺) is a universal second messenger that regulates a myriad of cellular processes such as cell division, cell proliferation and apoptosis. Multiple patterns of Ca²⁺ signaling including single-cell spikes, multicellular Ca²⁺ transients, large-scale Ca²⁺ waves, and global “fluttering” have been observed in epithelial systems during organ development. Key molecular players and biophysical mechanisms involved in formation of these patterns during organ development are not well understood. In this work, we developed a generalized multicellular model of Ca²⁺ that captures all the key categories of Ca²⁺ activity as a function of key hormonal signals. Integration of model predictions and experiments reveals two subclasses of cell populations and demonstrates that Ca²⁺ signaling activity at the organ scale is defined by a general decrease in gap junction communication as an organ grows. Our experiments also reveal that a “goldilocks zone” of optimal Ca²⁺ activity is required to achieve optimal growth at the organ level.

Hydrodynamic stress and phenotypic plasticity of the zebrafish regenerating fin

Understanding how extrinsic factors modulate genetically encoded information to produce a specific phenotype is of prime scientific interest. In particular, the feedback mechanism between abiotic forces and locomotory organs during morphogenesis to achieve efficient movement is a highly relevant example of such modulation. The study of this developmental process can provide unique insights on the transduction of cues at the interface between physics and biology. Here, we take advantage of the natural ability of adult zebrafish to regenerate their amputated fins to assess its morphogenic plasticity upon external modulations. Using a variety of surgical and chemical treatments, we could induce phenotypic responses to the structure of the fin. Through the ablation of specific rays in regenerating caudal fins, we generated artificially narrowed appendages in which the fin cleft depth and the positioning of rays bifurcations were perturbed compared with normal regenerates. To dissect the role of mechanotransduction in this process, we investigated the patterns of hydrodynamic forces acting on the surface of a zebrafish fin during regeneration by using particle tracking velocimetry on a range of biomimetic hydrofoils. This experimental approach enabled us to quantitatively compare hydrodynamic stress distributions over flapping fins of varying sizes and shapes. As a result, viscous shear stress acting on the distal margin of regenerating fins and the resulting internal tension are proposed as suitable signals for guiding the regulation of ray growth dynamics and branching pattern. Our findings suggest that mechanical forces are involved in the fine-tuning of the locomotory organ during fin morphogenesis.

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.

Growth control in the Drosophila wing disk

The regulation of size and shape is a fundamental requirement of biological development and has been a subject of scientific study for centuries, but we still lack an understanding of how organisms know when to stop growing. Imaginal wing disks of the fruit fly Drosophila melanogaster, which are precursors of the adult wings, are an archetypal tissue for studying growth control. The growth of the disks is dependent on many inter- and intra-organ factors such as morphogens, mechanical forces, nutrient levels, and hormones that influence gene expression and cell growth. Extracellular signals are transduced into gene-control signals via complex signal transduction networks, and since cells typically receive many different signals, a mechanism for integrating the signals is needed. Our understanding of the effect of morphogens on tissue-level growth regulation via individual pathways has increased significantly in the last half century, but our understanding of how multiple biochemical and mechanical signals are integrated to determine whether or not a cell decides to divide is still rudimentary. Numerous fundamental questions are involved in understanding the decision-making process, and here we review the major biochemical and mechanical pathways involved in disk development with a view toward providing a basis for beginning to understand how multiple signals can be integrated at the cell level, and how this translates into growth control at the level of the imaginal disk. This article is categorized under: Analytical and Computational Methods > Computational Methods Biological Mechanisms > Cell Signaling Models of Systems Properties and Processes > Cellular Models.

Gene expression atlas of a developing tissue by single cell expression correlation analysis

The Drosophila wing disc has been a fundamental model system for the discovery of key signaling pathways and for our understanding of developmental processes. However, a complete map of gene expression in this tissue is lacking. To obtain a complete gene expression atlas in the wing disc, we employed single-cell sequencing (scRNA-seq) and developed a new method for analyzing scRNA-seq data based on gene expression correlations rather than cell mapping. This enables us to compute expression maps for all detected genes in the wing disc and to discover 824 genes with spatially restricted expression patterns. This approach identifies both known and new clusters of genes with similar expression patterns and functional relevance. As proof of concept, we characterize the previously unstudied gene CG5151 and show that it regulates Wnt signaling. This novel method will enable the leveraging of scRNA-seq data for generating expression atlases of undifferentiated tissues during development.

Retinal stem cells modulate proliferative parameters to coordinate post-embryonic morphogenesis in the eye of fish

Combining clonal analysis with a computational agent based model, we investigate how tissue-specific stem cells for neural retina (NR) and retinal pigmented epithelium (RPE) of the teleost medaka (Oryzias latipes) coordinate their growth rates. NR cell division timing is less variable, consistent with an upstream role as growth inducer. RPE cells divide with greater variability, consistent with a downstream role responding to inductive signals. Strikingly, the arrangement of the retinal ciliary marginal zone niche results in a spatially biased random lineage loss, where stem- and progenitor cell domains emerge spontaneously. Further, our data indicate that NR cells orient division axes to regulate organ shape and retinal topology. We highlight an unappreciated mechanism for growth coordination, where one tissue integrates cues to synchronize growth of nearby tissues. This strategy may enable evolution to modulate cell proliferation parameters in one tissue to adapt whole-organ morphogenesis in a complex vertebrate organ.

By the time babies reach adulthood, they have grown many times larger than they were at birth. This development is driven by an increase in the number and size of cells in the body. In particular, special types of cells, called stem cells, act as a reservoir for tissues: they divide to create new cells that will mature into various specialized structures.

The retina is the light-sensitive part of the eye. It consists of the neural retina, a tissue that contains light-detecting cells, which is supported by the retinal pigment epithelium or RPE. In fish, the RPE and neural retina are replenished by distinct groups of stem cells that do not mix, despite the tissues being close together.

Unlike humans, fish grow throughout adulthood, and their eyes must then keep pace with the body. This means that the different tissues in the retina must somehow coordinate to expand at the same rate: otherwise, the retina would get wrinkled and not work properly. Tsingos et al. therefore wanted to determine how stem cells in the neural retina and RPE co-operated to produce the right number of new cells at the right time.

First, stem cells in the eyes of newly hatched fish were labelled with a visible marker so that their divisions could be tracked over time to build cell family trees. This showed that stem cells behaved differently in the neural retina and the RPE. Computer simulations of the growing retina explained this behavior: stem cells in the neural retina were telling the RPE stem cells when it was time to divide. Combining results from the simulations with data from the experiments revealed that a stem cell decided to keep up dividing partly because of its position in the tissue, and partly because of random chance.

To be healthy, the body needs to fine-tune the number of cells it produces: creating too few cells may make it difficult to heal after injury, but making too many could lead to diseases such as cancer. Understanding how tissues normally agree to grow together could therefore open new avenues of treatment for these conditions.

Statistics of noisy growth with mechanical feedback in elastic tissues

Tissue growth is a fundamental aspect of development and is intrinsically noisy. Stochasticity has important implications for morphogenesis, precise control of organ size, and regulation of tissue composition and heterogeneity. However, the basic statistical properties of growing tissues, particularly when growth induces mechanical stresses that can in turn affect growth rates, have received little attention. Here, we study the noisy growth of elastic sheets subject to mechanical feedback. Considering both isotropic and anisotropic growth, we find that the density-density correlation function shows power law scaling. We also consider the dynamics of marked, neutral clones of cells. We find that the areas (but not the shapes) of two clones are always statistically independent, even when they are adjacent. For anisotropic growth, we show that clone size variance scales like the average area squared and that the mode amplitudes characterizing clone shape show a slow [Formula: see text] decay, where n is the mode index. This is in stark contrast to the isotropic case, where relative variations in clone size and shape vanish at long times. The high variability in clone statistics observed in anisotropic growth is due to the presence of two soft modes-growth modes that generate no stress. Our results lay the groundwork for more in-depth explorations of the properties of noisy tissue growth in specific biological contexts.

Microfluidics on the fly: Inexpensive rapid fabrication of thermally laminated microfluidic devices for live imaging and multimodal perturbations of multicellular systems

Microfluidic devices provide a platform for analyzing both natural and synthetic multicellular systems. Currently, substantial capital investment and expertise are required for creating microfluidic devices using standard soft-lithography. These requirements present barriers to entry for many nontraditional users of microfluidics, including developmental biology laboratories. Therefore, fabrication methodologies that enable rapid device iteration and work "out-of-the-box" can accelerate the integration of microfluidics with developmental biology. Here, we have created and characterized low-cost hybrid polyethylene terephthalate laminate (PETL) microfluidic devices that are suitable for cell and micro-organ culture assays. These devices were validated with mammalian cell lines and the Drosophila wing imaginal disc as a model micro-organ. First, we developed and tested PETLs that are compatible with both long-term cultures and high-resolution imaging of cells and organs. Further, we achieved spatiotemporal control of chemical gradients across the wing discs with a multilayered microfluidic device. Finally, we created a multilayered device that enables controllable mechanical loading of micro-organs. This mechanical actuation assay was used to characterize the response of larval wing discs at different developmental stages. Interestingly, increased deformation of the older wing discs for the same mechanical loading suggests that the compliance of the organ is increased in preparation for subsequent morphogenesis. Together, these results demonstrate the applicability of hybrid PETL devices for biochemical and mechanobiology studies on micro-organs and provide new insights into the mechanics of organ development.

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.

Systemic signalling and local effectors in developmental stability, body symmetry, and size

Symmetric growth and the origins of fluctuating asymmetry are unresolved phenomena of biology. Small, and sometimes noticeable, deviations from perfect bilateral symmetry reflect the vulnerability of development to perturbations. The degree of asymmetry is related to the magnitude of the perturbations and the ability of an individual to cope with them. As the left and right sides of an individual were presumed to be genetically identical, deviations of symmetry were traditionally attributed to non-genetic effects such as environmental and developmental noise. In this review, we draw attention to other possible sources of variability, especially to somatic mutations and transposons. Mutations are a major source of phenotypic variability and recent genomic data have highlighted somatic mutations as ubiquitous, even in phenotypically normal individuals. We discuss the importance of factors that are responsible for buffering and stabilizing the genome and for maintaining size robustness and quality through elimination of less-fit or damaged cells. However, the important question that arises from these studies is whether this self-correcting capacity and intrinsic organ size controls are sufficient to explain how symmetric structures can reach an identical size and shape. Indeed, recent discoveries in the fruit fly have uncovered a conserved hormone of the insulin/IGF/relaxin family, Dilp8, that is responsible for stabilizing body size and symmetry in the face of growth perturbations. Dilp8 alarm signals periphery growth status to the brain, where it acts on its receptor Lgr3. Loss of Dilp8-Lgr3 signaling renders flies incapable of detecting growth perturbations and thus maintaining a stable size and symmetry. These findings help to understand how size and symmetry of somatic tissues remain undeterred in noisy environments, after injury or illnesses, and in the presence of accumulated somatic mutations.

AIP1 and cofilin ensure a resistance to tissue tension and promote directional cell rearrangement

In order to understand how tissue mechanics shapes animal body, it is critical to clarify how cells respond to and resist tissue stress when undergoing morphogenetic processes, such as cell rearrangement. Here, we address the question in the Drosophila wing epithelium, where anisotropic tissue tension orients cell rearrangements. We found that anisotropic tissue tension localizes actin interacting protein 1 (AIP1), a cofactor of cofilin, on the remodeling junction via cooperative binding of cofilin to F-actin. AIP1 and cofilin promote actin turnover and locally regulate the Canoe-mediated linkage between actomyosin and the junction. This mechanism is essential for cells to resist the mechanical load imposed on the remodeling junction perpendicular to the direction of tissue stretching. Thus, the present study delineates how AIP1 and cofilin achieve an optimal balance between resistance to tissue tension and morphogenesis.

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.

Basement Membrane Manipulation in Drosophila Wing Discs Affects Dpp Retention but Not Growth Mechanoregulation

Basement membranes (BMs) are extracellular matrix polymers basally underlying epithelia, where they regulate cell signaling and tissue mechanics. Constriction by the BM shapes Drosophila wing discs, a well-characterized model of tissue growth. Recently, the hypothesis that mechanical factors govern wing growth has received much attention, but it has not been definitively tested. In this study, we manipulated BM composition to cause dramatic changes in tissue tension. We found that increased tissue compression when perlecan was knocked down did not affect adult wing size. BM elimination, decreasing compression, reduced wing size but did not visibly affect Hippo signaling, widely postulated to mediate growth mechanoregulation. BM elimination, in contrast, attenuated signaling by bone morphogenetic protein/transforming growth factor β ligand Dpp, which was not efficiently retained within the tissue and escaped to the body cavity. Our results challenge mechanoregulation of wing growth, while uncovering a function of BMs in preserving a growth-promoting tissue environment.

Release of Applied Mechanical Loading Stimulates Intercellular Calcium Waves in Drosophila Wing Discs

Mechanical forces are critical but poorly understood inputs for organogenesis and wound healing. Calcium ions (Ca2+) are critical second messengers in cells for integrating environmental and mechanical cues, but the regulation of Ca2+ signaling is poorly understood in developing epithelial tissues. Here we report a chip-based regulated environment for microorgans that enables systematic investigations of the crosstalk between an organ's mechanical stress environment and biochemical signaling under genetic and chemical perturbations. This method enabled us to define the essential conditions for generating organ-scale intercellular Ca2+ waves in Drosophila wing discs that are also observed in vivo during organ development. We discovered that mechanically induced intercellular Ca2+ waves require fly extract growth serum as a chemical stimulus. Using the chip-based regulated environment for microorgans, we demonstrate that not the initial application but instead the release of mechanical loading is sufficient, but not necessary, to initiate intercellular Ca2+ waves. The Ca2+ response depends on the prestress intercellular Ca2+ activity and not on the magnitude or duration of the mechanical stimulation applied. Mechanically induced intercellular Ca2+ waves rely on IP3R-mediated Ca2+-induced Ca2+ release and propagation through gap junctions. Thus, intercellular Ca2+ waves in developing epithelia may be a consequence of stress dissipation during organ growth.

Oxygenation and adenosine deaminase support growth and proliferation of ex vivo cultured Drosophila wing imaginal discs

The Drosophila wing imaginal disc has been an important model system over past decades for discovering novel biology related to development, signaling and epithelial morphogenesis. Novel experimental approaches have been enabled using a culture setup that allows ex vivo cultures of wing discs. Current setups, however, are not able to sustain both growth and cell-cycle progression of wing discs ex vivo We discover here a setup that requires both oxygenation of the tissue and adenosine deaminase activity in the medium, and supports both growth and proliferation of wing discs for 9 h. Nonetheless, further work will be required to extend the duration of the culturing and to enable live imaging of the cultured discs in the future.

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

Mitotic rounding during cell division is critical for preventing daughter cells from inheriting an abnormal number of chromosomes, a condition that occurs frequently in cancer cells. Cells must significantly expand their apical area and transition from a polygonal to circular apical shape to achieve robust mitotic rounding in epithelial tissues, which is where most cancers initiate. However, how cells mechanically regulate robust mitotic rounding within packed tissues is unknown. Here, we analyze mitotic rounding using a newly developed multi-scale subcellular element computational model that is calibrated using experimental data. Novel biologically relevant features of the model include separate representations of the sub-cellular components including the apical membrane and cytoplasm of the cell at the tissue scale level as well as detailed description of cell properties during mitotic rounding. Regression analysis of predictive model simulation results reveals the relative contributions of osmotic pressure, cell-cell adhesion and cortical stiffness to mitotic rounding. Mitotic area expansion is largely driven by regulation of cytoplasmic pressure. Surprisingly, mitotic shape roundness within physiological ranges is most sensitive to variation in cell-cell adhesivity and stiffness. An understanding of how perturbed mechanical properties impact mitotic rounding has important potential implications on, amongst others, how tumors progressively become more genetically unstable due to increased chromosomal aneuploidy and more aggressive.

A full lifespan model of vertebrate lens growth

The mathematical determinants of vertebrate organ growth have yet to be elucidated fully. Here, we utilized empirical measurements and a dynamic branching process-based model to examine the growth of a simple organ system, the mouse lens, from E14.5 until the end of life. Our stochastic model used difference equations to model immigration and emigration between zones of the lens epithelium and included some deterministic elements, such as cellular footprint area. We found that the epithelial cell cycle was shortened significantly in the embryo, facilitating the rapid growth that marks early lens development. As development progressed, epithelial cell division becomes non-uniform and four zones, each with a characteristic proliferation rate, could be discerned. Adjustment of two model parameters, proliferation rate and rate of change in cellular footprint area, was sufficient to specify all growth trajectories. Modelling suggested that the direction of cellular migration across zonal boundaries was sensitive to footprint area, a phenomenon that may isolate specific cell populations. Model runs consisted of more than 1000 iterations, in each of which the stochastic behaviour of thousands of cells was followed. Nevertheless, sequential runs were almost superimposable. This remarkable degree of precision was attributed, in part, to the presence of non-mitotic flanking regions, which constituted a path by which epithelial cells could escape the growth process. Spatial modelling suggested that clonal clusters of about 50 cells are produced during migration and that transit times lengthen significantly at later stages, findings with implications for the formation of certain types of cataract.

Development of a synthetic gene network to modulate gene expression by mechanical forces

The majority of (mammalian) cells in our body are sensitive to mechanical forces, but little work has been done to develop assays to monitor mechanosensor activity. Furthermore, it is currently impossible to use mechanosensor activity to drive gene expression. To address these needs, we developed the first mammalian mechanosensitive synthetic gene network to monitor endothelial cell shear stress levels and directly modulate expression of an atheroprotective transcription factor by shear stress. The technique is highly modular, easily scalable and allows graded control of gene expression by mechanical stimuli in hard-to-transfect mammalian cells. We call this new approach mechanosyngenetics. To insert the gene network into a high proportion of cells, a hybrid transfection procedure was developed that involves electroporation, plasmids replication in mammalian cells, mammalian antibiotic selection, a second electroporation and gene network activation. This procedure takes 1 week and yielded over 60% of cells with a functional gene network. To test gene network functionality, we developed a flow setup that exposes cells to linearly increasing shear stress along the length of the flow channel floor. Activation of the gene network varied logarithmically as a function of shear stress magnitude.

Feedback Regulation in a Cancer Stem Cell Model can Cause an Allee Effect

The exact mechanisms of spontaneous tumor remission or complete response to treatment are phenomena in oncology that are not completely understood. We use a concept from ecology, the Allee effect, to help explain tumor extinction in a model of tumor growth that incorporates feedback regulation of stem cell dynamics, which occurs in many tumor types where certain signaling molecules, such as Wnts, are upregulated. Due to feedback and the Allee effect, a tumor may become extinct spontaneously or after therapy even when the entire tumor has not been eradicated by the end of therapy. We quantify the Allee effect using an 'Allee index' that approximates the area of the basin of attraction for tumor extinction. We show that effectiveness of combination therapy in cancer treatment may occur due to the increased probability that the system will be in the Allee region after combination treatment versus monotherapy. We identify therapies that can attenuate stem cell self-renewal, alter the Allee region and increase its size. We also show that decreased response of tumor cells to growth inhibitors can reduce the size of the Allee region and increase stem cell densities, which may help to explain why this phenomenon is a hallmark of cancer.

Quantitative Measurement of Histone Tail Acetylation Reveals Stage-Specific Regulation and Response to Environmental Changes during Drosophila Development

Histone modification plays a major role in regulating gene transcription and ensuring the healthy development of an organism. Numerous studies have suggested that histones are dynamically modified during developmental events to control gene expression levels in a temporal and spatial manner. However, the study of histone acetylation dynamics using currently available techniques is hindered by the difficulty of simultaneously measuring acetylation of the numerous potential sites of modification present in histones. Here, we present a methodology that allows us to combine mass spectrometry-based histone analysis with Drosophila developmental genetics. Using this system, we characterized histone acetylation patterns during multiple developmental stages of the fly. Additionally, we utilized this analysis to characterize how treatments with pharmacological agents or environmental changes such as γ-irradiation altered histone acetylation patterns. Strikingly, γ-irradiation dramatically increased the level of acetylation at H3K18, a site linked to DNA repair via nonhomologous end joining. In mutant fly strains deficient in DNA repair proteins, however, this increase in the level of H3K18 acetylation was lost. These results demonstrate the efficacy of our combined mass spectrometry system with a Drosophila model system and provide interesting insight into the changes in histone acetylation during development, as well as the effects of both pharmacological and environmental agents on global histone acetylation.

Capabilities and Limitations of Tissue Size Control through Passive Mechanical Forces

Embryogenesis is an extraordinarily robust process, exhibiting the ability to control tissue size and repair patterning defects in the face of environmental and genetic perturbations. The size and shape of a developing tissue is a function of the number and size of its constituent cells as well as their geometric packing. How these cellular properties are coordinated at the tissue level to ensure developmental robustness remains a mystery; understanding this process requires studying multiple concurrent processes that make up morphogenesis, including the spatial patterning of cell fates and apoptosis, as well as cell intercalations. In this work, we develop a computational model that aims to understand aspects of the robust pattern repair mechanisms of the Drosophila embryonic epidermal tissues. Size control in this system has previously been shown to rely on the regulation of apoptosis rather than proliferation; however, to date little work has been done to understand the role of cellular mechanics in this process. We employ a vertex model of an embryonic segment to test hypotheses about the emergence of this size control. Comparing the model to previously published data across wild type and genetic perturbations, we show that passive mechanical forces suffice to explain the observed size control in the posterior (P) compartment of a segment. However, observed asymmetries in cell death frequencies across the segment are demonstrated to require patterning of cellular properties in the model. Finally, we show that distinct forms of mechanical regulation in the model may be distinguished by differences in cell shapes in the P compartment, as quantified through experimentally accessible summary statistics, as well as by the tissue recoil after laser ablation experiments.

Patterning of wound-induced intercellular Ca(2+) flashes in a developing epithelium

Differential mechanical force distributions are increasingly recognized to provide important feedback into the control of an organ's final size and shape. As a second messenger that integrates and relays mechanical information to the cell, calcium ions (Ca(2+)) are a prime candidate for providing important information on both the overall mechanical state of the tissue and resulting behavior at the individual-cell level during development. Still, how the spatiotemporal properties of Ca(2+) transients reflect the underlying mechanical characteristics of tissues is still poorly understood. Here we use an established model system of an epithelial tissue, the Drosophila wing imaginal disc, to investigate how tissue properties impact the propagation of Ca(2+) transients induced by laser ablation. The resulting intercellular Ca(2+) flash is found to be mediated by inositol 1,4,5-trisphosphate and depends on gap junction communication. Further, we find that intercellular Ca(2+) transients show spatially non-uniform characteristics across the proximal-distal axis of the larval wing imaginal disc, which exhibit a gradient in cell size and anisotropy. A computational model of Ca(2+) transients is employed to identify the principle factors explaining the spatiotemporal patterning dynamics of intercellular Ca(2+) flashes. The relative Ca(2+) flash anisotropy is principally explained by local cell shape anisotropy. Further, Ca(2+) velocities are relatively uniform throughout the wing disc, irrespective of cell size or anisotropy. This can be explained by the opposing effects of cell diameter and cell elongation on intercellular Ca(2+) propagation. Thus, intercellular Ca(2+) transients follow lines of mechanical tension at velocities that are largely independent of tissue heterogeneity and reflect the mechanical state of the underlying tissue.

Oriented cell division shapes carnivorous pitcher leaves of Sarracenia purpurea

Complex morphology is an evolutionary outcome of phenotypic diversification. In some carnivorous plants, the ancestral planar leaf has been modified to form a pitcher shape. However, how leaf development was altered during evolution remains unknown. Here we show that the pitcher leaves of Sarracenia purpurea develop through cell division patterns of adaxial tissues that are distinct from those in bifacial and peltate leaves, subsequent to standard expression of adaxial and abaxial marker genes. Differences in the orientation of cell divisions in the adaxial domain cause bifacial growth in the distal region and adaxial ridge protrusion in the middle region. These different growth patterns establish pitcher morphology. A computer simulation suggests that the cell division plane is critical for the pitcher morphogenesis. Our results imply that tissue-specific changes in the orientation of cell division underlie the development of a morphologically complex leaf.

The pitcher-shaped leaf of the carnivorous plant Sarracenia purpurea acts as a pitfall trap to capture small animals. Here, Fukushima et al. analyse pitcher leaf development and propose that this unusual shape evolved from ancestral planar leaves through changes in the orientation of cell division.