Standardized Evaluation of Cerebral Arteriovenous Malformations Using Flow Distribution Network Graphs and Dual-venc 4D Flow MRI

Macrovascular blood flow and microvascular cerebrovascular reactivity are regionally coupled in adolescence

Cerebrovascular imaging assessments are particularly challenging in adolescent cohorts, where not all modalities are appropriate, and rapid brain maturation alters hemodynamics at both macro- and microvascular scales. In a preliminary sample of healthy adolescents (n=12, 8-25 years), we investigated relationships between 4D flow MRI-derived blood velocity and blood flow in bilateral anterior, middle, and posterior cerebral arteries and BOLD cerebrovascular reactivity in associated vascular territories. As hypothesized, higher velocities in large arteries are associated with an earlier response to a vasodilatory stimulus (cerebrovascular reactivity delay) in the downstream territory. Higher blood flow through these arteries is associated with a larger BOLD response to a vasodilatory stimulus (cerebrovascular reactivity amplitude) in the associated territory. These trends are consistent in a case study of adult moyamoya disease. In our small adolescent cohort, macrovascular-microvascular relationships for velocity/delay and flow/CVR change with age, though underlying mechanisms are unclear. Our work emphasizes the need to better characterize this key stage of human brain development, when cerebrovascular hemodynamics are changing, and standard imaging methods offer limited insight into these processes. We provide important normative data for future comparisons in pathology, where combining macro- and microvascular assessments may better help us prevent, stratify, and treat cerebrovascular disease.

Comparing Vascular Morphology and Hemodynamics in Patients with Vein of Galen Malformations Using Intracranial 4D Flow MRI

BACKGROUND AND PURPOSE

Vein of Galen malformation (VOGM) is the most common congenital cerebrovascular malformation, and many patients have high mortality rates and poor cognitive outcomes. Quantitative diagnostic tools are needed to improve clinical outcomes, and the purpose of this study was to characterize intracranial blood flow in VOGM using quantitative 4D flow MRI.

MATERIALS AND METHODS

A prospective study of children with VOGM was conducted by acquiring 4D flow MRI to quantify total blood inflow to the brain, flow in the pathologic falcine sinus, and flow in the superior sagittal sinus. Linear regression was used to test the relationships between these flows and age, clinical status, and the mediolateral diameter of the outflow tract of the lesion through the falcine or straight sinus diameter, which is a known morphologic prognostic metric.

RESULTS

In all 11 subjects (mean age, 22 [SD,17 ] weeks), total blood flow to the brain always exceeded normal levels (mean, 1063 [SD, 403] mL/minute). Significant correlations were observed between falcine sinus flow and the mediolateral diameter of the straight or falcine sinus, the posterior cerebral artery/MCA flow ratio and age at scanning, and superior sagittal sinus flow proximal to malformation inflow and age at scanning.

CONCLUSIONS

Using 4D flow MRI, we established the hemodynamic underpinnings of the mediolateral diameter of the straight or falcine sinus and investigated metrics representing parenchymal venous drainage that could be used to monitor the normalization of hemodynamics during embolization therapy.

Unraveling diurnal and technical variability in cerebral hemodynamics from neurovascular 4D-Flow MRI

Neurovascular 4D-Flow MRI enables non-invasive evaluation of cerebral hemodynamics including measures of cerebral blood flow (CBF), vessel pulsatility index (PI), and cerebral pulse wave velocity (PWV). 4D-Flow measures have been linked to various neurovascular disorders including small vessel disease and Alzheimer's disease; however, physiological and technical sources of variability are not well established. Here, we characterized sources of diurnal physiological and technical variability in cerebral hemodynamics using 4D-Flow in a retrospective study of cognitively unimpaired older adults (N = 750) and a prospective study of younger adults (N = 10). Younger participants underwent repeated MRI sessions at 7am, 4 pm, and 10 pm. In the older cohort, having an MRI earlier on the day was significantly associated with higher CBF and lower PI. In prospective experiments, time of day significantly explained variability in CBF and PI; however, not in PWV. Test-retest experiments showed high CBF intra-session repeatability (repeatability coefficient (RPC) =7.2%), compared to lower diurnal repeatability (RPC = 40%). PI and PWV displayed similar intra-session and diurnal variability (PI intra-session RPC = 22%, RPC = 24% 7am vs 4 pm; PWV intra-session RPC = 17%, RPC = 21% 7am vs 4 pm). Overall, CBF measures showed low technical variability, supporting diurnal variability is from physiology. PI and PWV showed higher technical variability but less diurnal variability.

Automated intracranial vessel segmentation of 4D flow MRI data in patients with atherosclerotic stenosis using a convolutional neural network

