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.

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.

No abstract