Spatial Configurations of 3D Extracellular Matrix Collagen Density and Anisotropy Simultaneously Guide Angiogenesis.

Extracellular matrix (ECM) collagen density and fibril anisotropy are thought to affect the development of new vasculatures during pathologic and homeostatic angiogenesis. Computational simulation is emerging as a tool to investigate the role of matrix structural configurations on cell guidance. However, prior computational models have only considered the orientation of collagen as a model input. Recent experimental evidence indicates that cell guidance is simultaneously influenced by the direction and intensity of alignment (i.e., degree of anisotropy) as well as the local collagen density. The objective of this study was to explore the role of ECM collagen anisotropy and density during sprouting angiogenesis through simulation in the AngioFE and FEBio modeling frameworks. AngioFE is a plugin for FEBio (Finite Elements for Biomechanics) that simulates cell-matrix interactions during sprouting angiogenesis. We extended AngioFE to represent ECM collagen as deformable 3D ellipsoidal fibril distributions (EFDs). The rate and direction of microvessel growth were modified to depend simultaneously on the ECM collagen anisotropy (orientation and degree of anisotropy) and density. The sensitivity of growing neovessels to these stimuli was adjusted so that AngioFE could reproduce the growth and guidance observed in experiments where microvessels were cultured in collagen gels of varying anisotropy and density. We then compared outcomes from simulations using EFDs to simulations that used AngioFE's prior vector field representation of collagen anisotropy. We found that EFD simulations were more accurate than vector field simulations in predicting experimentally observed microvessel guidance. Predictive simulations demonstrated the ability of anisotropy gradients to recruit microvessels across short and long distances relevant to wound healing. Further, simulations predicted that collagen alignment could enable microvessels to overcome dense tissue interfaces such as tumor-associated collagen structures (TACS) found in desmoplasia and tumor-stroma interfaces. This approach can be generalized to other mechanobiological relationships during cell guidance phenomena in computational settings.

Dynamic Biophysical Cues Near the Tip Cell Microenvironment Provide Distinct Guidance Signals to Angiogenic Neovessels.

The formation of new vascular networks via angiogenesis is a crucial biological mechanism to balance tissue metabolic needs, yet the coordination of factors that influence the guidance of growing neovessels remain unclear. This study investigated the influence of extracellular cues within the immediate environment of sprouting tips over multiple hours and obtained quantitative relationships describing their effects on the growth trajectories of angiogenic neovessels. Three distinct microenvironmental cues-fibril tracks, ECM density, and the presence of nearby cell bodies-were extracted from 3D time series image data. The prominence of each cue was quantified along potential sprout trajectories to predict the response to multiple microenvironmental factors simultaneously. Sprout trajectories significantly correlated with the identified microenvironmental cues. Specifically, ECM density and nearby cellular bodies were the strongest predictors of the trajectories taken by neovessels (p < 0.001 and p = 0.016). Notwithstanding, direction changing trajectories, deviating from the initial neovessel orientation, were significantly correlated with fibril tracks (p = 0.003). Direction changes also occurred more frequently with strong microenvironmental cues. This provides evidence for the first time that local matrix fibril alignment influences changes in sprout trajectories but does not materially contribute to persistent sprouting. Together, our results suggest the microenvironmental cues significantly contribute to guidance of sprouting trajectories. Further, the presented methods quantitatively distinguish the influence of individual microenvironmental stimuli during guidance.

Matrix anisotropy promotes angiogenesis in a density-dependent manner.

Angiogenesis is necessary for wound healing, tumorigenesis, implant inosculation, and homeostasis. In each situation, matrix structure and mechanics play a role in determining whether new vasculatures can establish transport to new or hypoxic tissues. Neovessel growth and directional guidance are sensitive to three-dimensional (3-D) matrix anisotropy and density, although the individual and integrated roles of these matrix features have not been fully recapitulated in vitro. We developed a tension-based method to align 3-D collagen constructs seeded with microvessel fragments in matrices of three levels of collagen fibril anisotropy and two levels of collagen density. The extent and direction of neovessel growth from the parent microvessel fragments increased with matrix anisotropy and decreased with density. The proangiogenic effects of anisotropy were attenuated at higher matrix densities. We also examined the impact of matrix anisotropy in an experimental model of neovessel invasion across a tissue interface. Matrix density was found to dictate the success of interface crossing, whereas interface curvature and fibril alignment were found to control directional guidance. Our findings indicate that complex configurations of matrix density and alignment can facilitate or complicate the establishment or maintenance of vascular networks in pathological and homeostatic angiogenesis. Furthermore, we extend preexisting methods for tuning collagen anisotropy in thick constructs. This approach addresses gaps in tissue engineering and cell culture by supporting the inclusion of large multicellular structures in prealigned constructs.NEW & NOTEWORTHY Matrix anisotropy and density have a considerable effect on angiogenic vessel growth and directional guidance. However, the current literature relies on 2-D and simplified models of angiogenesis (e.g., tubulogenesis and vasculogenesis). We present a method to align 3-D collagen scaffolds embedded with microvessel fragments to different levels of anisotropy. Neovessel growth increases with anisotropy and decreases with density, which may guide angiogenic neovessels across tissue interfaces such as during implant inosculation and tumorigenesis.

