Optogenetic investigation of BMP target gene expression diversity

  1. Katherine W Rogers
  2. Mohammad ElGamacy
  3. Benjamin M Jordan
  4. Patrick Müller  Is a corresponding author
  1. Systems Biology of Development Group, Friedrich Miescher Laboratory of the Max Planck Society, Germany
  2. Modeling Tumorigenesis Group, Translational Oncology Division, Eberhard Karls University Tübingen, Germany
  3. Heliopolis Biotechnology Ltd, United Kingdom
  4. Department of Organismic and Evolutionary Biology, Harvard University, United States
9 figures, 3 tables and 2 additional files

Figures

Figure 1 with 1 supplement
BMP target genes have diverse spatial expression patterns at shield stage.

(A) The gradient threshold model states that a signaling gradient activates genes (blue, white, red) at different thresholds (dashed gray lines). (B) BMP binding induces receptor complex formation, phosphorylation of Smad1/5/9, and activation of target genes. (C) Schematic of shield-stage zebrafish embryos with BMP signaling gradients (magenta) along the dorsal-ventral axis. (D-E) Representative images (D) of pSmad1/5/9 immunofluorescence in embryos at the indicated time post-fertilization and quantification (E). (F-O’) Fluorescence in situ hybridization (FISH) showing spatial expression of the indicated high-confidence BMP target genes at shield stage (~6.75 h post-fertilization (hpf)). (F,G,H,I,J,K,L,M,N,O) are animal views, dorsal to the right. (F’,G’,H’,I’,J’,K’,L’,M’,N’,O’) are ventral views. Vertical white bars indicate regions where expression is excluded from the margin. (P-Y) Quantification of FISH signal along the dorsal-ventral axis for the indicated BMP target genes at shield stage (ventral on the left, dorsal on the right as in (E)). Normalized intensities are shown; error bars represent standard error. The Gaussian function Ae-x-μ2ς was fitted to each profile (gray lines), and gene expression range was defined as r=μ+2ς/2 (gray bars). Some BMP target genes could not be reliably quantified due to weak FISH signal (bmp4, id2a, smad6a, smad7, and znfl2b) or inability to reliably identify the ventral side in all assays (crabp2b). See the Figure 1—source data 1 file for source data.

Figure 1—figure supplement 1
BMP target gene identification and spatial quantification.

(A) Transcript levels in uninjected shield-stage (~6.75 hpf) embryos and siblings injected with 10 pg bmp2b or 100 pg chordin (Chd, BMP inhibitor) mRNA were assessed by RNA-sequencing (Supplementary file 1). (B) 16 high-confidence BMP target genes were both upregulated in bmp-overexpressing embryos and downregulated in chordin-overexpressing embryos. All RNA-seq data shown here is classified as significantly differentially regulated by edgeR, DESeq, and Cuff diff. (C) Embryos treated as in (A) were fixed at shield stage and BMP target gene expression assessed by in situ hybridization. (D) Opto-BMP (Figure 3—figure supplement 1A) and uninjected embryos were exposed to DMSO or 50 μg/ml of the translation inhibitor cycloheximide (CHX) at 6.25 hpf, illuminated with blue light from 6.75 to 7.25 hpf, and fixed 20 min post-exposure. A subset of BMP target genes was assessed by in situ hybridization. (E) Embryos were imaged using selective plane illumination microscopy to assess BMP signaling (pSmad1/5/9 immunofluorescence) or BMP target gene expression (fluorescence in situ hybridization). Average intensities along vertical lines from the ventral to dorsal side were calculated. ROIs excluded the embryo center for quantification of margin-restricted genes. See the Figure 1—figure supplement 1—source data 1 file for source data.

BMP target genes have diverse temporal expression profiles.

(A–N) Embryos were collected every 30 min from 2.75 to 7.25 hpf, and transcript levels were quantified using NanoString technology. Error bars represent standard error. Temporal profiles were fit with the cumulative distribution function of the normal distribution (gray lines), and activation time (arrowheads) was defined as the average time point at which the curves reached about two mean average deviations (i.e., 1.5τ) from the inflection point ν (excluding the maternally deposited genes id2a [Chong et al., 2005] and smad6a [White et al., 2017]). NanoString probes for two high-confidence activated BMP target genes (apoc1l and znfl2b) were not functional. (O) Average gene expression spatial range is plotted against average activation time. See the Figure 2—source data 1 file for source data.

