Encoding of temporal signals by the TGF-β pathway and implications for embryonic patterning

Frequency-Dependent Premature Differentiation of Pheochromocytoma Cells Exhibits Band-Pass Filter Behavior Correlation with Intracellular Enzyme Activation Kinetics

Advances in microfluidics, optogenetics and electronics have enabled the study of dynamically controlled inputs on cellular fate. Here, we applied a microfluidic system to deliver periodic inputs of growth factors to pheochromocytoma cells and measured the extent of premature differentiation as a function of input frequency. Epidermal growth factor-triggered differentiation peaked at two cycles/hour, while nerve growth factor-triggered differentiation peaked at one cycle/hour. To interpret the results, we analyzed a published model that attributed pheochromocytoma cell differentiation to the linear combination of activated enzymes extracellular signal-regulated kinase (ERK), cAMP response element binding protein (CREB), protein kinase B (AKT) and c-Jun N-terminal kinase (JNK) at specific times after step input stimulation. Transfer functions for enzyme activation were derived from the published time-domain activation kinetics and these transfer functions were combined in a parallel architecture as a predictor of neurite outgrowth, as a function of input frequency. Qualitative agreement was observed between the model and the experiments.

Temporal dynamics of BMP/Nodal ratio drive tissue-specific gastrulation morphogenesis

Anteroposterior elongation of the vertebrate body plan is driven by convergence and extension (C&E) gastrulation movements in both the mesoderm and neuroectoderm, but how or whether molecular regulation of C&E differs between tissues remains an open question. Using a zebrafish explant model of anteroposterior axis extension, we show that C&E of the neuroectoderm and mesoderm can be uncoupled ex vivo, and that morphogenesis of individual tissues results from distinct morphogen signaling dynamics. Using precise temporal manipulation of BMP and Nodal signaling, we identify a critical developmental window during which high or low BMP/Nodal ratios induce neuroectoderm- or mesoderm-driven C&E, respectively. Increased BMP activity similarly enhances C&E specifically in the ectoderm of intact zebrafish gastrulae, highlighting the in vivo relevance of our findings. Together, these results demonstrate that temporal dynamics of BMP and Nodal morphogen signaling activate distinct morphogenetic programs governing C&E gastrulation movements within individual tissues.

Optogenetic control of Nodal signaling patterns

A crucial step in early embryogenesis is the establishment of spatial patterns of signaling activity. Tools to perturb morphogen signals with high resolution in space and time can help reveal how embryonic cells decode these signals to make appropriate fate decisions. Here, we present new optogenetic reagents and an experimental pipeline for creating designer Nodal signaling patterns in live zebrafish embryos. Nodal receptors were fused to the light-sensitive heterodimerizing pair Cry2/CIB1N, and the type II receptor was sequestered to the cytosol. The improved optoNodal2 reagents eliminate dark activity and improve response kinetics, without sacrificing dynamic range. We adapted an ultra-widefield microscopy platform for parallel light patterning in up to 36 embryos, and demonstrated precise spatial control over Nodal signaling activity and downstream gene expression. Patterned Nodal activation drove precisely controlled internalization of endodermal precursors. Furthermore, we used patterned illumination to generate synthetic signaling patterns in Nodal signaling mutants, rescuing several characteristic developmental defects. This study establishes an experimental toolkit for systematic exploration of Nodal signaling patterns in live embryos.

Bioinspired Nucleic Acid-Based Bandpass Filters and Their Concentration-Adaptive Functions

Natural signaling networks can act as bandpass filters to interpret external stimuli within defined concentration ranges for differential cellular activities. Replicating such a bandpass filtering mechanism by synthetic networks poses a significant challenge. Herein, we introduce a modular design of nucleic acid-based multilayer threshold-gated incoherent feedforward networks as multiband bandpass filters to produce mutually exclusive responses within defined input concentration ranges. In these networks, nucleic acids demonstrate triple functionality by acting as threshold-gated entities to discern input concentration levels, serving as network nodes to assemble incoherent feedforward loops for nonlinear signal processing, and functioning as signal transduction units for coupling downstream functional modules. These modular networks enable the fine-tuning of filtering performance in terms of band position, bandwidth, cascades, and responses. A mathematical simulation model allows us to predict the filtering behaviors under various conditions. Also, the networks are integrated with upstream signal conversion modules to process concentration information on molecules beyond nucleic acids, such as adenosine and its derivatives. Furthermore, connections to downstream functional modules allow the system to regulate various processes in a concentration bandpass manner, realizing concentration-adaptive DNAzyme biocatalysis, tristate logic operations, RNA transcription, and DNA condensate formation. These findings underscore the potential of enzyme-free DNA reaction networks in complex signal processing and lay a solid foundation for developing chemical and material systems with highly adaptive and autonomous behaviors.

Phenotypic consequences of logarithmic signaling in MAPK stress response

How cells respond to dynamic environmental changes is crucial for understanding fundamental biological processes and cell physiology. In this study, we developed an experimental and quantitative analytical framework to explore how dynamic stress gradients that change over time regulate cellular volume, signaling activation, and growth phenotypes. Our findings reveal that gradual stress conditions substantially enhance cell growth compared to conventional acute stress. This growth advantage correlates with a minimal reduction in cell volume dependent on the dynamic of stress. We explain the growth phenotype with our finding of a logarithmic signal transduction mechanism in the yeast mitogen-activated protein kinase (MAPK) osmotic stress response pathway. These insights into the interplay between gradual environments, cell volume change, dynamic cell signaling, and growth, advance our understanding of fundamental cellular processes in gradual stress environments.

Transcriptional regulators ensuring specific gene expression and decision-making at high TGFβ doses

TGFβ-signaling regulates cancer progression by controlling cell division, migration, and death. These outcomes are mediated by gene expression changes, but the mechanisms of decision-making toward specific fates remain unclear. Here, we combine SMAD transcription factor imaging, genome-wide RNA sequencing, and morphological assays to quantitatively link signaling, gene expression, and fate decisions in mammary epithelial cells. Fitting genome-wide kinetic models to our time-resolved data, we find that most of the TGFβ target genes can be explained as direct targets of SMAD transcription factors, whereas the remainder show signs of complex regulation, involving delayed regulation and strong amplification at high TGFβ doses. Knockdown experiments followed by global RNA sequencing revealed transcription factors interacting with SMADs in feedforward loops to control delayed and dose-discriminating target genes, thereby reinforcing the specific epithelial-to-mesenchymal transition at high TGFβ doses. We identified early repressors, preventing premature activation, and a late activator, boosting gene expression responses for a sufficiently strong TGFβ stimulus. Taken together, we present a global view of TGFβ-dependent gene regulation and describe specificity mechanisms reinforcing cellular decision-making.

