Towards the improved quantification of in vivo abnormal wall shear stresses in BAV-affected patients from 4D-flow imaging: Benchmarking and application to real data

Left ventricle diastolic vortex ring characterization in ischemic cardiomyopathy: insight into atrio-ventricular interplay

Diastolic vortex ring (VR) plays a key role in the blood-pumping function exerted by the left ventricle (LV), with altered VR structures being associated with LV dysfunction. Herein, we sought to characterize the VR diastolic alterations in ischemic cardiomyopathy (ICM) patients with systo-diastolic LV dysfunction, as compared to healthy controls, in order to provide a more comprehensive understanding of LV diastolic function. 4D Flow MRI data were acquired in ICM patients (n = 15) and healthy controls (n = 15). The λ2 method was used to extract VRs during early and late diastolic filling. Geometrical VR features, e.g., circularity index (CI), orientation (α), and inclination with respect to the LV outflow tract (ß), were extracted. Kinetic energy (KE), rate of viscous energy loss ( EL ˙ ), vorticity (W), and volume (V) were computed for each VR; the ratios with the respective quantities computed for the entire LV were derived. At peak E-wave, the VR was less circular (p = 0.032), formed a smaller α with the LV long-axis (p = 0.003) and a greater ß (p = 0.002) in ICM patients as compared to controls. At peak A-wave, CI was significantly increased (p = 0.034), while α was significantly smaller (p = 0.016) and β was significantly increased (p = 0.036) in ICM as compared to controls. At both peak E-wave and peak A-wave,EL ˙ VR / EL ˙ LV , WVR/WLV, and VVR/VLV significantly decreased in ICM patients vs. healthy controls. KEVR/VVR showed a significant decrease in ICM patients with respect to controls at peak E-wave, while VVR remained comparable between normal and pathologic conditions. In the analyzed ICM patients, the diastolic VRs showed alterations in terms of geometry and energetics. These derangements might be attributed to both structural and functional alterations affecting the infarcted wall region and the remote myocardium.

Does the AVNeo valve reduce wall stress on the aortic wall? A cardiac magnetic resonance analysis with 4D-flow for the evaluation of aortic valve replacement with the Ozaki technique

OBJECTIVES

Aortic valve neocuspidalization aims to replace the 3 aortic cusps with autologous pericardium pre-treated with glutaraldehyde, and it is a surgical alternative to the classical aortic valve replacement (AVR). Image-based patient-specific computational fluid dynamics allows the derivation of shear stress on the aortic wall [wall shear stress (WSS)]. Previous studies support a potential link between increased WSS and histological alterations of the aortic wall. The aim of this study is to compare the WSS of the ascending aorta in patients undergoing aortic valve neocuspidalization versus AVR with biological prostheses.

METHODS

This is a prospective nonrandomized clinical trial. Each patient underwent a 4D-flow cardiac magnetic resonance scan after surgery, which informed patient-specific computational fluid dynamics models to evaluate WSS at the ascending aortic wall. The adjusted variables were calculated by summing the residuals obtained from a multivariate linear model (with ejection fraction and left ventricle outflow tract-aorta angle as covariates) to the mean of the variables.

RESULTS

Ten patients treated with aortic valve neocuspidalization were enrolled and compared with 10 AVR patients. The aortic valve neocuspidalization group showed a significantly lower WSS in the outer curvature segments of the proximal and distal ascending aorta as compared to AVR patients (P = 0.0179 and 0.0412, respectively). WSS levels remained significantly lower along the outer curvature of the proximal aorta in the aortic valve neocuspidalization population, even after adjusting the WSS for the ejection fraction and the left ventricle outflow tract-aorta angle [2.44 Pa (2.17-3.01) vs 1.94 Pa (1.72-2.01), P = 0.02].

CONCLUSIONS

Aortic valve neocuspidalization hemodynamical features are potentially associated with a lower WSS in the ascending aorta as compared to commercially available bioprosthetic valves.

Fluid-structure interactions of peripheral arteries using a coupled in silico and in vitro approach

