Spatial correlations between MRI-derived wall shear stress and vessel wall thickness in the carotid bifurcation

The Correlation Between Wall Shear Stress and Plaque Composition in Advanced Human Carotid Atherosclerosis

The role of wall shear stress (WSS) in atherosclerotic plaque development is evident, but the relation between WSS and plaque composition in advanced atherosclerosis, potentially resulting in plaque destabilization, is a topic of discussion. Using our previously developed image registration pipeline, we investigated the relation between two WSS metrics, time-averaged WSS (TAWSS) and the oscillatory shear index (OSI), and the local histologically determined plaque composition in a set of advanced human carotid plaques. Our dataset of 11 carotid endarterectomy samples yielded 87 histological cross-sections, which yielded 511 radial bins for analysis. Both TAWSS and OSI values were subdivided into patient-specific low, mid, and high tertiles. This cross-sectional study shows that necrotic core (NC) size and macrophage area are significantly larger in areas exposed to high TAWSS or low OSI. Local TAWSS and OSI tertile values were generally inversely related, as described in the literature, but other combinations were also found. Investigating the relation between plaque vulnerability features and different combinations of TAWSS and OSI tertile values revealed a significantly larger cap thickness in areas exposed to both low TAWSS and low OSI. In conclusion, our study confirmed previous findings, correlating high TAWSS to larger macrophage areas and necrotic core sizes. In addition, our study demonstrated new relations, correlating low OSI to larger macrophage areas, and a combination of low TAWSS and low OSI to larger cap thickness.

Exploring the Relationships Between Hemodynamic Stresses in the Carotid Arteries

Background: Atherosclerosis manifests as a focal disease, often affecting areas with complex hemodynamics such as the carotid bifurcation. The magnitude and regularity of the hemodynamic shear stresses acting on the vessel wall are thought to generate risk patterns unique to each patient and play a role in the pathogenesis of atherosclerosis. The involvement of different expressions of shear stress in the pathogenesis of carotid atherosclerosis highlights the need to characterize and compare the differential impact of the various expressions of shear stress in the atherosclerotic carotid bifurcation. Therefore, the aim of this study is to characterize and compare hemodynamic wall shear stresses (WSS) in the carotid arteries of subjects with asymptomatic atherosclerotic plaques. Shear stresses were also compared against vessel diameter and bifurcation angle to examine the relationships with the geometry of the carotid bifurcation. Methods: 4D Flow MRI and contrast-enhanced MRA data were acquired for 245 subjects with atherosclerotic plaques of at least 2.7 mm in conjunction with the Swedish CArdioPulmonary bioImage Study (SCAPIS). Following automatic segmentation and geometric analysis, time-resolved WSS and near-wall turbulent kinetic energy (nwTKE) were derived from the 4D Flow data. Whole-cycle parameters including time-averaged WSS and nwTKE, and the oscillatory shear index (OSI) were calculated. Pairwise Spearman rank-correlation analyses were used to investigate relationships among the hemodynamic as well as geometric parameters. Results: One hundred and seventy nine subjects were successfully segmented using automated tools and subsequently geometric and hemodynamic analyses were performed. Temporally resolved WSS and nwTKE were strongly correlated, ρ = 0.64. Cycle-averaged WSS and nwTKE were moderately correlated, ρ = 0.57. Cycle-average nwTKE was weakly correlated to OSI (ρ = -0.273), revealing that nwTKE provides information about disturbed flow on the vessel wall that OSI does not. In this cohort, there was large inter-individual variation for both WSS and nwTKE. Both WSS and nwTKE varied most within the external carotid artery. WSS, nwTKE, and OSI were weakly correlated to vessel diameter and bifurcation angle. Conclusion: The turbulent and mean component of WSS were examined together in vivo for the first time, and a strong correlation was found between them. nwTKE presents the opportunity to quantify turbulent wall stresses in vivo and gain insight into the effects of disturbed flow on the vessel wall. Neither vessel diameter nor bifurcation angle were found to be strongly correlated to the turbulent or mean component of WSS in this cohort.

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.

Carotid geometry is an independent predictor of wall thickness - a 3D cardiovascular magnetic resonance study in patients with high cardiovascular risk

BACKGROUND

The posterior wall of the proximal internal carotid artery (ICA) is the predilection site for the development of stenosis. To optimally prevent stroke, identification of new risk factors for plaque progression is of high interest. Therefore, we studied the impact of carotid geometry and wall shear stress on cardiovascular magnetic resonance (CMR)-depicted wall thickness in the ICA of patients with high cardiovascular disease risk.

METHODS

