Single-cell-resolution imaging of the impact of Notch signaling and mitosis on segmentation clock dynamics

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.

DeltaC and DeltaD ligands play different roles in the segmentation clock dynamics

The vertebrate segmentation clock drives periodic somite segmentation during embryonic development. Her1 and Her7 clock proteins generate oscillatory expression of their own genes as well as that of deltaC in zebrafish. In turn, DeltaC and DeltaD ligands activate Notch signaling, which then activates transcription of clock genes in neighboring cells. While DeltaC and DeltaD proteins form homo- and heterodimers, only DeltaC-containing oscillatory dimers were expected to be functional. To investigate the contributions of DeltaC and DeltaD proteins on the transcription of her1 and her7 segmentation clock genes, we counted their transcripts by performing single molecule fluorescent in situ hybridization imaging in different genetic backgrounds of zebrafish embryos. Surprisingly, we found that DeltaD homodimers are also functional. We further found that Notch signaling promotes transcription of both deltaC and deltaD genes, thereby creating a previously unnoticed positive feedback loop. Our computational model highlighted the intriguing differential roles of DeltaC and DeltaD dimers on the clock synchronization and transcript numbers, respectively. We anticipate that a mechanistic understanding of the Notch signaling pathway will not only shed light on the mechanism driving robust somite segmentation but also inspire similar quantitative studies in other tissues and organs.

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.

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.

Substrate Rigidity Modulates Segmentation Clock Dynamics in Isolated Presomitic Mesoderm Cells

The segmentation clock, a genetic oscillator in the presomitic mesoderm (PSM), is known to be influenced by biochemical signals, yet its potential regulation by mechanical cues remains unclear. The complex PSM microenvironment has made it challenging to isolate the effects of mechanical perturbations on clock behavior. Here we investigated how mechanical stimuli affect clock oscillations by culturing zebrafish PSM cells on PDMS micropost arrays with tunable rigidities (0.6-1200 kPa). We observed an inverse sigmoidal relationship between surface rigidity and both the percentage of oscillating cells and the number of oscillation cycles, with a switching threshold between 3-6 kPa. The periods of oscillating cells showed a consistently broad distribution across rigidity changes. Moreover, these cells exhibited distinct biophysical properties, such as reduced motility, contractility, and sustained circularity. These findings highlight the crucial role of cell-substrate interactions in regulating segmentation clock behavior, providing insights into the mechanobiology of somitogenesis.

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.

Cell coupling compensates for changes in single-cell Her6 dynamics and provides phenotypic robustness

This paper investigates the effect of altering the protein expression dynamics of the bHLH transcription factor Her6 at the single-cell level in the embryonic zebrafish telencephalon. Using a homozygote endogenous Her6:Venus reporter and 4D single-cell tracking, we show that Her6 oscillates in neural telencephalic progenitors and that the fusion of protein destabilisation (PEST) domain alters its expression dynamics, causing most cells to downregulate Her6 prematurely. However, counterintuitively, oscillatory cells increase, with some expressing Her6 at high levels, resulting in increased heterogeneity of Her6 expression in the population. These tissue-level changes appear to be an emergent property of coupling between single-cells, as revealed by experimentally disrupting Notch signalling and by computationally modelling alterations in Her6 protein stability. Despite the profound differences in the single-cell Her6 dynamics, the size of the telencephalon is only transiently altered and differentiation markers do not exhibit significant differences early on; however, a small increase is observed at later developmental stages. Our study suggests that cell coupling provides a compensation strategy, whereby an almost normal phenotype is maintained even though single-cell gene expression dynamics are abnormal, granting phenotypic robustness.

Oscillatory control of embryonic development

