Development of a model of a multi-lymphangion lymphatic vessel incorporating realistic and measured parameter values

The Lymphatic Vascular System: Does Nonuniform Lymphangion Length Limit Flow-Rate?

A previously developed model of a lymphatic vessel as a chain of lymphangions was investigated to determine whether lymphangions of unequal length reduce pumping relative to a similar chain of equal-length ones. The model incorporates passive elastic and active contractile properties taken from ex vivo measurements, and intravascular lymphatic valves as transvalvular pressure-dependent resistances to flow with hysteresis and transmural pressure-dependent bias to the open state as observed experimentally. Coordination of lymphangion contractions is managed by marrying an autonomous transmural pressure-dependent pacemaker for each lymphangion with bidirectional transmission of activation signals between lymphangions, qualitatively matching empirical observations. With eight lymphangions as used here and many nonlinear constraints, the model is capable of complex outcomes. The expected flow-rate advantage conferred by longer lymphangions everywhere was confirmed. However, the anticipated advantage of uniform lymphangions over those of unequal length, compared in chains of equal overall length, was not found. A wide variety of dynamical outcomes was observed, with the most powerful determinant being the adverse pressure difference, rather than the arrangement of long and short lymphangions. This work suggests that the wide variation in lymphangion length which is commonly observed in collecting lymphatic vessels does not confer disadvantage in pumping lymph.

Potentiating glymphatic drainage minimizes post-traumatic cerebral oedema

Cerebral oedema is associated with morbidity and mortality after traumatic brain injury (TBI)1. Noradrenaline levels are increased after TBI2-4, and the amplitude of the increase in noradrenaline predicts both the extent of injury5 and the likelihood of mortality6. Glymphatic impairment is both a feature of and a contributor to brain injury7,8, but its relationship with the injury-associated surge in noradrenaline is unclear. Here we report that acute post-traumatic oedema results from a suppression of glymphatic and lymphatic fluid flow that occurs in response to excessive systemic release of noradrenaline. This post-TBI adrenergic storm was associated with reduced contractility of cervical lymphatic vessels, consistent with diminished return of glymphatic and lymphatic fluid to the systemic circulation. Accordingly, pan-adrenergic receptor inhibition normalized central venous pressure and partly restored glymphatic and cervical lymphatic flow in a mouse model of TBI, and these actions led to substantially reduced brain oedema and improved functional outcomes. Furthermore, post-traumatic inhibition of adrenergic signalling boosted lymphatic export of cellular debris from the traumatic lesion, substantially reducing secondary inflammation and accumulation of phosphorylated tau. These observations suggest that targeting the noradrenergic control of central glymphatic flow may offer a therapeutic approach for treating acute TBI.

An Enhanced 3D Model of Intravascular Lymphatic Valves to Assess Leaflet Apposition and Transvalvular Differences in Wall Distensibility

Lymphatic valves operate in a fluid-dynamically viscous environment that has little in common with that of cardiac valves, and accordingly have a different, axially lengthened, shape. A previously developed 3D fluid/structure interaction model of a lymphatic valve was extended to allow the simulation of stages of valve closure after the leaflets come together. This required that the numerical leaflet be prevented from passing into space occupied by the similar other leaflet. The resulting large deflections of the leaflet and lesser deflections of the rest of the valve were mapped as functions of the transvalvular pressure. In a second new development, the model was reconstructed to allow the vessel wall to have different material properties on either side of where the valve leaflet inserts into the wall. As part of this, a new pre-processing scheme was devised which allows easier construction of models with modified valve dimensions, and techniques for successfully interfacing the CAD software to the FE software are described. A two-fold change in wall properties either side of the leaflet made relatively little difference to valve operation apart from affecting the degree of sinus distension during valve closure. However, the numerically permitted strains were modest (<14%), and did not allow examination of the large-scale highly nonlinear elastic properties exhibited by real lymphatic vessels. A small series of murine popliteal, mesenteric, and inguinal-axillary lymphatic vessel segments containing a valve were experimentally investigated ex vivo. The pressure–diameter curves measured just upstream and just downstream of the valve were parameterised by computing the difference in tubular distensibility at three values of transmural pressure. In the popliteal and mesenteric segments, it was found that the distensibility was usually greater just downstream, i.e., in the sinus region, than upstream, at low and intermediate transmural pressure. However, there was wide variation in the extent of difference, and possible reasons for this are discussed.

A multiscale sliding filament model of lymphatic muscle pumping

The lymphatics maintain fluid balance by returning interstitial fluid to veins via contraction/compression of vessel segments with check valves. Disruption of lymphatic pumping can result in a condition called lymphedema with interstitial fluid accumulation. Lymphedema treatments are often ineffective, which is partially attributable to insufficient understanding of specialized lymphatic muscle lining the vessels. This muscle exhibits cardiac-like phasic contractions and smooth muscle-like tonic contractions to generate and regulate flow. To understand the relationship between this sub-cellular contractile machinery and organ-level pumping, we have developed a multiscale computational model of phasic and tonic contractions in lymphatic muscle and coupled it to a lymphangion pumping model. Our model uses the sliding filament model (Huxley in Prog Biophys Biophys Chem 7:255-318, 1957) and its adaptation for smooth muscle (Mijailovich in Biophys J 79(5):2667-2681, 2000). Multiple structural arrangements of contractile components and viscoelastic elements were trialed but only one provided physiologic results. We then coupled this model with our previous lumped parameter model of the lymphangion to relate results to experiments. We show that the model produces similar pressure, diameter, and flow tracings to experiments on rat mesenteric lymphatics. This model provides the first estimates of lymphatic muscle contraction energetics and the ability to assess the potential effects of sub-cellular level phenomena such as calcium oscillations on lymphangion outflow. The maximum efficiency value predicted (40%) is at the upper end of estimates for other muscle types. Spontaneous calcium oscillations during diastole were found to increase outflow up to approximately 50% in the range of frequencies and amplitudes tested.

