Resilience of three-dimensional sinusoidal networks in liver tissue

Can three-dimensional, microvasculature networks still ensure blood supply if individual links fail? We address this question in the sinusoidal network, a plexus-like microvasculature network, which transports nutrient-rich blood to every hepatocyte in liver tissue, by building on recent advances in high-resolution imaging and digital reconstruction of adult mice liver tissue. We find that the topology of the three-dimensional sinusoidal network reflects its two design requirements of a space-filling network that connects all hepatocytes, while using shortest transport routes: sinusoidal networks are sub-graphs of the Delaunay graph of their set of branching points, and also contain the corresponding minimum spanning tree, both to good approximation. To overcome the spatial limitations of experimental samples and generate arbitrarily-sized networks, we developed a network generation algorithm that reproduces the statistical features of 0.3-mm-sized samples of sinusoidal networks, using multi-objective optimization for node degree and edge length distribution. Nematic order in these simulated networks implies anisotropic transport properties, characterized by an empirical linear relation between a nematic order parameter and the anisotropy of the permeability tensor. Under the assumption that all sinusoid tubes have a constant and equal flow resistance, we predict that the distribution of currents in the network is very inhomogeneous, with a small number of edges carrying a substantial part of the flow—a feature known for hierarchical networks, but unexpected for plexus-like networks. We quantify network resilience in terms of a permeability-at-risk, i.e., permeability as function of the fraction of removed edges. We find that sinusoidal networks are resilient to random removal of edges, but vulnerable to the removal of high-current edges. Our findings suggest the existence of a mechanism counteracting flow inhomogeneity to balance metabolic load on the liver.

[1]  Nancy K. Drew,et al.  Metrics for Assessing Cytoskeletal Orientational Correlations and Consistency , 2015, PLoS Comput. Biol..

[2]  Felix Kramer,et al.  How to Pare a Pair: Topology Control and Pruning in Intertwined Complex Networks , 2019 .

[3]  Roux,et al.  Force Distributions in Dense Two-Dimensional Granular Systems. , 1996, Physical review letters.

[4]  Tobias Preußer,et al.  Spatio-Temporal Simulation of First Pass Drug Perfusion in the Liver , 2014, PLoS Comput. Biol..

[5]  A. Pries,et al.  Blood flow in microvascular networks. Experiments and simulation. , 1990, Circulation research.

[6]  D. R. Anderson,et al.  Effect of CO2 on intracellular pH and contraction of retinal capillary pericytes. , 1997, Investigative ophthalmology & visual science.

[7]  T. Lubensky,et al.  Principles of condensed matter physics , 1995 .

[8]  P. Arratia,et al.  Rheology of human blood plasma: viscoelastic versus Newtonian behavior. , 2013, Physical review letters.

[9]  3D spatially-resolved geometrical and functional models of human liver tissue reveal new aspects of NAFLD progression , 2019, Nature Medicine.

[10]  Henrik Ronellenfitsch,et al.  Limits of multifunctionality in tunable networks , 2018, Proceedings of the National Academy of Sciences.

[11]  E. Merrill,et al.  Viscosity of human blood: transition from Newtonian to non-Newtonian. , 1967, Journal of applied physiology.

[12]  Felix Kramer,et al.  How to Pare a Pair: Topology Control and Pruning in Intertwined Complex Networks , 2019 .

[13]  Physics of Lumen growth , 2018, bioRxiv.

[14]  J. Happel,et al.  Low Reynolds number hydrodynamics: with special applications to particulate media , 1973 .

[15]  J. S. Andrade,et al.  Modeling river delta formation , 2007, Proceedings of the National Academy of Sciences.

[16]  Karen Alim,et al.  Flow rate of transport network controls uniform metabolite supply to tissue , 2017, Journal of The Royal Society Interface.

[17]  R. Glenny,et al.  Pulmonary blood flow remains fractal down to the level of gas exchange. , 2000, Journal of applied physiology.

[18]  J. Tuszynski,et al.  Building a 3D Virtual Liver: Methods for Simulating Blood Flow and Hepatic Clearance on 3D Structures , 2016, PloS one.

[19]  Y. Kalaidzidis,et al.  A versatile pipeline for the multi-scale digital reconstruction and quantitative analysis of 3D tissue architecture , 2015, eLife.

[20]  E. Katifori,et al.  The Phenotypes of Fluctuating Flow: Development of Distribution Networks in Biology and the Trade-off between Efficiency, Cost, and Resilience , 2017 .

[21]  Sidney Redner,et al.  Fractal and Multifractal Scaling of Electrical Conduction in Random Resistor Networks , 2008 .

[22]  M. Mosharaf-Dehkordi.  A fully coupled porous media and channels flow approach for simulation of blood and bile flow through the liver lobules , 2019, Computer methods in biomechanics and biomedical engineering.

[23]  Thomas M A Fink,et al.  Easily repairable networks: reconnecting nodes after damage. , 2014, Physical review letters.

[24]  Rémi Louf,et al.  Balancing building and maintenance costs in growing transport networks. , 2016, Physical review. E.