Depolarization induces calcium-dependent BMP4 release from mouse embryonic palate mesenchymal cells

Bone Morphogenetic Protein (BMP) signaling is essential for craniofacial development, though little is known about the mechanisms that govern BMP secretion. We show that depolarization induces calcium-dependent BMP4 release from mouse embryonic palate mesenchyme. We show endogenous transient changes in intracellular calcium occur in cranial neural crest cells, the cells from which embryonic palate mesenchyme derives. Waves of transient changes in intracellular calcium suggest that these cells are electrically coupled and may temporally coordinate BMP release. These transient changes in intracellular calcium persist in palate mesenchyme cells from embryonic day 9.5 to 13.5 mice. Disruption of a potassium channel called Kcnj2 significantly decreases the amplitude of calcium transients and the ability of cells to secrete BMP. Kcnj2 knockout mice have cleft palate and reduced BMP signaling. Our data suggest that temporal control of developmental cues is regulated by ion channels, depolarization, and intracellular calcium for mammalian craniofacial morphogenesis.

Depolarization induces calcium-dependent BMP4 release from mouse embryonic palate mesenchyme

Ion channels are essential for proper morphogenesis of the craniofacial skeleton. However, the molecular mechanisms underlying this phenomenon are unknown. Loss of the Kcnj2 potassium channel disrupts Bone Morphogenetic Protein (BMP) signaling within the developing palate. BMP signaling is essential for the correct development of several skeletal structures, including the palate, though little is known about the mechanisms that govern BMP secretion. We introduce a tool to image the release of bone morphogenetic protein 4 (BMP4) from mammalian cells. Using this tool, we show that depolarization induces BMP4 release from mouse embryonic palate mesenchyme cells in a calcium-dependent manner. We show native transient changes in intracellular calcium occur in cranial neural crest cells, the cells from which embryonic palate mesenchyme derives. Waves of transient changes in intracellular calcium suggest that these cells are electrically coupled and may temporally coordinate BMP release. These transient changes in intracellular calcium persist in palate mesenchyme cells from embryonic day (E) 9.5 to 13.5 mice. Disruption of Kcnj2 significantly decreases the amplitude of calcium transients and the ability of cells to secrete BMP. Together, these data suggest that temporal control of developmental cues is regulated by ion channels, depolarization, and changes in intracellular calcium for mammalian craniofacial morphogenesis.

SUMMARY

We show that embryonic palate mesenchyme cells undergo transient changes in intracellular calcium. Depolarization of these cells induces BMP4 release suggesting that ion channels are a node in BMP4 signaling.

Liebig's law of the minimum in the TGF-β/SMAD pathway

Cells use signaling pathways to sense and respond to their environments. The transforming growth factor-β (TGF-β) pathway produces context-specific responses. Here, we combined modeling and experimental analysis to study the dependence of the output of the TGF-β pathway on the abundance of signaling molecules in the pathway. We showed that the TGF-β pathway processes the variation of TGF-β receptor abundance using Liebig's law of the minimum, meaning that the output-modifying factor is the signaling protein that is most limited, to determine signaling responses across cell types and in single cells. We found that the abundance of either the type I (TGFBR1) or type II (TGFBR2) TGF-β receptor determined the responses of cancer cell lines, such that the receptor with relatively low abundance dictates the response. Furthermore, nuclear SMAD2 signaling correlated with the abundance of TGF-β receptor in single cells depending on the relative expression levels of TGFBR1 and TGFBR2. A similar control principle could govern the heterogeneity of signaling responses in other signaling pathways.

Optogenetic control of Nodal signaling patterns

A crucial step in early embryogenesis is the establishment of spatial patterns of signaling activity. Tools to perturb morphogen signals with high resolution in space and time can help reveal how embryonic cells decode these signals to make appropriate fate decisions. Here, we present new optogenetic reagents and an experimental pipeline for creaHng designer Nodal signaling patterns in live zebrafish embryos. Nodal receptors were fused to the light-sensitive heterodimerizing pair Cry2/CIB1N, and the Type II receptor was sequestered to the cytosol. The improved optoNodal2 reagents eliminate dark activity and improve response kinetics, without sacrificing dynamic range. We adapted an ultra-widefield microscopy platform for parallel light patterning in up to 36 embryos and demonstrated precise spatial control over Nodal signaling activity and downstream gene expression. Patterned Nodal activation drove precisely controlled internalization of endodermal precursors. Further, we used patterned illumination to generate synthetic signaling patterns in Nodal signaling mutants, rescuing several characteristic developmental defects. This study establishes an experimental toolkit for systematic exploration of Nodal signaling patterns in live embryos.

Temporal dynamics of BMP/Nodal ratio drive tissue-specific gastrulation morphogenesis

Anteroposterior (AP) elongation of the vertebrate body plan is driven by convergence and extension (C&E) gastrulation movements in both the mesoderm and neuroectoderm, but how or whether molecular regulation of C&E differs between tissues remains an open question. Using a zebrafish explant model of AP axis extension, we show that C&E of the neuroectoderm and mesoderm can be uncoupled ex vivo, and that morphogenesis of individual tissues results from distinct morphogen signaling dynamics. Using precise temporal manipulation of BMP and Nodal signaling, we identify a critical developmental window during which high or low BMP/Nodal ratios induce neuroectoderm- or mesoderm-driven C&E, respectively. Increased BMP activity similarly enhances C&E specifically in the ectoderm of intact zebrafish gastrulae, highlighting the in vivo relevance of our findings. Together, these results demonstrate that temporal dynamics of BMP and Nodal morphogen signaling activate distinct morphogenetic programs governing C&E gastrulation movements within individual tissues.

Phenotypic consequences of logarithmic signaling in MAPK stress response