Pump efficacy in a two-dimensional, fluid-structure interaction model of a chain of contracting lymphangions

The transport of lymph through the lymphatic vasculature is the mechanism for returning excess interstitial fluid to the circulatory system, and it is essential for fluid homeostasis. Collecting lymphatic vessels comprise a significant portion of the lymphatic vasculature and are divided by valves into contractile segments known as lymphangions. Despite its importance, lymphatic transport in collecting vessels is not well understood. We present a computational model to study lymph flow through chains of valved, contracting lymphangions. We used the Navier-Stokes equations to model the fluid flow and the immersed boundary method to handle the two-way, fluid-structure interaction in 2D, non-axisymmetric simulations. We used our model to evaluate the effects of chain length, contraction style, and adverse axial pressure difference (AAPD) on cycle-mean flow rates (CMFRs). In the model, longer lymphangion chains generally yield larger CMFRs, and they fail to generate positive CMFRs at higher AAPDs than shorter chains. Simultaneously contracting pumps generate the largest CMFRs at nearly every AAPD and for every chain length. Due to the contraction timing and valve dynamics, non-simultaneous pumps generate lower CMFRs than the simultaneous pumps; the discrepancy diminishes as the AAPD increases. Valve dynamics vary with the contraction style and exhibit hysteretic opening and closing behaviors. Our model provides insight into how contraction propagation affects flow rates and transport through a lymphangion chain.

Fluid pumping of peristaltic vessel fitted with elastic valves

Using numerical simulations, we probe the fluid flow in an axisymmetric peristaltic vessel fitted with elastic bi-leaflet valves. In this biomimetic system that mimics the flow generated in lymphatic vessels, we investigate the effects of the valve and vessel properties on pumping performance of the valved peristaltic vessel. The results indicate that valves significantly increase pumping by reducing backflow. The presence of valves, however, increases the viscous resistance therefore requiring greater work compared to valveless vessels. The benefit of the valves is the most significant when the fluid is pumped against an adverse pressure gradient and for low vessel contraction wave speeds. We identify the optimum vessel and valve parameters leading to the maximum pumping efficiency. We show that the optimum valve elasticity maximizes the pumping flow rate by allowing the valve to block more effectively the backflow while maintaining low resistance during the forward flow. We also examine the pumping in vessels where the vessel contraction amplitude is a function of the adverse pressure gradient as found in lymphatic vessels. We find that in this case the flow is limited by the work generated by the contracting vessel, suggesting that the pumping in lymphatic vessels is constrained by the performance of lymphatic muscle. Given the regional heterogeneity of valve morphology observed throughout the lymphatic vasculature, these results provide insight into how these variations might facilitate efficient lymphatic transport in the vessel's local physiologic context.

A semi-automated finite difference mesh creation method for use with immersed boundary software IB2d and IBAMR

Numerous fluid-structure interaction problems in biology have been investigated using the immersed boundary method. The advantage of this method is that complex geometries, e.g., internal or external morphology, can easily be handled without the need to generate matching grids for both the fluid and the structure. Consequently, the difficulty of modeling the structure lies often in discretizing the boundary of the complex geometry (morphology). Both commercial and open source mesh generators for finite element methods have long been established; however, the traditional immersed boundary method is based on a finite difference discretization of the structure. Here we present a software library for obtaining finite difference discretizations of boundaries for direct use in the 2D immersed boundary method. This library provides tools for extracting such boundaries as discrete mesh points from digital images. We give several examples of how the method can be applied that include passing flow through the veins of insect wings, within lymphatic capillaries, and around starfish using open-source immersed boundary software.

Characterization of rat tail lymphatic contractility and biomechanics: incorporating nitric oxide-mediated vasoregulation

The lymphatic system transports lymph from the interstitial space back to the great veins via a series of orchestrated contractions of chains of lymphangions. Biomechanical models of lymph transport, validated with ex vivo or in vivo experimental results, have proved useful in revealing novel insight into lymphatic pumping; however, a need remains to characterize the contributions of vasoregulatory compounds in these modelling tools. Nitric oxide (NO) is a key mediator of lymphatic pumping. We quantified the active contractile and passive biaxial biomechanical response of rat tail collecting lymphatics and changes in the contractile response to the exogenous NO administration and integrated these findings into a biomechanical model. The passive mechanical response was characterized with a three-fibre family model. Nonlinear regression and non-parametric bootstrapping were used to identify best-fit material parameters to passive cylindrical biaxial mechanical data, assessing uniqueness and parameter confidence intervals; this model yielded a good fit (R² = 0.90). Exogenous delivery of NO via sodium nitroprusside (SNP) elicited a dose-dependent suppression of contractions; the amplitude of contractions decreased by 30% and the contraction frequency decreased by 70%. Contractile function was characterized with a modified Rachev-Hayashi model, introducing a parameter that is related to SNP concentration; the model provided a good fit (R² = 0.89) to changes in contractile responses to varying concentrations of SNP. These results demonstrated the significant role of NO in lymphatic pumping and provide a predictive biomechanical model to integrate the combined effect of mechanical loading and NO on lymphatic contractility and mechanical response.

The mechanical responses of advecting cells in confined flow

