Strength of interactions in the Notch gene regulatory network determines patterning and fate in the notochord

  1. Héctor Sánchez-Iranzo  Is a corresponding author
  2. Aliaksandr Halavatyi
  3. Alba Diz-Muñoz  Is a corresponding author
  1. Cell Biology and Biophysics Unit, European Molecular Biology Laboratory, Germany
  2. Advanced Light Microscopy Facility, European Molecular Biology Laboratory, Germany
7 figures, 2 tables and 10 additional files

Figures

Figure 1 with 1 supplement
Jag1a generates a lateral inhibition pattern.

(A) Schematic representation of notochord development. At 18–19 hpf most of the notochord is composed of coin-shaped precursor cells. During the following 8 hours, progressively, in an antero-posterior order, coin-shaped precursor cells begin their differentiation into sheath cells and vacuolated cells. (B) Airyscan confocal section of a zebrafish notochord at 19 hpf using the rcn3:lyn-mNeonGreen transgenic line. (C) Transmission electron microscopy of a zebrafish notochord at 19 hpf. (D) Magnification of boxed area in (C). (E) (Top) Schematic representation of the model for a Lateral Induction Network shows a pair of cells where the ligand in one cell activates NICD release in the neighboring cell. NICD activates ligand expression in its own cell. (Bottom) Representative simulation of this network applied to an array of cells unidimensionally arranged. (F) (Top) Schematic representation of the model for a Lateral Inhibition Network shows a pair of cells where the ligand in one cell activates NICD release in the neighboring cell. NICD activates the expression of the repressor, which in turn inhibits ligand expression. (Bottom) Simulation of this network applied to an array of cells unidimensionally arranged. (G–L) Maximal intensity projection of Airyscan confocal sections of a zebrafish tail at 22 hpf. (J–L) Magnification of boxed area in (G). n = 10 fish. (L) Intensity profile across a horizontal line in panel (L). (M) jag1a:mScarlet and tp1:GFP expression levels across a 1 μm thick horizontal line on a single plane of the image shown in J. Scale bars, 1 μm (D) 10 μm (B, C, J), 100 μm (G).

Figure 1—figure supplement 1
Lateral induction and lateral inhibition simulations, and mRNA expression pattern of jag1a and jag1b.

(A and B) Representative simulations of NICD and Ligand molecules in each cell for the lateral induction (A) and lateral inhibition (B) models. Each line represents a different cell. Related to models and simulations shown in Figure 1E and F. (C) Confocal projections of 24-hpf zebrafish stained with in situ HCR probes against jag1a (magenta) and jag1b (green). (D–F) Airyscan confocal projections at a higher magnification of boxed region in (C). n = 5 fish. (G) jag1a:mScarlet and tp1:GFP expression levels of each individual cell. Shown values are the mean intensity in a manually selected region of the cytoplasm of each independent cell on a single plane of the image shown in Figure 1J. Cells that are not labeled with jag1a:mScarlet or tp1:GFP cannot be distinguished, and were not included in the analysis. (H) Scatter plot of the intensities shown in Figure 1M. Pearson correlation and p-value of the correlation is shown in the plot. Scale bars, 200 μm (C) and 50 μm (D).

Figure 2 with 4 supplements
Jag1a/jag1b inhibition disrupts lateral inhibition pattern, and this pattern correlates with fate.

(A–D) jag1a:mScarlet tp1:GFP 21 hpf fish embryos injected with control MO (A), jag1a + control MO (B), jag1b + control MO (C), or control MO (D), together with lyn-miRFP mRNA to visualize membranes. (E) Quantification of the average number of jag1a-positive cells directly adjacent to each jag1a-positive cell. Two-tailed p-value is shown in the plot. (F) Cell fate of cells expressing jag1a or the tp1:GFP at the coin-shape stage. Quantifications from images as shown in G (standard deviation jag1a = 2.696, tp1 = 2.631; n = 5 fish). (G) Time lapse of optical sections of notochord cells using the tp1:GFP; jag1a:mScarlet double transgenic line. First frame corresponds to 24 hpf. Scale bars, 20μm.

