Biomechanical characterization of murine pulmonary arteries

Remodeling of Murine Branch Pulmonary Arteries Under Chronic Hypoxia and Short-Term Normoxic Recovery

Chronic hypoxia plays a central role in diverse pulmonary pathologies, but its effects on longitudinal changes in the biomechanical behavior of proximal pulmonary arteries remain poorly understood. Similarly, effects of normoxic recovery have not been well studied. Here, we report hypoxia-induced changes in composition, vasoactivity, and passive biaxial mechanics in the main branch pulmonary artery of male C57BL/6J mice exposed to 10% FiO2 for 1, 2, or 3 weeks. We observed significant changes in extracellular matrix, and consequently wall mechanics, as early as 1 week of hypoxia. While circumferential stress and stiffness returned toward normal values by 2-3 weeks of hypoxia, area fractions of cytoplasm and thin collagen fibers did not return toward normal until after 1 week of normoxic recovery. By contrast, elastic energy storage and overall distensibility remained reduced after 3 weeks of hypoxia as well as following 1 week of normoxic recovery. While smooth muscle and endothelial cell responses were attenuated under hypoxia, smooth muscle but not endothelial cell responses recovered following 1 week of subsequent normoxia. Collectively, these data suggest that homeostatic processes were unable to preserve or restore overall function, at least over a brief period of normoxic recovery. Longitudinal changes are critical in understanding large pulmonary artery remodeling under hypoxia, and its reversal, and will inform predictive models of vascular adaptation.

Guidelines for mechanistic modeling and analysis in cardiovascular research

Computational, or in silico, models are an effective, noninvasive tool for investigating cardiovascular function. These models can be used in the analysis of experimental and clinical data to identify possible mechanisms of (ab)normal cardiovascular physiology. Recent advances in computing power and data management have led to innovative and complex modeling frameworks that simulate cardiovascular function across multiple scales. While commonly used in multiple disciplines, there is a lack of concise guidelines for the implementation of computer models in cardiovascular research. In line with recent calls for more reproducible research, it is imperative that scientists adhere to credible practices when developing and applying computational models to their research. The goal of this manuscript is to provide a consensus document that identifies best practices for in silico computational modeling in cardiovascular research. These guidelines provide the necessary methods for mechanistic model development, model analysis, and formal model calibration using fundamentals from statistics. We outline rigorous practices for computational, mechanistic modeling in cardiovascular research and discuss its synergistic value to experimental and clinical data.

Hemodynamics and Wall Mechanics of Vascular Graft Failure

Blood vessels are subjected to complex biomechanical loads, primarily from pressure-driven blood flow. Abnormal loading associated with vascular grafts, arising from altered hemodynamics or wall mechanics, can cause acute and progressive vascular failure and end-organ dysfunction. Perturbations to mechanobiological stimuli experienced by vascular cells contribute to remodeling of the vascular wall via activation of mechanosensitive signaling pathways and subsequent changes in gene expression and associated turnover of cells and extracellular matrix. In this review, we outline experimental and computational tools used to quantify metrics of biomechanical loading in vascular grafts and highlight those that show potential in predicting graft failure for diverse disease contexts. We include metrics derived from both fluid and solid mechanics that drive feedback loops between mechanobiological processes and changes in the biomechanical state that govern the natural history of vascular grafts. As illustrative examples, we consider application-specific coronary artery bypass grafts, peripheral vascular grafts, and tissue-engineered vascular grafts for congenital heart surgery as each of these involves unique circulatory environments, loading magnitudes, and graft materials.

A Systematic Comparison of Normal Structure and Function of the Greater Thoracic Vessels

