Advanced flow MRI: emerging techniques and applications. Clin Radiol. 2016;71:779-795. [PMID: 26944696 DOI: 10.1016/j.crad.2016.01.011]

Effect of MRI acquisition parameters on accuracy and precision of phase-contrast measurements in a small-lumen vessel phantom

BACKGROUND

Phase-contrast magnetic resonance imaging (PC-MRI) quantifies blood flow and velocity noninvasively. Challenges arise in neurovascular disorders due to small vessels. We evaluated the impact of voxel size, number of signal averages (NSA), and velocity encoding (VENC) on PC-MRI measurement accuracy and precision in a small-lumen vessel phantom.

METHODS

We constructed an in vitro model with a constant flow rate using a 2.2-mm inner diameter plastic tube. A reservoir with a weight scale and timer was used as standard reference. Gradient-echo T1 weighted PC-MRI sequence was performed on a 3-T scanner with varying voxel size (2.5, 5.0, 7.5 mm3), NSA (1, 2, 3), and VENC (200, 300, 400 cm/s). We repeated measurements nine times per setting, calculating mean flow rate, maximum velocity, and least detectable difference (LDD).

RESULTS

PC-MRI flow measurements were higher than standard reference values (mean ranging from 7.3 to 9.5 mL/s compared with 6.6 mL/s). Decreased voxel size improved accuracy, reducing flow rate measurements from 9.5 to 7.3 mL/s. The LDD for flow rate and velocity varied between 1 and 5%. The LDD for flow rate decreased with increased voxel size and NSA (p = 0.033 and 0.042). The LDD for velocity decreased with increased voxel size (p < 10-16). No change was observed when VENC varied.

CONCLUSIONS

PC-MRI overestimated flow. However, it has high precision in a small-vessel phantom with constant flow rate. Improved accuracy was obtained with increasing spatial resolution (smaller voxels). Improved precision was obtained with increasing signal-to-noise ratio (larger voxels and/or higher NSA).

RELEVANCE STATEMENT

Phase-contrast MRI is clinically used in large vessels. To further investigate the possibility of using phase-contrast MRI for smaller intracranial vessels in neurovascular disorders, we need to understand how acquisition parameters affect phase-contrast MRI-measured flow rate and velocity in small vessels.

KEY POINTS

• PC-MRI measures flow and velocity in a small lumen phantom with high precision but overestimates flow rate. • The precision of PC-MRI measurements matches the precision of standard reference for flow rate measurements. • Optimizing PC-MRI settings can enhance accuracy and precision in flow rate and velocity measurements.

Patient-Specific Haemodynamic Analysis of Virtual Grafting Strategies in Type-B Aortic Dissection: Impact of Compliance Mismatch

INTRODUCTION

Compliance mismatch between the aortic wall and Dacron Grafts is a clinical problem concerning aortic haemodynamics and morphological degeneration. The aortic stiffness introduced by grafts can lead to an increased left ventricular (LV) afterload. This study quantifies the impact of compliance mismatch by virtually testing different Type-B aortic dissection (TBAD) surgical grafting strategies in patient-specific, compliant computational fluid dynamics (CFD) simulations.

MATERIALS AND METHODS

A post-operative case of TBAD was segmented from computed tomography angiography data. Three virtual surgeries were generated using different grafts; two additional cases with compliant grafts were assessed. Compliant CFD simulations were performed using a patient-specific inlet flow rate and three-element Windkessel outlet boundary conditions informed by 2D-Flow MRI data. The wall compliance was calibrated using Cine-MRI images. Pressure, wall shear stress (WSS) indices and energy loss (EL) were computed.

RESULTS

Increased aortic stiffness and longer grafts increased aortic pressure and EL. Implementing a compliant graft matching the aortic compliance of the patient reduced the pulse pressure by 11% and EL by 4%. The endothelial cell activation potential (ECAP) differed the most within the aneurysm, where the maximum percentage difference between the reference case and the mid (MDA) and complete (CDA) descending aorta replacements increased by 16% and 20%, respectively.

CONCLUSION

This study suggests that by minimising graft length and matching its compliance to the native aorta whilst aligning with surgical requirements, the risk of LV hypertrophy may be reduced. This provides evidence that compliance-matching grafts may enhance patient outcomes.

Evaluating a Phase-Specific Approach to Aortic Flow: A 4D Flow MRI Study

BACKGROUND

Aortic flow parameters can be quantified using 4D flow MRI. However, data are sparse on how different methods of analysis influence these parameters and how these parameters evolve during systole.

PURPOSE

To assess multiphase segmentations and multiphase quantification of flow-related parameters in aortic 4D flow MRI.

STUDY TYPE

Prospective.

POPULATION

40 healthy volunteers (50% male, 28.9 ± 5.0 years) and 10 patients with thoracic aortic aneurysm (80% male, 54 ± 8 years).

FIELD STRENGTH/SEQUENCE

4D flow MRI with a velocity encoded turbo field echo sequence at 3 T.

ASSESSMENT

Phase-specific segmentations were obtained for the aortic root and the ascending aorta. The whole aorta was segmented in peak systole. In all aortic segments, time to peak (TTP; for flow velocity, vorticity, helicity, kinetic energy, and viscous energy loss) and peak and time-averaged values (for velocity and vorticity) were calculated.

STATISTICAL TESTS

Static vs. phase-specific models were assessed using Bland-Altman plots. Other analyses were performed using phase-specific segmentations for aortic root and ascending aorta. TTP for all parameters was compared to TTP of flow rate using paired t-tests. Time-averaged and peak values were assessed using Pearson correlation coefficient. P < 0.05 was considered statistically significant.

RESULTS

In the combined group, velocity in static vs. phase-specific segmentations differed by 0.8 cm/sec for the aortic root, and 0.1 cm/sec (P = 0.214) for the ascending aorta. Vorticity differed by 167 sec-1  mL-1 (P = 0.468) for the aortic root, and by 59 sec-1  mL-1 (P = 0.481) for the ascending aorta. Vorticity, helicity, and energy loss in the ascending aorta, aortic arch, and descending aorta peaked significantly later than flow rate. Time-averaged velocity and vorticity values correlated significantly in all segments.

DATA CONCLUSION

Static 4D flow MRI segmentation yields comparable results as multiphase segmentation for flow-related parameters, eliminating the need for time-consuming multiple segmentations. However, multiphase quantification is necessary for assessing peak values of aortic flow-related parameters.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY STAGE: 3.

A characterization of cardiac-induced noise in R2 * maps of the brain

PURPOSE

Cardiac pulsation increases the noise level in brain maps of the transverse relaxation rate R2 *. Cardiac-induced noise is challenging to mitigate during the acquisition of R2 * mapping data because its characteristics are unknown. In this work, we aim to characterize cardiac-induced noise in brain maps of the MRI parameter R2 *.

METHODS

We designed a sampling strategy to acquire multi-echo 3D data in 12 intervals of the cardiac cycle, monitored with a fingertip pulse-oximeter. We measured the amplitude of cardiac-induced noise in this data and assessed the effect of cardiac pulsation on R2 * maps computed across echoes. The area of k-space that contains most of the cardiac-induced noise in R2 * maps was then identified. Based on these characteristics, we introduced a tentative sampling strategy that aims to mitigate cardiac-induced noise in R2 * maps of the brain.

RESULTS

In inferior brain regions, cardiac pulsation accounts for R2 * variations of up to 3 s-1 across the cardiac cycle (i.e., ∼35% of the overall variability). Cardiac-induced fluctuations occur throughout the cardiac cycle, with a reduced intensity during the first quarter of the cycle. A total of 50% to 60% of the overall cardiac-induced noise is localized near the k-space center (k < 0.074 mm-1 ). The tentative cardiac noise mitigation strategy reduced the variability of R2 * maps across repetitions by 11% in the brainstem and 6% across the whole brain.

CONCLUSION

We provide a characterization of cardiac-induced noise in brain R2 * maps that can be used as a basis for the design of mitigation strategies during data acquisition.

Deep learning methods for blood flow reconstruction in a vessel with contrast enhanced x-ray computed tomography

