Traffic Patterns of the Migrating Endothelium: How Force Transmission Regulates Vascular Malformation and Functional Shunting During Angiogenic Remodelling.

Angiogenesis occurs in distinct phases: initial spouting is followed by remodelling in which endothelial cells (ECs) composing blood vessels rearrange by migrating against the direction of flow. Abnormal remodelling can result in vascular malformation. Such is the case in mutation of the Alk1 receptor within the mouse retina which disrupts flow-migration coupling, creating mixed populations of ECs polarised with/against flow which aggregate into arteriovenous malformations (AVMs). The lack of live imaging options in vivo means that the collective EC dynamics that drive AVM and the consequences of mixed populations of polarity remain a mystery. Therefore, our goal is to present a novel agent-based model to provide theoretical insight into EC force transmission and collective dynamics during angiogenic remodelling. Force transmission between neighbouring agents consists of extrusive forces which maintain spacing and cohesive forces which maintain the collective. We performed migration simulations within uniformly polarised populations (against flow) and mixed polarity (with/against flow). Within uniformly polarised populations, extrusive forces stabilised the plexus by facilitating EC intercalation which ensures that cells remained evenly distributed. Excess cohesion disrupts intercalation, resulting in aggregations of cells and functional shunting. Excess cohesion between ECs prevents them from resolving diameter balances within the plexus, leading to prolonged flow reversals which exert a critical behaviour change within the system as they switch the direction of cell migration and traffic patterns at bifurcations. Introducing mixtures of cell polarity dramatically changed the role of extrusive forces within the system. At low extrusion, opposing ECs were able to move past each other; however, at high extrusion the pushing between cells resulted in migration speeds close to zero, forming traffic jams and disrupting migration. In our study, we produced vascular malformations and functional shunting with either excess cohesion between ECs or mixtures of cell polarity. At the centre of both these mechanisms are cell-cell adherens junctions, which are involved in flow sensing/polarity and must remodelling dynamically to allow rearrangements of cells during vascular patterning. Thus, our findings implicate junctional dysfunction as a new target in the treatment and prevention of vascular disease and AVMs.

WASp controls oriented migration of endothelial cells to achieve functional vascular patterning.

Endothelial cell migration and proliferation are essential for the establishment of a hierarchical organization of blood vessels and optimal distribution of blood. However, how these cellular processes are quantitatively coordinated to drive vascular network morphogenesis remains unknown. Here, using the zebrafish vasculature as a model system, we demonstrate that the balanced distribution of endothelial cells, as well as the resulting regularity of vessel calibre, is a result of cell migration from veins towards arteries and cell proliferation in veins. We identify the Wiskott-Aldrich Syndrome protein (WASp) as an important molecular regulator of this process and show that loss of coordinated migration from veins to arteries upon wasb depletion results in aberrant vessel morphology and the formation of persistent arteriovenous shunts. We demonstrate that WASp achieves its function through the coordination of junctional actin assembly and PECAM1 recruitment and provide evidence that this is conserved in humans. Overall, we demonstrate that functional vascular patterning in the zebrafish trunk is established through differential cell migration regulated by junctional actin, and that interruption of differential migration may represent a pathomechanism in vascular malformations.

Lymph node swelling combined with temporary effector T cell retention aids T cell response in a model of adaptive immunity.

Swelling of lymph nodes (LNs) is commonly observed during the adaptive immune response, yet the impact on T cell (TC) trafficking and subsequent immune response is not well known. To better understand the effect of macro-scale alterations, we developed an agent-based model of the LN paracortex, describing the TC proliferative response to antigen-presenting dendritic cells alongside inflammation-driven and swelling-induced changes in TC recruitment and egress, while also incorporating regulation of the expression of egress-modulating TC receptor sphingosine-1-phosphate receptor-1. Analysis of the effector TC response under varying swelling conditions showed that swelling consistently aided TC activation. However, subsequent effector CD8+ TC production was reduced in scenarios where swelling occurred too early in the TC proliferative phase or when TC cognate frequency was low due to increased opportunity for TC exit. Temporarily extending retention of newly differentiated effector TCs, mediated by sphingosine-1-phosphate receptor-1 expression, mitigated any negative effects of swelling by allowing facilitation of activation to outweigh increased access to exit areas. These results suggest that targeting temporary effector TC retention and egress associated with swelling offers new ways to modulate effector TC responses in, for example, immuno-suppressed patients and to optimize of vaccine design.

