Scholarly article on topic 'Faster embryonic segmentation through elevated Delta-Notch signalling'

Faster embryonic segmentation through elevated Delta-Notch signalling Academic research paper on "Biological sciences"

0
0
Share paper
Academic journal
Nature Communications
OECD Field of science
Keywords
{""}

Academic research paper on topic "Faster embryonic segmentation through elevated Delta-Notch signalling"

COMMUNICATIONS

ARTICLE

Received 23 Nov 2015 | Accepted 6 May 2016 | Published 15 Jun 2016^^BDOl7ia1o38/ncomms1186lM OPEN

Faster embryonic segmentation through elevated Delta-Notch signalling

Bo-Kai Liao1'2, David J. Jorg3'1 & Andrew C. Oates1'2'4

An important step in understanding biological rhythms is the control of period. A multicellular, rhythmic patterning system termed the segmentation clock is thought to govern the sequential production of the vertebrate embryo's body segments, the somites. Several genetic loss-of-function conditions, including the Delta-Notch intercellular signalling mutants, result in slower segmentation. Here, we generate DeltaD transgenic zebrafish lines with a range of copy numbers and correspondingly increased signalling levels, and observe faster segmentation. The highest-expressing line shows an altered oscillating gene expression wave pattern and shortened segmentation period, producing embryos with more, shorter body segments. Our results reveal surprising differences in how Notch signalling strength is quantitatively interpreted in different organ systems, and suggest a role for intercellular communication in regulating the output period of the segmentation clock by altering its spatial pattern.

1 Max Planck Institute of Molecular Cell Biology and Genetics, Pfotenhauerstr. 108, Dresden 01037, Germany. 2 Francis Crick Institute, Mill Hill Laboratory, The Ridgeway, London NW7 1AA, UK. 3 Max Planck Institute for the Physics of Complex Systems, Nöthnitzer Str. 38, Dresden 01187, Germany. 4 Department of

Cell and Developmental Biology, University College London, Gower Street, London WC1E 6BT, UK. f Present addresses: Cavendish Laboratory, Department of Physics, University of Cambridge, JJ Thomson Avenue, Cambridge CB3 0HE, UK and The Wellcome Trust/Cancer Research UK Gurdon Institute, University of

Cambridge, Tennis Court Road, Cambridge CB2 1QN, UK. Correspondence and requests for materials should be addressed to A.C.O. (email:

Andrew.Oates@crick.ac.uk).

The temporal coordination of events is an important, yet poorly understood aspect of development. A valuable model to investigate developmental timing is the rhythmic and sequential segmentation of the embryonic body axis into a series of multicellular units that anticipate the adult form. In vertebrates, these are called somites and give rise to segmented anatomy of the adult, such as vertebra, ribs and the associated muscles and skin. The rhythmic formation of embryonic segments is thought to be governed by a multicellular system of oscillating gene expression, termed as the 'segmentation clock', present in the posterior unsegmented tissue, the tailbud and presomitic mesoderm (PSM)1. Waves of gene expression appear in the posterior of the tissue and travel to the anterior where their arrest prefigures the time and position of each newly forming segment. These waves emerge at the tissue level from the coordinated output of cellular genetic oscillators, and their wavelength decreases as they move anteriorly. The period with which new segments are formed depends on the time scales of the genetic oscillations, and recent findings in zebrafish show that the period of segmentation can also be influenced by features of the wave pattern2. Despite this emerging picture, the genetic mechanisms that control the period at these various levels are not well-understood.

Although loss-of-function mutants in several genes increase the period of segmentation3, manipulating the system to run faster has been difficult. Removing introns from the mouse Hes7 gene leads to faster segmentation of the anterior axis, but the altered segmentation clock is unstable and segmentation is defective posterior to the neck4. Faster and stable segmentation of the entire body axis has not been engineered.

A central genetic player in segmentation is the Delta-Notch intercellular signalling pathway. The pathway has important roles in patterning and differentiation of many organ systems across the animal kingdom, and was associated with the rhythms of segmentation because of the waves of gene expression of several pathway members in the PSM5. In the zebrafish, the Delta-Notch pathway has been implicated in the synchronization of neighbouring oscillating PSM cells, leading to a coherent rhythm at the tissue level6-12. There are two delta genes expressed in the zebrafish PSM. Although deltaC shows oscillatory waves of expression and deltaD does not8, the formation of heterodimeric complexes between DeltaD and DeltaC13 and their joint regulation by the ubiquitin ligase Mind bomb14,15 suggest that both ligands contribute to rhythmic signalling between neighbouring PSM cells. deltaD mutants, known as after eight (aei)16,17, form their first 10 segments ~20% slower than wildtype (WT) and the length of these segments is correspondingly larger; other loss-of-function mutations and pharmacological treatments affecting this pathway show a similar phenotype18. Early during this developmental interval the wave pattern has a longer wavelength in the anterior PSM in mind bomb mutants than in WT siblings. In all loss-of-function cases, the coherent wave pattern and segment boundary formation is finally disrupted by the loss of synchrony leading to segmentation defects in the posterior trunk and the tail8,10,11,18.

These results raise the possibility that an increase in Delta-Notch signalling may shorten the period of segmentation. However, previous over-expression studies using mRNA injection or ubiquitous heat-shock induction of Delta and Notch resulted in segmentation defects10,12,19-21. The lack of the endogenous spatiotemporal expression pattern in these over-expression approaches may have masked a potential activity on the period of segmentation. In this report, we ask whether an appropriately directed increase in Delta-Notch signalling causes a change in segmentation period. We approach this by engineering zebrafish with extra copies of a transgene containing the deltaD locus

within its full genomic regulatory region to drive the correct spatiotemporal expression pattern. We assay protein levels, neurogenesis and segmentation in a set of DeltaD transgenic lines with different copy number and find that neurogenesis is sensitive to small changes in DeltaD levels, whereas only the highest level of DeltaD expression results in a shorter segmentation period. This shorter segmentation period is explained by elevated coupling between oscillators leading to changes in the synchronized wave pattern. Thus, beyond its previously described role in synchronizing oscillators, these findings suggest that Delta-Notch signalling can also affect the wavelength of the tissue's spatial pattern and thereby alter its output period.

Results

Transgenic lines with elevated accurate DeltaD expression.

To obtain accurate spatiotemporal expression, we used BAC recombineering to replace the stop codon of the deltaD gene with venus yellow fluorescent protein (YFP) in a 14-Kbp piece of genomic DNA that encompassed the deltaD exons22 (Supplementary Fig. 1). The venus-YFP tag was inserted to facilitate cellular-level imaging of DeltaD protein in vivo, which will be reported elsewhere; here we focus on the phenotypic effects of this transgene. We generated several deltaD transgenic lines using I-SceI meganuclease-mediated transgenesis23 to get single insertion-site, variable copy number lines in the aei/deltaD mutant background.

We analyzed two lines that we named Dover and Damascus. The mRNA patterns from the transgenic loci are indistinguishable from WT deltaD expression (Fig. 1a-f). Staining intensity from the homozygous Dover was higher than the heterozygote (Fig. 1d,e). We measured the copy number of deltaD-venus transgenes by Q-PCR and found that Dover has 7 copies and Damascus has approximately 100 (Fig. 1g-i). We examined DeltaD protein levels by quantitation of whole-mount immuno-histochemistry using monoclonal DeltaD antibodies13. We were able to quantitatively distinguish between WT, heterozygote and homozygote aei/deltaD mutants (Fig. 1j, inset of k). We observed an increased level of DeltaD in the tailbud and PSM in the transgenic lines ranging from a twofold increase in a Dover heterozygote up to a 50-fold increase in the Damascus heterozygote (Fig. 1k). We conclude that we made transgenic zebrafish with a range up to 50-fold excess of DeltaD, and with apparently normal expression patterns.

