Four-Dimensional Flow Magnetic Resonance Imaging for Assessment of Velocity Magnitudes and Flow Patterns in The Human Carotid Artery Bifurcation: Comparison with Computational Fluid Dynamics

Four-dimensional flow MRI for quantitative assessment of cerebrospinal fluid dynamics: Status and opportunities

Neurological disorders can manifest with altered neurofluid dynamics in different compartments of the central nervous system. These include alterations in cerebral blood flow, cerebrospinal fluid (CSF) flow, and tissue biomechanics. Noninvasive quantitative assessment of neurofluid flow and tissue motion is feasible with phase contrast magnetic resonance imaging (PC MRI). While two-dimensional (2D) PC MRI is routinely utilized in research and clinical settings to assess flow dynamics through a single imaging slice, comprehensive neurofluid dynamic assessment can be limited or impractical. Recently, four-dimensional (4D) flow MRI (or time-resolved three-dimensional PC with three-directional velocity encoding) has emerged as a powerful extension of 2D PC, allowing for large volumetric coverage of fluid velocities at high spatiotemporal resolution within clinically reasonable scan times. Yet, most 4D flow studies have focused on blood flow imaging. Characterizing CSF flow dynamics with 4D flow (i.e., 4D CSF flow) is of high interest to understand normal brain and spine physiology, but also to study neurological disorders such as dysfunctional brain metabolite waste clearance, where CSF dynamics appear to play an important role. However, 4D CSF flow imaging is challenged by the long T1 time of CSF and slower velocities compared with blood flow, which can result in longer scan times from low flip angles and extended motion-sensitive gradients, hindering clinical adoption. In this work, we review the state of 4D CSF flow MRI including challenges, novel solutions from current research and ongoing needs, examples of clinical and research applications, and discuss an outlook on the future of 4D CSF flow.

Integrating PointNet-Based Model and Machine Learning Algorithms for Classification of Rupture Status of IAs

BACKGROUND

The rupture of intracranial aneurysms (IAs) would result in subarachnoid hemorrhage with high mortality and disability. Predicting the risk of IAs rupture remains a challenge.

METHODS

This paper proposed an effective method for classifying IAs rupture status by integrating a PointNet-based model and machine learning algorithms. First, medical image segmentation and reconstruction algorithms were applied to 3D Digital Subtraction Angiography (DSA) imaging data to construct three-dimensional IAs geometric models. Geometrical parameters of IAs were then acquired using Geomagic, followed by the computation of hemodynamic clouds and hemodynamic parameters using Computational Fluid Dynamics (CFD). A PointNet-based model was developed to extract different dimensional hemodynamic cloud features. Finally, five types of machine learning algorithms were applied on geometrical parameters, hemodynamic parameters, and hemodynamic cloud features to classify and recognize IAs rupture status. The classification performance of different dimensional hemodynamic cloud features was also compared.

RESULTS

The 16-, 32-, 64-, and 1024-dimensional hemodynamic cloud features were extracted with the PointNet-based model, respectively, and the four types of cloud features in combination with the geometrical parameters and hemodynamic parameters were respectively applied to classify the rupture status of IAs. The best classification outcomes were achieved in the case of 16-dimensional hemodynamic cloud features, the accuracy of XGBoost, CatBoost, SVM, LightGBM, and LR algorithms was 0.887, 0.857, 0.854, 0.857, and 0.908, respectively, and the AUCs were 0.917, 0.934, 0.946, 0.920, and 0.944. In contrast, when only utilizing geometrical parameters and hemodynamic parameters, the accuracies were 0.836, 0.816, 0.826, 0.832, and 0.885, respectively, with AUC values of 0.908, 0.922, 0.930, 0.884, and 0.921.

CONCLUSION

In this paper, classification models for IAs rupture status were constructed by integrating a PointNet-based model and machine learning algorithms. Experiments demonstrated that hemodynamic cloud features had a certain contribution weight to the classification of IAs rupture status. When 16-dimensional hemodynamic cloud features were added to the morphological and hemodynamic features, the models achieved the highest classification accuracies and AUCs. Our models and algorithms would provide valuable insights for the clinical diagnosis and treatment of IAs.

Subjects with carotid webs demonstrate pro-thrombotic hemodynamics compared to subjects with carotid atherosclerosis