Association between erythrocyte dynamics and vessel remodelling in developmental vascular networks.

Sprouting angiogenesis is an essential vascularization mechanism consisting of sprouting and remodelling. The remodelling phase is driven by rearrangements of endothelial cells (ECs) within the post-sprouting vascular plexus. Prior work has uncovered how ECs polarize and migrate in response to flow-induced wall shear stress (WSS). However, the question of how the presence of erythrocytes (widely known as red blood cells (RBCs)) and their impact on haemodynamics affect vascular remodelling remains unanswered. Here, we devise a computational framework to model cellular blood flow in developmental mouse retina. We demonstrate a previously unreported highly heterogeneous distribution of RBCs in primitive vasculature. Furthermore, we report a strong association between vessel regression and RBC hypoperfusion, and identify plasma skimming as the driving mechanism. Live imaging in a developmental zebrafish model confirms this association. Taken together, our results indicate that RBC dynamics are fundamental to establishing the regional WSS differences driving vascular remodelling via their ability to modulate effective viscosity.

Reduced biomechanical models for precision-cut lung-slice stretching experiments.

Precision-cut lung-slices (PCLS), in which viable airways embedded within lung parenchyma are stretched or induced to contract, are a widely used ex vivo assay to investigate bronchoconstriction and, more recently, mechanical activation of pro-remodelling cytokines in asthmatic airways. We develop a nonlinear fibre-reinforced biomechanical model accounting for smooth muscle contraction and extracellular matrix strain-stiffening. Through numerical simulation, we describe the stresses and contractile responses of an airway within a PCLS of finite thickness, exposing the importance of smooth muscle contraction on the local stress state within the airway. We then consider two simplifying limits of the model (a membrane representation and an asymptotic reduction in the thin-PCLS-limit), that permit analytical progress. Comparison against numerical solution of the full problem shows that the asymptotic reduction successfully captures the key elements of the full model behaviour. The more tractable reduced model that we develop is suitable to be employed in investigations to elucidate the time-dependent feedback mechanisms linking airway mechanics and cytokine activation in asthma.

On the preservation of vessel bifurcations during flow-mediated angiogenic remodelling.

During developmental angiogenesis, endothelial cells respond to shear stress by migrating and remodelling the initially hyperbranched plexus, removing certain vessels whilst maintaining others. In this study, we argue that the key regulator of vessel preservation is cell decision behaviour at bifurcations. At flow-convergent bifurcations where migration paths diverge, cells must finely tune migration along both possible paths if the bifurcation is to persist. Experiments have demonstrated that disrupting the cells' ability to sense shear or the junction forces transmitted between cells impacts the preservation of bifurcations during the remodelling process. However, how these migratory cues integrate during cell decision making remains poorly understood. Therefore, we present the first agent-based model of endothelial cell flow-mediated migration suitable for interrogating the mechanisms behind bifurcation stability. The model simulates flow in a bifurcated vessel network composed of agents representing endothelial cells arranged into a lumen which migrate against flow. Upon approaching a bifurcation where more than one migration path exists, agents refer to a stochastic bifurcation rule which models the decision cells make as a combination of flow-based and collective-based migratory cues. With this rule, cells favour branches with relatively larger shear stress or cell number. We found that cells must integrate both cues nearly equally to maximise bifurcation stability. In simulations with stable bifurcations, we found competitive oscillations between flow and collective cues, and simulations that lost the bifurcation were unable to maintain these oscillations. The competition between these two cues is haemodynamic in origin, and demonstrates that a natural defence against bifurcation loss during remodelling exists: as vessel lumens narrow due to cell efflux, resistance to flow and shear stress increases, attracting new cells to enter and rescue the vessel from regression. Our work provides theoretical insight into the role of junction force transmission has in stabilising vasculature during remodelling and as an emergent mechanism to avoid functional shunting.

Validation of markerless strain-field optical tracking approach for soft tissue mechanical assessment.