Proper embryonic development depends on the timely progression of a genetic program. One of the key mechanisms for achieving precise control of developmental timing is to use gene expression oscillations. In this Review, we examine how gene expression oscillations encode temporal information during vertebrate embryonic development by discussing the gene expression oscillations occurring during somitogenesis, neurogenesis, myogenesis and pancreas development. These oscillations play important but varied physiological functions in different contexts. Oscillations control the period of somite formation during somitogenesis, whereas they regulate the proliferation-to-differentiation switch of stem cells and progenitor cells during neurogenesis, myogenesis and pancreas development. We describe the similarities and differences of the expression pattern in space (i.e. whether oscillations are synchronous or asynchronous across neighboring cells) and in time (i.e. different time scales) of mammalian Hes/zebrafish Her genes and their targets in different tissues. We further summarize experimental evidence for the functional role of their oscillations. Finally, we discuss the outstanding questions for future research.

Microtubule polarity determines the lineage of embryonic neural precursor in zebrafish spinal cord

The phenomenal diversity of neuronal types in the central nervous system is achieved in part by the asymmetric division of neural precursors. In zebrafish neural precursors, asymmetric dispatch of Sara endosomes (with its Notch signaling cargo) functions as fate determinant which mediates asymmetric division. Here, we found two distinct pools of neural precursors based on Sara endosome inheritance and spindle-microtubule enrichment. Symmetric or asymmetric levels of spindle-microtubules drive differently Sara endosomes inheritance and predict neural precursor lineage. We uncover that CAMSAP2a/CAMSAP3a and KIF16Ba govern microtubule asymmetry and endosome motility, unveiling the heterogeneity of neural precursors. Using a plethora of physical and cell biological assays, we determined the physical parameters and molecular mechanisms behind microtubule asymmetries and biased endosome motility. Evolutionarily, the values of those parameters explain why all sensory organ precursor cells are asymmetric in flies while, in zebrafish spinal cord, two populations of neural precursors (symmetric vs asymmetric) are possible.

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.

Spatiotemporal control of pattern formation during somitogenesis

Spatiotemporal patterns widely occur in biological, chemical, and physical systems. Particularly, embryonic development displays a diverse gamut of repetitive patterns established in many tissues and organs. Branching treelike structures in lungs, kidneys, livers, pancreases, and mammary glands as well as digits and bones in appendages, teeth, and palates are just a few examples. A fascinating instance of repetitive patterning is the sequential segmentation of the primary body axis, which is conserved in all vertebrates and many arthropods and annelids. In these species, the body axis elongates at the posterior end of the embryo containing an unsegmented tissue. Meanwhile, segments sequentially bud off from the anterior end of the unsegmented tissue, laying down an exquisite repetitive pattern and creating a segmented body plan. In vertebrates, the paraxial mesoderm is sequentially divided into somites. In this review, we will discuss the most prominent models, the most puzzling experimental data, and outstanding questions in vertebrate somite segmentation.

Single-cell transcriptional dynamics in a living vertebrate

The ability to quantify transcriptional dynamics in individual cells via live imaging has revolutionized our understanding of gene regulation. However, such measurements are lacking in the context of vertebrate embryos. We addressed this deficit by applying MS2-MCP mRNA labeling to the quantification of transcription in zebrafish, a model vertebrate. We developed a platform of transgenic organisms, light sheet fluorescence microscopy, and optimized image analysis that enables visualization and quantification of MS2 reporters. We used these tools to obtain the first single-cell, real-time measurements of transcriptional dynamics of the segmentation clock. Our measurements challenge the traditional view of smooth clock oscillations and instead suggest a model of discrete transcriptional bursts that are organized in space and time. Together, these results highlight how measuring single-cell transcriptional activity can reveal unexpected features of gene regulation and how this data can fuel the dialogue between theory and experiment.

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.

Stochastic gene expression and environmental stressors trigger variable somite segmentation phenotypes

