Improved SNR in phase contrast velocimetry with five-point balanced flow encoding

CSF dynamics throughout the ventricular system using 4D flow MRI: associations to arterial pulsatility, ventricular volumes, and age

BACKGROUND

Cerebrospinal fluid (CSF) dynamics are increasingly studied in aging and neurological disorders. Models of CSF-mediated waste clearance suggest that altered CSF dynamics could play a role in the accumulation of toxic waste in the CNS, with implications for Alzheimer's disease and other proteinopathies. Therefore, approaches that enable quantitative and volumetric assessment of CSF flow velocities could be of value. In this study we demonstrate the feasibility of 4D flow MRI for simultaneous assessment of CSF dynamics throughout the ventricular system, and evaluate associations to arterial pulsatility, ventricular volumes, and age.

METHODS

In a cognitively unimpaired cohort (N = 43; age 41-83 years), cardiac-resolved 4D flow MRI CSF velocities were obtained in the lateral ventricles (LV), foramens of Monro, third and fourth ventricles (V3 and V4), the cerebral aqueduct (CA) and the spinal canal (SC), using a velocity encoding (venc) of 5 cm/s. Cerebral blood flow pulsatility was also assessed with 4D flow (venc = 80 cm/s), and CSF volumes were obtained from T1- and T2-weighted MRI. Multiple linear regression was used to assess effects of age, ventricular volumes, and arterial pulsatility on CSF velocities.

RESULTS

Cardiac-driven CSF dynamics were observed in all CSF spaces, with region-averaged velocity range and root-mean-square (RMS) velocity encompassing from very low in the LVs (RMS 0.25 ± 0.08; range 0.85 ± 0.28 mm/s) to relatively high in the CA (RMS 6.29 ± 2.87; range 18.6 ± 15.2 mm/s). In the regression models, CSF velocity was significantly related to age in 5/6 regions, to CSF space volume in 2/3 regions, and to arterial pulsatility in 3/6 regions. Group-averaged waveforms indicated distinct CSF flow propagation delays throughout CSF spaces, particularly between the SC and LVs.

CONCLUSIONS

Our findings show that 4D flow MRI enables assessment of CSF dynamics throughout the ventricular system, and captures independent effects of age, CSF space morphology, and arterial pulsatility on CSF motion.

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.

Blood-brain barrier integrity is linked to cognitive function, but not to cerebral arterial pulsatility, among elderly

Blood-brain barrier (BBB) disruption may contribute to cognitive decline, but questions remain whether this association is more pronounced for certain brain regions, such as the hippocampus, or represents a whole-brain mechanism. Further, whether human BBB leakage is triggered by excessive vascular pulsatility, as suggested by animal studies, remains unknown. In a prospective cohort (N = 50; 68-84 years), we used contrast-enhanced MRI to estimate the permeability-surface area product (PS) and fractional plasma volume ( v p ), and 4D flow MRI to assess cerebral arterial pulsatility. Cognition was assessed by the Montreal Cognitive Assessment (MoCA) score. We hypothesized that high PS would be associated with high arterial pulsatility, and that links to cognition would be specific to hippocampal PS. For 15 brain regions, PS ranged from 0.38 to 0.85 (·10-3 min-1) and v p from 0.79 to 1.78%. Cognition was related to PS (·10-3 min-1) in hippocampus (β = - 2.9; p = 0.006), basal ganglia (β = - 2.3; p = 0.04), white matter (β = - 2.6; p = 0.04), whole-brain (β = - 2.7; p = 0.04) and borderline-related for cortex (β = - 2.7; p = 0.076). Pulsatility was unrelated to PS for all regions (p > 0.19). Our findings suggest PS-cognition links mainly reflect a whole-brain phenomenon with only slightly more pronounced links for the hippocampus, and provide no evidence of excessive pulsatility as a trigger of BBB disruption.

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.

Establishing the distribution of cerebrovascular resistance using computational fluid dynamics and 4D flow MRI

Cerebrovascular resistance (CVR) regulates blood flow in the brain, but little is known about the vascular resistances of the individual cerebral territories. We present a method to calculate these resistances and investigate how CVR varies in the hemodynamically disturbed brain. We included 48 patients with stroke/TIA (29 with symptomatic carotid stenosis). By combining flow rate (4D flow MRI) and structural computed tomography angiography (CTA) data with computational fluid dynamics (CFD) we computed the perfusion pressures out from the circle of Willis, with which CVR of the MCA, ACA, and PCA territories was estimated. 56 controls were included for comparison of total CVR (tCVR). CVR were 33.8 ± 10.5, 59.0 ± 30.6, and 77.8 ± 21.3 mmHg s/ml for the MCA, ACA, and PCA territories. We found no differences in tCVR between patients, 9.3 ± 1.9 mmHg s/ml, and controls, 9.3 ± 2.0 mmHg s/ml (p = 0.88), nor in territorial CVR in the carotid stenosis patients between ipsilateral and contralateral hemispheres. Territorial resistance associated inversely to territorial brain volume (p < 0.001). These resistances may work as reference values when modelling blood flow in the circle of Willis, and the method can be used when there is need for subject-specific analysis.

Phase-contrast magnetic resonance imaging of intracranial and extracranial blood flow in carotid near-occlusion

PURPOSE

Compare extracranial internal carotid artery flow rates and intracranial collateral use between conventional ≥ 50% carotid stenosis and carotid near-occlusion, and between symptomatic and asymptomatic carotid near-occlusion.

METHODS

We included patients with ≥ 50% carotid stenosis. Degree of stenosis was diagnosed on CTA. Mean blood flow rates were assessed with four-dimensional phase-contrast MRI.

RESULTS

We included 110 patients of which 83% were symptomatic, and 38% had near-occlusion. Near-occlusions had lower mean internal carotid artery flow (70 ml/min) than conventional ≥ 50% stenoses (203 ml/min, P < .001). Definite use of ≥ 1 collateral was found in 83% (35/42) of near-occlusions and 10% (7/68) of conventional stenoses (P < .001). However, there were no differences in total cerebral blood flow (514 ml/min vs. 519 ml/min, P = .78) or ipsilateral hemispheric blood flow (234 vs. 227 ml/min, P = .52), between near-occlusions and conventional ≥ 50% stenoses, based on phase-contrast MRI flow rates. There were no differences in total cerebral or hemispheric blood flow, or collateral use, between symptomatic and asymptomatic near-occlusions.

CONCLUSION

Near-occlusions have lower internal carotid artery flow rates and more collateral use, but similar total cerebral blood flow and hemispheric blood flow, compared to conventional ≥ 50% carotid stenosis.

Diagnostic separation of conventional ⩾50% carotid stenosis and near-occlusion with phase-contrast MRI

INTRODUCTION

The aim of this study was to assess sensitivity, specificity and interrater reliability of phase-contrast MRI (PC-MRI) for diagnosing carotid near-occlusion.

PATIENTS AND METHODS

Prospective cross-sectional study conducted between 2018 and 2021. We included participants with suspected 50%-100% carotid stenosis on at least one side, all were examined with CT angiography (CTA) and PC-MRI and both ICAs were analyzed. Degree of stenosis on CTA was the reference test. PC-MRI-based blood flow rates in extracranial ICA and intracranial cerebral arteries were assessed. ICA-cerebral blood flow (CBF) ratio was defined as ICA divided by sum of both ICAs and Basilar artery.

RESULTS

We included 136 participants. The ICAs were 102 < 50% stenosis, 88 conventional ⩾50% stenosis (31 with ⩾70%), 49 near-occlusion, 12 occlusions, 20 unclear cause of small distal ICA on CTA and one excluded. For separation of near-occlusion and conventional stenoses, ICA flow rate and ICA-CBF ratio had the highest area under the curve (AUC; 0.98-0.99) for near-occlusion. ICA-CBF ratio ⩽0.225 was 90% (45/49) sensitive and 99% (188/190) specific for near-occlusion. Inter-rater reliability for this threshold was excellent (kappa 0.98). Specificity was 94% (29/31) for cases with ⩾70% stenosis. PC-MRI had modest performance for separating <50% and conventional ⩾50% stenosis (highest AUC 0.74), and eight (16%) of near-occlusions were not distinguishable from occlusion (no visible flow).

CONCLUSION

ICA-CBF ratio ⩽0.225 on PC-MRI is an accurate and reliable method to separate conventional ⩾50% stenosis and near-occlusion that is feasible for routine use. PC-MRI should be considered further as a potential standard method for near-occlusion detection, to be used side-by-side with established modalities as PC-MRI cannot separate other degrees of stenosis.

Fat mitigation strategies to improve image quality of radial 4D flow MRI in obese subjects

PURPOSE

This study addresses the challenges in obtaining abdominal 4D flow MRI of obese patients. We aimed to evaluate spectral saturation and inner volume excitation as methods to mitigating artifacts originating from adipose signals, with the goal of enhancing image quality and improving quantification.

