Network Modelling of Murine Lymphatic System

Preprint

Article

Network Modelling of Murine Lymphatic System

Altmetrics

Downloads

162

Views

Comments

A peer-reviewed article of this preprint also exists.

Dmitry Grebennikov

Rostislav Savinkov,Ekaterina Zelenova,

Gennady Lobov

Gennady Bocharov^*

Dmitry Grebennikov

Rostislav Savinkov,Ekaterina Zelenova,

Gennady Lobov

Gennady Bocharov^*

This version is not peer-reviewed

Submitted:

17 February 2023

Posted:

21 February 2023

You are already at the latest version

Alerts

Abstract

Animal models of diseases in particular, the mice, are considered to be the cornerstone for translational research in immunology. The aim of the present study is to model the geometry and analyze the network structure of murine lymphatic system (LS). The algorithm for building the graph model of the LS makes use of anatomical data. To identify the edge directions of the graph model, a mass balance approach to lymph dynamics based on Hagen-Poiseuille equation is applied. It is the first study in which a geometric model of the murine LS has been developed and characterized in terms of its structural organization and the lymph transfer function. Our study meets the demand for quantitative mechanistic approaches in the growing field of immunoengineering to utilize or exploit the lymphatic system for immunotherapy.

Keywords:

Subject: Computer Science and Mathematics - Applied Mathematics

1. Introduction

Networks of various nature, e.g., structural, functional, spatial, underlie the dynamics and mechanisms of regulation in live systems ranging from cells to physiological organs and to whole organisms [1]. Network concepts are increasingly used to describe the structure and function of the immune system [2]. The immune system represents an example of a highly complex network of interacting and migration cell populations embedded into a spatially distributed components of the lymphatic system (LS) [3].

Animal models of diseases in particular, the mice, are considered to be the cornerstone for translational research in immunology [4,5]. The research with laboratory mice enabled invaluable insight into mammalian immune systems [6]. Despite numerous advances in understanding the immune system from mouse studies, there exist fundamental differences between mouse and human immune systems [5]. The structural organization of the lymphatic system represents a straightforward example. However, a mathematical network-based characterization of the LS in mice is still unavailable.

The application of graph theory methods to describe the spatial organization of the human LS has been addressed in a number of recent studies (see for a review [7]). In [7], we developed a computational algorithm for representing the anatomy-based and rule-based graphs of the LS in humans. The graph models enabled the analysis of different metrics of complexity of human LS such as spectral radius, clustering coefficient, average path length, number of separators. A similar analysis for the mouse LS remains to be done.

The aim of the present study is to model and analyze the network structure of the murine lymphatic system. The algorithm for building the graph model makes use of anatomical data. To identify the edge directions of the graph model, a mass balance approach to lymph dynamics based on Hagen-Poiseuille equation is applied. Various matrix forms for graph representation are specified. The lymph transfer times between various nodes are estimated. We summarize the properties of the graph model of the murine LS using metrics similar to the human LS graph thus providing a quantitative basis for understanding essential structural differences of the LS between the mice and humans.

2. Anatomy and physiology of murine lymphatic system

Available anatomical and physiological data provide the empirical basis for specifying the network structure of the lymphatic system in mice [8,9]. There are some variations in the number of lymph nodes, i.e. ranging from 22 to 36 as indicated in Table 1. A generalized graph of the murine LS, consisting of 88 nodes and 87 edges in shown in Figure 1. It is developed using anatomical descriptions from [8,9]. The vertices of the graph refer to either the lymph nodes, outlet vertices with out-degree

d e g^{+} = 0

corresponding to the sink into jugular veins, the confluences of lymphatic vessels, or inlet vertices with in-degree

d e g^{-} = 0

corresponding to the collecting lymphatics of various body tissues.

The geometric characteristics of the lymphatic vessels and the baseline parameters of the lymph flow through various parts of the LS network are detailed in Table 1. To set the pressure at the sink nodes

p_{o u t}

, we used the estimate of the murine central venous pressure: 7.4 (5.9-8.9) cm H₂O [15,16]. Lymph viscosity is taken to be equal to 1.81 mPa·s.

The anatomy data enable to specify a simple undirected graph of the murine LS. The adjacency matrix A of the LS graph is shown in Figure 2 (B).