Elevated DeltaD shortens segmentation period. Before we assessed the transgenes for an effect on the period of the segmentation clock, we first tested them for activity by their ability to rescue the fully penetrant aei/deltaD ~/ _ segmentation phenotype. Using a high-contrast stain for xirp2a mRNA to detect the segmental boundaries, both Dover and Damascus rescued (Fig. 2a,b). Damascus had a low frequency (<3%) of segment defects scattered along the axis. We counted the total number of segments made in the body (Fig. 2c). Although the total duration of segmentation was the same (Fig. 2e), Damascus made 7.6% more total-body segments (WT 34.1 ± 0.7 mean ± s.d., n = 34; Damascus 36.7 ±0.7, n = 25; t-test P<0.001), and made 8.2% more segments in the trunk than the WT (WT 17.0 ± 0.3, n = 34; Damascus 18.4 ±0.5, n = 25, t-test P< 0.001), indicating that the extra segments were not added to the posterior tail. The other transgenes had no significant effect (t-test P>0.05).

We measured the segmentation period directly by multiple embryo time-lapse microscopy (Fig. 2d). Damascus made its (trunk) somites faster than WT siblings, with segmentation

deltaD

dld copy # j

m-5b ii i—ni m to-

a aeT'-m b aei+/_ ■ c

ZF deltaD-venus, 14 Kb

+/- +/+

tg copy #

200 150 150 50 0

tg/+ tg/tg 0

. 50 -

^ 30 -

. S 20 -

y = 0.5063x + 0.0633 R2 = 0.9963

Damascus ,g/+ 49.5 ± 7.5

9.21 ± 1.6

t WT AB • 1.00 ± 0.2

Dover ,g/+ 1.96 ± 0.3

Copy number

Figure 1 | A range of transgenic lines with elevated spatiotemporally accurate DeltaD expression. (a-f) deltaD-venus recapitulates WT deltaD expression pattern. Whole-mount ISH, 9 somite-stage embryos lateral view with either deltaD or egfp anti-sense riboprobe in deltaD mutant after eight (aei)+ /_ (a), WT AB (b) or deltaD-venus transgenic lines (d-f). (c) Schematic representation of deltaD-venus construct for generating transgenic lines. See Supplementary Fig. 1. (a) Scale bar, 300 mm. (f-h) deltaD transgenic copy numbers determined by quantitative real-time PCR from genomic DNA. Copy number is insensitive to aei point mutation allele (g). Transgenic copies (h,i) are total deltaD copies minus 2 endogenous copies. Empty bars, heterozygous genotypes; black bars, homozygous genotypes. Data pooled from >3 independent experiments, mean±s.d. (j) DeltaD protein expression visualized and quantitated by immunostaining. (ja-jc) PSM of flat-mounted 9 somite-stage embryos, circle shows region used for intensity measurement. Scale bar, 50 mm. (k) DeltaD protein expression levels versus deltaD gene copy number. Expression level in aei~/_ was defined as 0. Grey dashed line, linear fit as the formula and R2 value at upper left corner. Inset shows magnification of data points at origin. Mean±s.d. (n>5).

0.47 ± 0.1

periods of 23.1 ±0.8 and 24.7 ± 0.6 min, respectively (wildtype n = 26; Damascus n = 12; t-test P< 0.001). The magnitude of the decrease in period, 6.5%, is consistent with the increase in trunk segment number and total-body segment numbers. The effect on period was lost by incubating Damascus embryos in N-[N-(3,5-difluorophenacetyl)-L-alanyl]-S-phenylglycine t-butyl ester (DAPT), which blocks the g-secretase cleavage of the Notch intracellular domain24. This suggests that an elevated Notch signalling level is responsible for the shorter segmentation period. We conclude that the Damascus transgenic line has a decreased segmentation period and makes more, shorter segments.

Altered gene expression patterns linked to period change.

Recent work suggests that the wave pattern of oscillating gene expression in the PSM contributes to the period of segmentation in a shortening tissue due to a Doppler effect2. This Doppler effect arises because, as the tissue shortens, the anterior end of the PSM where the new segments are generated moves into the waves of oscillating gene expression and hence acts like an observer moving towards a wave-emitting source. This movement locally decreases the observed period with which waves arrive at the anterior end relative to the period of genetic oscillations in the posterior PSM. Since each wave arriving at the anterior end corresponds to one newly forming segment, this Doppler effect mediates a decrease in the period of segment formation. The magnitude of the Doppler effect is determined by the velocity of tissue shortening and the wavelength at the anterior end2,25 (Supplementary Note 1).

