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

Measuring global cerebrovascular pulsatility transmission using 4D flow MRI

Pulse wave encephalopathy (PWE) is hypothesised to initiate many forms of dementia, motivating its identification and risk assessment. As candidate pulsatility based biomarkers for PWE, pulsatility index and pulsatility damping have been studied and, currently, do not adequately stratify risk due to variability in pulsatility and spatial bias. Here, we propose a locus-independent pulsatility transmission coefficient computed by spatially tracking pulsatility along vessels to characterise the brain pulse dynamics at a whole-organ level. Our preliminary analyses in a cohort of 20 subjects indicate that this measurement agrees with clinical observations relating blood pulsatility with age, heart rate, and sex, making it a suitable candidate to study the risk of PWE. We identified transmission differences between vascular regions perfused by the basilar and internal carotid arteries attributed to the identified dependence on cerebral blood flow, and some participants presented differences between the internal carotid perfused regions that were not related to flow or pulsatility burden, suggesting underlying mechanical differences. Large populational studies would benefit from retrospective pulsatility transmission analyses, providing a new comprehensive arterial description of the hemodynamic state in the brain. We provide a publicly available implementation of our tools to derive this coefficient, built into pre-existing open-source software.

Cerebral Blood Flow Patterns in Patients With Low-Flow Carotid Artery Stenosis, a 4D-PCMRI Assessment

BACKGROUND

Compromised cerebral blood flow can contribute to future ischemic events in patients with symptomatic carotid artery disease. However, there is limited knowledge of the effects on cerebral hemodynamics resulting from a reduced internal carotid artery (ICA) blood flow rate (BFR).

PURPOSE

Investigate how reduced ICA-BFR, relates to BFR in the cerebral arteries.

STUDY TYPE

Prospective.

SUBJECTS

Thirty-eight patients, age 72 ± 6 years (11 female).

FIELD STRENGTH/SEQUENCE

3-Tesla, four-dimensional phase-contrast magnetic resonance imaging (4D-PCMRI).

ASSESSMENT

Patients with ischemic stroke or transient ischemic attack were evaluated regarding the degree of stenosis. 4D-PCMRI was used to measure cerebral BFR in 38 patients with symptomatic carotid stenosis (≥50%). BFR in the cerebral arteries was assessed in two subgroups based on symptomatic ICA-BFR: reduced ICA-flow (<160 mL/minutes) and preserved ICA-flow (≥160 mL/minutes). BFR laterality was defined as a difference in the paired ipsilateral-contralateral arteries.

STATISTICAL TESTS

Patients were grouped based on ICA-BFR (reduced vs. preserved). Statistical tests (independent sample t-test/paired t-test) were used to compare groups and hemispheres. Significance was determined at P < 0.05.

RESULTS

The degree of stenosis was not significantly different, 80% (95% confidence interval [CI] = 73%-87%) in the reduced ICA-flow vs. 72% (CI = 66%-76%) in the preserved ICA-flow; P = 0.09. In the reduced ICA-flow group, a significantly reduced BFR was found in the ipsilateral middle cerebral artery and anterior cerebral artery (A1), while significantly increased in the contralateral A1. Retrograde BFR was found in the posterior communicating artery and ophthalmic artery. Significant BFR laterality was present in all paired arteries in the reduced ICA-flow group, contrasting the preserved ICA-flow group (P = 0.14-0.93).

DATA CONCLUSIONS

4D-PCMRI revealed compromised cerebral BFR due to carotid stenosis, not possible to detect by solely analyzing the degree of stenosis. In patients with reduced ICA-flow, collaterals were not sufficient to maintain symmetrical BFR distribution to the two hemispheres.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 3.

Automatic 4D Flow MRI Segmentation Using the Standardized Difference of Means Velocity