Fluid dynamics have long influenced cells in suspension. Red blood cells and white blood cells are advected through biological microchannels in both the cardiovascular and lymphatic systems and, as a result, are subject to a wide variety of complex fluidic forces as they pass through. In vivo, microfluidic forces influence different biological processes such as the spreading of infection, cancer metastasis, and cell viability, highlighting the importance of fluid dynamics in the blood and lymphatic vessels. This suggests that in vitro devices carrying cell suspensions may influence the viability and functionality of cells. Lab-on-a-chip, flow cytometry, and cell therapies involve cell suspensions flowing through microchannels of approximately 100-800  μ m. This review begins by examining the current fundamental theories and techniques behind the fluidic forces and inertial focusing acting on cells in suspension, before exploring studies that have investigated how these fluidic forces affect the reactions of suspended cells. In light of these studies' findings, both in vivo and in vitro fluidic cell microenvironments shall also be discussed before concluding with recommendations for the field.

Modelling secondary lymphatic valves with a flexible vessel wall: how geometry and material properties combine to provide function

A three-dimensional finite-element fluid/structure interaction model of an intravascular lymphatic valve was constructed, and its properties were investigated under both favourable and adverse pressure differences, simulating valve opening and valve closure, respectively. The shear modulus of the neo-Hookean material of both vascular wall and valve leaflet was varied, as was the degree of valve opening at rest. Also investigated was how the valve characteristics were affected by prior application of pressure inflating the whole valve. The characteristics were parameterised by the volume flow rate through the valve, the hydraulic resistance to flow, and the maximum sinus radius and inter-leaflet-tip gap on the plane of symmetry bisecting the leaflet, all as functions of the applied pressure difference. Maximum sinus radius on the leaflet-bisection plane increased with increasing pressure applied to either end of the valve segment, but also reflected the non-circular deformation of the sinus cross section caused by the leaflet, such that it passed through a minimum at small favourable pressure differences. When the wall was stiff, the inter-leaflet gap increased sigmoidally during valve opening; when it was as flexible as the leaflet, the gap increased more linearly. Less pressure difference was required both to open and to close the valve when either the wall or the leaflet material was more flexible. The degree of bias of the valve characteristics to the open position increased with the inter-leaflet gap in the resting position and with valve inflation pressure. The characteristics of the simulated valve were compared with those specified in an existing lumped-parameter model of one or more collecting lymphangions and used to estimate a revised value for the constant in that model which controls the rate of valve opening/closure with variation in applied pressure difference. The effects of the revised value on the lymph pumping efficacy predicted by the lumped-parameter model were evaluated.

Lymphatic remodelling in response to lymphatic injury in the hind limbs of sheep

Contractile activity in the lymphatic vasculature is essential for maintaining fluid balance within organs and tissues. However, the mechanisms by which collecting lymphatics adapt to changes in fluid load and how these adaptations influence lymphatic contractile activity are unknown. Here we report a model of lymphatic injury based on the ligation of one of two parallel lymphatic vessels in the hind limb of sheep and the evaluation of structural and functional changes in the intact, remodelling lymphatic vessel over a 42-day period. We show that the remodelled lymphatic vessel displayed increasing intrinsic contractile frequency, force generation and vessel compliance, as well as decreasing flow-mediated contractile inhibition via the enzyme endothelial nitric oxide synthase. A computational model of a chain of lymphatic contractile segments incorporating these adaptations predicted increases in the flow-generation capacity of the remodelled vessel at the expense of normal mitochondrial function and elevated oxidative stress within the lymphatic muscle. Our findings may inform interventions for mitigating lymphatic muscle fatigue in patients with dysfunctional lymphatics.

A computational model of a network of initial lymphatics and pre-collectors with permeable interstitium

Initial lymphatic vessels are made up of overlapped endothelial cells that act as unidirectional valves enabling one-way drainage of tissue fluid into the lumen of the initial lymphatics when there is a favourable pressure gradient. Initial lymphatics subsequently drain this fluid into the collecting lymphatics. This paper describes a computational model for a network of passive rat mesenteric lymphatic vessels with sparse secondary valves. The network was simulated with the secondary valves both operational and non-operational. The effects on the cycle-mean outflow-rate from the network of both inflammation and the resistance of the surrounding interstitium were considered. The cycle-mean outflow-rate is sensitive to vessel stiffness. If the influence of primary-valve resistance is reduced relative to that of interstitial resistance and intravascular resistance, there is no absolute advantage of extrinsic pumping, since maximum outflow-rate occurs when vessels are rigid. However, there is relative advantage, in that the outflow-rate at intermediate stiffness is higher with the secondary valves functioning than when they are deactivated. If primary-valve resistance dominates, then extrinsic pumping of non-rigid vessels provides absolute advantage. The nonlinear relation between pressure drop and flow-rate of the endothelial primary valves, combined with downstream compliance and pulsatile external pressure, constitutes a separate mechanism of pumping. By enabling the consideration of interactions between multiple phenomena (primary valves, secondary valves, a real network geometry with multiple branches, deformable vessel walls, interstitial resistance and external pressures), the model offers a perspective for delineating physiological phenomena that have not yet been fully linked to the biomechanics of fluid flow through initial lymphatic networks.

The effects of valve leaflet mechanics on lymphatic pumping assessed using numerical simulations

The lymphatic system contains intraluminal leaflet valves that function to bias lymph flow back towards the heart. These valves are present in the collecting lymphatic vessels, which generally have lymphatic muscle cells and can spontaneously pump fluid. Recent studies have shown that the valves are open at rest, can allow some backflow, and are a source of nitric oxide (NO). To investigate how these valves function as a mechanical valve and source of vasoactive species to optimize throughput, we developed a mathematical model that explicitly includes Ca²⁺ -modulated contractions, NO production and valve structures. The 2D lattice Boltzmann model includes an initial lymphatic vessel and a collecting lymphangion embedded in a porous tissue. The lymphangion segment has mechanically-active vessel walls and is flanked by deformable valves. Vessel wall motion is passively affected by fluid pressure, while active contractions are driven by intracellular Ca²⁺ fluxes. The model reproduces NO and Ca²⁺ dynamics, valve motion and fluid drainage from tissue. We find that valve structural properties have dramatic effects on performance, and that valves with a stiffer base and flexible tips produce more stable cycling. In agreement with experimental observations, the valves are a major source of NO. Once initiated, the contractions are spontaneous and self-sustained, and the system exhibits interesting non-linear dynamics. For example, increased fluid pressure in the tissue or decreased lymph pressure at the outlet of the system produces high shear stress and high levels of NO, which inhibits contractions. On the other hand, a high outlet pressure opposes the flow, increasing the luminal pressure and the radius of the vessel, which results in strong contractions in response to mechanical stretch of the wall. We also find that the location of contraction initiation is affected by the extent of backflow through the valves.