Strain measurement during tissue deformation is crucial to elucidate relationships between mechanical loading and functional changes in biological tissues. When combined with specified loading conditions, assessment of strain fields can be used to craft models that accurately represent the mechanical behavior of soft tissue. Inhomogeneities in strain fields may be indicative of normal or pathological inhomogeneities in mechanical properties. In this study, we present the validation of a modified Demons registration algorithm for non-contact, marker-less strain measurement of tissue undergoing uniaxial loading. We validate the algorithm on a synthetic dataset composed of artificial deformation fields applied to a speckle image, as well as images of aortic sections of varying perceptual quality. Initial results indicate that Demons outperforms recent Optical Flow and Digital Image Correlation methods in terms of accuracy and robustness to low image quality, with similar runtimes. Demons achieves at least 8% lower maximal deviation from ground truth on 50% biaxial and shear strain applied to aortic images. To illustrate utility, we quantified strain fields of multiple human aortic specimens undergoing uniaxial tensile testing, noting the formation of strain concentrations in areas of rupture. The modified Demons algorithm captured a large range of strains (up to 50%) and provided spatially resolved strain fields that could be useful in the assessment of soft tissue pathologies.

A fully coupled fluid-structure interaction model of the secondary lymphatic valve.

The secondary lymphatic valve is a bi-leaflet structure frequent throughout collecting vessels that serves to prevent retrograde flow of lymph. Despite its vital function in lymph flow and apparent importance in disease development, the lymphatic valve and its associated fluid dynamics have been largely understudied. The goal of this work was to construct a physiologically relevant computational model of an idealized rat mesenteric lymphatic valve using fully coupled fluid-structure interactions to investigate the relationship between three-dimensional flow patterns and stress/deformation within the valve leaflets. The minimum valve resistance to flow, which has been shown to be an important parameter in effective lymphatic pumping, was computed as 268 g/mm4-s. Hysteretic behavior of the lymphatic valve was confirmed by comparing resistance values for a given transvalvular pressure drop during opening and closing. Furthermore, eddy structures were present within the sinus adjacent to the valve leaflets in what appear to be areas of vortical flow; the eddy structures were characterized by non-zero velocity values (up to ∼4 mm/s) in response to an applied unsteady transvalvular pressure. These modeling capabilities present a useful platform for investigating the complex interplay between soft tissue motion and fluid dynamics of lymphatic valves and contribute to the breadth of knowledge regarding the importance of biomechanics in lymphatic system function.

The passive biomechanics of human pelvic collecting lymphatic vessels.

The lymphatic system has a major significance in the metastatic pathways in women's cancers. Lymphatic pumping depends on both extrinsic and intrinsic mechanisms, and the mechanical behavior of lymphatic vessels regulates the function of the system. However, data on the mechanical properties and function of human lymphatics are lacking. Our aim is to characterize, for the first time, the passive biomechanical behavior of human collecting lymphatic vessels removed at pelvic lymph node dissection during primary debulking surgeries for epithelial ovarian cancer. Isolated vessels were cannulated and then pressurized at varying levels of applied axial stretch in a calcium-free Krebs buffer. Pressurized vessels were then imaged using multi-photon microscopy for collagen-elastin structural composition and fiber orientation. Both pressure-diameter and force-elongation responses were highly nonlinear, and axial stretching of the vessel served to decrease diameter at constant pressure. Pressure-diameter behavior for the human vessels is very similar to data from rat mesenteric vessels, though the human vessels were approximately 10× larger than those from rats. Multiphoton microscopy revealed the vessels to be composed of an inner layer of elastin with an outer layer of aligned collagen fibers. This is the first study that successfully described the passive biomechanical response and composition of human lymphatic vessels in patients with ovarian cancer. Future work should expand on this knowledge base with investigations of vessels from other anatomical locations, contractile behavior, and the implications on metastatic cell transport.

In Silico Investigation of Angiogenesis with Growth and Stress Generation Coupled to Local Extracellular Matrix Density.

