Vein fate determined by flow-based but time-delayed integration of network architecture

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 flow-based vascular networks.