Haemodynamics in the mouse aortic arch computed from MRI-derived velocities at the aortic root

Morphological characterisation of pediatric Turner syndrome aortae: Insights from a small cohort study

Cardiovascular disease is widespread in girls and women living with Turner syndrome (TS). Despite this prevalence, cardiovascular risk evaluation using the current guidelines has seen life-threatening aortic events occurring at dimensions classified within the normal threshold. In this study, we characterized the three-dimensional aortic geometries of Turner syndrome children and their age-matched healthy counterparts to evaluate various morphological parameters. Turner syndrome girls had overall greater values in ten out of fifteen parameters examined (p > 0.05), when compared to healthy children: the aortic arch height and width; the ascending aorta, aortic arch (2 locations), and descending aorta diameters; the ratio of the ascending to descending aorta diameter; average curvature; average torsion; and average curvature-torsion score. Additionally, significant associations were found in the TS group: body surface area and both arch height (p = 0.03) and arch height to width ratio (p = 0.05), and aortic arch diameter and both body surface area (p = 0.04) and weight (p = 0.04). The new information resulting from this small cohort study contributes to an improved understanding of the morphological parameters affecting the hemodynamic environment in TS, and the clinical assessment of the increased cardiovascular risk in this population.

Reconstruction and Validation of Arterial Geometries for Computational Fluid Dynamics Using Multiple Temporal Frames of 4D Flow-MRI Magnitude Images

PURPOSE

Segmentation and reconstruction of arterial blood vessels is a fundamental step in the translation of computational fluid dynamics (CFD) to the clinical practice. Four-dimensional flow magnetic resonance imaging (4D Flow-MRI) can provide detailed information of blood flow but processing this information to elucidate the underlying anatomical structures is challenging. In this study, we present a novel approach to create high-contrast anatomical images from retrospective 4D Flow-MRI data.

METHODS

For healthy and clinical cases, the 3D instantaneous velocities at multiple cardiac time steps were superimposed directly onto the 4D Flow-MRI magnitude images and combined into a single composite frame. This new Composite Phase-Contrast Magnetic Resonance Angiogram (CPC-MRA) resulted in enhanced and uniform contrast within the lumen. These images were subsequently segmented and reconstructed to generate 3D arterial models for CFD. Using the time-dependent, 3D incompressible Reynolds-averaged Navier-Stokes equations, the transient aortic haemodynamics was computed within a rigid wall model of patient geometries.

RESULTS

Validation of these models against the gold standard CT-based approach showed no statistically significant inter-modality difference regarding vessel radius or curvature (p > 0.05), and a similar Dice Similarity Coefficient and Hausdorff Distance. CFD-derived near-wall hemodynamics indicated a significant inter-modality difference (p > 0.05), though these absolute errors were small. When compared to the in vivo data, CFD-derived velocities were qualitatively similar.

CONCLUSION

This proof-of-concept study demonstrated that functional 4D Flow-MRI information can be utilized to retrospectively generate anatomical information for CFD models in the absence of standard imaging datasets and intravenous contrast.

Calibration of patient-specific boundary conditions for coupled CFD models of the aorta derived from 4D Flow-MRI

Introduction: Patient-specific computational fluid dynamics (CFD) models permit analysis of complex intra-aortic hemodynamics in patients with aortic dissection (AD), where vessel morphology and disease severity are highly individualized. The simulated blood flow regime within these models is sensitive to the prescribed boundary conditions (BCs), so accurate BC selection is fundamental to achieve clinically relevant results. Methods: This study presents a novel reduced-order computational framework for the iterative flow-based calibration of 3-Element Windkessel Model (3EWM) parameters to generate patient-specific BCs. These parameters were calibrated using time-resolved flow information derived from retrospective four-dimensional flow magnetic resonance imaging (4D Flow-MRI). For a healthy and dissected case, blood flow was then investigated numerically in a fully coupled zero dimensional-three dimensional (0D-3D) numerical framework, where the vessel geometries were reconstructed from medical images. Calibration of the 3EWM parameters was automated and required ~3.5 min per branch. Results: With prescription of the calibrated BCs, the computed near-wall hemodynamics (time-averaged wall shear stress, oscillatory shear index) and perfusion distribution were consistent with clinical measurements and previous literature, yielding physiologically relevant results. BC calibration was particularly important in the AD case, where the complex flow regime was captured only after BC calibration. Discussion: This calibration methodology can therefore be applied in clinical cases where branch flow rates are known, for example, via 4D Flow-MRI or ultrasound, to generate patient-specific BCs for CFD models. It is then possible to elucidate, on a case-by-case basis, the highly individualized hemodynamics which occur due to geometric variations in aortic pathology high spatiotemporal resolution through CFD.

Blood flow modeling reveals improved collateral artery performance during the regenerative period in mammalian hearts

Collateral arteries bridge opposing artery branches, forming a natural bypass that can deliver blood flow downstream of an occlusion. Inducing coronary collateral arteries could treat cardiac ischemia, but more knowledge on their developmental mechanisms and functional capabilities is required. Here we used whole-organ imaging and three-dimensional computational fluid dynamics modeling to define spatial architecture and predict blood flow through collaterals in neonate and adult mouse hearts. Neonate collaterals were more numerous, larger in diameter and more effective at restoring blood flow. Decreased blood flow restoration in adults arose because during postnatal growth coronary arteries expanded by adding branches rather than increasing diameters, altering pressure distributions. In humans, adult hearts with total coronary occlusions averaged 2 large collaterals, with predicted moderate function, while normal fetal hearts showed over 40 collaterals, likely too small to be functionally relevant. Thus, we quantify the functional impact of collateral arteries during heart regeneration and repair-a critical step toward realizing their therapeutic potential.

Experimental and Mouse-Specific Computational Models of the Fbln4SMKO Mouse to Identify Potential Biomarkers for Ascending Thoracic Aortic Aneurysm

PURPOSE

To use computational methods to explore geometric, mechanical, and fluidic biomarkers that could correlate with mouse lifespan in the Fbln4SMKO mouse. Mouse lifespan was used as a surrogate for risk of a severe cardiovascular event in cases of ascending thoracic aortic aneurysm.

METHODS

Image-based, mouse-specific fluid-structure-interaction models were developed for Fbln4SMKO mice (n = 10) at ages two and six months. The results of the simulations were used to quantify potential biofluidic biomarkers, complementing the geometrical biomarkers obtained directly from the images.

RESULTS

Comparing the different geometrical and biofluidic biomarkers to the mouse lifespan, it was found that mean oscillatory shear index (OSImin) and minimum time-averaged wall shear stress (TAWSSmin) at six months showed the largest correlation with lifespan (r2 = 0.70, 0.56), with both correlations being positive (i.e., mice with high OSImean and high TAWSSmin tended to live longer). When change between two and six months was considered, the change in TAWSSmin showed a much stronger correlation than OSImean (r2 = 0.75 vs. 0.24), and the correlation was negative (i.e., mice with increasing TAWSSmin over this period tended to live less long).

CONCLUSION

The results highlight potential biomarkers of ATAA outcomes that can be obtained through noninvasive imaging and computational simulations, and they illustrate the potential synergy between small-animal and computational models.

Characterization of Flow Dynamics in the Pulmonary Bifurcation of Patients With Repaired Tetralogy of Fallot: A Computational Approach

The hemodynamic environment of the pulmonary bifurcation is of great importance for adult patients with repaired tetralogy of Fallot (rTOF) due to possible complications in the pulmonary valve and narrowing of the left pulmonary artery (LPA). The aim of this study was to computationally investigate the effect of geometrical variability and flow split on blood flow characteristics in the pulmonary trunk of patient-specific models. Data from a cohort of seven patients was used retrospectively and the pulmonary hemodynamics was investigated using averaged and MRI-derived patient-specific boundary conditions on the individualized models, as well as a statistical mean geometry. Geometrical analysis showed that curvature and tortuosity are higher in the LPA branch, compared to the right pulmonary artery (RPA), resulting in complex flow patterns in the LPA. The computational analysis also demonstrated high time-averaged wall shear stress (TAWSS) at the outer wall of the LPA and the wall of the RPA proximal to the junction. Similar TAWSS patterns were observed for averaged boundary conditions, except for a significantly modified flow split assigned at the outlets. Overall, this study enhances our understanding about the flow development in the pulmonary bifurcation of rTOF patients and associates some morphological characteristics with hemodynamic parameters, highlighting the importance of patient-specificity in the models. To confirm these findings, further studies are required with a bigger cohort of patients.

Effect of Subject-Specific, Spatially Reduced, and Idealized Boundary Conditions on the Predicted Hemodynamic Environment in the Murine Aorta

Mouse models of atherosclerosis have become effective resources to study atherogenesis, including the relationship between hemodynamics and lesion development. Computational methods aid the prediction of the in vivo hemodynamic environment in the mouse vasculature, but careful selection of inflow and outflow boundary conditions (BCs) is warranted to promote model accuracy. Herein, we investigated the impact of animal-specific versus reduced/idealized flow boundary conditions on predicted blood flow patterns in the mouse thoracic aorta. Blood velocities were measured in the aortic root, arch branch vessel, and descending aorta in ApoE^-/- mice using phase-contrast MRI. Computational geometries were derived from micro-CT imaging and combinations of high-fidelity or reduced/idealized MR-derived BCs were applied to predict the bulk flow field and hemodynamic metrics (e.g., wall shear stress, WSS; cross-flow index, CFI). Results demonstrate that pressure-free outlet BCs significantly overestimate outlet flow rates as compared to measured values. When compared to models that incorporate 3-component inlet velocity data [[Formula: see text]] and time-varying outlet mass flow splits [[Formula: see text]] (i.e., high-fidelity model), neglecting in-plane inlet velocity components (i.e., [Formula: see text])) leads to errors in WSS and CFI values ranging from 10 to 30% across the model domain whereas the application of a steady outlet mass flow splits results in negligible differences in these hemodynamics metrics. This investigation highlights that 3-component inlet velocity data and at least steady mass flow splits are required for accurate predictions of flow patterns in the mouse thoracic aorta.

The effect of coarctation degrees on wall shear stress indices

Coarctation of the aorta (CoA) is a congenital tightening of the proximal descending aorta. Flow quantification can be immensely valuable for an early and accurate diagnosis. However, there is a lack of appropriate diagnostic approaches for a variety of cardiovascular diseases, such as CoA. An accurate understanding of the disease depends on measurements of the global haemodynamics (criteria for heart function) and also the local haemodynamics (detailed data on the dynamics of blood flow). Playing a significant role in clinical processes, wall shear stress (WSS) cannot be measured clinically; thus, computation tools are needed to give an insight into this crucial haemodynamic parameter. In the present study, in order to enable the progress of non-invasive approaches that quantify global and local haemodynamics for different CoA severities, innovative computational blueprint simulations that include fluid-solid interaction models are developed. Since there is no clear approach for managing the CoA regarding its severity, this study proposes the use of WSS indices and pressure gradient to better establish a framework for treatment procedures in CoA patients with different severities. This provides a platform for improving CoA therapy on a patient-specific level, in which physicians can perform treatment methods based on WSS indices on top of using a mere experience. Results show how severe CoA affects the aorta in comparison to the milder cases, which can give the medical community valuable information before and after any intervention.

Numerical analysis of the hemodynamics of rat aorta based on magnetic resonance imaging and fluid-structure interaction

Murine models have been widely used to investigate the mechanobiology of aortic atherosclerosis and dissections, which develop preferably at different anatomic locations of aorta. Based MRI and finite element analysis with fluid-structure interaction, we numerically investigated factors that may affect the blood flow and structural mechanics of rat aorta. The results indicated that aortic root motion greatly increases time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), relative residence time (RRT), displacement of the aorta, and enhances helical flow pattern but has limited influence on effective stress, which is highly modulated by blood pressure. Moreover, the influence of the motion component on these indicators is different with axial motion more obvious than planar motion. Surrounding fixation of the intercostal arteries and the branch vessels on aortic arch would reduce the influence of aortic root motion. The compliance of the aorta has different influences at different regions, leading to decrease in TAWSS and helical flow, increase in OSI, RRT at the aortic arch, but has reversed effects on the branch vessels. When compared with the steady flow, the pulsatile blood flow would obviously increase the WSS, the displacement, and the effective stress in most regions. In conclusion, to accurately quantify the blood flow and structural mechanics of rat aorta, the motion of the aortic root, the compliance of aortic wall, and the pulsation of blood flow should be considered. However, when only focusing on the effective stress in rat aorta, the motion of the aortic root may be neglected.

Hemodynamic Abnormalities in the Aorta of Turner Syndrome Girls

Congenital abnormalities in girls and women with Turner syndrome (TS), alongside an underlying predisposition to obesity and hypertension, contribute to an increased risk of cardiovascular disease and ultimately reduced life expectancy. We observe that children with TS present a greater variance in aortic arch morphology than their healthy counterparts, and hypothesize that their hemodynamics is also different. In this study, computational fluid dynamic (CFD) simulations were performed for four TS girls, and three age-matched healthy girls, using patient-specific inlet boundary conditions, obtained from phase-contrast MRI data. The visualization of multidirectional blood flow revealed an increase in vortical flow in the arch, supra-aortic vessels, and descending aorta, and a correlation between the presence of aortic abnormalities and disturbed flow. Compared to the relatively homogeneous pattern of time-averaged wall shear stress (TAWSS) on the healthy aortae, a highly heterogeneous distribution with elevated TAWSS values was observed in the TS geometries. Visualization of further shear stress parameters, such as oscillatory shear index (OSI), normalized relative residence time (RRTn), and transverse WSS (transWSS), revealed dissimilar heterogeneity in the oscillatory and multidirectional nature of the aortic flow. Taking into account the young age of our TS cohort (average age 13 ± 2 years) and their obesity level (75% were obese or overweight), which is believed to accelerate the initiation and progression of endothelial dysfunction, these findings may be an indication of atherosclerotic disease manifesting earlier in life in TS patients. Age, obesity and aortic morphology may, therefore, play a key role in assessing cardiovascular risk in TS children.

Evaluation of Plaque Characteristics and Inflammation Using Magnetic Resonance Imaging

Atherosclerosis is an inflammatory disease of large and medium-sized arteries, characterized by the growth of atherosclerotic lesions (plaques). These plaques often develop at inner curvatures of arteries, branchpoints, and bifurcations, where the endothelial wall shear stress is low and oscillatory. In conjunction with other processes such as lipid deposition, biomechanical factors lead to local vascular inflammation and plaque growth. There is also evidence that low and oscillatory shear stress contribute to arterial remodeling, entailing a loss in arterial elasticity and, therefore, an increased pulse-wave velocity. Although altered shear stress profiles, elasticity and inflammation are closely intertwined and critical for plaque growth, preclinical and clinical investigations for atherosclerosis mostly focus on the investigation of one of these parameters only due to the experimental limitations. However, cardiovascular magnetic resonance imaging (MRI) has been demonstrated to be a potent tool which can be used to provide insights into a large range of biological parameters in one experimental session. It enables the evaluation of the dynamic process of atherosclerotic lesion formation without the need for harmful radiation. Flow-sensitive MRI provides the assessment of hemodynamic parameters such as wall shear stress and pulse wave velocity which may replace invasive and radiation-based techniques for imaging of the vascular function and the characterization of early plaque development. In combination with inflammation imaging, the analyses and correlations of these parameters could not only significantly advance basic preclinical investigations of atherosclerotic lesion formation and progression, but also the diagnostic clinical evaluation for early identification of high-risk plaques, which are prone to rupture. In this review, we summarize the key applications of magnetic resonance imaging for the evaluation of plaque characteristics through flow sensitive and morphological measurements. The simultaneous measurements of functional and structural parameters will further preclinical research on atherosclerosis and has the potential to fundamentally improve the detection of inflammation and vulnerable plaques in patients.

Changes in aortic arch geometry and the risk for Stanford B dissection

Background

The increase in aortic diameter is not closely associated with type B aortic dissection (TBAD); morphological risk factors other than aortic diameter may help to better identify patients at risk for TBAD. The purpose of this study was to investigate possible morphological factors associated with the occurrence of TBAD.

Methods

This study was a retrospective, multicenter, cross-sectional study. We collected 94 patients with TBAD as the TBAD group and 534 patients with healthy aortas as the healthy control group. Morphometric data were collected on three-dimensional models of the thoracic aorta. A propensity score matching (PSM) analysis was conducted to reduce the potential for confounding by baseline factors.

Results

The number of patients in the TBAD group was 75 after PSM. Longer lengths of the aortic arch (28.00±7.42 vs. 25.14±7.11 cm) were observed in patients with TBAD. The width (80.04±17.27 vs. 71.73±15.55 mm) and height (24.92±11.39 vs. 19.37±10.10 mm) of the aortic arch in patients with TBAD were both larger than those of healthy controls. The morphological changes associated with the occurrence of type B acute dissection were most pronounced in the geometry of the aortic arch.

Conclusions

This study demonstrates that TBAD was associated with longer lengths of aortic arch and with larger arch height and width.

3D confocal microscope imaging of macromolecule uptake in the intact brachiocephalic artery

BACKGROUND AND AIMS

Elevated uptake of plasma macromolecules by the arterial wall is an early event in atherogenesis. Existing optical techniques for detecting macromolecular tracers in the wall have poor depth penetration and hence require en face imaging of flattened arterial segments. Imaging uptake in undistorted curved and branched vessels would be useful in understanding disease development.

METHODS

Depth penetration was increased by applying optical clearing techniques. The rat aorto-brachiocephalic junction was imaged intact by confocal microscopy after it had been exposed to circulating rhodamine-labelled albumin in vivo, fixed in situ, excised and then cleared with benzyl alcohol/benzyl benzoate. Tracer uptake was mapped onto a 3D surface mesh of the arterial geometry.

RESULTS

Tracer fluorescence was detectable throughout the wall closest to the objective lens and, despite a vessel diameter of c. 1 mm, in the wall on the other side of the artery, across the lumen. By tile scanning, tracer concentrations were mapped in the aorta, the brachiocephalic artery and their junction without opening or flattening either vessel. Optical clearing was also shown to be compatible with immunofluorescent staining and imaging of experimental atherosclerosis.

CONCLUSIONS

The technique obviates the need for labour-intensive sample preparation associated with standard en face imaging. More importantly, it preserves arterial geometry, facilitating co-localisation of uptake maps with maps of biomechanical factors, which typically exist on 3D surface meshes. It will permit the correlation of haemodynamic wall shear stress with macromolecule permeability more accurately in regions of high curvature or branching, such as in the coronary arteries.

Blood flow simulations in the pulmonary bifurcation in relation to adult patients with repaired tetralogy of Fallot

Understanding the haemodynamic environment of the pulmonary bifurcation is important in adults with repaired conotruncal congenital heart disease. In these patients, dysfunction of the pulmonary valve and narrowing of the branch pulmonary arteries are common and can have serious clinical consequences. The aim of this study was to numerically investigate the underlying blood flow characteristics in the pulmonary trunk under a range of simplified conditions. For that, an in-depth analysis was conducted in idealised two-dimensional geometries that facilitate parametric investigation of healthy and abnormal conditions. Subtle variations in morphology influenced the haemodynamic environment and wall shear stress distribution. The pressure in the left pulmonary artery was generally higher than that in the right and main arteries, but was markedly reduced in the presence of a local stenosis. Different downstream pressure conditions altered the branch flow ratio, from 50:50% to more realistic 60:40% ratios in the right and left pulmonary artery, respectively. Despite some simplifications, this study highlights some previously undocumented aspects of the flow in bifurcating geometries, by clarifying the role of the stagnation point location on wall shear stress and differential branch pressures. In addition, measurements of the mean pressure ratios in the pulmonary bifurcation are discussed in the context of a new haemodynamic index which could potentially contribute to the assessment of left pulmonary artery stenosis in tetralogy of Fallot patients. Further studies are required to confirm the results in patient-specific models with personalised physiological flow conditions.

A key role for the novel coronary artery disease gene JCAD in atherosclerosis via shear stress mechanotransduction

AIMS

Genome-wide association studies (GWAS) have consistently identified an association between coronary artery disease (CAD) and a locus on chromosome 10 containing a single gene, JCAD (formerly KIAA1462). However, little is known about the mechanism by which JCAD could influence the development of atherosclerosis.

METHODS AND RESULTS

Vascular function was quantified in subjects with CAD by flow-mediated dilatation (FMD) and vasorelaxation responses in isolated blood vessel segments. The JCAD risk allele identified by GWAS was associated with reduced FMD and reduced endothelial-dependent relaxations. To study the impact of loss of Jcad on atherosclerosis, Jcad-/- mice were crossed to an ApoE-/- background and fed a high-fat diet from 6 to16 weeks of age. Loss of Jcad did not affect blood pressure or heart rate. However, Jcad-/-ApoE-/- mice developed significantly less atherosclerosis in the aortic root and the inner curvature of the aortic arch. En face analysis revealed a striking reduction in pro-inflammatory adhesion molecules at sites of disturbed flow on the endothelial cell layer of Jcad-/- mice. Loss of Jcad lead to a reduced recovery perfusion in response to hind limb ischaemia, a model of altered in vivo flow. Knock down of JCAD using siRNA in primary human aortic endothelial cells significantly reduced the response to acute onset of flow, as evidenced by reduced phosphorylation of NF-КB, eNOS, and Akt.

CONCLUSION

The novel CAD gene JCAD promotes atherosclerotic plaque formation via a role in the endothelial cell shear stress mechanotransduction pathway.

Computational study of aortic hemodynamics for patients with an abnormal aortic valve: The importance of secondary flow at the ascending aorta inlet

Blood flow in the aorta is helical, but most computational studies ignore the presence of secondary flow components at the ascending aorta (AAo) inlet. The aim of this study is to ascertain the importance of inlet boundary conditions (BCs) in computational analysis of flow patterns in the thoracic aorta based on patient-specific images, with a particular focus on patients with an abnormal aortic valve. Two cases were studied: one presenting a severe aortic valve stenosis and the other with a mechanical valve. For both aorta models, three inlet BCs were compared; these included the flat profile and 1D through-plane velocity and 3D phase-contrast magnetic resonance imaging derived velocity profiles, with the latter being used for benchmarking. Our results showed that peak and mean velocities at the proximal end of the ascending aorta were underestimated by up to 41% when the secondary flow components were neglected. The results for helical flow descriptors highlighted the strong influence of secondary velocities on the helical flow structure in the AAo. Differences in all wall shear stress (WSS)-derived indices were much more pronounced in the AAo and aortic arch (AA) than in the descending aorta (DAo). Overall, this study demonstrates that using 3D velocity profiles as inlet BC is essential for patient-specific analysis of hemodynamics and WSS in the AAo and AA in the presence of an abnormal aortic valve. However, predicted flow in the DAo is less sensitive to the secondary velocities imposed at the inlet; hence, the 1D through-plane profile could be a sufficient inlet BC for studies focusing on distal regions of the thoracic aorta.

Vascular Remodelling Relates to an Elevated Oscillatory Shear Index and Relative Residence Time in Spontaneously Hypertensive Rats

Haemodynamic disorders are common clinical findings in hypertension and lead to adverse cardiovascular events. However, the haemodynamic conditions in hypertension models are poorly understood. This study aimed to observe the characteristics of haemodynamics in spontaneously hypertensive rats (SHRs) and antihypertensive-treated SHRs. Twenty-four adult male SHRs and Wistar-Kyoto rats (WKYs) were randomly divided into four groups and treated for 7 days as follows: WKY-CON (WKYs + saline), WKY-NIF (WKYs + nifedipine, 50 mg/kg/day), SHR-CON (SHRs + saline), and SHR-NIF (SHRs + nifedipine). Aortic computational fluid dynamics (CFD) models were simulated to obtain the haemodynamic parameters. We found that in the hypertensive (SHR-CON) and blood pressure-controlled (SHR-NIF) groups, the oscillatory shear index (OSI) and relative residence time (RRT), which are key haemodynamics indices, were markedly elevated. Furthermore, there was a correlation between both the elevated OSI and RRT with the vascular wall thickening in regions near the inner wall of the aortic arch. Our research demonstrates that haemodynamics remains disturbed even if the blood pressure is normalized. In addition, vascular remodelling may play an important role in maintaining elevated OSI and RRT values.

Validation of numerical simulation methods in aortic arch using 4D Flow MRI

Computational fluid dynamics (CFD) are the gold standard in studying blood flow dynamics. However, CFD results are dependent on the boundary conditions and the computation model. The purpose of this study was to validate CFD methods using comparison with actual measurements of the blood flow vector obtained with four-dimensional (4D) flow magnetic resonance imaging (MRI). 4D Flow MRI was performed on a healthy adult and a child with double-aortic arch. The aortic lumen was segmented to visualize the blood flow. The CFD analyses were performed for the same geometries based on three turbulent models: laminar, large eddy simulation (LES), and the renormalization group k–ε model (RNG k–ε). The flow-velocity vector components, namely the wall shear stress (WSS) and flow energy loss (EL), of the MRI and CFD results were compared. The flow rate of the MRI results was underestimated in small vessels, including the neck vessels. Spiral flow in the ascending aorta caused by the left ventricular twist was observed by MRI. Secondary flow distal to the aortic arch was well realized in both CFD and MRI. The average correlation coefficients of the velocity vector components of MRI and CFD for the child were the highest for the RNG k–ε model (0.530 in ascending aorta, 0.768 in the aortic arch, 0.584 in the descending aorta). The WSS and EL values of MRI were less than half of those of CFD, but the WSS distribution patterns were quite similar. The WSS and EL estimates were higher in RNG k–ε and LES than in the laminar model because of eddy viscosity. The CFD computation realized accurate flow distal to the aortic arch, and the WSS distribution was well simulated compared to actual measurement using 4D Flow MRI. However, the helical flow was not simulated in the ascending aorta. The accuracy was enhanced by using the turbulence model, and the RNG k–ε model showed the highest correlation with 4D Flow MRI.

Experimental Approaches to Study Endothelial Responses to Shear Stress

SIGNIFICANCE

Shear stress controls multiple physiological processes in endothelial cells (ECs).

RECENT ADVANCES

The response of ECs to shear has been studied using a range of in vitro and in vivo models.

CRITICAL ISSUES

This article describes some of the experimental techniques that can be used to study endothelial responses to shear stress. It includes an appraisal of large animal, rodent, and zebrafish models of vascular mechanoresponsiveness. It also describes several bioreactors to apply flow to cells and physical methods to separate mechanoresponses from mass transport mechanisms.

FUTURE DIRECTIONS

We conclude that combining in vitro and in vivo approaches can provide a detailed mechanistic view of vascular responses to force and that high-throughput systems are required for unbiased assessment of the function of shear-induced molecules. Antioxid. Redox Signal. 25, 389-400.

Heart rate reduction with ivabradine promotes shear stress-dependent anti-inflammatory mechanisms in arteries

Blood flow generates wall shear stress (WSS) which alters endothelial cell (EC) function. Low WSS promotes vascular inflammation and atherosclerosis whereas high uniform WSS is protective. Ivabradine decreases heart rate leading to altered haemodynamics. Besides its cardio-protective effects, ivabradine protects arteries from inflammation and atherosclerosis via unknown mechanisms. We hypothesised that ivabradine protects arteries by increasing WSS to reduce vascular inflammation. Hypercholesterolaemic mice were treated with ivabradine for seven weeks in drinking water or remained untreated as a control. En face immunostaining demonstrated that treatment with ivabradine reduced the expression of pro-inflammatory VCAM-1 (p<0.01) and enhanced the expression of anti-inflammatory eNOS (p<0.01) at the inner curvature of the aorta. We concluded that ivabradine alters EC physiology indirectly via modulation of flow because treatment with ivabradine had no effect in ligated carotid arteries in vivo, and did not influence the basal or TNFα-induced expression of inflammatory (VCAM-1, MCP-1) or protective (eNOS, HMOX1, KLF2, KLF4) genes in cultured EC. We therefore considered whether ivabradine can alter WSS which is a regulator of EC inflammatory activation. Computational fluid dynamics demonstrated that ivabradine treatment reduced heart rate by 20 % and enhanced WSS in the aorta. In conclusion, ivabradine treatment altered haemodynamics in the murine aorta by increasing the magnitude of shear stress. This was accompanied by induction of eNOS and suppression of VCAM-1, whereas ivabradine did not alter EC that could not respond to flow. Thus ivabradine protects arteries by altering local mechanical conditions to trigger an anti-inflammatory response.

Plectin reinforces vascular integrity by mediating crosstalk between the vimentin and the actin networks

Mutations in the cytoskeletal linker protein plectin result in multisystemic diseases affecting skin and muscle with indications of additional vascular system involvement. To study the mechanisms underlying vascular disorders, we established plectin-deficient endothelial cell and mouse models. We show that apart from perturbing the vimentin cytoskeleton of endothelial cells, plectin deficiency leads to severe distortions of adherens junctions (AJs), as well as tight junctions, accompanied by an upregulation of actin stress fibres and increased cellular contractility. Plectin-deficient endothelial cell layers were more leaky and showed reduced mechanical resilience in fluid-shear stress and mechanical stretch experiments. We suggest that the distorted AJs and upregulated actin stress fibres in plectin-deficient cells are rooted in perturbations of the vimentin cytoskeleton, as similar phenotypes could be mimicked in wild-type cells by disruption of vimentin filaments. In vivo studies in endothelium-restricted conditional plectin-knockout mice revealed significant distortions of AJs in stress-prone aortic arch regions and increased pulmonary vascular leakage. Our study opens a new perspective on cytoskeleton-controlled vascular permeability, where a plectin-organized vimentin scaffold keeps actomyosin contractility ‘in-check’ and maintains AJ homeostasis.

Summary: Plectin-arranged vimentin scaffolds keep actomyosin contractility ‘in-check’ and maintain cell–cell junction homeostasis, providing a new perspective on cytoskeleton-controlled vascular permeability.

Inputs for Subject-Specific Computational Fluid Dynamics Simulation of Blood Flow in the Mouse Aorta

Mouse models are an important way for exploring relationships between blood hemodynamics and eventual plaque formation. We have developed a mouse model of aortic regurgitation (AR) that produces large changes in plaque burden with charges in hemodynamics [Zhou et al., 2010, "Aortic Regurgitation Dramatically Alters the Distribution of Atherosclerotic Lesions and Enhances Atherogenesis in Mice," Arterioscler. Thromb. Vasc. Biol., 30(6), pp. 1181–1188]. In this paper, we explore the amount of detail needed for realistic computational fluid dynamics (CFD) calculations in this experimental model. The CFD calculations use inputs based on experimental measurements from ultrasound (US), micro computed tomography (CT), and both anatomical magnetic resonance imaging (MRI) and phase contrast MRI (PC-MRI). The adequacy of five different levels of model complexity (a) subject-specific CT data from a single mouse; (b) subject-specific CT centerlines with radii from US; (c) same as (b) but with MRI derived centerlines; (d) average CT centerlines and averaged vessel radius and branching vessels; and (e) same as (d) but with averaged MRI centerlines) is evaluated by demonstrating their impact on relative residence time (RRT) outputs. The paper concludes by demonstrating the necessity of subject-specific geometry and recommends for inputs the use of CT or anatomical MRI for establishing the aortic centerlines, M-mode US for scaling the aortic diameters, and a combination of PC-MRI and Doppler US for estimating the spatial and temporal characteristics of the input wave forms.

