Probing the effect of morphology on lymphatic valve dynamic function

Numerical investigation on flow-induced wall shear stress variation of metastatic cancer cells in lymphatics with elastic valves

Hematogenous metastasis occurs when cancer cells detach from the extracellular matrix in the primary tumor into the bloodstream or lymphatic system. Elucidating the response of metastatic tumor cells in suspension to the flow conditions in lymphatics with valves from a mechanical/fluidic perspective is necessary. A physiologically relevant computational model of a lymphatic vessel with valves was constructed using fully coupled fluid-cell-vessel interactions to investigate the effects of lymphatic vessel contractility, valve properties, and cell size and stiffness on the variations in magnitude and gradient of the flow-induced wall shear stress (WSS) experienced by suspended tumor cells. Results indicated that the maximum WSSmax increased with the increments in cell diameter, vessel contraction amplitude, and valve stiffness. The decrease in vessel contraction period and valve aspect ratio also increased the maximum WSSmax. The influence of the properties of the valve on the WSS was more significant among the factors mentioned above. The maximum WSSmax acting on the cancer cell when the cell reversed the direction of its motion in the valve region increased by 0.5-1.4 times that before the cell entered the valve region. The maximum change in WSS was in the range of 0.004-0.028 Pa/µm depending on the factors studied. They slightly exceeded the values associated with breast cancer cell apoptosis. The results of this study provide biofluid mechanics-based support for mechanobiological research on the metastasis of metastatic cancer cells in suspension within the lymphatics.

Genetic causes of lymphatic disorders: recent updates on the clinical and molecular aspects of lymphatic disease

PURPOSE OF REVIEW

The lymphatic system facilitates several key functions that limit significant morbidity and mortality. Despite the impact and burden of lymphatic disorders, there are many remaining disorders whose genetic substrate remains unknown. The purpose of this review is to provide an update on the genetic causes of lymphatic disorders, while reporting on newly proposed clinical classifications of lymphatic disease.

RECENT FINDINGS

We reviewed several new mutations in genes that have been identified as potential causes of lymphatic disorders including: MDFIC, EPHB 4 , and ANGPT2. Furthermore, the traditional St. George's Classification system for primary lymphatic anomalies has been updated to reflect the use of genetic testing, both as a tool for the clinical identification of lymphatic disease and as a method through which new sub-classifications of lymphatic disorders have been established within this framework. Finally, we highlighted recent clinical studies that have explored the impact of therapies such as sirolimus, ketoprofen, and acebilustat on lymphatic disorders.

SUMMARY

Despite a growing body of evidence, current literature demonstrates a persistent gap in the number of known genes responsible for lymphatic disease entities. Recent clinical classification tools have been introduced in order to integrate traditional symptom- and time-based diagnostic approaches with modern genetic classifications, as highlighted in the updated St. George's classification system. With the introduction of this novel approach, clinicians may be better equipped to recognize established disease and, potentially, to identify novel causal mutations. Further research is needed to identify additional genetic causes of disease and to optimize current clinical tools for diagnosis and treatment.

The Role of Inflammation in Lymphedema: A Narrative Review of Pathogenesis and Opportunities for Therapeutic Intervention

Lymphedema is a chronic and progressive disease of the lymphatic system characterized by inflammation, increased adipose deposition, and tissue fibrosis. Despite early hypotheses identifying lymphedema as a disease of mechanical lymphatic disruption alone, the progressive inflammatory nature underlying this condition is now well-established. In this review, we provide an overview of the various inflammatory mechanisms that characterize lymphedema development and progression. These mechanisms contribute to the acute and chronic phases of lymphedema, which manifest clinically as inflammation, fibrosis, and adiposity. Furthermore, we highlight the interplay between current therapeutic modalities and the underlying inflammatory microenvironment, as well as opportunities for future therapeutic development.

Computational fluid dynamic modeling of the lymphatic system: a review of existing models and future directions

Historically, research into the lymphatic system has been overlooked due to both a lack of knowledge and limited recognition of its importance. In the last decade however, lymphatic research has gained substantial momentum and has included the development of a variety of computational models to aid understanding of this complex system. This article reviews existing computational fluid dynamic models of the lymphatics covering each structural component including the initial lymphatics, pre-collecting and collecting vessels, and lymph nodes. This is followed by a summary of limitations and gaps in existing computational models and reasons that development in this field has been hindered to date. Over the next decade, efforts to further characterize lymphatic anatomy and physiology are anticipated to provide key data to further inform and validate lymphatic fluid dynamic models. Development of more comprehensive multiscale- and multi-physics computational models has the potential to significantly enhance the understanding of lymphatic function in both health and disease.

Effect of valve spacing on peristaltic pumping