Figure 3 with 1 supplement
Reversible activation of BMP signaling using blue light-activated Opto-BMP.

(A) Schematic of Opto-BMP strategy. Blue light-dimerizable VfLOV domains were fused to zebrafish BMP receptor kinase domains. Blue light exposure activates BMP signaling. (B) Embryos injected with mRNA encoding Opto-BMP at the one-cell stage and their uninjected siblings were reared in the dark or exposed to blue light for 10 h starting 70–80 min post-fertilization. Ventralization phenotypes V1-V4 (indicating excess BMP signaling) were scored at 1 day post-fertilization. Number of embryos: uninjected dark = 59, Opto-BMP dark = 53, uninjected light = 55, Opto-BMP light = 60. (C-E) Uninjected and Opto-BMP-injected embryos were exposed to blue light (2300 lux) for 30 min starting at high stage (3.5 hpf) or shield stage (6.75 hpf) and fixed during and after exposure. pSmad1/5/9 immunofluorescence was quantified and plotted in (C) as Opto-BMP minus uninjected signal with piecewise linear interpolation between timepoints; error bars represent standard error (see Materials and methods for statistical analysis). Blue background represents light exposure. Representative embryos from the high-stage (D) and shield-stage (E) experiments quantified in (C). pSmad1/5/9 signal is shown in magenta, DAPI in cyan. See the Figure 3—source data 1 file for source data.

Figure 3—figure supplement 1
Opto-BMP characterization.

(A) Embryos injected with Opto-BMP-encoding mRNA at the one-cell stage were reared in the dark or exposed to white light, and phenotypes were scored at 1 day post-fertilization (dpf). Number of embryos analyzed is indicated on top of each bar. A combination of Opto-Alk3, Opto-Alk8, and Opto-BMPR2a (arrowheads) was used in subsequent experiments. Equimolar amounts of mRNA were injected: Opto-Alk3: 5.2 pg, Opto-Alk8: 5.2 pg, Opto-BMPR2a: 8.9 pg, Opto-BMPR2b: 5 pg. (B) Blue LEDs were embedded into the lid of a six-well dish. (C) LED intensity was measured three times with a luxmeter at different voltages. Averages are shown with piecewise linear interpolation; error bars represent standard error. (D) Circuit schematic for the LED array. A program-generated square wave on a single-board computer acts as a signal source (S1) controlling an amplifier to time the operation and intensity of an LED array (D1-6) run on a constant-voltage power supply (V1). S1: Square wave signal source. C1-3: 1000 μF, 10 V electrolytic capacitor, V1: Constant-voltage regulated power supply, D1-6: 470 nm blue light diode, R1: 200 Ω resistor, R2: 68 Ω, 2W choke resistor, R3: 70 Ω resistor, R4: 7 kΩ resistor, R5: 20 kΩ resistor. (E-H) Opto-BMP embryos were exposed to a 5–12 μm-thick sheet of 405 nm light at the indicated laser power for 3 min using SPIM. BMP signaling was assessed using pSmad1/5/9 immunofluorescence (animal views). Opto-BMP embryos exhibit localized BMP activation compared to uninjected (E) and mock-exposed Opto-BMP embryos (F). (I-J) Uninjected and Opto-BMP-injected embryos were exposed to blue light (2300 lux) for 30 min starting at high stage (3.5 hpf, I) or shield stage (6.75 hpf, J). pSmad1/5/9 immunofluorescence was quantified and plotted with piecewise linear interpolation between timepoints; error bars represent standard error (Materials and methods). Blue background represents light exposure. This data is shown in Figure 3C as Opto-BMP - uninjected intensity. (K) Uninjected and Opto-BMP embryos were reared in the dark or exposed transiently to blue light (2300 lux) for 30 min or 10 h at the indicated developmental stages. (L) Phenotypes were scored at 1 dpf; number of embryos analyzed is indicated to the right of each bar. See the Figure 3—figure supplement 1—source data 1 file for source data.

