Loss of anoctamin 1 reveals a subtle role for BK channels in lymphatic muscle action potentials

Ca2+ signalling plays a crucial role in determining lymphatic muscle cell excitability and contractility through its interaction with the Ca2+-activated Cl- channel anoctamin 1 (ANO1). In contrast, the large-conductance (BK) Ca2+-activated K+ channel (KCa) and other KCa channels have prominent vasodilatory actions by hyperpolarizing vascular smooth muscle cells. Here, we assessed the expression and contribution of the KCa family to mouse and rat lymphatic collecting vessel contractile function. The BK channel was the only KCa channel consistently expressed in fluorescence-activated cell sorting-purified mouse lymphatic muscle cell lymphatic muscle cells. We used a pharmacological inhibitor of BK channels, iberiotoxin, and small-conductance Ca2+-activated K+ channels, apamin, to inhibit KCa channels acutely in ex vivo isobaric myography experiments and intracellular membrane potential recordings. In basal conditions, BK channel inhibition had little to no effect on either mouse inguinal-axillary lymphatic vessel (MIALV) or rat mesenteric lymphatic vessel contractions or action potentials (APs). We also tested BK channel inhibition under loss of ANO1 either by genetic ablation (Myh11CreERT2-Ano1 fl/fl, Ano1ismKO) or by pharmacological inhibition with Ani9. In both Ano1ismKO MIALVs and Ani9-pretreated MIALVs, inhibition of BK channels increased contraction amplitude, increased peak AP and broadened the peak of the AP spike. In rat mesenteric lymphatic vessels, BK channel inhibition also abolished the characteristic post-spike notch, which was exaggerated with ANO1 inhibition, and significantly increased the peak potential and broadened the AP spike. We conclude that BK channels are present and functional on mouse and rat lymphatic muscle cells but are otherwise masked by the dominance of ANO1. KEY POINTS: Mouse and rat lymphatic muscle cells express functional BK channels. BK channels make little contribution to either rat or mouse lymphatic collecting vessel contractile function in basal conditions across a physiological pressure range. ANO1 limits the peak membrane potential achieved in the action potential and sets a plateau potential limiting the voltage-dependent activation of BK. BK channels are activated when ANO1 is absent or blocked and slightly impair contractile strength by reducing the peak membrane potential achieved in the action potential spike and accelerating the post-spike repolarization.

Dichotomous effects of in vivo and in vitro ionizing radiation exposure on lymphatic function

On the one hand, lymphatic dysfunction induces interstitial edema and inflammation. On the other hand, the formation of edema and inflammation induce lymphatic dysfunction. However, informed by the earlier reports of undetected apoptosis of irradiated lymphatic endothelial cells (LECs) in vivo, lymphatic vessels are commonly considered inconsequential to ionizing radiation (IR)-induced inflammatory injury to normal tissues. Primarily because of the lack of understanding of the acute effects of IR exposure on lymphatic function, acute edema and inflammation, common sequelae of IR exposure, have been ascribed solely to blood vessel damage. Therefore, in the present study, the lymphatic acute responses to IR exposure were quantified to evaluate the hypothesis that IR exposure impairs lymphatic pumping. Rat mesenteric lymphatic vessels were irradiated in vivo or in vitro, and changes in pumping were quantified in isolated vessels in vitro. Compared with sham-treated vessels, pumping was lowered in lymphatic vessels irradiated in vivo but increased in vessels irradiated in vitro. Furthermore, unlike in blood vessels, the acute effects of IR exposure in lymphatic vessels were not mediated by nitric oxide-dependent pathways in either in vivo or in vitro irradiated vessels. After cyclooxygenase blockade, pumping was partially restored in lymphatic vessels irradiated in vitro but not in vessels irradiated in vivo. Taken together, these findings demonstrated that lymphatic vessels are radiosensitive and LEC apoptosis alone may not account for all the effects of IR exposure on the lymphatic system.NEW & NOTEWORTHY Earlier studies leading to the common belief that lymphatic vessels are radioresistant either did not characterize lymphatic pumping, deemed necessary for the resolution of edema and inflammation, or did it in vivo. By characterizing pumping in vitro, the present study, for the first time, demonstrated that lymphatic pumping was impaired in vessels irradiated in vivo and enhanced in vessels irradiated in vitro. Furthermore, the pathways implicated in ionizing radiation-induced blood vessel damage did not mediate lymphatic responses.

Direct Interstitial Decongestion in an Animal Model of Acute-on-Chronic Ischemic Heart Failure