The greater thoracic vessels are central to a well-functioning circulatory system and are often targeted in congenital heart surgeries, yet the structure and function of these vessels have not been well studied. Here we use consistent methods to quantify and compare microstructural features and biaxial biomechanical properties of the following six greater thoracic vessels in wild-type mice: ascending thoracic aorta, descending thoracic aorta, right subclavian artery, right pulmonary artery, thoracic inferior vena cava, and superior vena cava. Specifically, we determine volume fractions and orientations of the structurally significant wall constituents (i.e., collagen, elastin, and cell nuclei) using multiphoton imaging, and we quantify vasoactive responses and mechanobiologically relevant mechanical quantities (e.g., stress, stiffness) using computer-controlled biaxial mechanical testing. Similarities and differences across systemic, pulmonary, and venous circulations highlight underlying design principles of the vascular system. Results from this study represent another step towards understanding growth and remodeling of greater thoracic vessels in health, disease, and surgical interventions by providing baseline information essential for developing and validating predictive computational models.

Application and reduction of a nonlinear hyperelastic wall model capturing ex vivo relationships between fluid pressure, area, and wall thickness in normal and hypertensive murine left pulmonary arteries

Pulmonary hypertension is a cardiovascular disorder manifested by elevated mean arterial blood pressure (>20 mmHg) together with vessel wall stiffening and thickening due to alterations in collagen, elastin, and smooth muscle cells. Hypoxia-induced (type 3) pulmonary hypertension can be studied in animals exposed to a low oxygen environment for prolonged time periods leading to biomechanical alterations in vessel wall structure. This study introduces a novel approach to formulating a reduced order nonlinear elastic structural wall model for a large pulmonary artery. The model relating blood pressure and area is calibrated using ex vivo measurements of vessel diameter and wall thickness changes, under controlled pressure conditions, in left pulmonary arteries isolated from control and hypertensive mice. A two-layer, hyperelastic, and anisotropic model incorporating residual stresses is formulated using the Holzapfel-Gasser-Ogden model. Complex relations predicting vessel area and wall thickness with increasing blood pressure are derived and calibrated using the data. Sensitivity analysis, parameter estimation, subset selection, and physical plausibility arguments are used to systematically reduce the 16-parameter model to one in which a much smaller subset of identifiable parameters is estimated via solution of an inverse problem. Our final reduced one layer model includes a single set of three elastic moduli. Estimated ranges of these parameters demonstrate that nonlinear stiffening is dominated by elastin in the control animals and by collagen in the hypertensive animals. The pressure-area relation developed in this novel manner has potential impact on one-dimensional fluids network models of vessel wall remodeling in the presence of cardiovascular disease.

Optimization of resting tension for wire myography in male rat pulmonary arteries

Wire myography to test vasomotor functions of blood vessels ex-vivo are well-established for the systemic circulation, however, there is no consensus on protocols for pulmonary arteries. We created a standardized wire myography protocol for healthy rat PAs and validated this in a pulmonary hypertension (PH) model. Vessels stretched to higher initial tensions (5.0, 7.5 and 10.0 mN) exhibited a uniform response to phenylephrine, a larger dynamic range, and lower EC50 values. The endothelium-mediated relaxation showed that moderate tensions (7.5 and 10.0 mN) produced robust responses with higher maximum relaxation and lower EC50 values. For endothelium independent responses, the higher initial tension groups had lower and more consistent EC50 values than the lower initial tension groups. Pulmonary arteries from rats with PH were more responsive to vasoactive drugs when subjected to a higher initial tension. Notably, vessels in the PH group subjected to 15.0 mN exhibited high dynamic ranges in contractile and relaxation responses without tearing. Lastly, we observed attenuated cholinergic responses in these vessels-consistent with endothelial dysfunction in PH. Therefore, a moderate initial tension of 7.5-10.0 mN is optimal for healthy rat pulmonary arteries and a higher initial tension of 15.0 mN is optimal for pulmonary arteries from animals with PH.

A computational growth and remodeling framework for adaptive and maladaptive pulmonary arterial hemodynamics