As the LS functions to transport the lymph from the drained tissues to the venous part of cardiovascular system, additional analysis is required to transfer the undirected graph representation to a physiologically meaningful oriented graph of the LS. To generate an oriented graph mode of the LS, we used a combination of experimental studies on fluid dynamics in various parts of the LS in mice ([10,12,13,14]) and computation of the lymph flow through the system in accordance to an overall mass balance using the Hagen-Poiseuille equations.

3. Oriented graph model of murine LS

The graph of the lymphatic system g can be divided in two disconnected subgraphs: collecting the lymph to the left common jugular vein (

g_{l}

) and to the right common jugular vein (

g_{r}

) (Figure 1,Figure 3). The thoracic duct belongs to the subgraph

g_{l}

3.1. Computing the direction of lymph flows

As we aim to reproduce the target flow through thoracic duct

q_{t d} = 10

mL/day (Table 1), we compute the flows in the left subgraph first. The pressure at the sink vertex is set to

p_{o u t} = 7.4

cm H₂O, and the outflow from the left LS subgraph to the left jugular vein

q_{o u t}^{(l)}

is calibrated so that the computed flow through thoracic duct (the edge incident to jugular vein) is equal to

q_{t d}

. The inflows in the collecting lymphatics (inlet vertices with zero in-degree) are given to be the same and equal to

q_{i n}^{(l)} = q_{o u t}^{(l)} / n_{i n}^{(l)}

, where

n_{i n}^{(l)}

is the number of inlet vertices. After obtaining the flows in the left subgraph

g_{l}

, we compute the flows in the right subgraph

g_{r}

by setting the same inflows

q_{i n}^{(r)} = q_{i n}^{(l)}

as in the left one, and the same output pressure

p_{o u t}

. The outlet outflow is given by

q_{o u t}^{(r)} = q_{i n}^{(r)} n_{i n}^{(r)}

The distribution of the steady flows in the graph

g (n, m)

with n vertices and m edges is considered to be governed

by the Hagen-Poiseuille equation:

$q_{i j} = g_{i j} (p_{i} - p_{j}) = \frac{π r_{i j}^{4}}{8 μ l_{i j}} (p_{i} - p_{j}),$

(1)

which links the flow $q_{i j}$ through the edge $e_{i j}$ with the drop of pressure from the tail i to the head j vertices ( $p_{i} - p_{j}$ ) by the conductances $g_{i j}$ that depend on lymph viscosity $μ$ and the radii and the lengths of the edges,
and by the balance of flow through the vertices due to mass conservation:

$\sum_{j \in A (i)} q_{i j} = \{\begin{matrix} q_{i n}, & if i is inlet vertex \\ - q_{o u t}, & if i is outlet vertex \\ 0, & otherwise \end{matrix}$

(2)

where $A (i)$ is a set of vertices adjacent to i.

Using the oriented incidence matrix

M \in R^{n \times m}

and the diagonal conductance matrix

G \in R^{m \times m}

(with elements indexed by edges rather than nodes), one can rewrite the equations (1-2) as

q = - G M^{T} p, M q = - \hat{q},

(3)

where

p \in R^{n}

are the nodal pressures,

\hat{q} \in R^{n}

are the net flows through the vertices and

q \in R^{m}

are the flows through the edges. Hence, we get the linear system to solve for the nodal pressures:

M G M^{T} p = L p = \hat{q},

(4)

where

L = M G M^{T}

is a symmetric weighted Laplacian matrix.

As at the outlet vertex the pressure is known (

p_{o u t}

), we substitute the vector

\hat{p} = {[0, \dots, p_{o u t}, \dots, 0]}^{T}

of zeros with the known pressure at the corresponding index in (4). By shifting

L \hat{p}

to the right hand side, we obtain the system

L_{r e c t} p_{u n k n o w n} = \hat{q} - L \hat{p},

(5)

where

L_{r e c t}

is the matrix L without the column corresponding to the index of the outlet vertex with known pressure. The pseudo-inverse of

L_{r e c t}

provides the vector of unknown pressures:

p_{u n k n o w n} = L_{r e c t}^{+} (\hat{q} - L \hat{p})

An oriented graph of the murine LS resulting from the analysis of the global lymphatic flow balance is shown in Figure 1. It is derived assuming constant diameter of vessels in the LS.

3.2. Matrix representation

To visualize the graph structure of the murine LS, we use the adjacency matrix, which indicates whether a pair of nodes are adjacent or not. For the directed graph, the adjacency matrix is shown in Figure 2 (A). The adjacency matrix of the undirected graph presented in Figure 2 (B) is symmetric. A complementary representation of the graph is provided by the incidence matrix (see Figure 2 (C)). The incidence matrix is different from an adjacency matrix, and it encodes the relation of node-vertex pairs. Finally, the graph Laplacian matrix is displayed in Figure 2 (D). It is related to the degree matrix D and the adjacency matrix A of the graph

L = D - A

, representing an edge-weighted graph.

4. Quantitative characterization of lymph flows through the LS

The estimated values of the lymph flow through various vessels of the murine LS are specified in Figure 3. The upper panel shows the flow intensity in the major section of the LS, draining the left and lower parts of the body. The baseline values vary from about 20 to 420

μ L / h r

. The lower panel characterizes the flow intensity in the minor section of the LS, draining the upper right part of the body. The baseline values vary from about 20 to only 100

μ L / h r

The radius of the thoracic duct is known to be around 300

μ

m, while the radius of the lymphatic vessels afferent to popliteal nodes is around 20-40

μ

m (Table 1). Due to the lack of detailed data on the diameters of all lymphatic vessels, we made three complementary assumptions on the distribution of the radii of the edges of the lymphatic graph, as specified in the following three scenarios:

Scenario 1. All radii in the graph are assumed to be the same, equal to 150 $μ$ m (half of the radius on the thoracic duct).
Scenario 2. Edge radii vary linearly with distance from the outlet vertices (jugular veins) to the inlet vertices. On the thoracic duct, the radius is assumed to be 300 $μ$ m, on the most distant edges (from the hind legs)– 41 $μ$ m. Therefore, on other edges from the inlet vertices, the value of the radius is equal to 41 $μ$ m and increases linearly when approaching the vein. On the subgraph collecting lymph into the right jugular vein, the radii are set symmetrically, equal to the radii in the left subgraph.
Scenario 3. Edge radii are distributed so that the cross-sectional area of incoming and outgoing vessels for each vertex of the graph is preserved. On all inlet edges, the radii are assumed to be the same and are estimated so that the radius on the thoracic duct would be equal to 300 $μ$ m.

The histograms of the vessel radii distribution for the above scenarios are presented in the left column of Figure 4. They clearly indicate that the median value of the vessel lumen decreases as we move from the 1st to the third scenario (from 150

μ m

(Scenario 1), to 115

μ m

(Scenario 2), to 60

μ m

(Scenario 3)).. Note, that the estimated lymph flows shown in Figure 3 refer to Scenario 1.

4.1. Lymph transfer rates between lymph nodes

The key characteristic of the LS function is the rate of transfer of fluid through the system. Using the developed graph model, we estimated the lymph flow rates and the transition times between the lymph nodes. The computational results for three scenarios are summarized in Figure 4. The central column provides the estimates of the flow velocity. The median values turn out to be the smallest one for the scenario of a uniform vessel diameters and the largest one under assumption of conservation of the cross-sectional area at the vessel junctions. In particular, it increases from about 66

μ m / s

(Scenario 1), to 242

μ m / s

(Scenario 2), to 409

μ m / s

(Scenario 3). In addition, it is predicted to by practically homogeneous across the LS in the third case. The respective median transfer times between neighboring edges increase from about 19 seconds to 95 second with the longest transfer times from 2 minutes 25 seconds to 15 minutes 15 seconds when we compare the third and the first scenarios of the vessel geometry.

4.2. Sensitivity to variations in vessels diameter

To analyze the sensitivity of the system to the diameter of the vessels, for all three scenarios, situations were simulated when the diameter of the vessel represented by the corresponding branch of the graph was varied by +10% and -10% from the initial one. The effect of change of the vessels diameters was expressed in terms of the histograms of the relative change of pressure in the vessels of the LS as shown in Figure 5. In all three scenarios, the reduction of the radii led expectedly to a pressure increase, while the increase of the vessel radii had an opposite effect. The degree of variation was smallest for Scenario 1 and was largest for Scenario 2. The 10% diameter variation led to les than 1% change in the pressure for most of the vessels.

5. Topological properties of the LS graph

Following our previous study of the human LS [7], we characterize the topological properties of the murine LS graph model using some fundamental metrics. Let

G = (n, m)

be the graph with n nodes and m edges, respectively. Consider the following characteristics:

the number of input nodes $N_{i n p}$ , i.e., the number of nodes with degree 1 and out-degree 0;
maximum degree of graph $Δ_{G}$ , i.e., the maximum degree of its vertices;
girth of the graph g, which is the length of the shortest (undirected) cycle in the graph;
diameter, i.e., the longest geodesic distance (in other terms, maximum eccentricity of any vertex),

$D = max_{v \in V} ϵ (v) = max_{v \in V} max_{u \in V} d (u, v),$

(6)

where $d (u, v)$ is the geodesic distance (shortest directed path connecting vertices u and v), $ϵ (v)$ is the eccentricity of vertex v;
radius of the graph (minimum eccentricity of any vertex),

$r = min_{v \in V} ϵ (v) = min_{v \in V} max_{u \in V} d (u, v) .$

(7)
average path length (mean geodesic distance),

$l_{G} = \frac{1}{n (n - 1)} \sum_{u, v \in V, u \neq v} d (u, v) .$

(8)
the energy and the spectral radius of the graph are defined as follows,

$E n (A) = \sum_{j = 1}^{n} ∣ λ_{j} ∣, ρ (A) = max {| λ_{j} |},$

(9)

where $λ_{j}$ stand for the eigenvalues of the adjacency matrix A of the graph;
edge density of the graph, i.e., the number of edges divided by the number of all possible edges,

$ρ_{d} = \frac{m}{n (n - 1)} .$

(10)
The clustering coefficient C (transitivity) measures the probability that two neighbors of a vertex are connected. It can be computed as function of adjacency matrix A:

$C (A) = \frac{\sum_{i = 1, j = 1, k = 1}^{n, n, n} a_{i j} \cdot a_{j k} \cdot a_{k i}}{\sum_{i = 1}^{n} ((\sum_{j = 1}^{n} a_{i j}) \cdot ((\sum_{j = 1}^{n} a_{i j}) - 1))} .$

(11)
Number of separators $n_{s e p}$ , i.e., the vertices removal of which disconnects the graph.
Topological diversity of the vertices as a function of the Shannon entropy associated with flow rates through the incident edges,

$D_{f l o w} (v_{i}) = \frac{H (v_{i})}{log (k)} = \frac{- \sum_{j = 1}^{k} p_{i j} log (p_{i j})}{log (k)}, p_{i j} = \frac{| Q_{i j} |}{\sum_{j = 1}^{k} | Q_{i j |}},$

(12)

where k is the number of $v_{i}$ ’s incident edges and $p_{i j}$ is the proportion of the flow between the adjacent $v_{i}$ and $v_{j}$ to the total flow through the edges involving $v_{i}$ . The flow diversity is defined similar to the definition of network diversity in [?].

To analyze the robustness of the mouse LS to damage, we sequentially removed individual nodes of the graph and checked how many source vertices remain connected to the sink vertex into the circulatory system. The subgraphs of the LS and whole graph were analyzed. Accordingly, the robustness of the graph was estimated as the arithmetic mean of the ratio of the number of source vertices that retained the connection to the sink to their total number in the graph/subgraphs.

The summary topological properties of the murine LS graph model are presented in Table 1. The characteristics of the whole LS graph and the two subgraphs representing the LS parts draining the right∪low and left parts of the body are specified.

6. Conclusions

In this work, we have developed a graph model of the LS network in mice. It is the first study in which a geometric model of the murine LS has been developed and characterized in terms of its structural organization and the lymph transfer function. The study complements our previous analysis of the human lymphatic system [7]. The developed graph model of the LS in mice provides a computational tool for studying the spatial aspects of the immune system functioning. It goes in line with recent advances in experimental techniques to characterize the whole-body dynamics of systemic infections in experimental mice [17].

The estimated parameters of the LS function in terms of the lymph flow rate and transfer time between various parts of the mouse body can be used on compartmental modelling for evaluation of the pharmacokinetic characteristics of drugs and adoptive cell therapies in advance of experiments. A remarkable example of the use of our recently developed graph model of human LS [7] is given in the study of how the topology of the lymphatic network affects the time for an immune search through the lymphatic network to be completed [18].

Further development of the presented graph model of the LS can be visioned to proceed along three lines:

considering the biomechanics of lymphatic pumping through a chain of lymphangions and lymph nodes,
coupling the LS model with the cardiovascular system,
integration with multi-physics models of the immune system.

Overall, our study meets the demand for quantitative rigorous approaches in the growing field of immunoengineering to utilize or exploit the lymphatic system for immunotherapy first in experimental animals and then to cure human immune-dependent diseases [19,20,21].

Author Contributions

Conceptualization, D.G., R.S., G.L., and G.B.; methodology, D.G., R.S., G.B.; software, D.G., R.S., E.Z.,; validation, D.G., R.S., E.Z., G.L., and G.B.; formal analysis, D.G., R.S., E.Z.,; investigation, D.G., R.S., E.Z., G.L., and G.B.; data curation, G.L. ; writing—original draft preparation, D.G., R.S., and G.B.; writing—review and editing, D.G., R.S., E.Z., G.L., and G.B.; visualization, D.G. and R.S.; supervision, G.B.; funding acquisition, G.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Russian Science Foundation grant number 18-11-00171. R.S., D.G. and G.B. were partly supported by the Russian Foundation for Basic Research grant number 20-01-00352 and by Moscow Center for Fundamental and Applied Mathematics (agreement with the Ministry of Education and Science of the Russian Federation No. 075-15-2022-286).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:

LS	Lymphatic system
LN	Lymph node

References

Ivanov, P. The New Field of Network Physiology: Building the Human Physiolome. Front. Netw. Physiol. 2021, 1. [Google Scholar] [CrossRef]
Hao, S.; Koon-Kiu, Y.; Liang, D.; Chenxi, Q.; Hongbo, C.; Jiyang, Y. Network Approaches for Dissecting the Immune System. iScience 2020, 23. [Google Scholar] [CrossRef]
Randolph, G.; Ivanov, S.; Zinselmeyer, B.; Scallan, J. The Lymphatic System: Integral Roles in Immunity. Annu Rev Immunol 2017, 35, 31–52. [Google Scholar] [CrossRef] [PubMed]
Ernst, P.; Carvunis, A. Of mice, men and immunity: a case for evolutionary systems biology. Nat Immunol 2018, 19, 421–425. [Google Scholar] [CrossRef] [PubMed]
Tao, L.; Reese, T. Making Mouse Models That Reflect Human Immune Responses. Trends Immunol 2017, 38, 181–193. [Google Scholar] [CrossRef] [PubMed]
Graham, A. Naturalizing mouse models for immunology. Nat Immunol 2021, 22, 111–117. [Google Scholar] [CrossRef] [PubMed]
Savinkov, R.; Grebennikov, D.; Puchkova, D.; Chereshnev, V.; Sazonov, I.; Bocharov, G. Graph Theory for Modeling and Analysis of the Human Lymphatic System. Mathematics 2020, 8, 2236. [Google Scholar] [CrossRef]
Yoshitsugu, K.; Makoto, S.; Yann-Ching, H.; Norio, K. THE LYMPH SYSTEM IN MICE. Japanese Journal of Veterinary Research 1964, 12, 69–78. [Google Scholar] [CrossRef]
Van den Broeck, W.; Derore, A.; Simoens, P. Anatomy and nomenclature of murine lymph nodes: Descriptive study and nomenclatory standardization in BALB/cAnNCrl mice. Journal of Immunological Methods 2006, 312, 12–19. [Google Scholar] [CrossRef] [PubMed]
Takeda, K.; Mori, S.; Kodama, T. Study of fluid dynamics reveals direct communications between lymphatic vessels and venous blood vessels at lymph nodes of mice. Journal of Immunological Methods 2017, 445, 1–9. [Google Scholar] [CrossRef] [PubMed]
Economopoulos, V.; Noad, J.C.; Krishnamoorthy, S.; Rutt, B.K.; Foster, P.J. Comparing the MRI Appearance of the Lymph Nodes and Spleen in Wild-Type and Immuno-Deficient Mouse Strains. PLoS ONE 2011, 6, e27508. [Google Scholar] [CrossRef]
Lindena, J.; Küpper, W.; Trautschold †, I. Catalytic Enzyme Activity Concentration in Thoracic Duct, Liver, and Intestinal Lymph of the Dog, the Rabbit, the Rat and the Mouse. Approach to a Quantitative Diagnostic Enzymology, II. Communication. Clinical Chemistry and Laboratory Medicine 1986, 24. [Google Scholar] [CrossRef]
Blatter, C.; Meijer, E.F.J.; Nam, A.S.; Jones, D.; Bouma, B.E.; Padera, T.P.; Vakoc, B.J. In vivo label-free measurement of lymph flow velocity and volumetric flow rates using Doppler optical coherence tomography. Scientific Reports 2016, 6, 29035. [Google Scholar] [CrossRef] [PubMed]
Bouta, E.M.; Wood, R.W.; Brown, E.B.; Rahimi, H.; Ritchlin, C.T.; Schwarz, E.M. In vivo quantification of lymph viscosity and pressure in lymphatic vessels and draining lymph nodes of arthritic joints in mice: Lymph viscosity and pressure in arthritic joints. The Journal of Physiology 2014, 592, 1213–1223. [Google Scholar] [CrossRef] [PubMed]
Nielsen, J.M.; Kristiansen, S.B.; Ringgaard, S.; Nielsen, T.T.; Flyvbjerg, A.; Redington, A.N.; Bøtker, H.E. Left ventricular volume measurement in mice by conductance catheter: evaluation and optimization of calibration. American Journal of Physiology-Heart and Circulatory Physiology 2007, 293, H534–H540. [Google Scholar] [CrossRef] [PubMed]
Xie, C.; Wei, W.; Zhang, T.; Dirsch, O.; Dahmen, U. Monitoring of Systemic and Hepatic Hemodynamic Parameters in Mice. Journal of Visualized Experiments, 2014; 5195. [Google Scholar] [CrossRef]
Wen, Y.; Xu, H.; Wan, W.; Shang, W.; Jin, R.; Zhou, F.; Mei, H.; Wang, J.; Xiao, G.; Chen, H.; Wu, X.; Zhang, L. Visualizing lymphocytic choriomeningitis virus infection in cells and living mice. iScience 2022, 25, 105090. [Google Scholar] [CrossRef] [PubMed]
Ferdous, J.; Fricke, G.; Moses, M. Modeling Immune Search Through the Lymphatic Network. In: Eds. Dorigo, Marco et al. Swarm Intelligence. Lecture Notes in Computer Science. ANTS 2022, 13491, 332–340. [Google Scholar] [CrossRef]
Thomas, S.; Rohner, N.; Edwards, E. Implications of Lymphatic Transport to Lymph Nodes in Immunity and Immunotherapy. Annu Rev Biomed Eng 2016, 18, 207–33. [Google Scholar] [CrossRef] [PubMed]
Maisel, K.; Stella Sasso, M.; Potin, L.; Swartz, M.A.
O’Melia, M.; Lund, A.; Thomas, S. The Biophysics of Lymphatic Transport: Engineering Tools and Immunological Consequences. iScience 2019, 22, 28–43. [Google Scholar] [CrossRef] [PubMed]

Figure 1. Directed graph of the murine lymphatic system based on the anatomic data with 88 vertices and 87 edges. The vertices of the graph belong to four groups as detailed in the legend box: (a) lymphatic nodes (large blue), (b) outlet vertices with out-degree

d e g^{+} = 0

corresponding to the sink into jugular veins (red), (c) connectors, i.e., the confluences of lymphatic vessels (light blue), (d) inlet vertices with in-degree

d e g^{-} = 0

corresponding to the collecting lymphatics of various body tissues (orange). The vertex ids and edge ids are enumerated arbitrarily for correspondence with the matrix representation of the graph in Figure 2.

d e g^{+} = 0

corresponding to the sink into jugular veins (red), (c) connectors, i.e., the confluences of lymphatic vessels (light blue), (d) inlet vertices with in-degree

d e g^{-} = 0

Figure 2. Matrix representation of the anatomy-based graph of murine lymphatic system presented on Figure 1. (A, B) Adjacency matrix A (for directed and undirected graph). (C) Incidence matrix, M. (D) Weighted Laplacian matrix,

L = M G M^{T}

, normalized by the maximum matrix element. The conductance matrix G for constant vessel diameters is used for the illustration.

L = M G M^{T}

, normalized by the maximum matrix element. The conductance matrix G for constant vessel diameters is used for the illustration.

Figure 3. The distribution of the flows in the lymphatic vessels for Scenario 1 (constant vessel diameters). Flows are shown in the subgraphs of the lymphatics collecting lymph into the left (A) and into the right (B) jugular vein. Colors and the arrows indicate the heads and the tails of the edges correspondingly. Edges are sorted by their distance from the sink.

Figure 4. The distribution of the velocities and transient transfer times of lymph flow in the vessels for the three scenarios of vessel radii distribution. At the bottom, the distribution of the edge lengths is shown which is the same for all scenarios.

Figure 5. The distributions of the changes of pressures in the lymphatic vessels of the LS (vertices of the graph) after variation of the length diameters (± 10%) for Scenario 1 (A,B), Scenario 2 (C,D) and Scenario 3 (E,F).

Table 1. Physiological and anatomical properties of the murine lymphatic system.

Property	Value (range)	References
Lymph nodes:
Number	22 28-36	(BALB/cAnNCr) [9] (DD/NIH) [8]
Diameter	1-2.3 mm 1-17.3 mm	(C57Bl/6J, Nude, CB-17 SCID) [10,11] (DD/NIH) [8]
Thoracic duct:
Radius	300 $μ$ m	[12]
Flow	417-1250 $μ$ L/h (10-30 mL/day)	(10 mL/day for immobilized mice, 30 mL/day after movements) [12]
Velocity	25-75 mm/min (410-1228 $μ$ m/s)	[12]
Vessels afferent to popliteal nodes:
Radius	20-40 $μ$ m	[12,13]
Flow	0.3-3.4 $μ$ L/h	(mean flow = 0.3 $μ$ L/h [13]) [12,13,14]
Velocity	37-186 $μ$ m/s	(mean velocity = 53±16 $μ$ m/s [13]) [12,13,14]
Collecting lymphatics in hind limbs:
Velocity	50-100 $μ$ m/s	[10]
Collecting lymphatics in ears:
Velocity	0-400 $μ$ m/s	[13]
Collecting lymphatics in the tail:
Velocity	4.2 $μ$ m/s	[14]

Table 2. Summary statistics for the anatomy-based graph of murine lymphatic system.

Property	Whole graph, g	Left subgraph, $g_{l}$	Right subgraph, $g_{r}$
$G (n, m)$	$g (88, 87)$	$g_{l} (61, 60)$	$g_{r} (27, 27)$
Number of inlet vertices	52	36	16
Maximum degree	5	5	5
Girth	3	0	3
Diameter, directed (undirected)	7 (11)	7 (11)	4 (7)
Radius	4	6	4
Average path length, dir. (undir.)	2.5 (5.3)	2.7 (5.5)	1.9 (3.9)
Energy	95.1	65.6	29.5
Spectral radius	2.93	2.9	2.93
Edge density	0.0114	0.0164	0.0385
Clustering coefficient	0.019	0	0.059
Number of separators	-	25	11
Robustness	0.917	0.863	0.825
Average topological flow diversity,
- scenario 1: - scenario 2: - scenario 3:	0.8252 0.8242 0.8242	0.8028 0.8028 0.8028	0.9039 0.9015 0.9023
Number of LNs	27	19	8

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer

Network Modelling of Murine Lymphatic System

Abstract

1. Introduction

2. Anatomy and physiology of murine lymphatic system

3. Oriented graph model of murine LS

3.1. Computing the direction of lymph flows

3.2. Matrix representation

4. Quantitative characterization of lymph flows through the LS

4.1. Lymph transfer rates between lymph nodes

4.2. Sensitivity to variations in vessels diameter

5. Topological properties of the LS graph

6. Conclusions

Author Contributions

Funding

Conflicts of Interest

Abbreviations

References

MDPI Initiatives

Important Links

Subscribe