Carotid artery webs (CaW) are non-atherosclerotic projections into the vascular lumen and have been linked to up to one-third of cryptogenic strokes in younger patients. Determining how CaW affects local hemodynamics is essential for understanding clot formation and stroke risk. Computational fluid dynamics simulations were used to investigate patient-specific hemodynamics in carotid artery bifurcations with CaW, bifurcations with atherosclerotic lesions having a similar degree of lumen narrowing, and with healthy carotid bifurcations. Simulations were conducted using segmented computed tomography angiography geometries with inlet boundary conditions extracted from 2D phase contrast MRI scans. The study included carotid bifurcations with CaW (n = 13), mild atherosclerosis (n = 7), and healthy bifurcation geometries (n = 6). Hemodynamic parameters associated with vascular dysfunction and clot formation, including shear rate, oscillatory shear index (OSI), low velocity, and flow stasis were calculated and compared between the subject groups. Patients with CaW had significantly larger regions containing low shear rate, high OSI, low velocity, and flow stasis in comparison to subjects with mild atherosclerosis or normal bifurcations. These abnormal hemodynamic metrics in patients with CaW are associated with clot formation and vascular dysfunction and suggest that hemodynamic assessment may be a tool to assess stroke risk in these patients.

Turbulent blood dynamics in the left heart in the presence of mitral regurgitation: a computational study based on multi-series cine-MRI

In this work, we performed a computational image-based study of blood dynamics in the whole left heart, both in a healthy subject and in a patient with mitral valve regurgitation. We elaborated multi-series cine-MRI with the aim of reconstructing the geometry and the corresponding motion of left ventricle, left atrium, mitral and aortic valves, and aortic root of the subjects. This allowed us to prescribe such motion to computational blood dynamics simulations where, for the first time, the whole left heart motion of the subject is considered, allowing us to obtain reliable subject-specific information. The final aim is to investigate and compare between the subjects the occurrence of turbulence and the risk of hemolysis and of thrombi formation. In particular, we modeled blood with the Navier-Stokes equations in the arbitrary Lagrangian-Eulerian framework, with a large eddy simulation model to describe the transition to turbulence and a resistive method to manage the valve dynamics, and we used a finite element discretization implemented in an in-house code for the numerical solution.

Optimization of 4D Flow MRI Spatial and Temporal Resolution for Examining Complex Hemodynamics in the Carotid Artery Bifurcation

BACKGROUND

Three-dimensional, ECG-gated, time-resolved, three-directional, velocity-encoded phase-contrast MRI (4D flow MRI) has been applied extensively to measure blood velocity in great vessels but has been much less used in diseased carotid arteries. Carotid artery webs (CaW) are non-inflammatory intraluminal shelf-like projections into the internal carotid artery (ICA) bulb that are associated with complex flow and cryptogenic stroke.

PURPOSE

Optimize 4D flow MRI for measuring the velocity field of complex flow in the carotid artery bifurcation model that contains a CaW.

METHODS

A 3D printed phantom model created from computed tomography angiography (CTA) of a subject with CaW was placed in a pulsatile flow loop within the MRI scanner. 4D Flow MRI images of the phantom were acquired with five different spatial resolutions (0.50-2.00  mm3) and four different temporal resolutions (23-96 ms) and compared to a computational fluid dynamics (CFD) solution of the flow field as a reference. We examined four planes perpendicular to the vessel centerline, one in the common carotid artery (CCA) and three in the internal carotid artery (ICA) where complex flow was expected. At these four planes pixel-by-pixel velocity values, flow, and time average wall shear stress (TAWSS) were compared between 4D flow MRI and CFD.

HYPOTHESIS

An optimized 4D flow MRI protocol will provide a good correlation with CFD velocity and TAWSS values in areas of complex flow within a clinically feasible scan time (~ 10 min).

RESULTS