Figure 2—figure supplement 1
In vivo imaging of notochord development.

(A–D) In vivo time-lapse imaging of zebrafish notochords starting acquisition at 22 hpf using the rcn3:mNeonGreen transgenic line. Acquisition was based on a feedback microscopy protocol, where low quality images were first acquired and then analyzed at the time of acquisition to perform high zoom tile scan imaging only in the notochord cells. (A, C) show maximum projection of Airyscan confocal notochord reconstructions. (B, D) show zoomed images of single Airyscan confocal optical sections magnified from (A, C), respectively. Scale bars, 100 μm (A, C), and 20 μm (B, D).

Figure 2—video 1
Time lapse optical section of notochord cells using the tp1:GFP; jag1a:mScarlet double transgenic line.
Figure 2—video 2
Maximal projection of the zebrafish notochord optical planes acquired using the feedback microscopy protocol to optimize quality of the region of interest.
Figure 2—video 3
Selected plane from Video S2 stabilizing and magnifying a specific region of the notochord.
Figure 3 with 2 supplements
her9 and her6 show a complementary pattern to jag1a.

(A) Heatmap showing the expression levels of the zebrafish HES/HEY family genes. Values represent average normalized UMIs in all notochord cells at 18 and 24 hpf. (B) Projection of confocal optical sections of 18 hpf zebrafish stained with in situ HCR probes against her6 (green) and jag1a (magenta). Transmitted light is shown in gray scale. (C) Maximal projection of confocal Airyscan optical sections of the boxed area in (B). (D), Intensity profile of her6 (green) and jag1a (magenta) along a 1 μm thick horizontal line on the in situ HCR shown in (C). (E–G) Magnified views of boxed area in (C), n = 8. (H) Scatter plot of the intensities shown in D. Each point represents her6 and jag1a intensity in a 1-pixel width times 1 um height square. Pearson correlation and p-value of the correlation is shown in the plot. (I–O) Analogous images to (B–H) based on the her9 probe instead of her6 probe, n = 9. Scale bars, 50 μm (B, I), 20 μm (C, E, J, L).

Figure 3—figure supplement 1
her12 expression is not detected in the notochord.

(A–F) UMAP plots showing notochord cells at 10, 14, 18, and 24 hpf. Cells are labeled depending on the developmental stage (A) or using a logarithmic color scale of of emilin3a, her6, her9, jag1a, and jag1b normalized expression (B–F). (G) Projection of confocal optical sections of 18 hpf zebrafish stained with in situ HCR probes against her12 (green) and jag1a (magenta). Transmitted light is shown in gray scale. (H) Maximal projection of confocal Airyscan optical sections of the boxed area in (G). (I) Intensity profile of her6 (green) and jag1a (magenta) along the notochord based on in situ HCR shown in (H). (J–M) Magnified views of boxed area in (H) n = 8. (M), Scatter plot of the intensities shown in (I). R represents Pearson’s correlation. (N) Guides used for her6 and her9 deletion and rbpj binding sites. Distances are represented as number of base pairs related taking the transcription start (TSS) site as cero. Only binding sites closer than 20 kb from the TSS are shown. Scale bars, 50 μm (G) 20 μm (H, J).

Figure 3—figure supplement 2
Expression of Notch receptors.

(A–D) UMAP plots showing notochord cells at 10, 14, 18, and 24 hpf. Cells are labeled based on a logarithmic color scale of notch1a (A), notch1b (B), notch2 (C) and notch3 (D) normalized expression. (E) Heatmap showing the expression levels of the Notch receptor genes. Values represent average normalized UMIs in all notochord cells at 18 and 24 hpf. (F) Projection of confocal optical sections of 18 hpf zebrafish stained with in situ HCR probe against notch2 (green). Transmitted light is shown in gray scale, n = 11. (G) Projection of confocal Airyscan optical sections of the boxed area in (F). Scale bars, 50 (F) 20 μm (G).