Peristaltic fluid pumping due to a periodically propagating contraction wave in a vessel fitted with one-way elastic valves is investigated numerically. It is concluded that the valve spacing within the vessel relative to the contraction wavelength plays a critical role in providing efficient pumping. When the valve spacing does not match the wavelength, the valves open asynchronously and the volume of the vessel segments bounded by two consecutive valves changes periodically, thereby inducing volumetric fluid pumping. The volumetric pumping leads to higher pumping flowrate and efficiency against an adverse pressure gradient. The optimum pumping occurs when the ratio of valve spacing to contraction wavelength is about2/3. This pumping regime is characterized by a longer period during which the valves are open. The results are useful for further understanding the pumping features of lymphatic system and provide insight into the design of biomimetic pumping devices.

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.

Computational simulations of the effects of gravity on lymphatic transport

Physical forces, including mechanical stretch, fluid pressure, and shear forces alter lymphatic vessel contractions and lymph flow. Gravitational forces can affect these forces, resulting in altered lymphatic transport, but the mechanisms involved have not been studied in detail. Here, we combine a lattice Boltzmann-based fluid dynamics computational model with known lymphatic mechanobiological mechanisms to investigate the movement of fluid through a lymphatic vessel under the effects of gravity that may either oppose or assist flow. Regularly spaced, mechanical bi-leaflet valves in the vessel enforce net positive flow as the vessel walls contract autonomously in response to calcium and nitric oxide (NO) levels regulated by vessel stretch and shear stress levels. We find that large gravitational forces opposing flow can stall the contractions, leading to no net flow, but transient mechanical perturbations can re-establish pumping. In the case of gravity strongly assisting flow, the contractions also cease due to high shear stress and NO production, which dilates the vessel to allow gravity-driven flow. In the intermediate range of oppositional gravity forces, the vessel actively contracts to offset nominal gravity levels or to modestly assist the favorable hydrostatic pressure gradients.

Imaging of fluorescent polymer dots in relation to channels and immune cells in the lymphatic system

Polymer dots (Pdots) have been applied to imaging lymph nodes (LNs) and lymphatic vessels (LVs) in living mice and rats. However, the mechanism of absorption, distribution, metabolism, and excretion of Pdots in LNs and LVs is still unclear. Therefore, the relationship between Pdots and immune cells, LVs and collagen fibers in lymphatics was studied by multiple in vivo and ex vivo microscopic imaging methods and detection techniques. Flow cytometry showed that Pdots could be phagocytosed by macrophages and monocytes, and had no relationship with B cells, T cells and dendric cells in LNs. Silver staining, immunofluorescence and two-photon microscope showed that Pdots gathered in collagen fibers and LVs of LNs. Furthermore, immunofluorescence imaging results verified that Pdots were distributed in the extracellular space of collecting LVs endothelial cells. In addition, Pdots in the collecting LVs were basically cleared by leaking into the surrounding tissue or draining LNs after 21 days of injection. During the long-time observation, Pdots also helped monitor the contraction frequency and variation range of LV. Our study lays a foundation on the research of Pdots as the carrier to study lymphatic structure and function in the future.

Mathematical Modeling of Lymph Node Drainage Function by Neural Network

The lymph node (LN) represents a key structural component of the lymphatic system network responsible for the fluid balance in tissues and the immune system functioning. Playing an important role in providing the immune defense of the host organism, LNs can also contribute to the progression of pathological processes, e.g., the spreading of cancer cells. To gain a deeper understanding of the transport function of LNs, experimental approaches are used. Mathematical modeling of the fluid transport through the LN represents a complementary tool for studying the LN functioning under broadly varying physiological conditions. We developed an artificial neural network (NN) model to describe the lymph node drainage function. The NN model predicts the flow characteristics through the LN, including the exchange with the blood vascular systems in relation to the boundary and lymphodynamic conditions, such as the afferent lymph flow, Darcy’s law constants and Starling’s equation parameters. The model is formulated as a feedforward NN with one hidden layer. The NN complements the computational physics-based model of a stationary fluid flow through the LN and the fluid transport across the blood vessel system of the LN. The physical model is specified as a system of boundary integral equations (IEs) equivalent to the original partial differential equations (PDEs; Darcy’s Law and Starling’s equation) formulations. The IE model has been used to generate the training dataset for identifying the NN model architecture and parameters. The computation of the output LN drainage function characteristics (the fluid flow parameters and the exchange with blood) with the trained NN model required about 1000-fold less central processing unit (CPU) time than computationally tracing the flow characteristics of interest with the physics-based IE model. The use of the presented computational models will allow for a more realistic description and prediction of the immune cell circulation, cytokine distribution and drug pharmacokinetics in humans under various health and disease states as well as assisting in the development of artificial LN-on-a-chip technologies.

An update on oral drug delivery via intestinal lymphatic transport