We present a method to automatically segment 4D flow magnetic resonance imaging (MRI) by identifying net flow effects using the standardized difference of means (SDM) velocity. The SDM velocity quantifies the ratio between the net flow and observed flow pulsatility in each voxel. Vessel segmentation is performed using an F-test, identifying voxels with significantly higher SDM velocity values than background voxels. We compare the SDM segmentation algorithm against pseudo-complex difference (PCD) intensity segmentation of 4D flow measurements in in vitro cerebral aneurysm models and 10 in vitro Circle of Willis (CoW) datasets. We also compared the SDM algorithm to convolutional neural network (CNN) segmentation in 5 thoracic vasculature datasets. The in vitro flow phantom geometry is known, while the ground truth geometries for the CoW and thoracic aortas are derived from high-resolution time-of-flight (TOF) magnetic resonance angiography and manual segmentation, respectively. The SDM algorithm demonstrates greater robustness than PCD and CNN approaches and can be applied to 4D flow data from other vascular territories. The SDM to PCD comparison demonstrated an approximate 48% increase in sensitivity in vitro and 70% increase in the CoW, respectively; the SDM and CNN sensitivities were similar. The vessel surface derived from the SDM method was 46% closer to the in vitro surfaces and 72% closer to the in vitro TOF surfaces than the PCD approach. The SDM and CNN approaches both accurately identify vessel surfaces. The SDM algorithm is a repeatable segmentation method, enabling reliable computation of hemodynamic metrics associated with cardiovascular disease.

Automated hemodynamic assessment for cranial 4D flow MRI

Cranial 4D flow MRI post-processing typically involves manual user interaction which is time-consuming and associated with poor repeatability. The primary goal of this study is to develop a robust quantitative velocity tool (QVT) that utilizes threshold-based segmentation techniques to improve segmentation quality over prior approaches based on centerline processing schemes (CPS) that utilize k-means clustering segmentation. This tool also includes an interactive 3D display designed for simplified vessel selection and automated hemodynamic visualization and quantification. The performances of QVT and CPS were compared in vitro in a flow phantom and in vivo in 10 healthy participants. Vessel segmentations were compared with ground-truth computed tomography in vitro (29 locations) and manual segmentation in vivo (13 locations) using linear regression. Additionally, QVT and CPS MRI flow rates were compared to perivascular ultrasound flow in vitro using linear regression. To assess internal consistency of flow measures in vivo, conservation of flow was assessed at vessel junctions using linear regression and consistency of flow along vessel segments was analyzed by fitting a Gaussian distribution to a histogram of normalized flow values. Post-processing times were compared between the QVT and CPS using paired t-tests. Vessel areas segmented in vitro (CPS: slope = 0.71, r = 0.95 and QVT: slope = 1.03, r = 0.95) and in vivo (CPS: slope = 0.61, r = 0.96 and QVT: slope = 0.93, r = 0.96) were strongly correlated with ground-truth area measurements. However, CPS (using k-means segmentation) consistently underestimated vessel areas. Strong correlations were observed between QVT and ultrasound flow (slope = 0.98, r = 0.96) as well as flow at junctions (slope = 1.05, r = 0.98). Mean and standard deviation of flow along vessel segments was 9.33e-16 ± 3.05%. Additionally, the QVT demonstrated excellent interobserver agreement and significantly reduced post-processing by nearly 10 min (p < 0.001). By completely automating post-processing and providing an easy-to-use 3D visualization interface for interactive vessel selection and hemodynamic quantification, the QVT offers an efficient, robust, and repeatable means to analyze cranial 4D flow MRI. This software is freely available at: https://github.com/uwmri/QVT.

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.

Four-dimensional flow magnetic resonance imaging for assessment of hemodynamic changes in the portal venous system before and after balloon-occluded retrograde transvenous obliteration: a pilot feasibility study

Background

The effectiveness of four-dimensional (4D) flow magnetic resonance imaging (MRI) for assessing hemodynamic changes before and after balloon-occluded retrograde transvenous obliteration (BRTO) remains unclear.

Purpose

To evaluate the feasibility of 4D flow MRI for assessing hemodynamic changes in the portal venous system before and after BRTO.

Material and Methods

We included 10 patients (7 men, 3 women; mean age = 67 years) with liver cirrhosis who had a high risk of gastric variceal bleeding or hepatic encephalopathy. Non-contrast 4D flow MRI of the upper abdomen was performed before and after BRTO. In addition, we compared the blood flow rates in the portal vein (PV), superior mesenteric vein (SMV), splenic vein (SV), left renal vein, and inferior vena cava before and after BRTO. Moreover, the flow directions of the SMV and SV before and after BRTO were assessed using both portography and 4D flow MRI.

