Delayed coupling theory of vertebrate segmentation

Modularity of the segmentation clock and morphogenesis

Vertebrates have evolved great diversity in the number of segments dividing the trunk body, however, the developmental origin of the evolvability of this trait is poorly understood. The number of segments is thought to be determined in embryogenesis as a product of morphogenesis of the pre-somitic mesoderm (PSM) and the periodicity of a molecular oscillator active within the PSM known as the segmentation clock. Here, we explore whether the clock and PSM morphogenesis exhibit developmental modularity, as independent evolution of these two processes may explain the high evolvability of segment number. Using a computational model of the clock and PSM parameterised for zebrafish, we find that the clock is broadly robust to variation in morphogenetic processes such as cell ingression, motility, compaction, and cell division. We show that this robustness is in part determined by the length of the PSM and the strength of phase coupling in the clock. As previous studies report no changes to morphogenesis upon perturbing the clock, we suggest that the clock and morphogenesis of the PSM exhibit developmental modularity.

Oscillatory DeltaC Expression in Neural Progenitors Primes the Prototype of Forebrain Development

Notch signaling plays a pivotal role in regulating various developmental processes, particularly in controlling the timing of neuronal production within the developing neocortex. Central to this regulatory mechanism is the oscillatory pattern of Delta, which functions as a developmental clock modulator. Its deficiency profoundly impairs mammalian brain formation, highlighting its fundamental role in brain development. However, zebrafish carrying a mutation in the functional ortholog DeltaC (dlc) within their functional ortholog exhibit an intact forebrain structure, implying evolutionary variations in Notch signaling within the forebrain. In this study, we unveil the distinct yet analogous expression profiles of Delta and Her genes in the developing vertebrate forebrain. Specifically, for the first time, we detected the oscillatory expression of the Delta gene dlc in the developing zebrafish forebrain. Although this oscillatory pattern appeared irregular and was not pervasive among the progenitor population, attenuation of the dlc-involved Notch pathway using a γ-secretase inhibitor impaired neuronal differentiation in the developing zebrafish forebrain, revealing the indispensable role of the dlc-involved Notch pathway in regulating early zebrafish neurogenesis. Taken together, our results demonstrate the foundational prototype of dlc-involved Notch signaling in the developing zebrafish forebrains, upon which the intricate patterns of the mammalian neocortex may have been sculpted.

Progress in understanding the vertebrate segmentation clock

The segmentation clock is a molecular oscillator that regulates the periodic formation of somites from the presomitic mesoderm during vertebrate embryogenesis. Synchronous oscillatory expression of a Hairy homologue or Hairy-related basic helix-loop-helix (bHLH) transcriptional repressor in presomitic mesoderm cells regulates periodic expression of downstream factors that control somite segmentation with a periodicity that varies across species. Although many of the key components of the clock have been identified and characterized, less is known about how the clock is synchronized across cells and how species-specific periodicity is achieved. Advances in live imaging, stem cell and organoid technologies, and synthetic approaches have started to uncover the detailed mechanisms underlying these aspects of somitogenesis, providing insight into how morphogenesis is coordinated in space and time during embryonic development.

Statistical description of mobile oscillators in embryonic pattern formation

Synchronization of mobile oscillators occurs in numerous contexts, including physical, chemical, biological, and engineered systems. In vertebrate embryonic development, a segmental body structure is generated by a population of mobile oscillators. Cells in this population produce autonomous gene expression rhythms and interact with their neighbors through local signaling. These cells form an extended tissue where frequency and cell mobility gradients coexist. Gene expression kinematic waves travel through this tissue and pattern the segment boundaries. It has been shown that oscillator mobility promotes global synchronization. However, in vertebrate segment formation, mobility may also introduce local fluctuations in kinematic waves and impair segment boundaries. Here, we derive a general framework for mobile oscillators that relates local mobility fluctuations to synchronization dynamics and pattern robustness. We formulate a statistical description of mobile phase oscillators in terms of probability density. We obtain and solve diffusion equations for the average phase and variance, revealing the relationship between local fluctuations and global synchronization in a homogeneous population of oscillators. Analysis of the probability density for large mobility identifies a mean-field onset, where locally coupled oscillators start behaving as if each oscillator was coupled with all the others. We extend the statistical description to inhomogeneous systems to address the gradients present in the vertebrate segmenting tissue. The theory relates pattern stability to mobility, coupling, and pattern wavelength. The general approach of the statistical description may be applied to mobile oscillators in other contexts, as well as to other patterning systems where mobility is present.

Cell-autonomous timing drives the vertebrate segmentation clock's wave pattern

Rhythmic and sequential segmentation of the growing vertebrate body relies on the segmentation clock, a multi-cellular oscillating genetic network. The clock is visible as tissue-level kinematic waves of gene expression that travel through the presomitic mesoderm (PSM) and arrest at the position of each forming segment. Here, we test how this hallmark wave pattern is driven by culturing single maturing PSM cells. We compare their cell-autonomous oscillatory and arrest dynamics to those we observe in the embryo at cellular resolution, finding similarity in the relative slowing of oscillations and arrest in concert with differentiation. This shows that cell-extrinsic signals are not required by the cells to instruct the developmental program underlying the wave pattern. We show that a cell-autonomous timing activity initiates during cell exit from the tailbud, then runs down in the anterior-ward cell flow in the PSM, thereby using elapsed time to provide positional information to the clock. Exogenous FGF lengthens the duration of the cell-intrinsic timer, indicating extrinsic factors in the embryo may regulate the segmentation clock via the timer. In sum, our work suggests that a noisy cell-autonomous, intrinsic timer drives the slowing and arrest of oscillations underlying the wave pattern, while extrinsic factors in the embryo tune this timer's duration and precision. This is a new insight into the balance of cell-intrinsic and -extrinsic mechanisms driving tissue patterning in development.

Dynamic processes of fate decision in inducible bistable systems

The process of biological fate decision regulated by gene regulatory networks involves numerous complex dynamical interactions among many components. Mathematical modeling typically employed ordinary differential equations and steady-state analysis, which has yielded valuable quantitative insights. However, stable states predicted by theoretical models often fail to capture transient or metastable phenomena that occur during most observation periods in experimental or real biological systems. We attribute this discrepancy to the omission of dynamic processes of various complex interactions. Here, we demonstrate the influence of delays in gene regulatory steps and the timescales of the external induction on the dynamic processes of the fate decision in inducible bistable systems. We propose that steady-state parameters determine the landscape of fate decision. However, during the dynamic evolution along the landscape, the unequal delays of biochemical interactions as well as the timescale of external induction cause deviations in the differentiation trajectories, leading to the formation of new transient distributions that persist long term. Our findings emphasize the importance of considering dynamic processes in fate decision instead of relying solely on steady-state analysis. We provide insights into the interpretation of experimental phenomena and offer valuable guidance for future efforts in dynamical modeling and synthetic biology design.