Mutations of several genes cause incomplete penetrance and variable expressivity of phenotypes, which are usually attributed to modifier genes or gene-environment interactions. Here, we show stochastic gene expression underlies the variability of somite segmentation defects in embryos mutant for segmentation clock genes her1 or her7. Phenotypic strength is further augmented by low temperature and hypoxia. By performing live imaging of the segmentation clock reporters, we further show that groups of cells with higher oscillation amplitudes successfully form somites while those with lower amplitudes fail to do so. In unfavorable environments, the number of cycles with high amplitude oscillations and the number of successful segmentations proportionally decrease. These results suggest that individual oscillation cycles stochastically fail to pass a threshold amplitude, resulting in segmentation defects in mutants. Our quantitative methodology is adaptable to investigate variable phenotypes of mutant genes in different tissues.

Nfinder: automatic inference of cell neighborhood in 2D and 3D using nuclear markers

BACKGROUND

In tissues and organisms, the coordination of neighboring cells is essential to maintain their properties and functions. Therefore, knowing which cells are adjacent is crucial to understand biological processes that involve physical interactions among them, e.g. cell migration and proliferation. In addition, some signaling pathways, such as Notch or extrinsic apoptosis, are highly dependent on cell-cell communication. While this is straightforward to obtain from membrane images, nuclei labelling is much more ubiquitous for technical reasons. However, there are no automatic and robust methods to find neighboring cells based only on nuclear markers.

RESULTS

In this work, we describe Nfinder, a method to assess the cell's local neighborhood from images with nuclei labeling. To achieve this goal, we approximate the cell-cell interaction graph by the Delaunay triangulation of nuclei centroids. Then, links are filtered by automatic thresholding in cell-cell distance (pairwise interaction) and the maximum angle that a pair of cells subtends with shared neighbors (non-pairwise interaction). We systematically characterized the detection performance by applying Nfinder to publicly available datasets from Drosophila melanogaster, Tribolium castaneum, Arabidopsis thaliana and C. elegans. In each case, the result of the algorithm was compared to a cell neighbor graph generated by manually annotating the original dataset. On average, our method detected 95% of true neighbors, with only 6% of false discoveries. Remarkably, our findings indicate that taking into account non-pairwise interactions might increase the Positive Predictive Value up to + 11.5%.

CONCLUSION

Nfinder is the first robust and automatic method for estimating neighboring cells in 2D and 3D based only on nuclear markers and without any free parameters. Using this tool, we found that taking non-pairwise interactions into account improves the detection performance significantly. We believe that using our method might improve the effectiveness of other workflows to study cell-cell interactions from microscopy images. Finally, we also provide a reference implementation in Python and an easy-to-use napari plugin.

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.

Controlling human organoid symmetry breaking reveals signaling gradients drive segmentation clock waves

Axial development of mammals involves coordinated morphogenetic events, including axial elongation, somitogenesis, and neural tube formation. To gain insight into the signals controlling the dynamics of human axial morphogenesis, we generated axially elongating organoids by inducing anteroposterior symmetry breaking of spatially coupled epithelial cysts derived from human pluripotent stem cells. Each organoid was composed of a neural tube flanked by presomitic mesoderm sequentially segmented into somites. Periodic activation of the somite differentiation gene MESP2 coincided in space and time with anteriorly traveling segmentation clock waves in the presomitic mesoderm of the organoids, recapitulating critical aspects of somitogenesis. Timed perturbations demonstrated that FGF and WNT signaling play distinct roles in axial elongation and somitogenesis, and that FGF signaling gradients drive segmentation clock waves. By generating and perturbing organoids that robustly recapitulate the architecture of multiple axial tissues in human embryos, this work offers a means to dissect mechanisms underlying human embryogenesis.

Orchestration of tissue shape changes and gene expression patterns in development

In development, tissue shape changes and gene expression patterns give rise to morphogenesis. Understanding tissue shape changes requires the analysis of mechanical properties of the tissue such as tissue rigidity, cell influx from neighboring tissues, cell shape changes and cell proliferation. Local and global gene expression patterns can be influenced by neighbor exchange and tissue shape changes. Here we review recent studies on the mechanisms for tissue elongation and its influences on dynamic gene expression patterns by focusing on vertebrate somitogenesis. We first introduce mechanical and biochemical properties of the segmenting tissue that drive tissue elongation. Then, we discuss patterning in the presence of cell mixing, scaling of signaling gradients, and dynamic phase waves of rhythmic gene expression under tissue shape changes. We also highlight the importance of theoretical approaches to address the relation between tissue shape changes and patterning.