Spatial resolution affected the velocity values, time average flow, and TAWSS measurements. Qualitatively, a spatial resolution of 0.50  mm3 resulted in higher noise, while a lower spatial resolution of 1.50-2.00  mm3 did not adequately resolve the velocity profile. Isotropic spatial resolutions of 0.50-1.00  mm3 showed no significant difference in total flow compared to CFD. Pixel-by-pixel velocity correlation coefficients between 4D flow MRI and CFD were > 0.75 for 0.50-1.00  mm3 but were < 0.5 for 1.50 and 2.00  mm3. Regional TAWSS values determined from 4D flow MRI were generally lower than CFD and decreased at lower spatial resolutions (larger pixel sizes). TAWSS differences between 4D flow and CFD were not statistically significant at spatial resolutions of 0.50-1.00  mm3 but were different at 1.50 and 2.00 mm3. Differences in temporal resolution only affected the flow values when temporal resolution was > 48.4 ms; temporal resolution did not affect TAWSS values.

CONCLUSION

A spatial resolution of 0.74-1.00  mm3 and a temporal resolution of 23-48 ms (1-2 k-space segments) provides a 4D flow MRI protocol capable of imaging velocity and TAWSS in regions of complex flow within the carotid bifurcation at a clinically acceptable scan time.

Blood Flow Turbulence Quantification of Carotid Artery With a High-Frame Rate Vector Flow Imaging

OBJECTIVES

To assess the feasibility and performance of Turbulence (Tur) index as a quantitative tool for carotid artery flow turbulence; to detect and compare the blood flow patterns of common carotid artery (CCA) and carotid bulb (CB) at different ages and cardiac phases in healthy adults, and thus interpret the evolvement of etiology difference between CCA and CB.

METHODS

Carotid flow characteristics of 40 healthy volunteers were evaluated quantitatively by a high-frame rate vector flow imaging. Three types of flow patterns were defined depending on the distributive range of complex flow during systole in CB. Comparison of mean Tur value in CCA and CB at different age groups and cardiac phases was performed. And the correlation between Tur value and the diameter ratio of proximal internal carotid artery to common carotid artery (DRpro-ica/cca) was tested.

RESULTS

Mean Tur values in CB were remarkably higher than that in CCA, whether during systole or diastole (P < .001). Meanwhile Tur values in CB during systole were significantly higher than that during diastole (P < .001). Flow complexity of CB showed variations among 40 participants especially in systole, whereas the flow pattern of CCA was relatively consistent. Mean Tur values were positively correlated with DRpro-ica/cca in CB (ρ = 0.69, P < .05).

CONCLUSIONS

V Flow imaging provided a reliable method-Tur, for quantitative analysis of carotid blood flow. It had potential to be further applied in distinguishing complex hemodynamic characteristics in high-risk people of carotid diseases for the risk stratification of cardiovascular events.

Hemodynamics in a three-dimensional printed aortic model: a comparison of four-dimensional phase-contrast magnetic resonance and image-based computational fluid dynamics

OBJECTIVE

This study aims to compare an electrocardiogram (ECG)-gated four-dimensional (4D) phase-contrast (PC) magnetic resonance imaging (MRI) technique and computational fluid dynamics (CFD) using variables controlled in a laboratory environment to minimize bias factors.

MATERIALS AND METHODS

Data from 4D PC-MRI were compared with computational fluid dynamics using steady and pulsatile flows at various inlet velocities. Anatomically realistic models for a normal aorta, a penetrating atherosclerotic ulcer, and an abdominal aortic aneurysm were constructed using a three-dimensional printer.

RESULTS

For the normal aorta model, the errors in the peak and the average velocities were within 5%. The peak velocities of the penetrating atherosclerotic ulcer and the abdominal aortic aneurysm models displayed a more extensive range of differences because of the high-speed and vortical fluid flows generated by the shape of the blood vessel. However, the average velocities revealed only relatively minor differences.

CONCLUSIONS

This study compared the characteristics of PC-MRI and CFD through a phantom study that only included controllable experimental parameters. Based on these results, 4D PC-MRI and CFD are powerful tools for analyzing blood flow patterns in vivo. However, there is room for future developments to improve velocity measurement accuracy.

Improving Blood Flow Visualization of Recirculation Regions at Carotid Bulb in 4D Flow MRI Using Semi-Automatic Segmentation with ITK-SNAP

