Synchronized oscillation of the segmentation clock gene in vertebrate development

Computational characterization of recombinase circuits for periodic behaviors

Recombinases are site-specific proteins found in nature that are capable of rearranging DNA. This function has made them promising gene editing tools in synthetic biology, as well as key elements in complex artificial gene circuits implementing Boolean logic. However, since DNA rearrangement is irreversible, it is still unclear how to use recombinases to build dynamic circuits like oscillators. In addition, this goal is challenging because a few molecules of recombinase are enough for promoter inversion, generating inherent stochasticity at low copy number. Here, we propose six different circuit designs for recombinase-based oscillators operating at a single copy number. We model them in a stochastic setting, leveraging the Gillespie algorithm for extensive simulations, and show that they can yield coherent periodic behaviors. Our results support the experimental realization of recombinase-based oscillators and, more generally, the use of recombinases to generate dynamic behaviors in synthetic biology.

Collective Oscillations in Coupled-Cell Systems

We investigate oscillations in coupled systems. The methodology is based on the Hopf bifurcation theorem and a condition extended from the Routh–Hurwitz criterion. Such a condition leads to locating the bifurcation values of the parameters. With such an approach, we analyze a single-cell system modeling the minimal genetic negative feedback loop and the coupled-cell system composed by these single-cell systems. We study the oscillatory properties for these systems and compare these properties between the model with Hill-type repression and the one with protein-sequestration-based repression. As the parameters move from the Hopf bifurcation value for single cells to the one for coupled cells, we compute the eigenvalues of the linearized systems to obtain the magnitude of the collective frequency when the periodic solution of the coupled-cell system is generated. Extending from this information on the parameter values, we further compute and compare the collective frequency for the coupled-cell system and the average frequency of the decoupled individual cells. To compare these scenarios with other biological oscillators, we perform parallel analysis and computations on a segmentation clock model.

Exploring the mechanism of pancreatic cell fate decisions via cell-cell communication

The endocrine and exocrine cells in pancreas originate initially from a group of apparently identical endoderm cells in the early gut. The endocrine and exocrine tissues are composed of islet/acinar and duct cells respectively. To explore the mechanism of pancreas cell fate decisions, we first construct a minimal mathematical model related to pancreatic regulations. The regulatory mechanism of acinar-to-islet cell conversion is revealed by bifurcation analysis of the model. In addition, Notch signaling is critical in determining the fate of endocrine and exocrine in the developing pancreas and it is a typical mediator of lateral inhibition which instructs adjacent cells to make different fate decisions. Next, we construct a multicellular model of cell-cell communication mediated by Notch signaling with trans-activation and cis-inhibition. The roles of Notch signaling in regulating fate decisions of endocrine and exocrine cells during the differentiation of pancreatic cells are explored. The results indicate that high (or low) level of Notch signaling drive cells to select the fate of exocrine (or endocrine) progenitor cells. The networks and the models presented here might be good candidates for providing qualitative mechanisms of pancreatic cell fate decisions. These results can also provide some insight on choosing perturbation strategies for further experimental analysis.

The kinetics in mathematical models on segmentation clock genes in zebrafish

Somitogenesis is the process for the development of somites in vertebrate embryos. This process is timely regulated by synchronous oscillatory expression of the segmentation clock genes. Mathematical models expressed by delay equations or ODEs have been proposed to depict the kinetics of these genes in interacting cells. Through mathematical analysis, we investigate the parameter regimes for synchronous oscillations and oscillation-arrested in an ODE model and a model with transcriptional and translational delays, both with Michaelis-Menten type degradations. Comparisons between these regimes for the two models are made. The delay model has larger capacity to accommodate synchronous oscillations. Based on the analysis and numerical computations extended from the analysis, we explore how the periods and amplitudes of the oscillations vary with the degradation rates, synthesis rates, and coupling strength. For typical parameter values, the period and amplitude increase as some synthesis rate or the coupling strength increases in the ODE model. Such variational properties of oscillations depend also on the magnitudes of time delays in delay model. We also illustrate the difference between the dynamics in systems modeled with linear degradation and the ones in systems with Michaelis-Menten type reactions for the degradation. The chief concerns are the connections between the dynamics in these models and the mechanism for the segmentation clocks, and the pertinence of mathematical modeling on somitogenesis in zebrafish.

Modelling asymmetric somitogenesis: Deciphering the mechanisms behind species differences