Keeping development on time: Insights into post-transcriptional mechanisms driving oscillatory gene expression during vertebrate segmentation

Biological time keeping, or the duration and tempo at which biological processes occur, is a phenomenon that drives dynamic molecular and morphological changes that manifest throughout many facets of life. In some cases, the molecular mechanisms regulating the timing of biological transitions are driven by genetic oscillations, or periodic increases and decreases in expression of genes described collectively as a "molecular clock." In vertebrate animals, molecular clocks play a crucial role in fundamental patterning and cell differentiation processes throughout development. For example, during early vertebrate embryogenesis, the segmentation clock regulates the patterning of the embryonic mesoderm into segmented blocks of tissue called somites, which later give rise to axial skeletal muscle and vertebrae. Segmentation clock oscillations are characterized by rapid cycles of mRNA and protein expression. For segmentation clock oscillations to persist, the transcript and protein molecules of clock genes must be short-lived. Faithful, rhythmic, genetic oscillations are sustained by precise regulation at many levels, including post-transcriptional regulation, and such mechanisms are essential for proper vertebrate development. This article is categorized under: RNA Export and Localization > RNA Localization RNA Turnover and Surveillance > Regulation of RNA Stability Translation > Regulation.

Periodic inhibition of Erk activity drives sequential somite segmentation

Sequential segmentation creates modular body plans of diverse metazoan embryos^1-4. Somitogenesis establishes the segmental pattern of the vertebrate body axis. A molecular segmentation clock in the presomitic mesoderm sets the pace of somite formation⁴. However, how cells are primed to form a segment boundary at a specific location remains unclear. Here we developed precise reporters for the clock and double-phosphorylated Erk (ppErk) gradient in zebrafish. We show that the Her1-Her7 oscillator drives segmental commitment by periodically lowering ppErk, therefore projecting its oscillation onto the ppErk gradient. Pulsatile inhibition of the ppErk gradient can fully substitute for the role of the clock, and kinematic clock waves are dispensable for sequential segmentation. The clock functions upstream of ppErk, which in turn enables neighbouring cells to discretely establish somite boundaries in zebrafish⁵. Molecularly divergent clocks and morphogen gradients were identified in sequentially segmenting species^3,4,6-8. Our findings imply that versatile clocks may establish sequential segmentation in diverse species provided that they inhibit gradients.

The vertebrate Embryo Clock: Common players dancing to a different beat

Vertebrate embryo somitogenesis is the earliest morphological manifestation of the characteristic patterned structure of the adult axial skeleton. Pairs of somites flanking the neural tube are formed periodically during early development, and the molecular mechanisms in temporal control of this early patterning event have been thoroughly studied. The discovery of a molecular Embryo Clock (EC) underlying the periodicity of somite formation shed light on the importance of gene expression dynamics for pattern formation. The EC is now known to be present in all vertebrate organisms studied and this mechanism was also described in limb development and stem cell differentiation. An outstanding question, however, remains unanswered: what sets the different EC paces observed in different organisms and tissues? This review aims to summarize the available knowledge regarding the pace of the EC, its regulation and experimental manipulation and to expose new questions that might help shed light on what is still to unveil.

Imaging the onset of oscillatory signaling dynamics during mouse embryo gastrulation