How cells respond to dynamic environmental changes is crucial for understanding fundamental biological processes and cell physiology. In this study, we developed an experimental and quantitative analytical framework to explore how dynamic stress gradients that change over time regulate cellular volume, signaling activation, and growth phenotypes. Our findings reveal that gradual stress conditions substantially enhance cell growth compared to conventional acute stress. This growth advantage correlates with a minimal reduction in cell volume dependent on the dynamic of stress. We explain the growth phenotype with our finding of a logarithmic signal transduction mechanism in the yeast Mitogen-Activated Protein Kinase (MAPK) osmotic stress response pathway. These insights into the interplay between gradual environments, cell volume change, dynamic cell signaling, and growth, advance our understanding of fundamental cellular processes in gradual stress environments.

Endoplasmic Reticulum Calcium Mediates Drosophila Wing Development

Background

The temporal dynamics of morphogen presentation impacts transcriptional responses and tissue patterning. However, the mechanisms controlling morphogen release are far from clear. We found that inwardly rectifying potassium (Irk) channels regulate endogenous transient increases in intracellular calcium and bone morphogenetic protein (BMP/Dpp) release for Drosophila wing development. Inhibition of Irk channels reduces BMP/Dpp signaling, and ultimately disrupts wing morphology. Ion channels impact development of several tissues and organisms in which BMP signaling is essential. In neurons and pancreatic beta cells, Irk channels modulate membrane potential to affect intracellular Ca++ to control secretion of neurotransmitters and insulin. Based on Irk activity in neurons, we hypothesized that electrical activity controls endoplasmic reticulum (ER) Ca++ release into the cytoplasm to regulate the release of BMP.

Materials and Methods

To test this hypothesis, we reduced expression of four proteins that control ER calcium, Stromal interaction molecule 1 (Stim), Calcium release-activated calcium channel protein 1 (Orai), SarcoEndoplasmic Reticulum Calcium ATPase (SERCA), small conductance calcium-activated potassium channel (SK), and Bestrophin 2 (Best2) using RNAi and documented wing phenotypes. We use live imaging to study calcium and Dpp release within pupal wings and larval wing discs. Additionally, we employed immunohistochemistry to characterize Small Mothers Against Decapentaplegic (SMAD) phosphorylation downstream of the BMP/Dpp pathway following RNAi knockdown.

Results

We found that reduced Stim and SERCA function decreases amplitude and frequency of endogenous calcium transients in the wing disc and reduced BMP/Dpp release.

Conclusion

Our results suggest control of ER calcium homeostasis is required for BMP/Dpp release, and Drosophila wing development.

State- and stimulus-specific dynamics of SMAD signaling determine fate decisions in individual cells

SMAD-mediated signaling regulates apoptosis, cell cycle arrest, and epithelial-to-mesenchymal transition to safeguard tissue homeostasis. However, it remains elusive how the relatively simple pathway can determine such a broad range of cell fate decisions and how it differentiates between varying ligands. Here, we systematically investigate how SMAD-mediated responses are modulated by various ligands of the transforming growth factor β (TGFβ) family and compare these ligand responses in quiescent and proliferating MCF10A cells. We find that the nature of the phenotypic response is mainly determined by the proliferation status, with migration and cell cycle arrest being dominant in proliferating cells for all tested TGFβ family ligands, whereas cell death is the major outcome in quiescent cells. In both quiescent and proliferating cells, the identity of the ligand modulates the strength of the phenotypic response proportional to the dynamics of induced SMAD nuclear-to-cytoplasmic translocation and, as a consequence, the corresponding gene expression changes. Interestingly, the proliferation state of a cell has little impact on the set of genes induced by SMAD signaling; instead, it modulates the relative cellular sensitivity to TGFβ superfamily members. Taken together, diversity of SMAD-mediated responses is mediated by differing cellular states, which determine ligand sensitivity and phenotypic effects, while the pathway itself merely serves as a quantitative relay from the cell membrane to the nucleus.

Modeling Cellular Signaling Variability Based on Single-Cell Data: The TGFβ-SMAD Signaling Pathway

Nongenetic heterogeneity is key to cellular decisions, as even genetically identical cells respond in very different ways to the same external stimulus, e.g., during cell differentiation or therapeutic treatment of disease. Strong heterogeneity is typically already observed at the level of signaling pathways that are the first sensors of external inputs and transmit information to the nucleus where decisions are made. Since heterogeneity arises from random fluctuations of cellular components, mathematical models are required to fully describe the phenomenon and to understand the dynamics of heterogeneous cell populations. Here, we review the experimental and theoretical literature on cellular signaling heterogeneity, with special focus on the TGFβ/SMAD signaling pathway.

Rate thresholds in cell signaling have functional and phenotypic consequences in non-linear time-dependent environments

All cells employ signal transduction pathways to respond to physiologically relevant extracellular cytokines, stressors, nutrient levels, hormones, morphogens, and other stimuli that vary in concentration and rate in healthy and diseased states. A central unsolved fundamental question in cell signaling is whether and how cells sense and integrate information conveyed by changes in the rate of extracellular stimuli concentrations, in addition to the absolute difference in concentration. We propose that different environmental changes over time influence cell behavior in addition to different signaling molecules or different genetic backgrounds. However, most current biomedical research focuses on acute environmental changes and does not consider how cells respond to environments that change slowly over time. As an example of such environmental change, we review cell sensitivity to environmental rate changes, including the novel mechanism of rate threshold. A rate threshold is defined as a threshold in the rate of change in the environment in which a rate value below the threshold does not activate signaling and a rate value above the threshold leads to signal activation. We reviewed p38/Hog1 osmotic stress signaling in yeast, chemotaxis and stress response in bacteria, cyclic adenosine monophosphate signaling in Amoebae, growth factors signaling in mammalian cells, morphogen dynamics during development, temporal dynamics of glucose and insulin signaling, and spatio-temproral stressors in the kidney. These reviewed examples from the literature indicate that rate thresholds are widespread and an underappreciated fundamental property of cell signaling. Finally, by studying cells in non-linear environments, we outline future directions to understand cell physiology better in normal and pathophysiological conditions.

TGF-β signaling and Creb5 cooperatively regulate Fgf18 to control pharyngeal muscle development