METHODS

Radial 4D flow MRI acquisitions with fat mitigation (inner volume excitation [IVE] and intermittent fat saturation [FS]) were compared to a standard slab selective excitation (SSE) in a test-retest study of 15 obese participants. IVE selectively excited a cylindrical region of interest, avoiding contamination from peripheral adipose tissue, while FS globally suppressed fat based on spectral selection. Acquisitions were evaluated qualitatively based on expert ratings and quantitatively based on conservation of mass, test-retest repeatability, and a divergence free quality metric. Errors were evaluated statistically using the absolute and relative errors, regression, and Bland-Altman analysis.

RESULTS

IVE demonstrated superior performance quantitatively in the conservation of mass analysis in the portal vein, with higher correlation and lower bias in regression analysis. IVE also produced flow fields with the lowest divergence error and was rated best in overall image quality, delineating small vessels, and producing the least streaking artifacts. Evaluation results did not differ significantly between FS and SSE. Test-retest reproducibility was similarly high for all sequences, with data suggesting biological variations dominate the technical variability.

CONCLUSION

IVE improved hemodynamic assessment of radial 4D flow MRI in the abdomen of obese participants while FS did not lead to significant improvements in image quality or flow metrics.

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.

Focused navigation for respiratory-motion-corrected free-running radial 4D flow MRI

PURPOSE

To validate a respiratory motion correction method called focused navigation (fNAV) for free-running radial whole-heart 4D flow MRI.

METHODS

Using fNAV, respiratory signals derived from radial readouts are converted into three orthogonal displacements, which are then used to correct respiratory motion in 4D flow datasets. Hundred 4D flow acquisitions were simulated with non-rigid respiratory motion and used for validation. The difference between generated and fNAV displacement coefficients was calculated. Vessel area and flow measurements from 4D flow reconstructions with (fNAV) and without (uncorrected) motion correction were compared to the motion-free ground-truth. In 25 patients, the same measurements were compared between fNAV 4D flow, 2D flow, navigator-gated Cartesian 4D flow, and uncorrected 4D flow datasets.

RESULTS

For simulated data, the average difference between generated and fNAV displacement coefficients was 0.04± $$ \pm $$ 0.32 mm and 0.31± $$ \pm $$ 0.35 mm in the x and y directions, respectively. In the z direction, this difference was region-dependent (0.02± $$ \pm $$ 0.51 mm up to 5.85± $$ \pm $$ 3.41 mm). For all measurements (vessel area, net volume, and peak flow), the average difference from ground truth was higher for uncorrected 4D flow datasets (0.32± $$ \pm $$ 0.11 cm2 , 11.1± $$ \pm $$ 3.5 mL, and 22.3± $$ \pm $$ 6.0 mL/s) than for fNAV 4D flow datasets (0.10± $$ \pm $$ 0.03 cm2 , 2.6± $$ \pm $$ 0.7 mL, and 5.1± 0 $$ \pm 0 $$ .9 mL/s, p < 0.05). In vivo, average vessel area measurements were 4.92± $$ \pm $$ 2.95 cm2 , 5.06± $$ \pm $$ 2.64 cm2 , 4.87± $$ \pm $$ 2.57 cm2 , 4.87± $$ \pm $$ 2.69 cm2 , for 2D flow and fNAV, navigator-gated and uncorrected 4D flow datasets, respectively. In the ascending aorta, all 4D flow datasets except for the fNAV reconstruction had significantly different vessel area measurements from 2D flow. Overall, 2D flow datasets demonstrated the strongest correlation to fNAV 4D flow for both net volume (r2  = 0.92) and peak flow (r2  = 0.94), followed by navigator-gated 4D flow (r2  = 0.83 and r2  = 0.86, respectively), and uncorrected 4D flow (r2  = 0.69 and r2  = 0.86, respectively).

CONCLUSION

fNAV corrected respiratory motion in vitro and in vivo, resulting in fNAV 4D flow measurements that are comparable to those derived from 2D flow and navigator-gated Cartesian 4D flow datasets, with improvements over those from uncorrected 4D flow.

Asynchronous magnetic resonance elastography: Shear wave speed reconstruction using noise correlation of incoherent waves

PURPOSE

The noninvasive measurement of biological tissue elasticity is an evolving technology that enables the robust characterization of soft tissue mechanics for a wide array of biomedical engineering and clinical applications. We propose, design, and implement here a new MRI technique termed asynchronous magnetic resonance elastography (aMRE) that pushes the measurement technology toward a driverless implementation. This technique can be added to clinical MRI scanners without any additional specialized hardware.

THEORY

Asynchronous MRE is founded on the theory of diffuse wavefields and noise correlation previously developed in ultrasound to reconstruct shear wave speeds using seemingly incoherent wavefields. Unlike conventional elastography methods that solve an inverse problem, aMRE directly reconstructs a pixel-wise mapping of wave speed using the spatial-temporal statistics of the measured wavefield.

METHODS

Incoherent finger tapping served as the wave-generating source for all aMRE measurements. Asynchronous MRE was performed on a phantom using a Siemens Prismafit as an experimental validation of the theory. It was further performed on thigh muscles as a proof-of-concept implementation of in vivo imaging using a Siemens Skyra scanner.

RESULTS

Numerical and phantom experiments show an accurate reconstruction of wave speeds from seemingly noisy wavefields. The proof-of-concept thigh experiments also show that the aMRE protocol can reconstruct a pixel-wise mapping of wave speeds.

CONCLUSION

Asynchronous MRE is shown to accurately reconstruct shear wave speeds in phantom experiments and remains at the proof-of-concept stage for in vivo imaging. After further validation and improvements, it has the potential to lower both the technical and monetary barriers of entry to measuring tissue elasticity.

Comparison of Improved Unidirectional Dual Velocity-Encoding MRI Methods

BACKGROUND

In phase-contrast (PC) MRI, several dual velocity encoding methods have been proposed to robustly increase velocity-to-noise ratio (VNR), including a standard dual-VENC (SDV), an optimal dual-VENC (ODV), and bi- and triconditional methods.

PURPOSE

To develop a correction method for the ODV approach and to perform a comparison between methods.

STUDY TYPE

Case-control study.

POPULATION

Twenty-six volunteers.

FIELD STRENGTH/SEQUENCE

1.5 T phase-contrast MRI with VENCs of 50, 75, and 150 cm/second.

ASSESSMENT

Since we acquired single-VENC protocols, we used the background phase from high-VENC (VENC_H ) to reconstruct the low-VENC (VENC_L ) phase. We implemented and compared the unwrapping methods for different noise levels and also developed a correction of the ODV method.

STATISTICAL TESTS

Shapiro-Wilk's normality test, two-way analysis of variance with homogeneity of variances was performed using Levene's test, and the significance level was adjusted by Tukey's multiple post hoc analysis with Bonferroni (P < 0.05).

RESULTS

Statistical analysis revealed no extreme outliers, normally distributed residuals, and homogeneous variances. We found statistically significant interaction between noise levels and the unwrapping methods. This implies that the number of non-unwrapped pixels increased with the noise level. We found that for β = VENC_L /VENC_H  = 1/2, unwrapping methods were more robust to noise. The post hoc test showed a significant difference between the ODV corrected and the other methods, offering the best results regarding the number of unwrapped pixels.

DATA CONCLUSIONS

All methods performed similarly without noise, but the ODV corrected method was more robust to noise at the price of a higher computational time.

LEVEL OF EVIDENCE

4 TECHNICAL EFFICACY STAGE: 1.

CSF Flow and Spinal Cord Motion in Patients With Spontaneous Intracranial Hypotension: A Phase Contrast MRI Study

BACKGROUND AND OBJECTIVES

Spontaneous intracranial hypotension (SIH) is characterized by loss of CSF volume. We hypothesize that in this situation of low volume, a larger CSF flow and spinal cord motion at the upper spine can be measured by noninvasive phase contrast MRI.

METHODS

A prospective, age-, sex-, and body mass index (BMI)-matched controlled cohort study on patients with SIH presenting with spinal longitudinal extradural fluid collection (SLEC) was conducted from October 2021 to February 2022. Cardiac-gated 2D phase contrast MRI sequences were acquired at segment C2/C3, and C5/C6 for CSF flow, and spinal cord motion analysis. Data processing was fully automated. CSF flow and spinal cord motion were analyzed by peak-to-peak amplitude and total displacement per segment and heartbeat, respectively. Clinical data included age, height, BMI, duration of symptoms, Bern score according to Dobrocky et al., and type of the spinal CSF leak according to Schievink et al. Groups were compared via the Mann-Whitney U test; multiple linear regression analysis was performed to address possible relations.

RESULTS