A Nonparametric Approach for Estimating Three-Dimensional Fiber Orientation Distribution Functions (ODFs) in Fibrous Materials.

Many biological tissues contain an underlying fibrous microstructure that is optimized to suit a physiological function. The fiber architecture dictates physical characteristics such as stiffness, diffusivity, and electrical conduction. Abnormal deviations of fiber architecture are often associated with disease. Thus, it is useful to characterize fiber network organization from image data in order to better understand pathological mechanisms. We devised a method to quantify distributions of fiber orientations based on the Fourier transform and the Qball algorithm from diffusion MRI. The Fourier transform was used to decompose images into directional components, while the Qball algorithm efficiently converted the directional data from the frequency domain to the orientation domain. The representation in the orientation domain does not require any particular functional representation, and thus the method is nonparametric. The algorithm was verified to demonstrate its reliability and used on datasets from microscopy to show its applicability. This method increases the ability to extract information of microstructural fiber organization from experimental data that will enhance our understanding of structure-function relationships and enable accurate representation of material anisotropy in biological tissues.

Ventricular-vascular coupling is predictive of adverse clinical outcome in paediatric pulmonary arterial hypertension.

AIMS: Ventricular-vascular coupling, the ratio between the right ventricle's contractile state (Ees) and its afterload (Ea), may be a useful metric in the management of paediatric pulmonary arterial hypertension (PAH). In this study we assess the prognostic capacity of the ventricular-vascular coupling ratio (Ees/Ea) derived using right ventricular (RV) pressure alone in children with PAH. METHODS: One hundred and thirty paediatric patients who were diagnosed with PAH via right heart catheterisation were retrospectively reviewed over a 10-year period. Maximum RV isovolumic pressure and end-systolic pressure were estimated using two single-beat methods from Takeuchi et al (Ees/Ea_(Takeuchi)) and from Kind et al (Ees/Ea_(Kind)) and used with an estimate of end-systolic pressure to compute ventricular-vascular coupling from pressure alone. Patients were identified as either idiopathic/hereditary PAH or associated PAH (IPAH/HPAH and APAH, respectively). Haemodynamic data, clinical functional class and clinical worsening outcomes-separated into soft (mild) and hard (severe) event categories-were assessed. Adverse soft events included functional class worsening, syncopal event, hospitalisation due to a proportional hazard-related event and haemoptysis. Hard events included death, transplantation, initiation of prostanoid therapy and hospitalisation for atrial septostomy and Pott's shunt. Cox proportional hazard modelling was used to assess whether Ees/Ea was predictive of time-to-event. RESULTS: In patients with IPAH/HPAH, Ees/Ea_(Kind) and Ees/Ea_(Takeuchi) were both independently associated with time to hard event (p=0.003 and p=0.001, respectively) and when adjusted for indexed pulmonary vascular resistance (p=0.032 and p=0.013, respectively). Neither Ees/Ea_(Kind) nor Ees/Ea_(Takeuchi) were associated with time to soft event. In patients with APAH, neither Ees/Ea_(Kind) nor Ees/Ea_(Takeuchi) were associated with time to hard event or soft event. CONCLUSIONS: Ees/Ea derived from pressure alone is a strong independent predictor of adverse outcome and could be a potential powerful prognostic tool for paediatric PAH.

Stromal Cells Promote Neovascular Invasion Across Tissue Interfaces.

Vascular connectivity between adjacent vessel beds within and between tissue compartments is essential to any successful neovascularization process. To establish new connections, growing neovessels must locate other vascular elements during angiogenesis, often crossing matrix and other tissue-associated boundaries and interfaces. How growing neovessels traverse any tissue interface, whether part of the native tissue structure or secondary to a regenerative procedure (e.g., an implant), is not known. In this study, we developed an experimental model of angiogenesis wherein growing neovessels must interact with a 3D interstitial collagen matrix interface that separates two distinct tissue compartments. Using this model, we determined that matrix interfaces act as a barrier to neovessel growth, deflecting growing neovessels parallel to the interface. Computational modeling of the neovessel/matrix biomechanical interactions at the interface demonstrated that differences in collagen fibril density near and at the interface are the likely mechanism of deflection, while fibril alignment guides deflected neovessels along the interface. Interestingly, stromal cells facilitated neovessel interface crossing during angiogenesis via a vascular endothelial growth factor (VEGF)-A dependent process. However, ubiquitous addition of VEGF-A in the absence of stromal cells did not promote interface invasion. Therefore, our findings demonstrate that vascularization of a tissue via angiogenesis involves stromal cells providing positional cues to the growing neovasculature and provides insight into how a microvasculature is organized within a tissue.