Mechanical interactions during angiogenesis, i.e., traction applied by neovessels to the extracellular matrix and the corresponding deformation, are important regulators of growth and neovascularization. We have previously designed, implemented, and validated a coupled model of angiogenesis in which a discrete microvessel growth model interacts with a continuous finite element mesh through the application of local remodeling sprout stresses (Edgar et al. in Biomech Model Mechanobiol, 2014). However, the initial implementation of this framework does not take matrix density into account when determined these remodeling stresses and is therefore insufficient for the study of angiogenesis within heterogeneous matrix environments such as those found in vivo. The objective of this study was to implement sensitivity to matrix density in the active stress generation within AngioFE in order to allow the study of angiogenic growth within a heterogeneous density environment. We accomplished this by scaling active sprout stresses relative to local matrix density using a scaling factor previously determined from experimental data. We then exercised the new functionality of the model by simulating angiogenesis within four different scenarios: homogeneous density, a narrow gap model, and matrix density gradient, and a construct subjected to repeated loading/unloading and preconditioning. These numerical experiments predicted heterogeneous matrix density in the initially homogeneous case, the closure and alignment of microvessels along a low-density gap, the formation of a unique cap-like structure during angiogenesis within a density gradient, and the alignment of microvessels in the absence of applied load due to preconditioning. The result of these in silico investigations demonstrate how matrix heterogeneity affects neovascularization and matrix deformation and provides a platform for studying angiogenesis in complicated and multi-faceted mechanical environments that microvessels experience in vivo.

Large-scale time series microscopy of neovessel growth during angiogenesis.

During angiogenesis, growing neovessels must effectively navigate through the tissue space as they elongate and subsequently integrate into a microvascular network. While time series microscopy has provided insight into the cell activities within single growing neovessel sprouts, less is known concerning neovascular dynamics within a large angiogenic tissue bed. Here, we developed a time-lapse imaging technique that allowed visualization and quantification of sprouting neovessels as they form and grow away from adult parent microvessels in three dimensions over cubic millimeters of matrix volume during the course of up to 5 days on the microscope. Using a new image acquisition procedure and novel morphometric analysis tools, we quantified the elongation dynamics of growing neovessels and found an episodic growth pattern accompanied by fluctuations in neovessel diameter. Average elongation rate was 5 µm/h for individual vessels, but we also observed considerable dynamic variability in growth character including retraction and complete regression of entire neovessels. We observed neovessel-to-neovessel directed growth over tens to hundreds of microns preceding tip-to-tip inosculation. As we have previously described via static 3D imaging at discrete time points, we identified different collagen fibril structures associated with the growing neovessel tip and stalk, and observed the coordinated alignment of growing neovessels in a deforming matrix. Overall analysis of the entire image volumes demonstrated that although individual neovessels exhibited episodic growth and regression, there was a monotonic increase in parameters associated with the entire vascular bed such as total network length and number of branch points. This new time-lapse imaging approach corroborated morphometric changes in individual neovessels described by us and others, as well as captured dynamic neovessel behaviors unique to days-long angiogenesis within the forming neovascular network.

A coupled model of neovessel growth and matrix mechanics describes and predicts angiogenesis in vitro.

During angiogenesis, sprouting microvessels interact with the extracellular matrix (ECM) by degrading and reorganizing the matrix, applying traction forces, and producing deformation. Morphometric features of the resulting microvascular network are affected by the interaction between the matrix and angiogenic microvessels. The objective of this study was to develop a continuous-discrete modeling approach to simulate mechanical interactions between growing neovessels and the deformation of the matrix in vitro. This was accomplished by coupling an existing angiogenesis growth model which uses properties of the ECM to regulate angiogenic growth with the nonlinear finite element software FEBio (www.febio.org). FEBio solves for the deformation and remodeling of the matrix caused by active stress generated by neovessel sprouts, and this deformation was used to update the ECM into the current configuration. After mesh resolution and parameter sensitivity studies, the model was used to accurately predict vascular alignment for various matrix boundary conditions. Alignment primarily arises passively as microvessels convect with the deformation of the matrix, but active alignment along collagen fibrils plays a role as well. Predictions of alignment were most sensitive to the range over which active stresses were applied and the viscoelastic time constant in the material model. The computational framework provides a flexible platform for interpreting in vitro investigations of vessel-matrix interactions, predicting new experiments, and simulating conditions that are outside current experimental capabilities.

Cell-generated traction forces and the resulting matrix deformation modulate microvascular alignment and growth during angiogenesis.