Inhibition of contraction strength and frequency by wall shear stress in a single-lymphangion model

The phasic contractions of collecting lymphatic vessels are reduced in strength and occur at diminished frequency when the favourable pressure difference and the resulting antegrade flow create large fluid shear stresses at the luminal surface. This paper describes a minimal phenomenological model of this mechanism, that is applied to a previously validated numerical model of a phasically contracting lymphangion. The parameters of the inhibition model are quantitatively matched to observations in isolated segments of rat lymphatic vessel, first for mesenteric lymphatics then for thoracic duct, and outcomes from the numerical model are then qualitatively compared with recent observations in isolated segments of rat thoracic duct.

Lymphatic Vessel Network Structure and Physiology

The lymphatic system is comprised of a network of vessels interrelated with lymphoid tissue, which has the holistic function to maintain the local physiologic environment for every cell in all tissues of the body. The lymphatic system maintains extracellular fluid homeostasis favorable for optimal tissue function, removing substances that arise due to metabolism or cell death, and optimizing immunity against bacteria, viruses, parasites, and other antigens. This article provides a comprehensive review of important findings over the past century along with recent advances in the understanding of the anatomy and physiology of lymphatic vessels, including tissue/organ specificity, development, mechanisms of lymph formation and transport, lymphangiogenesis, and the roles of lymphatics in disease. © 2019 American Physiological Society. Compr Physiol 9:207-299, 2019.

A one-dimensional mathematical model of collecting lymphatics coupled with an electro-fluid-mechanical contraction model and valve dynamics

We propose a one-dimensional model for collecting lymphatics coupled with a novel Electro-Fluid-Mechanical Contraction (EFMC) model for dynamical contractions, based on a modified FitzHugh-Nagumo model for action potentials. The one-dimensional model for a deformable lymphatic vessel is a nonlinear system of hyperbolic Partial Differential Equations (PDEs). The EFMC model combines the electrical activity of lymphangions (action potentials) with fluid-mechanical feedback (circumferential stretch of the lymphatic wall and wall shear stress) and lymphatic vessel wall contractions. The EFMC model is governed by four Ordinary Differential Equations (ODEs) and phenomenologically relies on: (1) environmental calcium influx, (2) stretch-activated calcium influx, and (3) contraction inhibitions induced by wall shear stresses. We carried out a stability analysis of the stationary state of the EFMC model. Contractions turn out to be triggered by the instability of the stationary state. Overall, the EFMC model allows emulating the influence of pressure and wall shear stress on the frequency of contractions observed experimentally. Lymphatic valves are modelled by extending an existing lumped-parameter model for blood vessels. Modern numerical methods are employed for the one-dimensional model (PDEs), for the EFMC model and valve dynamics (ODEs). Adopting the geometrical structure of collecting lymphatics from rat mesentery, we apply the full mathematical model to a carefully selected suite of test problems inspired by experiments. We analysed several indices of a single lymphangion for a wide range of upstream and downstream pressure combinations which included both favourable and adverse pressure gradients. The most influential model parameters were identified by performing two sensitivity analyses for favourable and adverse pressure gradients.

Contraction of collecting lymphatics: organization of pressure-dependent rate for multiple lymphangions

The paper describes the extension of a previously developed model of pressure-dependent contraction rate to the case of multiple lymphangions. Mechanical factors are key modulators of active lymphatic pumping. As part of the evolution of our lumped-parameter model to match experimental findings, we have designed an algorithm whereby the time until the next contraction depends on lymphangion transmural pressure in the contraction just completed. The functional dependence of frequency on pressure is quantitatively matched to isobaric contraction experiments on isolated lymphatic segments. When each of several lymphangions is given this ability, a scheme for their coordination must be instituted to match the observed synchronization. Accordingly, and in line with an experiment on an isolated lymphatic vessel segment in which we measured contraction sequence and conduction delay, we took the fundamental principle to be that local timing can be overridden by signals to initiate contraction that start in adjacent lymphangions, conducted with a short delay. The scheme leads to retrograde conduction when the lymphangion chain is pumping against an adverse pressure difference, but antegrade conduction when contractions occur with no or a favourable pressure difference. Abolition of these conducted signals leads to chaotic variation of cycle-mean flow-rate from the chain, diastolic duration in each lymphangion, and inter-lymphangion delays. Chaotic rhythm is also seen under other circumstances. Because the model responds to increasing adverse pressure difference by increasing the repetition rate of contractions, it maintains time-average output flow-rate better than one with fixed repetition rate.

Lymphatic System Flows

The supply of oxygen and nutrients to tissues is performed by the blood system, and involves a net leakage of fluid outward at the capillary level. One of the principal functions of the lymphatic system is to gather this fluid and return it to the blood system to maintain overall fluid balance. Fluid in the interstitial spaces is often at subatmospheric pressure, and the return points into the venous system are at pressures of approximately 20 cmH₂O. This adverse pressure difference is overcome by the active pumping of collecting lymphatic vessels, which feature closely spaced one-way valves and contractile muscle cells in their walls. Passive vessel squeezing causes further pumping. The dynamics of lymphatic pumping have been investigated experimentally and mathematically, revealing complex behaviours indicating that the system performance is robust against minor perturbations in pressure and flow. More serious disruptions can lead to incurable swelling of tissues called lymphœdema.