Figure 4 with 2 supplements
her6 and her9 inhibit jag1a expression.

(A – E) Airyscan confocal optical sections of live 22 hpf transgenic jag1a:mScarlet zebrafish injected with emilin3a:GFP (A and B), emilin3a:GFP-p2a-her6 (C and D) or emilin3a:GFP-p2a-her9 (E and F). DNA constructs were injected at the one-cell stage together with I-SceI protein. (B, D, F) show the boundary of GFP segmentation in A, C, and E, respectively, and manual outline of the notochord. (G), Quantification of jag1a:mScarlet intensity inside GFP-positive cells segmented as exemplified in (B, D, F). Values in the plot represent the intensity of jag1a:mScarlet inside segmented cells divided by the jag1a:mScarlet intensity inside the notochord outside of the segmented cells. Each point represents an individual fish (n = 12 GFP, n = 7 GFP-p2a-her6, n = 11 GFP-p2a-her9). Two-tailed p-values are shown in the plot. (H), Airyscan confocal sections of embryo at 22 hpf injected with Cas9 together with a control guide (H) or Cas9 together with her6/her9 gRNAs (I). (J) Quantification of the average number of jag1a-positive cells directly adjacent to each jag1a-positive cell. For each individual fish, we count how many jag1a-positive cells are adjacent to each jag1a-postive cell, and then calculate the average for that fish. This value would be equal to 2 in case all the cells are jag1a-positive, and zero if no jag1a-positive cell is adjacent to another jag1a-positive cell. Each individual point in the plot represents the average value for an independent fish (n = 7 control, n = 8 her6/her9 KOs). Two-tailed p-value is shown in the plot. Scale bars, 20 μm.

Figure 4—figure supplement 1
emilin3a-5kb promoter drives expression to the notochord.

(A and B) Boxplot representing emilin3a expression in different tissues at 18 (A) and 24 hpf (B) according to single-cell transcriptomics. (C) Maximum projection of optical confocal sections of an emilin3a:mScarlet transgenic line at 21 hpf. Transmitted light is shown in gray scale. (D – I) Airyscan confocal optical sections of 21 hpf double transgenic emilin3a:mScarlet/tp1:GFP (D–F) or emilin3a:mScarlet/jag1a:mNeonGreen. Although at different intensities, emilin3a expression can be detected in most of the cells. Arrowheads point to the tp1:GFP- and jag1a:mNG-positive cells to facilitate their identification in the emilin3a:mScarlet channel. (J) 98% ± 3% of the tp1:GFP 99% ± 2% of the jag1a:mNeonGreen cells are emilin3a:mScarlet positive. Cells with an emilin3a:mScarlet signal that is more intense than the background were considered emilin3a:mScarlet-positive. Each dot represents an individual fish (n = 10 tp1:GFP, n = 6 jag1a:mNG). (n. s.) non-significant (0.67) differences by two-tailed t-test. Scale bars, 200 μm (C), 10 μm (D, G).

Figure 4—figure supplement 2
jag1a mRNA expression upon her6 or her9 overexpression.

(A – F) Airyscan confocal optical sections of fixed 22 hpf transgenic injected with emilin3a:GFP (A and B), emilin3a:GFP-p2a-her6 (C and D) or emilin3a:GFP-p2a-her9 (E and F) constructs. GFP was detected by antibody staining and her6 and her9 mRNA by in situ HCR in whole mount embryos. (B, D, F) show the boundary of GFP segmentation in A, C, and E, respectively, and manual outline of the notochord. (G) Quantification of jag1a mRNA intensity inside GFP-positive cells segmented as exemplified in A-F. Each point represents an individual fish (n = 13 GFP, n = 12 GFP-her6, n = 10 GFP-her9). Scale bars, 20 μm.