A fundamental requirement for embryonic development is the coordination of signaling activities in space and time. A notable example in vertebrate embryos is found during somitogenesis, where gene expression oscillations linked to the segmentation clock are synchronized across cells in the presomitic mesoderm (PSM) and result in tissue-level wave patterns. To examine their onset during mouse embryo development, we studied the dynamics of the segmentation clock gene Lfng during gastrulation. To this end, we established an imaging setup using selective plane illumination microscopy (SPIM) that enables culture and simultaneous imaging of up to four embryos (‘SPIM- for-4’). Using SPIM-for-4, combined with genetically encoded signaling reporters, we detected the onset of Lfng oscillations within newly formed mesoderm at presomite stages. Functionally, we found that initial synchrony and the first ∼6-8 oscillation cycles occurred even when Notch signaling was impaired, revealing similarities to previous findings made in zebrafish embryos. Finally, we show that a spatial period gradient is present at the onset of oscillatory activity, providing a potential mechanism accounting for our observation that wave patterns build up gradually over the first oscillation cycles.

Summary: A versatile light-sheet imaging setup enabling simultaneous live imaging of multiple mouse embryos for 48 h, an approach that offers insight into the onset of oscillatory signaling dynamics and the segmentation clock.

Signaling oscillations in embryonic development

Tight spatiotemporal control of cellular behavior and cell fate decisions is paramount to the formation of multicellular organisms during embryonic development. Intercellular communication via signaling pathways mediates this control. Interestingly, these signaling pathways are not static, but dynamic and change in activity over time. Signaling oscillations as a specific type of dynamics are found in various signaling pathways and model systems. Functions of oscillations include the regulation of periodic events or the transmission of information by encoding signals in the dynamic properties of a signaling pathway. For instance, signaling oscillations in neural or pancreatic progenitor cells modulate their proliferation and differentiation. Oscillations between neighboring cells can also be synchronized, leading to the emergence of waves traveling through the tissue. Such population-wide signaling oscillations regulate for example the consecutive segmentation of vertebrate embryos, a process called somitogenesis. Here, we outline our current understanding of signaling oscillations in embryonic development, how signaling oscillations are generated, how they are studied and how they contribute to the regulation of embryonic development.

Morphogen-regulated contact-mediated signaling between cells can drive the transitions underlying body segmentation in vertebrates

We propose a unified mechanism that reproduces the sequence of dynamical transitions observed during somitogenesis, the process of body segmentation during embryonic development, that is invariant across all vertebrate species. This is achieved by combining inter-cellular interactions mediated via receptor-ligand coupling with global spatial heterogeneity introduced through a morphogen gradient known to occur along the anteroposterior axis. Our model reproduces synchronized oscillations in the gene expression in cells at the anterior of the presomitic mesoderm as it grows by adding new cells at its posterior, followed by travelling waves and subsequent arrest of activity, with the eventual appearance of somite-like patterns. This framework integrates a boundary-organized pattern formation mechanism, which uses positional information provided by a morphogen gradient, with the coupling-mediated self-organized emergence of collective dynamics, to explain the processes that lead to segmentation.

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.

Current approaches to fate mapping and lineage tracing using image data

Visualizing, tracking and reconstructing cell lineages in developing embryos has been an ongoing effort for well over a century. Recent advances in light microscopy, labelling strategies and computational methods to analyse complex image datasets have enabled detailed investigations into the fates of cells. Combined with powerful new advances in genomics and single-cell transcriptomics, the field of developmental biology is able to describe the formation of the embryo like never before. In this Review, we discuss some of the different strategies and applications to lineage tracing in live-imaging data and outline software methodologies that can be applied to various cell-tracking challenges.

Positional information encoded in the dynamic differences between neighboring oscillators during vertebrate segmentation

A central problem in developmental biology is to understand how cells interpret their positional information to give rise to spatial patterns, such as the process of periodic segmentation of the vertebrate embryo into somites. For decades, somite formation has been interpreted according to the clock-and-wavefront model. In this conceptual framework, molecular oscillators set the frequency of somite formation while the positional information is encoded in signaling gradients. Recent experiments using ex vivo explants have challenged this interpretation, suggesting that positional information is encoded in the properties of the oscillators, independent of long-range modulations such as signaling gradients. Here, we propose that positional information is encoded in the difference in the levels of neighboring oscillators. The differences gradually increase because both the amplitude and the period of the oscillators increase with time. When this difference exceeds a certain threshold, the segmentation program starts. Using this framework, we quantitatively fit experimental data from in vivo and ex vivo mouse segmentation, and propose mechanisms of somite scaling. Our results suggest a novel mechanism of spatial pattern formation based on the local interactions between dynamic molecular oscillators.

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.