Lymphatic System Flows

The supply of oxygen and nutrients to tissues is performed by the blood system, and involves a net leakage of fluid outward at the capillary level. One of the principal functions of the lymphatic system is to gather this fluid and return it to the blood system to maintain overall fluid balance. Fluid in the interstitial spaces is often at subatmospheric pressure, and the return points into the venous system are at pressures of approximately 20 cmH₂O. This adverse pressure difference is overcome by the active pumping of collecting lymphatic vessels, which feature closely spaced one-way valves and contractile muscle cells in their walls. Passive vessel squeezing causes further pumping. The dynamics of lymphatic pumping have been investigated experimentally and mathematically, revealing complex behaviours indicating that the system performance is robust against minor perturbations in pressure and flow. More serious disruptions can lead to incurable swelling of tissues called lymphœdema.

Characterization of Mouse Mesenteric Lymphatic Valve Structure and Function

Intraluminal valves of collecting lymphatic vessels ensure unidirectional lymph transport against hydrostatic pressure gradient. Mouse mesentery harbors up to 800 valves and represents a convenient model for lymphatic valve quantification, high resolution imaging of different stages of valve development as well as for analysis of valve function. The protocol describes embryonic and postnatal mesenteric lymphatic vessel preparation for whole-mount immunofluorescent staining and visualization of valve organization, quantification of main morphological parameters such as valve size and leaflet length, and the quantitative assessment of functional properties of adult valves using back-leak and closure tests.

Integrated geometric and mechanical analysis of an image-based lymphatic valve

Lymphatic valves facilitate the lymphatic system's role in maintaining fluid homeostasis. Malformed valves are found in several forms of primary lymphœdema, resulting in incurable swelling of the tissues and immune dysfunction. Their experimental study is complicated by their small size and operation in low pressure and low Reynolds number environments. Mathematical models of these structures can give insight and complement experimentation. In this work, we present the first valve geometry reconstructed from confocal imagery and used in the construction of a subject-specific model in a closing mode. A framework is proposed whereby an image is converted into a valve model. An FEA study was performed to identify the significance of the shear modulus, the consequences of smoothing the leaflet surface and the effect of wall motion on valve behaviour. Smoothing is inherent to any analysis from imagery. The nature of the image, segmentation and meshing all cause attenuation of high-frequency features. Smoothing not only causes loss of surface area but also the loss of high-frequency geometric features which may reduce stiffness. This work aimed to consider these effects and inform studies by taking a manual reconstruction and through manifold harmonic analysis, attenuating higher frequency features to replicate lower resolution images or lower degree-of-freedom reconstructions. In conclusion, two metrics were considered: trans-valvular pressure required to close the valve, ΔP_c, and the retrograde volume displacement after closure. The higher ΔP_c, the greater the volume of lymph that will pass through the valve during closure. Retrograde volume displacement after closure gives a metric of compliance of the valve and for the quality of the valve seal. In the case of the image-specific reconstructed valve, removing features with a wavelength longer than four μm caused changes in ΔP_c. Varying the shear modulus from 10 kPa to 60 kPa caused a 3.85-fold increase in the retrograde volume displaced. The inclusion of a non-rigid wall caused ΔP_c to increase from 1.56 to 2.52 cmH₂O.

Demonstration and Analysis of the Suction Effect for Pumping Lymph from Tissue Beds at Subatmospheric Pressure

Many tissues exhibit subatmospheric interstitial pressures under normal physiologic conditions. The mechanisms by which the lymphatic system extracts fluid from these tissues against the overall pressure gradient are unknown. We address this important physiologic issue by combining experimental measurements of contractile function and pressure generation with a previously validated mathematical model. We provide definitive evidence for the existence of 'suction pressure' in collecting lymphatic vessels, which manifests as a transient drop in pressure downstream of the inlet valve following contraction. This suction opens the inlet valve and is required for filling in the presence of low upstream pressure. Positive transmural pressure is required for this suction, providing the energy required to reopen the vessel. Alternatively, external vessel tethering can serve the same purpose when the transmural pressure is negative. Suction is transmitted upstream, allowing fluid to be drawn in through initial lymphatics. Because suction plays a major role in fluid entry to the lymphatics and is affected by interstitial pressure, our results introduce the phenomenon as another important factor to consider in the study of lymphoedema and its treatment.

Opportunities for Studying the Hydrodynamic Context for Breast Cancer Cell Spread Through Lymph Flow

The lymphatic system serves as the primary route for the metastatic spread of breast cancer cells (BCCs). A scarcity of information exists with regard to the advection of BCCs in lymph flow and a fundamental understanding of the response of BCCs to the forces in the lymphatics needs to be established. This review summarizes the flow environment metastatic BCCs are exposed to in the lymphatics. Special attention is paid to the behavior of cells/particles in microflows in an attempt to elucidate the behavior of BCCs under lymph flow conditions (Reynolds number <1).

Valve-related modes of pump failure in collecting lymphatics: numerical and experimental investigation