Hemodynamic loading is known to contribute to the development and progression of pulmonary arterial hypertension (PAH). This loading drives changes in mechanobiological stimuli that affect cellular phenotypes and lead to pulmonary vascular remodeling. Computational models have been used to simulate mechanobiological metrics of interest, such as wall shear stress, at single time points for PAH patients. However, there is a need for new approaches that simulate disease evolution to allow for prediction of long-term outcomes. In this work, we develop a framework that models the pulmonary arterial tree through adaptive and maladaptive responses to mechanical and biological perturbations. We coupled a constrained mixture theory-based growth and remodeling framework for the vessel wall with a morphometric tree representation of the pulmonary arterial vasculature. We show that non-uniform mechanical behavior is important to establish the homeostatic state of the pulmonary arterial tree, and that hemodynamic feedback is essential for simulating disease time courses. We also employed a series of maladaptive constitutive models, such as smooth muscle hyperproliferation and stiffening, to identify critical contributors to development of PAH phenotypes. Together, these simulations demonstrate an important step toward predicting changes in metrics of clinical interest for PAH patients and simulating potential treatment approaches.

Smooth muscle contribution to vaginal viscoelastic response

Smooth muscle cells contribute to the mechanical function of various soft tissues, however, their contribution to the viscoelastic response when subjected to multiaxial loading remains unknown. The vagina is a fibromuscular viscoelastic organ that is exposed to prolonged and increased pressures with daily activities and physiologic processes such as vaginal birth. The vagina changes in geometry over time under prolonged pressure, known as creep. Vaginal smooth muscle cells may contribute to creep. This may be critical for the function of vaginal and other soft tissues that experience fluctuations in their biomechanical environment. Therefore, the objective of this study was to develop methods to evaluate the contribution of smooth muscle to vaginal creep under multiaxial loading using extension - inflation tests. The vaginas from wildtype mice (C57BL/6 × 129SvEv; 3-6 months; n = 10) were stimulated with various concentrations of potassium chloride then subjected to the measured in vivo pressure (7 mmHg) for 100 s. In a different cohort of mice (n = 5), the vagina was stimulated with a single concentration of potassium chloride then subjected to 5 and 15 mmHg. A laser micrometer measured vaginal outer diameter in real-time. Immunofluorescence evaluated the expression of alpha-smooth muscle actin and myosin heavy chain in the vaginal muscularis (n = 6). When smooth muscle contraction was activated, vaginal creep behavior increased compared to the relaxed state. However, increased pressure decreased the active creep response. This study demonstrated that extension - inflation protocols can be used to evaluate smooth muscle contribution to the viscoelastic response of tubular soft tissues.

New Computational Approach to Shunt Design in Congenital Heart Palliation

Shunts are commonly used to redirect blood to pulmonary arteries in procedures that palliate congenital cardiovascular defects. Previous clinical studies and hemodynamic simulations reveal a critical role of shunt diameter in balancing flow to pulmonary versus systemic vessels, but the biomechanical process of creating the requisite anastomosis between the shunt and host vessel has received little attention. Here, we report a new Lagrange multiplier-based finite element approach that represents the shunt and host vessels as individual structures and predicts the anastomosis geometry and attachment force that result when the shunt is sutured at an incision in the host, followed by pressurization. Simulations suggest that anastomosis orifice opening increases markedly with increasing length of the host incision and moderately with increasing blood pressure. The host artery is further predicted to conform to common stiff synthetic shunts, whereas more compliant umbilical vessel shunts should conform to the host, with orifice area transitioning between these two extremes via a Hill-type function of shunt stiffness. Moreover, a direct relationship is expected between attachment forces and shunt stiffness. This new computational approach promises to aid in surgical planning for diverse vascular shunts by predicting in vivo pressurized geometries.

Ex vivo biomechanical characterization of umbilical vessels: Possible shunts in congenital heart palliation