A mathematical model for pancreatic cancer during intraepithelial neoplasia

Cancer is the result of complex interactions of intrinsic and extrinsic cell processes, which promote sustained proliferation, resistance to apoptosis, reprogramming and reorganization. The evolution of any type of cancer emerges from the role of the microenvironmental conditions and their impact of some molecular complexes on certain signalling pathways. The understanding of the early onset of cancer requires a multiscale analysis of the cellular microenvironment. In this paper, we analyse a qualitative multiscale model of pancreatic adenocarcinoma by modelling the cellular microenvironment through elastic cell interactions and their intercellular communication mechanisms, such as growth factors and cytokines. We focus on the low-grade dysplasia (PanIN 1) and moderate dysplasia (PanIN 2) stages of pancreatic adenocarcinoma. To this end, we propose a gene-regulatory network associated with the processes of proliferation and apoptosis of pancreatic cells and its kinetics in terms of delayed differential equations to mimic cell development. Likewise, we couple the cell cycle with the spatial distribution of cells and the transport of growth factors to show that the adenocarcinoma evolution is triggered by inflammatory processes. We show that the oncogene RAS may be an important target for developing anti-inflammatory strategies that limit the emergence of more aggressive adenocarcinomas.

Unravelling differential Hes1 dynamics during axis elongation of mouse embryos through single-cell tracking

The intricate dynamics of Hes expression across diverse cell types in the developing vertebrate embryonic tail have remained elusive. To address this, we have developed an endogenously tagged Hes1-Achilles mouse line, enabling precise quantification of dynamics at the single-cell resolution across various tissues. Our findings reveal striking disparities in Hes1 dynamics between presomitic mesoderm (PSM) and preneural tube (pre-NT) cells. While pre-NT cells display variable, low-amplitude oscillations, PSM cells exhibit synchronized, high-amplitude oscillations. Upon the induction of differentiation, the oscillation amplitude increases in pre-NT cells. Additionally, our study of Notch inhibition on Hes1 oscillations unveils distinct responses in PSM and pre-NT cells, corresponding to differential Notch ligand expression dynamics. These findings suggest the involvement of separate mechanisms driving Hes1 oscillations. Thus, Hes1 demonstrates dynamic behaviour across adjacent tissues of the embryonic tail, yet the varying oscillation parameters imply differences in the information that can be transmitted by these dynamics.

Nonreciprocal synchronization in embryonic oscillator ensembles

Synchronization of coupled oscillators is a universal phenomenon encountered across different scales and contexts, e.g., chemical wave patterns, superconductors, and the unison applause we witness in concert halls. The existence of common underlying coupling rules defines universality classes, revealing a fundamental sameness between seemingly distinct systems. Identifying rules of synchronization in any particular setting is hence of paramount relevance. Here, we address the coupling rules within an embryonic oscillator ensemble linked to vertebrate embryo body axis segmentation. In vertebrates, the periodic segmentation of the body axis involves synchronized signaling oscillations in cells within the presomitic mesoderm (PSM), from which somites, the prevertebrae, form. At the molecular level, it is known that intact Notch-signaling and cell-to-cell contact are required for synchronization between PSM cells. However, an understanding of the coupling rules is still lacking. To identify these, we develop an experimental assay that enables direct quantification of synchronization dynamics within mixtures of oscillating cell ensembles, for which the initial input frequency and phase distribution are known. Our results reveal a "winner-takes-it-all" synchronization outcome, i.e., the emerging collective rhythm matches one of the input rhythms. Using a combination of theory and experimental validation, we develop a coupling model, the "Rectified Kuramoto" (ReKu) model, characterized by a phase-dependent, nonreciprocal interaction in the coupling of oscillatory cells. Such nonreciprocal synchronization rules reveal fundamental similarities between embryonic oscillators and a class of collective behaviors seen in neurons and fireflies, where higher-level computations are performed and linked to nonreciprocal synchronization.

Cellular and molecular control of vertebrate somitogenesis

Segmentation is a fundamental feature of the vertebrate body plan. This metameric organization is first implemented by somitogenesis in the early embryo, when paired epithelial blocks called somites are rhythmically formed to flank the neural tube. Recent advances in in vitro models have offered new opportunities to elucidate the mechanisms that underlie somitogenesis. Notably, models derived from human pluripotent stem cells introduced an efficient proxy for studying this process during human development. In this Review, we summarize the current understanding of somitogenesis gained from both in vivo studies and in vitro studies. We deconstruct the spatiotemporal dynamics of somitogenesis into four distinct modules: dynamic events in the presomitic mesoderm, segmental determination, somite anteroposterior polarity patterning, and epithelial morphogenesis. We first focus on the segmentation clock, as well as signalling and metabolic gradients along the tissue, before discussing the clock and wavefront and other models that account for segmental determination. We then detail the molecular and cellular mechanisms of anteroposterior polarity patterning and somite epithelialization.

The Clock and Wavefront Self-Organizing model recreates the dynamics of mouse somitogenesis in vivo and in vitro

During mouse development, presomitic mesoderm cells synchronize Wnt and Notch oscillations, creating sequential phase waves that pattern somites. Traditional somitogenesis models attribute phase waves to a global modulation of the oscillation frequency. However, increasing evidence suggests that they could arise in a self-organizing manner. Here, we introduce the Sevilletor, a novel reaction-diffusion system that serves as a framework to compare different somitogenesis patterning hypotheses. Using this framework, we propose the Clock and Wavefront Self-Organizing model that considers an excitable self-organizing region where phase waves form independent of global frequency gradients. The model recapitulates the change in relative phase of Wnt and Notch observed during mouse somitogenesis and provides a theoretical basis for understanding the excitability of mouse presomitic mesoderm cells in vitro.

Species-specific roles of the Notch ligands, receptors, and targets orchestrating the signaling landscape of the segmentation clock