Results

There was a significant post-BRTO increase in the blood flow rate in the PV and SV ( P < 0.05). There was no significant post-BRTO change in the blood flow rates in the SMV, inferior vena cava, and left renal vein. In four patients, portography confirmed that hepatofugal flow in the SV and SMV changed to hepatopetal flow after BRTO. Moreover, 4D flow MRI correctly assessed the flow directions in the SMV and SV in 70%–100% of the patients.

Conclusion

4D flow MRI can be used to detect hemodynamic changes in the portal venous system before and after BRTO.

5-Year Associations among Cerebral Arterial Pulsatility, Perivascular Space Dilation, and White Matter Lesions

OBJECTIVE

High cerebral arterial pulsatility index (PI), white matter lesions (WMLs), enlarged perivascular spaces (PVSs), and lacunar infarcts are common findings in the elderly population, and considered indicators of small vessel disease (SVD). Here, we investigate the potential temporal ordering among these variables, with emphasis on determining whether high PI is an early or delayed manifestation of SVD.

METHODS

In a population-based cohort, 4D flow MRI data for cerebral arterial pulsatility was collected for 159 participants at baseline (age 64-68), and for 122 participants at follow-up 5 years later. Structural MRI was used for WML and PVS segmentation, and lacune identification. Linear mixed-effects (LME) models were used to model longitudinal changes testing for pairwise associations, and latent change score (LCS) models to model multiple relationships among variables simultaneously.

RESULTS

Longitudinal 5-year increases were found for WML, PVS, and PI. Cerebral arterial PI at baseline did not predict changes in WML or PVS volume. However, WML and PVS volume at baseline predicted 5-year increases in PI. This was shown for PI increases in relation to baseline WML and PVS volumes using LME models (R  ≥  0.24; p < 0.02 and R  ≥  0.23; p < 0.03, respectively) and LCS models ( β  = 0.28; p = 0.015 and β  = 0.28; p = 0.009, respectively). Lacunes at baseline were unrelated to PI.

INTERPRETATION

In healthy older adults, indicators of SVD are related in a lead-lag fashion, in which the expression of WML and PVS precedes increases in cerebral arterial PI. Hence, we propose that elevated PI is a relatively late manifestation, rather than a risk factor, for cerebral SVD. ANN NEUROL 2022;92:871-881.

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.

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.

Hemodynamic Disturbances in Posterior Circulation Stroke: 4D Flow Magnetic Resonance Imaging Added to Computed Tomography Angiography

Objective: A clinically feasible, non-invasive method to quantify blood flow, hemodynamics, and collateral flow in the vertebrobasilar arterial tree is missing. The objective of this study was to evaluate the feasibility of quantifying blood flow and blood flow patterns using 4D flow magnetic resonance imaging (MRI) in consecutive patients after an ischemic stroke in the posterior circulation. We also explore if 4D-flow, analyzed in conjunction with computed tomography angiography (CTA), has potential as a diagnostic tool in posterior circulation stroke. Methods: Twenty-five patients (mean age 62 years; eight women) with acute ischemic stroke in the posterior circulation were investigated. At admission, all patients were examined with CTA followed by MRI (4D flow MRI and diffusion-weighted sequences) at median 4 days after the presenting event. Based on the classification of Caplan, patients were divided into proximal/middle (n = 16) and distal territory infarcts (n = 9). Absolute and relative blood flow rates were calculated for internal carotid arteries (ICA), vertebral arteries (VA), basilar artery (BA), posterior cerebral arteries (P1 and P2), and the posterior communicating arteries (Pcom). In a control group consisting of healthy elderly, the 90th and 10th percentiles of flow were calculated in order to define normal, increased, or decreased blood flow in each artery. "Major hemodynamic disturbance" was defined as low BA flow and either low P2 flow or high Pcom flow. Various minor hemodynamic disturbances were also defined. Blood flow rates were compared between groups. In addition, a comprehensive analysis of each patient's blood flow profile was performed by assessing relative blood flow rates in each artery in conjunction with findings from CTA. Results: There was no difference in total cerebral blood flow between patients and controls [604 ± 117 ml/min vs. 587 ± 169 ml/min (mean ± SD), p = 0.39] or in total inflow to the posterior circulation (i.e., the sum of total VA and Pcom flows, 159 ± 63 ml/min vs. 164 ± 52 ml/min, p = 0.98). In individual arteries, there were no significant differences between patients and controls in absolute or relative flow. However, patients had larger interindividual relative flow variance in BA, P1, and P2 (p = 0.01, <0.01, and 0.02, respectively). Out of the 16 patients that had proximal/middle territory infarcts, nine had CTA findings in VA and/or BA generating five with major hemodynamic disturbance identified with 4D flow MRI. For those without CTA findings, seven had no or minor 4D flow MRI hemodynamic disturbance. Among nine patients with distal territory infarcts, one had major hemodynamic disturbances, while the remaining had minor disturbances. Conclusion: 4D flow MRI contributed to the identification of the patients who had major hemodynamic disturbances from the vascular pathologies revealed on CTA. We thus conclude that 4D flow MRI could add valuable hemodynamic information when used in conjunction with CTA.