•

In ADHF, elevated CVP opposes thoracic duct lymph flow and impairs decongestion of the interstitial space.

•

The use of a novel device for reducing CVP at the outflow of the thoracic duct was shown to be safe, well-tolerated, and effectively reduced EVLW, in an animal model of acute-on-chronic ischemic HF.

•

Similar results were observed when translating this therapy to a human case study.

•

Additional human studies to confirm these findings may establish device-based direct interstitial decongestion as a new treatment for ADHF.

Removal of excess fluid in acute decompensated heart failure (ADHF) targets the intravascular space, whereas most fluid resides in the interstitial space. The authors evaluated an approach to interstitial decongestion using a device to enhance lymph flow. The device was deployed in sheep with induced heart failure (HF) and acute volume overload to create a low-pressure zone at the thoracic duct outlet. Treatment decreased extravascular lung water (EVLW) volume (mL/kg) (-32% ± 9%, P = 0.029) compared to controls (+46% ± 9%, P = 0.003). Device-mediated thoracic duct decompression effectively reduced EVLW. Human studies may establish device-based interstitial decongestion as a new ADHF treatment.

Dichotomous effects on lymphatic transport with loss of caveolae in mice

AIM

Fluid and macromolecule transport from the interstitium into and through lymphatic vessels is necessary for tissue homeostasis. While lymphatic capillary structure suggests that passive, paracellular transport would be the predominant route of macromolecule entry, active caveolae-mediated transcellular transport has been identified in lymphatic endothelial cells (LECs) in vitro. Caveolae also mediate a wide array of endothelial cell processes, including nitric oxide regulation. Thus, how does the lack of caveolae impact "lymphatic function"?

METHODS

Various aspects of lymphatic transport were measured in mice constitutively lacking caveolin-1 ("CavKO"), the protein required for caveolae formation in endothelial cells, and in mice with a LEC-specific Cav1 gene deletion (Lyve1-Cre x Cav1^flox/flox ; "LyCav") and ex vivo in their vessels and cells.

RESULTS

In each model, lymphatic architecture was largely unchanged. The lymphatic conductance, or initial tissue uptake, was significantly higher in both CavKO mice and LyCav mice by quantitative microlymphangiography and the permeability to 70 kDa dextran was significantly increased in monolayers of LECs isolated from CavKO mice. Conversely, transport within the lymphatic system to the sentinel node was significantly reduced in anaesthetized CavKO and LyCav mice. Isolated, cannulated collecting vessel studies identified significantly reduced phasic contractility when lymphatic endothelium lacks caveolae. Inhibition of nitric oxide synthase was able to partially restore ex vivo vessel contractility.

CONCLUSION

Macromolecule transport across lymphatics is increased with loss of caveolae, yet phasic contractility reduced, resulting in reduced overall lymphatic transport function. These studies identify lymphatic caveolar biology as a key regulator of active lymphatic transport functions.

Adaptation of the hepatic transudation barrier to sinusoidal hypertension

The role of the hepatic transudation barrier in determining ascites volume and protein content in chronic liver disease is poorly understood. Therefore, the purpose of the present study was to characterize how chronic sinusoidal hypertension impacts hepatic transudation barrier properties and the transudation rate. The suprahepatic inferior vena cava was surgically constricted, and animals were exposed to either short-term (SVH; 2-3 wk) or long-term venous hypertension (LVH; 5-6 wk). Compared with SVH, LVH resulted in lower peritoneal fluid pressure, ascites volume, and ascites protein concentration. The transudation barrier protein reflection coefficient was significantly higher, and the transudation barrier hydraulic conductivity, transudation rate, and transudate-to-lymph protein concentration ratio were significantly lower in LVH animals compared with SVH animals. The sensitivity of transudation rates to acute changes in interstitial fluid pressures was also significantly lower in LVH animals compared with SVH animals. In contrast, there was no detectable difference in hepatic lymph flow rate or sensitivity of lymph flow to acute changes in interstitial fluid pressures between SVH and LVH animals. Taken together, these data suggest that decreased hepatic transudation barrier permeability to fluid and protein and increased reflection coefficient led to a decrease in the hepatic contribution to ascites volume. The present work, to the best of our knowledge, is the first to quantify an anti-ascites adaptation of the hepatic transudation barrier in response to chronic hepatic sinusoidal hypertension.

Hepatic transudation barrier properties

OBJECTIVE