One hundred twenty-one consecutive patients ≥50 years with hypertension, ≥1 additional cardiovascular risk factor and ICA plaque ≥1.5 mm thickness and < 50% stenosis were prospectively included. High-resolution 3D-multi-contrast (time of flight, T1, T2, proton density) and 4D flow CMR were performed for the assessment of morphological (bifurcation angle, ICA/common carotid artery (CCA) diameter ratio, tortuosity, and wall thickness) and hemodynamic parameters (absolute/systolic wall shear stress (WSS), oscillatory shear index (OSI)) in 242 carotid bifurcations.

RESULTS

We found lower absolute/systolic WSS, higher OSI and increased wall thickness in the posterior compared to the anterior wall of the ICA bulb (p < 0.001), whereas this correlation disappeared in ≥10% stenosis. Higher carotid tortuosity (regression coefficient = 0.764; p < 0.001) and lower ICA/CCA diameter ratio (regression coefficient = - 0.302; p < 0.001) were independent predictors of increased wall thickness even after adjustment for cardiovascular risk factors. This association was not found for bifurcation angle, WSS or OSI in multivariate regression analysis.

CONCLUSIONS

High carotid tortuosity and low ICA diameter were independent predictors for wall thickness of the ICA bulb in this cross-sectional study, whereas this association was not present for WSS or OSI. Thus, consideration of geometric parameters of the carotid bifurcation could be helpful to identify patients at increased risk of carotid plaque generation. However, this association and the potential benefit of WSS measurement need to be further explored in a longitudinal study.

4D Flow with MRI

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

Bilateral Multiband 4D Flow MRI of the Carotid Arteries at 7T

PURPOSE

Simultaneous multislab (SMSb) 4D flow MRI was developed and implemented at 7T for accelerated acquisition of the 3D blood velocity vector field in both carotid bifurcations.

METHODS

SMSb was applied to 4D flow to acquire blood velocities in both carotid bifurcations in sagittal orientation using a local transmit/receive coil at 7T. B 1 + transmit efficiency was optimized by B 1 + shimming. SMSb 4D flow was obtained in 8 healthy subjects in single-band (SB) and multiband (MB) fashion. Additionally, MB data were retrospectively undersampled to simulate GRAPPA R = 2 (MB2_GRAPPA2), and both SB datasets were added to form an artificial MB dataset (SumSB). The band separation performance was quantified by signal leakage. Peak velocity and total flow values were calculated and compared to SB via intraclass correlation analysis (ICC).

RESULTS

Clean slab separation was achieved yielding a mean signal leakage of 13% above the mean SB noise level. Mean total flow for MB2, SumSB, and MB_GRAPPA2 deviated less than 9% from the SB values. Peak velocities averaged over all vessels and subjects were 0.48 ± 0.11 m/s for SB, 0.47 ± 0.12 m/s for SumSB, 0.50 ± 0.13 m/s for MB2, and 0.53 ± 0.13 m/s for MB2_GRAPPA2. ICC revealed excellent absolute agreement and consistency of total flow for all methods compared to SB2. Peak velocity showed good to excellent agreement and consistency for SumSB and MB2 and MB2_GRAPPA2 method showed poor to excellent agreement and good to excellent consistency.

CONCLUSION

Simultaneous multislab 4D Flow MRI allows accurate quantification of total flow and peak velocity while reducing scan times.

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.

An MRI-based method to register patient-specific wall shear stress data to histology

Wall shear stress (WSS), the frictional force exerted on endothelial cells by blood flow, is hypothesised to influence atherosclerotic plaque growth and composition. We developed a methodology for image registration of MR and histology images of advanced human carotid plaques and corresponding WSS data, obtained by MRI and computational fluid dynamics. The image registration method requires four types of input images, in vivo MRI, ex vivo MRI, photographs of transversally sectioned plaque tissue and histology images. These images are transformed to a shared 3D image domain by applying a combination of rigid and non-rigid registration algorithms. Transformation matrices obtained from registration of these images are used to transform subject-specific WSS data to the shared 3D image domain as well. WSS values originating from the 3D WSS map are visualised in 2D on the corresponding lumen locations in the histological sections and divided into eight radial segments. In each radial segment, the correlation between WSS values and plaque composition based on histological parameters can be assessed. The registration method was successfully applied to two carotid endarterectomy specimens. The resulting matched contours from the imaging modalities had Hausdorff distances between 0.57 and 0.70 mm, which is in the order of magnitude of the in vivo MRI resolution. We simulated the effect of a mismatch in the rigid registration of imaging modalities on WSS results by relocating the WSS data with respect to the stack of histology images. A 0.6 mm relocation altered the mean WSS values projected on radial bins on average by 0.59 Pa, compared to the output of original registration. This mismatch of one image slice did not change the correlation between WSS and plaque thickness. In conclusion, we created a method to investigate correlations between WSS and plaque composition.