Reply to: Assessing the precision of morphogen gradients in neural tube development | Nature Communications

Blog

HomeHome / Blog / Reply to: Assessing the precision of morphogen gradients in neural tube development | Nature Communications

Oct 15, 2024

Reply to: Assessing the precision of morphogen gradients in neural tube development | Nature Communications

Nature Communications volume 15, Article number: 930 (2024) Cite this article 1009 Accesses 1 Altmetric Metrics details The Original Article was published on 01 February 2024 replying to M. Zagorski

Nature Communications volume 15, Article number: 930 (2024) Cite this article

1009 Accesses

1 Altmetric

Metrics details

The Original Article was published on 01 February 2024

replying to M. Zagorski et al. Nature Communications https://doi.org/10.1038/s41467-024-45148-8 (2024)

In a recent article1, we demonstrated that single morphogen gradients in the developing mouse neural tube (NT) can carry sufficient positional accuracy to explain the patterning precision of progenitor domain boundaries. Zagorski et al. had previously concluded otherwise2, based on methodological inconsistencies that we have revealed. The authors now comment on our work with a Matters Arising letter. We rebut their criticism point by point in the Supplement, and summarize the main aspects here.

The authors criticize our comparison of the different analysis methods (FitEPM, NumEPM, DEEM) using a set of exponential functions. It is a fundamental part of quantitative science to validate all methods used to process data—be it part of experimental, theoretical or numerical work—against known results. If a method fails to provide the correct results for known problems, it is unreasonable to apply it to similar problems, for which the answer is not known. The authors’ indirect approximation (FitEPM) leads to a vast overestimation of the positional error in case of exponential gradients, and is accurate only near the morphogen source1. The same limitations apply to noisy gradient data, and the challenges of background subtraction and smoothening apply to all three methods.

The authors insist that further away from the source the observed gradients are flat, such that our arguments would not apply. First, their approximation (FitEPM) will yield the wrong positional error whenever the mean deviates substantially from an exponential1 (i.e., also in flat parts of the gradient), as it makes use of the exponential shape explicitly. Second, the flat part of the gradients is not reliable biological data. We speculated in our paper that the switch to a flat gradient shape may be due to insufficient imaging depth, but the employed imaging depth remained unknown to us1. We have since received confirmation that Zagorski et al. indeed employed 8-bit imaging (A. Kicheva, personal communication; Supplementary Information). 8-bit imaging only permits the detection of a 28 = 256-fold intensity change. As such, it is technically impossible to detect an exponential gradient beyond 5.5 times its decay length (~110 µm) from the source, which coincides with the point where Zagorski et al. find the transition from an exponential to a flat gradient shape. The physical limits of their imaging setup and the mathematical limitation of their approximation of the positional error make it impossible to evaluate the positional error of the gradients at later time points or further away from the source.

A key conclusion of Zagorki et al.2 is that morphogen gradients are too noisy to specify cell fate in the NT beyond the first 30 h. The authors now argue that their analysis must be considered correct because it 1) provides an explanation for the sensitivity of the progenitor markers to both SHH and BMP, and 2) is consistent with their previous postulate that cell differentiation rather than morphogen gradients define the progenitor domain boundary positions at later stages3. However, one can imagine other roles than precision for this parallel SHH/BMP input4, and their own work3 and that by others5 showed that the key marker of the motor neuron domain, OLIG2, remains sensitive to SHH signaling also at later stages. The limits posed by 8-bit imaging and the inaccuracy of the chosen error approximation are rigorous mathematical facts. They cannot be challenged biologically.

Given that the experimental data cannot be used beyond the 8-bit limit, we developed computational approaches to estimate the gradient variability further away from the source based on available experimental measurements1. This necessitated assumptions, which the authors now question. For one, based on error propagation, we determined the expected gradient variability beyond the 8-bit detection limit for the case that gradients maintain their measured exponential shape across the domain. Secondly, we developed a cell-based simulation framework that allows us to estimate the positional error from measured molecular noise levels. Both the error propagation method and the cell-based simulations showed that gradients remain sufficiently precise such that gradients can, in principle, pattern the NT throughout its development.

Zagorski et al. claim that it is unrealistic to assume that gradients remain exponential. However, they use exponential gradients, rather than their measured gradients, in their own paper2, both in their error propagation method (FitEPM), and also as input when evaluating the potential of opposing gradients in NT patterning via their decoding map, as the reported gradients are too flat and noisy outside the 8-bit limit, also when both gradients are considered simultaneously. As we showed in follow-up work, our results also apply to non-exponential gradients6, and gradient precision is substantially higher still when taking into consideration that morphogens spread at least in 2D, rather than 1D7.

