Persistence, period and precision of autonomous cellular oscillators from the zebrafish segmentation clock

  1. Alexis B Webb
  2. Iván M Lengyel
  3. David J Jörg
  4. Guillaume Valentin
  5. Frank Jülicher
  6. Luis G Morelli
  7. Andrew C Oates  Is a corresponding author
  1. MRC-National Institute for Medical Research, United Kingdom
  2. Max Planck Institute of Molecular Cell Biology and Genetics, Germany
  3. CONICET, Argentina
  4. Max Planck Institute for the Physics of Complex Systems, Germany
  5. University College London, United Kingdom
3 figures, 2 videos and 1 additional file

Figures

Figure 1 with 9 supplements
Zebrafish segmentation clock cells oscillate autonomously in culture.

(A) Confocal section through the tailbud of a Looping zebrafish embryo in dorsal view where the dotted line indicates the anterior limit of tissue removed. Nuclei are shown in red and YFP expression …

https://doi.org/10.7554/eLife.08438.003
Figure 1—source data 1

Summary table of segmentation clock tissue and cellular oscillatory properties.

Summary of statistics of time series traces recorded and analyzed in vitro in tailbud explants or tailbud cells. Peaks were identified, and the period/amplitude of cycles was determined as described in Materials and methods. A maximum period is defined in the method at 140 min, approximately twice the mean. Serum only cells were from the same cell suspension as those that were then treated with Fgf8b for the experiments 280711 and 250112. Pooled data from N = 2 independent cultures, for a total of n = 52 serum only cells. Pooled data from N = 4 independent cultures, for a total of n = 149 Fgf8b treated cells. To culture fully isolated cells, a cell suspension was serially diluted in wells within a 96-well plate, producing wells with a single tailbud cell. These were also treated with Fgf8b. N = 5 independent 96-well plate experiments, with a total of n = 10 fully isolated cells in these experiments.

https://doi.org/10.7554/eLife.08438.004
Figure 1—source data 2

Summary table of low-density segmentation clock cell experiments.

Description of in vitro cultured tailbud cell population treated with Fgf8b (n = 547), using multiple donor embryos in each of 4 independent experimental replicates (N = 4), carried out on separate days. Across the 29 fields recorded, we observed cell divisions in both YFP-negative (30, 5% of total cells) and YFP-positive cells (13, 2% of total cells). We found a range in the number of cell divisions from 0 to 5 cells per field, with an average of 1.5 (±1 SD) divisions per field. The categories of disqualification list the first event in a recording that led to disqualification. For example, four divisions in YFP-positive cells occurred after the cell had been disqualified for another reason (movement in and out of field, touching another cell).

https://doi.org/10.7554/eLife.08438.005
Figure 1—source data 3

Time series data from low-density segmentation clock cells.

XLS file containing all time series data for each of the 147 low-density segmentation clock cells in the presence of Fgfb. The file contains 4 work-sheets corresponding to each of the 4 independent replicates and to the plots in Figure 1—figure supplement 5. In each sheet, each cell is described by 3 neighboring columns: average fluorescence, local background, and background subtracted signal. Cells are also listed by their field of view in the original microscopy files.

https://doi.org/10.7554/eLife.08438.006
Figure 1—figure supplement 1
Her1-YFP-expressing cells in the zebrafish tailbud.

A confocal section through the tailbud of a Her1-YFP and Histone 2A-mCherry expressing 8-somite stage Looping zebrafish embryo in both lateral and dorsal orientations. The approximate location of …

https://doi.org/10.7554/eLife.08438.007
Figure 1—figure supplement 2
Peak finding in time series to estimate period and amplitude.

(A) An example of peak finding from a representative low-density cell trace, showing the steps used to find and estimate inter-peak intervals and amplitude. Top: raw data (red line) and background …

https://doi.org/10.7554/eLife.08438.008
Figure 1—figure supplement 3
Persistent oscillations in explanted tailbud.

(A) Montage of brightfield and corresponding YFP images from representative explanted tailbud over ~7 hr recording. Brightfield image is a single z-plane, while YFP signal is shown as an average …

https://doi.org/10.7554/eLife.08438.009
Figure 1—figure supplement 4
Time series of low-density segmentation clock cells in serum-only culture.

(A) Individual tailbud cells from experiment 280711 in the presence of serum. Black time trace is the raw data after background subtraction, the background level is the black line, the smoothed …

https://doi.org/10.7554/eLife.08438.010
Figure 1—figure supplement 5
Time series of low-density segmentation clock cells.

Traces from each independent low-density culture replicate (serum + Fgf8b) are shown in separate panels (#070312: yellow, #012512: green, #251012: red, #280711: blue). Each raw trace (black) was …

https://doi.org/10.7554/eLife.08438.011
Figure 1—figure supplement 6
Characterization of Ntla and Tbx16 antibodies.

(A, D) Graphical representation of the full-length Ntla and Tbx16 proteins with exons depicted in different colors. Blue arrows show the predicted T-box domain from amino acid 35 to 212, and 31 to …

https://doi.org/10.7554/eLife.08438.012
Figure 1—figure supplement 7
Expression of tailbud markers in vivo and in low-density cultures of segmentation clock cells.

(A) z-stack projection showing the expression patterns of Ntla (green) and Tbx16 (red) protein in a 12-somite stage wild type embryo detected using immunohistochemistry with monoclonal antibodies …

https://doi.org/10.7554/eLife.08438.013
Figure 1—figure supplement 8
Analysis of low-density segmentation clock cell cultures.

(A) Histogram of the number of peaks observed in each cell in the presence of serum + Fgf8b (blue) compared to those from cells in the presence of serum alone (orange). (B) Histogram of periods from …

https://doi.org/10.7554/eLife.08438.014
Figure 1—figure supplement 9
Time series of fully isolated segmentation clock cells.

(Top row) An example 40x transmitted light field of a single, isolated tailbud cell in a 96-well plate well. Scale bar = 20 μm. The cell’s corresponding fluorescence time series at the right shows …

https://doi.org/10.7554/eLife.08438.015
Figure 2 with 3 supplements
Dispersed low-density cells show a variety of behaviors compatible with slow amplitude fluctuations.

(A) Representative background-subtracted traces displaying different oscillatory behaviors: persistent oscillations, oscillations that initiate, stop, or start and stop within the recording time of …

https://doi.org/10.7554/eLife.08438.018
Figure 2—figure supplement 1
Statistics of amplitude and period correlations.

(A) Definition of averages Ai* (blue dot) of a quantity A measured in consecutive cycles (red square), and the distance Di (green line) to the identity (grey line), employed in panels (B) and (C). (B

https://doi.org/10.7554/eLife.08438.019
Figure 2—figure supplement 2
Numerical simulations.

(A) Numerical simulation of equation (S30), see Supplementary ile 1, with μ = 1, b = 1, ω = 2π/78 min-1, q = 0, . Signal x(t) oscillates with a constant amplitude and a period T = 78 min. There is …

https://doi.org/10.7554/eLife.08438.020
Figure 2—figure supplement 3
Both additive noise and color noise are necessary for the theory to describe the observed fluctuations.

(A) Amplitude is not affected by additive noise. Data points show the median peak amplitude for 1000 stochastic simulations (S30), see Supplementary file 1, for σμ2=0. Error bars are the 68% confidence …

https://doi.org/10.7554/eLife.08438.021
Figure 3 with 1 supplement
Precision of persistent segmentation clock oscillators

(A) Quality factor workflow for time series analysis for an example persistent oscillator. Sub-panel 1: Background-subtracted intensity over time trace from a single tailbud cell (black) with phase …

https://doi.org/10.7554/eLife.08438.022
Figure 3—source data 1

Precision and period calculation for persistent segmentation clock oscillators.

Each set of panels shows, successively, the background-subtracted average YFP intensity levels over time from a single persistently oscillating cell in black; the cosine of the phase calculated from the wavelet transformation in blue; and the autocorrelation function in green. The dashed green curve shows the analytical fit of the autocorrelation. Both period and quality factor can be calculated from this procedure (see Supplementary file 1). This is the complete persistent cell data set, a sub-set of the low-density set, from which the plots of period andquality factor QP in Figure 3B and D are generated.

https://doi.org/10.7554/eLife.08438.023
Figure 3—source data 2

Precision and period calculation for the tissue-level segmentation clock in the zebrafish embryo.

As for data set supplement 3–1, each set of panels shows, successively, the background-subtracted average YFP intensity levels from a region of posterior PSM tissue in a Looping embryo in black; the cosine of the phase calculated from the wavelet transformation in blue; and the autocorrelation function in green. The dashed green curve shows the analytical fit of the autocorrelation. Both period and quality factor can be calculated from this procedure. The original intensity versus time data comes from Soroldoni et al. (2014). This is the complete dataset from time-lapse data of 24 embryos from which the plot of quality factor QEmbryo in Figure 3B is generated.

https://doi.org/10.7554/eLife.08438.024
Figure 3—source data 3

Precision and period calculation for persistent circadian clock oscillators.

As for data set supplement 3–1, each set of panels shows, successively, the background-subtracted intensity levels from a single persistently oscillating Per2-Lucifcerase-expressing fibroblast over time in black; the cosine of the phase calculated from the wavelet transformation in blue; and the autocorrelation function in green. The dashed green curve shows the analytical fit of the autocorrelation. Both period and quality factor can be calculated from this procedure. The original intensity versus time data comes from Leise et al. (2012). This is the complete fibroblast dataset from which the plot of quality factor QF in Figure 3D is generated.

https://doi.org/10.7554/eLife.08438.025
Figure 3—figure supplement 1
Quality factor value depends on length of time series.

Time series length is defined in terms of the number of cycles. The plot shows the quality factor from stochastic simulations for two parameter sets A and B that display QA = 4 and QB = 10 when the …

https://doi.org/10.7554/eLife.08438.026

Videos

Video 1
Low-density segmentation clock cells oscillate in vitro.

Field of view containing cell in Figure 1B–C, highlighted by the red arrow. This field contains 18 cells in total, with 9 expressing YFP. We observed 5 cell divisions, the highest number in any …

https://doi.org/10.7554/eLife.08438.016
Video 2
Isolated segmentation clock cells oscillate in vitro.

Field of view corresponding to cell in top row of Figure 1—figure supplement 9. Total duration = 6.2 hr; Time interval = 2.14 min, field size = 205 x 205 μm, Scale bar = 50 μm.

https://doi.org/10.7554/eLife.08438.017

Additional files

Supplementary file 1

Theory supplement including Stuart-Landau theory with slow fluctuations, quantification of noisy oscillations, and numerical methods.

https://doi.org/10.7554/eLife.08438.027

Download links