Most BMP target genes are competent to respond to an early BMP signaling pulse.

(A) Schematic of competence model: Late-activated genes should respond to a late (shield stage, solid line), but not early (high stage, dashed line) BMP signaling pulse. (B-O) High-confidence BMP target gene responses after an early (high stage,~3.5 hpf, dashed line) or late (shield stage,~6.75 hpf, solid line) BMP signaling pulse delivered by exposing uninjected and Opto-BMP-injected embryos to 30 min blue light (Figure 3C–E and Figure 3—figure supplement 1I,J). To assess induced transcription, NanoString transcript counts from uninjected embryos were subtracted from Opto-BMP transcript counts and are plotted here with piecewise linear interpolation between timepoints; error bars represent standard error (see Materials and methods for statistical analysis). See the Figure 4—source data 1 file for source data.

Figure 5 with 1 supplement
Differential expression kinetics do not fully explain BMP target gene spatiotemporal diversity.

(A) Different transcription kinetics may lead to differences in apparent activation times (arrowheads) based on assay detection thresholds (gray line). Assuming similar degradation kinetics, transcripts with faster induction rates should accumulate to higher levels in response to BMP. (B-E) Uninjected and Opto-BMP-injected embryos were exposed to blue light for 30 min at high (~3.5 hpf, B,C) or shield stage (~6.75 hpf, D,E), and target gene expression in response to the resulting BMP signaling pulses (Figure 3C–E and Figure 3—figure supplement 1I,J) was quantified using NanoString technology (Figure 4). Maximum average transcript counts were determined, and are plotted against activation time (B,D) (Figure 2) or spatial range (C,E). Error bars represent standard error, gray lines represent linear fits, ρs = Spearman correlation coefficient, ρp = Pearson correlation coefficient. crabp2b is not included due to lack of significant induction. (F-L) All three target gene response repeats were fitted with a model of induction and decay (Materials and methods). The average induction constant (σ) is plotted against activation time (F), spatial range (G), or maximum transcript count (J). The average decay constant (λ) is plotted against activation time (H), range (I), or maximum transcript count (K). Error bars represent standard error, ρs = Spearman correlation coefficient, ρp = Pearson correlation coefficient. crabp2b is not included due to lack of significant induction. pSmad1/5/9 immunofluorescence (Figure 3) was fitted with a polynomial (gray line, L) and used as signaling input. (M-Z) Individual fits of transcriptional responses (Figure 4); closed circles represent averages of three data points, open circles represent individual data points, and gray lines represent individual fits of each repeat. See the Figure 5—source data 1 source data file for source data.

Figure 5—figure supplement 1
Opto-BMP-induced BMP target gene responses and alternative fitting method.

(A–F) Opto-BMP and uninjected embryos were exposed to blue light (2300 lux) for 30 min at shield stage (~6.75 hpf), and BMP target gene expression was assessed using colorimetric in situ hybridization in response to the BMP signaling pulse (Figure 3C and Figure 3—figure supplement 1J). Animal pole views. (G-H) Transcript levels in embryos treated as described in A-F were quantified using NanoString technology (Figure 4). BMP target gene induction (σ) and decay constants (λ) were determined by fitting the data using two approaches (Materials and methods): In Method 1 (main text and Figure 5), individual Opto-BMP – uninjected subtracted data series were fitted with a model of induction and decay for each transcript; in Method 2, all raw data points from the Opto-BMP and uninjected data sets were fitted simultaneously for each transcript. Both approaches produced similar σ (G) and λ (H) values (gray lines represent linear fits, ρs = Spearman correlation coefficient, ρp = Pearson correlation coefficient, Method one median shown). (I-W) In Method 2, pSmad1/5/9 immunofluorescence (Figure 3C and Figure 3—figure supplement 1J) was fitted with a model of induction and decay (gray line, I) and used as signaling input to model the transcript time series. Fits of uninjected (black) and Opto-BMP (blue) transcript data using Method two are shown in J-W. See the Figure 5—figure supplement 1—source data 1 file for source data.