Twenty patients with SIH and 40 healthy controls were analyzed; each group consisted of 70% women. Eleven patients with SIH presented with type 1 leak, 8 with type 2, and 1 was indeterminate. CSF flow per heartbeat was increased at C2/C3 (peak-to-peak amplitude 65.68 ± 18.3 vs 42.50 ± 9.8 mm/s, total displacement 14.32 ± 3.5 vs 9.75 ± 2.7 mm, p < 0.001, respectively). Craniocaudal spinal cord motion per heartbeat was larger at segment C2/C3 (peak-to-peak amplitude 7.30 ± 2.4 vs 5.82 ± 2.0 mm/s, total displacement 1.01 ± 0.4 vs 0.74 ± 0.4 mm, p = 0.006, respectively) and at segment C5/C6 (total displacement 1.41 ± 0.7 vs 0.97 ± 0.4 mm, p = 0.021).

DISCUSSION

SLEC-positive patients with SIH show higher CSF flow and higher spinal cord motion at the upper cervical spine. This increased craniocaudal motion of the spinal cord per heartbeat might produce increased mechanical strain on neural tissue and adherent structures, which may be a mechanism leading to cranial nerve dysfunction, neck pain, and stiffness in SIH. Noninvasive phase contrast MRI of CSF flow and spinal cord motion is a promising diagnostic tool in SIH.

TRIAL REGISTRATION INFORMATION

German Clinical Trials Register, identification number: DRKS00017351.

CLASSIFICATION OF EVIDENCE

This study provides Class III evidence that noninvasive phase contrast MRI of the upper spine identifies differences in CSF flow and spinal cord motion in patients with SIH compared with healthy controls.

Accelerated 4D-flow MRI with 3-point encoding enabled by machine learning

PURPOSE

To investigate the acceleration of 4D-flow MRI using a convolutional neural network (CNN) that produces three directional velocities from three flow encodings, without requiring a fourth reference scan measuring background phase.

METHODS

A fully 3D CNN using a U-net architecture was trained in a block-wise fashion to take complex images from three flow encodings and to produce three real-valued images for each velocity component. Using neurovascular 4D-flow scans (n = 144), the CNN was trained to predict velocities computed from four flow encodings by standard reconstruction including correction for residual background phase offsets. Methods to optimize loss functions were investigated, including magnitude, complex difference, and uniform velocity weightings. Subsequently, 3-point encoding was evaluated using cross validation of pixelwise correlation, flow measurements in major arteries, and in experiments with data at differing acceleration rates than the training data.

RESULTS

The CNN-produced 3-point velocities showed excellent agreements with the 4-point velocities, both qualitatively in velocity images, in flow rate measures, and quantitatively in regression analysis (slope = 0.96, R²  = 0.992). Optimizing the training to focus on vessel velocities rather than the global velocity error and improved the correlation of velocity within vessels themselves. The lowest error was observed when the loss function used uniform velocity weighting, in which the magnitude-weighted inverse of the velocity frequency uniformly distributed weighting across all velocity ranges. When applied to highly accelerated data, the 3-point network maintained a high correlation with ground truth data and demonstrated a denoising effect.

CONCLUSION

The 4D-flow MRI can be accelerated using machine learning requiring only three flow encodings to produce three-directional velocity maps with small errors.

Spinal cord motion assessed by phase-contrast MRI - An inter-center pooled data analysis

BACKGROUND

Phase-contrast MRI of CSF and spinal cord dynamics has evolved among diseases caused by altered CSF volume (spontaneous intracranial hypotension, normal pressure hydrocephalus) and by altered CSF space (degenerative cervical myelopathy (DCM), Chiari malformation). While CSF seems to be an obvious target for possible diagnostic use, craniocaudal spinal cord motion analysis offers the benefit of fast and reliable assessments. It is driven by volume shifts between the intracranial and the intraspinal compartments (Monro-Kellie hypothesis). Despite promising initial reports, comparison of spinal cord motion data across different centers is challenged by reports of varying value, raising questions about the validity of the findings.

OBJECTIVE

To systematically investigate inter-center differences between phase-contrast MRI data.

METHODS

Age- and gender matched, retrospective, pooled-data analysis across two centers: cardiac-gated, sagittal phase-contrast MRI of the cervical spinal cord (segments C2/C3 to C7/T1) including healthy participants and DCM patients; comparison and analysis of different MRI sequences and processing techniques (manual versus fully automated).

RESULTS

A genuine craniocaudal spinal cord motion pattern and an increased focal spinal cord motion among DCM patients were depicted by both MRI sequences (p < 0.01). Higher time-resolution resolved steeper and larger peaks, causing inter-center differences (p < 0.01). Comparison of different processing methods showed a high level of rating reliability (ICC > 0.86 at segments C2/C3 to C6/C7).

DISCUSSION

Craniocaudal spinal cord motion is a genuine finding. Differences between values were attributed to time-resolution of the MRI sequences. Automated processing confers the benefit of unbiased and consistent analysis, while data did not reveal any superiority.

Longitudinal stability in working memory and frontal activity in relation to general brain maintenance

Cognitive functions are well-preserved for some older individuals, but the underlying brain mechanisms remain disputed. Here, 5-year longitudinal 3-back in-scanner and offline data classified individuals in a healthy older sample (baseline age = 64-68 years) into having stable or declining working-memory (WM). Consistent with a vital role of the prefrontal cortex (PFC), WM stability or decline was related to maintained or reduced longitudinal PFC functional responses. Subsequent analyses of imaging markers of general brain maintenance revealed higher levels in the stable WM group on measures of neurotransmission and vascular health. Also, categorical and continuous analyses showed that rate of WM decline was related to global (ventricles) and local (hippocampus) measures of neuronal integrity. Thus, our findings support a role of the PFC as well as general brain maintenance in explaining heterogeneity in longitudinal WM trajectories in aging.

Prediction of cerebral perfusion pressure during carotid surgery - A computational fluid dynamics approach

BACKGROUND

Maintaining cerebral perfusion pressure in the brain when a carotid artery is closed during vascular surgery is critical for avoiding intraoperative hypoperfusion and risk of ischemic stroke. Here we propose and evaluate a method based on computational fluid dynamics for predicting patient-specific cerebral perfusion pressures at carotid clamping during carotid endarterectomy.

METHODS

The study consisted of 22 patients with symptomatic carotid stenosis who underwent carotid endarterectomy (73 ± 5 years, 59-80 years, 17 men). The geometry of the circle of Willis was obtained preoperatively from computed tomography angiography and corresponding flow rates from four-dimensional flow magnetic resonance imaging. The patients were also classified as having a present or absent ipsilateral posterior communicating artery based on computed tomography angiography. The predicted mean stump pressures from computational fluid dynamics were compared with intraoperatively measured stump pressures from carotid endarterectomy.

FINDINGS

On group level, there was no difference between the predicted and measured stump pressures (-0.5 ± 13 mmHg, P = 0.86) and the pressures were correlated (r = 0.44, P = 0.039). Omitting two outliers, the correlation increased to r = 0.78 (P < 0.001) (-1.4 ± 8.0 mmHg, P = 0.45). Patients with a present ipsilateral posterior communicating artery (n = 8) had a higher measured stump pressure than those with an absent artery (n = 12) (P < 0.001).

INTERPRETATION

The stump pressure agreement indicates that the computational fluid dynamics approach was promising in predicting cerebral perfusion pressures during carotid clamping, which may prove useful in the preoperative planning of vascular interventions.

Alias-Free Arrays

Nonuniform array geometries provide freedom for increased aperture and reduced mutual coupling. A necessary and sufficient condition is given for an array of isotropic sensor elements to be unambiguous for any specified set of directions of arrival. The set of unambiguous spatial frequencies is shown to be a parallelepiped, admitting simple geometrical interpretation. Results are used in design of linear, planar, and 3D arrays.

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.

SRflow: Deep learning based super-resolution of 4D-flow MRI data

Exploiting 4D-flow magnetic resonance imaging (MRI) data to quantify hemodynamics requires an adequate spatio-temporal vector field resolution at a low noise level. To address this challenge, we provide a learned solution to super-resolve in vivo 4D-flow MRI data at a post-processing level. We propose a deep convolutional neural network (CNN) that learns the inter-scale relationship of the velocity vector map and leverages an efficient residual learning scheme to make it computationally feasible. A novel, direction-sensitive, and robust loss function is crucial to learning vector-field data. We present a detailed comparative study between the proposed super-resolution and the conventional cubic B-spline based vector-field super-resolution. Our method improves the peak-velocity to noise ratio of the flow field by 10 and 30% for in vivo cardiovascular and cerebrovascular data, respectively, for 4 × super-resolution over the state-of-the-art cubic B-spline. Significantly, our method offers 10x faster inference over the cubic B-spline. The proposed approach for super-resolution of 4D-flow data would potentially improve the subsequent calculation of hemodynamic quantities.

Characterization of mesenteric and portal hemodynamics using 4D flow MRI: the effects of meals and diurnal variation