To test the hypothesis that a change in the Doppler effect contributes to the change in segmentation period in Damascus, we measured tissue-length change and anterior wavelength of the gene expression pattern. The rate of extension of the body axis and rate of shortening of the PSM was not different in Damascus and WT (WT axial extension 1.13 ±0.05 mm min _ 1, Damascus 1.06 ±0.11 mm min "1, t-test P = 0.154; WT PSM shortening 0.80 ±0.07 mm min "Damascus 0.79 ±0.03 mm min " t-test P = 0.791; rate determined from the linear fit from 0 to 600 min) (Fig. 2f), ruling out a contribution from a change in the velocity of tissue shortening. In contrast, we observed a systematically shorter anterior wavelength in the wave pattern of endogenous herl expression in Damascus compared to WT across the developmental interval encompassing trunk segmentation (Fig. 3a,b; Supplementary Fig. 2). Features of the wave pattern were also changed more posteriorly in the PSM, as evidenced by an increase in the number of gene expression waves observed in Damascus (Fig. 3; Supplementary Fig. 2). Thus, elevated DeltaD expression creates a characteristic new wave pattern in Damascus.

To assess the magnitude of the Doppler effect, only the anterior-most stripe is relevant and from the contribution of this change in the wave pattern, we calculated the expected period in Damascus to be 23.3 ± 0.3, 5.7% shorter than WT (see Methods, Supplementary Note 2; Supplementary Table 1; Supplementary Fig. 3). This accounts for the majority of the observed 6.5% difference in segmentation period in this interval (Fig. 2d). Combined, these results fulfil the expectation of a 'segmentation clock' with a shorter output period in an otherwise normally growing embryo.

......34

% Defects

80 ■

60 ■

40 ■

20 ■

= aei +" = dovertg/+ = dovertg/tg Dover sib. Damascus tg/+ Damascus sib.

WT AB Dovertg/+ Dovert3/ts Damascustg/+

15 20 25

Somite number

H ■ H — f ■

10 15 20 25 30

Somite number

140% 130% 120% 110% 100% 90% 80%

/ ^ ? «r ^

^ 110% m

1,400 700

3 ii»*

t1 • ■ WT AB Damascus'91*

0 100 200 300 400 500 600

' i ......

0 100 200 300 400 500 600 Time (min)

Figure 2 | Elevated DeltaD speeds somitogenesis. (a) Representative embryos at 36 h post fertilization (h.p.f.) stained with segment boundary marker cb1045 (xirp2a). Rescue of segmentation defects in aei-/- genotype background by dover and damascus. Numbers above arrowheads indicate the posterior boundaries of the corresponding ordinal segments. Orange triangle, segment boundary aligned with proctodeum; purple triangle, last segment boundary observed; grey triangle, defective segment boundaries or onset of defects. Scale bar = 300 mm. (b) Distribution of somite defects showing deltaD-venus rescues aei -' - throughout somitogenesis. No defective boundaries observed in WT AB (n = 74), aei +/- (n = 76), Dovertg/ + (n = 95) and Dovertg/ tg (n = 60) (blue line). Damascus showed low penetrance of defects (2.6% from the 1st to 30th segments, n = 64) (green line). aei-/- siblings (sib) of dover (n = 29) or Damascus (n = 33) showed typical deltaD lost-of-function phenotypes, mean ALD around eight (red lines). (c) Damascus has more segments and shifted segmental position of proctodeum. Quantification from embryos as shown in a. Orange, level of proctodeum; magenta, total segment number; grey bars, extreme values. (d) Somitogenesis period of trunk region (s4-s19) measured at 28.3 ±0.3 °C by multiple embryo time-lapse imaging. Mean periods normalized to WT AB. Data pooled from >3 independent trials (n>12 embryos for each condition), except aei-/- siblings from 1 trial (n = 6). (e) Total segmentation time is the interval between the 4th segment to the last segment observed in time-lapse movies. Quantification of corresponding movies described in b. Data pooled from >3 independent trials (n>12 embryos for each condition). Mean durations normalized to WT AB (f) Axis elongation and PSM shortening are not altered in Damascus. Data analyzed from same time-lapse movies as d. Time 0 is defined at 5 ss. ***P<0.001, NS, not significant (P>0.05), Student's t-test. Error bars show the 5th and 95th percentiles of the data.

Why does elevated DeltaD expression lead to an altered wave pattern with a shorter anterior wavelength in Damascus? Elevated Delta expression would be anticipated to cause elevated Notch signalling, and this in turn would lead to an elevated coupling

between neighbouring oscillators8,11. Previous theoretical work has suggested that altered coupling can alter the spatial pattern and output period of a 'segmentation clock'18,26,27, and we return to examine this process in the last section of the results.

Interpreting the period phenotype requires a better understanding of the relationship between Delta expression levels and Notch signalling. We expect that in the various lines we have generated (Fig. 1j,k), increased levels of DeltaD protein expression give rise to correspondingly elevated levels of Delta-Notch coupling between neighbouring cells. However, the lack of an effect on the period in Dover raises the possibility that Notch

signalling was not elevated in this line. Alternatively, the Damascus allele may have a neomorphic activity. A third possibility is that the segmentation clock is strongly buffered against increases to Delta-Notch signalling. To investigate this further, we examined two indexes of Notch signalling in the embryo, lateral inhibition and synchronization, across the range of DeltaD expression levels.

Damascus tg/+

her1 antisense riboprobe

N= 17 21

Elevated DeltaD increases lateral inhibition in the CNS. Firstly, we tested whether the transgenic lines had an effect on primary neurogenesis. The role of Delta-Notch signalling in neurogenesis is better understood than in somitogenesis, and is thought to involve a process termed lateral inhibition that generates finegrained patterns in which neighbouring cells take on differing fates28-30. Elevated Delta levels, and thus increased Notch activation, yields a reduced number of cells selecting the primary neuronal fate in all species investigated31-34.

We examined the numbers of primary motor neurons (PMN), Rohan-Beard neurons (RB) and trigeminal neurons (TG) by isletl expression35. For each subtype, aei/deltaD homozygote embryos had more primary neurons than WT (Fig. 4a,g-i), as expected16,22,34. With higher copy number of transgenes, we found a graded decrease in the density of neurons for all categories (Fig. 4d-i). Primary motor neurons were the most sensitive to changes in deltaD copy number: the embryo distinguished between all copy number differences except between 1 and 2 (Fig. 4g). This shows that the Dover transgene produced elevated levels of active Delta and this translated into correspondingly elevated Notch signalling. It also indicates that Damascus acts as a strong hypermorphic allele of deltaD in the nervous system. These results suggest that, in contrast to Delta's effect on segmentation period, the lateral inhibition mechanism in the CNS converts the input level of Delta protein smoothly to the Notch signalling output.

Elevated DeltaD increases coupling between PSM oscillators.

Secondly, we tested whether Notch signalling was elevated in the segmenting tissue using a modified, quantitative re-synchronization assay9,11 (Fig. 5a). In this assay, Notch signalling is blocked before the onset of oscillations in the early embryo by incubation in saturating DAPT, which prevents proteolytic cleavage of the intracellular domain. Over time, the population of cellular oscillators desynchronizes at a rate set by the noise level of the individual oscillators, until falling below a critical level where penetrant segmentation boundary defects are observed posterior

Figure 3 | Damascus has shorter anterior wavelength. (a) Measurement of herí anterior wavelength from 4 to 20 somite-stage. In each representative embryo panel, left part: examples of WT TL or Damascustg/ + embryos in lateral view, anterior to left, after ISH with herí probe. The line of interest for image processing is visible in yellow, extending from the tip of the tail bud to the somites. Right part: (inverted) gene expression signal intensity (x axis) versus position along the anterior-posterior axis (y axis) plots for the embryos in the left-hand part. The tick marker on the y axis indicates the X0 position, and y axis is the scale bar (700 mm) to its representative embryo as well. Peak maxima is indicated with red triangles. The anterior wavelength, defined as the distance between the most anterior two gene expression maxima, is indicated by a red bracket. A, anterior; Int., intensity in arbitrary unit; P, posterior. (b) The most anterior wavelength of herí mRNA gene expression pattern. The average ratio (in percentage) between WT TL or Damascustg/ + embryos for each stage examined are labelled above the plots. Blue, WT; Green, Damascustg/ +. The central boxes of box-and-whisker plots cover the interquartile range with median as line within box. Whiskers are 5th and 95th percentiles, extreme values are grey bars. *P<0.05, **P<0.01, ***P< 0.001, Student's t-test.

/ ^ o^ J> J * 0° 0°

Figure 4 | Elevated DeltaD increases the strength of lateral inhibition in the CNS. (a-f) Representative embryos stained with isletl, an early marker of primary neuronal differentiation. In each box, 4 somite-stage embryos lateral view, anterior to left (upper panel); dorsal view with anterior to top (lower left panel), of head showing bilateral trigeminal ganglia (TG); mid-trunk showing PMN at the midline, and RB neurons at lateral margin of neural plate (lower right panel). (a) Scale bars, 250 mm. (g) Cell density of PMN. (h) Cell density of RB sensory neurons. (i) Area of the TG. Lower densities or smaller areas indicate higher Notch signalling strength. Data from aei~7~ (n = 33), aei+7" (n = 18), WT (n = 35), Dovertg/ + (n = 14), Dovertg/tg (n = 12) and Damascustg/ + (n = 22) pooled from >2 independent trials. The central boxes of the box-and-whisker plots cover the interquartile range with median as line within box. Whiskers are 5th and 95th percentiles, extreme values are small bars. *P<0.05, **P<0.01, ***P< 0.001, NS, not significant (P>0.05), Student's t-test.

to approximately segment number 5, the anterior limit of defects (ALD)8'11. When DAPT is washed out, Notch signalling begins to couple the neighbouring oscillators, and the population gradually resynchronizes. The position, and hence time, of the first segment boundary to form correctly again is recorded as a proxy for the synchronization state of the population having risen beyond the critical level. We term this the first recovered segment (FRS). The time taken to re-synchronize is expected to depend on the coupling strength11, and faster re-synchronization, that is, lower FRS, is expected from a higher Delta-Notch coupling strength.

We asked whether re-synchronization was sensitive to the difference between one and two functional copies of deltaD, and observed a higher FRS in aei/deltaD heterozygotes than WT (WT 19.1 ±2.2, n = 94; aei+1 ~ 21.3 ±2.0, n = 48; t-test Po 0.001), indicating a 20.5% slower re-synchronization time (Fig. 5b,c). We measured the re-synchronization time in the transgenic lines and found that Damascus had a lower FRS than WT and Dover (Damascus 17.8 ± 1.9, n = 78; t-test to WT P = 0.002), indicating that Damascus rescued segment boundaries 15.0% faster (Fig. 5b,c). Although we cannot rule out subtle changes to the intracellular oscillator circuit, for example, a potential increase in

the noise as evidenced by a decrease in ALD (Fig. 5c), this faster re-synchronization can be explained by a higher coupling strength in Damascus, consistent with elevated Notch signalling between oscillators in the population. Because the re-synchronization time did not change in Dover, these results suggest that the synchronization mechanism is strongly buffered across a range of Notch signalling.

Wave patterns in a simplified model of coupling with delays.

In systems of coupled oscillators, time delays in the coupling can give rise to a plethora of complex phenomena. In considering these effects, it is important to distinguish between the strength of coupling and its delay. Coupling strength can be thought of as the number of activated ligand-receptor pairs over time, which is expected to increase as a consequence of the over-expression of DeltaD. In contrast, the time taken to transfer this signal from one cell to another is termed the coupling delay. Because of the time required to synthesize and traffic Delta proteins, a non-negligible time delay in the coupling is expected for oscillators coupled by the Delta-Notch signalling13,18,36. Direct measurement of coupling strength and time delays is not

currently possible in any embryo, although values for WT and mutant zebrafish have been inferred18. To investigate whether an internally consistent description of the spatial features of the wave patterns observed in Damascus can be found in the presence of elevated coupling strength and realistic time delays, we turned to a simplified theoretical description, which ignores the effects of tissue shortening26.

We numerically calculated wave patterns of gene expression for the zebrafish system using a range of coupling strengths and time delays (Supplementary Note 3). We observed wave patterns with features corresponding to the WT embryo, as expected, as well as

wave patterns corresponding to Damascus. These latter resulted from elevated coupling strength and a plausible alteration of the coupling delay (Supplementary Fig. 4). Supporting the robustness of this altered state to parameter changes, we found multiple regions in parameter space in which elevated coupling strength gives rise to a stable time-periodic wave pattern with an increase in the number of waves and a shortened anterior wavelength (Supplementary Fig. 5) as observed experimentally in Damascus. These theoretical results provide an internally consistent description that qualitatively describes the differences between the phenotypes of WT and Damascus embryos, supporting the interpretation of an elevated coupling strength in Damascus.

w Zmin

+ DAPT

No DAPT

I I I I I I I I I

I I I I

5 h.p.f.

10 h.p.f.

15 h.p.f.

I I I I

20 h.p.f.

TDe-sync - ALD t0 TRe-sync - FRS - tW ~ e

I-► I-—

aei WT AB Dovertg/+ Dovertg/tg Damascustg/

Somite number

Discussion

In this report, we describe the first embryo genetically engineered to tick more quickly throughout development. This was achieved by the generation of a zebrafish line carrying multiple copies of a YFP-tagged deltaD transgene expressed from within its endogenous flanking genomic DNA. Direct measurements of a shorter segmentation period and a corresponding increase in the total number of segments in an otherwise normally-developing embryonic axis argue strongly for a 'segmentation clock' with a short output period.

A change in segmentation period in an animal with altered Delta/Notch signalling could potentially result from the effects of coupling on the overall period of a population of oscillators, if the coupling has time delays on the order of the period, as is thought to be the case in vertebrate embryos27. A lengthening of period due to this process has been inferred for embryos with a reduction in Delta-Notch signalling in the zebrafish18. In contrast, mouse segmentation appears to get faster with pharmacological blockade of Notch37, an outcome predicted if the length of the time delay in the coupling was shorter than half the period27. These descriptions of the system, as used in Supplementary Figs 4 and 5 to check for the effect of an elevated coupling strength, are based on a simplified scenario with a time-periodic wave pattern and constant tissue length.

Recent findings of the effects of waves of gene expression and shortening tissue on segmentation period raise the possibility that changes to the normal pattern of the waves could alter the segmentation period2. Our results show that Damascus embryos have a shorter gene expression wavelength in the anterior PSM, and the magnitude of this change accounts for most of the

Figure 5 | Elevated DeltaD increases the strength of coupling in the segmentation clock. (a) Schematic of DAPT pulse-chase synchrony rescue assay for estimating relative coupling strength. Presumptive synchrony level (Z) versus experimental time scale shown in upper panel. Z is an order parameter that measures the degree of synchrony between neighbouring oscillators, with 1 indicating perfect synchrony and 0 indicating complete lack of temporal correlation, as defined in ref. 11. Z0, initial synchrony level at shield stage (6 h.p.f.); Zc, critical level of synchrony for forming non-defective segment boundaries; Zmin, putative minimal synchrony at de-synchronized state; darkness of blue lines depict relative coupling strength. Time intervals (TRe-sync) between DAPT wash-out and resynchronized states reflect coupling strength (e). (b) Segment boundaries after DAPT pulse-chase in mutant and transgenic backgrounds. Numbers above arrowheads indicate the values of ALD (red triangle) or FRS (blue triangle). Scale bar, 300 mm. (c) Quantification of DAPT pulse-chase experiments. ALD (red) and FRS (blue) from left or right side independently, as shown in a. Data from aei+/_ (n = 48), WT (n = 94), Dovertg/ + (n = 43), Dovertg/tg (n = 36) and Damascustg/ + (n = 78) pooled from three independent trials. Data is mean value with error bars showing 5th and 95th percentiles. *P<0.05, ***P< 0.001, NS, not significant (P>0.05), Student's t-test.

observed change in segmentation period through a change in the contribution of the Doppler effect. Regulation of the wave patterns has been proposed to result from gradients of signalling that span the PSM, which are generated by elevated levels of Fgf and Wnt in the posterior27,38-42, but theoretical work shows that changes in the coupling could also alter the wave patterns18,26,27,43.

Our observations of elevated Delta-Notch signalling throughout the embryo, including a change in re-synchronization rate, imply that a change in cell-cell communication is responsible for the shorter anterior wavelength. These effects seem likely to be caused by an elevation in Delta-Notch signalling strength between neighbouring oscillators, but we cannot rule out a contribution from changes to signalling delays, in addition. Combining these lines of evidence, this finding supports the prediction that the changed wave pattern results in a change in the magnitude of the Doppler effect. Since one arriving wave at the anterior end is related to one newly formed segment, this wave effect corresponds to a shortening of the period of segment formation.

Direct observation of the genetic oscillations is not currently possible in Damascus because available reporters of the

segmentation clock in zebrafish express Venus-YFP2,6, as does Damascus. However, future experiments will investigate whether the observed changes to the anterior wavelength are accompanied by changes to the period of genetic oscillations in the tailbud, or longer-term changes in the wave pattern2,25. This will provide a fuller picture of the effects of elevated DeltaD signalling in zebrafish segmentation. Our results reinforce the idea that the period of somitogenesis is influenced by several different processes: by single-cell Hes/Her gene feedback circuits4,44,45, multicellular local Delta-Notch signalling18,37, and by the global patterns of cyclic gene expression waves and changing tissue length2,25.

Delta-Notch signalling is essential for multiple events in development, coordinating patterning and differentiation in space and time in various organ systems. It is not well-understood how similar signals are used to different effect in different tissues, for example, synchronization versus lateral inhibition, nor how such differences evolve. We observed a striking difference in how a change in gene copy number was interpreted in segmentation compared to neurogenesis. Segmentation is sensitive to changes in coupling strength in the range from 0 to 2 copies of WT Delta

Output V Segment intactness

Output II Lateral inhibition strength

Output I Expression level

Gene copies

Figure 6 | Robustness and flexibility of Delta-Notch signalling. Five experimentally measured system outputs (Assay, left-had axes) and their corresponding biological meaning (Implication, right-hand) are displayed as functions of deltaD copy numbers, shown on x axis in log2 scale, the 0.5 position is defined as 0 copy for aei~^~. DeltaD protein levels in posterior PSM from Fig. 1k shown in log2 scale, the 0.125 position on the axis defined as 0. PMN intervals from Fig. 4g. Re-synchronization (re-sync) rate calculated from the reciprocal of re-synchronization time from DAPT pulse-chase experiments (Fig. 5c), value of aei~^~ is defined as 0. Somitogenesis frequency from Fig. 2d. Proportion of non-defective somites from Fig. 2b. Lines are connections of data points except for DeltaD protein levels, which is linear fit. The five outputs were ordered from bottom to top based on presumptive increasing robustness of the system. The light green background marks the robust zone, where copy number does not affect system output, and the outer areas reflect the flexibility of output based on the change of WTand transgenic deltaD copy numbers. Data is mean value with error bars showing 5th and 95th percentiles. Green lines connect data points that are not significant (P>0.05) to WT in Student's t-test.

or Notch genes (this study and refs 9,11), but its properties appear to be well buffered for increased Notch signalling strength until much higher levels are attained (Fig. 6). In contrast, neurogenesis appears to be sensitive to changes in Notch signalling strength across the full range of endogenous and transgenic copy numbers in this study.

Given the frequency of gene duplication observed in metazoan lineages46,47, the range of copy numbers we assayed may mimic a relatively common mutational scenario. From our findings, we speculate that following the repeated duplication of a Delta gene locus, an evolutionary alteration of the nervous system might be possible without alteration of the segmented body plan. A difference in the sensitivity of developmental processes to gene copy number may enable the dissociation of evolutionary trajectories of different organ systems. Whether this would be of advantage to the mutant individual would of course depend on local selective pressures.

Do differences in levels of Notch signalling components, for example, Notch receptors, underlie the observed differences between the responses of neurogenesis and segmentation to elevated Delta? From our data, it does not appear that signalling component levels are limiting, since there is no evidence of saturation in either embryonic context; both assays show changes between Dover and Damascus (Fig. 6). Rather it appears that the tissues show a difference in the gain of their outputs across the range of Delta input. These differences may be due to tissue-specific signal transduction components that confer distinct properties of information processing in lateral inhibition and synchronization mechanisms. Answering these questions will allow multicellular biological systems to be better understood in development and disease states, and also to be engineered with a greater range of spatial and temporal characteristics in the future.

Methods

Animals and embryos. Zebrafish (Danio rerio) adult stocks in the fish facilities of the Max Planck Institute of Molecular Cell Biology and Genetics and University College London were kept in 28 0C fresh water under a 14:10 h light:dark photo-period. Embryos were collected within 30 min following fertilization and incubated in petri dishes with E3 media. Until the desired developmental stages48, embryos were incubated at 28.5 0C. For whole-mount in situ hybridization (ISH), PTU (1-phenyl 2-thiourea) at a final concentration of 0.003% was added before 12h.p.f. to prevent melanogenesis. All WTs were AB, unless otherwise specified, and the deltaD mutant allele used was aeitr233 16 crossed onto the same AB background. Zebrafish experimentation was carried out in strict accordance with the ethics and regulations of the Saxonian Ministry of the Environment and Agriculture in Germany under licence Az. 74-9165.40-9-2001, and the Home Office in the United Kingdom under project licence PPL No. 70/7675.

DAPT treatment and pulse-chase experiments. DAPT treatment was carried out as previously described11. In brief, 50 mM DAPT stock solution (Merck) was prepared in 100% DMSO (Sigma) and stored in a small volume at — 20 0C. Embryos in their chorions were transferred to 12-well plates at 2h.p.f. in 1.6 ml E3 medium with 20 embryos per well. A total of 50 mM DAPT in E3 medium was prepared immediately before the treatment. To prevent precipitation, the DAPT stock solution was added into E3 medium while vortexing, and then filtered by 0.22 mm PES syringe filter (Millipore). DAPT treatment was initiated at dome stage (4 h.p.f.) by replacing E3 medium with E3/DAPT medium. For pulse-chase experiments, DAPT was washed out at least twice with fresh E3 medium + 0.03% PTU. Embryos were dechorionated and fixed at 36 h.p.f. All experiments were incubated at 28.5 0C, except for short operations, for example, washing out, which were at room temperature.

Scoring segment defects and segment number. Embryos stained with the cb1045 riboprobe marking mRNA from the xirp2a gene were visually scored under an Olympus SZX-12 stereomicroscope and images were acquired with a QImaging Micropublisher 5.0 RTV camera. Defective segmentation boundaries were scored following11 with modification. In short, any abnormal boundary was defined as defective, and two measures were collected for each left-right side of the embryo: an ALD, that is, the position of the first defective boundary and the FRS, that is, the anterior boundary position of the first non-defective segment after the DAPT wash-out at 9 ss. In all scoring, the first two somites and those posterior to the 30th were defined as non-defective since the first two were clutch-dependent and the

posterior-most segments were too small for reliable analysis. Total segment number was counted in 36-39 h.p.f. embryos; embryos earlier than this stage had weaker xirp2a staining patterns. The re-synchronization rate is defined as the multiplicative inverse of the boundary difference between DAPT wash-out and FRS, that is, FRS minus 9.

BAC recombineering. Red/ET recombination was carried out as published (Supplementary Fig. 1, upper panel)49-51. Briefly, the recombination was performed in an E. coli HS996 strain carrying a toolbox plasmid, pRed-flp4. pRed-flp4 has a thermo-sensitive replication origin, pSC101, and contains anhydrotetracycline inducible Flippase (Flp) and an operon of homologous recombinases, Reda, Redft, Redg and RecA, controlled by a L-rhamnose-inducible promoter. The venus-kanamycin tagging cassette was generated by PCR from the plasmid template pR6K-venus in which the R6K origin provides background-free tagging because its replication cannot proceeded without a specific initiator52.

A commercially available BAC clone encompassing the deltaD locus (ID: CH1073-241D9) was first electroporated into HS996/pRed-flp4, then confirmed by colony PCR of overnight cultured bacterial colonies on LB-agar solid media with corresponding antibiotics. The positive clone of HS996/pRed-flp4/deltaD-BAC was cultured overnight in LB liquid media with antibiotics, and then 60 ml of the saturated culture was transferred to 1.8 ml fresh medium at 30 0C on Thermomixer shaker (Eppendorf) until its OD600 reached 0.2. Red expression was induced by adding 30 ml of 25% L-rhamnose to the liquid culture for 1 h at 37 0C. After the OD600 reached 0.4, 500 ng of the tagging cassette was eletroporated into HS996/pRed-flp4/deltaD-BAC, and then the culture was transferred into SOC medium without antibiotic for 1 h at 37 0C. After re-growth, the culture was centrifuged and transferred to LB medium or an LB-agar plate with 10mgml— 1 kanamycin for overnight culture. The positive clones were screening by colony PCR for the region across the deltaD stop codon, and the kanamycin flanked by two FRT sites was then removed by the induction of Flp with 200 nM anhydrotetracycline in the medium.

To subclone, the vector backbone for recombination was amplified by PCR from a plasmid template, pNPC2, which is a p15A origin plasmid containing two I-SceI recognition sites, a ClonNAT resistant gene, and ccdB toxic gene for reducing recombination background. Homology arms targeting the 5' upstream or the 3' downstream deltaD locus resulting in a 14-Kbp fragment were designed to the 5' ends of the subcloning primers. Recombination is similar to the tagging described above, with the exception that the final culture was plated on LB-Agar plates with 5 mgml — 1 ClonNAT (WERNER BioAgents). Positive clones were screened by colony PCR and plasmid DNA was prepared using a Qiagen Midi-Prep kit according to the manufacturer's instructions. The final construct, pNPC2-deltaD-venus, was verified by sequencing the tagging region and the conjunctions of the vector backbone to the subcloned fragment (MPI-CBG sequencing facility). The list of primer sequences and PCR conditions can be found as Supplementary Table 2.

Microinjection and I-SceI-mediated transgenesis. Meganuclease transgenesis was performed as the previously described (Supplementary Fig. 1, lower panel)23. In brief, injection needles were pulled in preparation from 1.0 mm O.D. x 0.58 mm I.D capillaries (Harvard Apparatus) using a Flaming/Brown micropipette puller (Sutter Instruments), and meganuclease I-SceI (Roche) was stored in — 80 0C as 2 ml aliquots.

Staged embryos from an incross of aei/deltaD homozygous adults were placed into a 1.5% agarose (Sigma) injection mould in 0.3 x Danieau's buffer and injected at the one cell stage. A 20 ml injection mixture was made immediately before injection from 14 ml of pNPC2-deltaD-venus (100 ng ml— 1), 2 ml of 10 x I-SceI buffer, 2 ml of 50 mM MgCl2, and 2 ml I-SceI (10 U ml— 1). A pneumatic Pico Pump PV 820 (World Precision Instruments) was used for injection. After loading the injection mixture into the needle with a microloader tip (Eppendorf), the injection volume was adjusted to ~ 1 nl by measuring the diameter of bolus injected into a drop of mineral oil (Sigma) on the objective micrometre. The injection mixture was directly injected into the cytoplasm of the embryos. After injection, embryos were transferred into E3 medium and screened for positive Venus-YFP fluorescence around 24 h.p.f. I-SceI meganuclease-mediated transgenesis is known for its high transgenesis efficiency and resulting in variation of transgene copy number with few genome insertions53-55. The potential founders were screened by whole-mount ISH staining of transgene expression and the rescue of their aei— / — background phenotype.

Genomic DNA extraction and quantitative real-time PCR. Genomic DNA was extracted from either fin-clips of adults or twenty 5-days post fertilization embryos using a DNeasy Blood & Tissue Kit (Qiagen, Cat No./ID 69504) with RNaseA treatment as described in the manufacturer's instructions. The concentrations of genomic DNA extracts were then measured by Nanodrop 2000 (Thermo scientific) and adjusted to 6ng ml— 1. Quantitative real-time PCR was carried out using the Mx3000 Real-time PCR system (Stratagene) with its built-in acquisition/analysis software. To prepare the PCR, two replicates of each sample were loaded into 96-well plates with a final volume of 10 ml in each well, containing 5 ml 2 x SYBR Green master mix (ThermoFisher Scientific), 30 nM ROX free dye (ThermoFisher Scientific), 70 nM primer pair, and 3 ml genomic DNA. The PCR programme was 95 OC 15 min for initiation, 40 cycles of 95 OC 30 s, 60 OC 30 s and 72 OC 30 s,

followed by the default dissociation analysis. Seven pairs of primers targeting the deltaD genome locus, deltaC genome locus or venus open reading frame were used (Supplementary Table 3), and the standard curve and PCR efficiency of each primer pair was determined by serial dilution of WT samples.

To analyze the samples, we first normalized the raw Ct. values from deltaC or deltaD primer pairs to the in-group WT values according to the PCR efficiency of each primer pair (Supplementary Table 3). The relative concentrations within each primer group were then divided by the averaged deltaC concentrations to calibrate the loading error. In the case of venus primer pairs, the raw Ct. values were normalized to one assigned transgenic sample. Thus, the copy numbers of each sample were averaged from the relative deltaD concentrations to WT samples and, if applicable, from the relative venus concentrations to each other. The copy numbers of dover and damascus were obtained from the average of at least three independent trials.

Immunohistochemistry and quantification protein level. Dechorionated embryos were fixed in 4% paraformaldehyde/PBS for ISH. For whole-mount immunohistochemistry, 9 somite-stage embryos were washed several times with PBS/0.1% Tween and then permeabilized by washing with PBS/0.2% Triton X-100 for 30 min at room temperature (RT). Blocking solution (Roche) prepared in PBS/0.2% Triton X-100 was used to block at RT for 2h, then the embryos were incubated with 1:100 PBS/0.2% Triton X-100 diluted mouse anti-DeltaD monoclonal antibody (Abcam) at 4 0C overnight. After primary antibody incubation, embryos were washed with PBS/0.2% Triton X-100 four times for 30 min, then incubated with 1:200 PBS/0.2% Triton X-100 diluted goat anti-mouse IgG conjugated with Alexa Fluor 488 (Invitrogen) for 2 h at RT. Secondary antibody was washed out with PBS/0.2% Triton X-100 at least four times for 30 min at RT. After staining, samples were gradually transferred into 87% glycerol/13% PBS and kept at 4 0C overnight before flat-mounting. ProLong Gold Antifade Reagent (Invirogen) was used as a mounting reagent.

To compare genotypes, embryos were stained in parallel in each experiment. Images from PSM of flat-mounted samples were acquired by a Zeiss LSM 510 META confocal laser-scanning microscope. Acquisition configurations were identical for all samples. Mean fluorescence intensity of the PSM was obtained from z-projection of three non-overlapping z-sections with 15 mm intervals, and the middle section was aligned to the centre of the notochord. The relative DeltaD levels in the posterior PSM were the average intensities obtained from circular areas of 68 mm diameter.

Time-lapse imaging and analysis. Somitogenesis period was determined by multiple embryo time-lapse microscopy as previously described56,57. In brief, 20-25-staged embryos were dechorionated by forceps immediately before imaging. Dechorionated embryos were carefully aligned laterally on the mould made by 2% low-melting point agarose (Sigma-Aldrich) on the cell culture dish with glass bottom (CELLview). To prevent the embryos moving during imaging, buffered Tricaine (ethyl 3-aminobenzoate methanesulfonic acid) (Sigma-Aldrich) in a final concentration of 0.016% in E3 medium was used for anaesthesia. Time-lapse images were acquired by an inverted Zeiss Axiovert 200M microscope with an Andor iXOn 888 EM-CCD camera and software controlled motor stage. Temperature during imaging was controlled by a Warner Thermal Cooling Module TCM-1 (Warner Instruments), and recorded during imaging by an electro-thermometer Voltcraft Plus K202. Imaging was started at bud stage (10h.p.f.) and lasted for 16.5 h; the acquisition and image post-processing were performed using Andor IQ2.0 software. To treat with DAPT, embryos at 4 h.p.f. were transferred to 50 mM DAPT or DMSO control in E3 medium. Four-chamber glass bottom dishes holding WT or transgenic embryos with or without DAPT were used for imaging.

To calculate the somitogenesis period, the timing of somitic furrow emergence in the trunk was assessed in each embryo from the time-lapse movies for segments 4-19. The slope of linear regression from time point versus somite number was calculated using Microsoft Excel. Somitogenesis period is the slope multiplied by time interval of the time-lapse movie. The somitogenesis period of each genotype was obtained from at least three independent trials at the same temperature setting.

The measurement of the axis elongation and PSM shortening was performed as previously published2,18,39. Briefly, in the time-lapse movie collected above, a measurement was made at each relevant time interval using the segmented line tool in Fiji. This line extended from the posterior tip of the tailbud to either the anterior end of the head or the most recently formed somite boundary, respectively, and followed the course of the notochord. The time 0 was defined as 5 somite-stage in the present study.

Measurement of anterior wavelength from herl whole-mount ISH. A total of 4 to 20 somite-stage embryos were stained with the her1 riboprobe according to published ISH protocol58 and imaged laterally on an Olympus SZX-10 stereomicroscope equipped with a QImaging Micropublisher 5.0 RTV camera. TIFF-formatted images were converted to inverted grayscale images from the red channel to give higher signal/noise ratio than the original RGB images.

To measure the intensity distribution along the axis, segment lines with 35 pixels (36.6 mm) width were drawn from anterior to the most anterior her1 stripe to the end of PSM (Fig. 3a, top-left panel). The anterior wavelengths were

calculated from the intervals between the first two peaks (Fig. 3a). The anterior border of the first peak, X0, is defined by the intersection of the background line with the linear regression line with the highest r2 value starting from the inflection point of the anterior side of the peak, extending anteriorly along the curve of the peak. her1 expression length, her1-L, is the length in pixels from the X0 to the end of line measurement.

Measurement of cell densities/areas in isletl whole-mount ISH. Four somite-stage embryos were hybridized with the islet1 riboprobe and imaged as described above. TIFF-formatted images were adjusted by auto contrast in Adobe Photoshop CS4 and converted to grayscale images from the red channel.

To measure the cell density of PMN and RB, we used a plug-in, Peak Finder (Johannes Schindelin, unpublished) in Fiji59. Along the PMN or RB in the anterior axis (Fig. 4a, lower right panel), we drew a line 5 pixels (3.78 mm) wide and 228.1 ± 40.2 mm (PMN) or 233.9 ± 40.4 mm (RB) long, then analyzed the number of peaks in pixel intensity along this line by Peak Finder. To measure TG areas, monochrome bitmap images were generated by watershedding images using Adobe Photoshop CS4 with various levels to match the black areas resembling the TG areas. TG areas were thus obtained from the number of black pixels.

Calculation of segmentation period from anterior wavelength. In Damascus,

the anterior wavelength is shorter than in WT, which we expect to affect the frequency contribution from the Doppler effect (Supplementary Note 1). Here, we tested how this altered Doppler contribution relates to the observed difference in the segmentation rate between WT and Damascus. As a reference for the relative contribution of the Doppler effect to the total rate of segmentation, we used the measurements on the Looping live reporter2. From this we estimated that the Doppler effect accounts for 22% of the observed rate of segmentation (Supplementary Note 2). The rate of segmentation in WT at 28.5 0C is R = T—1 = 0.0405 ± 0.0009 min — 1, where T is the segmentation period, which yielded an expected Doppler contribution of D = 22% • R = 0.0089 ± 0.0002 min — 1. Using the fact that the tissue shortening speed is the same in WT and Damascus and taking into account the effects of the altered wavelength on the Doppler contribution, we expect in Damascus a Doppler contribution of

D = ~ D (1)

where 1 is the anterior wavelength in WT and 1 is the anterior wavelength in Damascus (Supplementary Text 1). For the case that the altered Doppler contribution is the only change in the wave pattern affecting the segmentation rate, we obtained an estimate for the segmentation rate in Damascus as R = R + (D — D) and the segmentation period follows as ~ = (R) — . The measured relative wavelength 1/1, the inferred Doppler contribution D and resulting segmentation period T~ at each of the time-points measured in Fig. 3b is given in Supplementary Table 1.

Illustration of the effects of coupling on the wave pattern. We used a theory of coupled oscillators that takes coupling delays into account to demonstrate how the wavelength of the resulting wave pattern depends on coupling (Supplementary Note 3). We described the wave pattern in the PSM by a one-dimensional phase field along the elongating body axis in a reference frame co-moving with the tailbud tip (Supplementary Fig. 4a). In this reference frame, cells move in an anterior direction to become part of the formed segments, represented by the arrested wave pattern. Wave patterns emerge from a position-dependent frequency profile along the body axis with maximum frequency at the tailbud tip and oscillations gradually slow down as cells approach the anterior end (Supplementary Fig. 4b). For a constant PSM length, this system attains a state with a time-periodic wave pattern. The time-independent phase profile describes the spatial distribution of the waves along the PSM (Supplementary Fig. 4c-f). The wave pattern can be characterized by its local wavelength and the total number of waves in the PSM. Here we focused on the wavelength at the anterior end, which is experimentally evaluated in this study (Fig. 3b; Supplementary Fig. 2). Not all solutions for wave patterns were robust under small perturbations of the phase field. In certain parameter regions, small perturbations grow exponentially and prevent the system from forming a stable wave pattern. We performed a linear stability analysis that reveals the regions in which the pattern is stable (Supplementary Note 3). In certain parameter regions, the system exhibited multiple solutions for a given set of parameters. However, here we only considered parameter regions that enable a single stable solution, which corresponds to a robust pattern forming system (Supplementary Fig. 5A). In the regions with stable wave patterns, we simulated the system numerically and characterized the resulting wave patterns by their wavelength and the number of waves (Supplementary Fig. 5). Alterations in the strength and delay of coupling led to a variety of different wavelengths and number of waves in the PSM (Supplementary Figs 4C-F and 5C,D).

Statistics. Values were presented as mean ± s.d. unless otherwise described. Student's t-test with the parameters of two-tailed and unequal variance was performed in Microsoft Excel.

Data availability. The data that support the findings of this study are available

from the corresponding author upon request.

References

1. Palmeirim, I., Henrique, D., Ish-Horowicz, D. & Pourquie, O. Avian hairy gene expression identifies a molecular clock linked to vertebrate segmentation and somitogenesis. Cell 91, 639-648 (1997).

2. Soroldoni, D. et al. Genetic oscillations. A Doppler effect in embryonic pattern formation. Science 345, 222-225 (2014).

3. Oates, A. C., Morelli, L. G. & Ares, S. Patterning embryos with oscillations: structure, function and dynamics of the vertebrate segmentation clock. Development 139, 625-639 (2012).

4. Harima, Y., Takashima, Y., Ueda, Y., Ohtsuka, T. & Kageyama, R. Accelerating the tempo of the segmentation clock by reducing the number of introns in the Hes7 gene. Cell Rep. 3, 1-7 (2013).

5. Krol, A. J. et al. Evolutionary plasticity of segmentation clock networks. Development 138, 2783-2792 (2011).

6. Delaune, E. A., Francois, P., Shih, N. P. & Amacher, S. L. Single-cell-resolution imaging of the impact of Notch signalling and mitosis on segmentation clock dynamics. Dev. Cell 23, 995-1005 (2012).

7. Horikawa, K., Ishimatsu, K., Yoshimoto, E., Kondo, S. & Takeda, H. Noise-resistant and synchronized oscillation of the segmentation clock. Nature 441, 719-723 (2006).

8. Jiang, Y. J. et al. Notch signalling and the synchronization of the somite segmentation clock. Nature 408, 475-479 (2000).

9. Mara, A., Schroeder, J., Chalouni, C. & Holley, S. A. Priming, initiation and synchronization of the segmentation clock by deltaD and deltaC. Nat. Cell Biol. 9, 523-530 (2007).

10. Ozbudak, E. M. & Lewis, J. Notch signalling synchronizes the zebrafish segmentation clock but is not needed to create somite boundaries. PLoS Genet. 4, e15 (2008).

11. Riedel-Kruse, I. H., Muller, C. & Oates, A. C. Synchrony dynamics during initiation, failure, and rescue of the segmentation clock. Science 317, 1911-1915 (2007).

12. Soza-Ried, C., Ozturk, E., Ish-Horowicz, D. & Lewis, J. Pulses of Notch activation synchronise oscillating somite cells and entrain the zebrafish segmentation clock. Development 141, 1780-1788 (2014).

13. Wright, G. J. et al. DeltaC and DeltaD interact as Notch ligands in the zebrafish segmentation clock. Development 138, 2947-2956 (2011).

14. Itoh, M. et al. Mind bomb is a ubiquitin ligase that is essential for efficient activation of Notch signalling by Delta. Dev. Cell 4, 67-82 (2003).

15. Zhang, C., Li, Q., Lim, C. H., Qiu, X. & Jiang, Y. J. The characterization of zebrafish antimorphic mib alleles reveals that Mib and Mind bomb-2 (Mib2) function redundantly. Dev. Biol. 305, 14-27 (2007).

16. Holley, S. A., Geisler, R. & Nusslein-Volhard, C. Control of her1 expression during zebrafish somitogenesis by a delta-dependent oscillator and an independent wave-front activity. Genes Dev. 14, 1678-1690 (2000).

17. van Eeden, F. J. et al. Mutations affecting somite formation and patterning in the zebrafish, Danio rerio. Development 123, 153-164 (1996).

18. Herrgen, L. et al. Intercellular coupling regulates the period of the segmentation clock. Curr. Biol. 20, 1244-1253 (2010).

19. Jen, W. C., Wettstein, D., Turner, D., Chitnis, A. & Kintner, C. The Notch ligand, X-Delta-2, mediates segmentation of the paraxial mesoderm in Xenopus embryos. Development 124, 1169-1178 (1997).

20. Julich, D. et al. beamter/deltaC and the role of Notch ligands in the zebrafish somite segmentation, hindbrain neurogenesis and hypochord differentiation. Dev. Biol. 286, 391-404 (2005).

21. Takke, C. & Campos-Ortega, J. A. her1, a zebrafish pair-rule like gene, acts downstream of notch signalling to control somite development. Development 126, 3005-3014 (1999).

22. Hans, S. & Campos-Ortega, J. A. On the organisation of the regulatory region of the zebrafish deltaD gene. Development 129, 4773-4784 (2002).

23. Soroldoni, D., Hogan, B. M. & Oates, A. C. Simple and efficient transgenesis with meganuclease constructs in zebrafish. Methods Mol. Biol. 546, 117-130 (2009).

24. Geling, A., Steiner, H., Willem, M., Bally-Cuif, L. & Haass, C.

A gamma-secretase inhibitor blocks Notch signalling in vivo and causes a severe neurogenic phenotype in zebrafish. EMBO Rep. 3, 688-694 (2002).

25. Jorg, D. J., Morelli, L. G., Soroldoni, D., Oates, A. C. & Julicher, F. Continuum theory of gene expression waves during vertebrate segmentation. New J. Phys. 17, 093042 (2015).

26. Ares, S., Morelli, L. G., Jorg, D. J., Oates, A. C. & Julicher, F. Collective modes of coupled phase oscillators with delayed coupling. Phys. Rev. Lett. 108, 204101 (2012).

27. Morelli, L. G. etal. Delayed coupling theory of vertebrate segmentation. HFSP J. 3, 55-66 (2009).

28. Collier, J. R., Monk, N. A., Maini, P. K. & Lewis, J. H. Pattern formation by lateral inhibition with feedback: a mathematical model of delta-notch intercellular signalling. J. Theor. Biol. 183, 429-446 (1996).

29. Matsuda, M., Koga, M., Woltjen, K., Nishida, E. & Ebisuya, M. Synthetic lateral inhibition governs cell-type bifurcation with robust ratios. Nat. Commun. 6, 6195 (2015).

30. Sprinzak, D. et al. Cis-interactions between Notch and Delta generate mutually exclusive signalling states. Nature 465, 86-90 (2010).

31. Henrique, D. et al. Expression of a Delta homologue in prospective neurons in the chick. Nature 375, 787-790 (1995).

32. Chitnis, A. & Kintner, C. Sensitivity of proneural genes to lateral inhibition affects the pattern of primary neurons in Xenopus embryos. Development 122, 2295-2301 (1996).

33. Pan, D. & Rubin, G. M. Kuzbanian controls proteolytic processing of Notch and mediates lateral inhibition during Drosophila and vertebrate neurogenesis. Cell 90, 271-280 (1997).

34. Haddon, C. et al. Multiple delta genes and lateral inhibition in zebrafish primary neurogenesis. Development 125, 359-370 (1998).

35. Inoue, A., Takahashi, M., Hatta, K., Hotta, Y. & Okamoto, H. Developmental regulation of islet-1 mRNA expression during neuronal differentiation in embryonic zebrafish. Dev. Dyn. 199, 1-11 (1994).

36. Heuss, S. F., Ndiaye-Lobry, D., Six, E. M., Israel, A. & Logeat, F. The intracellular region of Notch ligands Dll1 and Dll3 regulates their trafficking and signalling activity. Proc. Natl Acad. Sci. USA 105, 11212-11217 (2008).

37. Kim, W. et al. The period of the somite segmentation clock is sensitive to Notch activity. Mol. Biol. Cell 22, 3541-3549 (2011).

38. Aulehla, A. et al. A beta-catenin gradient links the clock and wavefront systems in mouse embryo segmentation. Nat. Cell Biol. 10, 186-193 (2008).

39. Bajard, L. et al. Wnt-regulated dynamics of positional information in zebrafish somitogenesis. Development 141, 1381-1391 (2014).

40. Dubrulle, J., McGrew, M. J. & Pourquie, O. FGF signalling controls somite boundary position and regulates segmentation clock control of spatiotemporal Hox gene activation. Cell 106, 219-232 (2001).

41. Ishimatsu, K., Takamatsu, A. & Takeda, H. Emergence of traveling waves in the zebrafish segmentation clock. Development 137, 1595-1599 (2010).

42. Sawada, A. et al. Fgf/MAPK signalling is a crucial positional cue in somite boundary formation. Development 128, 4873-4880 (2001).

43. Murray, P. J., Maini, P. K. & Baker, R. E. The clock and wavefront model revisited. J. Theor. Biol. 283, 227-238 (2011).

44. Schroter, C. et al. Topology and dynamics of the zebrafish segmentation clock core circuit. PLoS Biol. 10, e1001364 (2012).

45. Schroter, C. & Oates, A. C. Segment number and axial identity in a segmentation clock period mutant. Curr. Biol. 20, 1254-1258 (2010).

46. Meyer, A. & Schartl, M. Gene and genome duplications in vertebrates: the one-to-four (-to-eight in fish) rule and the evolution of novel gene functions. Curr. Opin. Cell Biol. 11, 699-704 (1999).

47. Taylor, J. S., Van de Peer, Y., Braasch, I. & Meyer, A. Comparative genomics provides evidence for an ancient genome duplication event in fish. Philos. Trans. R. Soc. Lond. B Biol. Sci. 356, 1661-1679 (2001).

48. Kimmel, C. B., Ballard, W. W., Kimmel, S. R., Ullmann, B. & Schilling, T. F. Stages of embryonic development of the zebrafish. Dev. Dyn. 203, 253-310 (1995).

49. Sarov, M. et al. A recombineering pipeline for functional genomics applied to Caenorhabditis elegans. Nat. Methods 3, 839-844 (2006).

50. Fu, J., Teucher, M., Anastassiadis, K., Skarnes, W. & Stewart, A. F. A recombineering pipeline to make conditional targeting constructs. Methods Enzymol. 477, 125-144 (2010).

51. Ejsmont, R. K. et al. Recombination-mediated genetic engineering of large genomic DNA transgenes. Methods Mol. Biol. 772, 445-458 (2011).

52. Inuzuka, M. & Wada, Y. An initiator protein for plasmid R6K DNA replication. Mutations affecting the copy-number control. FEBS Lett. 228, 7-11 (1988).

53. Thermes, V. et al. I-SceI meganuclease mediates highly efficient transgenesis in fish. Mech. Dev. 118, 91-98 (2002).

54. Ogino, H., McConnell, W. B. & Grainger, R. M. High-throughput transgenesis in Xenopus using I-SceI meganuclease. Nat. Protoc. 1, 1703-1710 (2006).

55. Ogino, H., McConnell, W. B. & Grainger, R. M. Highly efficient transgenesis in Xenopus tropicalis using I-SceI meganuclease. Mech. Dev. 123, 103-113 (2006).

56. Herrgen, L., Schroter, C., Bajard, L. & Oates, A. C. Multiple embryo time-lapse imaging of zebrafish development. Methods Mol. Biol. 546, 243-254 (2009).

57. Schroter, C. et al. Dynamics of zebrafish somitogenesis. Dev. Dyn. 237, 545-553 (2008).

58. Thisse, C. & Thisse, B. High-resolution in situ hybridization to whole-mount zebrafish embryos. Nat. Protoc. 3, 59-69 (2008).

59. Schindelin, J. et al. Fiji: an open-source platform for biological-image analysis. Nat. Methods 9, 676-682 (2012).

Acknowledgements

Thanks to Pavel Tomancak, Laurel Rohde, Daniele Soroldoni, Alexis Webb and Saul Ares

for comments on the manuscript and to Luis Morelli and members of the Oates and

Julicher groups for discussion. This work was supported by the Francis Crick Institute,

which receives its core funding from the Cancer Research UK, the UK Medical Research

Council and the Wellcome Trust, as well as the Max Planck Society through the Max

Planck Institute for Molecular Cell Biology and Genetics. In addition, the research was supported by the European Research Council under the European Communities 7th Framework Programme (FP7/2007-2013)/(ERC grant 207634), the Wellcome Trust (WT098025MA) and the Medical Research Council (MC_UP_1202/3).

Author contributions

B.-K.L. and A.C.O. conceived the study; B.-K.L. carried out the experiments; all authors analyzed the data; D.J.J. made theoretical calculations, and all authors wrote, revised and edited the manuscript.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/ naturecommunications

Competing financial interests: The authors declare no competing financial interests.

Reprints and permission information is available online at http://npg.nature.com/ reprintsandpermissions/

How to cite this article: Liao, B.-K. et al. Faster embryonic segmentation through elevated Delta-Notch signalling. Nat. Commun. 7:11861 doi: 10.1038/ncomms11861 (2016).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/