Assessment of carotid bulb hemodynamics using four-dimensional (4D) flow magnetic resonance imaging (MRI) requires accurate segmentation of recirculation regions that is frequently hampered by limited resolution. This study aims to improve the accuracy of 4D flow MRI carotid bulb segmentation and subsequent recirculation regions analysis. Time-of-flight (TOF) MRI and 4D flow MRI were performed on bilateral carotid artery bifurcations in seven healthy volunteers. TOF-MRI data was segmented into 3D geometry for computational fluid dynamics (CFD) simulations. ITK-SNAP segmentation software was included in the workflow for the semi-automatic generation of 4D flow MRI angiographic data. This study compared the velocities calculated at the carotid bifurcations and the 3D blood flow visualization at the carotid bulbs obtained by 4D flow MRI and CFD. By applying ITK-SNAP segmentation software, an obvious improvement in the 4D flow MRI visualization of the recirculation regions was observed. The 4D flow MRI images of the recirculation flow characteristics of the carotid artery bulbs coincided with the CFD. A reasonable agreement was found in terms of velocity calculated at the carotid bifurcation between CFD and 4D flow MRI. However, the dispersion of velocity data points relative to the local errors of measurement in 4D flow MRI remains. Our proposed strategy showed the feasibility of improving recirculation regions segmentation and the potential for reliable blood flow visualization in 4D flow MRI. However, quantitative analysis of recirculation regions in 4D flow MRI with ITK-SNAP should be enhanced for use in clinical situations.

Evaluation of Plaque Vulnerability via Combination of Hemodynamic Analysis and Simultaneous Non-Contrast Angiography and Intraplaque Hemorrhage (SNAP) Sequence for Carotid Intraplaque Hemorrhage

The purpose of this study was to assess the vulnerability of plaque using a combination of simultaneous non-contrast angiography, intraplaque hemorrhage (SNAP) sequence, and local hemodynamic analysis in an intraplaque hemorrhage (IPH), and to evaluate the association between morphological and hemodynamic factors and IPH by comparing the IPH (presence of IPH) and non-IPH (plaque with absence of IPH) groups. In total, 27 IPH patients and 27 non-IPH patients were involved in this study, and baseline characteristics were collected. For morphological factors, diameters, and areas of the internal carotid artery (ICA), external carotid artery, and common carotid artery were measured, and bifurcation angle (α) and ICA angle (β) were also measured for comparison between the IPH group and non-IPH group. For hemodynamic factors, time-averaged wall shear stress (WSS), minimum WSS, maximum WSS, and oscillatory shear index were calculated using computational fluid dynamics (CFD) simulations. For the qualitative analysis, cross-sectional images with analyzed WSS and SNAP sequences were combined to precisely assess local hemodynamics. Bifurcation angle (α) was significantly different between the IPH and non-IPH groups (39.47 degrees vs. 47.60 degrees, p = 0.041). Significantly higher time-averaged WSS, minimum WSS, and maximum WSS were observed in the IPH group compared to the non-IPH group. In the IPH group, when using the combined analysis with SNAP sequences and WSS, the WSS of the region with IPH was significantly higher than the region without IPH (2.32 vs. 1.21 Pa, p = 0.005). A smaller bifurcation angle (α) and higher time-averaged WSS, minimum WSS, and maximum WSS were associated with IPH. The combined analysis of SNAP sequences and WSS might help to evaluate the risk of carotid IPH.

Patient-Specific 3D Bioprinted Models of Developing Human Heart

The heart is the first organ to develop in the human embryo through a series of complex chronological processes, many of which critically rely on the interplay between cells and the dynamic microenvironment. Tight spatiotemporal regulation of these interactions is key in heart development and diseases. Due to suboptimal experimental models, however, little is known about the role of microenvironmental cues in the heart development. This study investigates the use of 3D bioprinting and perfusion bioreactor technologies to create bioartificial constructs that can serve as high-fidelity models of the developing human heart. Bioprinted hydrogel-based, anatomically accurate models of the human embryonic heart tube (e-HT, day 22) and fetal left ventricle (f-LV, week 33) are perfused and analyzed both computationally and experimentally using ultrasound and magnetic resonance imaging. Results demonstrate comparable flow hemodynamic patterns within the 3D space. We demonstrate endothelial cell growth and function within the bioprinted e-HT and f-LV constructs, which varied significantly in varying cardiac geometries and flow. This study introduces the first generation of anatomically accurate, 3D functional models of developing human heart. This platform enables precise tuning of microenvironmental factors, such as flow and geometry, thus allowing the study of normal developmental processes and underlying diseases.