PURPOSE

To determine the variability of blood flow measurements using 4D flow MRI in the portal and mesenteric circulations and to characterize the effects of meal ingestion, time of day, and between-day (diurnal) variations on portal and mesenteric hemodynamics.

METHODS

In this IRB-approved and HIPAA-compliant study, 7 healthy and 7 portal hypertension patients imaged. MRI exams were conducted at 3 T using a 32-channel body coil with large volumetric coverage and 1.25-mm isotropic true spatial resolution. Blood flow was quantified (L/min) in the hepatic and splanchnic vasculature. The first MR scan was performed after at least 8 h of fasting. Subsequently, subjects ingested 574 mL EnSure Plus® orally. A second acquisition was started 20 min after the meal ingestion. A third scan was performed before lunch and a fourth acquisition took place 20 min after lunch. A fifth scan was performed around 4 pm. Finally, subjects returned one week later for a repeat morning visit, with identical conditions as the first visit.

RESULTS

In healthy controls significant increase in blood flow was seen in the PV, SMV, SMA, HA, and SCAo in response to breakfast but only the SCAo, SMA, SMV, and PV had a significant response to lunch. In general, patients with cirrhosis showed reduced response to meals compared to that in healthy controls. Additionally, PV flow in patients had the highest value in the afternoon.

CONCLUSION

Effects of meal ingestion, time of day, and between-day variations were characterized using Radial 4D flow MRI in patients with cirrhosis and healthy controls.

Non-invasive assessment of mesenteric hemodynamics in patients with suspected chronic mesenteric ischemia using 4D flow MRI

PURPOSE

Chronic mesenteric ischemia (CMI) is a rare disease with a particularly difficult diagnosis. In this study, 4D flow MRI is used to quantitatively evaluate mesenteric hemodynamics before and after a meal in patients suspected of having CMI and healthy individuals.

METHODS

Nineteen patients suspected of CMI and twenty control subjects were analyzed. Subjects were scanned using a radially undersampled 4D flow MR sequence (PC-VIPR). Flow rates were assessed in the supraceliac (SCAo) and infrarenal aorta, celiac artery, superior mesenteric artery (SMA), left and right renal arteries, superior mesenteric vein (SMV), splenic vein, and portal vein (PV) in a fasting state (preprandial) and 20 min after a 700-kcal meal (postprandial). Patients were subcategorized into positive diagnosis (CMI+, N = 6) and negative diagnosis (CMI-, N = 13) groups based on imaging and clinical findings. Preprandial, postprandial, and percent change in flow rates were compared between subgroups using a Welch t test.

RESULTS

In controls and CMI- patients, SCAo, SMA, SMV, and PV flow increased significantly after meal ingestion. No significant flow increases were observed in CMI+ patients. Percent changes in SMA, SMV, and PV flow were significantly greater in controls compared to CMI+ patients. Additionally, percent changes in flow in the SMV and PV were significantly greater in CMI- patients compared to CMI+ patients.

CONCLUSIONS

4D flow MRI with large volumetric coverage demonstrated significant differences in the redistribution of blood flow in SMA, SMV, and PV in CMI+ patients after a meal challenge. This approach may assist in the challenging diagnosis of CMI.

Deep learning-based velocity antialiasing of 4D-flow MRI

Purpose

To develop a convolutional neural network (CNN) for the robust and fast correction of velocity aliasing in 4D‐flow MRI.

Methods

This study included 667 adult subjects with aortic 4D‐flow MRI data with existing velocity aliasing (n = 362) and no velocity aliasing (n = 305). Additionally, 10 controls received back‐to‐back 4D‐flow scans with systemically varied velocity‐encoding sensitivity (vencs) at 60, 100, and 175 cm/s. The no‐aliasing data sets were used to simulate velocity aliasing by reducing the venc to 40%–70% of the original, alongside a ground truth locating all aliased voxels (153 training, 152 testing). The 152 simulated and 362 existing aliasing data sets were used for testing and compared with a conventional velocity antialiasing algorithm. Dice scores were calculated to quantify CNN performance. For controls, the venc 175‐cm/s scans were used as the ground truth and compared with the CNN‐corrected venc 60 and 100 cm/s data sets

Results

The CNN required 176 ± 30 s to perform compared with 162 ± 14 s for the conventional algorithm. The CNN showed excellent performance for the simulated data compared with the conventional algorithm (median range of Dice scores CNN: [0.89–0.99], conventional algorithm: [0.84–0.94], p < 0.001, across all simulated vencs) and detected more aliased voxels in existing velocity aliasing data sets (median detected CNN: 159 voxels [31–605], conventional algorithm: 65 [7–417], p < 0.001). For controls, the CNN showed Dice scores of 0.98 [0.95–0.99] and 0.96 [0.87–0.99] for venc = 60 cm/s and 100 cm/s, respectively, while flow comparisons showed moderate‐excellent agreement.

Conclusion

Deep learning enabled fast and robust velocity anti‐aliasing in 4D‐flow MRI.

MAXIMIZING UNAMBIGUOUS VELOCITY RANGE IN PHASE-CONTRAST MRI WITH MULTIPOINT ENCODING

In phase-contrast magnetic resonance imaging (PC-MRI), the velocity of spins at a voxel is encoded in the image phase. The strength of the velocity encoding gradient offers a trade-off between the velocity-to-noise ratio (VNR) and the extent of phase aliasing. Phase differences provide invariance to an unknown background phase. Existing literature proposes processing a reduced set of phase difference equations, simplifying the phase unwrapping problem at the expense of VNR or unaliased range of velocities, or both. Here, we demonstrate that the fullest unambiguous range of velocities is a parallelepiped, which can be accessed by jointly processing all phase differences. The joint processing also maximizes the velocity-to-noise ratio. The simple understanding of the unambiguous parallelepiped provides the potential for analyzing new multi-point acquisitions for an enhanced range of unaliased velocities; two examples are given.

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.

Virtual injections using 4D flow MRI with displacement corrections and constrained probabilistic streamlines

PURPOSE

Streamlines from 4D-flow MRI have been used clinically for intracranial blood-flow tracking. However, deterministic and stochastic errors degrade streamline quality. The purpose of this study is to integrate displacement corrections, probabilistic streamlines, and novel fluid constraints to improve selective blood-flow tracking and emulate "virtual bolus injections."

METHODS

Both displacement artifacts (deterministic) and velocity noise (stochastic) inherently occur during phase-contrast MRI acquisitions. Here, two displacement correction methods, single-step and iterative, were tested in silico with simulated displacements and were compared with ground-truth velocity fields. Next, the effects of combining displacement corrections and constrained probabilistic streamlines were performed in 10 healthy volunteers using time-averaged 4D-flow data. Measures of streamline length and depth into vasculature were then compared with streamlines generated with no corrections and displacement correction alone using one-way repeated-measures analysis of variance and Friedman's tests. Finally, virtual injections with improved streamlines were generated for three intracranial pathology cases.

RESULTS

Iterative displacement correction outperformed the single-step method in silico. In volunteers, the combination of displacement corrections and constrained probabilistic streamlines allowed for significant improvements in streamline length and increased the number of streamlines entering the circle of Willis relative to streamlines with no corrections and displacement correction alone. In the pathology cases, virtual injections with improved streamlines were qualitatively similar to dynamic arterial spin labeling images and allowed for forward/reverse selective flow tracking to characterize cerebrovascular malformations.

CONCLUSION

Virtual injections with improved streamlines from 4D-flow MRI allow for flexible, robust, intracranial flow tracking.

Identification of hemodynamic biomarkers for bicuspid aortic valve induced aortic dilation using machine learning

Recent advances in medical imaging have confirmed the presence of altered hemodynamics in bicuspid aortic valve (BAV) patients. Therefore, there is a need for new hemodynamic biomarkers to refine disease monitoring and improve patient risk stratification. This research aims to analyze and extract multiple correlation patterns of hemodynamic parameters from 4D Flow MRI data and find which parameters allow an accurate classification between healthy volunteers (HV) and BAV patients with dilated and non-dilated ascending aorta using machine learning. Sixteen hemodynamic parameters were calculated in the ascending aorta (AAo) and aortic arch (AArch) at peak systole from 4D Flow MRI. We used sequential forward selection (SFS) and principal component analysis (PCA) as feature selection algorithms. Then, eleven machine-learning classifiers were implemented to separate HV and BAV patients (non- and dilated ascending aorta). Multiple correlation patterns from hemodynamic parameters were extracted using hierarchical clustering. The linear discriminant analysis and random forest are the best performing classifiers, using five hemodynamic parameters selected with SFS (velocity angle, forward velocity, vorticity, and backward velocity in AAo; and helicity density in AArch) a 96.31 ± 1.76% and 96.00 ± 0.83% accuracy, respectively. Hierarchical clustering revealed three groups of correlated features. According to this analysis, we observed that features selected by SFS have a better performance than those selected by PCA because the five selected parameters were distributed according to 3 different clusters. Based on the proposed method, we concluded that the feature selection method found five potentially hemodynamic biomarkers related to this disease.