Introduction

Intracranial 4D flow MRI enables quantitative assessment of hemodynamics in patients with intracranial atherosclerotic disease (ICAD). However, quantitative assessments are still challenging due to the time-consuming vessel segmentation, especially in the presence of stenoses, which can often result in user variability. To improve the reproducibility and robustness as well as to accelerate data analysis, we developed an accurate, fully automated segmentation for stenosed intracranial vessels using deep learning.

Methods

154 dual-VENC 4D flow MRI scans (68 ICAD patients with stenosis, 86 healthy controls) were retrospectively selected. Manual segmentations were used as ground truth for training. For automated segmentation, deep learning was performed using a 3D U-Net. 20 randomly selected cases (10 controls, 10 patients) were separated and solely used for testing. Cross-sectional areas and flow parameters were determined in the Circle of Willis (CoW) and the sinuses. Furthermore, the flow conservation error was calculated. For statistical comparisons, Dice scores (DS), Hausdorff distance (HD), average symmetrical surface distance (ASSD), Bland-Altman analyses, and interclass correlations were computed using the manual segmentations from two independent observers as reference. Finally, three stenosis cases were analyzed in more detail by comparing the 4D flow-based segmentations with segmentations from black blood vessel wall imaging (VWI).

Results

Training of the network took approximately 10 h and the average automated segmentation time was 2.2 ± 1.0 s. No significant differences in segmentation performance relative to two independent observers were observed. For the controls, mean DS was 0.85 ± 0.03 for the CoW and 0.86 ± 0.06 for the sinuses. Mean HD was 7.2 ± 1.5 mm (CoW) and 6.6 ± 3.7 mm (sinuses). Mean ASSD was 0.15 ± 0.04 mm (CoW) and 0.22 ± 0.17 mm (sinuses). For the patients, the mean DS was 0.85 ± 0.04 (CoW) and 0.82 ± 0.07 (sinuses), the HD was 8.4 ± 3.1 mm (CoW) and 5.7 ± 1.9 mm (sinuses) and the mean ASSD was 0.22 ± 0.10 mm (CoW) and 0.22 ± 0.11 mm (sinuses). Small bias and limits of agreement were observed in both cohorts for the flow parameters. The assessment of the cross-sectional lumen areas in stenosed vessels revealed very good agreement (ICC: 0.93) with the VWI segmentation but a consistent overestimation (bias ± LOA: 28.1 ± 13.9%).

Discussion

Deep learning was successfully applied for fully automated segmentation of stenosed intracranial vasculatures using 4D flow MRI data. The statistical analysis of segmentation and flow metrics demonstrated very good agreement between the CNN and manual segmentation and good performance in stenosed vessels. To further improve the performance and generalization, more ICAD segmentations as well as other intracranial vascular pathologies will be considered in the future.

Generalized super-resolution 4D Flow MRI-using ensemble learning to extend across the cardiovascular system

4D Flow Magnetic Resonance Imaging (4D Flow MRI) is a non-invasive measurement technique capable of quantifying blood flow across the cardiovascular system. While practical use is limited by spatial resolution and image noise, incorporation of trained super-resolution (SR) networks has potential to enhance image quality post-scan. However, these efforts have predominantly been restricted to narrowly defined cardiovascular domains, with limited exploration of how SR performance extends across the cardiovascular system; a task aggravated by contrasting hemodynamic conditions apparent across the cardiovasculature. The aim of our study was to explore the generalizability of SR 4D Flow MRI using a combination of heterogeneous training sets and dedicated ensemble learning. With synthetic training data generated across three disparate domains (cardiac, aortic, cerebrovascular), varying convolutional base and ensemble learners were evaluated as a function of domain and architecture, quantifying performance on both in-silico and acquired in-vivo data from the same three domains. Results show that both bagging and stacking ensembling enhance SR performance across domains, accurately predicting high-resolution velocities from low-resolution input data in-silico. Likewise, optimized networks successfully recover native resolution velocities from downsampled in-vivo data, as well as show qualitative potential in generating denoised SR-images from clinicallevel input data. In conclusion, our work presents a viable approach for generalized SR 4D Flow MRI, with ensemble learning extending utility across various clinical areas of interest.

Rupture-Risk Stratifying Patients with Cerebral Arteriovenous Malformations Using Quantitative Hemodynamic Flow Measurements