Fluid and protein continuously transude from the surface of the liver. Despite a common understanding that transudation plays a critical role in hepatic interstitial and peritoneal fluid balance, transudation from the entire liver has not been studied. Therefore, the goal of the present work was to provide the first direct measurement of the hepatic transudation rate and transudation barrier properties.

METHODS

Transudation rates were determined by collecting transudate from the entire liver. Hydraulic conductivity, and fluid transudation and protein reflection coefficients of the transudation barrier (formed by the subscapular interstitial matrix, capsule, and peritoneum) were determined from changes in fluid and protein transudation rates in response to hepatic venous pressure elevation.

RESULTS

Following hepatic venous pressure elevation from 6.1 ± 0.9 to 11.1 ± 0.6 mm Hg, transudation rate increased from 0.13 ± 0.03 to 0.37 ± 0.03 mL/min·100 g. Transudation barrier hydraulic conductivity, fluid transudation and protein reflection coefficients (3.9 × 10^-4  ± 5.7 × 10^-5  mL/min·mm Hg·cm² , 0.36 ± 0.04 mL/min·mm Hg, and 0.09 ± 0.03, respectively) were comparable to those reported for hepatic sinusoids.

CONCLUSIONS

Taken together, these findings suggest that the hepatic transudation barrier is highly permeable at elevated sinusoidal pressures. These fundamental studies provide a better understanding of the hepatic transudation barrier properties and transudation under conditions that are physiologically and clinically relevant to ascites formation.

Mesenteric lymphatic vessels adapt to mesenteric venous hypertension by becoming weaker pumps

Lymphangions, the segments of lymphatic vessels between two adjacent lymphatic valves, actively pump lymph. Acute changes in transmural pressure and lymph flow have profound effects on lymphatic pump function in vitro. Chronic changes in pressure and flow in vivo have also been reported to lead to significant changes in lymphangion function. Because changes in pressure and flow are both cause and effect of adaptive processes, characterizing adaptation requires a more fundamental analysis of lymphatic muscle properties. Therefore, the purpose of the present work was to use an intact lymphangion isovolumetric preparation to evaluate changes in mesenteric lymphatic muscle mechanical properties and the intracellular Ca(2+) in response to sustained mesenteric venous hypertension. Bovine mesenteric veins were surgically occluded to create mesenteric venous hypertension. Postnodal mesenteric lymphatic vessels from mesenteric venous hypertension (MVH; n = 6) and sham surgery (Sham; n = 6) animals were isolated and evaluated 3 days after the surgery. Spontaneously contracting MVH vessels generated end-systolic active tension and end-diastolic active tension lower than the Sham vessels. Furthermore, steady-state active tension and intracellular Ca(2+) concentration levels in response to KCl stimulation were also significantly lower in MVH vessels compared with those of the Sham vessels. There was no significant difference in passive tension in lymphatic vessels from the two groups. Taken together, these results suggest that following 3 days of mesenteric venous hypertension, postnodal mesenteric lymphatic vessels adapt to become weaker pumps with decreased cytosolic Ca(2+) concentration.

Functional adaptation of bovine mesenteric lymphatic vessels to mesenteric venous hypertension

Lymph flow is the primary mechanism for returning interstitial fluid to the blood circulation. Currently, the adaptive response of lymphatic vessels to mesenteric venous hypertension is not known. This study sought to determine the functional responses of postnodal mesenteric lymphatic vessels. We surgically occluded bovine mesenteric veins to create mesenteric venous hypertension to elevate mesenteric lymph flow. Three days after surgery, postnodal mesenteric lymphatic vessels from mesenteric venous hypertension (MVH; n = 7) and sham surgery (Sham; n = 6) group animals were evaluated and compared. Contraction frequency (MVH: 2.98 ± 0.75 min(-1); Sham: 5.42 ± 0.81 min(-1)) and fractional pump flow (MVH: 1.14 ± 0.30 min(-1); Sham: 2.39 ± 0.32 min(-1)) were significantly lower in the venous occlusion group. These results indicate that postnodal mesenteric lymphatic vessels adapt to mesenteric venous hypertension by reducing intrinsic contractile activity.

Adaptation of mesenteric lymphatic vessels to prolonged changes in transmural pressure