Simultaneous 3D-TOF angiography and 4D-flow MRI with enhanced flow signal using multiple overlapping thin slab acquisition and magnetization transfer

PURPOSE

To investigate the fusion of 3D time-of-flight principles into 4D-flow MRI to enhance vessel contrast and signal without an exogenous contrast agent, enabling simultaneous in-flow based angiograms.

METHODS

A 4D-flow MRI technique was developed consisting of multiple overlapping slabs with intermittent magnetization transfer preparation. The scan time penalty associated with multiple slab acquisitions was mitigated by using undersampled distributed spiral trajectories and compressed sensing reconstruction. A flow phantom was used to characterize in-flow enhancement, velocity noise improvement, and flow rate measurements against the single-slab 4D-flow MRI. In a patient-volunteer cohort (n = 15), magnitude-based angiograms were radiologically evaluated against 3D time-of-flight, and velocity measurements were compared pixel-wise against single-slab and contrast-enhanced 4D-flow MRI.

RESULTS

Multiple-slab acquisitions, together with magnetization transfer preparation, substantially improved vessel signal, contrast, and vessel conspicuity in magnitude angiograms. Both clinical 3D time-of-flight and the proposed technique produced equivalent vessel depictions with no statistically significant difference (p < .1). Both techniques also produced clear depictions of brain aneurysms in all patients; however, very small vessels tended to show reduced conspicuity in the proposed technique. Velocity measurements agreed with contrast-enhanced and single-slab scans with high correlations (R² = 0.941-0.974) and agreements (slopes = 0.994-1.071). Slab boundary and magnetization transfer-related artifacts were not observed in velocity measurements, and velocity noise was reduced with in-flow enhancement over single-slab scans (phantom).

CONCLUSION

The vessel signal and contrast can be improved in 4D-flow MRI without exogenous contrast agents by utilizing in-flow enhancement, efficient sampling, and compressed sensing. The in-flow enhancement also enables simultaneous 3D time-of-flight angiograms useful for flow quantification and diagnosis.

Cerebral arterial pulsatility is linked to hippocampal microvascular function and episodic memory in healthy older adults

Microvascular damage in the hippocampus is emerging as a central cause of cognitive decline and dementia in aging. This could be a consequence of age-related decreases in vascular elasticity, exposing hippocampal capillaries to excessive cardiac-related pulsatile flow that disrupts the blood-brain barrier and the neurovascular unit. Previous studies have found altered intracranial hemodynamics in cognitive impairment and dementia, as well as negative associations between pulsatility and hippocampal volume. However, evidence linking features of the cerebral arterial flow waveform to hippocampal function is lacking. We used a high-resolution 4D flow MRI approach to estimate global representations of the time-resolved flow waveform in distal cortical arteries and in proximal arteries feeding the brain in healthy older adults. Waveform-based clustering revealed a group of individuals featuring steep systolic onset and high amplitude that had poorer hippocampus-sensitive episodic memory (p = 0.003), lower whole-brain perfusion (p = 0.001), and weaker microvascular low-frequency oscillations in the hippocampus (p = 0.035) and parahippocampal gyrus (p = 0.005), potentially indicating compromised neurovascular unit integrity. Our findings suggest that aberrant hemodynamic forces contribute to cerebral microvascular and hippocampal dysfunction in aging.

Sildenafil administration improves right ventricular function on 4D flow MRI in young adults born premature

Extreme preterm birth conveys an elevated risk of heart failure by young adulthood. Smaller biventricular chamber size, diastolic dysfunction, and pulmonary hypertension may contribute to reduced ventricular-vascular coupling. However, how hemodynamic manipulations may affect right ventricular (RV) function and coupling remains unknown. As a pilot study, 4D flow MRI was used to assess the effect of afterload reduction and heart rate reduction on cardiac hemodynamics and function. Young adults born premature were administered sildenafil (a pulmonary vasodilator) and metoprolol (a β blocker) on separate days, and MRI with 4D flow completed before and after each drug administration. Endpoints include cardiac index (CI), direct flow fractions, and ventricular kinetic energy including E/A wave kinetic energy ratio. Sildenafil resulted in a median CI increase of 0.24 L/min/m2 (P = 0.02), mediated through both an increase in heart rate (HR) and stroke volume. Although RV ejection fraction improved only modestly, there was a significant increase (4% of end diastolic volume) in RV direct flow fraction (P = 0.04), consistent with hemodynamic improvement. Metoprolol administration resulted in a 5-beats/min median decrease in HR (P = 0.01), a 0.37 L/min/m2 median decrease in CI (P = 0.04), and a reduction in time-averaged kinetic energy (KE) in both ventricles (P < 0.01), despite increased RV diastolic E/A KE ratio (P = 0.04). Despite reduced right atrial workload, metoprolol significantly depressed overall cardiac systolic function. Sildenafil, however, increased CI and improved RV function, as quantified by the direct flow fraction. The preterm heart appears dependent on HR but sensitive to RV afterload manipulations.NEW & NOTEWORTHY We assessed the effect of right ventricular afterload reduction with sildenafil and heart rate reduction with metoprolol on cardiac hemodynamics and function in young adults born premature using 4D flow MRI. Metoprolol depressed cardiac function, whereas sildenafil improved cardiac function including right ventricular direct flow fraction by 4D flow, consistent with hemodynamic improvement. This suggests that the preterm heart is dependent on heart rate and sensitive to right ventricular afterload changes.

Assessment of vascular stiffness in the internal carotid artery proximal to the carotid canal in Alzheimer's disease using pulse wave velocity from low rank reconstructed 4D flow MRI

Clinical evidence shows vascular factors may co-occur and complicate the expression of Alzheimer's disease (AD); yet, the pathologic mechanisms and involvement of different compartments of the vascular network are not well understood. Diseases such as arteriosclerosis diminish vascular compliance and will lead to arterial stiffness, a well-established risk factor for cardiovascular morbidity. Arterial stiffness can be assessed using pulse wave velocity (PWV); however, this is usually done from carotid-to-femoral artery ratios. To probe the brain vasculature, intracranial PWV measures would be ideal. In this study, high temporal resolution 4D flow MRI was used to assess transcranial PWV in 160 subjects including AD, mild cognitive impairment (MCI), healthy controls, and healthy subjects with apolipoprotein ɛ4 positivity (APOE4+) and parental history of AD dementia (FH+). High temporal resolution imaging was achieved by high temporal binning of retrospectively gated data using a local-low rank approach. Significantly higher transcranial PWV in AD dementia and MCI subjects was found when compared to old-age-matched controls (AD vs. old-age-matched controls: P <0.001, AD vs. MCI: P = 0.029, MCI vs. old-age-matched controls P = 0.013). Furthermore, vascular changes were found in clinically healthy middle-age adults with APOE4+ and FH+ indicating significantly higher transcranial PWV compared to controls (P <0.001).

Middle cerebral artery pressure laterality in patients with symptomatic ICA stenosis

An internal carotid artery (ICA) stenosis can potentially decrease the perfusion pressure to the brain. In this study, computational fluid dynamics (CFD) was used to study if there was a hemispheric pressure laterality between the contra- and ipsilateral middle cerebral artery (MCA) in patients with a symptomatic ICA stenosis. We further investigated if this MCA pressure laterality (ΔPMCA) was related to the hemispheric flow laterality (ΔQ) in the anterior circulation, i.e., ICA, proximal MCA and the proximal anterior cerebral artery (ACA). Twenty-eight patients (73±6 years, range 59-80 years, 21 men) with symptomatic ICA stenosis were included. Flow rates were measured using 4D flow MRI data (PC-VIPR) and vessel geometries were obtained from computed tomography angiography. The ΔPMCA was calculated from CFD, where patient-specific flow rates were applied at all input- and output boundaries. The ΔPMCA between the contra- and ipsilateral side was 6.4±8.3 mmHg (p<0.001) (median 3.9 mmHg, range -1.3 to 31.9 mmHg). There was a linear correlation between the ΔPMCA and ΔQICA (r = 0.85, p<0.001) and ΔQACA (r = 0.71, p<0.001), respectively. The correlation to ΔQMCA was weaker (r = 0.47, p = 0.011). In conclusion, the MCA pressure laterality obtained with CFD, is a promising physiological biomarker that can grade the hemodynamic disturbance in patients with a symptomatic ICA stenosis.

Modeling and Reducing the Effect of Geometric Uncertainties in Intracranial Aneurysms with Polynomial Chaos Expansion, Data Decomposition, and 4D-Flow MRI

PURPOSE