Arteriovenous malformations (AVMs) are high-pressure, low-resistance arterial-venous shunts without intervening capillaries. Up to 60% of AVMs present with an intracranial hemorrhage; however, noninvasive neuroimaging has increasingly diagnosed incidental AVMs. AVM management depends on weighing the lifetime rupture risk against the risks of intervention. Although AVM rupture risk relies primarily on angioarchitectural features, measuring hemodynamic flow is gaining traction. Accurate understanding of AVM hemodynamic flow parameters will help endovascular neurosurgeons and interventional neuroradiologists stratify patients by rupture risk and select treatment plans. This review examines various neuroimaging modalities and their capabilities to quantify AVM flow, as well as the relationship between AVM flow and rupture risk. Quantitative hemodynamic studies on the relationship between AVM flow and rupture risk have not reached a clear consensus; however, the preponderance of data suggests that higher arterial inflow and lower venous outflow in the AVM nidus contribute to increased hemorrhagic risk. Future studies should consider using larger sample sizes and standardized definitions of hemodynamic parameters to reach a consensus. In the meantime, classic angioarchitectural features may be more strongly correlated with AVM rupture than the amount of blood flow.

Understanding development of jugular bulb stenosis in vein of galen malformations: identifying metrics of complex flow dynamics in the cerebral venous vasculature of infants

Introduction: Computational fluid dynamics (CFD) assess biological systems based on specific boundary conditions. We propose modeling more advanced hemodynamic metrics, such as core line length (CL) and critical points which characterize complexity of flow in the context of cerebral vasculature, and specifically cerebral veins during the physiologically evolving early neonatal state of vein of Galen malformations (VOGM). CFD has not been applied to the study of arteriovenous shunting in Vein of Galen Malformations but could help illustrate the pathophysiology of this malformation. Methods: Three neonatal patients with VOGM at Boston Children's Hospital met inclusion criteria for this study. Structural MRI data was segmented to generate a mesh of the VOGM and venous outflow. Boundary condition flow velocity was derived from PC-MR sequences with arterial and venous dual velocity encoding. The mesh and boundary conditions were applied to model the cerebral venous flow. We computed flow variables including mean wall shear stress (WSSmean), mean OSI, CL, and the mean number of critical points (nCrPointsmean) for each patient specific model. A critical point is defined as the location where the shear stress vector field is zero (stationary point) and can be used to describe complexity of flow. Results: The division of flow into the left and right venous outflow was comparable between PC-MR and CFD modeling. A high complexity recirculating flow pattern observed on PC-MR was also identified on CFD modeling. Regions of similar WSSmean and OSImean (<1.3 fold) in the left and right venous outflow channels of a single patient have several-fold magnitude difference in higher order hemodynamic metrics (> 3.3 fold CL, > 1.7 fold nCrPointsmean). Specifically, the side which developed JBS in each model had greater nCrPointsmean compared to the jugular bulb with no stenosis (VOGM1: 4.49 vs. 2.53, VOGM2: 1.94 vs. 0, VOGM3: 1 vs. 0). Biologically, these regions had subsequently divergent development, with increased complexity of flow associating with venous stenosis. Discussion: Advanced metrics of flow complexity identified in computational models may reflect observed flow phenomena not fully characterized by primary or secondary hemodynamic parameters. These advanced metrics may indicate physiological states that impact development of jugular bulb stenosis in VOGM.

Cerebrovascular super-resolution 4D Flow MRI - Sequential combination of resolution enhancement by deep learning and physics-informed image processing to non-invasively quantify intracranial velocity, flow, and relative pressure

The development of cerebrovascular disease is tightly coupled to regional changes in intracranial flow and relative pressure. Image-based assessment using phase contrast magnetic resonance imaging has particular promise for non-invasive full-field mapping of cerebrovascular hemodynamics. However, estimations are complicated by the narrow and tortuous intracranial vasculature, with accurate image-based quantification directly dependent on sufficient spatial resolution. Further, extended scan times are required for high-resolution acquisitions, and most clinical acquisitions are performed at comparably low resolution (>1 mm) where biases have been observed with regard to the quantification of both flow and relative pressure. The aim of our study was to develop an approach for quantitative intracranial super-resolution 4D Flow MRI, with effective resolution enhancement achieved by a dedicated deep residual network, and with accurate quantification of functional relative pressures achieved by subsequent physics-informed image processing. To achieve this, our two-step approach was trained and validated in a patient-specific in-silico cohort, showing good accuracy in estimating velocity (relative error: 15.0 ± 0.1%, mean absolute error (MAE): 0.07 ± 0.06 m/s, and cosine similarity: 0.99 ± 0.06 at peak velocity) and flow (relative error: 6.6 ± 4.7%, root mean square error (RMSE): 0.56 mL/s at peak flow), and with the coupled physics-informed image analysis allowing for maintained recovery of functional relative pressure throughout the circle of Willis (relative error: 11.0 ± 7.3%, RMSE: 0.3 ± 0.2 mmHg). Furthermore, the quantitative super-resolution approach is applied to an in-vivo volunteer cohort, effectively generating intracranial flow images at <0.5 mm resolution and showing reduced low-resolution bias in relative pressure estimation. Our work thus presents a promising two-step approach to non-invasively quantify cerebrovascular hemodynamics, being applicable to dedicated clinical cohorts in the future.

New approaches for brain arteriovenous malformations-related epilepsy

BACKGROUND

The purpose of this review is to present the current literature and to highlight the most recent findings in brain arteriovenous malformations (bAVM)-related epilepsy research.

METHODS

We searched Medline, PubMed, Biblioinserm, Cochrane Central to study the latest research reports about the different factors that could be responsible for the genesis of bAVM-related epilepsy. We analyzed if epileptogenesis has any characteristics traits and its relation with the vascular malformation. The results of different treatments on epilepsy were considered. Typical errors that may lead towards incorrect or worse management of the seizures for these patients were also examined.

RESULTS

The development of bAVM results from multifactorial etiologies and bAVM-related epileptogenesis is likely specific for this pathology. Different types of evidence demonstrate a bidirectional relationship between bAVM and epilepsy. Currently, there is not enough published data to determine what may be the right management for these patients.

CONCLUSIONS

A better understanding of epileptogenesis in conjunction with knowledge of the complex alterations of structures and functions following bAVM-related seizures is necessary. Identification of biomarkers that can identify subgroups most likely to benefit from a specific intervention are needed to help guide clinical management. A new concept for the treatment of epilepsy related to an unruptured bAVM that cannot be treated invasively is proposed as well as new therapeutic perspectives. The next necessary step will be to propose additional algorithms to improve the development of future trials.

Quantified flow and angioarchitecture show similar associations with hemorrhagic presentation of brain arteriovenous malformations

INTRODUCTION

The purpose of our study was to elucidate the impact of brain arteriovenous malformation (BAVM) flow and wall shear stress (WSS) on angioarchitecture and to evaluate their association with hemorrhagic presentations.

MATERIALS AND METHODS

Forty-one patients with BAVMs were evaluated by phase-contrast MR angiography. Volume flow rate and WSS were quantified. Angioarchitectural features such as location, angiogenesis, venous stenosis, venous ectasia, venous phlebitis, venous rerouting, exclusive deep vein and venous sac were evaluated by two neuroradiologists. The correlation between BAVM flow and size was evaluated with Spearman correlation coefficients. Differences of size, flow, and WSS between the hemorrhagic and non-hemorrhagic groups, the seizure and non-seizure groups, and between the different groups based on angioarchitecture were evaluated with Mann-Whitney U tests. Accuracy in predicting hemorrhage was evaluated with receiver operating characteristic curves.

RESULT

BAVM flow was highly correlated with volume (ρ = 0.77). Higher flow was more commonly associated with angiogenesis, venous ectasia, venous rerouting, and venous phlebitis. Flow and angioarchitecture showed similar efficacy in differentiating hemorrhagic from non-hemorrhagic BAVMs. WSS did not demonstrate differences across any clinical groups.

CONCLUSION

Flow quantification and angioarchitecture analysis of BAVMs showed similar efficacy as evaluated by associations with hemorrhagic presentation. High flow affects both arterial and venous angioarchitecture, reflecting the nature of low vascular resistance in BAVMs.

Motion-corrected 4D-Flow MRI for neurovascular applications

Neurovascular 4D-Flow MRI has emerged as a powerful tool for comprehensive cerebrovascular hemodynamic characterization. Clinical studies in at risk populations such as aging adults indicate hemodynamic markers can be confounded by motion-induced bias. This study develops and characterizes a high fidelity 3D self-navigation approach for retrospective rigid motion correction of neurovascular 4D-Flow data. A 3D radial trajectory with pseudorandom ordering was combined with a multi-resolution low rank regularization approach to enable high spatiotemporal resolution self-navigators from extremely undersampled data. Phantom and volunteer experiments were performed at 3.0T to evaluate the ability to correct for different amounts of induced motions. In addition, the approach was applied to clinical-research exams from ongoing aging studies to characterize performance in the clinical setting. Simulations, phantom and volunteer experiments with motion correction produced images with increased vessel conspicuity, reduced image blurring, and decreased variability in quantitative measures. Clinical exams revealed significant changes in hemodynamic parameters including blood flow rates, flow pulsatility index, and lumen areas after motion correction in probed cerebral arteries (Flow: P<0.001 Lt ICA, P=0.002 Rt ICA, P=0.004 Lt MCA, P=0.004 Rt MCA; Area: P<0.001 Lt ICA, P<0.001 Rt ICA, P=0.004 Lt MCA, P=0.004 Rt MCA; flow pulsatility index: P=0.042 Rt ICA, P=0.002 Lt MCA). Motion induced bias can lead to significant overestimation of hemodynamic markers in cerebral arteries. The proposed method reduces measurement bias from rigid motion in neurovascular 4D-Flow MRI in challenging populations such as aging adults.