Figure 6 with 1 supplement
Differential sensitivity to BMP does not fully explain target gene expression diversity.

(A) The activation threshold model predicts that broadly expressed genes will be activated by lower amplitude signaling. (B) pSmad1/5/9 immunofluorescence in uninjected and Opto-BMP-injected embryos exposed to 10 (triangle, dashed line) or 20 (diamond, solid line) min of 70 (light pink) or 3900 (magenta) lux blue light starting at shield stage. Immunofluorescence was quantified and plotted as Opto-BMP signal - uninjected with piecewise linear interpolation between timepoints; error bars represent standard error (see Materials and methods for statistical analysis). Embryos for the transcriptional response experiment were collected 30 (x), 40 (square), or 50 (circle) min after the start of light exposure. (C-F) Transcriptional responses in Opto-BMP embryos exposed to conditions shown in (B) were quantified using NanoString technology and are plotted against spatial range (C,E) or activation time (D,F). Embryos were collected 30 (x), 40 (square), or 50 (circle) min after the start of light exposure. Responses that are not statistically significant are anchored to the x-axis (N.S.; see Materials and methods for statistical analysis). See the Figure 6—source data 1 file for source data.

Figure 6—figure supplement 1
Responses to different amplitudes and durations of BMP signaling.

(A–N) High-confidence BMP target gene responses were measured in zebrafish embryos by NanoString technology after exposure to optogenetically generated BMP signaling pulses with different amplitudes and durations (see Figure 6B for BMP signaling pulse profiles). Uninjected and Opto-BMP-injected embryos were exposed to 10 (triangle, dashed line) or 20 (diamond, solid line) min of 70 (lilac) or 3900 (dark purple) lux blue light starting at shield stage (~6.75 hpf), and collected 30, 40, and 50 min after the beginning of light exposure. To assess induced transcription, average NanoString transcript counts from uninjected embryos were subtracted from Opto-BMP transcript counts and are plotted here with piecewise linear interpolation between timepoints; large symbols represent averages of three data points, small symbols represent individual data points, and error bars represent standard error (see Materials and methods for statistical analysis). See the Figure 6—source data 1 file for source data.

Figure 7 with 3 supplements
FGF and Nodal contribute to the spatial diversity of BMP target genes.

(A) pSmad1/5/9 immunofluorescence in untreated embryos at the indicated times (data also shown in Figure 1E). (B-E) Embryos were treated with 10 μM FGF inhibitor SU-5402 (B), 50 μM Nodal inhibitor SB-505124 (C), or both (D) starting at 2 hpf, or injected with 0.5 pg bmp2b mRNA at the one-cell stage (E). The BMP signaling gradient was quantified along the dorsal-ventral axis at shield stage using pSmad1/5/9 immunofluorescence; error bars represent standard error. Note that the embryos in panels B-E came from different experiments and were processed and imaged on different days, but untreated controls were always siblings of treated embryos and processed and imaged simultaneously. (F-J) BMP target gene expression along the dorsal-ventral axis at shield stage in untreated (F), SU-5402-treated (G), SB-505124-treated (H), SB-505124 + SU-5402-treated (I), and bmp-overexpressing (J) embryos quantified using fluorescence in situ hybridization (untreated data from Figure 1). (K) FGF and Nodal block expression of a subset of BMP target genes at the margin, and restrict BMP signaling in part by activating the BMP inhibitor Chordin. (L) Ventral views of margin-excluded BMP target gene expression at shield stage assessed by FISH in untreated embryos (top row), bmp-overexpressing embryos (0.5 pg bmp2b mRNA, middle row), and embryos treated with SU-5402 + SB-505124 (bottom row). Vertical white bars indicate regions where expression is excluded from the margin. (M) Expression levels at the margin quantified by calculating the average normalized intensity from 5–10% embryo length in untreated versus bmp-overexpressing embryos (left) or untreated versus SU-5402+SB-505124-treated embryos (right). Lines connect treated and untreated conditions to visualize shifts, error bars represent standard error. Lower numbers indicate less expression at the margin. (N) Spatial coefficient of variation for the 10 BMP target genes assessed here in untreated (black), SU-5402-treated (yellow), SB-505124-treated (red), SU-5402+SB-505124-treated (salmon), and bmp-overexpressing (magenta) embryos. Lower numbers indicate less spatial diversity. See the Figure 7—source data 1 file for source data.

