Piezo regulates epithelial topology and promotes precision in organ size control.

Mechanosensitive Piezo channels regulate cell division, cell extrusion, and cell death. However, systems-level functions of Piezo in regulating organogenesis remain poorly understood. Here, we demonstrate that Piezo controls epithelial cell topology to ensure precise organ growth by integrating live-imaging experiments with pharmacological and genetic perturbations and computational modeling. Notably, the knockout or knockdown of Piezo increases bilateral asymmetry in wing size. Piezo's multifaceted functions can be deconstructed as either autonomous or non-autonomous based on a comparison between tissue-compartment-level perturbations or between genetic perturbation populations at the whole-tissue level. A computational model that posits cell proliferation and apoptosis regulation through modulation of the cutoff tension required for Piezo channel activation explains key cell and tissue phenotypes arising from perturbations of Piezo expression levels. Our findings demonstrate that Piezo promotes robustness in regulating epithelial topology and is necessary for precise organ size control.

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.

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.

Vestigial-dependent induction contributes to robust patterning but is not essential for wing-fate recruitment in Drosophila.

Cell recruitment is a process by which a differentiated cell induces neighboring cells to adopt its same cell fate. In Drosophila, cells expressing the protein encoded by the wing selector gene, vestigial (vg), drive a feed-forward recruitment signal that expands the Vg pattern as a wave front. However, previous studies on Vg pattern formation do not reveal these dynamics. Here, we use live imaging to show that multiple cells at the periphery of the wing disc simultaneously activate a fluorescent reporter of the recruitment signal, suggesting that cells may be recruited without the need for their contact neighbors be recruited in advance. In support of this observation, when Vg expression is inhibited either at the dorsal-ventral boundary or away from it, the activation of the recruitment signal still occurs at a distance, suggesting that Vg expression is not absolutely required to send or propagate the recruitment signal. However, the strength and extent of the recruitment signal is clearly compromised. We conclude that a feed-forward, contact-dependent cell recruitment process is not essential for Vg patterning, but it is necessary for robustness. Overall, our findings reveal a previously unidentified role of cell recruitment as a robustness-conferring cell differentiation mechanism.

Rational Design and Identification of Harmine-Inspired, N-Heterocyclic DYRK1A Inhibitors Employing a Functional Genomic In Vivo Drosophila Model System.

Deregulation of dual-specificity tyrosine phosphorylation-regulated kinase 1A (DYRK1A) plays a significant role in developmental brain defects, early-onset neurodegeneration, neuronal cell loss, dementia, and several types of cancer. Herein, we report the discovery of three new classes of N-heterocyclic DYRK1A inhibitors based on the potent, yet toxic kinase inhibitors, harmine and harmol. An initial in vitro evaluation of the small molecule library assembled revealed that the core heterocyclic motifs benzofuranones, oxindoles, and pyrrolones, showed statistically significant DYRK1A inhibition. Further, the utilization of a low cost, high-throughput functional genomic in vivo model system to identify small molecule inhibitors that normalize DYRK1A overexpression phenotypes is described. This in vivo assay substantiated the in vitro results, and the resulting correspondence validates generated classes as architectural motifs that serve as potential DYRK1A inhibitors. Further expansion and analysis of these core compound structures will allow discovery of safe, more effective chemical inhibitors of DYRK1A to ameliorate phenotypes caused by DYRK1A overexpression.

Rab11b-mediated integrin recycling promotes brain metastatic adaptation and outgrowth.

Breast cancer brain metastases (BCBM) have a 5-20 year latency and account for 30% of mortality; however, mechanisms governing adaptation to the brain microenvironment remain poorly defined. We combine time-course RNA-sequencing of BCBM development with a Drosophila melanogaster genetic screen, and identify Rab11b as a functional mediator of metastatic adaptation. Proteomic analysis reveals that Rab11b controls the cell surface proteome, recycling proteins required for successful interaction with the microenvironment, including integrin ß1. Rab11b-mediated control of integrin ß1 surface expression allows efficient engagement with the brain ECM, activating mechanotransduction signaling to promote survival. Lipophilic statins prevent membrane association and activity of Rab11b, and we provide proof-of principle that these drugs prevent breast cancer adaptation to the brain microenvironment. Our results identify Rab11b-mediated recycling of integrin ß1 as regulating BCBM, and suggest that the recycleome, recycling-based control of the cell surface proteome, is a previously unknown driver of metastatic adaptation and outgrowth.