Somitogenesis is a hallmark feature of all vertebrates and some invertebrate species that involves the periodic formation of block-like structures called somites. Somites are transient embryonic segments that eventually establish the entire vertebral column. A highly conserved molecular oscillator called the segmentation clock underlies this periodic event and the pace of this clock regulates the pace of somite formation. Although conserved signaling pathways govern the clock in most vertebrates, the mechanisms underlying the species-specific divergence in various clock characteristics remain elusive. For example, the segmentation clock in classical model species such as zebrafish, chick, and mouse embryos tick with a periodicity of ∼30, ∼90, and ∼120 min respectively. This enables them to form the species-specific number of vertebrae during their overall timespan of somitogenesis. Here, we perform a systematic review of the species-specific features of the segmentation clock with a keen focus on mouse embryos. We perform this review using three different perspectives: Notch-responsive clock genes, ligand-receptor dynamics, and synchronization between neighboring oscillators. We further review reports that use non-classical model organisms and in vitro model systems that complement our current understanding of the segmentation clock. Our review highlights the importance of comparative developmental biology to further our understanding of this essential developmental process.

Generation of patterns in the paraxial mesoderm

The Segmentation Clock is a tissue-level patterning system that enables the segmentation of the vertebral column precursors into transient multicellular blocks called somites. This patterning system comprises a set of elements that are essential for correct segmentation. Under the so-called "Clock and Wavefront" model, the system consists of two elements, a genetic oscillator that manifests itself as traveling waves of gene expression, and a regressing wavefront that transforms the temporally periodic signal encoded in the oscillations into a permanent spatially periodic pattern of somite boundaries. Over the last twenty years, every new discovery about the Segmentation Clock has been tightly linked to the nomenclature of the "Clock and Wavefront" model. This constrained allocation of discoveries into these two elements has generated long-standing debates in the field as what defines molecularly the wavefront and how and where the interaction between the two elements establishes the future somite boundaries. In this review, we propose an expansion of the "Clock and Wavefront" model into three elements, "Clock", "Wavefront" and signaling gradients. We first provide a detailed description of the components and regulatory mechanisms of each element, and we then examine how the spatiotemporal integration of the three elements leads to the establishment of the presumptive somite boundaries. To be as exhaustive as possible, we focus on the Segmentation Clock in zebrafish. Furthermore, we show how this three-element expansion of the model provides a better understanding of the somite formation process and we emphasize where our current understanding of this patterning system remains obscure.

Insect visuomotor delay adjustments in group flight support swarm cohesion

Flying insects routinely demonstrate coordinated flight in crowded assemblies despite strict communication and processing constraints. This study experimentally records multiple flying insects tracking a moving visual stimulus. System identification techniques are used to robustly identify the tracking dynamics, including a visuomotor delay. The population delay distributions are quantified for solo and group behaviors. An interconnected visual swarm model incorporating heterogeneous delays is developed, and bifurcation analysis and swarm simulation are applied to assess swarm stability under the delays. The experiment recorded 450 insect trajectories and quantified visual tracking delay variation. Solitary tasks showed a 30ms average delay and standard deviation of 50ms, while group behaviors show a 15ms average and 8ms standard deviation. Analysis and simulation indicate that the delay adjustments during group flight support swarm formation and center stability, and are robust to measurement noise. These results quantify the role of visuomotor delay heterogeneity in flying insects and their role in supporting swarm cohesion through implicit communication.

Ripply suppresses Tbx6 to induce dynamic-to-static conversion in somite segmentation

The metameric pattern of somites is created based on oscillatory expression of clock genes in presomitic mesoderm. However, the mechanism for converting the dynamic oscillation to a static pattern of somites is still unclear. Here, we provide evidence that Ripply/Tbx6 machinery is a key regulator of this conversion. Ripply1/Ripply2-mediated removal of Tbx6 protein defines somite boundary and also leads to cessation of clock gene expression in zebrafish embryos. On the other hand, activation of ripply1/ripply2 mRNA and protein expression is periodically regulated by clock oscillation in conjunction with an Erk signaling gradient. Whereas Ripply protein decreases rapidly in embryos, Ripply-triggered Tbx6 suppression persists long enough to complete somite boundary formation. Mathematical modeling shows that a molecular network based on results of this study can reproduce dynamic-to-static conversion in somitogenesis. Furthermore, simulations with this model suggest that sustained suppression of Tbx6 caused by Ripply is crucial in this conversion.

How clock heterogeneity affects synchronization and can enhance stability

The production process of integrated electronic circuitry inherently leads to large heterogeneities on the component level. For electronic clock networks this implies detuned intrinsic frequencies and differences in coupling strength and the characteristic time delays associated with signal transmission, processing, and feedback. Using a phase-model description, we study the effects of such component heterogeneity on the dynamical properties of synchronization in networks of mutually delay-coupled Kuramoto oscillators. We test the theory against experimental results and circuit-level simulations in a prototype system of mutually delay-coupled electronic clocks, so-called phase-locked loops. Interestingly, our results show that heterogeneity in the system can actually enhance the stability of synchronized states. That means that beyond the optimizations that can be achieved by tuning homogeneous coupling strengths, time delays, and loop-filter cut-off frequencies, heterogeneities in these system parameters enable much better optimization of perturbation decay rates, stabilization of synchronous states, and tuning of phase differences between the clocks. Our theory enables the design of custom-fit synchronization layers according to the specific requirements and properties of electronic systems, such as operational frequencies, phase relations, and, e.g., transmission delays. These results are not restricted to electronic systems, because signal transmission, processing, and feedback delays are common to networks of spatially distributed and coupled autonomous oscillators.

Heterogeneity-induced synchronization in delay-coupled electronic oscillators

We study synchronization in networks of delay-coupled electronic oscillators, so-called phase-locked loops (PLLs). Using a phase-model description, we study the collective dynamics of mutually coupled PLLs and report the phenomenon of heterogeneity-induced synchronization. This phenomenon refers to the observation that heterogeneity in the system's parameters can induce synchronization by stabilizing the states which are unstable without such heterogeneity. In systems where component heterogeneity can be tuned and controlled, we show how the complex collective self-organized dynamics can be guided towards synchronized states with specific operational frequencies and phase relations. This is of importance for the technical applicability of self-organized dynamics. In electrical engineering, for example, where components can be strongly heterogeneous, our theoretical framework can inform the design process for networks of spatially distributed PLLs. The results presented here are also useful in understanding the collective dynamics in ensembles of phase oscillators with time-delayed interactions, inertia, and heterogeneity.

Arnold tongue entrainment reveals dynamical principles of the embryonic segmentation clock

Living systems exhibit an unmatched complexity, due to countless, entangled interactions across scales. Here, we aim to understand a complex system, that is, segmentation timing in mouse embryos, without a reference to these detailed interactions. To this end, we develop a coarse-grained approach, in which theory guides the experimental identification of the segmentation clock entrainment responses. We demonstrate period- and phase-locking of the segmentation clock across a wide range of entrainment parameters, including higher-order coupling. These quantifications allow to derive the phase response curve (PRC) and Arnold tongues of the segmentation clock, revealing its essential dynamical properties. Our results indicate that the somite segmentation clock has characteristics reminiscent of a highly non-linear oscillator close to an infinite period bifurcation and suggests the presence of long-term feedbacks. Combined, this coarse-grained theoretical-experimental approach reveals how we can derive simple, essential features of a highly complex dynamical system, providing precise experimental control over the pace and rhythm of the somite segmentation clock.