Figure 7—figure supplement 1
Signaling and BMP target gene expression in bmp-overexpressing and FGF/Nodal-inhibited embryos.

(A–J’) BMP, FGF, and Nodal signaling were assessed at shield stage (~6.75 hpf) using pSmad1/5/9, pErk, or pSmad2/3 immunofluorescence, respectively, in untreated embryos and siblings injected with 0.5 pg bmp2b mRNA (A-C’), or exposed to 10 μM FGF inhibitor SU-5402 (D-E’), 50 μM Nodal inhibitor SB-505124 (F-G’), or both (H-J’) starting at 2 h post-fertilization (hpf). pSmad1/5/9: animal views, pErk: lateral views, pSmad2/3: dorsal views. Quantification is shown in Figure 7B–E. (K-O) Fluorescence in situ hybridization for the indicated 10 high-confidence BMP target genes in shield-stage embryos treated as described in A-J’. Untreated images are also shown in Figure 1; Left columns are animal views (dorsal to the right), and right columns are ventral views. Vertical white bars indicate regions where expression is excluded from the margin. (P-S) Expression levels at the margin quantified by calculating the average normalized intensity from 5–10% embryo length in embryos treated as described in A-J’. P and S are also shown in Figure 7. Lines connect treated and untreated conditions to visualize shifts; error bars represent standard error. See the Figure 7—figure supplement 1—source data 1 file for source data.

Figure 7—figure supplement 2
BMP signaling and temporal target gene expression in the absence of FGF/Nodal signaling.

Embryos were treated with 10 μM FGF inhibitor SU-5402 and 50 μM Nodal inhibitor SB-505124 starting at 2 hpf. (A-K) BMP signaling was assessed from 2.75 to 7.25 h post-fertilization (hpf) in treated and untreated embryos using pSmad1/5/9 immunofluorescence. Spatial BMP signaling gradients along the ventral-to-dorsal axis were quantified (B-K). Error bars represent standard error; data from untreated embryos also shown in Figures 1E and 7A. (L-ZA) BMP target gene transcript levels (L-Y) were assessed using NanoString technology in embryos treated as described above from 2.75 to 7.25 hpf. Temporal profiles were fit with the cumulative distribution function of the normal distribution (gray and light salmon lines), and activation time (black arrowhead = untreated, salmon arrowhead = drug-treated) was defined as the average time point at which the curves reached about two mean average deviations (i.e., 1.5τ) from the inflection point ν (excluding the maternally deposited genes id2a [Chong et al., 2005] and smad6a [White et al., 2017]). Expression of bmp2b (Z) and chordin (ZA) was also assessed. Error bars represent standard error (see Materials and methods for statistical analysis). See the Figure 7—figure supplement 2—source data 1 file for source data.

Figure 7—figure supplement 3
Individual fluorescence in situ hybridization profiles in bmp-overexpressing and FGF/Nodal-inhibited embryos.

Spatial quantification of fluorescence in situ hybridization signal along the dorsal-ventral axis for BMP target genes at shield stage in embryos treated as indicated, compared to untreated spatial distributions (black). Normalized intensities are shown, error bars represent standard error. (A) Embryos were injected with 0.5 pg bmp2b mRNA at the one-cell stage. (B-D) Embryos were treated with 10 μM FGF inhibitor SU-5402 (B), 50 μM Nodal inhibitor SB-505124 (C), or both (D) starting at 2 hpf. See the Figure 7—figure supplement 3—source data 1 file for source data.

Author response image 1
BMP target gene responses to different amplitudes and durations of BMP signalling.

Embryos injected with Opto-BMP and uninjected siblings were either left in the dark (A) or exposed to 70 or 3900 lux blue light for 10 (B) or 20 (C) min starting at shield stage (~6.75 hpf), and fixed 20 min post-exposure. BMP target gene expression in response to the resulting low- and high-amplitude BMP signaling pulses (Figure 6B) was assessed using colorimetric in situ hybridization.