Orally administered drug entities have to survive the harsh gastrointestinal environment, penetrate the enteric epithelia and circumvent hepatic metabolism before reaching the systemic circulation. Whereas the gastrointestinal stability can be well maintained by taking proper measures, hepatic metabolism presents as a formidable barrier to drugs suffering from first-pass metabolism. The pharmaceutical academia and industries are seeking alternative pathways for drug transport to circumvent problems associated with the portal pathway. Intestinal lymphatic transport is emerging as a promising pathway to this end. In this review, we intend to provide an updated overview on the rationale, strategies, factors and applications involved in intestinal lymphatic transport. There are mainly two pathways for peroral lymphatic transport-the chylomicron and the microfold cell pathways. The underlying mechanisms are being unraveled gradually and nowadays witness increasing research input and applications.

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.

Microfluidic valvular chips and a numerical lymphatic vessel model for the study of lymph transport characteristics

Lymph transport inside lymphatic vessels is highly complex and not yet fully understood. So far, a consensus has not been reached among existing analytical models on how spatiotemporal coordination of contracting adjacent lymphangions affects lymph transport. To understand complex lymph transport, we created a novel microfluidic valvular chip with flexible bicuspid valves and segmental pneumatic pumps based on a microfluidic device with an inside 3D structure made of hydrogels. Inside the chip, water moved unidirectionally when the microfluidic channel was locally compressed, with its velocity profile closely resembling the waveform of lymph observed in vivo. Furthermore, for a systematic and mechanistic study, we constructed a numerical model based on fluid-structure interaction and validated the model via demonstration of similarities in water transport characteristics between the model and the chip. Using this model, we examined various mechanical and time-dependent parameters, such as period, phase delay, sequence, and strength of contractions, valve compliance, fluid viscosity, and pressure differences, for their effects on water transport. Although our model is simplified, it enabled a parametric study that helped clarify the mechano-temporal correlations between compressions of adjacent chambers via transmissions of hydrodynamic forces, which regulate complex lymph transport. Moreover, our chip demonstrated technical advances that enable unidirectional discrete movement of fluid in the picoliter range by phenumatic pumping. The velocity profile is also similar to the pulse waveform of arteries under pathological conditions such as increased aortic stiffness, allowing our chip to be used for in vitro mechanobiology studies of endothelial cells.

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.

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.

A novel mouse tail lymphedema model for observing lymphatic pump failure during lymphedema development

It has been suggested that many forms of secondary lymphedema in humans are driven by a progressive loss of lymphatic pump function after an initial risk-inducing event. However, the link between pump failure and disease progression has remained elusive due to experimental challenges in the clinical setting and a lack of adequate animal models. Using a novel surgical model of lymphatic injury, we track the adaptation and functional decline of the lymphatic network in response to surgery. This model mimics the histological hallmarks of the typical mouse tail lymphedema model while leaving an intact collecting vessel for analysis of functional changes during disease progression. Lymphatic function in the intact collecting vessel negatively correlated with swelling, while a loss of pumping pressure generation remained even after resolution of swelling. By using this model to study the role of obesity in lymphedema development, we show that obesity exacerbates acquired lymphatic pump failure following lymphatic injury, suggesting one mechanism through which obesity may worsen lymphedema. This lymphatic injury model will allow for future studies investigating the molecular mechanisms leading to lymphedema development.

Renal Lymphatics: Anatomy, Physiology, and Clinical Implications

Renal lymphatics are abundant in the cortex of the normal kidney but have been largely neglected in discussions around renal diseases. They originate in the substance of the renal lobule as blind-ended initial capillaries, and can either follow the main arteries and veins toward the hilum, or penetrate the capsule to join capsular lymphatics. There are no valves present in interlobular lymphatics, which allows lymph formed in the cortex to exit the kidney in either direction. There are very few lymphatics present in the medulla. Lymph is formed from interstitial fluid in the cortex, and is largely composed of capillary filtrate, but also contains fluid reabsorbed from the tubules. The two main factors that contribute to renal lymph formation are interstitial fluid volume and intra-renal venous pressure. Renal lymphatic dysfunction, defined as a failure of renal lymphatics to adequately drain interstitial fluid, can occur by several mechanisms. Renal lymphatic inflow may be overwhelmed in the setting of raised venous pressure (e.g., cardiac failure) or increased capillary permeability (e.g., systemic inflammatory response syndrome). Similarly, renal lymphatic outflow, at the level of the terminal thoracic duct, may be impaired by raised central venous pressures. Renal lymphatic dysfunction, from any cause, results in renal interstitial edema. Beyond a certain point of edema, intra-renal collecting lymphatics may collapse, further impairing lymphatic drainage. Additionally, in an edematous, tense kidney, lymphatic vessels exiting the kidney via the capsule may become blocked at the exit point. The reciprocal negative influences between renal lymphatic dysfunction and renal interstitial edema are expected to decrease renal function due to pressure changes within the encapsulated kidney, and this mechanism may be important in several common renal conditions.