Lymph is transported along collecting lymphatic vessels by intrinsic and extrinsic pumping. The walls have muscle of a type intermediate between blood-vascular smooth muscle and myocardium; a contracting segment between two valves (a lymphangion) constitutes a pump. This intrinsic mechanism is investigated ex vivo in isolated, spontaneously contracting, perfused segments subjected to controlled external pressures. The reaction to varying afterload is probed by slowly ramping up the outlet pressure until pumping fails. Often the failure occurs when the contraction raises intra-lymphangion pressure insufficiently to overcome the outlet pressure, open the outlet valve and cause ejection, but many segments fail by other means, the mechanisms of which are not clear. We here elucidate those mechanisms by resort to a numerical model. Experimental observations are paired with comparable findings from computer simulations, using a lumped-parameter model that incorporates previously measured valve properties, plus new measurements of active contractile and passive elastic properties, and the dependence of contraction frequency on transmural pressure, all taken from isobaric twitch contraction experiments in the same vessel. Surprisingly, the model predicts seven different possible modes of pump failure, each defined by a different sequence of valve events, with their occurrence depending on the parameter values and boundary conditions. Some, but not all, modes were found experimentally. Further model investigation reveals routes by which a vessel exhibiting one mode of failure might under altered circumstances exhibit another.

A lumped parameter model of mechanically mediated acute and long-term adaptations of contractility and geometry in lymphatics for characterization of lymphedema

A primary purpose of the lymphatic system is to transport fluid from peripheral tissues to the central venous system in order to maintain tissue-fluid balance. Failure to perform this task results in lymphedema marked by swelling of the affected limb as well as geometric remodeling and reduced contractility of the affected lymphatic vessels. The mechanical environment has been implicated in the regulation of lymphatic contractility, but it is unknown how changes in the mechanical environment are related to loss of contractile function and remodeling of the tissue. The purpose of this paper was to introduce a new theoretical framework for acute and long-term adaptations of lymphatic vessels to changes in mechanical loading. This theoretical framework combines a simplified version of a published lumped parameter model for lymphangion function and lymph transport, a published microstructurally motivated constitutive model for the active and passive mechanical behavior of isolated rat thoracic ducts, and novel models for acute mechanically mediated vasoreactive adaptations and long-term volumetric growth to simulate changes in muscle contractility and geometry of a single isolated rat thoracic duct in response to a sustained elevation in afterload. The illustrative examples highlight the potential role of the mechanical environment in the acute maintenance of contractility and long-term geometric remodeling, presumably aimed at meeting fluid flow demands while also maintaining mechanical homeostasis. Results demonstrate that contractility may adapt in response to shear stress to meet fluid flow demands and show that pressure-induced long-term geometric remodeling may attenuate these adaptations and reduce fluid flow. The modeling framework and illustrative simulations help suggest relevant experiments that are necessary to accurately quantify and predict the acute and long-term adaptations of lymphangions to altered mechanical loading.

In vivo label-free measurement of lymph flow velocity and volumetric flow rates using Doppler optical coherence tomography

Direct in vivo imaging of lymph flow is key to understanding lymphatic system function in normal and disease states. Optical microscopy techniques provide the resolution required for these measurements, but existing optical techniques for measuring lymph flow require complex protocols and provide limited temporal resolution. Here, we describe a Doppler optical coherence tomography platform that allows direct, label-free quantification of lymph velocity and volumetric flow rates. We overcome the challenge of very low scattering by employing a Doppler algorithm that operates on low signal-to-noise measurements. We show that this technique can measure lymph velocity at sufficiently high temporal resolution to resolve the dynamic pulsatile flow in collecting lymphatic vessels.

Estimation of the Pressure Drop Required for Lymph Flow through Initial Lymphatic Networks

BACKGROUND

Lymphatic function is critical for maintaining interstitial fluid balance and is linked to multiple pathological conditions. While smooth muscle contractile mechanisms responsible for fluid flow through collecting lymphatic vessels are well studied, how fluid flows into and through initial lymphatic networks remains poorly understood. The objective of this study was to estimate the pressure difference needed for flow through an intact initial lymphatic network.

METHODS AND RESULTS

Pressure drops were computed for real and theoretical networks with varying branch orders using a segmental Poiseuille flow model. Vessel geometries per branch order were based on measurements from adult Wistar rat mesenteric initial lymphatic networks. For computational predications based on real network geometries and combinations of low or high output velocities (2 mm/s, 4 mm/s) and viscosities (1 cp, 1.5 cp), pressure drops were estimated to range 0.31-2.57 mmHg. The anatomical data for the real networks were also used to create a set of theoretical networks in order to identify possible minimum and maximum pressure drops. The pressure difference range for the theoretical networks was 0.16-3.16 mmHg.

CONCLUSIONS

The results support the possibility for suction pressures generated from cyclic smooth muscle contractions of upstream collecting lymphatics being sufficient for fluid flow through an initial lymphatic network.

Pump function curve shape for a model lymphatic vessel

The transport capacity of a contractile segment of lymphatic vessel is defined by its pump function curve relating mean flow-rate and adverse pressure difference. Numerous system characteristics affect curve shape and the magnitude of the generated flow-rates and pressures. Some cannot be varied experimentally, but their separate and interacting effects can be systematically revealed numerically. This paper explores variations in the rate of change of active tension and the form of the relation between active tension and muscle length, factors not known from experiment to functional precision. Whether the pump function curve bends toward or away from the origin depends partly on the curvature of the passive pressure-diameter relation near zero transmural pressure, but rather more on the form of the relation between active tension and muscle length. A pump function curve bending away from the origin defines a well-performing pump by maximum steady output power. This behaviour is favoured by a length/active-tension relationship which sustains tension at smaller lengths. Such a relationship also favours high peak mechanical efficiency, defined as output power divided by the input power obtained from the lymphangion diameter changes and active-tension time-course. The results highlight the need to pin down experimentally the form of the length/active-tension relationship.

Quantification of the passive and active biaxial mechanical behaviour and microstructural organization of rat thoracic ducts