Author response image 2
The activation timing and expression ranges for hypothetical genes with different activation thresholds were calculated based on quantified pSmad1/5/9/ activation kinetics (Figure 1E).

Tables

Key resources table
Reagent type (species)
or resource
DesignationSource or referenceIdentifiersAdditional information
Strain, strain background
(E. coli)
One Shot TOP10Life TechnologiesC4040Chemically competent
Strain, strain background
(Danio rerio)
TE zebrafishPomreinke et al., 2017
Donovan et al., 2017
Wild type
Antibodyanti-phospho-
Smad1/Smad5/Smad9
(Rabbit monoclonal)
Cell Signaling Technology13820,
RRID:AB_2493181
IF (1:100)
Antibodyanti-rabbit Alexa Fluor
488-conjugated
secondary (Goat polyclonal)
Life TechnologiesA11008,
RRID:AB_143165
IF (1:5000)
Antibodyanti-phospho-Smad2/Smad3
(Rabbit monoclonal)
Cell Signaling Technology8828,
RRID:AB_2631089
IF (1:5000)
Antibodyanti-rabbit horseradish
peroxidase (Goat polyclonal)
Jackson ImmunoResearch111-035-003,
RRID:AB_2313567
IF (1:500)
Antibodyanti-pErk (Mouse monoclonal)SigmaM8159,
RRID:AB_477245
IF (1:5000)
Antibodyanti-mouse horseradish
peroxidase (Donkey polyclonal)
Jackson ImmunoResearch715-035-150,
RRID:AB_2340770
IF (1:500)
Antibodyanti-digoxigenin horseradish
peroxidase Fab fragments
(Sheep polyclonal)
Roche11207733910,
RRID:AB_514500
FISH (1:150)
Recombinant
DNA reagent
pCS2-Opto-Alk3Generated in this study
Recombinant DNA reagentpCS2-Opto-Alk8Generated in this study
Recombinant DNA reagentpCS2-Opto-BMPR2aGenerated in this study
Recombinant DNA reagentpCS2-Opto-BMPR2bGenerated in this study
Chemical compound, drugTRIzol reagentInvitrogen5596026
Chemical compound, drugCo-Precipitant PinkBiolineBIO-37075
Chemical compound, drugCycloheximideSigmaC4859
Chemical compound, drugPronaseRoche11459643001
Chemical compound, drugDMSORothA994.2
Chemical compound, drugFBSBiochromS0415
Chemical compound, drugDAPILife TechnologiesD13061:5000
Chemical compound, drugBlocking reagentRoche11096176001
Chemical compound, drugLow melting
temperature agarose
Lonza50080
Chemical compound, drugNodal inhibitor SB-505124SigmaS4696-5MG50 μM
Chemical compound, drugFGF inhibitor SU-5402SigmaSML0443-5MG10 μM
Commercial assay or kitTSA plus cyanine three systemPerkin ElmerNEL744001KTFISH/IF (1:75)
Commercial assay or kitRNeasy kitQIAGEN74104
Commercial assay or kitWizard SV Gel and
PCR Clean-up System
PromegaA9282
Commercial assay or kitpCR-bluntII TOPO kitThermo Fisher Scientific450245
Commercial assay or kitSP6 mMessage
mMachine transcription kit
Thermo Fisher ScientificAM1340
Commercial assay or kitDIG RNA labeling mixSigma-Aldrich11277073910
Software, algorithmFijiSchindelin et al., 2012https://fiji.sc/
RRID:SCR_002285
Software, algorithmPrismGraphPad Softwarehttps://www.graphpad.com/scientific-software/prism
RRID:SCR_002798
Software, algorithmCOMSOL Multiphysics 3.5aCOMSOL, Inchttps://www.comsol.com/
RRID:SCR_014767
Software, algorithmMatlabMathworkshttp://mathworks.com
RRID:SCR_001622
Software, algorithmedgeR 3.2.3Robinson et al., 2010RRID:SCR_012802
Software, algorithmDESeq 1.12.0Anders and Huber, 2010RRID:SCR_000154
Software, algorithmCuff diff 2.1.1Trapnell et al., 2010https://github.com/cole-trapnell-lab/cufflinks
Software, algorithmPWM code for
controlling LED array
Generated in this study
Software, algorithmnSolver 4.0 softwareNanoStringRRID:SCR_003420
Software, algorithmExcelMicrosoftRRID:SCR_016137
Software, algorithmMapleWaterloo Maple IncRRID:SCR_014449
OtherRNA-sequencing dataGenerated in this studyGEO: GSE135100
OtherTIP122 complementary
power NPN Darlington
STMicroelectronics
OtherRegulated power supplyDisrelec Group AGRND 320-KD3000D
Other6-well platesGreiner Bio-One657160
OtherBlue LEDsNichiaNSPB510AS
OtherBlue LEDsEverlight1363-2SUBC/C470/S400-A4
OtherTemperature-
controlled incubator,
Heratherm IMC 18
ThermoScientific50125882
OtherRaspberry Pi model BRaspberry Pi Foundation
OtherLM37 luxmeterDOSTMANN
electronic GmbH
OtherWhite worklightREV Ritter GmbH90910
OtherRed color filtersRoscoE106 Primary Red
Author response table 1
Bootstrapping analysis for Figure 2O.
100 random bootstrapped Pearson correlation coefficients
-0.550.04-0.22-0.16-0.13-0.33-0.11-0.39-0.35-0.02
-0.29-0.38-0.13-0.38-0.09-0.16-0.24-0.10-0.38-0.02
-0.26-0.13-0.450.11-0.08-0.38-0.26-0.13-0.17-0.22
0.02-0.08-0.21-0.33-0.430.03-0.62-0.36-0.24-0.26
-0.71-0.20-0.01-0.66-0.05-0.260.11-0.250.15-0.02
-0.13-0.310.11-0.350.03-0.21-0.16-0.31-0.66-0.31
-0.33-0.36-0.08-0.33-0.66-0.25-0.14-0.39-0.33-0.04
-0.62-0.39-0.13-0.11-0.55-0.39-0.03-0.20-0.08-0.30
-0.22-0.24-0.04-0.250.15-0.14-0.13-0.68-0.13-0.20
-0.160.11-0.22-0.21-0.24-0.68-0.31-0.66-0.28-0.22
Author response table 2
Jackknife analysis for Figure 2O.
Pearson correlation coefficients with the indicated gene removed
bambiacdx4eve1foxi1gata2aklf2bsizzledtfap2cved
-0.21-0.65-0.17-0.35-0.10-0.04-0.13-0.23-0.24