Children born with congenital heart defects typically undergo staged palliative surgeries to reconstruct the circulation to improve transport of deoxygenated blood to the lungs. As part of the first surgery, a temporary shunt (Blalock-Thomas-Taussig) is often created in neonates to connect a systemic and a pulmonary artery. Standard-of-care shunts are synthetic, which can lead to thrombosis, and much stiffer than the two host vessels, which can cause adverse mechanobiological responses. Moreover, the neonatal vasculature can undergo significant changes in size and structure over a short period, thus constraining the use of a non-growing synthetic shunt. Recent studies suggest that autologous umbilical vessels could serve as improved shunts, but there has not been a detailed biomechanical characterization of the four primary vessels - subclavian artery, pulmonary artery, umbilical vein, and umbilical artery. Herein, we biomechanically phenotype umbilical veins and arteries from prenatal mice (E18.5) and compare them to subclavian and pulmonary arteries harvested at two critical postnatal developmental ages (P10, P21). Comparisons include age-specific physiological conditions and simulated 'surgical-like' shunt conditions. Results suggest that the intact umbilical vein is a better choice as a shunt than the umbilical artery due to concerns with lumen closure and constriction related intramural damage in the latter. Yet, decellularization of umbilical arteries may be a viable alternative, with the possibility of host cellular infiltration and subsequent remodeling. Given recent efforts using autologous umbilical vessels as Blalock-Thomas-Taussig shunts in a clinical trial, our findings highlight aspects of the associated biomechanics that deserve further investigation.

A fluid-structure interaction model accounting arterial vessels as a key part of the blood-flow engine for the analysis of cardiovascular diseases

According to the classical Windkessel model, the heart is the only power source for blood flow, while the arterial system is assumed to be an elastic chamber that acts as a channel and buffer for blood circulation. In this paper we show that in addition to the power provided by the heart for blood circulation, strain energy stored in deformed arterial vessels in vivo can be transformed into mechanical work to propel blood flow. A quantitative relationship between the strain energy increment and functional (systolic, diastolic, mean and pulse blood pressure) and structural (stiffness, diameter and wall thickness) parameters of the aorta is described. In addition, details of blood flow across the aorta remain unclear due to changes in functional and other physiological parameters. Based on the arterial strain energy and fluid-structure interaction theory, the relationship between physiological parameters and blood supply to organs was studied, and a corresponding mathematical model was developed. The findings provided a new understanding about blood-flow circulation, that is, cardiac output allows blood to enter the aorta at an initial rate, and then strain energy stored in the elastic arteries pushes blood toward distal organs and tissues. Organ blood supply is a key factor in cardio-cerebrovascular diseases (CCVD), which are caused by changes in blood supply in combination with multiple physiological parameters. Also, some physiological parameters are affected by changes in blood supply, and vice versa. The model can explain the pathophysiological mechanisms of chronic diseases such as CCVD and hypertension among others, and the results are in good agreement with epidemiological studies of CCVD.

Compromised Cardiopulmonary Function in Fibulin-5 Deficient Mice

Competent elastic fibers are critical to the function of the lung and right circulation. Murine models of elastopathies can aid in understanding the functional roles of the elastin and elastin-associated glycoproteins that constitute elastic fibers. Here, we quantify together lung and pulmonary arterial structure, function, and mechanics with right heart function in a mouse model deficient in the elastin-associated glycoprotein fibulin-5. Differences emerged as a function of genotype, sex, and arterial region. Specifically, functional studies revealed increased lung compliance in fibulin-5 deficiency consistent with a histologically observed increased alveolar disruption. Biaxial mechanical tests revealed that the primary branch pulmonary arteries exhibit decreased elastic energy storage capacity and wall stress despite only modest differences in circumferential and axial material stiffness in the fibulin-5 deficient mice. Histological quantifications confirm a lower elastic fiber content in the fibulin-5 deficient pulmonary arteries, with fragmented elastic laminae in the outer part of the wall - likely the reason for reduced energy storage. Ultrasound measurements confirm sex differences in compromised right ventricular function in the fibulin-5 deficient mice. These results reveal compromised right heart function, but opposite effects of elastic fiber dysfunction on the lung parenchyma (significantly increased compliance) and pulmonary arteries (trend toward decreased distensibility), and call for further probing of ventilation-perfusion relationships in pulmonary pathologies. Amongst many other models, fibulin-5 deficient mice can contribute to our understanding of the complex roles of elastin in pulmonary health and disease.