Figure 5 with 1 supplement
her6, her9, and jag1a determine cell fate in the zebrafish notochord.

(A–C, E–F) Confocal optical sections of 2 dpf live zebrafish that were injected with the emilin3a:GFP (A, E) emilin3a:GFP-p2a-her6 (B), emilin3a:GFP-p2a-her9 (C) or emilin3a:GFP-p2a-jag1a (F) constructs. DNA constructs were injected at the one-cell stage together with I-SceI protein. (D and G) Proportion of vacuolated cells at 2 dpf are shown. Proportion of vacuolated cells was calculated by counting the number of vacuolated cells divided by the sum of the number of sheath and vacuolated cells. Each point in D, G represents an independent fish quantified from on z-stack confocal planes (D, n = 9 GFP, n = 10 GFP-her6, n = 10 GFP-her9, G, n = 5 GFP, n = 5 GFP-jag1a). Two-tailed p-values are shown in D and G. Scale bars, 50 μm.

Figure 5—figure supplement 1
Jag1a intracellular domain does not have an effect on notochord cell fate.

(A–B) Optical sections of 2 dpf zebrafish that were injected with the emilin3a:GFP or emilin3a:GFP-p2a-JICD constructs. (C) Proportion of vacuolated cells at 2 dpf. This proportion was calculated by counting the number of vacuolated cells divided by the sum of the number of sheath and vacuolated cells. Each point represents an independent fish quantified from on z-stack confocal acquisitions (n = 9 GFP, n = 6 GFP-jag1aICD). (n.s.) p-value > 0.05 by two-tailed t-test. Scale bars, 50μm.

Figure 6 with 2 supplements
Modeling and experimental results of cis- and trans-interactions in the notochord.

(A) Two possible circuits may explain the effect of jag1a on fate of the cell where jag1a is overexpressed. Circuit 1 is based on the interaction of ligand and receptor in trans. Circuit 2 is based on a possible role of cis-inhibition of the Notch receptor by the ligand. Cells where we overexpress the ligand are represented as the cell on the left. Adjacent cells are represented on the right. OE, overexpression. (B – G), Airyscan confocal planes of fixed 22 hpf transgenic fish injected with emilin3a:GFP (B–D) or emilin3a:GFP-p2a-jag1a (E–G) constructs. GFP was detected by antibody staining and her6 and her9 mRNA by in situ HCR in whole mount embryos. (C and F) show the notochord outline manually selected and the outline of GFP-positive cells automatically segmented. (D and G) show the outline of the manually selected notochord and the neighborhood to the GFP-positive cells. On the right side of each panel, a magnified view of the boxed region is shown. (H, I) Quantification of her6 and her9 mRNA expression after GFP-based segmentation as shown in (C, F) or (D, G), respectively. Values of her6 and her9 expression levels inside the segmented area inside the notochord were divided by the expression levels of the same genes in the region outside the segmented area, also inside the notochord. Each point represents a different fish. Two-tailed p-values are shown in the plots. n.s., non-significant. (J) Escape rates from the homogeneous steady state (indicated by Maximum Lyapunov Exponents, or MLE) as a function of Kcis and Ktrans parameters. Positive MLE values (red) support patterning, while negative MLE values (blue) do not. Scale bars, 20 μm.

Figure 6—figure supplement 1
Model of lateral inhibition including cis-interactions.

A) Schematic representation of a lateral inhibition model with mutual inhibition. Ligand binding to Notch receptor releases NICD, that activates Repressor expression. The repressor inhibits Ligand expression. In addition, when Notch receptor and Ligand are present in the same cell, they mutually inhibit to each other. (B) Representative simulation of the model shown in (A).

Figure 6—figure supplement 2
Escape rates from the homogeneous steady state.