Vascular compliance is considered both a cause and a consequence of cardiovascular disease and a significant factor in the mid- and long-term patency of vascular grafts. However, the biomechanical effects of localised changes in compliance cannot be satisfactorily studied with the available medical imaging technologies or surgical simulation materials. To address this unmet need, we developed a coupled silico-vitro platform which allows for the validation of numerical fluid-structure interaction results as a numerical model and physical prototype. This numerical one-way and two-way fluid-structure interaction study is based on a three-dimensional computer model of an idealised femoral artery which is validated against patient measurements derived from the literature. The numerical results are then compared with experimental values collected from compliant arterial phantoms via direct pressurisation and ring tensile testing. Phantoms within a compliance range of 1.4-68.0%/100 mmHg were fabricated via additive manufacturing and silicone casting, then mechanically characterised via ring tensile testing and optical analysis under direct pressurisation with moderately statistically significant differences in measured compliance ranging between 10 and 20% for the two methods. One-way fluid-structure interaction coupling underestimated arterial wall compliance by up to 14.7% compared with two-way coupled models. Overall, Solaris™ (Smooth-On) matched the compliance range of the numerical and in vivo patient models most closely out of the tested silicone materials. Our approach is promising for vascular applications where mechanical compliance is especially important, such as the study of diseases which commonly affect arterial wall stiffness, such as atherosclerosis, and the model-based design, surgical training, and optimisation of vascular prostheses.

Segmentation of the Aorta and Pulmonary Arteries Based on 4D Flow MRI in the Pediatric Setting Using Fully Automated Multi-Site, Multi-Vendor, and Multi-Label Dense U-Net

BACKGROUND

Automated segmentation using convolutional neural networks (CNNs) have been developed using four-dimensional (4D) flow magnetic resonance imaging (MRI). To broaden usability for congenital heart disease (CHD), training with multi-institution data is necessary. However, the performance impact of heterogeneous multi-site and multi-vendor data on CNNs is unclear.

PURPOSE

To investigate multi-site CNN segmentation of 4D flow MRI for pediatric blood flow measurement.

STUDY TYPE

Retrospective.

POPULATION

A total of 174 subjects across two sites (female: 46%; N = 38 healthy controls, N = 136 CHD patients). Participants from site 1 (N = 100), site 2 (N = 74), and both sites (N = 174) were divided into subgroups to conduct 10-fold cross validation (10% for testing, 90% for training).

FIELD STRENGTH/SEQUENCE

3 T/1.5 T; retrospectively gated gradient recalled echo-based 4D flow MRI.

ASSESSMENT

Accuracy of the 3D CNN segmentations trained on data from single site (single-site CNNs) and data across both sites (multi-site CNN) were evaluated by geometrical similarity (Dice score, human segmentation as ground truth) and net flow quantification at the ascending aorta (Qs), main pulmonary artery (Qp), and their balance (Qp/Qs), between human observers, single-site and multi-site CNNs.

STATISTICAL TESTS

Kruskal-Wallis test, Wilcoxon rank-sum test, and Bland-Altman analysis. A P-value <0.05 was considered statistically significant.

RESULTS

No difference existed between single-site and multi-site CNNs for geometrical similarity in the aorta by Dice score (site 1: 0.916 vs. 0.915, P = 0.55; site 2: 0.906 vs. 0.904, P = 0.69) and for the pulmonary arteries (site 1: 0.894 vs. 0.895, P = 0.64; site 2: 0.870 vs. 0.869, P = 0.96). Qs site-1 medians were 51.0-51.3 mL/cycle (P = 0.81) and site-2 medians were 66.7-69.4 mL/cycle (P = 0.84). Qp site-1 medians were 46.8-48.0 mL/cycle (P = 0.97) and site-2 medians were 76.0-77.4 mL/cycle (P = 0.98). Qp/Qs site-1 medians were 0.87-0.88 (P = 0.97) and site-2 medians were 1.01-1.03 (P = 0.43). Bland-Altman analysis for flow quantification found equivalent performance.

DATA CONCLUSION

Multi-site CNN-based segmentation and blood flow measurement are feasible for pediatric 4D flow MRI and maintain performance of single-site CNNs.

LEVEL OF EVIDENCE

3 TECHNICAL EFFICACY: Stage 2.

A novel MRI-based data fusion methodology for efficient, personalised, compliant simulations of aortic haemodynamics

