Impacts of the feedback loop between sense-antisense RNAs in regulating circadian rhythms.

Antisense transcripts are a unique group of non-coding RNAs that are transcribed from the opposite strand of a sense coding gene in an antisense orientation. Even though they do not encode a protein, these transcripts play a regulatory role in a variety of biological processes, including circadian rhythms. We and others found an antisense transcript, Per2AS , that is transcribed from the strand opposite the sense transcript Period2 ( Per2 ) and exhibits a rhythmic and antiphasic expression pattern compared to Per2 in mouse. By assuming that Per2AS and Per2 mutually repress each other, our previous mathematical model predicted that Per2AS regulates the robustness and the amplitude of circadian rhythms. In this study, we revised our previous model and developed a new mathematical model that mechanistically described the mutually repressive relationship between Per2 and Per2AS via transcriptional interference. We found that the simulation results are largely consistent with experimental observations including the counterintuitive ones that could not be fully explained by our previous model. These results indicate that our revised model serves as a foundation to build more detailed models in the future to better understand the impact of Per2AS-Per2 interaction in the mammalian circadian clock. Our mechanistic description of Per2AS-Per2 interaction can also be extended to other mathematical models that involve sense-antisense RNA pairs that mutually repress each other.

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.

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.

Cell size homeostasis under the circadian regulation of cell division in cyanobacteria.

Bacterial cells maintain their characteristic cell size over many generations. Several rod-shaped bacteria, such as Escherichia coli and the cyanobacteria Synechococcus elongatus, divide after adding a constant length to their length at birth. Through this division control known as the adder mechanism, perturbation in cell length due to physiological fluctuation decays over generations at a rate of 2-1 per cell division. However, previous experiments have shown that the circadian clock in cyanobacteria reduces cell division frequency at a specific time of day under constant light. This circadian gating should modulate the division control by the adder mechanism, but its significance remains unknown. Here we address how the circadian gating affects cell length, doubling time, and cell length stability in cyanobacteria by using mathematical models. We show that a cell subject to circadian gating grows for a long time, and gives birth to elongated daughter cells. These elongated daughter cells grow faster than the previous generation, as elongation speed is proportional to cell length and divide in a short time before the next gating. Hence, the distributions of doubling time and cell length become bimodal, as observed in experimental data. Interestingly, the average doubling time over the population of cells is independent of gating because the extension of doubling time by gating is compensated by its reduction in the subsequent generation. On the other hand, average cell length is increased by gating, suggesting that the circadian clock controls cell length. We then show that the decay rate of perturbation in cell length depends on the ratio of delay in division by the gating τG to the average doubling time τ0 as [Formula: see text] . We estimated τG≈2.5, τ0≈13.6 hours, and τG/τ0≈0.18 from experimental data, indicating that a long doubling time in cyanobacteria maintains the decay rate similar to that of the adder mechanism. Thus, our analysis suggests that the acquisition of the circadian clock during evolution did not impose a constraint on cell size homeostasis in cyanobacteria.

Complementary phase responses via functional differentiation of dual negative feedback loops.

Multiple feedback loops are often found in gene regulations for various cellular functions. In mammalian circadian clocks, oscillations of Period1 (Per1) and Period2 (Per2) expression are caused by interacting negative feedback loops (NFLs) whose protein products with similar molecular functions repress each other. However, Per1 expression peaks earlier than Per2 in the pacemaker tissue, raising the question of whether the peak time difference reflects their different dynamical functions. Here, we address this question by analyzing phase responses of the circadian clock caused by light-induced transcription of both Per1 and Per2 mRNAs. Through mathematical analyses of dual NFLs, we show that phase advance is mainly driven by light inputs to the repressor with an earlier expression peak as Per1, whereas phase delay is driven by the other repressor with a later peak as Per2. Due to the complementary contributions to phase responses, the ratio of light-induced transcription rates between Per1 and Per2 determines the magnitude and direction of phase shifts at each time of day. Specifically, stronger Per1 light induction than Per2 results in a phase response curve (PRC) with a larger phase advance zone than delay zone as observed in rats and hamsters, whereas stronger Per2 induction causes a larger delay zone as observed in mice. Furthermore, the ratio of light-induced transcription rates required for entrainment is determined by the relation between the circadian and light-dark periods. Namely, if the autonomous period of a circadian clock is longer than the light-dark period, a larger light-induced transcription rate of Per1 than Per2 is required for entrainment, and vice versa. In short, the time difference between Per1 and Per2 expression peaks can differentiate their dynamical functions. The resultant complementary contributions to phase responses can determine entrainability of the circadian clock to the light-dark cycle.

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.

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.