Variations in the vessel radius of segmented surfaces of intracranial aneurysms significantly influence the fluid velocities given by computer simulations. It is important to generate models that capture the effect of these variations in order to have a better interpretation of the numerically predicted hemodynamics. Also, it is highly relevant to develop methods that combine experimental observations with uncertainty modeling to get a closer approximation to the blood flow behavior.

METHODS

This work applies polynomial chaos expansion to model the effect of geometric uncertainties on the simulated fluid velocities of intracranial aneurysms. The radius of the vessel is defined as the uncertainty variable. Proper orthogonal decomposition is applied to characterize the solution space of fluid velocities. Next, a process of projecting the 4D-Flow MRI velocities on the basis vectors followed by coefficient mapping using generalized dynamic mode decomposition enables the merging of 4D-Flow MRI with the uncertainty propagated fluid velocities.

RESULTS

Polynomial chaos expansion propagates the fluid velocities with an error of 2% in velocity magnitude relative to computer simulations. Also, the bifurcation region (or impingement location) shows a standard deviation of 0.17 m/s (since an available reported variance in the vessel radius is adopted to model the uncertainty, the expected standard deviation may be different). Numerical phantom experiments indicate that the proposed approach reconstructs the fluid velocities with 0.3% relative error in presence of geometric uncertainties.

CONCLUSION

Polynomial chaos expansion is an effective approach to propagate the effect of the uncertainty variable in the blood flow velocities of intracranial aneurysms. Merging 4D-Flow MRI and uncertainty propagated fluid velocities leads to more realistic flow trends relative to ignoring the uncertainty in the vessel radius.

A technique for intra-procedural blood velocity quantitation using time-resolved 2D digital subtraction angiography

BACKGROUND

2D digital subtraction angiography (DSA) is utilized qualitatively to assess blood velocity changes that occur during arterial interventions. Quantitative angiographic metrics, such as blood velocity, could be used to standardize endpoints during angiographic interventions.

PURPOSE

To assess the accuracy and precision of a quantitative 2D DSA (qDSA) technique and to determine its feasibility for in vivo measurements of blood velocity.

MATERIALS AND METHODS

A quantitative DSA technique was developed to calculate intra-procedural blood velocity. In vitro validation was performed by comparing velocities from the qDSA method and an ultrasonic flow probe in a bifurcation phantom. Parameters of interest included baseline flow rate, contrast injection rate, projection angle, and magnification. In vivo qDSA analysis was completed in five different branches of the abdominal aorta in two 50 kg swine and compared to 4D Flow MRI. Linear regression, Bland-Altman, Pearson's correlation coefficient and chi squared tests were used to assess the accuracy and precision of the technique.

RESULTS

In vitro validation showed strong correlation between qDSA and flow probe velocities over a range of contrast injection and baseline flow rates (slope = 1.012, 95% CI [0.989,1.035], Pearson's r = 0.996, p < .0001). The application of projection angle and magnification corrections decreased variance to less than 5% the average baseline velocity (p = 0.999 and p = 0.956, respectively). In vivo validation showed strong correlation with a small bias between qDSA and 4D Flow MRI velocities for all five abdominopelvic arterial vessels of interest (slope = 1.01, Pearson's r = 0.880, p = <.01, Bias = 0.117 cm/s).

CONCLUSION

The proposed method allows for accurate and precise calculation of blood velocities, in near real-time, from time resolved 2D DSAs.

Techniques for Imaging Vascular Supply of Peripheral Nerves

Few studies have been developed to map the vascular structures feeding peripheral nerves, with the majority using cadaveric models and inadequate sample sizes. Preliminary evidence, while limited, indicates that the mapping of these vessels may allow or preclude certain procedures in nerve reconstruction due to the location of essential arterial inflow to the vasa nervorum. This review evaluates the evidence regarding historical, current, and emerging techniques for visualizing these vascular structures in vivo and considers their potential application in peripheral nerve vasculature.

Additional value and new insights by four-dimensional flow magnetic resonance imaging in congenital heart disease: application in neonates and young children

Cardiovascular MRI has become an essential imaging modality in children with congenital heart disease (CHD) in the last 15-20 years. With use of appropriate sequences, it provides important information on cardiovascular anatomy, blood flow and function for initial diagnosis and post-surgical or -interventional monitoring in children. Although considered as more sophisticated and challenging than CT, in particular in neonates and infants, MRI is able to provide information on intra- and extracardiac haemodynamics, in contrast to CT. In recent years, four-dimensional (4-D) flow MRI has emerged as an additional MR technique for retrospective assessment and visualisation of blood flow within the heart and any vessel of interest within the acquired three-dimensional (3-D) volume. Its application in young children requires special adaptations for the smaller vessel size and faster heart rate compared to adolescents or adults. In this article, we provide an overview of 4-D flow MRI in various types of complex CHD in neonates and infants to demonstrate its potential indications and beneficial application for optimised individual cardiovascular assessment. We focus on its application in clinical routine cardiovascular workup and, in addition, show some examples with pathologies other than CHD to highlight that 4-D flow MRI yields new insights in disease understanding and therapy planning. We shortly review the essentials of 4-D flow data acquisition, pre- and post-processing techniques in neonates, infants and young children. Finally, we conclude with some details on accuracy, limitations and pitfalls of the technique.

Super-resolution and denoising of 4D-Flow MRI using physics-Informed deep neural nets

BACKGROUND AND OBJECTIVE

Time resolved three-dimensional phase contrast magnetic resonance imaging (4D-Flow MRI) has been used to non-invasively measure blood velocities in the human vascular system. However, issues such as low spatio-temporal resolution, acquisition noise, velocity aliasing, and phase-offset artifacts have hampered its clinical application. In this research, we developed a purely data-driven method for super-resolution and denoising of 4D-Flow MRI.

METHODS

The flow velocities, pressure, and the MRI image magnitude are modeled as a patient-specific deep neural net (DNN). For training, 4D-Flow MRI images in the complex Cartesian space are used to impose data-fidelity. Physics of fluid flow is imposed through regularization. Creative loss function terms have been introduced to handle noise and super-resolution. The trained patient-specific DNN can be sampled to generate noise-free high-resolution flow images. The proposed method has been implemented using the TensorFlow DNN library and tested on numerical phantoms and validated in-vitro using high-resolution particle image velocitmetry (PIV) and 4D-Flow MRI experiments on transparent models subjected to pulsatile flow conditions.

RESULTS

In case of numerical phantoms, we were able to increase spatial resolution by a factor of 100 and temporal resolution by a factor of 5 compared to the simulated 4D-Flow MRI. There is an order of magnitude reduction of velocity normalized root mean square error (vNRMSE). In case of the in-vitro validation tests with PIV as reference, there is similar improvement in spatio-temporal resolution. Although the vNRMSE is reduced by 50%, the method is unable to negate a systematic bias with respect to the reference PIV that is introduced by the 4D-Flow MRI measurement.

CONCLUSIONS

This work has demonstrated the feasibility of using the readily available machinery of deep learning to enhance 4D-Flow MRI using a purely data-driven method. Unlike current state-of-the-art methods, the proposed method is agnostic to geometry and boundary conditions and therefore eliminates the need for tedious tasks such as accurate image segmentation for geometry, image registration, and estimation of boundary flow conditions. Arbitrary regions of interest can be selected for processing. This work will lead to user-friendly analysis tools that will enable quantitative hemodynamic analysis of vascular diseases in a clinical setting.

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.

In Vitro Assessment of Flow Variability in an Intracranial Aneurysm Model Using 4D Flow MRI and Tomographic PIV

Aneurysm rupture has been suggested to be related to aneurysm geometry, morphology, and complex flow activity; therefore, understanding aneurysm-specific hemodynamics is crucial. 4D Flow MRI has been shown to be a feasible tool for assessing hemodynamics in intracranial aneurysms with high spatial resolution. However, it requires averaging over multiple heartbeats and cannot account for cycle-to-cycle hemodynamics variations. This study aimed to assess cycle-to-cycle flow dynamics variations in a patient-specific intracranial aneurysm model using tomographic particle image velocimetry (tomo-PIV) at a high image rate under pulsatile flow conditions. Time-resolved and time-averaged velocity flow fields within the aneurysm sac and estimations of wall shear stress (WSS) were compared with those from 4D Flow MRI. A one-way ANOVA showed a significant difference between cardiac cycles (p value < 0.0001); however, differences were not significant after PIV temporal and spatial resolution was matched to that of MRI (p value 0.9727). This comparison showed the spatial resolution to be the main contributor to assess cycle-to-cycle variability. Furthermore, the comparison with 4D Flow MRI between velocity components, streamlines, and estimated WSS showed good qualitative and quantitative agreement. This study showed the feasibility of patient-specific in-vitro experiments using tomo-PIV to assess 4D Flow MRI with high repeatability in the measurements.

Towards multi-modal data fusion for super-resolution and denoising of 4D-Flow MRI