Assessing cerebral arterial pulse wave velocity using 4D flow MRI

Intracranial arterial stiffening is a potential early marker of emerging cerebrovascular dysfunction and could be mechanistically involved in disease processes detrimental to brain function via several pathways. A prominent consequence of arterial wall stiffening is the increased velocity at which the systolic pressure pulse wave propagates through the vasculature. Previous non-invasive measurements of the pulse wave propagation have been performed on the aorta or extracranial arteries with results linking increased pulse wave velocity to brain pathology. However, there is a lack of intracranial "target-organ" measurements. Here we present a 4D flow MRI method to estimate pulse wave velocity in the intracranial vascular tree. The method utilizes the full detectable branching structure of the cerebral vascular tree in an optimization framework that exploits small temporal shifts that exists between waveforms sampled at varying depths in the vasculature. The method is shown to be stable in an internal consistency test, and of sufficient sensitivity to robustly detect age-related increases in intracranial pulse wave velocity.

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.

Blood Flow Velocity Pulsatility and Arterial Diameter Pulsatility Measurements of the Intracranial Arteries Using 4D PC-MRI

4D phase contrast magnetic resonance imaging (PC-MRI) allows for the visualization and quantification of the cerebral blood flow. A drawback of software that is used to quantify the cerebral blood flow is that it oftentimes assumes a static arterial luminal area over the cardiac cycle. Quantifying the lumen area pulsatility index (aPI), i.e. the change in lumen area due to an increase in distending pressure over the cardiac cycle, can provide insight in the stiffness of the arteries. Arterial stiffness has received increased attention as a predictor in the development of cerebrovascular disease. In this study, we introduce software that allows for measurement of the aPI as well as the blood flow velocity pulsatility index (vPI) from 4D PC-MRI. The internal carotid arteries of seven volunteers were imaged using 7 T MRI. The aPI and vPI measurements from 4D PC-MRI were validated against measurements from 2D PC-MRI at two levels of the internal carotid arteries (C3 and C7). The aPI and vPI computed from 4D PC-MRI were comparable to those measured from 2D PC-MRI (aPI: mean difference: 0.03 (limits of agreement: -0.14 - 0.23); vPI: 0.03 (-0.17-0.23)). The measured blood flow rate for the C3 and C7 segments was similar, indicating that our proposed software correctly captures the variation in arterial lumen area and blood flow velocity that exists along the distal end of the carotid artery. Our software may potentially aid in identifying changes in arterial stiffness of the intracranial arteries caused by pathological changes to the vessel wall.

Exercise-induced fluid shifts are distinct to exercise mode and intensity: a comparison of blood flow-restricted and free-flow resistance exercise