Measuring, reversing, and modeling the mechanical changes due to the absence of Fibulin-4 in mouse arteries

Mice with a smooth muscle cell (SMC)-specific deletion of Fibulin-4 (SMKO) show decreased expression of SMC contractile genes, decreased circumferential compliance, and develop aneurysms in the ascending aorta. Neonatal administration of drugs that inhibit the angiotensin II pathway encourages the expression of contractile genes and prevents aneurysm development, but does not increase compliance in SMKO aorta. We hypothesized that multidimensional mechanical changes in the aorta and/or other elastic arteries may contribute to aneurysm pathophysiology. We found that the SMKO ascending aorta and carotid artery showed mechanical changes in the axial direction. These changes were not reversed by angiotensin II inhibitors, hence reversing the axial changes is not required for aneurysm prevention. Mechanical changes in the circumferential direction were specific to the ascending aorta; therefore, mechanical changes in the carotid do not contribute to aortic aneurysm development. We also hypothesized that a published model of postnatal aortic growth and remodeling could be used to investigate mechanisms behind the changes in SMKO aorta and aneurysm development over time. Dimensions and mechanical behavior of adult SMKO aorta were reproduced by the model after modifying the initial component material constants and the aortic dilation with each postnatal time step. The model links biological observations to specific mechanical responses in aneurysm development and treatment.

The impact of MRI-based inflow for the hemodynamic evaluation of aortic coarctation

Aortic coarctation (CoA) accounting for 3-11% of congenital heart disease can be successfully treated. Long-term results, however, have revealed decreased life expectancy associated with abnormal hemodynamics. Accordingly, an assessment of hemodynamics is the key factor in treatment decisions and successful long-term results. In this study, 3D angiography whole heart (3DWH) and 4D phase-contrast magnetic resonance imaging (MRI) data were acquired. Geometries of the thoracic aorta with CoAs were reconstructed using ZIB-Amira software. X-ray angiograms were used to evaluate the post-treatment geometry. Computational fluid dynamics models in three patients were created to simulate pre- and post-treatment situations using the FLUENT program. The aim of the study was to investigate the impact of the inlet velocity profile (plug vs. MRI-based) with a focus on the peak systole pressure gradient and wall shear stress (WSS). Results show that helical flow at the aorta inlet can significantly affect the assessment of pressure drop and WSS. Simplified plug inlet velocity profiles significantly (p < 0.05) overestimate the pressure drop in pre- and post-treatment geometries and significantly (p < 0.05) underestimate surface-averaged WSS. We conclude that the use of the physiologically correct but time-expensive 4D MRI-based in vivo velocity profile in CFD studies may be an important step towards a patient-specific analysis of CoA hemodynamics.

Co-registration of the distribution of wall shear stress and 140 complex plaques of the aorta

Previous studies provide evidence that atherosclerosis develops in vascular regions exposed to low wall shear stress (WSS) and high oscillatory shear index (OSI). 4D flow MRI was analyzed in 70 stroke patients with complex plaques (≥4 mm thickness, ulcerated or superimposed thrombi) and in 12 young healthy volunteers. The segmental distribution of peak systolic WSSsystole and OSI was quantified in analysis planes positioned directly at the location of 140 complex plaques found in the 70 patients. In addition, WSSsystole and OSI were evaluated in 8 standard analysis planes distributed along the aorta. Complex plaques were predominantly found at the inner curvature of the aortic arch and of the descending aorta. High OSI was co-located with the segments mostly affected by complex plaque while WSSsystole demonstrated a homogenous distribution. In standard analysis planes, patients demonstrated significantly (p<0.01) altered distribution of wall parameters compared to volunteers (reduced WSSsystole in 91% of aortic wall segments, increased OSI in 48% of segments). OSI distribution was asymmetric with higher values at the inner curvature of the aorta. While WSS and OSI showed expected changes in patients compared to healthy controls, their distribution pattern at complex plaques indicated a more complex and heterogeneous relationship than previously anticipated.