The reconstruction of blood velocity in a vessel from contrast enhanced x-ray computed tomography projections is a complex inverse problem. It can be formulated as reconstruction problem with a partial differential equation constraint. A solution can be estimated with the a variational adjoint method and proper orthogonal decomposition (POD) basis. In this work, we investigate new inversion approaches based on PODs coupled with deep learning methods. The effectiveness of the reconstruction methods is shown with simulated realistic stationary blood flows in a vessel. The methods outperform the reduced adjoint method and show large speed-up at the online stage.

Modelling blood flow in coronary arteries: Newtonian or shear-thinning non-Newtonian rheology?

BACKGROUND

The combination of medical imaging and computational hemodynamics is a promising technology to diagnose/prognose coronary artery disease (CAD). However, the clinical translation of in silico hemodynamic models is still hampered by assumptions/idealizations that must be introduced in model-based strategies and that necessarily imply uncertainty. This study aims to provide a definite answer to the open question of how to properly model blood rheological properties in computational fluid dynamics (CFD) simulations of coronary hemodynamics.

METHODS

The geometry of the right coronary artery (RCA) of 144 hemodynamically stable patients with different stenosis degree were reconstructed from angiography. On them, unsteady-state CFD simulations were carried out. On each reconstructed RCA two different simulation strategies were applied to account for blood rheological properties, implementing (i) a Newtonian (N) and (ii) a shear-thinning non-Newtonian (non-N) rheological model. Their impact was evaluated in terms of wall shear stress (WSS magnitude, multidirectionality, topological skeleton) and helical flow (strength, topology) profiles. Additionally, luminal surface areas (SAs) exposed to shear disturbances were identified and the co-localization of paired N and non-N SAs was quantified in terms of similarity index (SI).

RESULTS

The comparison between paired N vs. shear-thinning non-N simulations revealed remarkably similar profiles of WSS-based and helicity-based quantities, independent of the adopted blood rheology model and of the degree of stenosis of the vessel. Statistically, for each paired N and non-N hemodynamic quantity emerged negligible bias from Bland-Altman plots, and strong positive linear correlation (r > 0.94 for almost all the WSS-based quantities, r > 0.99 for helicity-based quantities). Moreover, a remarkable co-localization of N vs. non-N luminal SAs exposed to disturbed shear clearly emerged (SI distribution 0.95 [0.93, 0.97]). Helical flow topology resulted to be unaffected by blood rheological properties.

CONCLUSIONS

This study, performed on 288 angio-based CFD simulations on 144 RCA models presenting with different degrees of stenosis, suggests that the assumptions on blood rheology have negligible impact both on WSS and helical flow profiles associated with CAD, thus definitively answering to the question "is Newtonian assumption for blood rheology adequate in coronary hemodynamics simulations?".

Valvular and ascending aortic hemodynamics of the On-X aortic valved conduit by same-day echocardiography and 4D flow MRI

This study aims to assess whether the On-X aortic valved conduit better restores normal valvular and ascending aortic hemodynamics than other commonly used bileaflet mechanical valved conduit prostheses from St. Jude Medical and Carbomedics by using same-day transthoracic echocardiography (TTE) and 4D flow magnetic resonance imaging (MRI) examinations. TTE and 4D flow MRI were performed back-to-back in 10 patients with On-X, six patients with St. Jude (two) and Carbomedics (four) prostheses, and 36 healthy volunteers. TTE evaluated valvular hemodynamic parameters: transvalvular peak velocity (TPV), mean and peak transvalvular pressure gradient (TPG), and effective orifice area (EOA). 4D flow MRI evaluated the peak systolic 3D viscous energy loss rate (VELR) density and mean vorticity magnitude in the ascending aorta (AAo). While higher TPV and mean and peak TPG were recorded in all patients compared to healthy subjects, the values in On-X patients were closer to those in healthy subjects (TPV 1.9 ± 0.3 vs. 2.2 ± 0.3 vs. 1.2 ± 0.2 m/s, mean TPG 7.4 ± 1.9 vs. 9.2 ± 2.3 vs. 3.1 ± 0.9 mmHg, peak TPG 15.3 ± 5.2 vs. 18.9 ± 5.2 vs. 6.1 ± 1.8 mmHg, p < 0.001). Likewise, while higher VELR density and mean vorticity magnitude were recorded in all patients than in healthy subjects, the values in On-X patients were closer to those in healthy subjects (VELR: 50.6 ± 20.1 vs. 89.8 ± 35.2 vs. 21.4 ± 9.2 W/m3, p < 0.001) and vorticity (147.6 ± 30.0 vs. 191.2 ± 26.0 vs. 84.6 ± 20.5 s-1, p < 0.001). This study demonstrates that the On-X aortic valved conduit may produce less aberrant hemodynamics in the AAo while maintaining similar valvular hemodynamics to St. Jude Medical and Carbomedics alternatives.

Energetics of Cardiac Blood Flow in Hypertrophic Cardiomyopathy through Individualized Computational Modeling

Hypertrophic cardiomyopathy (HCM) is a congenital heart disease characterized by thickening of the heart's left ventricle (LV) wall that can lead to cardiac dysfunction and heart failure. Ventricular wall thickening affects the motion of cardiac walls and blood flow within the heart. Because abnormal cardiac blood flow in turn could lead to detrimental remodeling of heart walls, aberrant ventricular flow patterns could exacerbate HCM progression. How blood flow patterns are affected by hypertrophy and inter-patient variability is not known. To address this gap in knowledge, we present here strategies to generate personalized computational fluid dynamics (CFD) models of the heart LV from patient cardiac magnetic resonance (cMR) images. We performed simulations of CFD LV models from three cases (one normal, two HCM). CFD computations solved for blood flow velocities, from which flow patterns and the energetics of flow within the LV were quantified. We found that, compared to a normal heart, HCM hearts exhibit anomalous flow patterns and a mismatch in the timing of energy transfer from the LV wall to blood flow, as well as changes in kinetic energy flow patterns. While our results are preliminary, our presented methodology holds promise for in-depth analysis of HCM patient hemodynamics in clinical practice.

4D Flow cardiovascular magnetic resonance consensus statement: 2023 update

Hemodynamic assessment is an integral part of the diagnosis and management of cardiovascular disease. Four-dimensional cardiovascular magnetic resonance flow imaging (4D Flow CMR) allows comprehensive and accurate assessment of flow in a single acquisition. This consensus paper is an update from the 2015 '4D Flow CMR Consensus Statement'. We elaborate on 4D Flow CMR sequence options and imaging considerations. The document aims to assist centers starting out with 4D Flow CMR of the heart and great vessels with advice on acquisition parameters, post-processing workflows and integration into clinical practice. Furthermore, we define minimum quality assurance and validation standards for clinical centers. We also address the challenges faced in quality assurance and validation in the research setting. We also include a checklist for recommended publication standards, specifically for 4D Flow CMR. Finally, we discuss the current limitations and the future of 4D Flow CMR. This updated consensus paper will further facilitate widespread adoption of 4D Flow CMR in the clinical workflow across the globe and aid consistently high-quality publication standards.

A clinician's guide to understanding aortic 4D flow MRI

Four-dimensional flow magnetic resonance imaging is an emerging technique which may play a role in diagnosis and risk-stratification of aortic disease. Some knowledge of flow dynamics and related parameters is necessary to understand and apply this technique in clinical workflows. The purpose of the current review is to provide a guide for clinicians to the basics of flow imaging, frequently used flow-related parameters, and their relevance in the context of aortic disease.Clinical relevance statement Understanding normal and abnormal aortic flow could improve clinical care in patients with aortic disease.

Near-wall hemodynamic parameters quantification in in vitro intracranial aneurysms with 7 T PC-MRI

OBJECTIVE

Wall shear stress (WSS) and its derived spatiotemporal parameters have proven to play a major role on intracranial aneurysms (IAs) growth and rupture. This study aims to demonstrate how ultra-high field (UHF) 7 T phase contrast magnetic resonance imaging (PC-MRI) coupled with advanced image acceleration techniques allows a highly resolved visualization of near-wall hemodynamic parameters patterns in in vitro IAs, paving the way for more robust risk assessment of their growth and rupture.

MATERIALS AND METHODS

We performed pulsatile flow measurements inside three in vitro models of patient-specific IAs using 7 T PC-MRI. To this end, we built an MRI-compatible test bench, which faithfully reproduced a typical physiological intracranial flow rate in the models.

RESULTS

The ultra-high field 7 T images revealed WSS patterns with high spatiotemporal resolution. Interestingly, the high oscillatory shear index values were found in the core of low WSS vortical structures and in flow stream intersecting regions. In contrast, maxima of WSS occurred around the impinging jet sites.