Accelerated dual-venc 4D flow MRI with variable high-venc spatial resolution for neurovascular applications

PURPOSE

Dual-velocity encoded (dual-venc or DV) 4D flow MRI achieves wide velocity dynamic range and velocity-to-noise ratio (VNR), enabling accurate neurovascular flow characterization. To reduce scan time, we present interleaved dual-venc 4D Flow with independently prescribed, prospectively undersampled spatial resolution of the high-venc (HV) acquisition: Variable Spatial Resolution Dual Venc (VSRDV).

METHODS

A prototype VSRDV sequence was developed based on a Cartesian acquisition with eight-point phase encoding, combining PEAK-GRAPPA acceleration with zero-filling in phase and partition directions for HV. The VSRDV approach was optimized by varying z, the zero-filling fraction of HV relative to low-venc, between 0%-80% in vitro (realistic neurovascular model with pulsatile flow) and in vivo (n = 10 volunteers). Antialiasing precision, mean and peak velocity quantification accuracy, and test-retest reproducibility were assessed relative to reference images with equal-resolution HV and low venc (z = 0%).

RESULTS

In vitro results for all z demonstrated an antialiasing true positive rate at least 95% for R PEAK - GRAPPA $$ {R}_{\mathrm{PEAK}-\mathrm{GRAPPA}} $$  = 2 and 5, with no linear relationship to z (p = 0.62 and 0.13, respectively). Bland-Altman analysis for z = 20%, 40%, 60%, or 80% versus z = 0% in vitro and in vivo demonstrated no bias >1% of venc in mean or peak velocity values at any R ZF $$ {R}_{\mathrm{ZF}} $$ . In vitro mean and peak velocity, and in vivo peak velocity, had limits of agreement within 15%.

CONCLUSION

VSRDV allows up to 34.8% scan time reduction compared to PEAK-GRAPPA accelerated DV 4D Flow MRI, enabling large spatial coverage and dynamic range while maintaining VNR and velocity measurement accuracy.

Challenges in Modeling Hemodynamics in Cerebral Aneurysms Related to Arteriovenous Malformations

PURPOSE

The significantly higher incidence of aneurysms in patients with arteriovenous malformations (AVMs) suggests a strong hemodynamic relationship between these lesions. The presence of an AVM alters hemodynamics in proximal vessels by drastically changing the distal resistance, thus affecting intra-aneurysmal flow. This study discusses the challenges associated with patient-specific modeling of aneurysms in the presence of AVMs.

METHODS

We explore how the presence of a generic distal AVM affects upstream aneurysms by examining the relationship between distal resistance and aneurysmal wall shear stress using physiologically realistic estimates for the influence of the AVM on hemodynamics. Using image-based computational models of aneurysms and surrounding vasculature, aneurysmal wall-shear stress is calculated for a range of distal resistances corresponding to the presence of AVMs of various sizes and compared with a control case representing the absence of an AVM.

RESULTS

In the patient cases considered, the alteration in aneurysmal wall shear stress due to the presence of an AVM is considerable, as much as 19 times the base case wall shear stress. Furthermore, the relationship between aneurysmal wall shear stress and distal resistance is shown to be highly geometry-dependent and nonlinear. In most cases, the range of physiologically realistic possibilities for AVM-related distal resistance are so large that patient-specific flow measurements are necessary for meaningful predictions of wall shear stress.

CONCLUSIONS

The presented work offers insight on the impact of distal AVMs on aneurysmal wall shear stress using physiologically realistic computational models. Patient-specific modeling of hemodynamics in aneurysms and associated AVMs has great potential for understanding lesion pathogenesis, surgical planning, and assessing the effect of treatment of one lesion relative to another. However, we show that modeling approaches cannot usually meaningfully quantify the impact of AVMs if based solely on imaging data from CT and X-ray angiography, currently used in clinical practice. Based on recent studies, it appears that 4D flow MRI is one promising approach to obtaining meaningful patient-specific flow boundary conditions that improve modeling fidelity.

Shortened cerebral circulation time correlates with seizures in brain arteriovenous malformation: a cross-sectional quantitative digital subtraction angiography study