The details of the mechanical factors that modulate angiogenesis remain poorly understood. Previous in vitro studies of angiogenesis using microvessel fragments cultured within collagen constructs demonstrated that neovessel alignment can be induced via mechanical constraint of the boundaries (i.e., boundary conditions). The objective of this study was to investigate the role of mechanical boundary conditions in the regulation of angiogenic alignment and growth in an in vitro model of angiogenesis. Angiogenic microvessels within three-dimensional constructs were subjected to different boundary conditions, thus producing different stress and strain fields during growth. Neovessel outgrowth and orientation were quantified from confocal image data after 6 days. Vascularity and branching decreased as the amount of constraint imposed on the culture increased. In long-axis constrained hexahedral constructs, microvessels aligned parallel to the constrained axis. In contrast, constructs that were constrained along the short axis had random microvessel orientation. Finite element models were used to simulate the contraction of gels under the various boundary conditions and to predict the local strain field experienced by microvessels. Results from the experiments and simulations demonstrated that microvessels aligned perpendicular to directions of compressive strain. Alignment was due to anisotropic deformation of the matrix from cell-generated traction forces interacting with the mechanical boundary conditions. These findings demonstrate that boundary conditions and thus the effective stiffness of the matrix regulate angiogenesis. This study offers a potential explanation for the oriented vascular beds that occur in native tissues and provides the basis for improved control of tissue vascularization in both native tissues and tissue-engineered constructs.

Extracellular matrix density regulates the rate of neovessel growth and branching in sprouting angiogenesis.

Angiogenesis is regulated by the local microenvironment, including the mechanical interactions between neovessel sprouts and the extracellular matrix (ECM). However, the mechanisms controlling the relationship of mechanical and biophysical properties of the ECM to neovessel growth during sprouting angiogenesis are just beginning to be understood. In this research, we characterized the relationship between matrix density and microvascular topology in an in vitro 3D organ culture model of sprouting angiogenesis. We used these results to design and calibrate a computational growth model to demonstrate how changes in individual neovessel behavior produce the changes in vascular topology that were observed experimentally. Vascularized gels with higher collagen densities produced neovasculatures with shorter vessel lengths, less branch points, and reduced network interconnectivity. The computational model was able to predict these experimental results by scaling the rates of neovessel growth and branching according to local matrix density. As a final demonstration of utility of the modeling framework, we used our growth model to predict several scenarios of practical interest that could not be investigated experimentally using the organ culture model. Increasing the density of the ECM significantly reduced angiogenesis and network formation within a 3D organ culture model of angiogenesis. Increasing the density of the matrix increases the stiffness of the ECM, changing how neovessels are able to deform and remodel their surroundings. The computational framework outlined in this study was capable of predicting this observed experimental behavior by adjusting neovessel growth rate and branching probability according to local ECM density, demonstrating that altering the stiffness of the ECM via increasing matrix density affects neovessel behavior, thereby regulated vascular topology during angiogenesis.

Mechanical interaction of angiogenic microvessels with the extracellular matrix.

Angiogenesis is the process by which new blood vessels sprout from existing blood vessels, enabling new vascular elements to be added to an existing vasculature. This review discusses our investigations into the role of cell-matrix mechanics in the mechanical regulation of angiogenesis. The experimental aspects of the research are based on in vitro experiments using an organ culture model of sprouting angiogenesis with the goal of developing new treatments and techniques to either promote or inhibit angiogenic outgrowth, depending on the application. Computational simulations were performed to simulate angiogenic growth coupled to matrix deformation, and live two-photon microscopy was used to obtain insight into the dynamic mechanical interaction between angiogenic neovessels and the extracellular matrix. In these studies, we characterized how angiogenic neovessels remodel the extracellular matrix (ECM) and how properties of the matrix such as density and boundary conditions influence vascular growth and alignment. Angiogenic neovessels extensively deform and remodel the matrix through a combination of applied traction, proteolytic activity, and generation of new cell-matrix adhesions. The angiogenic phenotype within endothelial cells is promoted by ECM deformation and remodeling. Sensitivity analysis using our finite element model of angiogenesis suggests that cell-generated traction during growth is the most important parameter controlling the deformation of the matrix and, therefore, angiogenic growth and remodeling. Live two-photon imaging has also revealed numerous neovessel behaviors during angiogenesis that are poorly understood such as episodic growth/regression, neovessel colocation, and anastomosis. Our research demonstrates that the topology of a resulting vascular network can be manipulated directly by modifying the mechanical interaction between angiogenic neovessels and the matrix.