In vitro studies have revealed that acute increases in transmural pressure increase lymphatic vessel contractile function. However, adaptive responses to prolonged changes in transmural pressure in vivo have not been reported. Therefore, we developed a novel bovine mesenteric lymphatic partial constriction model to test the hypothesis that lymphatic vessels exposed to higher transmural pressures adapt functionally to become stronger pumps than vessels exposed to lower transmural pressures. Postnodal mesenteric lymphatic vessels were partially constricted for 3 days. On postoperative day 3, constricted vessels were isolated, and divided into upstream (UP) and downstream (DN) segment groups, and instrumented in an isolated bath. Although there were no differences between the passive diameters of the two groups, both diastolic diameter and systolic diameter were significantly larger in the UP group than in the DN group. The pump index of the UP group was also higher than that in the DN group. In conclusion, this is the first work to report how lymphatic vessels adapt to prolonged changes in transmural pressure in vivo. Our results suggest that vessel segments upstream of the constriction adapt to become both better fluid conduits and lymphatic pumps than downstream segments.

Award article: Microcirculatory Society Award for Excellence in Lymphatic Research: time course of myocardial interstitial edema resolution and associated left ventricular dysfunction

OBJECTIVE

Although the causal relationship between acute myocardial edema and cardiac dysfunction has been established, resolution of myocardial edema and subsequent recovery of cardiac function have not been established. The time to resolve myocardial edema and the degree that cardiac function is depressed after edema resolves are not known. We therefore characterized temporal changes in cardiac function as acute myocardial edema formed and resolved.

METHODS

Acute myocardial edema was induced in the canine model by elevating coronary sinus pressure for three hours. Myocardial water content and cardiac function were determined before and during coronary sinus pressure elevation, and after coronary sinus pressure restoration.

RESULTS

Although no change in systolic properties was detected, accumulation of water in myocardial interstitium was associated with increased diastolic stiffness. When coronary sinus pressure was relieved, myocardial edema resolved within 180 minutes. Diastolic stiffness, however, remained significantly elevated compared with baseline values, and cardiac function remained compromised.

CONCLUSIONS

The present work suggests that the cardiac dysfunction caused by the formation of myocardial edema may persist after myocardial edema resolves. With the advent of new imaging techniques to quantify myocardial edema, this insight provides a new avenue for research to detect and treat a significant cause of cardiac dysfunction.

Blood flow augmentation by intrinsic venular contraction in vivo

Venomotion, spontaneous cyclic contractions of venules, was first observed in the bat wing 160 years ago. Of all the functional roles proposed since then, propulsion of blood by venomotion remains the most controversial. Common animal models that require anesthesia and surgery have failed to provide evidence for venular pumping of blood. To determine whether venomotion actively pumps blood in a minimally invasive, unanesthetized animal model, we reintroduced the batwing model. We evaluated the temporal and functional relationship between the venous contraction cycle and blood flow and luminal pressure. Furthermore, we determined the effect of inhibiting venomotion on blood flow. We found that the active venous contractions produced an increase in the blood flow and exhibited temporal vessel diameter-blood velocity and pressure relationships characteristic of a peristaltic pump. The presence of valves, a characteristic of reciprocating pumps, enhances the efficiency of the venular peristaltic pump by preventing retrograde flow. Instead of increasing blood flow by decreasing passive resistance, venular dilation with locally applied sodium nitroprusside decreased blood flow. Taken together, these observations provide evidence for active venular pumping of blood. Although strong venomotion may be unique to bats, venomotion has also been inferred from venous pressure oscillations in other animal models. The conventional paradigm of microvascular pressure and flow regulation assumes venules only act as passive resistors, a proposition that must be reevaluated in the presence of significant venomotion.

Evaluation of gravimetric techniques to estimate the microvascular filtration coefficient

Microvascular permeability to water is characterized by the microvascular filtration coefficient (K(f)). Conventional gravimetric techniques to estimate K(f) rely on data obtained from either transient or steady-state increases in organ weight in response to increases in microvascular pressure. Both techniques result in considerably different estimates and neither account for interstitial fluid storage and lymphatic return. We therefore developed a theoretical framework to evaluate K(f) estimation techniques by 1) comparing conventional techniques to a novel technique that includes effects of interstitial fluid storage and lymphatic return, 2) evaluating the ability of conventional techniques to reproduce K(f) from simulated gravimetric data generated by a realistic interstitial fluid balance model, 3) analyzing new data collected from rat intestine, and 4) analyzing previously reported data. These approaches revealed that the steady-state gravimetric technique yields estimates that are not directly related to K(f) and are in some cases directly proportional to interstitial compliance. However, the transient gravimetric technique yields accurate estimates in some organs, because the typical experimental duration minimizes the effects of interstitial fluid storage and lymphatic return. Furthermore, our analytical framework reveals that the supposed requirement of tying off all draining lymphatic vessels for the transient technique is unnecessary. Finally, our numerical simulations indicate that our comprehensive technique accurately reproduces the value of K(f) in all organs, is not confounded by interstitial storage and lymphatic return, and provides corroboration of the estimate from the transient technique.