OBJECTIVE

Seizure is the most common clinical presentation in patients with nonhemorrhagic brain arteriovenous malformations (BAVMs) and it influences their quality of life. Angioarchitectural analysis of the seizure risk for BAVMs is subjective and does not consider hemodynamics. This study aimed to investigate the angioarchitectural and hemodynamic factors that may be associated with seizure in patients with BAVMs.

METHODS

From 2011 to 2019, 104 patients with supratentorial BAVMs without previous hemorrhage or treatment were included and grouped according to the initial presentation of seizure. Their angiograms and MRI results were analyzed for morphological characteristics and quantitative digital subtraction angiography (QDSA) parameters. Modified cerebral circulation time (mCCT) was defined as the difference between the bolus arrival time of the ipsilateral cavernous internal carotid artery and the parietal vein on lateral DSA. Logistic regression analysis was performed to estimate the odds ratio (OR) for BAVMs presenting with seizure.

RESULTS

The seizure group had shorter mCCT (1.98 s vs. 2.44 s, p = 0.005) and more BAVMs with temporal location (45% vs. 30.8%, p = 0.013), neoangiogenesis (55% vs. 33%, p = 0.03), and long draining veins (95% vs. 72%, p = 0.004) than did the nonseizure group. Shorter mCCT (OR: 3.4, p = 0.02), temporal location (OR: 13.4, p < 0.001), and neoangiogenesis (OR: 4.7, p = 0.013) were independently associated with higher risks of seizure, after adjustments for age, gender, BAVM volume, and long draining vein.

CONCLUSIONS

Shorter mCCT, temporal location, and neoangiogenesis were associated with epileptic BAVMs. QDSA can objectively evaluate hemodynamic changes in epileptic BAVMs.

KEY POINTS

• Quantitative digital subtraction angiography may be used to evaluate the hemodynamic differences between brain arteriovenous malformations presenting with and without seizure. • BAVMs with temporal location, neoangiogenesis, and shortened cerebral circulation time were more likely to present with seizure.

Quality Control for 4D Flow MR Imaging

In recent years, 4D flow MRI has become increasingly important in clinical applications for the blood vessels in the whole body, heart, and cerebrospinal fluid. 4D flow MRI has advantages over 2D cine phase-contrast (PC) MRI in that any targeted area of interest can be analyzed post-hoc, but there are some factors to be considered, such as ensuring measurement accuracy, a long imaging time and post-processing complexity, and interobserver variability.Due to the partial volume phenomenon caused by low spatial and temporal resolutions, the accuracy of flow measurement in 4D flow MRI is reduced. For spatial resolution, it is recommended to include at least four voxels in the vessel of interest, and if possible, six voxels. In large vessels such as the aorta, large voxels can be secured and SNR can be maintained, but in small cerebral vessels, SNR is reduced, resulting in reduced accuracy. A temporal resolution of less than 40 ms is recommended. The velocity-to-noise ratio (VNR) of low-velocity blood flow is low, resulting in poor measurement accuracy. The use of dual velocity encoding (VENC) or multi-VENC is recommended to avoid velocity wrap around and to increase VNR. In order to maintain sufficient spatio-temporal resolution, a longer imaging time is required, leading to potential patient movement during examination and a corresponding decrease in measurement accuracy.For the clinical application of new technologies, including various acceleration techniques, in vitro and in vivo accuracy verification based on existing accuracy-validated 2D cine PC MRI and 4D flow MRI, as well as accuracy verification on the conservation of mass' principle, should be performed, and intraobserver repeatability, interobserver reproducibility, and test-retest reproducibility should be checked.

Influence of Spatial Resolution and Compressed SENSE Acceleration Factor on Flow Quantification with 4D Flow MRI at 3 Tesla