Somitogenesis is one of the major hallmarks of bilateral symmetry in vertebrates. This symmetry is lost when retinoic acid (RA) signalling is inhibited, allowing the left-right determination pathway to influence somitogenesis. In all three studied vertebrate model species, zebrafish, chicken and mouse, the frequency of somite formation becomes asymmetric, with slower gene expression oscillations driving somitogenesis on the right side. Still, intriguingly, the resulting left-right asymmetric phenotypes differ significantly between these model species. While somitogenesis is generally considered as functionally equivalent among different vertebrates, substantial differences exist in the subset of oscillating genes between different vertebrate species. Variation also appears to exist in the way oscillations cease and somite boundaries become patterned. In addition, in absence of RA, the FGF8 gradient thought to constitute the determination wavefront becomes asymmetric in zebrafish and mouse, extending more anteriorly to the right, while remaining symmetric in chicken. Here we use a computational modelling approach to decipher the causes underlying species differences in asymmetric somitogenesis. Specifically, we investigate to what extent differences can be explained from observed differences in FGF asymmetry and whether differences in somite determination dynamics may also be involved. We demonstrate that a simple clock-and-wavefront model incorporating the observed left-right differences in somitogenesis frequency readily reproduces asymmetric somitogenesis in chicken. However, incorporating asymmetry in FGF signalling was insufficient to robustly reproduce mouse or zebrafish asymmetry phenotypes. In order to explain these phenoptypes we needed to extend the basic model, incorporating species-specific details of the somitogenesis determination mechanism. Our results thus demonstrate that a combination of differences in FGF dynamics and somite determination cause species differences in asymmetric somitogenesis. In addition,they highlight the power of using computational models as well as studying left-right asymmetry to obtain more insight in somitogenesis.

A multi-cell, multi-scale model of vertebrate segmentation and somite formation

Somitogenesis, the formation of the body's primary segmental structure common to all vertebrate development, requires coordination between biological mechanisms at several scales. Explaining how these mechanisms interact across scales and how events are coordinated in space and time is necessary for a complete understanding of somitogenesis and its evolutionary flexibility. So far, mechanisms of somitogenesis have been studied independently. To test the consistency, integrability and combined explanatory power of current prevailing hypotheses, we built an integrated clock-and-wavefront model including submodels of the intracellular segmentation clock, intercellular segmentation-clock coupling via Delta/Notch signaling, an FGF8 determination front, delayed differentiation, clock-wavefront readout, and differential-cell-cell-adhesion-driven cell sorting. We identify inconsistencies between existing submodels and gaps in the current understanding of somitogenesis mechanisms, and propose novel submodels and extensions of existing submodels where necessary. For reasonable initial conditions, 2D simulations of our model robustly generate spatially and temporally regular somites, realistic dynamic morphologies and spontaneous emergence of anterior-traveling stripes of Lfng. We show that these traveling stripes are pseudo-waves rather than true propagating waves. Our model is flexible enough to generate interspecies-like variation in somite size in response to changes in the PSM growth rate and segmentation-clock period, and in the number and width of Lfng stripes in response to changes in the PSM growth rate, segmentation-clock period and PSM length.

Measurement of phase synchrony of coupled segmentation clocks

The temporal behavior of segmentation clock oscillations shows phase synchrony via mean field like coupling of delta protein restricting to nearest neighbors only, in a configuration of cells arranged in a regular three dimensional array. We found the increase of amplitudes of oscillating dynamical variables of the cells as the activation rate of delta-notch signaling is increased, however, the frequencies of oscillations are decreased correspondingly. Our results show the phase transition from desynchronized to synchronized behavior by identifying three regimes, namely, desynchronized, transition and synchronized regimes supported by various qualitative and quantitative measurements.

Intercellular coupling regulates the period of the segmentation clock

BACKGROUND

Coupled biological oscillators can tick with the same period. How this collective period is established is a key question in understanding biological clocks. We explore this question in the segmentation clock, a population of coupled cellular oscillators in the vertebrate embryo that sets the rhythm of somitogenesis, the morphological segmentation of the body axis. The oscillating cells of the zebrafish segmentation clock are thought to possess noisy autonomous periods, which are synchronized by intercellular coupling through the Delta-Notch pathway. Here we ask whether Delta-Notch coupling additionally influences the collective period of the segmentation clock.

RESULTS

Using multiple-embryo time-lapse microscopy, we show that disruption of Delta-Notch intercellular coupling increases the period of zebrafish somitogenesis. Embryonic segment length and the spatial wavelength of oscillating gene expression also increase correspondingly, indicating an increase in the segmentation clock's period. Using a theory based on phase oscillators in which the collective period self-organizes because of time delays in coupling, we estimate the cell-autonomous period, the coupling strength, and the coupling delay from our data. Further supporting the role of coupling delays in the clock, we predict and experimentally confirm an instability resulting from decreased coupling delay time.

CONCLUSIONS

Synchronization of cells by Delta-Notch coupling regulates the collective period of the segmentation clock. Our identification of the first segmentation clock period mutants is a critical step toward a molecular understanding of temporal control in this system. We propose that collective control of period via delayed coupling may be a general feature of biological clocks.