MRI can provide fundamental tools in decoding physiological stressors stimulated by training paradigms. Acute physiological changes induced by three diverse exercise protocols known to elicit similar levels of muscle hypertrophy were evaluated using muscle functional magnetic resonance imaging (mfMRI). The study was a cross-over study with participants (n = 10) performing three acute unilateral knee extensor exercise protocols to failure and a work matched control exercise protocol. Participants were scanned after each exercise protocol; 70% 1 repetition maximum (RM) (FF70); 20% 1RM (FF20); 20% 1RM with blood flow restriction (BFR20); free-flow (FF) control work matched to BFR20 (FF20_WM). Post exercise mfMRI scans were used to obtain interleaved measures of muscle R2 (indicator of edema), R2' (indicator of deoxyhemoglobin), muscle cross sectional area (CSA) blood flow, and diffusion. Both BFR20 and FF20 exercise resulted in a larger acute decrease in R2, decrease in R2', and expansion of the extracellular compartment with slower rates of recovery. BFR20 caused greater acute increases in muscle CSA than FF20_WM and FF70. Only BFR20 caused acute increases in intracellular volume. Postexercise muscle blood flow was higher after FF70 and FF20 exercise than BFR20. Acute changes in mean diffusivity were similar across all exercise protocols. This study was able to differentiate the acute physiological responses between anabolic exercise protocols. Low-load exercise protocols, known to have relatively higher energy contributions from glycolysis at task failure, elicited a higher mfMRI response. Noninvasive mfMRI represents a promising tool for decoding mechanisms of anabolic adaptation in muscle.NEW & NOTEWORTHY Using muscle functional MRI (mfMRI), this study was able to differentiate the acute physiological responses following three established hypertrophic resistance exercise strategies. Low-load exercise protocols performed to failure, with or without blood flow restriction, resulted in larger changes in R₂ (i.e. greater T₂-shifts) with a slow rate of return to baseline indicative of myocellular fluid shifts. These data were cross evaluated with interleaved measures of macrovascular blood flow, water diffusion, muscle cross sectional area (i.e. acute macroscopic muscle swelling), and intracellular water fraction measured using MRI.

Quantification and mapping of cerebral hemodynamics before and after carotid endarterectomy, using four-dimensional flow magnetic resonance imaging

BACKGROUND

Carotid stenosis can profoundly affect cerebral hemodynamics, which cannot simply be inferred from the degree of stenosis. We quantified and mapped the distribution of the blood flow rate (BFR) in the cerebral arteries before and after carotid endarterectomy using four-dimensional (4D) phase-contrast (PC) magnetic resonance imaging (MRI).

METHODS

Nineteen patients (age, 71 ± 6 years; 2 women) with symptomatic carotid stenosis (≥50%) undergoing carotid endarterectomy (CEA) were investigated using 4D PC-MRI before and after surgery. The BFR was measured in 17 cerebral arteries and the ophthalmic arteries. Collateral recruitment through the anterior and posterior communicating arteries, ophthalmic arteries, and leptomeningeal arteries was quantified. BFR laterality was significantly different between the paired contralateral and ipsilateral arteries. Subgroups were defined according to the presence of collateral recruitment.

RESULTS

The total cerebral blood flow had increased by 15% (P < .01) after CEA. Before CEA, laterality was seen in the internal carotid artery, anterior cerebral artery, and middle cerebral artery (MCA). On the ipsilateral side, an increased BFR was found after CEA in the internal carotid artery (246 ± 62 mL/min vs 135 ± 80 mL/min; P < .001), anterior cerebral artery (87 ± mL/min vs 38 ± 58 mL/min; P < .01), and MCA (149 ± 43 mL/min vs 119 ± 34 mL/min; P < .01), resulting in a postoperative BFR distribution without signs of laterality. In the nine patients with preoperatively recruited collaterals, BFR laterality was found in the MCA before, but not after, CEA (P < .01). This laterality was not found in the 10 patients without collateral recruitment (P = .2). The degree of stenosis did not differ between the groups with and without collateral recruitment (P = .85).

CONCLUSIONS