A saturated reaction in repressor synthesis creates a daytime dead zone in circadian clocks.

Negative feedback loops (NFLs) for circadian clocks include light-responsive reactions that allow the clocks to shift their phase depending on the timing of light signals. Phase response curves (PRCs) for light signals in various organisms include a time interval called a dead zone where light signals cause no phase shift during daytime. Although the importance of the dead zone for robust light entrainment is known, how the dead zone arises from the biochemical reactions in an NFL underlying circadian gene expression rhythms remains unclear. In addition, the observation that the light-responsive reactions in the NFL vary between organisms raises the question as to whether the mechanism for dead zone formation is common or distinct between different organisms. Here we reveal by mathematical modeling that the saturation of a biochemical reaction in repressor synthesis in an NFL is a common mechanism of daytime dead zone generation. If light signals increase the degradation of a repressor protein, as in Drosophila, the saturation of repressor mRNA transcription nullifies the effect of light signals, generating a dead zone. In contrast, if light signals induce the transcription of repressor mRNA, as in mammals, the saturation of repressor translation can generate a dead zone by cancelling the influence of excess amount of mRNA induced by light signals. Each of these saturated reactions is located next to the light-responsive reaction in the NFL, suggesting a design principle for daytime dead zone generation.

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.

Assessment of Instructions on Protection Against Food Contaminated with Radiocesium in Japan in 2011.

The Japan Ministry of Health, Labour and Welfare (MHLW) has published instructions for radiological protection against food after the Fukushima Daiichi nuclear power plant accident in 2011. Following the instructions, the export and consumption of food items identified as being contaminated were restricted for a certain period. We assessed the validity of the imposed restriction periods for two representative vegetables (spinach and cabbage) grown in Fukushima Prefecture from two perspectives: effectiveness for reducing dietary dose and economic efficiency. To assess effectiveness, we estimated the restriction period required to maintain consumers' dose below the guidance dose levels. To assess economic efficiency, we estimated the restriction period that maximizes the net benefit to taxpayers. All estimated restriction periods were shorter than the actual restriction periods imposed on spinach and cabbage from Fukushima in 2011, which indicates that the food restriction effectively maintained consumers' dietary dose below the guidance dose level, but in an economically inefficient manner. We also evaluated the response of the restriction period to the sample size for each weekly food safety test and the instructions for when to remove the restriction. Stringent MHLW instructions seemed to sufficiently reduce consumers' health risk even when the sample size for the weekly food safety test was small, but tended to increase the economic cost to taxpayers.

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.

Feedback loops interlocked at competitive binding sites amplify and facilitate genetic oscillations.

Positive and negative feedback loops are often present in regulatory networks for genetic oscillations. Relative time scales and integration of these feedback loops are key to robust oscillations in expression levels. Using examples from the circadian clock and synthetic genetic oscillators, we study positive and negative feedback loops interlocked at competitive binding sites. In the mammalian circadian clock, a key clock gene Bmal1 is regulated by the activator ROR and the repressor REV-ERB. Conversely, Bmal1 activates both of them, forming interlocked feedback loops. Previous experiments indicate that the activator and repressor compete for the same binding sites in the Bmal1 promoter. Transcription patterns predict that ROR peaks later than REV-ERB and, moreover, the peak phase difference between them is small. Using mathematical modeling we reveal an optimal ratio of dissociation constants of an activator and a repressor for the competitive binding sites to enhance the amplitude of Bmal1 oscillations. This optimal ratio arises only when the amplitude of the repressor is larger than that of the activator. Secondly, we reveal that the preference of binding sites for an activator and a repressor depends on their relative time scales. A previous study demonstrated that noncompetitive binding sites are preferable for synthetic genetic oscillators that comprise a fast activator and a slow repressor with a large time scale separation. Here we show that when their time scales are similar, competitive binding sites are more likely to generate oscillation than noncompetitive sites. In contrast, for a slow activator and a fast repressor with a small phase difference as in Bmal1 regulation, noncompetitive binding sites are advantageous for amplifying oscillations. Our results, therefore, predict that additional mechanisms are necessary to compensate the disadvantage of the Bmal1 promoter and further facilitate amplification under the regulation by ROR and REV-ERB.

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.

Correction: Object Segmentation and Ground Truth in 3D Embryonic Imaging.

[This corrects the article DOI: 10.1371/journal.pone.0150853.].

Object Segmentation and Ground Truth in 3D Embryonic Imaging.