Nicotine Affects Murine Aortic Stiffness and Fatigue Response During Supraphysiological Cycling

Nicotine exposure is a major risk factor for several cardiovascular diseases. Although the deleterious effects of nicotine on aortic remodeling processes have been studied to some extent, the biophysical consequences are not fully elucidated. In this investigation, we applied quasi-static and dynamic loading to quantify ways in which exposure to nicotine affects the mechanical behavior of murine arterial tissue. Segments of thoracic aortas from C57BL/6 mice exposed to 25 mg/kg/day of subcutaneous nicotine for 28 days were subjected to uniaxial tensile loading in an open-circumferential configuration. Comparing aorta segments from nicotine-treated mice relative to an equal number of control counterparts, stiffness in the circumferential direction was nearly twofold higher (377 kPa ± 165 kPa versus 191 kPa ± 65 kPa, n = 5, p = 0.03) at 50% strain. Using a degradative power-law fit to fatigue data at supraphysiological loading, we observed that nicotine-treated aortas exhibited significantly higher peak stress, greater loss of tension, and wider oscillation band than control aortas (p ≤ 0.01 for all three variables). Compared to simple stress relaxation tests, fatigue cycling is shown to be more sensitive and versatile in discerning nicotine-induced changes in mechanical behavior over many cycles. Supraphysiological fatigue cycling thus may have broader potential to reveal subtle changes in vascular mechanics caused by other exogenous toxins or pathological conditions.

Fully-Coupled FSI Computational Analyses in the Ascending Thoracic Aorta Using Patient-Specific Conditions and Anisotropic Material Properties

Computational hemodynamics has become increasingly important within the context of precision medicine, providing major insight in cardiovascular pathologies. However, finding appropriate compromise between speed and accuracy remains challenging in computational hemodynamics for an extensive use in decision making. For example, in the ascending thoracic aorta, interactions between the blood and the aortic wall must be taken into account for the sake of accuracy, but these fluid structure interactions (FSI) induce significant computational costs, especially when the tissue exhibits a hyperelastic and anisotropic response. The objective of the current study is to use the Small On Large (SOL) theory to linearize the anisotropic hyperelastic behavior in order to propose a reduced-order model for FSI simulations of the aorta. The SOL method is tested for fully-coupled FSI simulations in a patient-specific aortic geometry presenting an Ascending Thoracic Aortic Aneurysm (aTAA). The same model is also simulated with a fully-coupled FSI with non-linear material behavior, without SOL linearization. Eventually, the results and computational times with and without the SOL are compared. The SOL approach is demonstrated to provide a significant reduction of computational costs for FSI analysis in the aTAA, and the results in terms of stress state distribution are comparable. The method is implemented in ANSYS and will be further evaluated for clinical applications.

Role of fibulin-5 insufficiency and prolapse progression on murine vaginal biomechanical function

The vagina plays a critical role in supporting the pelvic organs and loss of support leads to pelvic organ prolapse. It is unknown what microstructural changes influence prolapse progression nor how decreased elastic fibers contributes to vaginal remodeling and smooth muscle contractility. The objective for this study was to evaluate the effect of fibulin-5 haploinsufficiency, and deficiency with progressive prolapse on the biaxial contractile and biomechanical function of the murine vagina. Vaginas from wildtype (n = 13), haploinsufficient (n = 13), and deficient mice with grade 1 (n = 9) and grade 2 or 3 (n = 9) prolapse were explanted for biaxial contractile and biomechanical testing. Multiaxial histology (n = 3/group) evaluated elastic and collagen fiber microstructure. Western blotting quantified protein expression (n = 6/group). A one-way ANOVA or Kruskal-Wallis test evaluated statistical significance. Pearson's or Spearman's test determined correlations with prolapse grade. Axial contractility decreased with fibulin-5 deficiency and POP (p < 0.001), negatively correlated with prolapse grade (ρ = - 0.80; p < 0.001), and positively correlated with muscularis elastin area fraction (ρ = - 0.78; p = 0.004). Circumferential (ρ = 0.71; p < 0.001) and axial (ρ = 0.69; p < 0.001) vaginal wall stresses positively correlated with prolapse grade. These findings demonstrated that fibulin-5 deficiency and prolapse progression decreased vaginal contractility and increased vaginal wall stress. Future work is needed to better understand the processes that contribute to prolapse progression in order to guide diagnostic, preventative, and treatment strategies.