Combined Scaffold Evaluation and Systems-Level Transcriptome-Based Analysis for Accelerated Lead Optimization Reveals Ribosomal Targeting Spirooxindole Cyclopropanes.

With evolutionary drug resistance impacting efforts to treat disease, the need for small molecules that exhibit novel molecular mechanisms of action is paramount. In this study, we combined scaffold-directed synthesis with a hybrid experimental and transcriptome analysis to identify bis-spirooxindole cyclopropanes that inhibit cancer cell proliferation through disruption of ribosomal function. These findings demonstrate the value of an integrated, biologically inspired synthesis and assay strategy for the accelerated identification of first-in-class cancer therapeutic candidates.

Potassium channel activity controls breast cancer metastasis by affecting ß-catenin signaling.

Potassium ion channels are critical in the regulation of cell motility. The acquisition of cell motility is an essential parameter of cancer metastasis. However, the role of K+ channels in cancer metastasis has been poorly studied. High expression of the hG1 gene, which encodes for Kv11.1 channel associates with good prognosis in estrogen receptor-negative breast cancer (BC). We evaluated the efficacy of the Kv11.1 activator NS1643 in arresting metastasis in a triple negative breast cancer (TNBC) mouse model. NS1643 significantly reduces the metastatic spread of breast tumors in vivo by inhibiting cell motility, reprogramming epithelial-mesenchymal transition via attenuation of Wnt/ß-catenin signaling and suppressing cancer cell stemness. Our findings provide important information regarding the clinical relevance of potassium ion channel expression in breast tumors and the mechanisms by which potassium channel activity can modulate tumor biology. Findings suggest that Kv11.1 activators may represent a novel therapeutic approach for the treatment of metastatic estrogen receptor-negative BC. Ion channels are critical factor for cell motility but little is known about their role in metastasis. Stimulation of the Kv11.1 channel suppress the metastatic phenotype in TNBC. This work could represent a paradigm-shifting approach to reducing mortality by targeting a pathway that is central to the development of metastases.

Calcium as a signal integrator in developing epithelial tissues.

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.

Preclinical study of a Kv11.1 potassium channel activator as antineoplastic approach for breast cancer.

Potassium ion (K+) channels have been recently found to play a critical role in cancer biology. Despite that pharmacologic manipulation of ion channels is recognized as an important therapeutic approach, very little is known about the effects of targeting of K+ channels in cancer. In this study, we demonstrate that use of the Kv11.1 K+ channel activator NS1643 inhibits tumor growth in an in vivo model of breast cancer. Tumors exposed to NS1643 had reduced levels of proliferation markers, high expression levels of senescence markers, increased production of ROS and DNA damage compared to tumors of untreated mice. Importantly, mice treated with NS1643 did not exhibit significant cardiac dysfunction. In conclusion, pharmacological stimulation of Kv11.1 activity produced arrested TNBC-derived tumor growth by generating DNA damage and senescence without significant side effects. We propose that use of Kv11.1 channels activators could be considered as a possible pharmacological strategy against breast tumors.

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.

Robust cell tracking in epithelial tissues through identification of maximum common subgraphs.

Tracking of cells in live-imaging microscopy videos of epithelial sheets is a powerful tool for investigating fundamental processes in embryonic development. Characterizing cell growth, proliferation, intercalation and apoptosis in epithelia helps us to understand how morphogenetic processes such as tissue invagination and extension are locally regulated and controlled. Accurate cell tracking requires correctly resolving cells entering or leaving the field of view between frames, cell neighbour exchanges, cell removals and cell divisions. However, current tracking methods for epithelial sheets are not robust to large morphogenetic deformations and require significant manual interventions. Here, we present a novel algorithm for epithelial cell tracking, exploiting the graph-theoretic concept of a 'maximum common subgraph' to track cells between frames of a video. Our algorithm does not require the adjustment of tissue-specific parameters, and scales in sub-quadratic time with tissue size. It does not rely on precise positional information, permitting large cell movements between frames and enabling tracking in datasets acquired at low temporal resolution due to experimental constraints such as phototoxicity. To demonstrate the method, we perform tracking on the Drosophila embryonic epidermis and compare cell-cell rearrangements to previous studies in other tissues. Our implementation is open source and generally applicable to epithelial tissues.

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.