The communication between myogenic cells and their surrounding connective tissues is indispensable for muscle morphogenesis. During late embryonic development in mice, myogenic progenitors migrate to discrete sites to form individual muscles. The detailed mechanism of this process remains unclear. Using mouse levator veli palatini (LVP) development as a model, we systematically investigated how a distinct connective tissue subpopulation, perimysial fibroblasts, communicates with myogenic cells to regulate mouse pharyngeal myogenesis. Using single-cell RNAseq data analysis, we identified that TGF-β signaling is a key regulator for the perimysial fibroblasts. Loss of TGF-β signaling in the neural crest-derived palatal mesenchyme leads to defects in perimysial fibroblasts and muscle malformation in the soft palate in Osr2Cre;Tgfbr1fl/fl mice. In particular, Creb5, a transcription factor expressed in the perimysial fibroblasts, cooperates with TGF-β signaling to activate expression of Fgf18. Moreover, Fgf18 supports pharyngeal muscle development in vivo and exogenous Fgf18 can partially rescue myogenic cell numbers in Osr2Cre;Tgfbr1fl/fl samples, illustrating that TGF-β-regulated Fgf18 signaling is required for LVP development. Collectively, our findings reveal the mechanism by which TGF-β signaling achieves its functional specificity in defining the perimysial-to-myogenic signals for pharyngeal myogenesis.

Roles of TGF-β in cancer hallmarks and emerging onco-therapeutic design

Transforming growth factor-beta (TGF-β) is a double-edged sword in cancer treatment because of its pivotal yet complex and roles played during cancer initiation/development. Current anti-cancer strategies involving TGF-β largely view TGF-β as an onco-therapeutic target that not only substantially hinders its full utilisation for cancer control, but also considerably restricts innovations in this field. Thereby, how to take advantages of therapeutically favourable properties of TGF-β for cancer management represents an interesting and less investigated problem. Here, by categorising cancer hallmarks into four critical transition events and one enabling characteristic controlling cancer initiation and progression, and delineating TGF-β complexities according to these cancer traits, we identify the suppressive role of TGF-β in tumour initiation and early-stage progression and its promotive functionalities in cancer metastasis as well as other cancer hallmarks. We also propose the feasibility and possible scenarios of combining cold atmospheric plasma (CAP) with onco-therapeutics utilising TGF-β for cancer control given the intrinsic properties of CAP against cancer hallmarks.

Overcoming therapeutic resistance to platinum-based drugs by targeting Epithelial–Mesenchymal transition

Platinum-based drugs (PBDs), including cisplatin, carboplatin, and oxaliplatin, have been widely used in clinical practice as mainstay treatments for various types of cancer. Although there is firm evidence of notable achievements with PBDs in the management of cancers, the acquisition of resistance to these agents is still a major challenge to efforts at cure. The introduction of the epithelial-mesenchymal transition (EMT) concept, a critical process during embryonic morphogenesis and carcinoma progression, has offered a mechanistic explanation for the phenotypic switch of cancer cells upon PBD exposure. Accumulating evidence has suggested that carcinoma cells can enter a resistant state via induction of the EMT. In this review, we discussed the underlying mechanism of PBD-induced EMT and the current understanding of its role in cancer drug resistance, with emphasis on how this novel knowledge can be exploited to overcome PBD resistance via EMT-targeted compounds, especially those under clinical trials.

An automated do-it-yourself system for dynamic stem cell and organoid culture in standard multi-well plates

We present a low-cost, do-it-yourself system for complex mammalian cell culture under dynamically changing medium formulations by integrating conventional multi-well tissue culture plates with simple microfluidic control and system automation. We demonstrate the generation of complex concentration profiles, enabling the investigation of sophisticated input-response relations. We further apply our automated cell-culturing platform to the dynamic stimulation of two widely employed stem-cell-based in vitro models for early mammalian development: the conversion of naive mouse embryonic stem cells into epiblast-like cells and mouse 3D gastruloids. Performing automated medium-switch experiments, we systematically investigate cell fate commitment along the developmental trajectory toward mouse epiblast fate and examine symmetry-breaking, germ layer formation, and cardiac differentiation in mouse 3D gastruloids as a function of time-varying Wnt pathway activation. With these proof-of-principle examples, we demonstrate a highly versatile and scalable tool that can be adapted to specific research questions, experimental demands, and model systems.

Data-based stochastic modeling reveals sources of activity bursts in single-cell TGF-β signaling

Cells sense their surrounding by employing intracellular signaling pathways that transmit hormonal signals from the cell membrane to the nucleus. TGF-β/SMAD signaling encodes various cell fates, controls tissue homeostasis and is deregulated in diseases such as cancer. The pathway shows strong heterogeneity at the single-cell level, but quantitative insights into mechanisms underlying fluctuations at various time scales are still missing, partly due to inefficiency in the calibration of stochastic models that mechanistically describe signaling processes. In this work we analyze single-cell TGF-β/SMAD signaling and show that it exhibits temporal stochastic bursts which are dose-dependent and whose number and magnitude correlate with cell migration. We propose a stochastic modeling approach to mechanistically describe these pathway fluctuations with high computational efficiency. Employing high-order numerical integration and fitting to burst statistics we enable efficient quantitative parameter estimation and discriminate models that assume noise in different reactions at the receptor level. This modeling approach suggests that stochasticity in the internalization of TGF-β receptors into endosomes plays a key role in the observed temporal bursting. Further, the model predicts the single-cell dynamics of TGF-β/SMAD signaling in untested conditions, e.g., successfully reflects memory effects of signaling noise and cellular sensitivity towards repeated stimulation. Taken together, our computational framework based on burst analysis, noise modeling and path computation scheme is a suitable tool for the data-based modeling of complex signaling pathways, capable of identifying the source of temporal noise.

Dynamic Visualization of TGF-β/SMAD3 Transcriptional Responses in Single Living Cells

Simple Summary

How a single cytokine can induce a variety of cellular responses in the same cell or in different cells is a longstanding question. Transforming growth factor β (TGF-β) is a prototypical multifunctional cytokine of which biological responses are highly dependent on in a cellular context. TGF-β signals via intracellular SMAD transcription factors, and the duration and intensity of SMAD activation are key determinants for the responses that are elicited by TGF-β. To visualize the TGF-β signaling kinetics, we developed a dynamic TGF-β/SMAD3 transcriptional reporter using a quickly folded and highly unstable green florescent protein. We demonstrate the specificity and sensitivity of this reporter and its wide application to monitor dynamic TGF-β-induced responses in cells cultured on plastic dishes, and in living animals. This tool allows for the analysis of TGF-β signaling at a single living cell level, and allows for the discovery of dynamic TGF-β SMAD- induced transcriptional responses in multi-step biological processes.

Abstract

Transforming growth factor-β (TGF-β) signaling is tightly controlled in duration and intensity during embryonic development and in the adult to maintain tissue homeostasis. To visualize the TGF-β/SMAD3 signaling kinetics, we developed a dynamic TGF-β/SMAD3 transcriptional fluorescent reporter using multimerized SMAD3/4 binding elements driving the expression of a quickly folded and highly unstable GFP protein. We demonstrate the specificity and sensitivity of this reporter and its wide application to monitor dynamic TGF-β/SMAD3 transcriptional responses in both 2D and 3D systems in vitro, as well as in vivo, using live-cell and intravital imaging. Using this reporter in B16F10 cells, we observed single cell heterogeneity in response to TGF-β challenge, which can be categorized into early, late, and non-responders. Because of its broad application potential, this reporter allows for new discoveries into how TGF-β/SMAD3-dependent transcriptional dynamics are affected during multistep and reversible biological processes.

Using Microfluidics and Live Cell Reporters to Dissect the Dynamics of TGF-β Signaling in Mouse Embryonic Stem Cells

The TGF-β pathway is known to behave as a classical morphogen, meaning that it can dictate cell fate decisions in a dose-dependent manner. Recent observations however showed that in addition to the absolute value of morphogen concentration, cells could also extract information from its temporal variations. In the present article we describe how to use automated microfluidics cell culture to stimulate cells with precisely defined temporal profiles of morphogens and how to engineer mouse embryonic stem cells with fluorescent reporters of pathway activity to record in real time their response to the applied stimulations. The combination of automated cell culture and of live cell reporter provides a complete toolbox to study how cells encode the information carried by time-varying TGF-β signals.

Geometry of gene regulatory dynamics

Embryonic development leads to the reproducible and ordered appearance of complexity from egg to adult. The successive differentiation of different cell types that elaborate this complexity results from the activity of gene networks and was likened by Waddington to a flow through a landscape in which valleys represent alternative fates. Geometric methods allow the formal representation of such landscapes and codify the types of behaviors that result from systems of differential equations. Results from Smale and coworkers imply that systems encompassing gene network models can be represented as potential gradients with a Riemann metric, justifying the Waddington metaphor. Here, we extend this representation to include parameter dependence and enumerate all three-way cellular decisions realizable by tuning at most two parameters, which can be generalized to include spatial coordinates in a tissue. All diagrams of cell states vs. model parameters are thereby enumerated. We unify a number of standard models for spatial pattern formation by expressing them in potential form (i.e., as topographic elevation). Turing systems appear nonpotential, yet in suitable variables the dynamics are low dimensional and potential. A time-independent embedding recovers the original variables. Lateral inhibition is described by a saddle point with many unstable directions. A model for the patterning of the Drosophila eye appears as relaxation in a bistable potential. Geometric reasoning provides intuitive dynamic models for development that are well adapted to fit time-lapse data.

Desensitisation of Notch signalling through dynamic adaptation in the nucleus

During embryonic development, signalling pathways orchestrate organogenesis by controlling tissue‐specific gene expression programmes and differentiation. Although the molecular components of many common developmental signalling systems are known, our current understanding of how signalling inputs are translated into gene expression outputs in real‐time is limited. Here we employ optogenetics to control the activation of Notch signalling during Drosophila embryogenesis with minute accuracy and follow target gene expression by quantitative live imaging. Light‐induced nuclear translocation of the Notch Intracellular Domain (NICD) causes a rapid activation of target mRNA expression. However, target gene transcription gradually decays over time despite continuous photo‐activation and nuclear NICD accumulation, indicating dynamic adaptation to the signalling input. Using mathematical modelling and molecular perturbations, we show that this adaptive transcriptional response fits to known motifs capable of generating near‐perfect adaptation and can be best explained by state‐dependent inactivation at the target cis‐regulatory region. Taken together, our results reveal dynamic nuclear adaptation as a novel mechanism controlling Notch signalling output during tissue differentiation.

Growth across scales: Dynamic signaling impacts tissue size and shape

Organ and tissue growth result from an integration of biophysical communication across biological scales, both in time and space. In this review, we highlight new insight into the dynamic connections between control mechanisms operating at different length scales. First, we consider how the dynamics of chemical and electrical signaling in the shape of gradients or waves affect spatiotemporal signal interpretation. Then, we discuss the mechanics underlying dynamic cell behavior during oriented tissue growth, followed by the connections between signaling at the tissue and organismal levels.

Differential compartmentalization of BMP4/NOGGIN requires NOGGIN trans-epithelial transport

Using self-organizing human models of gastrulation, we previously showed that (1) BMP4 initiates the cascade of events leading to gastrulation, (2) BMP4 signal reception is restricted to the basolateral domain, and (3) in a human-specific manner, BMP4 directly induces the expression of NOGGIN. Here, we report the surprising discovery that in human epiblasts, NOGGIN and BMP4 were secreted into opposite extracellular spaces. Interestingly, apically presented NOGGIN could inhibit basally delivered BMP4. Apically imposed microfluidic flow demonstrated that NOGGIN traveled in the apical extracellular space. Our co-localization analysis detailed the endocytotic route that trafficked NOGGIN from the apical space to the basolateral intercellular space where BMP4 receptors were located. This apical-basal transcytosis was indispensable for NOGGIN inhibition. Taken together, the segregation of activator/inhibitor into distinct extracellular spaces challenges classical views of morphogen movement. We propose that the transport of morphogen inhibitors regulates the spatial availability of morphogens during embryogenesis.

Self-organized signaling in stem cell models of embryos

Mammalian embryonic development is a complex process driven by self-organization. Understanding how a fertilized egg develops into an embryo composed of more than 200 cell types in precise spatial patterns remains one of the fundamental challenges in biology. Pluripotent stem cells have been used as in vitro models for investigating mammalian development, and represent promising building blocks for regenerative therapies. Recently, sophisticated stem cell-based models that recapitulate early embryonic fate patterning and morphogenesis have been developed. In this article, we review recent advances in stem cell models of embryos in particular focusing on signaling activities underpinning cell fate decisions in space and time.

Ubiquitin-dependent regulation of transcription in development and disease

Transcription is an elaborate process that is required to establish and maintain the identity of the more than two hundred cell types of a metazoan organism. Strict regulation of gene expression is therefore vital for tissue formation and homeostasis. An accumulating body of work found that ubiquitylation of histones, transcription factors, or RNA polymerase II is crucial for ensuring that transcription occurs at the right time and place during development. Here, we will review principles of ubiquitin-dependent control of gene expression and discuss how breakdown of these regulatory circuits leads to a wide array of human diseases.

Kinetics of osmotic stress regulate a cell fate switch of cell survival

Exposure of cells to diverse types of stressful environments differentially regulates cell fate. Although many types of stresses causing this differential regulation are known, it is unknown how changes over time of the same stressor regulate cell fate. Changes in extracellular osmolarity are critically involved in physiological and pathophysiological processes in several tissues. We observe that human cells survive gradual but not acute hyperosmotic stress. We find that stress, caspase, and apoptosis signaling do not activate during gradual stress in contrast to acute treatments. Contrary to the current paradigm, we see a substantial accumulation of proline in cells treated with gradual but not acute stresses. We show that proline can protect cells from hyperosmotic stress similar to the osmoprotection in plants and bacteria. Our studies found a cell fate switch that enables cells to survive gradually changing stress environments by preventing caspase activation and protect cells through proline accumulation.

A rate threshold mechanism regulates MAPK stress signaling and survival

Cells are exposed to changes in extracellular stimulus concentration that vary as a function of rate. However, how cells integrate information conveyed from stimulation rate along with concentration remains poorly understood. Here, we examined how varying the rate of stress application alters budding yeast mitogen-activated protein kinase (MAPK) signaling and cell behavior at the single-cell level. We show that signaling depends on a rate threshold that operates in conjunction with stimulus concentration to determine the timing of MAPK signaling during rate-varying stimulus treatments. We also discovered that the stimulation rate threshold and stimulation rate-dependent cell survival are sensitive to changes in the expression levels of the Ptp2 phosphatase, but not of another phosphatase that similarly regulates osmostress signaling during switch-like treatments. Our results demonstrate that stimulation rate is a regulated determinant of cell behavior and provide a paradigm to guide the dissection of major stimulation rate dependent mechanisms in other systems.

Combined transcriptomic and phosphoproteomic analysis of BMP4 signaling in human embryonic stem cells

Human embryonic stem cells (hESCs) are an invaluable tool in the fields of embryology and regenerative medicine. Activin A and BMP4 are well-characterised growth factors implicated in pluripotency and differentiation. In the current study, hESCs are cultured in a modified version of mTeSR1, where low concentrations of ActivinA substitute for TGFβ. This culture system is further used to investigate the changes induced by BMP4 on hESCs by employing a combination of transcriptomic and phosphoproteomic approaches. Results indicate that in a pluripotent state, hESCs maintain WNT signaling under negative regulation by expressing pathway inhibitors. Initial stages of differentiation are characterized by upregulation of WNT pathway ligands, TGFβ pathway inhibitors which have been shown in Xenopus to expand the BMP signaling range essential for embryonic patterning, and mesendodermal transcripts. Moreover, BMP4 enhances the phosphorylation of proteins associated with migration and transcriptional regulation. Results further indicate the vital regulatory role of Activin A and BMP4 in crucial fate decisions in hESCs.

Development of a microfluidic approach for the real-time analysis of extrinsic TGF-β signalling

Autocrine and paracrine signalling are traditionally difficult to study due to the sub-micromolar concentrations involved. This has proven to be especially limiting in the study of embryonic stem cells that rely on such signalling for viability, self-renewal, and proliferation. Microfluidics allows to achieve local concentrations of ligands representative of the in vivo stem cell niche, gaining more precise control over the cell microenvironment, as well as to manipulate ligands availability with high temporal resolution and minimal amount of reagents. Here we developed a microfluidics-based system for monitoring the dynamics of TGF-β pathway activity by means of a SMAD2/3-dependent luciferase reporter. We first validated our system by showing dose-dependent transcriptional activation. We then tested the effects of pulsatile stimulation and delayed inhibition of TGF-β activity on signalling dynamics. Finally, we show that our microfluidic system, unlike conventional culture systems, can detect TGF-β ligands secreted in the conditioned medium from hESCs.

Diverse Cell Stimulation Kinetics Identify Predictive Signal Transduction Models

Computationally understanding the molecular mechanisms that give rise to cell signaling responses upon different environmental, chemical, and genetic perturbations is a long-standing challenge that requires models that fit and predict quantitative responses for new biological conditions. Overcoming this challenge depends not only on good models and detailed experimental data but also on the rigorous integration of both. We propose a quantitative framework to perturb and model generic signaling networks using multiple and diverse changing environments (hereafter "kinetic stimulations") resulting in distinct pathway activation dynamics. We demonstrate that utilizing multiple diverse kinetic stimulations better constrains model parameters and enables predictions of signaling dynamics that would be impossible using traditional dose-response or individual kinetic stimulations. To demonstrate our approach, we use experimentally identified models to predict signaling dynamics in normal, mutated, and drug-treated conditions upon multitudes of kinetic stimulations and quantify which proteins and reaction rates are most sensitive to which extracellular stimulations.

Information transmission from NFkB signaling dynamics to gene expression

The dynamic signal encoding paradigm suggests that information flows from the extracellular environment into specific signaling patterns (encoding) that are then read by downstream effectors to control cellular behavior. Previous work empirically quantified the information content of dynamic signaling patterns. However, whether this information can be faithfully transmitted to the gene expression level is unclear. Here we used NFkB signaling as a model to understand the accuracy of information transmission from signaling dynamics into gene expression. Using a detailed mathematical model, we simulated realistic NFkB signaling patterns with different degrees of variability. The NFkB patterns were used as an input to a simple gene expression model. Analysis of information transmission between ligand and NFkB and ligand and gene expression allows us to determine information loss in transmission between receptors to dynamic signaling patterns and between signaling dynamics to gene expression. Information loss could occur due to biochemical noise or due to a lack of specificity. We found that noise-free gene expression has very little information loss suggesting that gene expression can preserve specificity in NFkB patterns. As expected, the addition of noise to the gene expression model results in information loss. Interestingly, this effect can be mitigated by a specific choice of parameters that can substantially reduce information loss due to biochemical noise during gene expression. Overall our results show that the cellular capacity for information transmission from dynamic signaling patterns to gene expression can be high enough to preserve ligand specificity and thereby the accuracy of cellular response to environmental cues.

On-Target Anti-TGF-β Therapies Are Not Succeeding in Clinical Cancer Treatments: What Are Remaining Challenges?