Imaging the Dynamic Interaction Between Sprouting Microvessels and the Extracellular Matrix.

Thorough understanding of growth and evolution of tissue vasculature is fundamental to many fields of medicine including cancer therapy, wound healing, and tissue engineering. Angiogenesis, the growth of new vessels from existing ones, is dynamically influenced by a variety of environmental factors, including mechanical and biophysical factors, chemotactic factors, proteolysis, and interaction with stromal cells. Yet, dynamic interactions between neovessels and their environment are difficult to study with traditional fixed time imaging techniques. Advancements in imaging technologies permit time-series and volumetric imaging, affording the ability to visualize microvessel growth over 3D space and time. Time-lapse imaging has led to more informative investigations of angiogenesis. The environmental factors implicated in angiogenesis span a wide range of signals. Neovessels advance through stromal matrices by forming attachments and pulling and pushing on their microenvironment, reorganizing matrix fibers, and inducing large deformations of the surrounding stroma. Concurrently, neovessels secrete proteolytic enzymes to degrade their basement membrane, create space for new vessels to grow, and release matrix-bound cytokines. Growing neovessels also respond to a host of soluble and matrix-bound growth factors, and display preferential growth along a cytokine gradient. Lastly, stromal cells such as macrophages and mesenchymal stem cells (MSCs) interact directly with neovessels and their surrounding matrix to facilitate sprouting, vessel fusion, and tissue remodeling. This review highlights how time-lapse imaging techniques advanced our understanding of the interaction of blood vessels with their environment during sprouting angiogenesis. The technology provides means to characterize the evolution of microvessel behavior, providing new insights and holding great promise for further research on the process of angiogenesis.

Chronic kidney disease and aging differentially diminish bone material and microarchitecture in C57Bl/6 mice.

Chronic kidney disease (CKD) is a common disease of aging and increases fracture risk over advanced age alone. Aging and CKD differently impair bone turnover and mineralization. We thus hypothesize that the loss of bone quality would be greatest with the combination of advanced age and CKD. We evaluated bone from young adult (6 mo.), middle-age (18 mo.), and old (24 mo.) male C57Bl/6 mice three months following either 5/6th nephrectomy, to induce CKD, or Sham procedures. CKD exacerbated losses of cortical and trabecular microarchitecture associated with aging. Aging and CKD each resulted in thinner, more porous cortices and fewer and thinner trabeculae. Bone material quality was also reduced with CKD, and these changes to bone material were distinct from those due to age. Aging reduced whole-bone flexural strength and modulus, micrometer-scale nanoindentation modulus, and nanometer-scale tissue and collagen strain (small-angle x-ray scattering [SAXS]. By contrast, CKD reduced work to fracture and variation in bone tissue modulus and composition (Raman spectroscopy), and increased percent collagen strain. The increased collagen strain burden was associated with loss of toughness in CKD. In addition, osteocyte lacunae became smaller, sparser, and more disordered with age for Sham mice, yet these age-related changes were not clearly observed in CKD. However, for CKD, larger lacunae positively correlated with increased serum phosphate levels, suggesting that osteocytes play a role in systemic mineral homeostasis. This work demonstrates that CKD reduces bone quality, including microarchitecture and bone material properties, and that loss of bone quality with age is compounded by CKD. These findings may help reconcile why bone mass does not consistently predict fracture in the CKD population, as well as why older individuals with CKD are at high risk of fragility.

A new open-source tool for measuring 3D osteocyte lacunar geometries from confocal laser scanning microscopy reveals age-related changes to lacunar size and shape in cortical mouse bone.

Osteocytes can participate in systemic mineral homeostasis through perilacunar maintenance and remodeling, where changes to osteocyte lacunar morphology may affect bone structural integrity, tissue strains, and osteocyte mechanosensitivity. Though aging is associated with both decreased bone quality and altered mineral metabolism, it is not known if osteocyte lacunae undergo age-related changes in geometry. In order to survey lacunar changes with age, we developed an open-source program whereby 3D osteocyte lacunae are automatically segmented and then subsequently reconstructed from confocal laser scanning microscopy (CLSM) depth stacks for quantitative analysis of geometry and orientation. This approach takes advantage of the availability and speed of CLSM while avoiding time-consuming and bias-prone manual segmentation. Unlike conventional approaches used to quantify osteocyte lacunar morphology, CLSM enables facile analysis in three-dimensions with clear identification of osteocyte lacunae. We report that 3D osteocyte lacunae measured by CLSM become smaller, more spherical, more oblate, more spatially disorganized, and more sparsely populated with increased age in C57Bl/6 mouse cortical bone in groups spanning 6-24 months old. Critically, these age-related changes are in large part not observed in 2D analyses from the same samples. These results (1) demonstrate proof-of-concept of an efficient method to quantitatively assess osteocyte lacunae in 3D for application to a wide range of studies and (2) motivate further inquiry into how changes to osteocyte lacunar geometries and perilacunar material contribute to diminished bone quality in aging.