Pattern formation by a moving morphogen source.

During Drosophila melanogaster oogenesis, the follicular epithelium that envelops the germline cyst gives rise to an elaborate eggshell, which houses the future embryo and mediates its interaction with the environment. A prominent feature of the eggshell is a pair of dorsal appendages, which are needed for embryo respiration. Morphogenesis of this structure depends on broad, a zinc-finger transcription factor, regulated by the EGFR pathway. While much has been learned about the mechanisms of broad regulation by EGFR, current understanding of processes that shape the spatial pattern of broad expression is incomplete. We propose that this pattern is defined by two different phases of EGFR activation: an early, posterior-to-anterior gradient of EGFR signaling sets the posterior boundary of broad expression, while the anterior boundary is set by a later phase of EGFR signaling, distributed in a dorsoventral gradient. This model can explain the wild-type pattern of broad in D. melanogaster, predicts how this pattern responds to genetic perturbations, and provides insight into the mechanisms driving diversification of eggshell patterning. The proposed model of the broad expression pattern can be used as a starting point for the quantitative analysis of a large number of gene expression patterns in Drosophila oogenesis.

Feedback control of the EGFR signaling gradient: superposition of domain-splitting events in Drosophila oogenesis.

The morphogenesis of structures with repeated functional units, such as body segments and appendages, depends on multi-domain patterns of cell signaling and gene expression. We demonstrate that during Drosophila oogenesis, the two-domain expression pattern of Broad, a transcription factor essential for the formation of the two respiratory eggshell appendages, is established by a single gradient of EGFR activation that induces both Broad and Pointed, which mediates repression of Broad. Two negative-feedback loops provided by the intracellular inhibitors of EGFR signaling, Kekkon-1 and Sprouty, control the number and position of Broad-expressing cells and in this way influence eggshell morphology. Later in oogenesis, the gradient of EGFR activation is split into two smaller domains in a process that depends on Argos, a secreted antagonist of EGFR signaling. In contrast to the previously proposed model of eggshell patterning, we show that the two-domain pattern of EGFR signaling is not essential for specifying the number of appendages. Thus, the processes that define the two-domain patterns of Broad and EGFR activation are distinct; their actions are separated in time and have different effects on eggshell morphology.

Bistability coordinates activation of the EGFR and DPP pathways in Drosophila vein differentiation.

Cell differentiation in developing tissues is controlled by a small set of signaling pathways, which must coordinate the timing and levels of activation to ensure robust and precise outcomes. Highly coordinated activation of signaling pathways can result from cross-regulatory interactions in multi-pathway networks. Here we explore the dynamics and function of pathway coordination between the EGFR and DPP pathways during Drosophila wing-vein differentiation. We show that simultaneous activation of both the EGFR and DPP pathways must be maintained for vein cell differentiation and that above-threshold ectopic activation of either pathway is sufficient to drive vein cell differentiation outside the proveins. The joint activation of the EGFR and DPP signaling systems is ensured by a positive feedback loop, in which the two pathways stimulate each other at the level of ligand production.

Cad74A is regulated by BR and is required for robust dorsal appendage formation in Drosophila oogenesis.

Drosophila egg development is an established model for studying epithelial patterning and morphogenesis, but the connection between signaling pathways and egg morphology is still incompletely understood. We have identified a non-classical cadherin, Cad74A, as a putative adhesion gene that bridges epithelial patterning and morphogenesis in the follicle cells. Starting in mid-oogenesis, Cad74A is expressed in the follicle cells that contact the oocyte, including the border cells and most of the columnar follicle cells. However, Cad74A is repressed in two dorsolateral patches of follicle cells, which participate in the formation of tubular respiratory appendages. We show genetically that Cad74A is downstream of the EGFR and BMP signaling pathways and is repressed by the Zn-finger transcription factor Broad. The correlation of Cad74A repression in the cells that bend out of the plane of the follicular epithelium is preserved across Drosophila species and mutant backgrounds exhibiting a range of eggshell phenotypes. Complete removal of Cad74A from the follicle cells causes defects in dorsal appendage formation. Ectopic expression of Cad74A in the roof cells results in shortened, flattened appendages due to the hindered migration of the roof cells. Based on these results, we propose that Cad74A is part of the adhesive machinery that enables robust dorsal appendage formation, and as such provides a link between the patterning of the follicle cells and eggshell morphogenesis.