Four-dimensional (4D) flow MRI allows quantifying flow in blood vessels–non invasively and in vivo. The clinical use of 4D flow MRI in small vessels, however, is hampered by long examination times and limited spatial resolution. Compressed SENSE (CS-SENSE) is a technique that can accelerate 4D flow dramatically. Here, we investigated the effect of spatial resolution and CS acceleration on flow measurements by using 4D flow MRI in small vessels in vitro at 3 T. We compared the flow in silicon tubes (inner diameters of 2, 3, 4, and 5 mm) measured with 4D flow MRI, accelerated with four CS factors (CS = 2.5, 4.5, 6.5, and 13) and three voxel sizes (0.5, 1, and 1.5 mm³) to 2D flow MRI and a flow sensor. Additionally, the velocity field in an aneurysm model acquired with 4D flow MRI was compared to the one simulated with computational fluid dynamics (CFD). A strong correlation was observed between flow sensor, 2D flow MRI, and 4D flow MRI (rho > 0.94). The use of fewer than seven voxels per vessel diameter (nROI) resulted in an overestimation of flow in more than 5% of flow measured with 2D flow MRI. A negative correlation (rho = −0.81) between flow error and nROI were found for CS = 2.5 and 4.5. No statistically significant impact of CS factor on differences in flow rates was observed. However, a trend of increased flow error with increased CS factor was observed. In an aneurysm model, the peak velocity and stagnation zone were detected by CFD and all 4D flow MRI variants. The velocity difference error in the aneurysm sac did not exceed 11% for CS = 4.5 in comparison to CS = 2.5 for all spatial resolutions. Therefore, CS factors from 2.5–4.5 can appear suitable to improve spatial or temporal resolution for accurate quantification of flow rate and velocity. We encourage reporting the number of voxels per vessel diameter to standardize 4D flow MRI protocols.

4D flow MRI hemodynamic biomarkers for cerebrovascular diseases

Alterations in cerebral blood flow are common in several neurological diseases among the elderly including stroke, cerebral small vessel disease, vascular dementia, and Alzheimer's disease. 4D flow magnetic resonance imaging (MRI) is a relatively new technique to investigate cerebrovascular disease, and makes it possible to obtain time-resolved blood flow measurements of the entire cerebral arterial venous vasculature and can be used to derive a repertoire of hemodynamic biomarkers indicative of cerebrovascular health. The information that can be obtained from one single 4D flow MRI scan allows both the investigation of aberrant flow patterns at a focal location in the vasculature as well as estimations of brain-wide disturbances in blood flow. Such focal and global hemodynamic biomarkers show the potential of being sensitive to impending cerebrovascular disease and disease progression and can also become useful during planning and follow-up of interventions aiming to restore a normal cerebral circulation. Here, we describe 4D flow MRI approaches for analyzing the cerebral vasculature. We then survey key hemodynamic biomarkers that can be reliably assessed using the technique. Finally, we highlight cerebrovascular diseases where one or multiple hemodynamic biomarkers are of central interest.

Noninvasive quantification of cerebrovascular pressure changes using 4D Flow MRI

PURPOSE

Hemodynamic alterations are indicative of cerebrovascular disease. However, the narrow and tortuous cerebrovasculature complicates image-based assessment, especially when quantifying relative pressure. Here, we present a systematic evaluation of image-based cerebrovascular relative pressure mapping, investigating the accuracy of the routinely used reduced Bernoulli (RB), the extended unsteady Bernoulli (UB), and the full-field virtual work-energy relative pressure ( ν WERP) method.

METHODS

Patient-specific in silico models were used to generate synthetic cerebrovascular 4D Flow MRI, with RB, UB, and ν WERP performance quantified as a function of spatiotemporal sampling and image noise. Cerebrovascular relative pressures were also derived in 4D Flow MRI from healthy volunteers ( n = 8 ), acquired at two spatial resolutions (dx = 1.1 and 0.8 mm).

RESULTS

The in silico analysis indicate that accurate relative pressure estimations are inherently coupled to spatial sampling: at dx = 1.0 mm high errors are reported for all methods; at dx = 0.5 mm ν WERP recovers relative pressures at a mean error of 0.02 ± 0.25 mm Hg, while errors remain higher for RB and UB (mean error of -2.18 ± 1.91 and -2.18 ± 1.87 mm Hg, respectively). The dependence on spatial sampling is also indicated in vivo, albeit with higher correlative dependence between resolutions using ν WERP (k = 0.64, R2 = 0.81 for dx = 1.1 vs. 0.8 mm) than with RB or UB (k = 0.04, R2 = 0.03, and k = 0.07, R2 = 0.07, respectively).

CONCLUSION

Image-based full-field methods such as ν WERP enable cerebrovascular relative pressure mapping; however, accuracy is directly dependent on utilized spatial resolution.

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.

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

BACKGROUND

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

PURPOSE

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

STUDY TYPE

Retrospective.

POPULATION

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

FIELD STRENGTH/SEQUENCE

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

ASSESSMENT

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

STATISTICAL TESTS

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

RESULTS

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

DATA CONCLUSION

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

LEVEL OF EVIDENCE

3 TECHNICAL EFFICACY: Stage 2.

Hemodynamic Analysis of Cerebral AVMs with 3D Phase-Contrast MR Imaging

BACKGROUND AND PURPOSE

The hemodynamics associated with cerebral AVMs have a significant impact on their clinical presentation. This study aimed to evaluate the hemodynamic features of AVMs using 3D phase-contrast MR imaging with dual velocity-encodings.

