Vein fate determined by flow-based but time-delayed integration of network architecture
Abstract
Veins in vascular networks, such as in blood vasculature or leaf networks, continuously reorganize, grow or shrink, to minimize energy dissipation. Flow shear stress on vein walls has been set forth as the local driver for a vein's continuous adaptation. Yet, shear feedback alone cannot account for the observed diversity of vein dynamics - a puzzle made harder by scarce spatiotemporal data. Here, we resolve network-wide vein dynamics and shear rate during spontaneous reorganization in the prototypical vascular networks of Physarum polycephalum. Our experiments reveal a plethora of vein dynamics (stable, growing, shrinking) where the role of shear is ambiguous. Quantitative analysis of our data reveals that (a) shear rate indeed feeds back on vein radius, yet, with a time delay of 1-3 min. Further, we reconcile the experimentally observed disparate vein fates by developing a model for vein adaptation within a network and accounting for the observed time delay. The model reveals that (b) vein fate is determined by parameters - local pressure or relative vein resistance - which integrate the entire network's architecture, as they result from global conservation of fluid volume. Finally, we observe avalanches of network reorganization events that cause entire clusters of veins to vanish. Such avalanches are consistent with network architecture integrating parameters governing vein fate as vein connections continuously change. As the network architecture integrating parameters intrinsically arise from laminar fluid flow in veins, we expect our findings to play a role across ow-based vascular networks.
Data availability
Original microscopic images of all the specimens used for this study are available as movies in MP4 format.Data sharing plan:All data used to generate the figures and the custom written matlab codes will be available on the open access data repository platform mediaTUM if the paper is accepted and which will correspond to the final versions of the figures.
Article and author information
Author details
Funding
MRSEC Program of the National Science Foundation (Award Number DMR- 1420073)
- Sophie Marbach
Max Planck Society
- Karen Alim
European Research Council under the European Union's Horizon 2020 research and innovation programme (grant agreement No. 947630,FlowMem)
- Karen Alim
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Copyright
© 2023, Marbach et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 916
- views
-
- 189
- downloads
-
- 14
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
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)
Further reading
-
- Physics of Living Systems
Bacterial biofilms are communities of bacteria usually attached to solid strata and often differentiated into complex structures. Communication across biofilms has been shown to involve chemical signaling and, more recently, electrical signaling in Gram-positive biofilms. We report for the first time, community-level synchronized membrane potential dynamics in three-dimensional Escherichia coli biofilms. Two hyperpolarization events are observed in response to light stress. The first requires mechanically sensitive ion channels (MscK, MscL, and MscS) and the second needs the Kch-potassium channel. The channels mediated both local spiking of single E. coli biofilms and long-range coordinated electrical signaling in E. coli biofilms. The electrical phenomena are explained using Hodgkin-Huxley and 3D fire-diffuse-fire agent-based models. These data demonstrate that electrical wavefronts based on potassium ions are a mechanism by which signaling occurs in Gram-negative biofilms and as such may represent a conserved mechanism for communication across biofilms.
-
- Cell Biology
- Physics of Living Systems
The Gillespie algorithm is commonly used to simulate and analyze complex chemical reaction networks. Here, we leverage recent breakthroughs in deep learning to develop a fully differentiable variant of the Gillespie algorithm. The differentiable Gillespie algorithm (DGA) approximates discontinuous operations in the exact Gillespie algorithm using smooth functions, allowing for the calculation of gradients using backpropagation. The DGA can be used to quickly and accurately learn kinetic parameters using gradient descent and design biochemical networks with desired properties. As an illustration, we apply the DGA to study stochastic models of gene promoters. We show that the DGA can be used to: (1) successfully learn kinetic parameters from experimental measurements of mRNA expression levels from two distinct Escherichia coli promoters and (2) design nonequilibrium promoter architectures with desired input–output relationships. These examples illustrate the utility of the DGA for analyzing stochastic chemical kinetics, including a wide variety of problems of interest to synthetic and systems biology.