Metastasis is the leading cause of death for cancer patients. During cancer progression, the initial detachment of cells from the primary tumor and the later colonization of a secondary organ are characterized as limiting steps for metastasis. Epithelial-mesenchymal transition (EMT) and mesenchymal-epithelial transition (MET) are opposite dynamic multistep processes that enable these critical events in metastasis by altering the phenotype of cancer cells and improving their ability to migrate, invade and seed at distant organs. Among the molecular pathways that promote tumorigenesis in late-stage cancers, transforming growth factor-β (TGF-β) is described as an EMT master inducer by controlling different genes and proteins related to cytoskeleton assembly, cell-cell attachment and extracellular matrix remodeling. Still, despite the successful outcomes of different TGF-β pharmacological inhibitors in cell culture (in vitro) and animal models (in vivo), results in cancer clinical trials are poor or inconsistent at least, highlighting the existence of crucial components in human cancers that have not been properly explored. Here we review most recent findings to provide perspectives bridging the gap between on-target anti-TGF-β therapies in vitro and in pre-clinical models and the poor clinical outcomes in treating cancer patients. Specifically, we focus on (i) the dual roles of TGF-β signaling in cancer metastasis; (ii) dynamic signaling; (iii) functional differences of TGF-β free in solution vs. in exosomes; (iv) the regulatory effects of tumor microenvironment (TME) – particularly by cancer-associated fibroblasts – on TGF-β signaling pathway. Clearly identifying and establishing those missing links may provide strategies to revitalize and clinically improve the efficacy of TGF-β targeted therapies.

A review on microfluidics manipulation of the extracellular chemical microenvironment and its emerging application to cell analysis

Spatiotemporal manipulation of extracellular chemical environments with simultaneous monitoring of cellular responses plays an essential role in exploring fundamental biological processes and expands our understanding of underlying mechanisms. Despite the rapid progress and promising successes in manipulation strategies, many challenges remain due to the small size of cells and the rapid diffusion of chemical molecules. Fortunately, emerging microfluidic technology has become a powerful approach for precisely controlling the extracellular chemical microenvironment, which benefits from its integration capacity, automation, and high-throughput capability, as well as its high resolution down to submicron. Here, we summarize recent advances in microfluidics manipulation of the extracellular chemical microenvironment, including the following aspects: i) Spatial manipulation of chemical microenvironments realized by convection flow-, diffusion-, and droplet-based microfluidics, and surface chemical modification; ii) Temporal manipulation of chemical microenvironments enabled by flow switching/shifting, moving/flowing cells across laminar flows, integrated microvalves/pumps, and droplet manipulation; iii) Spatiotemporal manipulation of chemical microenvironments implemented by a coupling strategy and open-space microfluidics; and iv) High-throughput manipulation of chemical microenvironments. Finally, we briefly present typical applications of the above-mentioned technical advances in cell-based analyses including cell migration, cell signaling, cell differentiation, multicellular analysis, and drug screening. We further discuss the future improvement of microfluidics manipulation of extracellular chemical microenvironments to fulfill the needs of biological and biomedical research and applications.

BMP Signaling Gradient Scaling in the Zebrafish Pectoral Fin

Secreted growth factors can act as morphogens that form spatial concentration gradients in developing organs, thereby controlling growth and patterning. For some morphogens, adaptation of the gradients to tissue size allows morphological patterns to remain proportioned as the organs grow. In the zebrafish pectoral fin, we found that BMP signaling forms a two-dimensional gradient. The length of the gradient scales with tissue length and its amplitude increases with fin size according to a power-law. Gradient scaling and amplitude power-laws are signatures of growth control by time derivatives of morphogenetic signaling: cell division correlates with the fold change over time of the cellular signaling levels. We show that Smoc1 regulates BMP gradient scaling and growth in the fin. Smoc1 scales the gradient by means of a feedback loop: Smoc1 is a BMP agonist and BMP signaling represses Smoc1 expression. Our work uncovers a layer of morphogen regulation during vertebrate appendage development.

Interplay between morphogen-directed positional information systems and physiological signaling

The development of an organism from an undifferentiated single cell into a spatially complex structure requires spatial patterning of cell fates across tissues. Positional information, proposed by Lewis Wolpert in 1969, has led to the characterization of many components involved in regulating morphogen signaling activity. However, how morphogen gradients are established, maintained, and interpreted by cells still is not fully understood. Quantitative and systems-based approaches are increasingly needed to define general biological design rules that govern positional information systems in developing organisms. This short review highlights a selective set of studies that have investigated the roles of physiological signaling in modulating and mediating morphogen-based pattern formation. Similarities between neural transmission and morphogen-based pattern formation mechanisms suggest underlying shared principles of active cell-based communication. Within larger tissues, neural networks provide directed information, via physiological signaling, that supplements positional information through diffusion. Further, mounting evidence demonstrates that physiological signaling plays a role in ensuring robustness of morphogen-based signaling. We conclude by highlighting several outstanding questions regarding the role of physiological signaling in morphogen-based pattern formation. Elucidating how physiological signaling impacts positional information is critical for understanding the close coupling of developmental and cellular processes in the context of development, disease, and regeneration.

Insights into mammalian morphogen dynamics from embryonic stem cell systems

Morphogens play an essential role in cell fate specification and patterning including in laying out the mammalian body plan during gastrulation. In vivo studies have shed light on the signaling pathways involved in this process and the phenotypes associated with their disruption, however, several important open questions remain regarding how morphogens function in space and time. Self-organized patterning systems based on embryonic stem cells have emerged as a powerful platform for beginning to address these questions that is complementary to in vivo approaches. Here we review recent progress in understanding morphogen signaling dynamics and patterning in early mammalian development by taking advantage of cutting-edge embryonic stem cell technology.

Cracking the context-specific PI3K signaling code

Specificity in signal transduction is determined by the ability of cells to "encode" and subsequently "decode" different environmental signals. Akin to computer software, this "signaling code" governs context-dependent execution of cellular programs through modulation of signaling dynamics and can be corrupted by disease-causing mutations. Class IA phosphoinositide 3-kinase (PI3K) signaling is critical for normal growth and development and is dysregulated in human disorders such as benign overgrowth syndromes, cancer, primary immune deficiency, and metabolic syndrome. Despite decades of PI3K research, understanding of context-dependent regulation of the PI3K pathway and of the underlying signaling code remains rudimentary. Here, we review current knowledge on context-specific PI3K signaling and how technological advances now make it possible to move from a qualitative to quantitative understanding of this pathway. Insight into how cellular PI3K signaling is encoded or decoded may open new avenues for rational pharmacological targeting of PI3K-associated diseases. The principles of PI3K context-dependent signal encoding and decoding described here are likely applicable to most, if not all, major cell signaling pathways.