MATERIALS AND METHODS

Thirty-two patients with supratentorial AVMs who had not received any previous treatment and had undergone 3D phase-contrast MR imaging were included in this study. The nidus diameter and volume were measured for classification of AVMs (small, medium, or large). Flow parameters measured included apparent AVM inflow, AVM inflow index, apparent AVM outflow, AVM outflow index, and the apparent AVM inflow-to-outflow ratio. Correlation coefficients between the nidus volume and each flow were calculated. The flow parameters between small and other AVMs as well as between nonhemorrhagic and hemorrhagic AVMs were compared.

RESULTS

Patients were divided into hemorrhagic (n = 8) and nonhemorrhagic (n = 24) groups. The correlation coefficient between the nidus volume and the apparent AVM inflow and outflow was .83. The apparent AVM inflow and outflow in small AVMs were significantly smaller than in medium AVMs (P < .001 for both groups). The apparent AVM inflow-to-outflow ratio was significantly larger in the hemorrhagic AVMs than in the nonhemorrhagic AVMs (P = .02).

CONCLUSIONS

The apparent AVM inflow-to-outflow ratio was the only significant parameter that differed between nonhemorrhagic and hemorrhagic AVMs, suggesting that a poor drainage system may increase AVM pressure, potentially causing cerebral hemorrhage.

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.

Quantification of hemodynamics of cerebral arteriovenous malformations after stereotactic radiosurgery using 4D flow magnetic resonance imaging

Stereotactic radiosurgery (SRS) is used to treat cerebral arteriovenous malformations (AVMs). However, early evaluation of efficacy is difficult as structural magnetic resonance imaging (MRI)/magnetic resonance angiography (MRA) often does not demonstrate appreciable changes within the first 6 months. The aim of this study was to evaluate the use of four-dimensional (4D) flow MRI to quantify hemodynamic changes after SRS as early as 2 months. This was a retrospective observational study, which included 14 patients with both pre-SRS and post-SRS imaging obtained at multiple time points from 1 to 27 months after SRS. A 3T MRI Scanner was used to obtain T₂ single-shot fast spin echo, time-of-flight MRA, and postcontrast 4D flow with three-dimensional velocity encoding between 150 and 200 cm/s. Post-hoc two-dimensional cross-sectional flow was measured for the dominant feeding artery, the draining vein, and the corresponding contralateral artery as a control. Measurements were performed by two independent observers, and reproducibility was assessed. Wilcoxon signed-rank tests were used to compare differences in flow, circumference, and pulsatility between the feeding artery and the contralateral artery both before and after SRS; and differences in nidus size and flow and circumference of the feeding artery and draining vein before and after SRS. Arterial flow (L/min) decreased in the primary feeding artery (mean: 0.1 ± 0.07 vs. 0.3 ± 0.2; p < 0.05) and normalized in comparison to the contralateral artery (mean: 0.1 ± 0.07 vs. 0.1 ± 0.07; p = 0.068). Flow decreased in the draining vein (mean: 0.1 ± 0.2 vs. 0.2 ± 0.2; p < 0.05), and the circumference of the draining vein also decreased (mean: 16.1 ± 8.3 vs. 15.7 ± 6.7; p < 0.05). AVM volume decreased after SRS (mean: 45.3 ± 84.8 vs. 38.1 ± 78.7; p < 0.05). However, circumference (mm) of the primary feeding artery remained similar after SRS (mean: 15.7 ± 2.7 vs. 16.1 ± 3.1; p = 0.600). 4D flow may be able to demonstrate early hemodynamic changes in AVMs treated with radiosurgery, and these changes appear to be more pronounced and occur earlier than the structural changes on standard MRI/MRA. Level of Evidence: 4 Technical Efficacy Stage: 1.

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.

4D Flow with MRI

Magnetic resonance imaging (MRI) has become an important tool for the clinical evaluation of patients with cardiac and vascular diseases. Since its introduction in the late 1980s, quantitative flow imaging with MRI has become a routine part of standard-of-care cardiothoracic and vascular MRI for the assessment of pathological changes in blood flow in patients with cardiovascular disease. More recently, time-resolved flow imaging with velocity encoding along all three flow directions and three-dimensional (3D) anatomic coverage (4D flow MRI) has been developed and applied to enable comprehensive 3D visualization and quantification of hemodynamics throughout the human circulatory system. This article provides an overview of the use of 4D flow applications in different cardiac and vascular regions in the human circulatory system, with a focus on using 4D flow MRI in cardiothoracic and cerebrovascular diseases.