Using 4D PC-MRI, we have presented a comprehensive and noninvasive method to evaluate the cerebral hemodynamics due to carotid stenosis before and after CEA. MCA laterality, seen in the patients with collateral recruitment before CEA, pointed toward a hemodynamic disturbance in MCA territory for those patients. This methodologic advancement provides an insight into the pathophysiology of cerebral hemodynamics in patients with carotid stenosis.

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.

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 (ΔP_MCA) 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 ΔP_MCA was calculated from CFD, where patient-specific flow rates were applied at all input- and output boundaries. The ΔP_MCA 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 ΔP_MCA and ΔQ_ICA (r = 0.85, p<0.001) and ΔQ_ACA (r = 0.71, p<0.001), respectively. The correlation to ΔQ_MCA 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.

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.

Semiautomatic cerebrovascular territory mapping based on dynamic ASL MR angiography without vessel-encoded labeling

PURPOSE

Characterizing vessel territories can provide crucial information for evaluation of cerebrovascular disorders. In this study, we present a novel postprocessing pipeline for vascular territorial imaging of cerebral arteries based on a noncontrast enhanced time-resolved 4D magnetic resonance angiography (MRA).

METHODS

Eight healthy participants, 1 Moyamoya patient, and 1 arteriovenous malformations patient were recruited. Territorial segmentation and relative blood flow rate calculations of cerebral arteries including left and right middle cerebral arteries and left and right posterior cerebral arteries were carried out based on the 4D MRA-derived arterial arrival time maps of intracranial vessels.

RESULTS

Among healthy young subjects, the average relative blood flow rate values corresponding to left and right middle cerebral arteries and left and right posterior cerebral arteries were 35.9 ± 5.9%, 32.9 ± 7.5%, 15.4 ± 3.8%, and 15.9 ± 2.5%, respectively. Excellent agreement was observed between relative blood flow rate values obtained from the proposed 4D MRA-based method and reference 2D phase contrast MRI. Abnormal cerebral circulations were visualized and quantified on both patients using the developed technique.

CONCLUSION

The vascular territorial imaging technique developed in this study allowed for the generation of territorial maps with user-defined level of details within a clinically feasible scan time, and as such may provide useful information to assess cerebral circulation balance in different pathologies.

Characterizing pulsatility in distal cerebral arteries using 4D flow MRI

Recent reports have suggested that age-related arterial stiffening and excessive cerebral arterial pulsatility cause blood-brain barrier breakdown, brain atrophy and cognitive decline. This has spurred interest in developing non-invasive methods to measure pulsatility in distal vessels, closer to the cerebral microcirculation. Here, we report a method based on four-dimensional (4D) flow MRI to estimate a global composite flow waveform of distal cerebral arteries. The method is based on finding and sampling arterial waveforms from thousands of cross sections in numerous small vessels of the brain, originating from cerebral cortical arteries. We demonstrate agreement with internal and external reference methods and show the ability to capture significant increases in distal cerebral arterial pulsatility as a function of age. The proposed approach can be used to advance our understanding regarding excessive arterial pulsatility as a potential trigger of cognitive decline and dementia.

Physiological responses of human skeletal muscle to acute blood flow restricted exercise assessed by multimodal MRI

Important physiological quantities for investigating muscle hypertrophy include blood oxygenation, cell swelling, and changes in blood flow. The purpose of this study was to compare the acute changes of these parameters in human skeletal muscle induced by low-load (20% 1-RM) blood flow-restricted (BFR-20) knee extensor exercise compared with free-flow work-matched (FF-20_WM) and free-flow 50% 1-RM (FF-50) knee extensor exercise using multimodal magnetic resonance imaging (MRI). Subjects (n = 11) completed acute exercise sessions for each exercise mode in an MRI scanner, where interleaved measures of muscle R₂ (indicator of edema), [Formula: see text] (indicator of deoxyhemoglobin), macrovascular blood flow, and diffusion were performed before, between sets, and after the final set for each exercise protocol. BFR-20 exercise resulted in larger acute decreases in R₂ and greater increases in cross-sectional area than FF-20_WM and FF-50 (P < 0.01). Blood oxygenation decreased between sets during BFR-20, as indicated by a 13.6% increase in [Formula: see text] values (P < 0.01)), whereas they remained unchanged for FF-20_WM and decreased during FF-50 exercise. Quadriceps blood flow between sets was highest for the heavier load (FF-50), averaging 305 mL/min, and lowest for BFR-20 at 123 ± 73 mL/min until post-exercise cuff release, where blood flow rates in BFR-20 exceeded both FF protocols (P < 0.01). Acute changes in diffusion rates were similar for all exercise protocols. This study was able to differentiate the acute exercise response of selected physiological factors associated with skeletal muscle hypertrophy. Marked differences in these parameters were found to exist between BFR and FF exercise conditions, which contribute to explain the anabolic potential of low-load blood flow restricted muscle exercise.NEW & NOTEWORTHY Acute changes in blood flow, diffusion, blood oxygenation, cross-sectional area, and the "T₂ shift" are evaluated in human skeletal muscle in response to blood flow-restricted (BFR) and conventional free-flow knee extensor exercise performed in an MRI scanner. The acute physiological response to exercise was dependent on the magnitude of load and the application of BFR. Physiological variables changed markedly and established a steady state rapidly after the first of four exercise sets.

Blood Flow Lateralization and Collateral Compensatory Mechanisms in Patients With Carotid Artery Stenosis

Background and Purpose- Four-dimensional phase-contrast magnetic resonance imaging enables quantification of blood flow rate (BFR; mL/min) in multiple cerebral arteries simultaneously, making it a promising technique for hemodynamic investigation in patients with stroke. The aim of this study was to quantify the hemodynamic disturbance and the compensatory pattern of collateral flow in patients with symptomatic carotid stenosis. Methods- Thirty-eight patients (mean, 72 years; 27 men) with symptomatic carotid stenosis (≥50%) or occlusion were investigated using 4-dimensional phase-contrast magnetic resonance imaging. For each patient, BFR was measured in 19 arteries/locations. The ipsilateral side to the symptomatic carotid stenosis was compared with the contralateral side. Results- Internal carotid artery BFR was lower on the ipsilateral side (134±87 versus 261±95 mL/min; P<0.001). BFR in anterior cerebral artery (A1 segment) was lower on ipsilateral side (35±58 versus 119±72 mL/min; P<0.001). Anterior cerebral artery territory bilaterally was primarily supplied by contralateral internal carotid artery. The ipsilateral internal carotid artery mainly supplied the ipsilateral middle cerebral artery (MCA) territory. MCA was also supplied by a reversed BFR found in the ophthalmic and the posterior communicating artery routes on the ipsilateral side (-5±28 versus 10±28 mL/min, P=0.001, and -2±12 versus 6±6 mL/min, P=0.03, respectively). Despite these compensations, BFR in MCA was lower on the ipsilateral side, and this laterality was more pronounced in patients with severe carotid stenosis (≥70%). Although comparing ipsilateral MCA BFR between stenosis groups (<70% and ≥70%), there was no difference ( P=0.95). Conclusions- With a novel approach using 4-dimensional phase-contrast magnetic resonance imaging, we could simultaneously quantify and rank the importance of collateral routes in patients with carotid stenosis. An important observation was that contralateral internal carotid artery mainly secured the bilateral anterior cerebral artery territory. Because of the collateral recruitment, compromised BFR in MCA is not necessarily related to the degree of carotid stenosis. These findings highlight the importance of simultaneous investigation of the hemodynamics of the entire cerebral arterial tree.

Cerebral arterial pulsatility is associated with features of small vessel disease in patients with acute stroke and TIA: a 4D flow MRI study

Cerebral small vessel disease (SVD) is a major cause of stroke and cognitive impairment. However, the underlying mechanisms behind SVD are still poorly understood. High cerebral arterial pulsatility has been suggested as a possible cause of SVD. In population studies, arterial pulsatility has been linked to white matter hyperintensities (WMH), cerebral atrophy, and cognitive impairment, all features of SVD. In stroke, pulsatility data are scarce and contradictory. The aim of this study was to investigate the relationship between arterial pulsatility and SVD in stroke patients. With a cross-sectional design, 89 patients with acute ischemic stroke or TIA were examined with MRI. A neuropsychological assessment was performed 1 year later. Using 4D flow MRI, pulsatile indices (PI) were calculated for the internal carotid artery (ICA) and middle cerebral artery (M1, M3). Flow volume pulsatility (FVP), a measure corresponding to the cyclic expansion of the arterial tree, was calculated for the same locations. These parameters were assessed for associations with WMH volume, brain volume and cognitive function. ICA-FVP was associated with WMH volume (β = 1.67, 95% CI: [0.1, 3.24], p = 0.037). M1-PI and M1-FVP were associated with decreasing cognitive function (β = − 4.4, 95% CI: [− 7.7, − 1.1], p = 0.009 and β = − 13.15, 95% CI: [− 24.26, − 2.04], p = 0.02 respectively). In summary, this supports an association between arterial pulsatility and SVD in stroke patients, and provides a potential target for further research and preventative treatment. FVP may become a useful biomarker for assessing pulsatile stress with PCMRI and 4D flow MRI.

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.