CONCLUSIONS

We showed that the elevated signal-to-noise ratio arising from 7 T PC-MRI enabled to resolve high and low WSS patterns with a high degree of detail.

Wall Shear Stress Differences Between Arterial and Venous Coronary Artery Bypass Grafts One Month After Surgery

Although coronary artery bypass graft (CABG) surgery is a well-established intervention, graft failure can occur, and the underlying mechanisms remain incompletely understood. The purpose of this prospective study is to utilize computational fluid dynamics (CFD) to investigate how graft hemodynamics one month post surgery may vary among graft types, which have different long-term patency rates. Twenty-four grafts from 10 participants (64.6 ± 8.5 years, 9 men) were scanned with coronary CT angiography and 4D flow MRI one month after CABG surgery. Grafts included 10 left internal mammary arteries (LIMA), 3 radial arteries (RA), and 11 saphenous vein grafts (SVG). Image-guided CFD was used to quantify blood flow rate and wall area exposed to abnormal wall shear stress (WSS). Arterial grafts had a lower abnormal WSS area than venous grafts (17.9% vs. 70.1%; p = 0.001), and a similar trend was observed for LIMA vs. SVG (13.8% vs. 70.1%; p = 0.001). Abnormal WSS area correlated positively to lumen diameter (p < 0.001) and negatively to flow rate (p = 0.001). This CFD study is the first of its kind to prospectively reveal differences in abnormal WSS area 1 month post surgery among CABG types, suggesting that WSS may influence the differential long-term graft failure rates observed among these groups.

Highly Accelerated Compressed-Sensing 4D Flow for Intracardiac Flow Assessment

BACKGROUND

Four-dimensional (4D) flow MRI allows for the quantification of complex flow patterns; however, its clinical use is limited by its inherently long acquisition time. Compressed sensing (CS) is an acceleration technique that provides substantial reduction in acquisition time.

PURPOSE

To compare intracardiac flow measurements between conventional and CS-based highly accelerated 4D flow acquisitions.

STUDY TYPE

Prospective.

SUBJECTS

Fifty healthy volunteers (28.0 ± 7.1 years, 24 males).

FIELD STRENGTH/SEQUENCE

Whole heart time-resolved 3D gradient echo with three-directional velocity encoding (4D flow) with conventional parallel imaging (factor 3) as well as CS (factor 7.7) acceleration at 3 T.

ASSESSMENT