We present a novel, cost-efficient methodology to simulate aortic haemodynamics in a patient-specific, compliant aorta using an MRI data fusion process. Based on a previously-developed Moving Boundary Method, this technique circumvents the high computational cost and numerous structural modelling assumptions required by traditional Fluid-Structure Interaction techniques. Without the need for Computed Tomography (CT) data, the MRI images required to construct the simulation can be obtained during a single imaging session. Black Blood MR Angiography and 2D Cine-MRI data were used to reconstruct the luminal geometry and calibrate wall movement specifically to each region of the aorta. 4D-Flow MRI and non-invasive pressure measurements informed patient-specific inlet and outlet boundary conditions. Luminal area closely matched 2D Cine-MRI measurements with a mean error of less than 4.6% across the cardiac cycle, while physiological pressure and flow distributions were simulated to within 3.3% of patient-specific targets. Moderate agreement with 4D-Flow MRI velocity data was observed. Despite lower peak velocity, an equivalent rigid-wall simulation predicted a mean Time-Averaged Wall Shear Stress (TAWSS) 13% higher than the compliant simulation. The agreement observed between compliant simulation results and MRI data is testament to the accuracy and efficiency of this MRI-based simulation technique.

Significance of Hemodynamics Biomarkers, Tissue Biomechanics and Numerical Simulations in the Pathogenesis of Ascending Thoracic Aortic Aneurysms

Guidelines for the treatment of aortic wall diseases are based on measurements of maximum aortic diameter. However, aortic rupture or dissections do occur for small aortic diameters. Growing scientific evidence underlines the importance of biomechanics and hemodynamics in aortic disease development and progression. Wall shear stress (WWS) is an important hemodynamics marker that depends on aortic wall morphology and on the aortic valve function. WSS could be helpful to interpret aortic wall remodeling and define personalized risk criteria. The complementarity of Computational Fluid Dynamics and 4D Magnetic Resonance Imaging as tools for WSS assessment is a promising reality. The potentiality of these innovative technologies will provide maps or atlases of hemodynamics biomarkers to predict aortic tissue dysfunction. Ongoing efforts should focus on the correlation between these non-invasive imaging biomarkers and clinico-pathologic situations for the implementation of personalized medicine in current clinical practice.

Molecular and Mechanical Mechanisms of Calcification Pathology Induced by Bicuspid Aortic Valve Abnormalities

Bicuspid aortic valve (BAV) is a congenital defect affecting 1-2% of the general population that is distinguished from the normal tricuspid aortic valve (TAV) by the existence of two, rather than three, functional leaflets (or cusps). BAV presents in different morphologic phenotypes based on the configuration of cusp fusion. The most common phenotypes are Type 1 (containing one raphe), where fusion between right coronary and left coronary cusps (BAV R/L) is the most common configuration followed by fusion between right coronary and non-coronary cusps (BAV R/NC). While anatomically different, BAV R/L and BAV R/NC configurations are both associated with abnormal hemodynamic and biomechanical environments. The natural history of BAV has shown that it is not necessarily the primary structural malformation that enforces the need for treatment in young adults, but the secondary onset of premature calcification in ~50% of BAV patients, that can lead to aortic stenosis. While an underlying genetic basis is a major pathogenic contributor of the structural malformation, recent studies have implemented computational models, cardiac imaging studies, and bench-top methods to reveal BAV-associated hemodynamic and biomechanical alterations that likely contribute to secondary complications. Contributions to the field, however, lack support for a direct link between the external valvular environment and calcific aortic valve disease in the setting of BAV R/L and R/NC BAV. Here we review the literature of BAV hemodynamics and biomechanics and discuss its previously proposed contribution to calcification. We also offer means to improve upon previous studies in order to further characterize BAV and its secondary complications.

Parametric Sequential Method for MRI-Based Wall Shear Stress Quantification

Wall shear stress (WSS) has been suggested as a potential biomarker in various cardiovascular diseases and it can be estimated from phase-contrast Magnetic Resonance Imaging (PC-MRI) velocity measurements. We present a parametric sequential method for MRI-based WSS quantification consisting of a geometry identification and a subsequent approximation of the velocity field. This work focuses on its validation, investigating well controlled high-resolution in vitro measurements of turbulent stationary flows and physiological pulsatile flows in phantoms. Initial tests for in vivo 2D PC-MRI data of the ascending aorta of three volunteers demonstrate basic applicability of the method to in vivo.

An LDV based method to quantify the error of PC-MRI derived Wall Shear Stress measurement

Wall Shear Stress (WSS) has been demonstrated to be a biomarker of the development of atherosclerosis. In vivo assessment of WSS is still challenging, but 4D Flow MRI represents a promising tool to provide 3D velocity data from which WSS can be calculated. In this study, a system based on Laser Doppler Velocimetry (LDV) was developed to validate new improvements of 4D Flow MRI acquisitions and derived WSS computing. A hydraulic circuit was manufactured to allow both 4D Flow MRI and LDV velocity measurements. WSS profiles were calculated with one 2D and one 3D method. Results indicated an excellent agreement between MRI and LDV velocity data, and thus the set-up enabled the evaluation of the improved performances of 3D with respect to the 2D-WSS computation method. To provide a concrete example of the efficacy of this method, the influence of the spatial resolution of MRI data on derived 3D-WSS profiles was investigated. This investigation showed that, with acquisition times compatible with standard clinical conditions, a refined MRI resolution does not improve WSS assessment, if the impact of noise is unreduced. This study represents a reliable basis to validate with LDV WSS calculation methods based on 4D Flow MRI.

4D flow evaluation of blood non-Newtonian behavior in left ventricle flow analysis

Blood is generally modeled as a Newtonian fluid, assuming a standard and constant viscosity; however, this assumption may not hold for the highly pulsatile and recirculating intracavitary flow in the left ventricle (LV), hampering the quantification of fluid dynamic indices of potential clinical relevance. Herein, we investigated the effect of three viscosity models on the patient-specific quantification of LV blood energetics, namely on viscous energy loss (EL), from 4D Flow magnetic resonance imaging: I) Newtonian with standard viscosity (3.7 cP), II) Newtonian with subject-specific hematocrit-dependent viscosity, III) non-Newtonian accounting for the effect of hematocrit and shear rate. Analyses were performed on 5 controls and 5 patients with cardiac light-chain amyloidosis. In Model II, viscosity ranged between 3.0 (-19%) and 4.3 cP (+16%), mildly deviating from the standard value. In the non-Newtonian model, this effect was emphasized: viscosity ranged from 3.2 to 6.0 cP, deviating maximally from the standard value in low shear rate (i.e., <100 s^-1) regions. This effect reflected on EL quantifications: in particular, as compared to Model I, Model III yielded markedly higher EL values (up to +40%) or markedly lower (down to -21%) for subjects with hematocrit higher than 39.5% and lower than 30%, respectively. Accounting for non-Newtonian blood behavior on a patient-specific basis may enhance the accuracy of intracardiac energetics assessment by 4D Flow, which may be explored as non-invasive index to discriminate between healthy and pathologic LV.

Towards quantitative evaluation of wall shear stress from 4D flow imaging

Wall shear stress (WSS) is a relevant hemodynamic indicator of the local stress applied on the endothelium surface. More specifically, its spatiotemporal distribution reveals crucial in the evolution of many pathologies such as aneurysm, stenosis, and atherosclerosis. This paper introduces a new solution, called PaLMA, to quantify the WSS from 4D Flow MRI data. It relies on a two-step local parametric model, to accurately describe the vessel wall and the velocity-vector field in the neighborhood of a given point of interest. Extensive validations have been performed on synthetic 4D Flow MRI data, including four datasets generated from patient specific computational fluid dynamics simulations on carotids. The validation tests are focused on the impact of the noise component, of the resolution level, and of the segmentation accuracy concerning the vessel position in the context of complex flow patterns. In simulated cases aimed to reproduce clinical acquisition conditions, the WSS quantification performance reached by PaLMA is significantly higher (with a gain in RMSE of 12 to 27%) than the reference one obtained using the smoothing B-spline method proposed by Potters et al. (2015) method, while the computation time is equivalent for both WSS quantification methods.

Four-dimensional Flow Magnetic Resonance Imaging Quantification of Blood Flow in Bicuspid Aortic Valve

BACKGROUND

Four-dimensional (D) flow magnetic resonance imaging (MRI) is limited by time-consuming and nonstandardized data analysis. We aimed to test the efficiency and interobserver reproducibility of a dedicated 4D flow MRI analysis workflow.

MATERIALS AND METHODS

Thirty retrospectively identified patients with bicuspid aortic valve (BAV, age=47.8±11.8 y, 9 male) and 30 healthy controls (age=48.8±12.5 y, 21 male) underwent Aortic 4D flow MRI using 1.5 and 3 T MRI systems. Two independent readers performed 4D flow analysis on a dedicated workstation including preprocessing, aorta segmentation, and placement of four 2D planes throughout the aorta for quantification of net flow, peak velocity, and regurgitant fraction. 3D flow visualization using streamlines was used to grade aortic valve outflow jets and extent of helical flow.

RESULTS

4D flow analysis workflow time for both observers: 5.0±1.4 minutes per case (range=3 to 10 min). Valve outflow jets and flow derangement was visible in all 30 BAV patients (both observers). Net flow, peak velocity, and regurgitant fraction was significantly elevated in BAV patients compared with controls except for regurgitant fraction in plane 4 (91.1±29.7 vs. 62.6±19.6 mL/s, 37.1% difference; 121.7±49.7 vs. 90.9±26.4 cm/s, 28.9% difference; 9.3±10.1% vs. 2.0±3.4%, 128.0% difference, respectively; P<0.001). Excellent intraclass correlation coefficient agreement for net flow: 0.979, peak velocity: 0.931, and regurgitant fraction: 0.928.

CONCLUSION

Our study demonstrates the potential of an efficient data analysis workflow to perform standardized 4D flow MRI processing in under 10 minutes and with good-to-excellent reproducibility for flow and velocity quantification in the thoracic aorta.

Evaluation of 4D flow MRI-based non-invasive pressure assessment in aortic coarctations

Severity of aortic coarctation (CoA) is currently assessed by estimating trans-coarctation pressure drops through cardiac catheterization or echocardiography. In principle, more detailed information could be obtained non-invasively based on space- and time-resolved magnetic resonance imaging (4D flow) data. Yet the limitations of this imaging technique require testing the accuracy of 4D flow-derived hemodynamic quantities against other methodologies. With the objective of assessing the feasibility and accuracy of this non-invasive method to support the clinical diagnosis of CoA, we developed an algorithm (4DF-FEPPE) to obtain relative pressure distributions from 4D flow data by solving the Poisson pressure equation. 4DF-FEPPE was tested against results from a patient-specific fluid-structure interaction (FSI) simulation, whose patient-specific boundary conditions were prescribed based on 4D flow data. Since numerical simulations provide noise-free pressure fields on fine spatial and temporal scales, our analysis allowed to assess the uncertainties related to 4D flow noise and limited resolution. 4DF-FEPPE and FSI results were compared on a series of cross-sections along the aorta. Bland-Altman analysis revealed very good agreement between the two methodologies in terms of instantaneous data at peak systole, end-diastole and time-averaged values: biases (means of differences) were +0.4 mmHg, -1.1 mmHg and +0.6 mmHg, respectively. Limits of agreement (2 SD) were ±0.978 mmHg, ±1.06 mmHg and ±1.97 mmHg, respectively. Peak-to-peak and maximum trans-coarctation pressure drops obtained with 4DF-FEPPE differed from FSI results by 0.75 mmHg and -1.34 mmHg respectively. The present study considers important validation aspects of non-invasive pressure difference estimation based on 4D flow MRI, showing the potential of this technology to be more broadly applied to the clinical practice.

Novel insights by 4D Flow imaging on aortic flow physiology after valve-sparing root replacement with or without neosinuses

OBJECTIVES

This study was undertaken to evaluate the flow dynamics in the aortic root after valve-sparing root replacement with and without neosinuses of Valsalva reconstruction, by exploiting the capability of 4D Flow imaging to measure in vivo blood velocity fields and 3D geometric flow patterns.

METHODS

Ten patients who underwent valve-sparing root replacement utilizing grafts with neosinuses or straight tube grafts (5 cases each) were evaluated by 4D Flow imaging at a mean of 46.5 months after surgery. We used in-house processing tools to quantify relevant bulk flow variables (flow rate, stroke volume, peak velocity and mean velocity), wall shear stresses and the amount of flow rotation characterizing the region enclosed by the graft and the aortic valve leaflets.

RESULTS

Despite bulk flows with similar peak velocities, flow rates and stroke volumes (P = 0.31-1.00), the neosinuses graft was associated with a lower mean velocity (P < 0.03) and magnitude of wall shear stress along the axial direction of the vessel wall (P < 0.05) at the proximal root level but remained comparable along the circumferential direction (P = 0.22-1.0) to the straight tube graft. Flow rotation was evidently and systematically higher in the neosinuses grafts, characterized by streamline rotations higher than 270°, nearly triple that of tubular grafts (10.3 ÷ 14.0% of all aortic streamline vs 2.2 ÷ 5.7%, P = 0.008).

CONCLUSIONS

Recreation of the sinuses of Valsalva during valve-sparing root replacement is associated with significantly lower wall shear stress and organized vortical flows at the level of the sinus that are not evident using the straight tube graft. These findings need confirmation in larger studies and could have important implications in terms of aortic valve durability.

New imaging techniques project the cellular and molecular alterations underlying bicuspid aortic valve development

Bicuspid aortic valve (BAV) disease is the most common congenital cardiac malformation associated with an increased lifetime risk and a high rate of surgically-relevant valve deterioration and aortic dilatation. Genomic data revealed that different genes are associated with BAV. A dominant genetic factor for the recent past was the basis to the recommendation for a more extensive aortic intervention. However very recent evidence that hemodynamic stressors and alterations of wall shear stress play an important role independent from the genetic trait led to more conservative treatment recommendations. Therefore, there is a current need to improve the ability to risk stratify BAV patients in order to obtain an early detection of valvulopathy and aortopathy while also to predict valve dysfunction and/or aortic disease development. Imaging studies based on new cutting-edge technologies, such us 4-dimensional (4D) flow magnetic resonance imaging (MRI), two-dimensional (2D) or three-dimensional (3D) speckle-tracking imaging (STI) and computation fluid dynamics, combined with studies demonstrating new gene mutations, specific signal pathways alterations, hemodynamic influences, circulating biomarkers modifications, endothelial progenitor cell impairment and immune/inflammatory response, all detected BAV valvulopathy progression and aortic wall abnormality. Overall, the main purpose of this review article is to merge the evidences of imaging and basic science studies in a coherent hypothesis that underlies and thus projects the development of both BAV during embryogenesis and BAV-associated aortopathy and its complications in the adult life, with the final goal to identifying aneurysm formation/rupture susceptibility to improve diagnosis and management of patients with BAV-related aortopathy.

Characterization and estimation of turbulence-related wall shear stress in patient-specific pulsatile blood flow

Disturbed, turbulent-like blood flow promotes chaotic wall shear stress (WSS) environments, impairing essential endothelial functions and increasing the susceptibility and progression of vascular diseases. These flow characteristics are today frequently detected at various anatomical, lesion and intervention-related sites, while their role as a pathological determinant is less understood. To present-day, numerous WSS-based descriptors have been proposed to characterize the spatiotemporal nature of the WSS disturbances, however, without differentiation between physiological laminar oscillations and turbulence-related WSS (tWSS) fluctuations. Also, much attention has been focused on magnetic resonance (MR) WSS estimations, so far with limited success; promoting the need of a near-wall surrogate marker. In this study, a new approach is explored to characterize the tWSS, by taking advantage of the tensor characteristics of the fluctuating WSS correlations, providing both a magnitude and an anisotropy measure of the disturbances. These parameters were studied in two patient-specific coarctation models (sever and mild), using large eddy simulations, and correlated against near-wall reciprocal Reynolds stress parameters. Collectively, results showed distinct regions of differing tWSS characteristics, features which were sensitive to changes in flow conditions. Generally, the post-stenotic tWSS was governed by near axisymmetric fluctuations, findings that where not consistent with conventional WSS disturbance predictors. At the 2-3 mm wall-offset range, a strong linear correlation was found between tWSS magnitude and near-wall turbulence kinetic energy (TKE), in contrast to the anisotropy indices, suggesting that MR-measured TKE can be used to assess elevated tWSS regions while tWSS anisotropy estimates request well-resolved simulation methods.

Wall shear stress estimation in the aorta: Impact of wall motion, spatiotemporal resolution, and phase noise

BACKGROUND

Wall shear stress (WSS) presents an important parameter for assessing blood flow characteristics and evaluating flow-mediated lesions in the aorta.

PURPOSE

To investigate the robustness of WSS and oscillatory shear index (OSI) estimation based on 4D flow MRI against vessel wall motion, spatiotemporal resolution, and velocity encoding (VENC).

STUDY TYPE

Simulated and prospective.

POPULATION

Synthetic 4D flow MRI data of the aorta, simulated using the Lattice-Boltzmann method; in vivo 4D flow MRI data of the aorta from healthy volunteers (n = 11) and patients with congenital heart defects (n = 17).

FIELD STRENGTH/SEQUENCE

1.5T; 4D flow MRI with PEAK-GRAPPA acceleration and prospective electrocardiogram triggering.

ASSESSMENT

Predicated upon 3D cubic B-splines interpolation of the image velocity field, WSS was estimated in mid-systole, early-diastole, and late-diastole and OSI was derived. We assessed the impact of spatiotemporal resolution and phase noise, and compared results based on tracked-using deformable registration-and static vessel wall location.

STATISTICAL TESTS

Bland-Altman analysis to assess WSS/OSI differences; Hausdorff distance (HD) to assess wall motion; and Pearson's correlation coefficient (PCC) to assess correlation of HD with WSS.

RESULTS

Synthetic data results show systematic over-/underestimation of WSS when different spatial resolution (mean ± 1.96 SD up to -0.24 ± 0.40 N/m2 and 0.5 ± 1.38 N/m2 for 8-fold and 27-fold voxel size, respectively) and VENC-depending phase noise (mean ± 1.96 SD up to 0.31 ± 0.12 N/m2 and 0.94 ± 0.28 N/m2 for 2-fold and 4-fold VENC increase, respectively) are given. Neglecting wall motion when defining the vessel wall perturbs WSS estimates to a considerable extent (1.96 SD up to 1.21 N/m2 ) without systematic over-/underestimation (Bland-Altman mean range -0.06 to 0.05).

DATA CONCLUSION

In addition to sufficient spatial resolution and velocity to noise ratio, accurate tracking of the vessel wall is essential for reliable image-based WSS estimation and should not be neglected if wall motion is present.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2018.

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.

Experimental quantification of the fluid dynamics in blood-processing devices through 4D-flow imaging: A pilot study on a real oxygenator/heat-exchanger module

The performance of blood-processing devices largely depends on the associated fluid dynamics, which hence represents a key aspect in their design and optimization. To this aim, two approaches are currently adopted: computational fluid-dynamics, which yields highly resolved three-dimensional data but relies on simplifying assumptions, and in vitro experiments, which typically involve the direct video-acquisition of the flow field and provide 2D data only. We propose a novel method that exploits space- and time-resolved magnetic resonance imaging (4D-flow) to quantify the complex 3D flow field in blood-processing devices and to overcome these limitations. We tested our method on a real device that integrates an oxygenator and a heat exchanger. A dedicated mock loop was implemented, and novel 4D-flow sequences with sub-millimetric spatial resolution and region-dependent velocity encodings were defined. Automated in house software was developed to quantify the complex 3D flow field within the different regions of the device: region-dependent flow rates, pressure drops, paths of the working fluid and wall shear stresses were computed. Our analysis highlighted the effects of fine geometrical features of the device on the local fluid-dynamics, which would be unlikely observed by current in vitro approaches. Also, the effects of non-idealities on the flow field distribution were captured, thanks to the absence of the simplifying assumptions that typically characterize numerical models. To the best of our knowledge, our approach is the first of its kind and could be extended to the analysis of a broad range of clinically relevant devices.

4D Flow Analysis of BAV-Related Fluid-Dynamic Alterations: Evidences of Wall Shear Stress Alterations in Absence of Clinically-Relevant Aortic Anatomical Remodeling

Bicuspid aortic valve (BAV) is the most common congenital cardiac disease and is a foremost risk factor for aortopathies. Despite the genetic basis of BAV and of the associated aortopathies, BAV-related alterations in aortic fluid-dynamics, and particularly in wall shear stresses (WSSs), likely play a role in the progression of aortopathy, and may contribute to its pathogenesis. To test whether WSS may trigger aortopathy, in this study we used 4D Flow sequences of phase-contrast cardiac magnetic resonance imaging (CMR) to quantitatively compare the in vivo fluid dynamics in the thoracic aorta of two groups of subjects: (i) five prospectively enrolled young patients with normo-functional BAV and with no aortic dilation and (ii) ten age-matched healthy volunteers. Through the semi-automated processing of 4D Flow data, the aortic bulk flow at peak systole was quantified, and WSSs acting on the endothelium of the ascending aorta were characterized throughout the systolic phase in terms of magnitude and time-dependency through a method recently developed by our group. Variables computed for each BAV patient were compared vs. the corresponding distribution of values obtained for healthy controls. In BAV patients, ascending aorta diameter was measured on cine-CMR images at baseline and at 3-year follow-up. As compared to controls, normo-functional BAV patients were characterized by minor bulk flow disturbances at peak systole. However, they were characterized by evident alterations of WSS distribution and peak values in the ascending aorta. In particular, in four BAV patients, who were characterized by right-left leaflet fusion, WSS peak values exceeded by 27–46% the 90th percentile of the distribution obtained for healthy volunteers. Only in the BAV patient with right-non-coronary leaflet fusion the same threshold was exceeded by 132%. Also, evident alterations in the time-dependency of WSS magnitude and direction were observed. Despite, these fluid-dynamic alterations, no clinically relevant anatomical remodeling was observed in the BAV patients at 3-year follow-up. In light of previous evidence from the literature, our results suggest that WSS alterations may precede the onset of aortopathy and may contribute to its triggering, but WSS-driven anatomical remodeling, if any, is a very slow process.