Segmentation-Based Blood Flow Parameter Refinement in Cerebrovascular Structures Using 4-D Arterial Spin Labeling MRA

OBJECTIVE

Cerebrovascular diseases are one of the main global causes of death and disability in the adult population. The preferred imaging modality for the diagnostic routine is digital subtraction angiography, an invasive modality. Time-resolved three-dimensional arterial spin labeling magnetic resonance angiography (4D ASL MRA) is an alternative non-invasive modality, which captures morphological and blood flow data of the cerebrovascular system, with high spatial and temporal resolution. This work proposes advanced medical image processing methods that extract the anatomical and hemodynamic information contained in 4D ASL MRA datasets.

METHODS

A previously published segmentation method, which uses blood flow data to improve its accuracy, is extended to estimate blood flow parameters by fitting a mathematical model to the measured vascular signal. The estimated values are then refined using regression techniques within the cerebrovascular segmentation. The proposed method was evaluated using fifteen 4D ASL MRA phantoms, with ground-truth morphological and hemodynamic data, fifteen 4D ASL MRA datasets acquired from healthy volunteers, and two 4D ASL MRA datasets from patients with a stenosis.

RESULTS

The proposed method reached an average Dice similarity coefficient of 0.957 and 0.938 in the phantom and real dataset segmentation evaluations, respectively. The estimated blood flow parameter values are more similar to the ground-truth values after the refinement step, when using phantoms. A qualitative analysis showed that the refined blood flow estimation is more realistic compared to the raw hemodynamic parameters.

CONCLUSION

The proposed method can provide accurate segmentations and blood flow parameter estimations in the cerebrovascular system using 4D ASL MRA datasets.

SIGNIFICANCE

The information obtained with the proposed method can help clinicians and researchers to study the cerebrovascular system non-invasively.

Composite MRA: statistical approach to generate an MR angiogram from multiple contrasts

PURPOSE

To develop a method to use information from multiple MRI contrasts to produce a composite angiogram with reduced sequence-specific artifacts and improved vessel depiction.

METHODS

Bayesian posterior vessel probability was determined as a function of black blood (BB), contrast enhanced angiography (CE-MRA), and phase-contrast MRA (PC-MRA) intensities from training subjects (N = 4). To generate composite angiogram in evaluation subjects (N = 12), the voxel-wise vessel probabilities were weighted with a confidence measure and combined as a weighted product to yield angiogram intensity. For 23 internal carotid artery (ICA) segments (N = 23) from evaluation subjects, segmentation accuracy of composite MRA was evaluated and compared against CE-MRA using dice similarity coefficient (DSC).

RESULTS

The composite MRA suppressed venous contaminations in CE-MRA, reduced flow artifacts, and velocity aliasing seen in PC-MRA and removed signal ambiguities in BB images. For ICA segmentations, the composite MRA improved segmentation over CE-MRA per DSC (0.908 ± 0.037 vs. 0.765 ± 0.079). Compared with CE-MRA, the composite MRA showed conservative changes in vessel appearance to small threshold changes. However, small vessels that are sensitive to registration errors or visible only weakly in CE-MRA were susceptible to poor depiction in composite MRA.

CONCLUSION

By dynamically weighting vessel information from multiple contrasts and extracting their complementary information, the composite MRA produces reduced sequence-specific artifacts and improved vessel contrast. It is a promising technique for semi-automatic segmentation of vessels that are hard to segment because of artifacts.