Mechanical loading conditions are likely to play a key role in passive and active (contractile) behaviour of lymphatic vessels. The development of a microstructurally motivated model of lymphatic tissue is necessary for quantification of mechanically mediated maladaptive remodelling in the lymphatic vasculature. Towards this end, we performed cylindrical biaxial testing of Sprague-Dawley rat thoracic ducts (n = 6) and constitutive modelling to characterize their mechanical behaviour. Spontaneous contraction was quantified at transmural pressures of 3, 6 and 9 cmH2O. Cyclic inflation in calcium-free saline was performed at fixed axial stretches between 1.30 and 1.60, while recording pressure, outer diameter and axial force. A microstructurally motivated four-fibre family constitutive model originally proposed by Holzapfel et al. (Holzapfel et al. 2000 J. Elast. 61, 1-48. (doi:10.1023/A:1010835316564)) was used to quantify the passive mechanical response, and the model of Rachev and Hayashi was used to quantify the active (contractile) mechanical response. The average error between data and theory was 8.9 ± 0.8% for passive data and 6.6 ± 2.6% and 6.8 ± 3.4% for the systolic and basal conditions, respectively, for active data. Multi-photon microscopy was performed to quantify vessel wall thickness (32.2 ± 1.60 µm) and elastin and collagen organization for three loading conditions. Elastin exhibited structural 'fibre families' oriented nearly circumferentially and axially. Sample-to-sample variation was observed in collagen fibre distributions, which were often non-axisymmetric, suggesting material asymmetry. In closure, this paper presents a microstructurally motivated model that accurately captures the biaxial active and passive mechanical behaviour in lymphatics and offers potential for future research to identify parameters contributing to mechanically mediated disease development.

Network Scale Modeling of Lymph Transport and Its Effective Pumping Parameters

The lymphatic system is an open-ended network of vessels that run in parallel to the blood circulation system. These vessels are present in almost all of the tissues of the body to remove excess fluid. Similar to blood vessels, lymphatic vessels are found in branched arrangements. Due to the complexity of experiments on lymphatic networks and the difficulty to control the important functional parameters in these setups, computational modeling becomes an effective and essential means of understanding lymphatic network pumping dynamics. Here we aimed to determine the effect of pumping coordination in branched network structures on the regulation of lymph flow. Lymphatic vessel networks were created by building upon our previous lumped-parameter model of lymphangions in series. In our network model, each vessel is itself divided into multiple lymphangions by lymphatic valves that help maintain forward flow. Vessel junctions are modeled by equating the pressures and balancing mass flows. Our results demonstrated that a 1.5 s rest-period between contractions optimizes the flow rate. A time delay between contractions of lymphangions at the junction of branches provided an advantage over synchronous pumping, but additional time delays within individual vessels only increased the flow rate for adverse pressure differences greater than 10.5 cmH2O. Additionally, we quantified the pumping capability of the system under increasing levels of steady transmural pressure and outflow pressure for different network sizes. We observed that peak flow rates normally occurred under transmural pressures between 2 to 4 cmH2O (for multiple pressure differences and network sizes). Networks with 10 lymphangions per vessel had the highest pumping capability under a wide range of adverse pressure differences. For favorable pressure differences, pumping was more efficient with fewer lymphangions. These findings are valuable for translating experimental measurements from the single lymphangion level to tissue and organ scales.

Consequences of intravascular lymphatic valve properties: a study of contraction timing in a multi-lymphangion model

The observed properties of valves in collecting lymphatic vessels include transmural pressure-dependent bias to the open state and hysteresis. The bias may reduce resistance to flow when the vessel is functioning as a conduit. However, lymphatic pumping implies a streamwise increase in mean pressure across each valve, suggesting that the bias is then potentially unhelpful. Lymph pumping by a model of several collecting lymphatic vessel segments (lymphangions) in series, which incorporated these properties, was investigated under conditions of adverse pressure difference while varying the refractory period between active muscular contractions and the inter-lymphangion contraction delay. It was found that many combinations of the timing parameters and the adverse pressure difference led to one or more intermediate valves remaining open instead of switching between open and closed states during repetitive contraction cycles. Cyclic valve switching was reliably indicated if the mean pressure in a lymphangion over a cycle was higher than that in the lymphangion upstream, but either lack of or very brief valve closure could cause mean pressure to be lower downstream. Widely separated combinations of refractory period and delay time were found to produce the greatest flow-rate for a given pressure difference. The efficiency of pumping was always maximized by a long refractory period and lymphangion contraction starting when the contraction of the lymphangion immediately upstream was peaking. By means of an ex vivo experiment, it was verified that intermediate valves in a chain of pumping lymphangions can remain open, while the lymphangions on either side of the open valve continue to execute contractions.

Determining the combined effect of the lymphatic valve leaflets and sinus on resistance to forward flow

The lymphatic system is vital to a proper maintenance of fluid and solute homeostasis. Collecting lymphatics are composed of actively contracting tubular vessels segmented by bulbous sinus regions that encapsulate bi-leaflet check valves. Valve resistance to forward flow strongly influences pumping performance. However, because of the sub-millimeter size of the vessels with flow rates typically <1 ml/h and pressures of a few cmH2O, resistance is difficult to measure experimentally. Using a newly defined idealized geometry, we employed an uncoupled approach where the solid leaflet deflections of the open valve were computed and lymph flow calculations were subsequently performed. We sought to understand: 1) the effect of sinus and leaflet size on the resulting deflections experienced by the valve leaflets and 2) the effects on valve resistance to forward flow of the fully open valve. For geometries with sinus-to-root diameter ratios >1.39, the average resistance to forward flow was 0.95×10(6)[g/(cm4 s)]. Compared to the viscous pressure drop that would occur in a straight tube the same diameter as the upstream lymphangion, valve leaflets alone increase the pressure drop up to 35%. However, the presence of the sinus reduces viscous losses, with the net effect that when combined with leaflets the overall resistance is less than that of the equivalent continuing straight tube. Accurately quantifying resistance to forward flow will add to the knowledge used to develop therapeutics for treating lymphatic disorders and may eventually lead to understanding some forms of primary lymphedema.