4D flow MRI data were postprocessed by applying a valve tracking algorithm. Acquisition times, flow volumes (mL/cycle) and diastolic function parameters (ratio of early to late diastolic left ventricular peak velocities [E/A] and ratio of early mitral inflow velocity to mitral annular early diastolic velocity [E/e']) were quantified by two readers.

STATISTICAL TESTS

Paired-samples t-test and Wilcoxon rank sum test to compare measurements. Pearson correlation coefficient (r), Bland-Altman-analysis (BA) and intraclass correlation coefficient (ICC) to evaluate agreement between techniques and readers. A P value < 0.05 was considered statistically significant.

RESULTS

A significant improvement in acquisition time was observed using CS vs. conventional accelerated acquisition (6.7 ± 1.3 vs. 12.0 ± 1.3 min). Net forward flow measurements for all valves showed good correlation (r > 0.81) and agreement (ICCs > 0.89) between conventional and CS acceleration, with 3.3%-8.3% underestimation by the CS technique. Evaluation of diastolic function showed 3.2%-17.6% error: E/A 2.2 [1.9-2.4] (conventional) vs. 2.3 [2.0-2.6] (CS), BA bias 0.08 [-0.81-0.96], ICC 0.82; and E/e' 4.6 [3.9-5.4] (conventional) vs. 3.8 [3.4-4.3] (CS), BA bias -0.90 [-2.31-0.50], ICC 0.89.

DATA CONCLUSION

Analysis of intracardiac flow patterns and evaluation of diastolic function using a highly accelerated 4D flow sequence prototype is feasible, but it shows underestimation of flow measurements by approximately 10%.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 1.

Full-volume three-component intraventricular vector flow mapping by triplane color Doppler

Objective. Intraventricular vector flow mapping (iVFM) is a velocimetric technique for retrieving two-dimensional velocity vector fields of blood flow in the left ventricular cavity. This method is based on conventional color Doppler imaging, which makes iVFM compatible with the clinical setting. We have generalized the iVFM for a three-dimensional reconstruction (3D-iVFM). Approach. 3D-iVFM is able to recover three-component velocity vector fields in a full intraventricular volume by using a clinical echocardiographic triplane mode. The 3D-iVFM problem was written in the spherical (radial, polar, azimuthal) coordinate system associated to the six half-planes produced by the triplane mode. As with the 2D version, the method is based on the mass conservation, and free-slip boundary conditions on the endocardial wall. These mechanical constraints were imposed in a least-squares minimization problem that was solved through the method of Lagrange multipliers. We validated 3D-iVFM in silico in a patient-specific CFD (computational fluid dynamics) model of cardiac flow and tested its clinical feasibility in vivo in patients and in one volunteer. Main results. The radial and polar components of the velocity were recovered satisfactorily in the CFD setup (correlation coefficients, r  = 0.99 and 0.78). The azimuthal components were estimated with larger errors ( r  = 0.57) as only six samples were available in this direction. In both in silico and in vivo investigations, the dynamics of the intraventricular vortex that forms during diastole was deciphered by 3D-iVFM. In particular, the CFD results showed that the mean vorticity can be estimated accurately by 3D-iVFM. Significance. Our results tend to indicate that 3D-iVFM could provide full-volume echocardiographic information on left intraventricular hemodynamics from the clinical modality of triplane color Doppler.

Phase contrast cardiac magnetic resonance imaging versus transoesophageal echocardiography for the evaluation of feasibility for transcatheter closure of atrial septal defects

Background

This prospective study was aimed at comparing phase contrast cardiac magnetic resonance imaging (PC-CMR) with 2D transoesophageal echocardiography (TEE) for determining potential candidature for transcatheter closure in ostium secundum ASD (OS-ASD) patients. We included consecutive adult patients with OS-ASD for the evaluation of feasibility for transcatheter closure using 2D-TEE and PC-CMR over a period of 2 years. Patients who fulfilled the conventional criteria for transcatheter closure, i.e. maximum ASD diameter ≤ 34 mm, adequate rims (≥ 5 mm, except for anterosuperior rim), and normal pulmonary venous drainage on both imaging modalities, were taken for device closure. In patients where there was discrepancy in the measurements of ASD diameter or rim size, making them eligible for device closure on one imaging modality and ineligible on the other hand, provisional device closure was attempted. All patients who underwent transcatheter closure were followed up to 6 months to rule out any complications.

Results

A total of 58 patients (mean age 35.93 ± 10.59 years) were enrolled in the study. Overall, there was significant positive correlation between 2D-TEE and CMR measurements of maximal ASD diameter and rim size (p < 0.001). However, TEE significantly underestimated maximal ASD diameter and posteroinferior rim size in comparison with CMR (p = 0.013 and p = 0.023, respectively). 46 (79.3%) patients were suitable for transcatheter closure on CMR, while 44 (75.9%) were eligible on TEE. Transcatheter closure was attempted in 48 patients based on imaging findings and was successful in 46 (95.8%) patients. Device closure was unsuccessful in 2 patients with defect size < 34 mm on TEE but > 34 mm on CMR. Among 7 patients with deficient posteroinferior rim on TEE, 5 had sufficient rim on CMR and underwent successful transcatheter closure. CMR detected anomalous pulmonary venous drainage in one patient which was missed on TEE, hence excluding the patient from transcatheter closure. Mean device size was 28.3 ± 7.4 mm and correlated more strongly with the defect dimensions on PC-CMR (r = 0.85, p < 0.001) compared to TEE (r = 0.71, p = 0.02).

Conclusions

PC-CMR may to be superior to 2D-TEE for the preprocedural planning and feasibility assessment for transcatheter closure in adult patients with ostium secundum ASD.

Rotational Vertebrobasilar Insufficiency: Is There a Physiological Spectrum? Phase-Contrast Magnetic Resonance Imaging Quantification in Healthy Volunteers

BACKGROUND

Some cases of cerebral ischemia have been attributed to dynamic flow limitation in neck vessels. It however remains unknown whether this represents the extreme end of a physiological response.

METHODS

Eighteen healthy volunteers were recruited to this prospective study. Cervical blood flow (ml/min/m²) was assessed using phase-contrast MRI, and cerebral perfusion ratios were assessed using arterial spin labeling perfusion at neutral position, predefined head rotations, as well as flexion and extension. Inter-reader agreements were assessed using intraclass correlation coefficient.

RESULTS

The mean age was 38.6 ± 10.8 (range = 22-56) years, for five male participants and 13 females. The means for height and weight were 168 cm and 73.2 kg, respectively. There were no significant differences in individual arterial blood flow with change in head position (P > 0.05). Similarly, the repeated-measures analysis of variance test demonstrated no significant difference in perfusion ratios in relation to head position movement (P > 0.05). Inter-reader agreement was excellent (intraclass correlation coefficient = 0.97).

CONCLUSIONS

There is neither significant change in either individual cervical arterial blood flow nor cerebral perfusion within the normal physiological/anatomical range of motion in healthy individuals. It is therefore reasonable to conclude that any such hemodynamic change identified in a patient with ischemic stroke be considered causative.

Measuring Arterial Pulsatility With Dynamic Inflow Magnitude Contrast

The pulsatility of blood flow through cerebral arteries is clinically important, as it is intrinsically associated with cerebrovascular health. In this study we outline a new MRI approach to measuring the real-time pulsatile flow in cerebral arteries, which is based on the inflow phenomenon associated with fast gradient-recalled-echo acquisitions. Unlike traditional phase-contrast techniques, this new method, which we dub dynamic inflow magnitude contrast (DIMAC), does not require velocity-encoding gradients as sensitivity to flow velocity is derived purely from the inflow effect. We achieved this using a highly accelerated single slice EPI acquisition with a very short TR (15 ms) and a 90° flip angle, thus maximizing inflow contrast. We simulate the spoiled GRE signal in the presence of large arteries and perform a sensitivity analysis. The sensitivity analysis demonstrates that in the regime of high inflow contrast, DIMAC shows much greater sensitivity to flow velocity over blood volume changes. We support this theoretical prediction with in-vivo data collected in two separate experiments designed to demonstrate the utility of the DIMAC signal contrast. We perform a hypercapnia challenge experiment in order to experimentally modulate arterial tone within subjects, and thus modulate the arterial pulsatile flow waveform. We also perform a thigh-cuff release challenge, designed to induce a transient drop in blood pressure, and demonstrate that the continuous DIMAC signal captures the complex transient change in the pulsatile and non-pulsatile components of flow. In summary, this study proposes a new role for a well-established source of MR image contrast and demonstrates its potential for measuring both steady-state and dynamic changes in arterial tone.

A multi-modality approach for enhancing 4D flow magnetic resonance imaging via sparse representation

This work evaluates and applies a multi-modality approach to enhance blood flow measurements and haemodynamic analysis with phase-contrast magnetic resonance imaging (4D flow MRI) in cerebral aneurysms (CAs). Using a library of high-resolution velocity fields from patient-specific computational fluid dynamic simulations and in vitro particle tracking velocimetry measurements, the flow field of 4D flow MRI data is reconstructed as the sparse representation of the library. The method was evaluated with synthetic 4D flow MRI data in two CAs. The reconstruction enhanced the spatial resolution and velocity accuracy of the synthetic MRI data, leading to reliable pressure and wall shear stress (WSS) evaluation. The method was applied on in vivo 4D flow MRI data acquired in the same CAs. The reconstruction increased the velocity and WSS by 6-13% and 39-61%, respectively, suggesting that the accuracy of these quantities was improved since the raw MRI data underestimated the velocity and WSS by 10-20% and 40-50%, respectively. The computed pressure fields from the reconstructed data were consistent with the observed flow structures. The results suggest that using the sparse representation flow reconstruction with in vivo 4D flow MRI enhances blood flow measurement and haemodynamic analysis.

Divergence-Free Constrained Phase Unwrapping and Denoising for 4D Flow MRI Using Weighted Least-Squares

A novel divergence-free constrained phase unwrapping method was proposed and evaluated for 4D flow MRI. The unwrapped phase field was obtained by integrating the phase variations estimated from the wrapped phase data using weighted least-squares. The divergence-free constraint for incompressible blood flow was incorporated to regulate and denoise the resulting phase field. The proposed method was tested on synthetic phase data of left ventricular flow and in vitro 4D flow measurement of Poiseuille flow. The method was additionally applied to in vivo 4D flow measurements in the thoracic aorta from 30 human subjects. The performance of the proposed method was compared to the state-of-the-art 4D single-step Laplacian algorithm. The synthetic phase data were completely unwrapped by the proposed method for all the cases with velocity encoding (venc) as low as 20% of the maximum velocity and signal-to-noise ratio as low as 5. The in vitro Poiseuille flow data were completely unwrapped with a 60% increase in the velocity-to-noise ratio. For the in-vivo aortic datasets with venc ratio less than 0.4, the proposed method significantly improved the success rate by as much as 40% and reduced the velocity error levels by a factor of 10 compared to the state-of-the-art method. The divergence-free constrained method exhibits reliability and robustness on phase unwrapping and shows improved accuracy of velocity and hemodynamic quantities by unwrapping the low-venc 4D flow MRI data.

Cross-modality and in-vivo validation of 4D flow MRI evaluation of uterine artery blood flow in human pregnancy

OBJECTIVES

Clinical assessment of uterine artery (UtA) hemodynamics is currently limited to Doppler ultrasound (US) velocimetry. We have demonstrated previously the feasibility of applying four-dimensional (4D) flow magnetic resonance imaging (MRI) to evaluate UtA hemodynamics during pregnancy, allowing flow quantification of the entire course of the vessel. In this study, we sought to further validate the physiological relevance of 4D flow MRI measurement of UtA blood flow by exploring its association with pregnancy outcome relative to US-based metrics.

METHODS

Recruited into this prospective, cross-sectional study were 87 women with a singleton pregnancy who underwent 4D flow MRI between May 2016 and April 2019 to measure the UtA pulsatility index (MRI-PI) and blood flow rate (MRI-flow, in mL/min). UtA-PI was also measured using US (US-PI). The primary outcome was a composite (COMP) of pre-eclampsia (PE) and/or small-for-gestational-age (SGA) neonate, and secondary outcomes were PE and SGA neonate individually. We assessed the ability of MRI-flow, MRI-PI and US-PI to distinguish between outcomes, and evaluated whether MRI-flow changed as gestation progressed.

RESULTS

Following 4D flow postprocessing and exclusions from the analysis, 74 women had 4D flow MRI data analyzed for both UtAs. Of these, 18 developed a COMP outcome: three developed PE only, 11 had a SGA neonate only and four had both. A comparison of the COMP group vs the no-COMP group found no differences in maternal age, body mass index, nulliparity, gravidity or race. For 66 of the 74 subjects, US data were also available. In these subjects, both median MRI-PI (0.95 vs 0.70; P < 0.01) and median US-PI (0.95 vs 0.73; P < 0.01) were significantly increased in subjects in the COMP group compared with those in the no-COMP group. The UtA blood-flow rate, as measured by MRI, did not increase significantly from the second to the third trimester (median flow (interquartile range (IQR)), 543 (419-698) vs 575 (440-746) mL/min; P = 0.77), but it was significantly lower overall in the COMP compared with the no-COMP group (median flow (IQR), 486 (366-598) vs 624 (457-749) mL/min; P = 0.04). The areas under the receiver-operating-characteristics curves for MRI-flow, MRI-PI and US-PI in predicting COMP were not significantly different (0.694, 0.737 and 0.731, respectively; P = 0.87).

CONCLUSIONS

4D flow MRI can yield physiological measures of UtA blood-flow rate and PI that are associated with adverse pregnancy outcome. This may open up new avenues in the future to expand the potential of this technique as a robust tool with which to evaluate UtA hemodynamics in pregnancy. © 2020 International Society of Ultrasound in Obstetrics and Gynecology.

Qualitative and Quantitative Assessments of Blood Flow on Tears in Type B Aortic Dissection With Different Morphologies

Objective: The interactions between aortic morphology and hemodynamics play a key role in determining type B aortic dissection (TBAD) progression and remodeling. The study aimed to provide qualitative and quantitative hemodynamic assessment in four different TBAD morphologies based on 4D flow MRI analysis. Materials and Methods: Four patients with different TBAD morphologies underwent CT and 4D flow MRI scans. Qualitative blood flow evaluation was performed by visualizing velocity streamlines and flow directionality near the tears. Quantitative analysis included flow rate, velocity and reverse flow index (RFI) measurements. Statistical analysis was performed to evaluate hemodynamic differences between the true lumen (TL) and false lumen (FL) of patients. Results: Qualitative analysis revealed blood flow splitting near the primary entry tears (PETs), often causing the formation of vortices in the FL. All patients exhibited clear hemodynamic differences between TL and FL, with the TL generally showing higher velocities and flow rates, and lower RFIs. Average velocity magnitude measurements were significantly different for Patient 1 (t = 5.61, p = 0.001), Patient 2 (t = 3.09, p = 0.02) and Patient 4 (t = 2.81, p = 0.03). At follow-up, Patient three suffered from left renal ischemia because of FL collapse. This patient presented a complex morphology with two FLs and marked flow differences between TL and FLs. In Patient 4, left renal artery malperfusion was observed at the 32-months follow-up, due to FL thrombosis growing after PET repair. Conclusion: The study demonstrates the clinical feasibility of using 4D flow MRI in the context of TBAD. Detailed patient-specific hemodynamics assessment before treatment may provide useful insights to better understand this pathology in the future.

The Influence of Respiration on Blood Flow in the Fontan Circulation: Insights for Imaging-Based Clinical Evaluation of the Total Cavopulmonary Connection

Congenital heart disease is the most common birth defect and functionally univentricular heart defects represent the most severe end of this spectrum. The Fontan circulation provides an unique solution for single ventricle patients, by connecting both caval veins directly to the pulmonary arteries. As a result, the pulmonary circulation in Fontan palliated patients is characterized by a passive, low-energy circulation that depends on increased systemic venous pressure to drive blood toward the lungs. The absence of a subpulmonary ventricle led to the widely believed concept that respiration, by sucking blood to the pulmonary circulation during inspiration, is of great importance as a driving force for antegrade blood flow in Fontan patients. However, recent studies show that respiration influences pulsatility, but has a limited effect on net forward flow in the Fontan circulation. Importantly, since MRI examination is recommended every 2 years in Fontan patients, clinicians should be aware that most conventional MRI flow sequences do not capture the pulsatility of the blood flow as a result of the respiration. In this review, the unique flow dynamics influenced by the cardiac and respiratory cycle at multiple locations within the Fontan circulation is discussed. The impact of (not) incorporating respiration in different MRI flow sequences on the interpretation of clinical flow parameters will be covered. Finally, the influence of incorporating respiration in advanced computational fluid dynamic modeling will be outlined.

Quantification of peak blood flow velocity at the cardiac valve and great thoracic vessels by four-dimensional flow and two-dimensional phase-contrast MRI compared with echocardiography: a systematic review and meta-analysis

AIM

To objectively examine the agreement and correlation between four-dimensional (4D) flow magnetic resonance imaging (MRI) and traditional two-dimensional (2D) phase-contrast (PC) MRI with the reference standard of Doppler echocardiography for measuring peak blood velocity at the cardiac valve and great arteries, and to assess if 4D flow MRI offers an advantage over the traditional 2D method.

MATERIALS AND METHODS

The literature was searched systematically for studies that evaluate the degree of correlation and agreement between 4D flow MRI or 2D PC MRI and Doppler retrieved from PubMed, EMBASE, and the Cochrane Library. A meta-analysis was conducted to determine the peak velocity pooled bias with 95% limits of agreement (LoA) and correlation coefficient (r) for 4D flow MRI and 2D PC MRI compared with Doppler.

RESULTS

Ten studies that compared 4D flow MRI with Doppler and 12 studies that compared 2D PC MRI with Doppler were included. 4D flow MRI showed an underestimation with bias and 95% LoA of -0.09 (-0.41, 0.24) m/s (p=0.079) while 2D PC MRI showed a poorer agreement with a bias and 95% LoA of -0.25 (-0.53, 0.03), p=0.596. 4D flow MRI and 2D PC MRI showed a strong correlation with R=0.80 (95% CI 0.75, 0.84; p<0.001) and R=0.83 (95% CI 0.79, 0.87; p<0.001), respectively.

CONCLUSION

In this meta-analysis, 4D flow MRI provides improved assessment of peak velocity when compared with traditional 2D PC MRI. 4D flow MRI can be considered an important complement or substitute to Doppler echocardiography for peak velocity assessment.

Recommendation for Cardiac Magnetic Resonance Imaging-Based Phenotypic Study: Imaging Part

Cardiac magnetic resonance (CMR) imaging provides important biomarkers for the early diagnosis of many cardiovascular diseases and has been reported to reveal phenome-wide associations of cardiac/aortic structure and functionality in population studies. Nevertheless, due to the complexity of operation and variations among manufactural vendors, magnetic field strengths, coils, sequences, scan parameters, and image analysis approaches, CMR is rarely used in large cohort studies. Existing guidelines mainly focused on the diagnosis of cardiovascular diseases, which did not aim to basic research. The purpose of this study was to propose a recommendation for CMR based phenotype measurements for cohort study. We classify the imaging sequences of CMR into three categories according to the importance and universality of corresponding measurable phenotypes. The acquisition time and repeatability of the phenotypic measurement were also taken into consideration during the categorization. Unlike other guidelines, this recommendation focused on quantitative measurement of large amount of phenotypes from CMR.

Early Atherosclerotic Changes in Coronary Arteries are Associated with Endothelium Shear Stress Contraction/Expansion Variability

Although unphysiological wall shear stress (WSS) has become the consensus hemodynamic mechanism for coronary atherosclerosis, the complex biomechanical stimulus affecting atherosclerosis evolution is still undetermined. This has motivated the interest on the contraction/expansion action exerted by WSS on the endothelium, obtained through the WSS topological skeleton analysis. This study tests the ability of this WSS feature, alone or combined with WSS magnitude, to predict coronary wall thickness (WT) longitudinal changes. Nine coronary arteries of hypercholesterolemic minipigs underwent imaging with local WT measurement at three time points: baseline (T1), after 5.6 ± 0.9 (T2), and 7.6 ± 2.5 (T3) months. Individualized computational hemodynamic simulations were performed at T1 and T2. The variability of the WSS contraction/expansion action along the cardiac cycle was quantified using the WSS topological shear variation index (TSVI). Alone or combined, high TSVI and low WSS significantly co-localized with high WT at the same time points and were significant predictors of thickening at later time points. TSVI and WSS magnitude values in a physiological range appeared to play an atheroprotective role. Both the variability of the WSS contraction/expansion action and WSS magnitude, accounting for different hemodynamic effects on the endothelium, (1) are linked to WT changes and (2) concur to identify WSS features leading to coronary atherosclerosis.

Clinical Applications of 4D Flow MR Imaging in Aortic Valvular and Congenital Heart Disease

4D flow MRI allows time-resolved 3D velocity-encoded phase-contrast imaging for 3D visualization and quantification of aortic and intracardiac flow. Radiologists should be familiar with the principles of 4D flow MRI and methods for evaluating blood flow qualitatively and quantitatively. The most substantial benefits of 4D flow MRI are that it enables the simultaneous comprehensive assessment of different vessels, and that retrospective analysis can be achieved in all vessels in any direction in the field of view, which is especially beneficial for patients with complicated congenital heart disease (CHD). For aortic valvular diseases, new parameters such as wall shear stress and energy loss may provide new prognostic values for 4D flow MRI. In this review, we introduce the clinical applications of 4D flow MRI for the visualization of blood flow and quantification of hemodynamic metrics in the setting of aortic valvular disease and CHD, including intracardiac shunt and coronary artery anomaly.

Haben medizinische Kompressionsstrümpfe im Liegen einen hämodynamischen Effekt?

Medizinische Kompressionsstrümpfe (MKS) bilden die Grundlage in der Therapie chronischer Venenerkrankungen. Aufgrund mangelnder Studiendaten über deren Auswirkung auf die Hämodynamik im Liegen kann für ihr Tragen in horizontaler Körperlage bisher jedoch keine klare Empfehlung gegeben werden. In folgender Kasuistik wurde der hämodynamische Effekt von MKS der Klasse 1 bei einer Patientin mit venentypischen Beschwerden im Stadium C4 über die Erhebung des veno-arteriellen Flow-Index (VAFI) quantitativ untersucht. Unter Kompression mit MKS zeigten sowohl die duplexsonografischen als auch die mittels Phasenkontrast-MRT akquirierten Messdaten eine deutliche Verbesserung der hämodynamischen Verhältnisse in Rückenlage. Dieser Krankheitsfall zeigt, dass das Tragen von MKS im Liegen insbesondere bei Immobilität einen therapeutischen Mehrwert bringen kann.

Numerical analysis of stenoses severity and aortic wall mechanics in patients with supravalvular aortic stenosis

Supravalvular aortic stenosis (SVAS) is an aortic malformation characterized by a narrowing of the ascending aorta, resulting in abnormal hemodynamics and pressure drop across the stenosed region. It has been observed that the pressure drops measured from Doppler ultrasound exams often tend to be higher than those obtained from invasive cardiac catheterization. These misleadingly elevated pressure measurements may drive the decision to refer patients for surgical treatment prematurely. Considering this strong clinical association, the purpose of this work is to develop a computational modeling approach using a two-way coupled fluid-structure interaction methodology to determine an accurate prediction of trans-stenotic pressure drop and to further highlight the discrepancy between the SVAS assessment methods. Blood is modeled using Navier-Stokes equations while the aortic wall is simulated by a composite poroelastic structure to represent the three main layers of the arterial wall. The relationship between aortic wall elasticity and the blood flow conditions is examined in varying levels of stenosis, ranging from mild to severe degrees of vessel diameter narrowing. A substantial overestimation of the traditional Doppler pressure drop measurement is observed, especially for severe stenosis levels. The simulation results indicate that elasticity of the aortic wall has a relatively little effect on trans-stenotic pressure drop for the range of mild to moderate SVAS cases, but predicted to have a profound effect for severe SVAS cases. Moreover, significant sensitivity to the pressure drop across the SVAS region from stenosis severity is observed.

DiSpect: Displacement spectrum imaging of flow and tissue perfusion using spin-labeling and stimulated echoes

PURPOSE

We propose a new method, displacement spectrum (DiSpect) imaging, for probing in vivo complex tissue dynamics such as motion, flow, diffusion, and perfusion. Based on stimulated echoes and image phase, our flexible approach enables observations of the spin dynamics over short (milliseconds) to long (seconds) evolution times.

METHODS

The DiSpect method is a Fourier-encoded variant of displacement encoding with stimulated echoes, which encodes bulk displacement of spins that occurs between tagging and imaging in the image phase. However, this method fails to capture partial volume effects as well as blood flow. The DiSpect variant mitigates this by performing multiple scans with increasing displacement-encoding steps. Fourier analysis can then resolve the multidimensional spectrum of displacements that spins exhibit over the mixing time. In addition, repeated imaging following tagging can capture dynamic displacement spectra with increasing mixing times.

RESULTS

We demonstrate properties of DiSpect MRI using flow phantom experiments as well as in vivo brain scans. Specifically, the ability of DiSpect to perform retrospective vessel-selective perfusion imaging at multiple mixing times is highlighted.

CONCLUSION

The DiSpect variant is a new tool in the arsenal of MRI techniques for probing complex tissue dynamics. The flexibility and the rich information it provides open the possibility of alternative ways to quantitatively measure numerous complex spin dynamics, such as flow and perfusion within a single exam.

Phase-Contrast MRI: Physics, Techniques, and Clinical Applications

With phase-contrast imaging, the MRI signal is used to visualize and quantify velocity. This imaging modality relies on phase data, which are intrinsic to all MRI signals. With use of bipolar gradients, degrees of phase shift are encoded and in turn correlated directly with the velocity of protons. The acquisition of diagnostic-quality images requires selection of the correct imaging plane to ensure accurate measurement and selection of the encoding velocity and thus prevent aliasing and achieve the highest signal-to-noise ratio. Multiple applications of phase-contrast imaging are actively used in clinical practice. One of the most common clinical uses is in cardiac valvular flow imaging, at which the data are used to assess the severity of valvular disease and quantify the shunt fraction. In neurologic imaging, phase-contrast imaging can be used to measure the flow of cerebrospinal fluid. This measurement can aid in the diagnosis and direct management of normal pressure hydrocephalus or be used to evaluate the severity of stenosis, such as that in Chiari I malformations. At vascular analysis, phase-contrast imaging can be used to visualize arterial and venous flow, and this application is used most commonly in the brain. Three-dimensional imaging can yield highly detailed flow data in a technique referred to as four-dimensional flow. A more recently identified application is in MR elastography. Shear waves created by using an impulse device can be velocity encoded, and this velocity is directly proportional to the stiffness of the organ, or the shear modulus. This imaging modality is most commonly used in the liver for evaluation of cirrhosis and steatosis, although research on the assessment of other organs is being performed. Phase-contrast imaging is an important tool in the arsenal of MRI examinations and has many applications. Proper use of phase-contrast imaging requires an understanding of the many practical and technical factors and unique physics principles underlying the technique.^©RSNA, 2020.

Overview of MRI for pulmonary functional imaging

Morphological evaluation of the lung is important in the clinical evaluation of pulmonary diseases. However, the disease process, especially in its early phases, may primarily result in changes in pulmonary function without changing the pulmonary structure. In such cases, the traditional imaging approaches to pulmonary morphology may not provide sufficient insight into the underlying pathophysiology. Pulmonary imaging community has therefore tried to assess pulmonary diseases and functions utilizing not only nuclear medicine, but also CT and MR imaging with various technical approaches. In this review, we overview state-of-the art MR methods and the future direction of: (1) ventilation imaging, (2) perfusion imaging and (3) biomechanical evaluation for pulmonary functional imaging.

4D flow MRI for non-invasive measurement of blood flow in the brain: A systematic review

The brain's vasculature is essential for brain health and its dysfunction contributes to the onset and development of many dementias and neurological disorders. While numerous in vivo imaging techniques exist to investigate cerebral haemodynamics in humans, phase-contrast magnetic resonance imaging (MRI) has emerged as a reliable, non-invasive method of quantifying blood flow within intracranial vessels. In recent years, an advanced form of this method, known as 4D flow, has been developed and utilised in patient studies, where its ability to capture complex blood flow dynamics within any major vessel across the acquired volume has proved effective in collecting large amounts of information in a single scan. While extremely promising as a method of examining the vascular system's role in brain-related diseases, the collection of 4D data can be time-consuming, meaning data quality has to be traded off against the acquisition time. Here, we review the available literature to examine 4D flow's capabilities in assessing physiological and pathological features of the cerebrovascular system. Emerging techniques such as dynamic velocity-encoding and advanced undersampling methods, combined with increasingly high-field MRI scanners, are likely to bring 4D flow to the forefront of cerebrovascular imaging studies in the years to come.

Complicated Double-Orifice Mitral Regurgitation: Combined Hemodynamic Assessment Using Echocardiography and Four-Dimensional Flow Magnetic Resonance Imaging

•

Angle correction may be needed when PISA is applied to a near commissural MR jet.

•

Four-dimensional flow MRI enabled MR volume quantification with 3D jet visualization.

•

Exercise-induced VT may have complex structural, genetic, and ECG etiologies.

Fully self-gated whole-heart 4D flow imaging from a 5-minute scan

PURPOSE

To develop and validate an acquisition and processing technique that enables fully self-gated 4D flow imaging with whole-heart coverage in a fixed 5-minute scan.

THEORY AND METHODS

The data are acquired continuously using Cartesian sampling and sorted into respiratory and cardiac bins using the self-gating signal. The reconstruction is performed using a recently proposed Bayesian method called ReVEAL4D. ReVEAL4D is validated using data from 8 healthy volunteers and 2 patients and compared with compressed sensing technique, L1-SENSE.

RESULTS

Healthy subjects-Compared with 2D phase-contrast MRI (2D-PC), flow quantification from ReVEAL4D shows no significant bias. In contrast, the peak velocity and peak flow rate for L1-SENSE are significantly underestimated. Compared with traditional parallel MRI-based 4D flow imaging, ReVEAL4D demonstrates small but significant biases in net flow and peak flow rate, with no significant bias in peak velocity. All 3 indices are significantly and more markedly underestimated by L1-SENSE. Patients-Flow quantification from ReVEAL4D agrees well with the 2D-PC reference. In contrast, L1-SENSE markedly underestimated peak velocity.

CONCLUSIONS

The combination of highly accelerated 5-minute Cartesian acquisition, self-gating, and ReVEAL4D enables whole-heart 4D flow imaging with accurate flow quantification.

Haemodynamic assessment of bicuspid aortic valve aortopathy: a systematic review of the current literature

Both genetic and haemodynamic theories explain the aetiology, progression and optimal management of bicuspid aortic valve aortopathy. In recent years, the haemodynamic theory has been explored with the help of magnetic resonance imaging and computational fluid dynamics. The objective of this review was to summarize the findings of these investigations with focus on the blood flow pattern and associated variables, including flow eccentricity, helicity, flow displacement, cusp opening angle, systolic flow angle, wall shear stress (WSS) and oscillatory shear index. A structured literature review was performed from January 1990 to January 2018 and revealed the following 3 main findings: (i) the bicuspid aortic valve is associated with flow eccentricity and helicity in the ascending aorta compared to healthy and diseased tricuspid aortic valve, (ii) flow displacement is easier to obtain than WSS and has been shown to correlate with valve morphology and type of aortopathy and (iii) the stenotic bicuspid aortic valve is associated with elevated WSS along the greater curvature of the ascending aorta, where aortic dilatation and aortic wall thinning are commonly found. We conclude that new haemodynamic variables should complement ascending aorta diameter as an indicator for disease progression and the type and timing of intervention. WSS describes the force that blood flow exerts on the vessel wall as a function of viscosity and geometry of the vessel, making it a potentially more reliable marker of disease progression.

Strategies for rapid reconstruction in 3D MRI with radial data acquisition: 3D fast Fourier transform vs two-step 2D filtered back-projection

For 3D radial data reconstruction in magnetic resonance imaging (MRI), fast Fourier transform via gridding (gFFT) is widely used for its fast processing and flexibility. In comparison, conventional 3D filtered back projection (cFBP), while more robust against common radial k-space centering errors, suffers from long computation times and is less frequently used. In this study, we revisit another back-projection reconstruction strategy, namely two-step 2D filtered back-projection (tsFBP), as an alternative 3D radial MRI reconstruction method that combines computational efficiency and certain error tolerance. In order to compare the three methods (gFFT, cFBP, and tsFBP), theoretical analysis was performed to evaluate the number of computational steps involved in each method. Actual reconstruction times were also measured and compared using 3D radial-MRI data of a phantom and a human brain. Additionally, the sensitivity of tsFBP to artifacts caused by radial k-space centering errors was compared with the other methods. Compared to cFBP, tsFBP dramatically improved the reconstruction speed while retaining the benefit of tolerance to the radial k-space errors. Our study therefore suggests that tsFBP can be a promising alternative to the conventional back projection method for 3D radial MRI reconstruction.

Does the inflow velocity profile influence physiologically relevant flow patterns in computational hemodynamic models of left anterior descending coronary artery?

Patient-specific computational fluid dynamics is a powerful tool for investigating the hemodynamic risk in coronary arteries. Proper setting of flow boundary conditions in computational hemodynamic models of coronary arteries is one of the sources of uncertainty weakening the findings of in silico experiments, in consequence of the challenging task of obtaining in vivo 3D flow measurements within the clinical framework. Accordingly, in this study we evaluated the influence of assumptions on inflow velocity profile shape on coronary artery hemodynamics. To do that, (1) ten left anterior descending coronary artery (LAD) geometries were reconstructed from clinical angiography, and (2) eleven velocity profiles with realistic 3D features such as eccentricity and differently shaped (single- and double-vortex) secondary flows were generated analytically and imposed as inflow boundary conditions. Wall shear stress and helicity-based descriptors obtained prescribing the commonly used parabolic velocity profile were compared with those obtained with the other velocity profiles. Our findings indicated that the imposition of idealized velocity profiles as inflow boundary condition is acceptable as long the results of the proximal vessel segment are not considered, in LAD coronary arteries. As a pragmatic rule of thumb, a conservative estimation of the length of influence of the shape of the inflow velocity profile on LAD local hemodynamics can be given by the theoretical entrance length for cylindrical conduits in laminar flow conditions.

Validation of four-dimensional flow cardiovascular magnetic resonance for aortic stenosis assessment

The management of patients with aortic stenosis (AS) crucially depends on accurate diagnosis. The main aim of this study were to validate the four-dimensional flow (4D flow) cardiovascular magnetic resonance (CMR) methods for AS assessment. Eighteen patients with clinically severe AS were recruited. All patients had pre-valve intervention 6MWT, echocardiography and CMR with 4D flow. Of these, ten patients had a surgical valve replacement, and eight patients had successful transcatheter aortic valve implantation (TAVI). TAVI patients had invasive pressure gradient assessments. A repeat assessment was performed at 3–4 months to assess the remodelling response. The peak pressure gradient by 4D flow was comparable to an invasive pressure gradient (54 ± 26 mmHG vs 50 ± 34 mmHg, P = 0.67). However, Doppler yielded significantly higher pressure gradient compared to invasive assessment (61 ± 32 mmHG vs 50 ± 34 mmHg, P = 0.0002). 6MWT was associated with 4D flow CMR derived pressure gradient (r = −0.45, P = 0.01) and EOA (r = 0.54, P < 0.01) but only with Doppler EOA (r = 0.45, P = 0.01). Left ventricular mass regression was better associated with 4D flow derived pressure gradient change (r = 0.64, P = 0.04). 4D flow CMR offers an alternative method for non-invasive assessment of AS. In addition, 4D flow derived valve metrics have a superior association to prognostically relevant 6MWT and LV mass regression than echocardiography.

Four-dimensional-flow Magnetic Resonance Imaging of the Aortic Valve and Thoracic Aorta

Blood flow through the heart and great vessels is sensitive to time and multiple velocity directions. The assessment of its three-dimensional nature has been limited. Recent advances in magnetic resonance imaging (MRI) allow the comprehensive visualization and quantification of in vivo flow dynamics using four-dimensional (4D)-flow MRI. In addition, the technique provides the opportunity to obtain advanced hemodynamic measures. This article introduces 4D-flow MRI as it is currently used for blood flow visualization and quantification of cardiac hemodynamic parameters. It discusses its advantages relative to other flow MRI techniques and describes its potential clinical applications.

Improved acceleration of phase-contrast flow imaging with magnitude difference regularization

PURPOSE

To develop a regularized image reconstruction algorithm for improved scan acceleration of phase-contrast (PC) flow MRI.

METHODS

Based on the magnitude similarity between bipolar-encoded k-space data, magnitude-difference regularization was incorporated into the conventional compressed sensing (CS) reconstruction. The gradient of the magnitude regularization was derived so the reconstruction problem can be solved using non-linear conjugate gradient with backtracking line search. Phase contrast flow data obtained in the peripheral arteries of healthy and patient subjects were retrospectively undersampled for testing the proposed reconstruction method. Three-dimensional velocity-encoded PC flow MRI was performed with prospective 4-fold undersampling for measuring arotic flow velocity in a healthy volunteer.

RESULTS

In the femoral arteries of healthy volunteers, the root-mean-square (RMS) errors of mean velocities were 0.56 ± 0.09 cm/s with CS-only reconstruction and 0.46 ± 0.08 cm/s with addition of magnitude regularization for three-fold acceleration; 1.34 ± 0.17 cm/s (CS only) and 1.08 ± 0.15 cm/s (magnitude regularized) for four-fold acceleration. In the iliac arteries of the patient, the RMS errors of mean velocities were 0.72 ± 0.12 cm/s and 0.56 ± 0.10 for three-fold acceleration, and 1.75 ± 0.21 and 1.24 ± 0.19 cm/s for four-fold acceleration (in the order of CS-only and magnitude regularized reconstructions). In the popliteal arteries, the RMS errors were 0.61 ± 0.10 cm/s and 0.42 ± 0.11 for three-fold acceleration, and 1.41 ± 0.19 and 1.12 ± 0.17 cm/s for four-fold acceleration. The maximum through-plane mean flow velocities were measured as 63.2 cm/s and 84.5 cm/s in ascending and descending aortas, respectively.

CONCLUSION

The addition of magnitude-difference regularization into conventional CS reconstruction improves the accuracy of image reconstruction using highly undersampled phase-contrast flow MR data.

On the limitations of echo planar 4D flow MRI

PURPOSE

To compare EPI and GRE readout in high-flow velocity regimes and evaluate their impact on measurement accuracy in silico and in vitro.

THEORY AND METHODS

Phase-contrast sequences for EPI and GRE were simulated using CFD velocity data to assess displacement artifacts as well as effective spatial resolution. In silico findings were validated experimentally using a steady flow phantom.

RESULTS

For EPI factor 5 and simulated stenotic flow with peak velocity of 2.2 m s - 1 , displacement artifacts resulted in misregistration of 7.3 mm at echo time and the effective resolution was locally reduced by factors 5 and 8 compared to GRE for flow along phase and frequency encoding directions, respectively. In vitro, a maximum velocity difference between EPI factor 5 and GRE of 0.97 m s - 1 was found.

CONCLUSIONS

Four-dimensional flow MRI using EPI readout results not only in considerable velocity misregistration but also in spatially varying degradation of resolution. The proposed work indicates that EPI is inferior to standard GRE for 4D flow MRI.

Magnetic Resonance Angiography of the Thoracic Vasculature: Technique and Applications

Magnetic resonance angiography (MRA) is a powerful clinical tool for evaluation of the thoracic vasculature. MRA can be performed on nearly any magnetic resonance imaging (MRI) scanner, and provides images of high diagnostic quality without the use of ionizing radiation. While computed tomographic angiography (CTA) is preferred in the evaluation of hemodynamically unstable patients, MRA represents an important tool for evaluation of the thoracic vasculature in stable patients. Contrast-enhanced MRA is generally performed unless there is a specific contraindication, as it shortens the duration of the exam and provides images of higher diagnostic quality than noncontrast MRA. However, intravenous contrast is often not required to obtain a diagnostic evaluation for most clinical indications. Indeed, a variety of noncontrast MRA techniques are used for thoracic imaging, often in conjunction with contrast-enhanced MRA, each of which has a differing degree of reliance on flowing blood to produce the desired vascular signal. In this article we review contrast-enhanced MRA, with a focus on contrast agents, methods of bolus timing, and considerations in imaging acquisition. Next, we cover the mechanism of contrast, strengths, and weaknesses of various noncontrast MRA techniques. Finally, we present an approach to protocol development and review representative protocols used at our institution for a variety of thoracic applications. Further attention will be devoted to additional techniques employed to address specific clinical questions, such as delayed contrast-enhanced imaging, provocative maneuvers, electrocardiogram and respiratory gating, and phase-contrast imaging. The purpose of this article is to review basic techniques and methodology in thoracic MRA, discuss an approach to protocol development, and illustrate commonly encountered pathology on thoracic MRA examinations. Level of Evidence 5 Technical Efficacy Stage 3.

The impact of helical flow on coronary atherosclerotic plaque development

BACKGROUND AND AIMS

Atherosclerosis has been associated with near-wall hemodynamics and wall shear stress (WSS). However, the role of coronary intravascular hemodynamics, in particular of the helical flow (HF) patterns that physiologically develop in those arteries, is rarely considered. The purpose of this study was to assess how HF affects coronary plaque initiation and progression, definitively demonstrating its atheroprotective nature.

METHODS

The three main coronary arteries of five adult hypercholesterolemic mini-pigs on a high fat diet were imaged by computed coronary tomography angiography (CCTA) and intravascular ultrasound (IVUS) at 3 (T1, baseline) and 9.4 ± 1.9 (T2) months follow-up. The baseline geometries of imaged coronary arteries (n = 15) were reconstructed, and combined with pig-specific boundary conditions (based on in vivo Doppler blood flow measurements) to perform computational fluid dynamic simulations. Local wall thickness (WT) was measured on IVUS images at T1 and T2, and its temporal changes were assessed. Descriptors of HF and WSS nature were computed for each model, and statistically compared to WT data.

RESULTS

HF intensity was strongly positively associated with WSS magnitude (p < 0.001). Overall, coronary segments exposed to high baseline levels of HF intensity exhibited a significantly lower WT growth (p < 0.05), compared to regions with either mid or low HF intensity.

CONCLUSIONS

This study confirms the physiological significance of HF in coronary arteries, revealing its protective role against atherosclerotic WT growth and its potential in predicting regions undergoing WT development. These findings support future in vivo measurement of coronary HF as atherosclerotic risk marker, overcoming current limitations of in vivo WSS assessment.

Segmentation of carotid arterial walls using neural networks

BACKGROUND

Automated, accurate, objective, and quantitative medical image segmentation has remained a challenging goal in computer science since its inception. This study applies the technique of convolutional neural networks (CNNs) to the task of segmenting carotid arteries to aid in the assessment of pathology.

AIM

To investigate CNN’s utility as an ancillary tool for researchers who require accurate segmentation of carotid vessels.

METHODS

An expert reader delineated vessel wall boundaries on 4422 axial T2-weighted magnetic resonance images of bilateral carotid arteries from 189 subjects with clinically evident atherosclerotic disease. A portion of this dataset was used to train two CNNs (one to segment the vessel lumen and the other to segment the vessel wall) with the remaining portion used to test the algorithm’s efficacy by comparing CNN segmented images with those of an expert reader.

RESULTS

Overall quantitative assessment between automated and manual segmentations was determined by computing the DICE coefficient for each pair of segmented images in the test dataset for each CNN applied. The average DICE coefficient for the test dataset (CNN segmentations compared to expert’s segmentations) was 0.96 for the lumen and 0.87 for the vessel wall. Pearson correlation values and the intra-class correlation coefficient (ICC) were computed for the lumen (Pearson = 0.98, ICC = 0.98) and vessel wall (Pearson = 0.88, ICC = 0.86) segmentations. Bland-Altman plots of area measurements for the CNN and expert readers indicate good agreement with a mean bias of 1%-8%.

CONCLUSION

Although the technique produces reasonable results that are on par with expert human assessments, our application requires human supervision and monitoring to ensure consistent results. We intend to deploy this algorithm as part of a software platform to lessen researchers’ workload to more quickly obtain reliable results.

Comprehensive MR Analysis of Cardiac Function, Aortic Hemodynamics and Left Ventricular Strain in Pediatric Cohort with Isolated Bicuspid Aortic Valve

Bicuspid aortic valve (BAV) disease demonstrates a range of clinical presentations and complications. We aim to use cardiac MRI (CMR) to evaluate left ventricular (LV) parameters, myocardial strain and aortic hemodynamics in pediatric BAV patients with and without aortic stenosis (AS) or regurgitation (AR) compared to tricuspid aortic valve (TAV) controls. We identified 58 pediatric BAV patients without additional cardiovascular pathology and 25 healthy TAV controls (15.3 ± 2.2 years) who underwent CMR with 4D flow. BAV cohort included subgroups with no valvulopathy (n = 13, 14.3 ± 4.7 years), isolated AS (n = 19, 14.5 ± 4.0 years), mixed valve disease (AS + AR) (n = 13, 17.1 ± 3.2 years), and prior valvotomy/valvuloplasty (n = 13, 13.9 ± 3.2 years). CMR data included LV volumetric and mass indices, myocardial strain and aortic hemodynamics. BAV patients with no valvulopathy or isolated AS had similar LV parameters to controls excepting cardiac output (p < 0.05). AS + AR and post-surgical patients had abnormal LV volumetric and mass indices (p < 0.01). Post-surgical patients had decreased global longitudinal strain (p = 0.02); other subgroups had comparable strain to controls. Patients with valvulopathy demonstrated elevated velocity and wall shear stress (WSS) in the ascending aorta (AAo) and arch (p < 0.01), while those without valve dysfunction had only elevated AAo velocity (p = 0.03). Across the cohort, elevated AAo velocity and WSS correlated to higher LV mass (p < 0.01), and abnormal hemodynamics correlated to decreased strain rates (p < 0.045). Pediatric BAV patients demonstrate abnormalities in LV parameters as a function of valvular dysfunction, most significantly in children with AS + AR or prior valvotomy/valvuloplasty. Correlations between aortic hemodynamics, LV mass and strain suggest valvular dysfunction could drive LV remodeling. Multiparametric CMR assessment in pediatric BAV may help stratify risk for cardiac remodeling and dysfunction.