On-Target Anti-TGF-β Therapies Are Not Succeeding in Clinical Cancer Treatments: What Are Remaining Challenges?

Metastasis is the leading cause of death for cancer patients. During cancer progression, the initial detachment of cells from the primary tumor and the later colonization of a secondary organ are characterized as limiting steps for metastasis. Epithelial-mesenchymal transition (EMT) and mesenchymal-epithelial transition (MET) are opposite dynamic multistep processes that enable these critical events in metastasis by altering the phenotype of cancer cells and improving their ability to migrate, invade and seed at distant organs. Among the molecular pathways that promote tumorigenesis in late-stage cancers, transforming growth factor-β (TGF-β) is described as an EMT master inducer by controlling different genes and proteins related to cytoskeleton assembly, cell-cell attachment and extracellular matrix remodeling. Still, despite the successful outcomes of different TGF-β pharmacological inhibitors in cell culture (in vitro) and animal models (in vivo), results in cancer clinical trials are poor or inconsistent at least, highlighting the existence of crucial components in human cancers that have not been properly explored. Here we review most recent findings to provide perspectives bridging the gap between on-target anti-TGF-β therapies in vitro and in pre-clinical models and the poor clinical outcomes in treating cancer patients. Specifically, we focus on (i) the dual roles of TGF-β signaling in cancer metastasis; (ii) dynamic signaling; (iii) functional differences of TGF-β free in solution vs. in exosomes; (iv) the regulatory effects of tumor microenvironment (TME) - particularly by cancer-associated fibroblasts - on TGF-β signaling pathway. Clearly identifying and establishing those missing links may provide strategies to revitalize and clinically improve the efficacy of TGF-β targeted therapies.

A matter of time: Formation and interpretation of the Bicoid morphogen gradient

Spatially distributed signaling molecules, known as morphogens, provide spatial information during development. A host of different morphogens have now been identified, from subcellular gradients through to morphogens that act across a whole embryo. These gradients form over a wide-range of timescales, from seconds to hours, and their time windows for interpretation are also highly variable; the processes of morphogen gradient formation and interpretation are highly dynamic. The morphogen Bicoid (Bcd), present in the early Drosophila embryo, is essential for setting up the future Drosophila body segments. Due to its accessibility for both genetic perturbations and imaging, this system has provided key insights into how precise patterning can occur within a highly dynamic system. Here, we review the temporal scales of Bcd gradient formation and interpretation. In particular, we discuss the quantitative evidence for different models of Bcd gradient formation, outline the time windows for Bcd interpretation, and describe how Bcd temporally adapts its own ability to be interpreted. The utilization of temporal information in morphogen readout may provide crucial inputs to ensure precise spatial patterning, particularly in rapidly developing systems.

Optogenetic approaches to investigate spatiotemporal signaling during development

Embryogenesis is coordinated by signaling pathways that pattern the developing organism. Many aspects of this process are not fully understood, including how signaling molecules spread through embryonic tissues, how signaling amplitude and dynamics are decoded, and how multiple signaling pathways cooperate to pattern the body plan. Optogenetic approaches can be used to address these questions by providing precise experimental control over a variety of biological processes. Here, we review how these strategies have provided new insights into developmental signaling and discuss how they could contribute to future investigations.

Signaling regulation during gastrulation: Insights from mouse embryos and in vitro systems

Gastrulation is the process whereby cells exit pluripotency and concomitantly acquire and pattern distinct cell fates. This is driven by the convergence of WNT, BMP, Nodal and FGF signals, which are tightly spatially and temporally controlled, resulting in regional and stage-specific signaling environments. The combination, level and duration of signals that a cell is exposed to, according its position within the embryo and the developmental time window, dictates the fate it will adopt. The key pathways driving gastrulation exhibit complex interactions, which are difficult to disentangle in vivo due to the complexity of manipulating multiple signals in parallel with high spatiotemporal resolution. Thus, our current understanding of the signaling dynamics regulating gastrulation is limited. In vitro stem cell models have been established, which undergo organized cellular differentiation and patterning. These provide amenable, simplified, deconstructed and scalable models of gastrulation. While the foundation of our understanding of gastrulation stems from experiments in embryos, in vitro systems are now beginning to reveal the intricate details of signaling regulation. Here we discuss the current state of knowledge of the role, regulation and dynamic interaction of signaling pathways that drive mouse gastrulation.

The calcium channel subunit α2δ-3 organizes synapses via an activity-dependent and autocrine BMP signaling pathway

Synapses are highly specialized for neurotransmitter signaling, yet activity-dependent growth factor release also plays critical roles at synapses. While efficient neurotransmitter signaling relies on precise apposition of release sites and neurotransmitter receptors, molecular mechanisms enabling high-fidelity growth factor signaling within the synaptic microenvironment remain obscure. Here we show that the auxiliary calcium channel subunit α2δ-3 promotes the function of an activity-dependent autocrine Bone Morphogenetic Protein (BMP) signaling pathway at the Drosophila neuromuscular junction (NMJ). α2δ proteins have conserved synaptogenic activity, although how they execute this function has remained elusive. We find that α2δ-3 provides an extracellular scaffold for an autocrine BMP signal, suggesting a mechanistic framework for understanding α2δ's conserved role in synapse organization. We further establish a transcriptional requirement for activity-dependent, autocrine BMP signaling in determining synapse density, structure, and function. We propose that activity-dependent, autocrine signals provide neurons with continuous feedback on their activity state for modulating both synapse structure and function.

Moran Model of Spatial Alignment in Microbial Colonies

We describe a spatial Moran model that captures mechanical interactions and directional growth in spatially extended populations. The model is analytically tractable and completely solvable under a mean-field approximation and can elucidate the mechanisms that drive the formation of population-level patterns. As an example we model a population of E. coli growing in a rectangular microfluidic trap. We show that spatial patterns can arise as a result of a tug-of-war between boundary effects and growth rate modulations due to cell-cell interactions: Cells align parallel to the long side of the trap when boundary effects dominate. However, when cell-cell interactions exceed a critical value, cells align orthogonally to the trap's long side. This modeling approach and analysis can be extended to directionally-growing cells in a variety of domains to provide insight into how local and global interactions shape collective behavior.