A computational model of in vitro angiogenesis based on extracellular matrix fibre orientation.

Recent interest in the process of vascularisation within the biomedical community has motivated numerous new research efforts focusing on the process of angiogenesis. Although the role of chemical factors during angiogenesis has been well documented, the role of mechanical factors, such as the interaction between angiogenic vessels and the extracellular matrix, remains poorly understood. In vitro methods for studying angiogenesis exist; however, measurements available using such techniques often suffer from limited spatial and temporal resolutions. For this reason, computational models have been extensively employed to investigate various aspects of angiogenesis. This paper outlines the formulation and validation of a simple and robust computational model developed to accurately simulate angiogenesis based on length, branching and orientation morphometrics collected from vascularised tissue constructs. Microvessels were represented as a series of connected line segments. The morphology of the vessels was determined by a linear combination of the collagen fibre orientation, the vessel density gradient and a random walk component. Excellent agreement was observed between computational and experimental morphometric data over time. Computational predictions of microvessel orientation within an anisotropic matrix correlated well with experimental data. The accuracy of this modelling approach makes it a valuable platform for investigating the role of mechanical interactions during angiogenesis.

Determinants of microvascular network topologies in implanted neovasculatures.

OBJECTIVE: During neovascularization, the end result is a new functional microcirculation composed of a network of mature microvessels with specific topologies. Although much is known concerning the mechanisms underlying the initiation of angiogenesis, it remains unclear how the final architecture of microcirculatory beds is regulated. To begin to address this, we determined the impact of angiogenic neovessel prepatterning on the final microvascular network topology using a model of implant neovascularization. METHODS AND RESULTS: We used 3D direct-write bioprinting or physical constraints in a manner permitting postangiogenesis vascular remodeling and adaptation to pattern angiogenic microvascular precursors (neovessels formed from isolated microvessel segments) in 3D collagen gels before implantation and subsequent network formation. Neovasculatures prepatterned into parallel arrays formed functional networks after 4 weeks postimplantation but lost the prepatterned architecture. However, maintenance of uniaxial physical constraints during postangiogenesis remodeling of the implanted neovasculatures produced networks with aligned microvessels, as well as an altered proportional distribution of arterioles, capillaries, and venules. CONCLUSIONS: Here we show that network topology resulting from implanted microvessel precursors is independent from prepatterning of precursors but can be influenced by a patterning stimulus involving tissue deformation during postangiogenesis remodeling and maturation.

An in vitro assay of collagen fiber alignment by acupuncture needle rotation.

BACKGROUND: During traditional acupuncture therapy, soft tissues attach to and wind around the acupuncture needle. To study this phenomenon in a controlled and quantitative setting, we performed acupuncture needling in vitro. METHODS: Acupuncture was simulated in vitro in three-dimensional, type I collagen gels prepared at 1.5 mg/ml, 2.0 mg/ml, and 2.5 mg/ml collagen, and either crosslinked with formalin or left untreated. Acupuncture needles were inserted into the gels and rotated via a computer-controlled motor at 0.3 rev/sec for up to 10 revolutions while capturing the evolution of birefringence under cross-polarization. RESULTS: Simulated acupuncture produced circumferential alignment of collagen fibers close to the needle that evolved into radial alignment as the distance from the needle increased, which generally matched observations from published tissue explant studies. All gels failed prior to 10 revolutions, and the location of failure was near the transition between circumferential and radial alignment. Crosslinked collagen failed at a significantly lower number of revolutions than untreated collagen, whereas collagen concentration had no effect on gel failure. The strength of the alignment field increased with increasing collagen concentration and decreased with crosslinking. Separate studies were performed in which the gel thickness and depth of needle insertion were varied. As gel thickness increased, gels failed at fewer needle revolutions. For the same depth of insertion, alignment was greater in thinner gels. Alignment increased as the depth of insertion increased. CONCLUSION: These results indicate that the mechanostructural properties of soft connective tissues may affect their response to acupuncture therapy. The in vitro model provides a platform to study mechanotransduction during acupuncture in a highly controlled and quantitative setting.