4D-Flow magnetic resonance imaging (MRI) has enabled in vivo time-resolved measurement of three-dimensional blood flow velocities in the human vascular system. However, its clinical use has been hampered by two main issues, namely, low spatio-temporal resolution and acquisition noise. While patient-specific computational fluid dynamics (CFD) simulations can address the resolution and noise issues, its fidelity is impacted by accuracy of estimation of boundary conditions, model parameters, vascular geometry, and flow model assumptions. In this paper a scheme to address limitations of both modalities through data-fusion is presented. The solutions of the patient-specific CFD simulation are characterized using proper orthogonal decomposition (POD). Next, a process of projecting the 4D-Flow MRI data onto the POD basis and projection coefficient mapping using generalized dynamic mode decomposition (DMD) enables simultaneous super-resolution and denoising of 4D-Flow MRI. The method has been tested using numerical phantoms derived from patient-specific aneurysmal geometries and applied to in vivo 4D-Flow MRI data.

Highly accelerated 4D flow cardiovascular magnetic resonance using a pseudo-spiral Cartesian acquisition and compressed sensing reconstruction for carotid flow and wall shear stress

BACKGROUND

4D flow cardiovascular magnetic resonance (CMR) enables visualization of complex blood flow and quantification of biomarkers for vessel wall disease, such as wall shear stress (WSS). Because of the inherently long acquisition times, many efforts have been made to accelerate 4D flow acquisitions, however, no detailed analysis has been made on the effect of Cartesian compressed sensing accelerated 4D flow CMR at different undersampling rates on quantitative flow parameters and WSS.

METHODS

We implemented a retrospectively triggered 4D flow CMR acquisition with pseudo-spiral Cartesian k-space filling, which results in incoherent undersampling of k-t space. Additionally, this strategy leads to small jumps in k-space thereby minimizing eddy current related artifacts. The pseudo-spirals were rotated in a tiny golden-angle fashion, which provides optimal incoherence and a variable density sampling pattern with a fully sampled center. We evaluated this 4D flow protocol in a carotid flow phantom with accelerations of R = 2-20, as well as in carotids of 7 healthy subjects (27 ± 2 years, 4 male) for R = 10-30. Fully sampled 2D flow CMR served as a flow reference. Arteries were manually segmented and registered to enable voxel-wise comparisons of both velocity and WSS using a Bland-Altman analysis.

RESULTS

Magnitude images, velocity images, and pathline reconstructions from phantom and in vivo scans were similar for all accelerations. For the phantom data, mean differences at peak systole for the entire vessel volume in comparison to R = 2 ranged from - 2.3 to - 5.3% (WSS) and - 2.4 to - 2.2% (velocity) for acceleration factors R = 4-20. For the in vivo data, mean differences for the entire vessel volume at peak systole in comparison to R = 10 were - 9.9, - 13.4, and - 16.9% (WSS) and - 8.4, - 10.8, and - 14.0% (velocity), for R = 20, 25, and 30, respectively. Compared to single slice 2D flow CMR acquisitions, peak systolic flow rates of the phantom showed no differences, whereas peak systolic flow rates in the carotid artery in vivo became increasingly underestimated with increasing acceleration.

CONCLUSION

Acquisition of 4D flow CMR of the carotid arteries can be highly accelerated by pseudo-spiral k-space sampling and compressed sensing reconstruction, with consistent data quality facilitating velocity pathline reconstructions, as well as quantitative flow rate and WSS estimations. At an acceleration factor of R = 20 the underestimation of peak velocity and peak WSS was acceptable (< 10%) in comparison to an R = 10 accelerated 4D flow CMR reference scan. Peak flow rates were underestimated in comparison with 2D flow CMR and decreased systematically with higher acceleration factors.

Low and Oscillatory Wall Shear Stress Is Not Related to Aortic Dilation in Patients With Bicuspid Aortic Valve: A Time-Resolved 3-Dimensional Phase-Contrast Magnetic Resonance Imaging Study

Supplemental Digital Content is available in the text.

Objective:

To assess the relationship between regional wall shear stress (WSS) and oscillatory shear index (OSI) and aortic dilation in patients with bicuspid aortic valve (BAV).

Approach and Results:

Forty-six consecutive patients with BAV (63% with right-left-coronary-cusp fusion, aortic diameter ≤ 45 mm and no severe valvular disease) and 44 healthy volunteers were studied by time-resolved 3-dimensional phase-contrast magnetic resonance imaging. WSS and OSI were quantified at different levels of the ascending aorta and the aortic arch, and regional WSS and OSI maps were obtained. Seventy percent of BAV had ascending aorta dilation. Compared with healthy volunteers, patients with BAV had increased WSS and decreased OSI in most of the ascending aorta and the aortic arch. In both BAV and healthy volunteers, regions of high WSS matched regions of low OSI and vice versa. No regions of both low WSS and high OSI were identified in BAV compared with healthy volunteers. Patients with BAV with dilated compared with nondilated aorta presented low and oscillatory WSS in the aortic arch, but not in the ascending aorta where dilation is more prevalent. Furthermore, no regions of concomitant low WSS and high OSI were identified when BAV were compared according to leaflet fusion pattern, despite the well-known differences in regional dilation prevalence.

Conclusions:

Regions with low WSS and high OSI do not match those with the highest prevalence of dilation in patients with BAV, thus providing no evidence to support the low and oscillatory shear stress theory in the pathogenesis of proximal aorta dilation in the presence of BAV.

Assessment of Cerebral Blood Flow Pulsatility and Cerebral Arterial Compliance With 4D Flow MRI

BACKGROUND

Four-dimensional flow magnetic resonance imaging (4D flow MRI) enables efficient investigation of cerebral blood flow pulsatility in the cerebral arteries. This is important for exploring hemodynamic mechanisms behind vascular diseases associated with arterial pulsations.

PURPOSE

To investigate the feasibility of pulsatility assessments with 4D flow MRI, its agreement with reference two-dimensional phase-contrast MRI (2D PC-MRI) measurements, and to demonstrate how 4D flow MRI can be used to assess cerebral arterial compliance and cerebrovascular resistance in major cerebral arteries.

STUDY TYPE

Prospective.

SUBJECTS

Thirty-five subjects (20 women, 79 ± 5 years, range 70-91 years).

FIELD STRENGTH/SEQUENCE

4D flow MRI (PC-VIPR) and 2D PC-MRI acquired with a 3T scanner.

ASSESSMENT

Time-resolved flow was assessed in nine cerebral arteries. From the pulsatile flow waveform in each artery, amplitude (ΔQ), volume load (ΔV), and pulsatility index (PI) were calculated. To reduce high-frequency noise in the 4D flow MRI data, the flow waveforms were low-pass filtered. From the total cerebral blood flow, total PI (PI_tot ), total volume load (ΔV_tot ), cerebral arterial compliance (C), and cerebrovascular resistance (R) were calculated.

STATISTICAL TESTS

Two-tailed paired t-test, intraclass correlation (ICC).

RESULTS

There was no difference in ΔQ between 4D flow MRI and the reference (0.00 ± 0.022 ml/s, mean ± SEM, P = 0.97, ICC = 0.95, n = 310) with a cutoff frequency of 1.9 Hz and 15 cut plane long arterial segments. For ΔV, the difference was -0.006 ± 0.003 ml (mean ± SEM, P = 0.07, ICC = 0.93, n = 310) without filtering. Total R was 11.4 ± 2.41 mmHg/(ml/s) (mean ± SD) and C was 0.021 ± 0.009 ml/mmHg (mean ± SD). ΔV_tot was 1.21 ± 0.29 ml (mean ± SD) with an ICC of 0.82 compared with the reference. PI_tot was 1.08 ± 0.21 (mean ± SD).

DATA CONCLUSION

We successfully assessed 4D flow MRI cerebral arterial pulsatility, cerebral arterial compliance, and cerebrovascular resistance. Averaging of multiple cut planes and low-pass filtering was necessary to assess accurate peak-to-peak features in the flow rate waveforms.

LEVEL OF EVIDENCE

2 Technical Efficacy Stage: 2 J. Magn. Reson. Imaging 2020;51:1516-1525.

Assessment of spinal cord motion as a new diagnostic MRI-parameter in cervical spinal canal stenosis: study protocol on a prospective longitudinal trial

Background

Increased spinal cord motion has been proven to be a relevant finding within spinal canal stenosis disclosed by phase-contrast MRI (PC-MRI). Adapted PC-MRI is a suitable and reliable method within the well deliberated setting. As the decision between conservative and operative treatment can be challenging in some cases, further diagnostic marker would facilitate the diagnostic process. We hypothesize that increased spinal cord motion will correlate to clinical course and functional impairment and will contribute as a new diagnostic marker.

Methods