Bilateral Feedback in Oscillator Model Is Required to Explain the Coupling Dynamics of Hes1 with the Cell Cycle

Biological processes are governed by the expression of proteins, and for some proteins, their level of expression can fluctuate periodically over time (i.e., they oscillate). Many oscillatory proteins (e.g., cell cycle proteins and those from the HES family of transcription factors) are connected in complex ways, often within large networks. This complexity can be elucidated by developing intuitive mathematical models that describe the underlying critical aspects of the relationships between these processes. Here, we provide a mathematical explanation of a recently discovered biological phenomenon: the phasic position of the gene Hes1’s oscillatory expression at the beginning of the cell cycle of an individual human breast cancer stem cell can have a predictive value on how long that cell will take to complete a cell cycle. We use a two-component model of coupled oscillators to represent Hes1 and the cell cycle in the same cell with minimal assumptions. Inputting only the initial phase angles, we show that this model is capable of predicting the dynamic mitosis to mitosis behaviour of Hes1 and predicting cell cycle length patterns as found in real-world experimental data. Moreover, we discover that bidirectional coupling between Hes1 and the cell cycle is critical within the system for the data to be reproduced and that nonfixed asymmetry in the interactions between the oscillators is required. The phase dynamics we present here capture the complex interplay between Hes1 and the cell cycle, helping to explain nongenetic cell cycle variability, which has critical implications in cancer treatment contexts.

Patterning with clocks and genetic cascades: Segmentation and regionalization of vertebrate versus insect body plans

Oscillatory and sequential processes have been implicated in the spatial patterning of many embryonic tissues. For example, molecular clocks delimit segmental boundaries in vertebrates and insects and mediate lateral root formation in plants, whereas sequential gene activities are involved in the specification of regional identities of insect neuroblasts, vertebrate neural tube, vertebrate limb, and insect and vertebrate body axes. These processes take place in various tissues and organisms, and, hence, raise the question of what common themes and strategies they share. In this article, we review 2 processes that rely on the spatial regulation of periodic and sequential gene activities: segmentation and regionalization of the anterior-posterior (AP) axis of animal body plans. We study these processes in species that belong to 2 different phyla: vertebrates and insects. By contrasting 2 different processes (segmentation and regionalization) in species that belong to 2 distantly related phyla (arthropods and vertebrates), we elucidate the deep logic of patterning by oscillatory and sequential gene activities. Furthermore, in some of these organisms (e.g., the fruit fly Drosophila), a mode of AP patterning has evolved that seems not to overtly rely on oscillations or sequential gene activities, providing an opportunity to study the evolution of pattern formation mechanisms.

Towards a physical understanding of developmental patterning

The temporal coordination of events at cellular and tissue scales is essential for the proper development of organisms, and involves cell-intrinsic processes that can be coupled by local cellular signalling and instructed by global signalling, thereby creating spatial patterns of cellular states that change over time. The timing and structure of these patterns determine how an organism develops. Traditional developmental genetic methods have revealed the complex molecular circuits regulating these processes but are limited in their ability to predict and understand the emergent spatio-temporal dynamics. Increasingly, approaches from physics are now being used to help capture the dynamics of the system by providing simplified, generic descriptions. Combined with advances in imaging and computational power, such approaches aim to provide insight into timing and patterning in developing systems.

Time and space in segmentation

Arthropod segmentation and vertebrate somitogenesis are leading fields in the experimental and theoretical interrogation of developmental patterning. However, despite the sophistication of current research, basic conceptual issues remain unresolved. These include: (i) the mechanistic origins of spatial organization within the segment addition zone (SAZ); (ii) the mechanistic origins of segment polarization; (iii) the mechanistic origins of axial variation; and (iv) the evolutionary origins of simultaneous patterning. Here, I explore these problems using coarse-grained models of cross-regulating dynamical processes. In the morphogenetic framework of a row of cells undergoing axial elongation, I simulate interactions between an 'oscillator', a 'switch' and up to three 'timers', successfully reproducing essential patterning behaviours of segmenting systems. By comparing the output of these largely cell-autonomous models to variants that incorporate positional information, I find that scaling relationships, wave patterns and patterning dynamics all depend on whether the SAZ is regulated by temporal or spatial information. I also identify three mechanisms for polarizing oscillator output, all of which functionally implicate the oscillator frequency profile. Finally, I demonstrate significant dynamical and regulatory continuity between sequential and simultaneous modes of segmentation. I discuss these results in the context of the experimental literature.

Hybrid reaction-diffusion and clock-and-wavefront model for the arrest of oscillations in the somitogenesis segmentation clock

The clock and wavefront paradigm is arguably the most widely accepted model for explaining the embryonic process of somitogenesis. According to this model, somitogenesis is based upon the interaction between a genetic oscillator, known as segmentation clock, and a differentiation wavefront, which provides the positional information indicating where each pair of somites is formed. Shortly after the clock and wavefront paradigm was introduced, Meinhardt presented a conceptually different mathematical model for morphogenesis in general, and somitogenesis in particular. Recently, Cotterell et al. [A local, self-organizing reaction-diffusion model can explain somite patterning in embryos, Cell Syst. 1, 257-269 (2015)] rediscovered an equivalent model by systematically enumerating and studying small networks performing segmentation. Cotterell et al. called it a progressive oscillatory reaction-diffusion (PORD) model. In the Meinhardt-PORD model, somitogenesis is driven by short-range interactions and the posterior movement of the front is a local, emergent phenomenon, which is not controlled by global positional information. With this model, it is possible to explain some experimental observations that are incompatible with the clock and wavefront model. However, the Meinhardt-PORD model has some important disadvantages of its own. Namely, it is quite sensitive to fluctuations and depends on very specific initial conditions (which are not biologically realistic). In this work, we propose an equivalent Meinhardt-PORD model and then amend it to couple it with a wavefront consisting of a receding morphogen gradient. By doing so, we get a hybrid model between the Meinhardt-PORD and the clock-and-wavefront ones, which overcomes most of the deficiencies of the two originating models.

A dynamic, spatially periodic, micro-pattern of HES5 underlies neurogenesis in the mouse spinal cord

Ultradian oscillations of HES Transcription Factors (TFs) at the single-cell level enable cell state transitions. However, the tissue-level organisation of HES5 dynamics in neurogenesis is unknown. Here, we analyse the expression of HES5 ex vivo in the developing mouse ventral spinal cord and identify microclusters of 4-6 cells with positively correlated HES5 level and ultradian dynamics. These microclusters are spatially periodic along the dorsoventral axis and temporally dynamic, alternating between high and low expression with a supra-ultradian persistence time. We show that Notch signalling is required for temporal dynamics but not the spatial periodicity of HES5. Few Neurogenin 2 cells are observed per cluster, irrespective of high or low state, suggesting that the microcluster organisation of HES5 enables the stable selection of differentiating cells. Computational modelling predicts that different cell coupling strengths underlie the HES5 spatial patterns and rate of differentiation, which is consistent with comparison between the motoneuron and interneuron progenitor domains. Our work shows a previously unrecognised spatiotemporal organisation of neurogenesis, emergent at the tissue level from the synthesis of single-cell dynamics.

Nonlinear delay differential equations and their application to modeling biological network motifs

Biological regulatory systems, such as cell signaling networks, nervous systems and ecological webs, consist of complex dynamical interactions among many components. Network motif models focus on small sub-networks to provide quantitative insight into overall behavior. However, such models often overlook time delays either inherent to biological processes or associated with multi-step interactions. Here we systematically examine explicit-delay versions of the most common network motifs via delay differential equation (DDE) models, both analytically and numerically. We find many broadly applicable results, including parameter reduction versus canonical ordinary differential equation (ODE) models, analytical relations for converting between ODE and DDE models, criteria for when delays may be ignored, a complete phase space for autoregulation, universal behaviors of feedforward loops, a unified Hill-function logic framework, and conditions for oscillations and chaos. We conclude that explicit-delay modeling simplifies the phenomenology of many biological networks and may aid in discovering new functional motifs.

Network motif models focus on small sub-networks in biological systems to quantitatively describe overall behavior but they often overlook time delays. Here, the authors systematically examine the most common network motifs via delay differential equations (DDE), often leading to more concise descriptions.

From local resynchronization to global pattern recovery in the zebrafish segmentation clock

Integrity of rhythmic spatial gene expression patterns in the vertebrate segmentation clock requires local synchronization between neighboring cells by Delta-Notch signaling and its inhibition causes defective segment boundaries. Whether deformation of the oscillating tissue complements local synchronization during patterning and segment formation is not understood. We combine theory and experiment to investigate this question in the zebrafish segmentation clock. We remove a Notch inhibitor, allowing resynchronization, and analyze embryonic segment recovery. We observe unexpected intermingling of normal and defective segments, and capture this with a new model combining coupled oscillators and tissue mechanics. Intermingled segments are explained in the theory by advection of persistent phase vortices of oscillators. Experimentally observed changes in recovery patterns are predicted in the theory by temporal changes in tissue length and cell advection pattern. Thus, segmental pattern recovery occurs at two length and time scales: rapid local synchronization between neighboring cells, and the slower transport of the resulting patterns across the tissue through morphogenesis.

Geometric models for robust encoding of dynamical information into embryonic patterns

During development, cells gradually assume specialized fates via changes of transcriptional dynamics, sometimes even within the same developmental stage. For anterior-posterior (AP) patterning in metazoans, it has been suggested that the gradual transition from a dynamic genetic regime to a static one is encoded by different transcriptional modules. In that case, the static regime has an essential role in pattern formation in addition to its maintenance function. In this work, we introduce a geometric approach to study such transition. We exhibit two types of genetic regime transitions arising through local or global bifurcations, respectively. We find that the global bifurcation type is more generic, more robust, and better preserves dynamical information. This could parsimoniously explain common features of metazoan segmentation, such as changes of periods leading to waves of gene expressions, ‘speed/frequency-gradient’ dynamics, and changes of wave patterns. Geometric approaches appear as possible alternatives to gene regulatory networks to understand development.

Spatial correlations in a finite-range Kuramoto model

We study spatial correlations between oscillator phases in the steady state of a Kuramoto model, in which phase oscillators that are randomly distributed in space interact with constant strength but within a limited range. Such a model could be relevant, for example, in the synchronization of gene expression oscillations in cells, where only oscillations of neighboring cells are coupled through cell-cell contacts. We analytically infer spatial phase-phase correlation functions from the known steady-state distribution of oscillators for the case of homogenous frequencies and show that these can contain information about the range and strength of interactions, provided the noise in the system can be estimated. We suggest a method for the latter, and also explore when correlations appear to be ergodic in this model, which would enable an experimental measurement of correlation functions to utilize temporal averages. Simulations show that our techniques also give qualitative results for the model with heterogenous frequencies. We illustrate our results by comparison with experimental data on genetic oscillations in the segmentation clock of zebrafish embryos.

Waiting on the Fringe: cell autonomy and signaling delays in segmentation clocks

The rhythmic and sequential segmentation of the vertebrate body axis into somites during embryogenesis is governed by a multicellular, oscillatory patterning system called the segmentation clock. Despite many overt similarities between vertebrates, differences in genetic and dynamic regulation have been reported, raising intriguing questions about the evolution and conservation of this fundamental patterning process. Recent studies have brought insights into two important and related issues: (1) whether individual cells of segmentation clocks are autonomous oscillators or require cell-cell communication for their rhythm; and (2) the role of delays in the cell-cell communication that synchronizes the population of genetic oscillators. Although molecular details differ between species, conservation may exist at the level of the dynamics, hinting at rules for evolutionary trajectories in the system.

Segmentation clock dynamics is strongly synchronized in the forming somite

During vertebrate somitogenesis an inherent segmentation clock coordinates the spatiotemporal signaling to generate segmented structures that pattern the body axis. Using our experimental and quantitative approach, we study the cell movements and the genetic oscillations of her1 expression level at single-cell resolution simultaneously and scale up to the entire pre-somitic mesoderm (PSM) tissue. From the experimentally determined phases of PSM cellular oscillators, we deduced an in vivo frequency profile gradient along the anterior-posterior PSM axis and inferred precise mathematical relations between spatial cell-level period and tissue-level somitogenesis period. We also confirmed a gradient in the relative velocities of cellular oscillators along the axis. The phase order parameter within an ensemble of oscillators revealed the degree of synchronization in the tailbud and the posterior PSM being only partial, whereas synchronization can be almost complete in the presumptive somite region but with temporal oscillations. Collectively, the degree of synchronization itself, possibly regulated by cell movement and the synchronized temporal phase of the transiently expressed clock protein Her1, can be an additional control mechanism for making precise somite boundaries.

Coupling delay controls synchronized oscillation in the segmentation clock