They are represented by Maximum Lyapunov Exponents, or MLE, as a function of Kcis and Ktrans parameters. Positive MLE values (red) support patterning, while negative MLE values (blue) do not. (A, C, E, G) Non-cooperativity. (B, D, E, H) Cooperativity. βN, production rate of Notch receptor; βL, production rate of Ligand, βR, production rate of receptor. Panel H has been duplicated from Figure 5J for easier comparison with the rest of plots.

Author response image 1
Validation of her6/her9 transient knock-out generation.

(A, B), Confirmation by PCR of her6 (A) and her9 (B) gene deletion after the injection of the Cas9 protein together with guides targeting her6 and her9 genes. Expected PCR product after deletion: 519 bp for her6 and 225 bp for her9. In some cases, due to repair not exactly as expected, different PCR product size can be detected. Only embryos were a her6 and her9 bands were detected were included in the analysis. (C) Original plot including all the imaged fish. (D) New plot including only the samples labeled in green.

Tables

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
AntibodyGFP-Booster Alexa Fluor 488.
anti-GFP Alpaca/recombinant VHH domain, monoclonal.
Chromotekgb2AF488Used at (1:500) dilution.
Genetic reagent (Danio rerio)jag1a:mScarletThis paper“Animal handling and generation of transgenic lines” Methods section
Genetic reagent (Danio rerio)jag1a:mNeonGreenThis paper“Animal handling and generation of transgenic lines” Methods section
Genetic reagent (Danio rerio)emilin3a:mScarletThis paper“Animal handling and generation of transgenic lines” Methods section
Genetic reagent (Danio rerio)rcn3:lyn-mNeonGreenThis paper“Animal handling and generation of transgenic lines” Methods section
Genetic reagent (Danio rerio)tp1:GFPParsons et al., 2009“Animal handling and generation of transgenic lines” Methods section
Sequence-based reagentPrimer 1: pTarBAC_HA1_iTol2_FThis paperPrimergcgtaagcggggcacatttcattacctctttctccgcacccgacatagatCCCTGCTCGAGCCGGGCCCAAGTG
Sequence-based reagentPrimer 2: pTarBAC_HA2_iTol2_RThis paperPrimergcggggcatgactattggcgcgccggatcgatccttaattaagtctactaATTATGATCCTCTAGATCAGATC
Sequence-based reagentPrimer 3: jag1a_HA1_mScarlet_FThis paperPrimergaggcgtgtggcggctgaagtggtagttttcacagcgacagacacacagacagacaaaccACCATGGTGAGCAAGGGC
Sequence-based reagentPrimer 4: jag1a_HA2_FRT_RThis paperPrimeragcagcacgtgagcggacagcgccgcaaaagttgagctcggtctgagaatGGAGGCTACCATGGAGAAG
Sequence-based reagentPrimer 5: jag1a_HA1_mNG_FThis paperPrimergaggcgtgtggcggctgaagtggtagttttcacagcgacagacacacagacagacaaaccACCATGGTGAGCAAGGGC
Sequence-based reagentPrimer 6: ScaffoldShah et al., 2015PrimerGATCCGCACCGACTCGGTGCCACTTTTTCAAGTTGATAACGGACTAGCCTTATTTTAACTT GCTATTTCTAGCTCTAAAAC
Sequence-based reagentPrimer 7: her6_guide1This paperPrimertaatacgactcactataGGTGGTCGGCGCCCCTCCATgttttagagctagaa
Sequence-based reagentPrimer 8: her6_guide2This paperPrimertaatacgactcactataGGGTGGCCATTCTTTGAAGGgttttagagctagaa
Sequence-based reagentPrimer 9: her9_guide1This paperPrimertaatacgactcactataGGGTGACTGACAGCCCGCGGgttttagagctagaa
Sequence-based reagentPrimer 10: her9_guide2This paperPrimertaatacgactcactataGGGGGAAACCCTGCGGCCGTgttttagagctagaa
Sequence-based reagentPrimer 11: univ_guideWierson et al., 2020PrimertaatacgactcactataGGGAGGCGTTCGGGCCACAGgttttagagctagaa
Sequence-based reagentPrimer 12: her6_FThis paperPrimerGTTTGCTGTTTCTGAGCGGAG
Sequence-based reagentPrimer 13: her6_RThis paperPrimerGGGAAGCACGTCTGAGTCTG
Sequence-based reagentPrimer 14: her9_FThis paperPrimerCCGCGCAGTATGTGAATGC
Sequence-based reagentPrimer 15: her9_RThis paperPrimerACCTTCACAGGCTACAGAACC
Sequence-based reagentControl MOYamamoto et al., 2010MorpholinoCCTCTTACCTCAGTTACAATTTATA
Sequence-based reagentjag1a_splMOYamamoto et al., 2010MorpholinoAAGCCAAACCCGCACATACCCGCAT
Sequence-based reagentjag1b_atgMOYamamoto et al., 2010MorpholinoCTGAACTCCGTCGCAGAATCATGCC
Recombinant DNA reagentmScarlet FRT kan FRTThis paperSequence available in Source data 1 file.
Recombinant DNA reagentmNG FRT kan FRTThis paperSequence available in Source data 2 file.
Recombinant DNA reagentemilin3a mScarletThis paperSequence available in Source data 3 file.
Recombinant DNA reagentrcn3 lyn mNeonGreenThis paperSequence available in Source data 4 file.
Recombinant DNA reagentSP6 lyn-miRFP-pAThis paperSequence available in Source data 5 file.
Recombinant DNA reagentemilin3a GFPThis paperSequence available in Source data 6 file.
Recombinant DNA reagentemilin3a GFP-p2a-her6This paperSequence available in Source data 7 file.
Recombinant DNA reagentemilin3a GFP-p2a-her9This paperSequence available in Source data 8 file.
Recombinant DNA reagentemilin3a mScarlet-p2a-JICDThis paperSequence available in Source data 9 file.
Commercial assay or kitmMESSAGEmMACHINESP6 Transcription KitThermo Fisher ScientificCat#:AM1340
Commercial assay or kitTricaine (MESAB)Sigma-AldrichCat#:A5040
Software, algorithmCode for image analysis and mathematical modelingThis paperAvailable in github: https://github.com/hsancheziranzo/notochord-lateral-inhibition (copy archived at swh:1:rev:2e5c5fe15e30ea6bacdc0282e1506b44b05415af, Sánchez-Iranzo, 2022)
Table 1
Genes with a highest score for specificity and expression levels in the notochord at 18 hpf.