Cell Tracking for Organoids: Lessons From Developmental Biology

Organoids have emerged as powerful model systems to study organ development and regeneration at the cellular level. Recently developed microscopy techniques that track individual cells through space and time hold great promise to elucidate the organizational principles of organs and organoids. Applied extensively in the past decade to embryo development and 2D cell cultures, cell tracking can reveal the cellular lineage trees, proliferation rates, and their spatial distributions, while fluorescent markers indicate differentiation events and other cellular processes. Here, we review a number of recent studies that exemplify the power of this approach, and illustrate its potential to organoid research. We will discuss promising future routes, and the key technical challenges that need to be overcome to apply cell tracking techniques to organoid biology.

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.

Imaging and manipulating the segmentation clock

During embryogenesis, the processes that control how cells differentiate and interact to form particular tissues and organs with precise timing and shape are of fundamental importance. One prominent example of such processes is vertebrate somitogenesis, which is governed by a molecular oscillator called the segmentation clock. The segmentation clock system is initiated in the presomitic mesoderm in which a set of genes and signaling pathways exhibit coordinated spatiotemporal dynamics to establish regularly spaced boundaries along the body axis; these boundaries provide a blueprint for the development of segment-like structures such as spines and skeletal muscles. The highly complex and dynamic nature of this in vivo event and the design principles and their regulation in both normal and abnormal embryogenesis are not fully understood. Recently, live-imaging has been used to quantitatively analyze the dynamics of selected components of the circuit, particularly in combination with well-designed experiments to perturb the system. Here, we review recent progress from studies using live imaging and manipulation, including attempts to recapitulate the segmentation clock in vitro. In combination with mathematical modeling, these techniques have become essential for disclosing novel aspects of the clock.

Pairing of segmentation clock genes drives robust pattern formation

Gene expression is an inherently stochastic process ^1,2; however, organismal development and homeostasis require that cells spatiotemporally coordinate the expression of large sets of genes. Coexpressed gene pairs in metazoans often reside in the same chromosomal neighborhood, with gene pairs representing 10% - 50% of all genes depending on species ^3–6. As shared upstream regulators can ensure correlated gene expression, the selective advantage of maintaining adjacent gene pairs remains unknown ⁶. Here, using two linked zebrafish segmentation clock genes, her1 and her7, and combining single-cell transcript counting, genetic engineering, real-time imaging and computational modeling, we reveal that gene pairing boosts correlated transcription and provides phenotypic robustness for developmental pattern formation. Our results demonstrate that disrupting gene pairing disrupts oscillations and segmentation, identifying the selective pressure retaining correlated transcription to sustain a robust and rapid developmental clock. We anticipate that these findings will inspire investigating advantages of gene pairing in other systems and engineering precise synthetic clocks in embryos and organoids.

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.

Understanding the Principles of Pattern Formation Driven by Notch Signaling by Integrating Experiments and Theoretical Models

Notch signaling is an evolutionary conserved cell-cell communication pathway. Besides regulating cell-fate decisions at an individual cell level, Notch signaling coordinates the emergent spatiotemporal patterning in a tissue through ligand-receptor interactions among transmembrane molecules of neighboring cells, as seen in embryonic development, angiogenesis, or wound healing. Due to its ubiquitous nature, Notch signaling is also implicated in several aspects of cancer progression, including tumor angiogenesis, stemness of cancer cells and cellular invasion. Here, we review experimental and computational models that help understand the operating principles of cell patterning driven by Notch signaling. First, we discuss the basic mechanisms of spatial patterning via canonical lateral inhibition and lateral induction mechanisms, including examples from angiogenesis, inner ear development and cancer metastasis. Next, we analyze additional layers of complexity in the Notch pathway, including the effect of varying cell sizes and shapes, ligand-receptor binding within the same cell, variable binding affinity of different ligand/receptor subtypes, and filopodia. Finally, we discuss some recent evidence of mechanosensitivity in the Notch pathway in driving collective epithelial cell migration and cardiovascular morphogenesis.