Additional files

Supplementary file 1

RNA-sequencing read counts and differential expression analyses.

Wild type TE zebrafish embryos were injected at the one-cell stage with 10 pg mRNA encoding zebrafish Bmp2b, 100 pg mRNA encoding zebrafish Chordin, or left uninjected. When uninjected siblings reached shield stage (~6.75 hpf), embryos were snap-frozen in liquid nitrogen. 10 embryos were collected per sample, three samples per condition. Total RNA was sequenced and aligned against the reference genome Danio rerio GRCz10 with STAR 2.4.1b. Differential expression analyses comparing BMP-injected versus uninjected (‘+BMP vs. uninj.’ tabs) and Chordin-injected versus uninjected (‘+Chd vs. uninj.’ tabs) were carried out with edgeR 3.2.3, DESeq 1.12.0, and Cuff diff 2.1.1. The p-value threshold for differentially expressed genes was set to 0.05. Read counts and DE analyses for all genes (‘all’ tabs) and for only the significantly differentially expressed genes (‘DE’ tabs) are available here.

https://cdn.elifesciences.org/articles/58641/elife-58641-supp1-v2.xlsx
Transparent reporting form
https://cdn.elifesciences.org/articles/58641/elife-58641-transrepform-v2.docx

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. Katherine W Rogers
  2. Mohammad ElGamacy
  3. Benjamin M Jordan
  4. Patrick Müller
(2020)
Optogenetic investigation of BMP target gene expression diversity
eLife 9:e58641.
https://doi.org/10.7554/eLife.58641