Expression: Average expression in Notochord cells (normalized UMIs per million). Enrichment: Average expression in notochord cells divided by average expression in the rest of the cells in the fish at 18 hpf. Score: Expression multiplied by enrichment (equivalent to the equation described above).

GeneExpressionEnrichmentScore
emilin3a5055.751118.165653.16
ntd55929.62134.03794.74
col2a1a9083.5067.75633.59
cmn1536.94330.67508.22
loxl5b1264.29330.84418.28
col9a1b1511.20226.72342.62
ta1443.23162.65234.74
LOC100333762642.86303.86195.34
lgals1l12163.8886.59187.38
col9a21331.86103.37137.68
si:ch211-125-g7.4501.98272.60136.84
si:dkey-12l12.1408.15304.23127.17
col9a3754.32140.39105.90
LOC100334188393.62205.6880.96
pmp22b1001.8068.6468.76
si:dkey-99l1.9499.03131.4565.60
si:ch73-23l24.1317.78205.6065.34
lgals2a1018.1764.0865.24
twist2306.83167.2051.60
lect1616.9779.2448.89

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Héctor Sánchez-Iranzo
  2. Aliaksandr Halavatyi
  3. Alba Diz-Muñoz
(2022)
Strength of interactions in the Notch gene regulatory network determines patterning and fate in the notochord
eLife 11:e75429.
https://doi.org/10.7554/eLife.75429