Dynamic properties of noise and Her6 levels are optimized by miR-9, allowing the decoding of the Her6 oscillator

Noise is prevalent in biology and has been widely quantified using snapshot measurements. This static view obscures our understanding of dynamic noise properties and how these affect gene expression and cell state transitions. Using a CRISPR/Cas9 Zebrafish her6::Venus reporter combined with mathematical and in vivo experimentation, we explore how noise affects the protein dynamics of Her6, a basic helix-loop-helix transcriptional repressor. During neurogenesis, Her6 expression transitions from fluctuating to oscillatory at single-cell level. We identify that absence of miR-9 input generates high-frequency noise in Her6 traces, inhibits the transition to oscillatory protein expression and prevents the downregulation of Her6. Together, these impair the upregulation of downstream targets and cells accumulate in a normally transitory state where progenitor and early differentiation markers are co-expressed. Computational modelling and double smFISH of her6 and the early neurogenesis marker, elavl3, suggest that the change in Her6 dynamics precedes the downregulation in Her6 levels. This sheds light onto the order of events at the moment of cell state transition and how this is influenced by the dynamic properties of noise. Our results suggest that Her/Hes oscillations, facilitated by dynamic noise optimization by miR-9, endow progenitor cells with the ability to make a cell state transition.

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.

Pumilio response and AU-rich elements drive rapid decay of Pnrc2-regulated cyclic gene transcripts

Vertebrate segmentation is regulated by the segmentation clock, a biological oscillator that controls periodic formation of somites, or embryonic segments, which give rise to many mesodermal tissue types. This molecular oscillator generates cyclic gene expression with the same periodicity as somite formation in the presomitic mesoderm (PSM), an area of mesenchymal cells that give rise to mature somites. Molecular components of the clock include the Hes/her family of genes that encode transcriptional repressors, but additional genes cycle. Cyclic gene transcripts are cleared rapidly, and clearance depends upon the pnrc2 (proline-rich nuclear receptor co-activator 2) gene that encodes an mRNA decay adaptor. Previously, we showed that the her1 3'UTR confers instability to otherwise stable transcripts in a Pnrc2-dependent manner, however, the molecular mechanism(s) by which cyclic gene transcripts are cleared remained largely unknown. To identify features of the her1 3'UTR that are critical for Pnrc2-mediated decay, we developed an array of transgenic inducible reporter lines carrying different regions of the 3'UTR. We find that the terminal 179 nucleotides (nts) of the her1 3'UTR are necessary and sufficient to confer rapid instability. Additionally, we show that the 3'UTR of another cyclic gene, deltaC (dlc), also confers Pnrc2-dependent instability. Motif analysis reveals that both her1 and dlc 3'UTRs contain terminally-located Pumilio response elements (PREs) and AU-rich elements (AREs), and we show that the PRE and ARE in the last 179 ​nts of the her1 3'UTR drive rapid turnover of reporter mRNA. Finally, we show that mutation of Pnrc2 residues and domains that are known to facilitate interaction of human PNRC2 with decay factors DCP1A and UPF1 reduce the ability of Pnrc2 to restore normal cyclic gene expression in pnrc2 mutant embryos. Our findings suggest that Pnrc2 interacts with decay machinery components and cooperates with Pumilio (Pum) proteins and ARE-binding proteins to promote rapid turnover of cyclic gene transcripts during somitogenesis.

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.

Somite boundary determination in normal and clock-less vertebrate embryos