Mechanisms of Hypoxia-Induced Pulmonary Arterial Stiffening in Mice Revealed by a Functional Genetics Assay of Structural, Functional, and Transcriptomic Data

Hypoxia adversely affects the pulmonary circulation of mammals, including vasoconstriction leading to elevated pulmonary arterial pressures. The clinical importance of changes in the structure and function of the large, elastic pulmonary arteries is gaining increased attention, particularly regarding impact in multiple chronic cardiopulmonary conditions. We establish a multi-disciplinary workflow to understand better transcriptional, microstructural, and functional changes of the pulmonary artery in response to sustained hypoxia and how these changes inter-relate. We exposed adult male C57BL/6J mice to normoxic or hypoxic (FiO₂ 10%) conditions. Excised pulmonary arteries were profiled transcriptionally using single cell RNA sequencing, imaged with multiphoton microscopy to determine microstructural features under in vivo relevant multiaxial loading, and phenotyped biomechanically to quantify associated changes in material stiffness and vasoactive capacity. Pulmonary arteries of hypoxic mice exhibited an increased material stiffness that was likely due to collagen remodeling rather than excessive deposition (fibrosis), a change in smooth muscle cell phenotype reflected by decreased contractility and altered orientation aligning these cells in the same direction as the remodeled collagen fibers, endothelial proliferation likely representing endothelial-to-mesenchymal transitioning, and a network of cell-type specific transcriptomic changes that drove these changes. These many changes resulted in a system-level increase in pulmonary arterial pulse wave velocity, which may drive a positive feedback loop exacerbating all changes. These findings demonstrate the power of a multi-scale genetic-functional assay. They also highlight the need for systems-level analyses to determine which of the many changes are clinically significant and may be potential therapeutic targets.

Computational simulation-derived hemodynamic and biomechanical properties of the pulmonary arterial tree early in the course of ventricular septal defects

Untreated ventricular septal defects (VSDs) can lead to pulmonary arterial hypertension (PAH) characterized by elevated pulmonary artery (PA) pressure and vascular remodeling, known as PAH associated with congenital heart disease (PAH-CHD). Though previous studies have investigated hemodynamic effects on vascular mechanobiology in late-stage PAH, hemodynamics leading to PAH-CHD initiation have not been fully quantified. We hypothesize that abnormal hemodynamics from left-to-right shunting in early stage VSDs affects PA biomechanical properties leading to PAH initiation. To model PA hemodynamics in healthy, small, moderate, and large VSD conditions prior to the onset of vascular remodeling, computational fluid dynamics simulations were performed using a 3D finite element model of a healthy 1-year-old's proximal PAs and a body-surface-area-scaled 0D distal PA tree. VSD conditions were modeled with increased pulmonary blood flow to represent degrees of left-to-right shunting. In the proximal PAs, pressure, flow, strain, and wall shear stress (WSS) increased with increasing VSD size; oscillatory shear index decreased with increasing VSD size in the larger PA vessels. WSS was higher in smaller diameter vessels and increased with VSD size, with the large VSD condition exhibiting WSS >100 dyn/cm[Formula: see text], well above values typically used to study dysfunctional mechanotransduction pathways in PAH. This study is the first to estimate hemodynamic and biomechanical metrics in the entire pediatric PA tree with VSD severity at the stage leading to PAH initiation and has implications for future studies assessing effects of abnormal mechanical stimuli on endothelial cells and vascular wall mechanics that occur during PAH-CHD initiation and progression.