Myocardial microvascular permeability, interstitial oedema, and compromised cardiac function

The heart, perhaps more than any other organ, is exquisitely sensitive to increases in microvascular permeability and the accumulation of myocardial interstitial oedema fluid. Whereas some organs can cope with profound increases in the interstitial fluid volume or oedema formation without a compromise in function, heart function is significantly compromised with only a few percent increase in the interstitial fluid volume. This would be of little consequence if myocardial oedema were an uncommon pathology. On the contrary, myocardial oedema forms in response to many disease states as well as clinical interventions such as cardiopulmonary bypass and cardioplegic arrest common to many cardiothoracic surgical procedures. The heart's inability to function effectively in the presence of myocardial oedema is further confounded by the perplexing fact that the resolution of myocardial oedema does not restore normal cardiac function. We will attempt to provide some insight as to how microvascular permeability and myocardial oedema formation compromise cardiac function and discuss the acute changes that might take place in the myocardium to perpetuate compromised cardiac function following oedema resolution. We will also discuss compensatory changes in the interstitial matrix of the heart in response to chronic myocardial oedema and the role they play to optimize myocardial function during chronic oedemagenic disease.

Balance point characterization of interstitial fluid volume regulation

The individual processes involved in interstitial fluid volume and protein regulation (microvascular filtration, lymphatic return, and interstitial storage) are relatively simple, yet their interaction is exceedingly complex. There is a notable lack of a first-order, algebraic formula that relates interstitial fluid pressure and protein to critical parameters commonly used to characterize the movement of interstitial fluid and protein. Therefore, the purpose of the present study is to develop a simple, transparent, and general algebraic approach that predicts interstitial fluid pressure (P(i)) and protein concentrations (C(i)) that takes into consideration all three processes. Eight standard equations characterizing fluid and protein flux were solved simultaneously to yield algebraic equations for P(i) and C(i) as functions of parameters characterizing microvascular, interstitial, and lymphatic function. Equilibrium values of P(i) and C(i) arise as balance points from the graphical intersection of transmicrovascular and lymph flows (analogous to Guyton's classical cardiac output-venous return curves). This approach goes beyond describing interstitial fluid balance in terms of conservation of mass by introducing the concept of inflow and outflow resistances. Algebraic solutions demonstrate that P(i) and C(i) result from a ratio of the microvascular filtration coefficient (1/inflow resistance) and effective lymphatic resistance (outflow resistance), and P(i) is unaffected by interstitial compliance. These simple algebraic solutions predict P(i) and C(i) that are consistent with reported measurements. The present work therefore presents a simple, transparent, and general balance point characterization of interstitial fluid balance resulting from the interaction of microvascular, interstitial, and lymphatic function.

Venomotion modulates lymphatic pumping in the bat wing

In skin, it is believed that lymph must be pumped by intrinsic contraction of lymphatic muscle, since investigators have not considered that cyclical dilation of venules could compress adjacent lymphatic microvessels. Because lymphatic vessels are sensitive to stretch, we hypothesized that venomotion not only can cause extrinsic pumping of lymph in nearby lymphatic vessels, but also can stimulate intrinsic contractions. Bat wing venules have pronounced venomotion and are in close proximity to lymphatic microvessels, and can be studied noninvasively without the confounding effects of anesthesia, surgical trauma, or contrast agents. Therefore, the interaction between venules and their paired lymphatic vessels in unanesthetized Pallid bats (n = 8) was evaluated by recording the diameters of both vessels. Four sets of observations suggested that lymphatic and venous contractions were partially coupled. First, venous dilation and contraction produced a significant change in lymphatic microvascular cross-sectional area. Second, lymphatic microvascular contractions were immediately preceded by a change in venular diameter. Third, venular and lymphatic vessel contraction frequencies were positively correlated (r = 0.75). Fourth, time delays between peak venular systole and onset of lymphatic microvascular contraction were negatively correlated with venomotion magnitude (r = -0.55) and velocity (r = -0.64). In a separate experiment, inhibiting venomotion resulted in a 54.3 +/- 20.0% (n = 8) decrease in lymphatic contraction frequency. Furthermore, 85.7% (n = 56) of lymphatic vessels switch sides and lie adjacent to arterioles when venules were too small to exhibit venomotion. These results are consistent with both extrinsic pumping of lymph and stretch-induced lymphatic contraction and imply that intrinsic and extrinsic pumping can be coupled.

First-order approximation for the pressure-flow relationship of spontaneously contracting lymphangions

To return lymph to the great veins of the neck, it must be actively pumped against a pressure gradient. Mean lymph flow in a portion of a lymphatic network has been characterized by an empirical relationship (Pin − Pout = −Pp + RLQL), where Pin − Pout is the axial pressure gradient and QL is mean lymph flow. RL and Pp are empirical parameters characterizing the effective lymphatic resistance and pump pressure, respectively. The relation of these global empirical parameters to the properties of lymphangions, the segments of a lymphatic vessel bounded by valves, has been problematic. Lymphangions have a structure like blood vessels but cyclically contract like cardiac ventricles; they are characterized by a contraction frequency ( f) and the slopes of the end-diastolic pressure-volume relationship [minimum value of resulting elastance ( Emin)] and end-systolic pressure-volume relationship [maximum value of resulting elastance ( Emax)]. Poiseuille's law provides a first-order approximation relating the pressure-flow relationship to the fundamental properties of a blood vessel. No analogous formula exists for a pumping lymphangion. We therefore derived an algebraic formula predicting lymphangion flow from fundamental physical principles and known lymphangion properties. Quantitative analysis revealed that lymph inertia and resistance to lymph flow are negligible and that lymphangions act like a series of interconnected ventricles. For a single lymphangion, Pp = Pin ( Emax − Emin)/ Emin and RL = Emax/ f. The formula was tested against a validated, realistic mathematical model of a lymphangion and found to be accurate. Predicted flows were within the range of flows measured in vitro. The present work therefore provides a general solution that makes it possible to relate fundamental lymphangion properties to lymphatic system function.

Edemagenic gain and interstitial fluid volume regulation

Under physiological conditions, interstitial fluid volume is tightly regulated by balancing microvascular filtration and lymphatic return to the central venous circulation. Even though microvascular filtration and lymphatic return are governed by conservation of mass, their interaction can result in exceedingly complex behavior. Without making simplifying assumptions, investigators must solve the fluid balance equations numerically, which limits the generality of the results. We thus made critical simplifying assumptions to develop a simple solution to the standard fluid balance equations that is expressed as an algebraic formula. Using a classical approach to describe systems with negative feedback, we formulated our solution as a "gain" relating the change in interstitial fluid volume to a change in effective microvascular driving pressure. The resulting "edemagenic gain" is a function of microvascular filtration coefficient (K(f)), effective lymphatic resistance (R(L)), and interstitial compliance (C). This formulation suggests two types of gain: "multivariate" dependent on C, R(L), and K(f), and "compliance-dominated" approximately equal to C. The latter forms a basis of a novel method to estimate C without measuring interstitial fluid pressure. Data from ovine experiments illustrate how edemagenic gain is altered with pulmonary edema induced by venous hypertension, histamine, and endotoxin. Reformulation of the classical equations governing fluid balance in terms of edemagenic gain thus yields new insight into the factors affecting an organ's susceptibility to edema.

Lymphangion coordination minimally affects mean flow in lymphatic vessels

The lymphatic system returns interstitial fluid to the central venous circulation, in part, by the cyclical contraction of a series of “lymphangion pumps” in a lymphatic vessel. The dynamics of individual lymphangions have been well characterized in vitro; their frequencies and strengths of contraction are sensitive to both preload and afterload. However, lymphangion interaction within a lymphatic vessel has been poorly characterized because it is difficult to experimentally alter properties of individual lymphangions and because the afterload of one lymphangion is coupled to the preload of another. To determine the effects of lymphangion interaction on lymph flow, we adapted an existing mathematical model of a lymphangion (characterizing lymphangion contractility, lymph viscosity, and inertia) to create a new lymphatic vessel model consisting of several lymphangions in series. The lymphatic vessel model was validated with focused experiments on bovine mesenteric lymphatic vessels in vitro. The model was then used to predict changes in lymph flow with different time delays between onset of contraction of adjacent lymphangions (coordinated case) and with different relative lymphangion contraction frequencies (noncoordinated case). Coordination of contraction had little impact on mean flow. Furthermore, orthograde and retrograde propagations of contractile waves had similar effects on flow. Model results explain why neither retrograde propagation of contractile waves nor the lack of electrical continuity between lymphangions adversely impacts flow. Because lymphangion coordination minimally affects mean flow in lymphatic vessels, lymphangions have flexibility to independently adapt to local conditions.