The authors further claim that our statistical approach was flawed because the gradient amplitude and decay length of exponential gradients would be correlated. As the authors only provide a correlation coefficient but no data, it remains unclear what data they base their statement on. Consistent with data of the SHH gradient in the mouse NT8 that we based our analysis on (Pearson’s R = −0.0061, Kendall’s τ = 0.0561), our cell-based model shows that such a correlation, while expected in a deterministic setting, is negligibly small with independent cell-to-cell molecular noise (Fig. 1).

Data from n = 1000 simulated gradients obtained with independent noise between individual cells in all three kinetic parameters, CVp,d,D = 0.3. See1 for method details. Pearson’s R = −0.082 (p = 0.010), Kendall’s τ = −0.053 (p = 0.012).

As we analyzed the processed gradient data that we received from the authors, we noticed that the gradients had been binned from five somite stages and scaled to the same domain length before determining the positional error. This introduces an artificial positional error, which can explain the remaining small deviation in the positional error of gradients and readouts. The authors claim that we overestimate this artificial positional error as the errors may not be fully additive. With the unprocessed data inaccessible to us, the degree of additivity of this artificial error in the authors’ methodology on the biological positional error remains unknown. However, as we show in our article, the positional error of the gradients and the readout are comparable even if the errors are only partially additive, as the difference is rather small.

Zagorski et al. criticize that we did not also correct the readout data. The method section of their Science paper2 mentions the scaling neither for the gradients nor for the readout, and the raw data is not published. We therefore measured the positional error of the dorsal NKX6.1 boundary ourselves. We could reproduce their results only if we bin, but do not scale the domains, from which we concluded that the authors’ readout dataset was likely not scaled. During the writing of our response to their Matters Arising letter, the authors informed us that they did not scale the readout data (A. Kicheva, personal communication; Supplementary Information). With the data and scripts that we have access to we cannot check this. While we remain interested in settling the point, we consider domain scaling a minor issue, given the considerable challenges in detecting gradients and aligning them with their readouts, in particular in pseudostratified epithelia, where nuclei are not perfectly aligned with their apical surface9.

In conclusion, none of the points raised by Zagorski et al. bear relevance to our conclusions. The remaining uncertainties, whilst in principle addressable, cannot be resolved with the data and method details that we have access to. See Supplementary Information for the detailed rebuttal.

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Vetter, R. & Iber, D. Precision of morphogen gradients in neural tube development. Nat. Commun. 13, 1145 (2022).

Article ADS CAS PubMed PubMed Central Google Scholar

Zagorski, M. et al. Decoding of position in the developing neural tube from antiparallel morphogen gradients. Science 356, 1379–1383 (2017).

Article ADS CAS PubMed PubMed Central Google Scholar

Kicheva, A. et al. Coordination of progenitor specification and growth in mouse and chick spinal cord. Science 345, 1254927 (2014).

Article PubMed PubMed Central Google Scholar

Vetter, R. & Iber, D. Morphogen gradients can convey position and time in growing tissues. BioRxiv. https://doi.org/10.1101/2023.08.28.555133 (2023).

Yu, K., McGlynn, S. & Matise, M. P. Floor plate-derived sonic hedgehog regulates glial and ependymal cell fates in the developing spinal cord. Development 140, 1594–1604 (2013).

Article CAS PubMed PubMed Central Google Scholar

Adelmann, J. A., Vetter, R. & Iber, D. Patterning precision under non-linear morphogen decay and molecular noise. eLife 12, e84757 (2023).

Article PubMed PubMed Central Google Scholar

Long, Y., Vetter, R. & Iber, D. 2D effects enhance precision of gradient-based tissue patterning. iScience 26, 107880 (2023).

Article ADS PubMed PubMed Central Google Scholar

Cohen, M. et al. Ptch1 and Gli regulate Shh signalling dynamics via multiple mechanisms. Nat. Commun. 6, 6709 (2015).

Article ADS CAS PubMed Google Scholar

Iber, D. & Vetter, R. Relationship between epithelial organization and morphogen interpretation. Curr. Opin. Genet. Dev. 75, 101916 (2022).

Article CAS PubMed Google Scholar

Download references

Department of Biosystems Science and Engineering, ETH Zürich, Schanzenstrasse 44, 4056, Basel, Switzerland

Roman Vetter & Dagmar Iber

Swiss Institute of Bioinformatics, Schanzenstrasse 44, 4056, Basel, Switzerland

Roman Vetter & Dagmar Iber

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

R.V. and D.I. jointly wrote this manuscript.

Correspondence to Dagmar Iber.

The authors declare no competing interests.

Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

Vetter, R., Iber, D. Reply to: Assessing the precision of morphogen gradients in neural tube development. Nat Commun 15, 930 (2024). https://doi.org/10.1038/s41467-024-45149-7

Download citation

Received: 23 May 2022

Accepted: 15 January 2024

Published: 01 February 2024

DOI: https://doi.org/10.1038/s41467-024-45149-7

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative