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 Gaq-mediated signaling. In vivo calcium signaling dynamics depend on both receptor tyrosine kinase/phospholipase C g and G protein-coupled receptor/phospholipase C b 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.

Decoding how tissue properties emerge across multiple spatial and temporal scales from the integration of local signals is a grand challenge in quantitative biology. For example, the collective behavior of epithelial cells is critical for shaping developing embryos. Understanding how epithelial cells interpret a diverse range of local signals to coordinate tissue-level processes requires a systems-level understanding of development. Integration of multiple signaling pathways that specify cell signaling information requires second messengers such as calcium ions. Increasingly, specific roles have been uncovered for calcium signaling throughout development. Calcium signaling regulates many processes including division, migration, death, and differentiation. However, the pleiotropic and ubiquitous nature of calcium signaling implies that many additional functions remain to be discovered. Here we review a selection of recent studies to highlight important insights into how multiple signals are transduced by calcium transients in developing epithelial tissues. Quantitative imaging and computational modeling have provided important insights into how calcium signaling integration occurs. Reverse-engineering the conserved features of signal integration mediated by calcium signaling will enable novel approaches in regenerative medicine and synthetic control of morphogenesis.

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 (Ca2+) 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 Ca2+ 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 Ca2+ transients induced by laser ablation. The resulting intercellular Ca2+ flash is found to be mediated by inositol 1,4,5-trisphosphate (IP3) and depends on gap junction communication. Further, we find that intercellular Ca2+ transients show spatially nonuniform characteristics across the proximal-distal (PD) axis of the larval wing imaginal disc, which exhibit a gradient in cell size and anisotropy. A computational model of Ca2+ transients is employed to identify the principle factors explaining the spatiotemporal patterning dynamics of intercellular Ca2+ flashes. The relative Ca2+ flash anisotropy is principally explained by local cell shape anisotropy. Further, Ca2+ 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 Ca2+ propagation. Thus, intercellular Ca2+ transients follow lines of mechanical tension at velocities that are largely independent of tissue heterogeneity and reflect the mechanical state of the underlying tissue.

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.

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.