Four-Dimensional Flow MRI of Abdominal Veins: A Systematic Review

The aim of this systematic review is to provide an overview of the use of Four-Dimensional Magnetic Resonance Imaging of vector blood flow (4D Flow MRI) in the abdominal veins. This study was composed according to the PRISMA guidelines 2009. The literature search was conducted in MEDLINE, Cochrane Library, EMBASE, and Web of Science. Quality assessment of the included studies was performed using the QUADAS-2 tool. The initial search yielded 781 studies and 21 studies were included. All studies successfully applied 4D Flow MRI in abdominal veins. Four-Dimensional Flow MRI was capable of discerning between healthy subjects and patients with cirrhosis and/or portal hypertension. The visual quality and inter-observer agreement of 4D Flow MRI were rated as excellent and good to excellent, respectively, and the studies utilized several different MRI data sampling strategies. By applying spiral sampling with compressed sensing to 4D Flow MRI, the blood flow of several abdominal veins could be imaged simultaneously in 18-25 s, without a significant loss of visual quality. Four-Dimensional Flow MRI might be a useful alternative to Doppler sonography for the diagnosis of cirrhosis and portal hypertension. Further clinical studies need to establish consensus regarding MRI sampling strategies in patients and healthy subjects.

Comparison of Hemodynamic Visualization in Cerebral Arteries: Can Magnetic Resonance Imaging Replace Computational Fluid Dynamics?

A multimodality approach was applied using four-dimensional flow magnetic resonance imaging (4D flow MRI), time-of-flight magnetic resonance angiography (TOF-MRA) signal intensity gradient (SIG), and computational fluid dynamics (CFD) to investigate the 3D blood flow characteristics and wall shear stress (WSS) of the cerebral arteries. TOF-MRA and 4D flow MRI were performed on the major cerebral arteries in 16 healthy volunteers (mean age 34.7 ± 7.6 years). The flow rate measured with 4D flow MRI in the internal carotid artery, middle cerebral artery, and anterior cerebral artery were 3.8, 2.5, and 1.2 mL/s, respectively. The 3D blood flow pattern obtained through CFD and 4D flow MRI on the cerebral arteries showed reasonable consensus. CFD delivered much greater resolution than 4D flow MRI. TOF-MRA SIG and CFD WSS of the major cerebral arteries showed reasonable consensus with the locations where the WSS was relatively high. However, the visualizations were very different between TOF-MRA SIG and CFD WSS at the internal carotid artery bifurcations, the anterior cerebral arteries, and the anterior communicating arteries. 4D flow MRI, TOF-MRA SIG, and CFD are complementary methods that can provide additional insight into the hemodynamics of the human cerebral artery.

From 2D to 4D Phase-Contrast MRI in the Neurovascular System: Will It Be a Quantum Jump or a Fancy Decoration?

Considering the crosstalk between the flow and vessel wall, hemodynamic assessment of the neurovascular system may offer a well-integrated solution for both diagnosis and management by adding prognostic significance to the standard CT/MR angiography. 4D flow MRI or time-resolved 3D velocity-encoded phase-contrast MRI has long been promising for the hemodynamic evaluation of the great vessels, but challenged in clinical studies for assessing intracranial vessels with small diameter due to long scan times and low spatiotemporal resolution. Current accelerated MRI techniques, including parallel imaging with compressed sensing and radial k-space undersampling acquisitions, have decreased scan times dramatically while preserving spatial resolution. 4D flow MRI visualized and measured 3D complex flow of neurovascular diseases such as aneurysm, arteriovenous shunts, and atherosclerotic stenosis using parameters including flow volume, velocity vector, pressure gradients, and wall shear stress. In addition to the noninvasiveness of the phase contrast technique and retrospective flow measurement through the wanted windows of the analysis plane, 4D flow MRI has shown several advantages over Doppler ultrasound or computational fluid dynamics. The evaluation of the flow status and vessel wall can be performed simultaneously in the same imaging modality. This article is an overview of the recent advances in neurovascular 4D flow MRI techniques and their potential clinical applications in neurovascular disease. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY STAGE: 3.