Many questions in developmental biology depend on measuring the position and movement of individual cells within developing embryos. Yet, tools that provide this data are often challenged by high cell density and their accuracy is difficult to measure. Here, we present a three-step procedure to address this problem. Step one is a novel segmentation algorithm based on image derivatives that, in combination with selective post-processing, reliably and automatically segments cell nuclei from images of densely packed tissue. Step two is a quantitative validation using synthetic images to ascertain the efficiency of the algorithm with respect to signal-to-noise ratio and object density. Finally, we propose an original method to generate reliable and experimentally faithful ground truth datasets: Sparse-dense dual-labeled embryo chimeras are used to unambiguously measure segmentation errors within experimental data. Together, the three steps outlined here establish a robust, iterative procedure to fine-tune image analysis algorithms and microscopy settings associated with embryonic 3D image data sets.

Genetic oscillators in development.

In development, morphogenetic processes are strictly coordinated in time. Cells in a developing tissue would need mechanisms for time-keeping. One such time-keeping mechanism is to use oscillations of gene expression. Oscillatory gene expression can be generated by transcriptional/translational feedback loops, usually referred to as a genetic oscillator. In this review article, we discuss genetic oscillators in the presence of developmental processes such as cell division, cell movement and cell differentiation. We first introduce the gene regulatory network for generating a rhythm of gene expression. We then discuss how developmental processes influence genetic oscillators. Examples include vertebrate somitogenesis and neural progenitor cell differentiation, as well as the circadian clock for comparison. To understand the behaviors of genetic oscillators in development, it is necessary to consider both gene expression dynamics and cellular behaviors simultaneously. Theoretical modeling combined with live imaging at single-cell resolution will be a powerful tool to analyze genetic oscillators in development.

Sample size allocation for food item radiation monitoring and safety inspection.

The objective of this study is to identify a procedure for determining sample size allocation for food radiation inspections of more than one food item to minimize the potential risk to consumers of internal radiation exposure. We consider a simplified case of food radiation monitoring and safety inspection in which a risk manager is required to monitor two food items, milk and spinach, in a contaminated area. Three protocols for food radiation monitoring with different sample size allocations were assessed by simulating random sampling and inspections of milk and spinach in a conceptual monitoring site. Distributions of (131)I and radiocesium concentrations were determined in reference to (131)I and radiocesium concentrations detected in Fukushima prefecture, Japan, for March and April 2011. The results of the simulations suggested that a protocol that allocates sample size to milk and spinach based on the estimation of (131)I and radiocesium concentrations using the apparent decay rate constants sequentially calculated from past monitoring data can most effectively minimize the potential risks of internal radiation exposure.

Collective cell movement promotes synchronization of coupled genetic oscillators.

Collective cell movement is a crucial component of embryonic development. Intercellular interactions regulate collective cell movement by allowing cells to transfer information. A key question is how collective cell movement itself influences information flow produced in tissues by intercellular interactions. Here, we study the effect of collective cell movement on the synchronization of locally coupled genetic oscillators. This study is motivated by the segmentation clock in zebrafish somitogenesis, where short-range correlated movement of cells has been observed. We describe the segmentation clock tissue by a Voronoi diagram, cell movement by the force balance of self-propelled and repulsive forces between cells, the dynamics of the direction of self-propelled motion, and the synchronization of genetic oscillators by locally coupled phase oscillators. We find that movement with a correlation length of about 2 ∼ 3 cell diameters is optimal for the synchronization of coupled oscillators. Quantification of cell mixing reveals that this short-range correlation of cell movement allows cells to exchange neighbors most efficiently. Moreover, short-range correlated movement strongly destabilizes nonuniform spatial phase patterns, further promoting global synchronization. Our theoretical results suggest that collective cell movement may enhance the synchronization of the segmentation clock in zebrafish somitogenesis. More generally, collective cell movement may promote information flow in tissues by enhancing cell mixing and destabilizing spurious patterns.

Interplay between intercellular signaling and cell movement in development.

Cell movement and local intercellular signaling are crucial components of morphogenesis during animal development. Intercellular signaling regulates the collective movement of a cell population via direct cell-cell contact. Cell movement, conversely, can influence local intercellular signaling by rearranging neighboring cells. Here, we first discuss theoretical models that address how intercellular signaling regulates collective cell movement during development. Examples include neural crest cell migration, convergent extension, and cell movement during vertebrate axis elongation. Second, we review theoretical studies on how cell movement may affect intercellular signaling, using the segmentation clock in zebrafish as an example. We propose that interplay between cell movement and intercellular signaling must be considered when studying morphogenesis in embryonic development.