[25]  P. Segers,et al.  Perfusion characteristics of the human hepatic microcirculation based on three-dimensional reconstructions and computational fluid dynamic analysis. , 2012, Journal of biomechanical engineering.

[26]  Henrik Ronellenfitsch,et al.  Global Optimization, Local Adaptation, and the Role of Growth in Distribution Networks. , 2016, Physical review letters.

[27]  P. Gennes,et al.  The physics of liquid crystals , 1974 .

[28]  A. Perna,et al.  Animal transportation networks , 2014, Journal of The Royal Society Interface.

[29]  Edward A. Bender,et al.  The Asymptotic Number of Labeled Graphs with Given Degree Sequences , 1978, J. Comb. Theory A.

[30]  Z. Kmieć,et al.  Introduction — Morphology of the Liver Lobule , 2001 .

[31]  Eleni Katifori,et al.  Extracting Hidden Hierarchies in 3D Distribution Networks , 2014, 1410.3951.

[32]  Chwee Teck Lim,et al.  Topological defects in epithelia govern cell death and extrusion , 2017, Nature.

[33]  A. Bonfiglio,et al.  Mathematical modeling of the circulation in the liver lobule. , 2010, Journal of biomechanical engineering.

[34]  Jörn Dunkel,et al.  Stochastic cycle selection in active flow networks , 2016, Proceedings of the National Academy of Sciences.

[35]  A. Pries,et al.  Red cell distribution at microvascular bifurcations. , 1989, Microvascular research.

[36]  Eleni Katifori,et al.  Damage and fluctuations induce loops in optimal transport networks. , 2009, Physical review letters.

[37]  H. Holzhütter,et al.  Modeling function–perfusion behavior in liver lobules including tissue, blood, glucose, lactate and glycogen by use of a coupled two-scale PDE–ODE approach , 2014, Biomechanics and Modeling in Mechanobiology.

[38]  Dan Hu,et al.  Adaptation and optimization of biological transport networks. , 2013, Physical review letters.

[39]  Marco Zimmerling,et al.  Renormalization group theory for percolation in time-varying networks , 2018, Scientific Reports.

[40]  Shyr-Shea Chang,et al.  Microvscular networks with uniform flow. , 2019, Journal of theoretical biology.

[41]  L Mahadevan,et al.  Pressure-driven occlusive flow of a confined red blood cell. , 2016, Soft matter.

[42]  M. Brenner,et al.  Robust Increase in Supply by Vessel Dilation in Globally Coupled Microvasculature. , 2019, Physical review letters.

[43]  A. Tero,et al.  Rules for Biologically Inspired Adaptive Network Design , 2010, Science.

[44]  Thomas G. Fai,et al.  Image-based model of the spectrin cytoskeleton for red blood cell simulation , 2017, PLoS Comput. Biol..

[45]  S. Safran,et al.  How cells feel their substrate: spontaneous symmetry breaking of active surface stresses , 2012 .

[46]  Franck Nicoud,et al.  Red cells’ dynamic morphologies govern blood shear thinning under microcirculatory flow conditions , 2016, Proceedings of the National Academy of Sciences.

[47]  Aric Hagberg,et al.  Exploring Network Structure, Dynamics, and Function using NetworkX , 2008 .

[48]  M. Brenner,et al.  Random network peristalsis in Physarum polycephalum organizes fluid flows across an individual , 2013, Proceedings of the National Academy of Sciences.

[49]  Hernan Morales-Navarrete,et al.  Automatic recognition and characterization of different non-parenchymal cells in liver tissue , 2016, 2016 IEEE 13th International Symposium on Biomedical Imaging (ISBI).

[50]  Hernan Morales-Navarrete,et al.  A Predictive 3D Multi-Scale Model of Biliary Fluid Dynamics in the Liver Lobule. , 2017, Cell systems.

[51]  J. Boyer,et al.  Bile formation and secretion. , 2013, Comprehensive Physiology.

[52]  Gou-Jen Wang,et al.  A lab-on-a-chip capillary network for red blood cell hydrodynamics , 2010 .

[53]  Y. Kalaidzidis,et al.  Liquid-crystal organization of liver tissue , 2018, bioRxiv.

[54]  Walter de Back,et al.  Quantification of nematic cell polarity in three-dimensional tissues , 2019, PLoS Comput. Biol..

[55]  J B Bassingthwaighte,et al.  Regional Distribution of Diffusible Tracers and Carbonized Microspheres in the Left Ventricle of Isolated Dog Hearts , 1973, Circulation research.

[56]  S. Whitaker.  Flow in porous media I: A theoretical derivation of Darcy's law , 1986 .

[57]  G. Dubini,et al.  Microcirculation in the murine liver: a computational fluid dynamic model based on 3D reconstruction from in vivo microscopy. , 2017, Journal of biomechanics.

[58]  Robin Fåhræus,et al.  THE VISCOSITY OF THE BLOOD IN NARROW CAPILLARY TUBES , 1931 .

[59]  S. Kirkpatrick.  Percolation and Conduction , 1973 .

[60]  R. Guthke,et al.  Integrated metabolic spatial‐temporal model for the prediction of ammonia detoxification during liver damage and regeneration , 2014, Hepatology.