Individual cellular activities fluctuate but are constantly coordinated at the population level via cell-cell coupling. A notable example is the somite segmentation clock, in which the expression of clock genes (such as Hes7) oscillates in synchrony between the cells that comprise the presomitic mesoderm (PSM)^1,2. This synchronization depends on the Notch signalling pathway; inhibiting this pathway desynchronizes oscillations, leading to somite fusion^3-7. However, how Notch signalling regulates the synchronicity of HES7 oscillations is unknown. Here we establish a live-imaging system using a new fluorescent reporter (Achilles), which we fuse with HES7 to monitor synchronous oscillations in HES7 expression in the mouse PSM at a single-cell resolution. Wild-type cells can rapidly correct for phase fluctuations in HES7 oscillations, whereas the absence of the Notch modulator gene lunatic fringe (Lfng) leads to a loss of synchrony between PSM cells. Furthermore, HES7 oscillations are severely dampened in individual cells of Lfng-null PSM. However, when Lfng-null PSM cells were completely dissociated, the amplitude and periodicity of HES7 oscillations were almost normal, which suggests that LFNG is involved mostly in cell-cell coupling. Mixed cultures of control and Lfng-null PSM cells, and an optogenetic Notch signalling reporter assay, revealed that LFNG delays the signal-sending process of intercellular Notch signalling transmission. These results-together with mathematical modelling-raised the possibility that Lfng-null PSM cells shorten the coupling delay, thereby approaching a condition known as the oscillation or amplitude death of coupled oscillators⁸. Indeed, a small compound that lengthens the coupling delay partially rescues the amplitude and synchrony of HES7 oscillations in Lfng-null PSM cells. Our study reveals a delay control mechanism of the oscillatory networks involved in somite segmentation, and indicates that intercellular coupling with the correct delay is essential for synchronized oscillation.

Intrinsic noise, Delta-Notch signalling and delayed reactions promote sustained, coherent, synchronized oscillations in the presomitic mesoderm

Using a stochastic individual-based modelling approach, we examine the role that Delta-Notch signalling plays in the regulation of a robust and reliable somite segmentation clock. We find that not only can Delta-Notch signalling synchronize noisy cycles of gene expression in adjacent cells in the presomitic mesoderm (as is known), but it can also amplify and increase the coherence of these cycles. We examine some of the shortcomings of deterministic approaches to modelling these cycles and demonstrate how intrinsic noise can play an active role in promoting sustained oscillations, giving rise to noise-induced quasi-cycles. Finally, we explore how translational/transcriptional delays can result in the cycles in neighbouring cells oscillating in anti-phase and we study how this effect relates to the propagation of noise-induced stochastic waves.

Cell cycle regulation of oscillations yields coupling of growth and form in a computational model of the presomitic mesoderm

A striking example of coupling between growth and form occurs during the segmentation of the vertebrate embryo. During segmentation, pairs of segments, one on either side of the anterior-posterior axis, bud off from the presomitic mesoderm (PSM) at regular intervals in time. In the clock and wavefront model, a multicellular oscillator regulates the time at which the next pair of segments form whilst a wavefront regulates their spatial location. In most mathematical models of segmentation, it is assumed that cells in the PSM are oscillators that have a constant natural frequency. Based on recent experimental findings, here we propose a model in which the natural oscillation frequency of each PSM cell is a function of its position in the cell cycle. Given adequate oscillator coupling and that cells in the PSM are randomly distributed in the cell cycle, we find that the emergent oscillator frequency is a weighted average of the constituent oscillator frequencies with the weightings dependent on the fraction of cells in a given cell cycle state. Here, we show that such a model can allow for coupling between pattern formation and growth rate in PSM tissue.

Arthropod segmentation

There is now compelling evidence that many arthropods pattern their segments using a clock-and-wavefront mechanism, analogous to that operating during vertebrate somitogenesis. In this Review, we discuss how the arthropod segmentation clock generates a repeating sequence of pair-rule gene expression, and how this is converted into a segment-polarity pattern by ‘timing factor’ wavefronts associated with axial extension. We argue that the gene regulatory network that patterns segments may be relatively conserved, although the timing of segmentation varies widely, and double-segment periodicity appears to have evolved at least twice. Finally, we describe how the repeated evolution of a simultaneous (Drosophila-like) mode of segmentation within holometabolan insects can be explained by heterochronic shifts in timing factor expression plus extensive pre-patterning of the pair-rule genes.

Synchronization dynamics of mobile oscillators in the presence of coupling delays

Individual biological oscillators can synchronize to generate a collective rhythm. During vertebrate development, mobile cells exchange signals to synchronize a rhythmic pattern generator that makes the embryonic segments. Previous theoretical works have shown that cell mobility can enhance synchronization of coupled oscillators when signal exchange is instantaneous. However, in vertebrate segmentation, the exchange of signals is thought to comprise delays from signal sending and processing, which could alter the effect of mobility on synchronization. Here, we study synchronization dynamics of mobile phase oscillators in the presence of coupling delays. We find that mobility can speed up synchronization when coupling delays are present. We derive an analytical expression for the characteristic time of synchronization dynamics, which is in very good agreement with numerical simulations. This analytical expression suggests a subdivision of the mobility range into different dynamical regimes and reveals that, with delayed coupling, synchronization is enhanced at a lower mobility rate than with instantaneous coupling. We argue that these results may be relevant to the synchronization of mobile oscillators in vertebrate segmentation.

Around the clock: gradient shape and noise impact the evolution of oscillatory segmentation dynamics

BACKGROUND

Segmentation, the subdivision of the major body axis into repeated elements, is considered one of the major evolutionary innovations in bilaterian animals. In all three segmented animal clades, the predominant segmentation mechanism is sequential segmentation, where segments are generated one by one in anterior-posterior order from a posterior undifferentiated zone. In vertebrates and arthropods, sequential segmentation is thought to arise from a clock-and-wavefront-type mechanism, where oscillations in the posterior growth zone are transformed into a segmental prepattern in the anterior by a receding wavefront. Previous evo-devo simulation studies have demonstrated that this segmentation type repeatedly arises, supporting the idea of parallel evolutionary origins in these animal clades. Sequential segmentation has been studied most extensively in vertebrates, where travelling waves have been observed that reflect the slowing down of oscillations prior to their cessation and where these oscillations involve a highly complex regulatory network. It is currently unclear under which conditions this oscillator complexity and slowing should be expected to evolve, how they are related and to what extent similar properties should be expected for sequential segmentation in other animal species.

RESULTS