Mechanobiology of lymphatic contractions

The lymphatic system is responsible for controlling tissue fluid pressure by facilitating flow of lymph (i.e. the plasma and cells that enter the lymphatic system). Because lymph contains cells of the immune system, its transport is not only important for fluid homeostasis, but also immune function. Lymph drainage can occur via passive flow or active pumping, and much research has identified the key biochemical and mechanical factors that affect output. Although many studies and reviews have addressed how tissue properties and fluid mechanics (i.e. pressure gradients) affect lymph transport [1-3] there is less known about lymphatic mechanobiology. As opposed to passive mechanical properties, mechanobiology describes the active coupling of mechanical signals and biochemical pathways. Lymphatic vasomotion is the result of a fascinating system affected by mechanical forces exerted by the flowing lymph, including pressure-induced vessel stretch and flow-induced shear stresses. These forces can trigger or modulate biochemical pathways important for controlling the lymphatic contractions. Here, I review the current understanding of lymphatic vessel function, focusing on vessel mechanobiology, and summarize the prospects for a comprehensive understanding that integrates the mechanical and biomechanical control mechanisms in the lymphatic system.

Mechanical forces and lymphatic transport

This review examines the current understanding of how the lymphatic vessel network can optimize lymph flow in response to various mechanical forces. Lymphatics are organized as a vascular tree, with blind-ended initial lymphatics, precollectors, prenodal collecting lymphatics, lymph nodes, postnodal collecting lymphatics and the larger trunks (thoracic duct and right lymph duct) that connect to the subclavian veins. The formation of lymph from interstitial fluid depends heavily on oscillating pressure gradients to drive fluid into initial lymphatics. Collecting lymphatics are segmented vessels with unidirectional valves, with each segment, called a lymphangion, possessing an intrinsic pumping mechanism. The lymphangions propel lymph forward against a hydrostatic pressure gradient. Fluid is returned to the central circulation both at lymph nodes and via the larger lymphatic trunks. Several recent developments are discussed, including evidence for the active role of endothelial cells in lymph formation; recent developments on how inflow pressure, outflow pressure, and shear stress affect the pump function of the lymphangion; lymphatic valve gating mechanisms; collecting lymphatic permeability; and current interpretations of the molecular mechanisms within lymphatic endothelial cells and smooth muscle. An improved understanding of the physiological mechanisms by which lymphatic vessels sense mechanical stimuli, integrate the information, and generate the appropriate response is key for determining the pathogenesis of lymphatic insufficiency and developing treatments for lymphedema.

Primary and secondary lymphatic valve development: molecular, functional and mechanical insights

Fluid homeostasis in vertebrates critically relies on the lymphatic system forming a hierarchical network of lymphatic capillaries and collecting lymphatics, for the efficient drainage and transport of extravasated fluid back to the cardiovascular system. Blind-ended lymphatic capillaries employ specialized junctions and anchoring filaments to encourage a unidirectional flow of the interstitial fluid into the initial lymphatic vessels, whereas collecting lymphatics are responsible for the active propulsion of the lymph to the venous circulation via the combined action of lymphatic muscle cells and intraluminal valves. Here we describe recent findings on molecular and physical factors regulating the development and maturation of these two types of valves and examine their role in tissue-fluid homeostasis.

Parameter sensitivity analysis of a lumped-parameter model of a chain of lymphangions in series

Any disruption of the lymphatic system due to trauma or injury can lead to edema. There is no effective cure for lymphedema, partly because predictive knowledge of lymphatic system reactions to interventions is lacking. A well-developed model of the system could greatly improve our understanding of its function. Lymphangions, defined as the vessel segment between two valves, are the individual pumping units. Based on our previous lumped-parameter model of a chain of lymphangions, this study aimed to identify the parameters that affect the system output the most using a sensitivity analysis. The system was highly sensitive to minimum valve resistance, such that variations in this parameter caused an order-of-magnitude change in time-average flow rate for certain values of imposed pressure difference. Average flow rate doubled when contraction frequency was increased within its physiological range. Optimum lymphangion length was found to be some 13-14.5 diameters. A peak of time-average flow rate occurred when transmural pressure was such that the pressure-diameter loop for active contractions was centered near maximum passive vessel compliance. Increasing the number of lymphangions in the chain improved the pumping in the presence of larger adverse pressure differences. For a given pressure difference, the optimal number of lymphangions increased with the total vessel length. These results indicate that further experiments to estimate valve resistance more accurately are necessary. The existence of an optimal value of transmural pressure may provide additional guidelines for increasing pumping in areas affected by edema.

Incorporating measured valve properties into a numerical model of a lymphatic vessel

An existing lumped-parameter model of multiple lymphangions (lymphatic vascular segments) in series is adapted for the incorporation of recent physiological measurements of lymphatic vascular properties. The new data show very marked nonlinearity of the passive pressure-diameter relation during distension, relative to comparable blood vessels, and complex valve behaviour. Since lymph is transported as a result of either the active contraction or the passive squeezing of vascular segments situated between two one-way valves, the performance of these valves is of primary importance. The valves display hysteresis (the opening and closing pressure drop thresholds differ), a bias to staying open (both state changes occur when the trans-valve pressure drop is adverse) and pressure-drop threshold dependence on transmural pressure. These properties, in combination with the strong nonlinearity that valve operation represents, have in turn caused intriguing numerical problems in the model, and we describe numerical stratagems by which we have overcome the problems. The principal problem is also generalised into a relatively simple mathematical example, for which solution detail is provided using two different solvers.