A monocentric, prospective longitudinal observational trial on cervical spinal canal stenosis will be conducted at the University Medical Center Freiburg. Patients (n = 130) with relevant cervical spinal canal stenosis, being defined by at least contact to the spinal cord, will be included. Also, we will examine a control group of healthy volunteers (n = 20) as proof-of-principle. We will observe two openly assigned branches of participants undergoing conservative and surgical decompressive treatment (based on current German Guidelines) over a time course of 12 month, including a total of 4 visits. We will conduct a broad assessment of clinical parameters, standard scores and gradings, electrophysiological measurements, standard MRI, and adapted functional PC-MRI of spinal cord motion. Primary endpoint is the evaluation of an expected negative correlation of absolute spinal cord displacement to clinical impairment. Secondary endpoints are the evaluation of positive correlation of increased absolute spinal cord displacement to prolonged evoked potentials, prediction of clinical course by absolute spinal cord displacement, and demonstration of normalized spinal cord motion after decompressive surgery.

Discussion

With the use of adapted, non-invasive PC-MRI as a quantitative method for assessment of spinal cord motion, further objective diagnostic information can be gained, that might improve the therapeutic decision-making process. This study will offer the needed data in order to establish PC-MRI on spinal cord motion within the diagnostic work-up of patients suffering from spinal canal stenosis.

Trial registration

German Clinical Trials Register, ID: DRKS00012962, Register date 2018/01/17

Efficient triple-VENC phase-contrast MRI for improved velocity dynamic range

PURPOSE

To evaluate the utility of an efficient triple velocity-encoding (VENC) 4D flow MRI implementation to improve velocity unwrapping of 4D flow MRI data with the same scan time as an interleaved dual-VENC acquisition.

METHODS

A balanced 7-point acquisition was used to derive 3 sets of 4D flow images corresponding to 3 different VENCs. These 3 datasets were then used to unwrap the aliased lowest VENC into a minimally aliased, triple-VENC dataset. Triple-VENC MRI was evaluated and compared with dual-VENC MRI over 3 different VENC ranges (50-150, 60-150, and 60-180 cm/s) in vitro in a steadily rotating phantom as well as in a pulsatile flow phantom. In vivo, triple-VENC data of the thoracic aorta were also evaluated in 3 healthy volunteers (2 males, 26-44 years old) with VENC = 50/75/150 cm/s. Two triple-VENC (triconditional and biconditional) and 1 dual-VENC unwrapping algorithms were quantitatively assessed through comparison to a reference, unaliased, single-VENC scan.

RESULTS

Triple-VENC 4D flow constant rotation phantom results showed high correlation with the analytical solution (intraclass correlation coefficient = 0.984-0.995, P < .001) and up to a 61% reduction in velocity noise compared with the corresponding single-VENC scans (VENC = 150, 180 cm/s). Pulsatile flow phantom experiments demonstrated good agreement between triple-VENC and single-VENC acquisitions (peak flow < 0.8% difference; peak velocity < 11.7% difference). Triconditional triple-VENC unwrapping consistently outperformed dual-VENC unwrapping, correctly unwrapping more than 83% and 46%-66% more voxels in vitro and in vivo, respectively.

CONCLUSION

Triple-VENC 4D flow MRI adds no additional scan time to dual-VENC MRI and has the potential for improved unwrapping to extend the velocity dynamic range beyond dual-VENC methods.

Increased rotational flow in the proximal aortic arch is associated with its dilation in bicuspid aortic valve disease

Aims

Aortic dilation in bicuspid aortic valve (BAV) might extend to the proximal arch. Arch flow dynamics and their relationship with this segment dilation are still unexplored. Using 4D-flow cardiovascular magnetic resonance, we analysed flow dynamics in the arch for each BAV morphotype and their association with this segment dilation.

Methods and results

One hundred and eleven BAV patients (aortic diameters ≤55 mm, non-severe valvular disease), 21 age-matched tricuspid aortic valve (TAV) patients with dilated arch and 24 healthy volunteers (HV) underwent 4D-flow. BAV were classified per fusion morphotype: 75% right-left (RL-BAV), and per arch dilation: 57% dilated, mainly affecting the right-noncoronary (RN) BAV (86% dilated vs. 47% in RL-BAV). Peak velocity, jet angle, normalized displacement, in-plane rotational flow (IRF), wall shear stress, and systolic flow reversal ratio (SFRR) were calculated along the thoracic aorta. ANCOVA and multivariate linear regression analyses were used to identify correlates of arch dilation. BAV had higher rotational flow and eccentricity than TAV in the proximal arch. Dilated compared with non-dilated BAV had higher IRF being more pronounced in the RN-morphotype. RN-BAV, IRF, and SFRR were independently associated with arch dilation. Aortic stenosis and male sex were independently associated with arch dilation in RL-BAV. Flow parameters associated with dilation converged to the values found in HV in the distal arch.

Conclusion

Increased rotational flow could explain dilation of the proximal arch in RN-BAV and in RL-BAV patients of male sex and with valvular stenosis. These patients may benefit from a closer follow-up with cardiac magnetic resonance or computed tomography.

Accuracy of blood flow assessment in cerebral arteries with 4D flow MRI: Evaluation with three segmentation methods

Background

Accelerated 4D flow MRI allows for high‐resolution velocity measurements with whole‐brain coverage. Such scans are increasingly used to calculate flow rates of individual arteries in the vascular tree, but detailed information about the accuracy and precision in relation to different postprocessing options is lacking.

Purpose

To evaluate and optimize three proposed segmentation methods and determine the accuracy of in vivo 4D flow MRI blood flow rate assessments in major cerebral arteries, with high‐resolution 2D PCMRI as a reference.

Study Type

Prospective.

Subjects

Thirty‐five subjects (20 women, 79 ± 5 years, range 70–91 years).

Field Strength/Sequence

4D flow MRI with PC‐VIPR and 2D PCMRI acquired with a 3 T scanner.

Assessment

We compared blood flow rates measured with 4D flow MRI, to the reference, in nine main cerebral arteries. Lumen segmentation in the 4D flow MRI was performed with k‐means clustering using four different input datasets, and with two types of thresholding methods. The threshold was defined as a percentage of the maximum intensity value in the complex difference image. Local and global thresholding approaches were used, with evaluated thresholds from 6–26%.

Statistical Tests

Paired t‐test, F‐test, linear correlation (P < 0.05 was considered significant) along with intraclass correlation (ICC).

Results

With the thresholding methods, the lowest average flow difference was obtained for 20% local (0.02 ± 15.0 ml/min, ICC = 0.97, n = 310) or 10% global (0.08 ± 17.3 ml/min, ICC = 0.97, n = 310) thresholding with a significant lower standard deviation for local (F‐test, P = 0.01). For all clustering methods, we found a large systematic underestimation of flow compared with 2D PCMRI (16.1–22.3 ml/min).

Data Conclusion

A locally adapted threshold value gives a more stable result compared with a globally fixed threshold. 4D flow with the proposed segmentation method has the potential to become a useful reliable clinical tool for assessment of blood flow in the major cerebral arteries.

Level of Evidence: 2

Technical Efficacy: Stage 2

J. Magn. Reson. Imaging 2019;50:511–518.

Four-dimensional Flow MRI as a Marker for Risk Stratification of Gastroesophageal Varices in Patients with Liver Cirrhosis

Purpose To assess the feasibility of four-dimensional (4D) flow MRI as a noninvasive imaging marker for stratifying the risk of variceal bleeding in patients with liver cirrhosis. Materials and Methods This study recruited participants scheduled for both liver MRI and gastroesophageal endoscopy. Risk of variceal bleeding was assessed at endoscopy by using a three-point scale: no varices, low risk, and high risk requiring treatment. Four-dimensional flow MRI was used to create angiograms for evaluating visibility of varices and to measure flow volumes in main portal vein (PV), superior mesenteric vein, splenic vein (SV), and azygos vein. Fractional flow changes in PV and SV were calculated to quantify shunting (outflow) from PV and SV into varices. Logistic analysis was used to identify the independent indicator of high-risk varices. Results There were 23 participants (mean age, 52.3 years; age range, 25-75 years), including 14 men (mean age, 51.7 years; age range, 25-75 years) and nine women (mean age, 53.2 years; age range, 31-72 years) with no varices (n = 8), low-risk varices (n = 8), and high-risk varices (n = 7) determined at endoscopy. Four-dimensional flow MRI-based angiography helped radiologists to view varices in four of 15 participants with varices. Independent indicators of high-risk varices were flow volume in the azygos vein greater than 0.1 L/min (P = .034; 100% sensitivity [seven of seven] and 62% specificity [10 of 16]) and fractional flow change in PV of less than 0 (P < .001; 100% sensitivity [seven of seven] and 94% specificity [15 of 16]). Conclusion Azygos flow greater than 0.1 L/min and portal venous flow less than the sum of splenic and superior mesenteric vein flow are useful markers to stratify the risk of gastroesophageal varices bleeding in patients with liver cirrhosis. © RSNA, 2018 Online supplemental material is available for this article.