To investigate these questions, we extend a previously developed computational model for the evolution of segmentation. We vary the slope of the posterior morphogen gradient and the strength of gene expression noise. We find that compared to a shallow gradient, a steep morphogen gradient allows for faster evolution and evolved oscillator networks are simpler. Furthermore, under steep gradients, damped oscillators often evolve, whereas shallow gradients appear to require persistent oscillators which are regularly accompanied by travelling waves, indicative of a frequency gradient. We show that gene expression noise increases the likelihood of evolving persistent oscillators under steep gradients and of evolving frequency gradients under shallow gradients. Surprisingly, we find that the evolutions of oscillator complexity and travelling waves are not correlated, suggesting that these properties may have evolved separately.

CONCLUSIONS

Based on our findings, we suggest that travelling waves may have evolved in response to shallow morphogen gradients and gene expression noise. These two factors may thus also be responsible for the observed differences between different species within both the arthropod and chordate phyla.

Information flow in the presence of cell mixing and signaling delays during embryonic development

Embryonic morphogenesis is organized by an interplay between intercellular signaling and cell movements. Both intercellular signaling and cell movement involve multiple timescales. A key timescale for signaling is the time delay caused by preparation of signaling molecules and integration of received signals into cells' internal state. Movement of cells relative to their neighbors may introduce exchange of positions between cells during signaling. When cells change their relative positions in a tissue, the impact of signaling delays on intercellular signaling increases because the delayed information that cells receive may significantly differ from the present state of the tissue. The time it takes to perform a neighbor exchange sets a timescale of cell mixing that may be important for the outcome of signaling. Here we review recent theoretical work on the interplay of timescales between cell mixing and signaling delays adopting the zebrafish segmentation clock as a model system. We discuss how this interplay can lead to spatial patterns of gene expression that could disrupt the normal formation of segment boundaries in the embryo. The effect of cell mixing and signaling delays highlights the importance of theoretical and experimental frameworks to understand collective cellular behaviors arising from the interplay of multiple timescales in embryonic developmental processes.

Size-reduced embryos reveal a gradient scaling-based mechanism for zebrafish somite formation

Little is known about how the sizes of animal tissues are controlled. A prominent example is somite size, which varies widely both within an individual and across species. Despite intense study of the segmentation clock governing the timing of somite generation, how it relates to somite size is poorly understood. Here, we examine somite scaling and find that somite size at specification scales with the length of the presomitic mesoderm (PSM) despite considerable variation in PSM length across developmental stages and in surgically size-reduced embryos. Measurement of clock period, axis elongation speed and clock gene expression patterns demonstrate that existing models fail to explain scaling. We posit a 'clock and scaled gradient' model, in which somite boundaries are set by a dynamically scaling signaling gradient across the PSM. Our model not only explains existing data, but also makes a unique prediction that we confirm experimentally - the formation of periodic 'echoes' in somite size following perturbation of the size of one somite. Our findings demonstrate that gradient scaling plays a central role in both progression and size control of somitogenesis.

Travelling waves in somitogenesis: Collective cellular properties emerge from time-delayed juxtacrine oscillation coupling

The sculpturing of the vertebrate body plan into segments begins with the sequential formation of somites in the presomitic mesoderm (PSM). The rhythmicity of this process is controlled by travelling waves of gene expression. These kinetic waves emerge from coupled cellular oscillators and sweep across the PSM. In zebrafish, the oscillations are driven by autorepression of her genes and are synchronized via Notch signalling. Mathematical modelling has played an important role in explaining how collective properties emerge from the molecular interactions. Increasingly more quantitative experimental data permits the validation of those mathematical models, yet leads to increasingly more complex model formulations that hamper an intuitive understanding of the underlying mechanisms. Here, we review previous efforts, and design a mechanistic model of the her1 oscillator, which represents the experimentally viable her7;hes6 double mutant. This genetically simplified system is ideally suited to conceptually recapitulate oscillatory entrainment and travelling wave formation, and to highlight open questions. It shows that three key parameters, the autorepression delay, the juxtacrine coupling delay, and the coupling strength, are sufficient to understand the emergence of the collective period, the collective amplitude, and the synchronization of neighbouring Her1 oscillators. Moreover, two spatiotemporal time delay gradients, in the autorepression and in the juxtacrine signalling, are required to explain the collective oscillatory dynamics and synchrony of PSM cells. The highlighted developmental principles likely apply more generally to other developmental processes, including neurogenesis and angiogenesis.

Chemical event chain model of coupled genetic oscillators

We introduce a stochastic model of coupled genetic oscillators in which chains of chemical events involved in gene regulation and expression are represented as sequences of Poisson processes. We characterize steady states by their frequency, their quality factor, and their synchrony by the oscillator cross correlation. The steady state is determined by coupling and exhibits stochastic transitions between different modes. The interplay of stochasticity and nonlinearity leads to isolated regions in parameter space in which the coupled system works best as a biological pacemaker. Key features of the stochastic oscillations can be captured by an effective model for phase oscillators that are coupled by signals with distributed delays.

A damped oscillator imposes temporal order on posterior gap gene expression in Drosophila

Insects determine their body segments in two different ways. Short-germband insects, such as the flour beetle Tribolium castaneum, use a molecular clock to establish segments sequentially. In contrast, long-germband insects, such as the vinegar fly Drosophila melanogaster, determine all segments simultaneously through a hierarchical cascade of gene regulation. Gap genes constitute the first layer of the Drosophila segmentation gene hierarchy, downstream of maternal gradients such as that of Caudal (Cad). We use data-driven mathematical modelling and phase space analysis to show that shifting gap domains in the posterior half of the Drosophila embryo are an emergent property of a robust damped oscillator mechanism, suggesting that the regulatory dynamics underlying long- and short-germband segmentation are much more similar than previously thought. In Tribolium, Cad has been proposed to modulate the frequency of the segmentation oscillator. Surprisingly, our simulations and experiments show that the shift rate of posterior gap domains is independent of maternal Cad levels in Drosophila. Our results suggest a novel evolutionary scenario for the short- to long-germband transition and help explain why this transition occurred convergently multiple times during the radiation of the holometabolan insects.

Different insect species exhibit one of two distinct modes of determining their body segments (known as segmentation) during development: they either use a molecular oscillator to position segments sequentially, or they generate segments simultaneously through a hierarchical gene-regulatory cascade. The sequential mode is ancestral, while the simultaneous mode has been derived from it independently several times during evolution. In this paper, we present evidence suggesting that simultaneous segmentation also involves an oscillator in the posterior end of the embryo of the vinegar fly, Drosophila melanogaster. This surprising result indicates that both modes of segment determination are much more similar than previously thought. Such similarity provides an important step towards our understanding of the frequent evolutionary transitions observed between sequential and simultaneous segmentation.

Recent advances in understanding vertebrate segmentation

Segmentation is the partitioning of the body axis into a series of repeating units or segments. This widespread body plan is found in annelids, arthropods, and chordates, showing it to be a successful developmental strategy for growing and generating diverse morphology and anatomy. Segmentation has been extensively studied over the years. Forty years ago, Cooke and Zeeman published the Clock and Wavefront model, creating a theoretical framework of how developing cells could acquire and keep temporal and spatial information in order to generate a segmented pattern. Twenty years later, in 1997, Palmeirim and co-workers found the first clock gene whose oscillatory expression pattern fitted within Cooke and Zeeman's model. Currently, in 2017, new experimental techniques, such as new ex vivo experimental models, real-time imaging of gene expression, live single cell tracking, and simplified transgenics approaches, are revealing some of the fine details of the molecular processes underlying the inner workings of the segmentation mechanisms, bringing new insights into this fundamental process. Here we review and discuss new emerging views that further our understanding of the vertebrate segmentation clock, with a particular emphasis on recent publications that challenge and/or complement the currently accepted Clock and Wavefront model.

Mobility-induced persistent chimera states

We study the dynamics of mobile, locally coupled identical oscillators in the presence of coupling delays. We find different kinds of chimera states in which coherent in-phase and antiphase domains coexist with incoherent domains. These chimera states are dynamic and can persist for long times for intermediate mobility values. We discuss the mechanisms leading to the formation of these chimera states in different mobility regimes. This finding could be relevant for natural and technological systems composed of mobile communicating agents.

Excitable Dynamics and Yap-Dependent Mechanical Cues Drive the Segmentation Clock

The periodic segmentation of the vertebrate body axis into somites, and later vertebrae, relies on a genetic oscillator (the segmentation clock) driving the rhythmic activity of signaling pathways in the presomitic mesoderm (PSM). To understand whether oscillations are an intrinsic property of individual cells or represent a population-level phenomenon, we established culture conditions for stable oscillations at the cellular level. This system was used to demonstrate that oscillations are a collective property of PSM cells that can be actively triggered in vitro by a dynamical quorum sensing signal involving Yap and Notch signaling. Manipulation of Yap-dependent mechanical cues is sufficient to predictably switch isolated PSM cells from a quiescent to an oscillatory state in vitro, a behavior reminiscent of excitability in other systems. Together, our work argues that the segmentation clock behaves as an excitable system, introducing a broader paradigm to study such dynamics in vertebrate morphogenesis.

Speed regulation of genetic cascades allows for evolvability in the body plan specification of insects

During the anterior-posterior fate specification of insects, anterior fates arise in a nonelongating tissue (called the "blastoderm"), and posterior fates arise in an elongating tissue (called the "germband"). However, insects differ widely in the extent to which anterior-posterior fates are specified in the blastoderm versus the germband. Here we present a model in which patterning in both the blastoderm and germband of the beetle Tribolium castaneum is based on the same flexible mechanism: a gradient that modulates the speed of a genetic cascade of gap genes, resulting in the induction of sequential kinematic waves of gap gene expression. The mechanism is flexible and capable of patterning both elongating and nonelongating tissues, and hence converting blastodermal to germband fates and vice versa. Using RNAi perturbations, we found that blastodermal fates could be shifted to the germband, and germband fates could be generated in a blastoderm-like morphology. We also suggest a molecular mechanism underlying our model, in which gradient levels regulate the switch between two enhancers: One enhancer is responsible for sequential gene activation, and the other is responsible for freezing temporal rhythms into spatial patterns. This model is consistent with findings in Drosophila melanogaster, where gap genes were found to be regulated by two nonredundant "shadow" enhancers.

Stochastic Kuramoto oscillators with discrete phase states

We present a generalization of the Kuramoto phase oscillator model in which phases advance in discrete phase increments through Poisson processes, rendering both intrinsic oscillations and coupling inherently stochastic. We study the effects of phase discretization on the synchronization and precision properties of the coupled system both analytically and numerically. Remarkably, many key observables such as the steady-state synchrony and the quality of oscillations show distinct extrema while converging to the classical Kuramoto model in the limit of a continuous phase. The phase-discretized model provides a general framework for coupled oscillations in a Markov chain setting.

A framework for quantification and physical modeling of cell mixing applied to oscillator synchronization in vertebrate somitogenesis

In development and disease, cells move as they exchange signals. One example is found in vertebrate development, during which the timing of segment formation is set by a 'segmentation clock', in which oscillating gene expression is synchronized across a population of cells by Delta-Notch signaling. Delta-Notch signaling requires local cell-cell contact, but in the zebrafish embryonic tailbud, oscillating cells move rapidly, exchanging neighbors. Previous theoretical studies proposed that this relative movement or cell mixing might alter signaling and thereby enhance synchronization. However, it remains unclear whether the mixing timescale in the tissue is in the right range for this effect, because a framework to reliably measure the mixing timescale and compare it with signaling timescale is lacking. Here, we develop such a framework using a quantitative description of cell mixing without the need for an external reference frame and constructing a physical model of cell movement based on the data. Numerical simulations show that mixing with experimentally observed statistics enhances synchronization of coupled phase oscillators, suggesting that mixing in the tailbud is fast enough to affect the coherence of rhythmic gene expression. Our approach will find general application in analyzing the relative movements of communicating cells during development and disease.

Determining the impact of cell mixing on signaling during development

Cell movement and intercellular signaling occur simultaneously to organize morphogenesis during embryonic development. Cell movement can cause relative positional changes between neighboring cells. When intercellular signals are local such cell mixing may affect signaling, changing the flow of information in developing tissues. Little is known about the effect of cell mixing on intercellular signaling in collective cellular behaviors and methods to quantify its impact are lacking. Here we discuss how to determine the impact of cell mixing on cell signaling drawing an example from vertebrate embryogenesis: the segmentation clock, a collective rhythm of interacting genetic oscillators. We argue that comparing cell mixing and signaling timescales is key to determining the influence of mixing. A signaling timescale can be estimated by combining theoretical models with cell signaling perturbation experiments. A mixing timescale can be obtained by analysis of cell trajectories from live imaging. After comparing cell movement analyses in different experimental settings, we highlight challenges in quantifying cell mixing from embryonic timelapse experiments, especially a reference frame problem due to embryonic motions and shape changes. We propose statistical observables characterizing cell mixing that do not depend on the choice of reference frames. Finally, we consider situations in which both cell mixing and signaling involve multiple timescales, precluding a direct comparison between single characteristic timescales. In such situations, physical models based on observables of cell mixing and signaling can simulate the flow of information in tissues and reveal the impact of observed cell mixing on signaling.