Arterial Wall Stiffening in Caveolin-1 Deficiency-Induced Pulmonary Artery Hypertension in Mice

BACKGROUND

Pulmonary artery hypertension (PAH) is a complex disorder that can lead to right heart failure. The generation of caveolin-1 deficient mice (CAV-1-/-) has provided an alternative genetic model to study the mechanisms of pulmonary hypertension. However, the vascular adaptations in these mice have not been characterized.

OBJECTIVE

To determine the histological and functional changes in the pulmonary and carotid arteries in CAV-1-/- induced PAH.

METHODS

Pulmonary and carotid arteries of young (4-6 months old) and mature (9-12 months old) CAV-1-/- mice were tested and compared to normal wild type mice.

RESULTS

Artery stiffness increases in CAV-1-/- mice, especially the circumferential stiffness of the pulmonary arteries. Increases in stiffness were quantified by a decrease in circumferential stretch and transition strain, increases in elastic moduli, and an increase in total strain energy at physiologic strains. Changes in mechanical properties for the pulmonary artery correlated with increased collagen content while carotid artery mechanical properties correlated with decreased elastin content.

CONCLUSIONS

We demonstrated that an increase in artery stiffness is associated with CAV-1 deficiency-induced pulmonary hypertension. These results improve our understanding of artery remodeling in PAH.

Vascular adaptation in the presence of external support - A modeling study

Vascular grafts have long been used to replace damaged or diseased vessels with considerable success, but a new approach is emerging where native vessels are merely supported, not replaced. Although external supports have been evaluated in diverse situations - ranging from aneurysmal disease to vein grafts or the Ross operation - optimal supports and procedures remain wanting. In this paper, we present a novel application of a growth and remodeling model well suited for parametrically exploring multiple designs of external supports while accounting for mechanobiological and immunobiological responses of the supported native vessel. These results suggest that a load bearing external support can reduce vessel thickening in response to pressure elevation. Results also suggest that the final adaptive state of the vessel depends on the structural stiffness of the support via a mechano-driven adaptation, although luminal encroachment may be a complication in the presence of chronic inflammation. Finally, the supported vessel can stiffen (structurally and materially) along circumferential and axial directions, which could have implications on overall hemodynamics and thus subsequent vascular remodeling. The proposed framework can provide valuable insights into vascular adaptation in the presence of external support, accelerate rational design, and aid translation of this emerging approach.

Robust latent-variable interpretation of in vivo regression models by nested resampling

Simple multilinear methods, such as partial least squares regression (PLSR), are effective at interrelating dynamic, multivariate datasets of cell-molecular biology through high-dimensional arrays. However, data collected in vivo are more difficult, because animal-to-animal variability is often high, and each time-point measured is usually a terminal endpoint for that animal. Observations are further complicated by the nesting of cells within tissues or tissue sections, which themselves are nested within animals. Here, we introduce principled resampling strategies that preserve the tissue-animal hierarchy of individual replicates and compute the uncertainty of multidimensional decompositions applied to global averages. Using molecular-phenotypic data from the mouse aorta and colon, we find that interpretation of decomposed latent variables (LVs) changes when PLSR models are resampled. Lagging LVs, which statistically improve global-average models, are unstable in resampled iterations that preserve nesting relationships, arguing that these LVs should not be mined for biological insight. Interestingly, resampling is less discriminatory for multidimensional regressions of in vitro data, where replicate-to-replicate variance is sufficiently low. Our work illustrates the challenges and opportunities in translating systems-biology approaches from cultured cells to living organisms. Nested resampling adds a straightforward quality-control step for interpreting the robustness of in vivo regression models.