Vertebrate segments called somites are generated by periodic segmentation of the presomitic mesoderm (PSM). In the most accepted theoretical model for somite segmentation, the clock and wavefront (CW) model, a clock that ticks to determine particular timings and a wavefront that moves posteriorly are presented in the PSM, and somite positions are determined when the clock meets the posteriorly moving wavefront somewhere in the PSM. Over the last two decades, it has been revealed that the molecular mechanism of the clock and wavefront in vertebrates is based on clock genes including Hes family transcription factors and Notch effectors that oscillate within the PSM to determine particular timings and fibroblast growth factor (FGF) gradients, acting as the posteriorly moving wavefront to determine the position of somite segmentation. A clock-less condition in the CW model was predicted to form no somites; however, irregularly sized somites were still formed in mice and zebrafish, suggesting that this was one of the limitations of the CW model. Recently, we performed interdisciplinary research of experimental and theoretical biological studies and revealed the mechanisms of somite boundary determination in normal and clock-less conditions by characterization of the FGF/extracellular signal-regulated kinase (ERK) activity dynamics. Since features of the molecular clock have already been described in-depth in several reviews, we summarized recent findings regarding the role of FGF/ERK signaling in somite boundary formation and described our current understanding of how FGF/ERK signaling contributes to somitogenesis in normal and clock-less conditions in this review.

In vitro characterization of the human segmentation clock

The segmental organization of the vertebral column is established early in embryogenesis when pairs of somites are rhythmically produced by the presomitic mesoderm (PSM). The tempo of somite formation is controlled by a molecular oscillator known as the segmentation clock^1,2. While this oscillator has been well-characterized in model organisms^1,2, whether a similar oscillator exists in humans remains unknown. Genetic analysis of patients with severe spine segmentation defects have implicated several human orthologs of cyclic genes associated with the mouse segmentation clock, suggesting that this oscillator might be conserved in humans³. Here we show that in vitro-derived human as well as mouse PSM cells⁴ recapitulate oscillations of the segmentation clock. Human PSM cells oscillate twice slower than mouse cells (5-hours vs. 2.5 hours), but are similarly regulated by FGF, Wnt, Notch and YAP⁵. Single cell RNA-sequencing reveals that mouse and human PSM cells in vitro follow a similar developmental trajectory to mouse PSM in vivo. Furthermore, we demonstrate that FGF signaling controls the phase and period of oscillations, expanding the role of this pathway beyond its classical interpretation in “Clock and Wavefront” models. Overall, our work identifying the human segmentation clock represents an important breakthrough for human developmental biology.

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.

Patterning via local cell-cell interactions in developing systems

Spatial patterning during embryonic development emerges from the differentiation of progenitor cells that share the same genetic program. One of the main challenges in systems biology is to understand the relationship between gene network and patterning, especially how the cells communicate to coordinate their differentiation. This review aims to describe the principles of pattern formation from local cell-cell interactions mediated by the Notch signalling pathway. Notch mediates signalling via direct cell-cell contact and regulates cell fate decisions in many tissues during embryonic development. Here, I will describe the patterning mechanisms via different Notch ligands and the critical role of Notch oscillations during the segmentation of the vertebrate body, brain development, and blood vessel formation.

Noise in the Vertebrate Segmentation Clock Is Boosted by Time Delays but Tamed by Notch Signaling

Taming cell-to-cell variability in gene expression is critical for precise pattern formation during embryonic development. To investigate the source and buffering mechanism of expression variability, we studied a biological clock, the vertebrate segmentation clock, controlling the precise spatiotemporal patterning of the vertebral column. By counting single transcripts of segmentation clock genes in zebrafish, we show that clock genes have low RNA amplitudes and expression variability is primarily driven by gene extrinsic sources, which is suppressed by Notch signaling. We further show that expression noise surprisingly increases from the posterior progenitor zone to the anterior segmentation and differentiation zone. Our computational model reproduces the spatial noise profile by incorporating spatially increasing time delays in gene expression. Our results, suggesting that expression variability is controlled by the balance of time delays and cell signaling in a vertebrate tissue, will shed light on the accuracy of natural clocks in multi-cellular systems and inspire engineering of robust synthetic oscillators.

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.