Association between flow skewness and aortic dilatation in patients with aortic stenosis

We investigated association between hemodynamic characteristics and aortic dilatation in patients with severe aortic stenosis (AS). Eighty patients with severe AS (mean age, 67.2 ± 12.5 years) who underwent multi-detector computed tomography and phase-contrast magnetic resonance imaging at the ascending aorta were retrospectively analyzed. Patients with an ascending aorta diameter >4 cm had a significantly higher forward flow rate at systole (28.5 ± 6.0 vs. 36.2 ± 8.6 L min, P < 0.001), and retrograde flow rate at systole (11.3 ± 4.2 vs. 18.8 ± 5.8 L min, P < 0.001), fractional reverse ratio (a ratio of retrograde flow rate to forward flow rate; 34.1 ± 11.9% vs. 43.5 ± 18.0%, P = 0.014), flow skewness R_skewness (a ratio of sum of forward and retrograde systole flow to net systole flow rate; 2.4 ± 0.7 vs. 3.2 ± 1.0, P < 0.001). The presence of bicuspid aortic valve (BAV; odds ratio [OR] 72.01, 95% confidence interval [CI] 10.57-490.46, P < 0.001), Left ventricular mass index (LVMI; OR 1.02 /g/m²; CI 1.00-1.04, P = 0.043) and R_skewness (OR 5.6 per 1, 95% CI 1.8-17.1, P = 0.001) were associated with aortic dilatation. BAV, LVMI, and increased R_skewness in the ascending aorta are associated with aortic dilatation in patients with AS.

Morphotype-Dependent Flow Characteristics in Bicuspid Aortic Valve Ascending Aortas: A Benchtop Particle Image Velocimetry Study

The bicuspid aortic valve (BAV) is a major risk factor for secondary aortopathy such as aortic dilation. The heterogeneous BAV morphotypes [left-right-coronary cusp fusion (LR), right-non-coronary cusp fusion (RN), and left-non-coronary cusp fusion (LN)] are associated with different dilation patterns, suggesting a role for hemodynamics in BAV aortopathogenesis. However, assessment of this theory is still hampered by the limited knowledge of the hemodynamic abnormalities generated by the distinct BAV morphotypes. The objective of this study was to compare experimentally the hemodynamics of a normal (i.e., non-dilated) ascending aorta (AA) subjected to tricuspid aortic valve (TAV), LR-BAV, RN-BAV, and NL-BAV flow. Tissue BAVs reconstructed from porcine TAVs were subjected to physiologic pulsatile flow conditions in a left-heart simulator featuring a realistic aortic root and compliant aorta. Phase-locked particle image velocimetry experiments were carried out to characterize the flow in the aortic root and in the tubular AA in terms of jet skewness and displacement, as well as mean velocity, viscous shear stress and Reynolds shear stress fields. While all three BAVs generated skewed and asymmetrical orifice jets (up to 1.7- and 4.0-fold increase in flow angle and displacement, respectively, relative to the TAV at the sinotubular junction), the RN-BAV jet was out of the plane of observation. The LR- and NL-BAV exhibited a 71% increase in peak-systolic orifice jet velocity relative to the TAV, suggesting an inherent degree of stenosis in BAVs. While these two BAV morphotypes subjected the convexity of the aortic wall to viscous shear stress overloads (1.7-fold increase in maximum peak-systolic viscous shear stress relative to the TAV-AA), the affected sites were morphotype-dependent (LR-BAV: proximal AA, NL-BAV: distal AA). Lastly, the LR- and NL-BAV generated high degrees of turbulence in the AA (up to 2.3-fold increase in peak-systolic Reynolds shear stress relative to the TAV) that were sustained from peak systole throughout the deceleration phase. This in vitro study reveals substantial flow abnormalities (increased jet skewness, asymmetry, jet velocity, turbulence, and shear stress overloads) in non-dilated BAV aortas, which differ from those observed in dilated aortas but still coincide with aortic wall regions prone to dilation.