 Technical notes
 Open Access
 Published:
In vivo cardiovascular magnetic resonance of 2D vessel wall diffusion anisotropy in carotid arteries
Journal of Cardiovascular Magnetic Resonance volume 18, Article number: 81 (2017)
Abstract
Background
Diffusion weighted (DW) cardiovascular magnetic resonance (CMR) has shown great potential to discriminate between healthy and diseased vessel tissue by evaluating the apparent diffusion coefficient (ADC) along the arterial axis. Recently, ex vivo studies on porcine arteries utilizing diffusion tensor imaging (DTI) revealed a circumferential fiber orientation rather than an organization in axial direction, suggesting dominant diffusion perpendicular to the slice direction. In the present study, we propose a method to access tangential and radial diffusion of carotids in vivo by utilizing a pulse sequence that enables high resolution DW imaging in combination with a twodimensional (2D) diffusion gradient direction sampling scheme perpendicular to the longitudinal axis of the artery.
Methods
High resolution DTI of 12 healthy male volunteers (age: 25–60 years) was performed on one selected axial slice using a readout segmented EPI (rsEPI) sequence on a 3T MR scanner.
Results
It was found consistently for all 12 volunteers, that the tangential component as the principle direction of diffusion. Mean vessel wall fractional anisotropy (FA) values ranged from 0.7 for the youngest to 0.56 for the oldest participant. Linear regression analysis between the FA values and volunteers age revealed a highly significant (P < 0.01) linear relationship with an adjusted R^{2} of 0.52. In addition, a linear trend (P < 0.1) could be observed between radial diffusivity (RD) and age.
Conclusion
These results point to FA being a sensitive parameter able to capture changes in the vascular architecture with age. In detail, the data demonstrate a decrease in FA with advancing age indicating possible alterations of tissue microstructural integrity. Moreover, analyzing 2D diffusion tensor directions is sufficient and applicable in a clinical setup concerning the overall scan time.
Background
Cardiovascular diseases (CVDs) remain the global killer number one. Estimations from the world health organization (WHO) revealed, that about one third of global deaths were attributable to this life threatening diseases in 2012 [1]. Of these events, three quarters are associated with stroke and coronary heart disease. The most severe pathological condition is atherosclerosis (AS), affecting blood and heart vessels. AS is a chronic immunemediated inflammatory state of large and medium sized arteries that manifests in the first decade of life, with silent progression until a clinical event takes place [2]. The detection and diagnosis of AS at an early stage plays a key role in health care allowing for personalized interventions adjusted to the needs of patients.
Imaging modalities, such as ultrasound, computed tomography (CT), cardiovascular magnetic resonance (CMR) can assist in the detection and classification of morphological pathologies. However, these imaging techniques are applied mainly to determine the degree of luminal narrowing, which is an indispensable information, but a limited biomarker to determine the true impact on future cardiovascular events [3]. A more advanced imaging technique is multicontrast CMR which became a wellrecognized tool to study the complex composition of atherosclerotic plaques [4, 5]. The intrinsic soft tissue contrast and high spatial resolution makes CMR a well suited tool to delineate key constituents of plaque vulnerability.
Recently, diffusionweighted (DW) CMR has shown great potential to discriminate between normal and pathological vessel tissue by evaluating the apparent diffusion coefficient (ADC) [6–8] for a certain vascular segment. Given the fact that the spatial resolution is limited in DWI as a consequence of a dramatically increase in echo time when using conventional diffusion measurements, the diffusion perpendicular to an artery’s axis has never been imaged in vivo so far. The specific architecture of vessels allows for the investigation of diffusion anisotropy in the vessel wall by diffusion tensor imaging (DTI), which has already been shown in ex vivo studies where the fibrous structure of porcine aortas has successfully been visualized [9, 10].
Two major constraints prevent the application of DTI to visualize diffusion in the carotid artery wall in vivo. Firstly, the low spatial resolution when using a conventional singleshot EPI (ssEPI) and a significant loss of SNR at higher spatial resolution as a result of the prolonged echo time. Secondly, the scan time that dramatically increases with the number of diffusion sensitizing gradient directions.
In this work we propose a method to overcome these limitations in order to measure tangential and radial diffusion, evaluated from a twodimensional diffusion tensor, in the carotid artery wall in vivo. High spatial resolution was achieved by utilizing a readout segmented EPI (rsEPI) pulse sequence [11]. Due to the specific properties of this pulsesequence, sufficient SNR can be gained to allow for DTI analysis. Clinically acceptable scan time was assured by using a novel gradient direction scheme. Motivated by results from ex vivo DTI studies [9, 10] that revealed a tangential orientated principle diffusion direction and further supported by our ex vivo measurements, the usually used threedimensional gradient direction scheme was reduced to a twodimensional acquisition scheme that significantly decreased the number of needed gradient directions.
In this proof of concept we demonstrate the feasibility of measuring diffusion in the vessel wall of carotid arteries in vivo in a nondestructive and noninvasive fashion and provide initial results from 12 healthy male volunteers.
Methods
Twodimensional diffusion acquisition scheme
The basic idea to reduce the amount of diffusion sensitizing gradients is the reduction of a threedimensional diffusion measurement to a twodimensional diffusion measurement perpendicular to the artery’s axis (Fig. 1). This strategy is justified by the work from Flamini et al. [10] revealing a tangential fiber orientation that indicates the principle direction of diffusion. Given the fact that diffusion in zaxis is low compared to tangential and radial diffusion, this simplification provides the basis for in vivo measurements due to the significant reduction of diffusion sensitizing gradients needed for DTI. To demonstrate the validity of our approach 3D and 2D diffusion measurements on a porcine aorta were performed.
Preliminary work on porcine aorta
Ex vivo DWI of a slaughterhouse harvested porcine aorta was performed on one selected axial slice using rsEPI and both a threedimensional gradient direction schema with 30 directions (Multi Directional Diffusion Weighting, MDDW; acquisition time = ~50 min) and a twodimensional schema with 18 directions (acquisition time = ~30 min). Sequence parameters were as follows: FOV = 111 × 111 mm^{2}, matrix = 308 × 308, slice thickness = 10 mm, TR = 2500 ms, TE = 81 ms, GRAPPA = 2, number of readout segments = 31, bvalues = 0 and 1000 s/mm^{2}. To extract angle distributions of the principle 3DDTI tangential vector with respect to the longitudinal axis of the artery Eq. (1) was used.
To optimize imaging conditions, ex vivo specimens were immersed in a fluid (1,1,2Trichloro1,2,2trifluoroethane, Genetron® 113, SigmaAldrich) having no MR signal and to minimize magnetic susceptibility effects at airtissue boundaries [12, 13].
In vivo MR imaging
The 2D diffusion acquisition scheme was used for in vivo measurements requiring that the diffusion sensitizing directions, lying within a plane in space, are oriented perpendicular to the artery’s axis. This was achieved by defining the diffusion directions in the (p,r,s) coordinate system that is linked to the orientation of the field of view.
High resolution diffusion weighted imaging (DW) of 12 healthy male volunteers (25–60 years) was performed on one selected axial slice using rsEPI on a 3T whole body MR scanner (Prisma, Siemens Medical Solutions, Erlangen, Germany) with a 2 × 4 channel multifunctional coil (NORAS MRI products GmbH, Germany) and FOV = 189 × 189 mm^{2}, matrix = 346 × 346, slice thickness = 10 mm, TR = 2RR intervals, trigger delay = 50 ms to acquire images in the diastolic phase, TE = 93 ms, GRAPPA = 2, number of readout segments = 9, bvalues = 0, 200, 400, 600 s/mm^{2}, acquisition time = ~12 min. Eighteen gradient directions were defined on a hemicycle within a plane oriented perpendicular to the longitudinal axis of the carotid artery. To minimize artifacts, introduced by pulsatile vessel motion, peripheral pulsetriggering was performed using the systems’ internal physiological monitoring unit. To avoid chemical shift artifacts caused by periadventitial fat, a spectral fat saturation pulse was applied. A series of twodimensional (2D) angiographic axial vessel scout images (TurboFLASH: TR = 500 ms, TI = 280 ms, TE = 4.65 ms, flip angle = 15 °, FOV = 150 × 150 mm^{2}, matrix = 256 × 102, slice thickness = 5 mm, acquisition time = ~1 min) preceded the main task to locate the carotids and orient the slice perpendicular to the long axis of the artery. Measured DWI data with an inplane resolution of 0.55 × 0.55 mm^{2} were interpolated to 0.2 × 0.2 mm^{2} using zerofilling.
Registration of the individual DW images was achieved by applying publicly available elastix software [14, 15]. In detail, registration of the individual DW images was applied to the b0 image. The registration process included a rigid registration (inplane translation only) to correct for motion introduced by the relaxing back muscles during the investigation. In a next step the wall of the individual DW images was manually segmented to prevent misregistration to the dominant spine/neck structures. As a final processing step in order to correct for distortions introduced by more complex motion (eddy currents, breathing, swallowing and residual vessel motion) a Bspline transform was applied to register the delineated wall DW images [16, 17].
The study was approved by the institutional review board (IRB) as well as followed Declaration of Helsinki recommendations and informed consent was signed by each participant prior to the MR measurements.
Reproducibility of in vivo FA measurements
To test the robustness of repeated FA measurements using the above introduced in vivo high resolution DW sequence, the carotid artery of four male volunteers was repeatedly (in total four visits) measured and the coefficient of variation (CV), defined as the ratio of the standard deviation (σ) to the mean (x̅), was calculated.
Multiple bvalue selection for the in vivo DWI imaging experiment
The afore mentioned in vivo high resolution diffusion weighted imaging (DWI) parameter setup was used to investigate optimal bvalues on one selected axial carotid slice of a male volunteer. With all other sequence parameters held constant, the TE was set to 97 ms in order to allow for measurements with increasing bvalues from 0 to 1000 s/mm^{2} (in steps of 200 s/mm^{2}) for a certain diffusion gradient direction perpendicular to the slice direction.
Statistical analysis
To investigate the relationship between fractional anisotropy (FA), radial diffusivity (RD) and primary eigenvalue (λ_{1}) data with volunteer’s age a linear regression analysis was used to estimate parameters (coefficient of determination (R^{2}), slope and intercept) and assess linearity. In addition, 95 % confidence and prediction intervals for the linear regression were calculated. An Ftest for linear models was employed to compare variances and a probability value of P < 0.05 was considered as significant. Normal distribution was tested based on ShapiroWilk Normality Test. All analyses were performed using freely available Rproject software [18].
Diffusion tensor mathematics, FA, RD, MD and primary eigenvalue calculation
The 3D diffusion tensor described by a symmetric 3 × 3 matrix and its derived quantities can easily be reduced to the 2D case where the tensor is defined by a 2 × 2 matrix. Mathematical aspects are discussed in the following and are detailed described by Kingsley et al. [19–21]. A LevenbergMarquardt nonlinear least squares fitting algorithm was applied based on Eq. (2) to calculate diffusion coefficients (D, measured apparent diffusion coefficient) for every voxel for all measured gradient directions. Attenuation of the MR signal S(b) is introduced by varying the b value (200, 400, 600 s/mm^{2}) sensitizing it to water diffusion. S_{0} represents the reference signal without diffusion weighting (no diffusion gradient).
The 2D gradient direction sampling scheme enables the evaluation of the diffusion tensor d by solving Eq. (3):
where Y represents the calculated apparent diffusion coefficients (18 element row vector), H is derived from the normalized gradient components Eq. (4) and d is the tensor given as a threeelement column vector Eq. (5).
From tensor d corresponding eigenvalues (λ) and eigenvectors (ε) can be calculated on a pixel by pixel basis using eigenvalue decomposition. To study tangential and radial diffusion of the vessel wall 2D FA values were obtained using Eq. (6).
FA values close to zero mean no directed diffusion whereas anisotropy is given if one of the eigenvalues will be higher than the other. Calculation included the primary eigenvalue (λ_{1}), RD (λ_{2} for the twodimensional tensor) and the mean diffusivity (MD) which is the average of the eigenvalues.
Results
In vivo vessel wall diffusion weighted imaging
Table 1 summarizes the calculated mean diffusion components as well as the mean diffusivity (MD) and FA for the high resolution 2DDTI in vivo case of the 12 male volunteers. The tangential and radial diffusion components correspond to the eigenvalues λ_{1} and λ_{2}, respectively. Figure 1 outlines the anatomical region of interest (carotid artery) and the orientation of the diffusion gradient directions on a hemicycle within a plane perpendicular to the longitudinal axis of the carotid artery. The concept of tangential and radial diffusion as well as the twodimensional diffusion ellipsoid with eigenvalues (λ) and eigenvectors (ε) is delineated. In addition, Fig. 1 shows representative in vivo DWI from the carotid artery at a bvalue of 600 s/mm^{2} for the predefined set of diffusion directions (headless white arrows). It can be observed that the signal is suppressed in the direction of the diffusion sensitizing gradient, marked by the white arrowhead. Calculating the 2D diffusion tensor revealed consistently for all 12 volunteers that the largest eigenvalue is oriented tangential to the vessel wall while the second component therefore shows a radial orientation. This fact provides clear evidence for the tangential direction as the principle contribution for diffusion within the vessel wall as shown in Fig. 2a. Based on the eigenvalues a twodimensional fractional anisotropy map (FA) was generated (Fig. 2b). Additionally overlain is the direction of the principal DTI vectors denoted by the red lines.
Robustness of in vivo FA measurements
Table 2 summarizes repeated measurements (#1 to #4) of individual mean FA values of four male volunteers (R1 to R2). The robustness of the measurement is demonstrated by the coefficient of variation (CV) with values ranging from 2.5 to 5.4 %.
3D and 2D diffusion comparison in porcine
To confirm the tangential diffusion as the direction of the largest eigenvalue ex vivo high resolution 3DDTI and 2DDTI of a slaughterhouse harvested porcine aorta was performed using rsEPI. Figure 2c and d represents the 3D diffusion tensor as ellipsoids in a three dimensional color space (RGB), and the 2D diffusion tensor in a twodimensional color space (RG), respectively. Calculating the angle between the principle diffusion direction and the artery’s axis for the 3D measurement, the mean value of 1.8 ± 7.2 ° confirms the dominant orientation of diffusion perpendicular to the artery’s axis. The low diffusion in zdirection justifies the 2D diffusion measurement perpendicular to the artery’s axis allowing for in vivo measurements in a clinically feasible time due to the lower number of diffusion sensitizing gradient directions. In addition, Table 3 summarizes the calculated mean diffusion components for the high resolution 2DDTI (tangential and radial) and 3DDTI (tangential, radial and longitudinal) ex vivo case of the slaughterhouse porcine aorta. For the 2D case the diffusion components correspond to the eigenvalues λ_{1} and λ_{2}. For the 3D case the tangential component corresponds to the largest eigenvalue (λ_{1}). The radial and longitudinal components can either be λ_{2} or λ_{3} because both eigenvalues are in a comparable range of magnitude (Fig. 2c).
bvalue estimation
To estimate optimal bvalues and study the tissues’ anisotropic properties, an in vivo vessel wall ADC map was generated from multiple measurements with increasing bvalues (0–1000 s/mm^{2}) for a certain diffusion direction perpendicular to the slice direction (Fig. 3a). Given that a signal drop can be observed in the direction of the applied diffusion gradient direction and a signal enhancement perpendicular to it, a diffusion process perpendicular to the vessels’ longitudinal axis is assumable as exemplified in Fig. 3b. As shown in Fig. 3c, two pools of ADC values could be identified (mean values: 2.11e3 ± 0.23e3 and 1.27e3 ± 0.29e3 mm^{2}/s). In addition, outlined in Fig. 3d is the final selected set of four equidistant bvalue increments (0, 200, 400, 600 s/mm^{2}) and the maximum difference of both populations found at a bvalue of about 605 s/mm^{2}. In addition to the diffusion decay factor e ^{(−bD)}, signal loss is given by the vessel’s tissue T2 relaxation time and predefined echo time as expressed by the decay constant e ^{(−TE/T2)}. All these effects contribute to a signal drop and justify a distribution of selected bvalues between zero and the found maximum.
FA, RD and primary eigenvalue vs. age relationship
In Fig. 4a, b and c, mean vessel wall FA, RD and primary eigenvalue (λ_{1}) data for 12 male volunteers ranging from 0.7 for the youngest to 0.56 for the oldest are plotted as a function of age. The linear regression for the FA and age relationship was highly significant P < 0.01 with an adjusted R^{2} of 0.52 and the linear equation y = −0.0049x + 0.79. For the relationship between mean RD and age a linear trend P < 0.1 with an adjusted R^{2} of 0.18 was observed. No linear relationship P > 0.6 between the mean primary eigenvalues (λ_{1}) and age was found. This analysis indicates that FA, measured within the vessel wall, is a sensitive parameter able to capture changes in the vascular architecture due to aging.
Discussion
In this work, we demonstrated the feasibility of measuring diffusion anisotropy in human carotids in vivo in a nondestructive and noninvasive fashion. In contrast to previous works that aim for quantitative diffusion imaging in vivo (DW imaging, ADC values) along the slice direction (zdirection) our results suggests that tangential and radial diffusion is accessible by using a pulse sequence that enables high resolution DW imaging measurements in combination with a 2D gradient direction samplingscheme orientated perpendicular to the vessel’s longitudinal axis. Furthermore, our in vivo and ex vivo results corroborate the fact that the principle diffusion direction is aligned perpendicular to the long axis of the artery as demonstrated by ex vivo studies on porcine aortas [9, 10]. Moreover, fiber orientation organizes more circumferential and coherent with increasing transmural pressure as was shown in an ex vivo investigation on human arteries [22]. Although our ex vivo setup was pressure free, the circumferential orientation of diffusion remained highly preserved (1.8 ± 7.2 °) another clear evidence for a dominant diffusion process perpendicular to the long axis of the vessel in vivo. Quantitative diffusion imaging along zdirection by evaluating the ADC may spare major information to ascertain arterial microstructure properties or changes in case of pathological tissue remodeling. Moreover, ADC depends on subject’s orientation relative to the magnetic field gradients introducing errors like scanrescan variations or possible confound comparisons between individuals [23]. In contrast, DTI allows overcoming these issues by determining rotationally invariant parameters (e.g. eigenvalues λ) which values do not depend on subject’s orientation in the scanner reference frame. Hence, our framework implies that the combination of high resolution DW imaging with a special set of diffusion gradients aligned perpendicular to the long axis of the artery is essential to uncover important structural vessel wall information. Moreover, it was proofed that analyzing twodimensional diffusion tensor directions is sufficient and also applicable in a clinical setup concerning the overall scan time. The excellent agreement of the tangential and radial components from the 2D and 3DDTI analysis (Table 3) of the ex vivo porcine aorta further supports the reduction of the threedimensional tensor to a twodimensional one and the feasibility of our proposed method. Differences in the magnitude of mean FA values between 2D and 3D (Table 3) occur due to the missing third component (longitudinal diffusion). The reduction to a twodimensional tensor causes an increase of the calculated FA values.
Recent investigations of ADC values along the zdirection (slice direction) were done on healthy volunteers and patients but the relevance of aging was not considered [6–8]. In contrast, our data reveal a highly significant relationship between FA vessel wall values and volunteer’s age. In detail, the data demonstrate a decrease in FA with advancing age indicating possible alterations of tissue microstructural integrity. Furthermore, RD values tend to increase with age indicating a possible decrease in cellularity and a conversion process to more extracellular space. In contrast, the primary eigenvalue λ_{1} indicate a rather constant behavior with advancing age. These indicators taken together may provide complex information about ongoing processes on the vascular architecture. However, the reported aging effect is notable but lifestyle may also influence vascular wall structure and was not covered by this proof of concept. To investigate this process more closely, further experiments are needed including more subjects. Another interesting application of the proposed method might be to study the influence of plaque components on the diffusion process. Future work on vessel wall diffusion imaging may help to understand the formation of such a complex disease.
To achieve accurate estimates of the selfdiffusion tensor the sensitivity to motioninduced phase errors has to be reduced to a minimum. Motion is a severe factor in diffusion imaging in general and specifically when trying to access diffusion in the carotid wall [6, 8]. Several strategies used in this work directly address this problem. Firstly the suggested RESOLVE sequence inherently corrects for motion artifacts using a 2D navigator echo (correction for nonlinear phase errors introduced by nonrigid body motion) [11, 24–26]. Secondly, pulse triggering accounts for the pulsatile motion of the arterial wall and thirdly, a coregistration procedure including translation followed bspline based registration minimizes the distance of subsequent images. With this strategy no subject had to be excluded from the analysis due to motion. The repeatability of mean FA measurements (Table 2; CV: 2.5–5.4 %) clearly indicates the feasibility and robustness to evaluate the twodimensional tensor in vivo.
Conventional singleshot EPI’s are fast sequences that enable the acquisition of a slice within a single repetition time but suffer from susceptibility artifacts and T2* blurring due to the prolonged time required to fill the kspace, low resolution and significant loss of SNR at higher resolution. Diffusion tensor imaging is inherently a technique suffering from low SNR. High spatial resolution cannot be achieved using such conventional ssEPI due to the signal loss as a result of the prolonged echo time. Hence the use of rsEPI is a prerequisite for high spatial resolution. It can be seen from the representative diffusion weighted images of Fig. 1 that the vessel wall can clearly be delineated and the signal is sufficient to calculate the diffusion tensor. The specific properties of the suggested RESOLVE sequence (short EPI echo spacing; 2D navigator based phase correction; parallel imaging using GRAPPA [27]) reduces the sensitivity to image distortions and motion induced artifacts and therefore make the sequence a well suited tool to analyze the twodimensional tensor in vivo with sufficient SNR.
Conclusions
In conclusion, we present a novel method to access tangential and radial diffusion of carotids in vivo by utilizing a pulse sequence that enables high resolution DW imaging in combination with a 2D diffusion gradient direction sampling scheme perpendicular to the longitudinal axis of the artery. In contrast to standard DW methods evaluating the ADC for a certain vascular segment along the artery’s axis our method allows to image the water diffusion tensor within the vessel wall with a high spatial resolution. Ex vivo studies on porcine arteries suggested a dominant circumferential fiber orientation rather than an organization in axial direction supporting our 2D approach. Furthermore, our results demonstrate a decrease in fractional anisotropy (FA) with advancing age indicating FA as a sensitive parameter able to capture changes in the vascular architecture with age. Moreover, analyzing 2D diffusion tensor directions is sufficient and applicable in a clinical setup concerning the overall scan time. Our data demonstrate for the first time the feasibility of measuring diffusion anisotropy in human carotids in vivo in a nondestructive and noninvasive fashion. Future work will investigate the influence of plaque components on the diffusion process, which may help to understand the formation of a complex disease.
Abbreviations
 2D:

Twodimensional
 3D:

Threedimensional
 ADC:

Apparent diffusion coefficient
 AS:

Atherosclerosis
 CMR:

Cardiovascular magnetic resonance
 CVDs:

Cardiovascular diseases
 DTI:

Diffusion tensor imaging
 DW:

Diffusion weighted
 FA:

Fractional anisotropy
 FLASH:

Fast low angle shot
 FOV:

Field of view
 GRAPPA:

Generalized autocalibrating partially parallel acquisition
 MDDW:

Multi directional diffusion weighting
 RD:

Radial diffusivity
 rsEPI:

Readout segmented echo planar imaging
 ssEPI:

Singleshot echo planar imaging
 TE:

Echo time
 TI:

Inversion time
 TR:

Repetition time
References
 1.
Health Organization World WHO: Global Status Report on Noncommunicable Diseases 2014. World Health Organization; 2015. ISBN Number: 9789241564854.
 2.
Lusis AJ. Atherosclerosis. Nature. 2000;407:233–41.
 3.
Saam T, Underhill HR, Chu B, Takaya N, Cai J, Polissar NL, Yuan C, Hatsukami TS. Prevalence of American Heart Association type VI carotid atherosclerotic lesions identified by magnetic resonance imaging for different levels of stenosis as measured by duplex ultrasound. J Am Coll Cardiol. 2008;51:1014–21.
 4.
Cai JM, Hatsukami TS, Ferguson MS, Small R, Polissar NL, Yuan C. Classification of human carotid atherosclerotic lesions with in vivo multicontrast magnetic resonance imaging. Circulation. 2002;106:1368–73.
 5.
Sun J, Zhao XQ, Balu N, Hippe DS, Hatsukami TS, Isquith DA, Yamada K, Neradilek MB, Cantón G, Xue Y, Fleg JL, DesvigneNickens P, Klimas MT, Padley RJ, Vassileva MT, Wyman BT, Yuan C. Carotid magnetic resonance imaging for monitoring atherosclerotic plaque progression: a multicenter reproducibility study. Int J Cardiovasc Imaging. 2015;31:95–103.
 6.
Kim SE, Jeong EK, Shi XF, Morrell G, Treiman GS, Parker DL. Diffusionweighted imaging of human carotid artery using 2D singleshot interleaved multislice inner volume diffusionweighted echo planar imaging (2D ssIMIVDWEPI) at 3 T: diffusion measurement in atherosclerotic plaque. J Magn Reson Imaging. 2009;30:1068–77.
 7.
Kim SE, Treiman GS, Roberts JA, Jeong EK, Shi X, Hadley JR, Parker DL. In vivo and ex vivo measurements of the mean ADC values of lipid necrotic core and hemorrhage obtained from diffusion weighted imaging in human atherosclerotic plaques. J Magn Reson Imaging. 2011;34:1167–75.
 8.
Xie Y, Yu W, Fan Z, Nguyen C, Bi X, An J, Zhang T, Zhang Z, Li D. High resolution 3D diffusion cardiovascular magnetic resonance of carotid vessel wall to detect lipid core without contrast media. J Cardiovasc Magn Reson. 2014;16:67.
 9.
Ghazanfari S, DriessenMol A, Strijkers GJ, Kanters FMW, Baaijens FPT, Bouten CVC. A comparative analysis of the collagen architecture in the carotid artery: second harmonic generation versus diffusion tensor imaging. Biochem Biophys Res Commun. 2012;426:54–8.
 10.
Flamini V, Kerskens C, Simms C, Lally C. Fibre orientation of fresh and frozen porcine aorta determined noninvasively using diffusion tensor imaging. Med Eng Phys. 2013;35:765–76.
 11.
Porter DA, Heidemann RM. High resolution diffusionweighted imaging using readoutsegmented echoplanar imaging, parallel imaging and a twodimensional navigatorbased reacquisition. Magn Reson Med. 2009;62:468–75.
 12.
Clarke SE, Hammond RR, Mitchell JR, Rutt BK. Quantitative assessment of carotid plaque composition using multicontrast MRI and registered histology. Magn Reson Med. 2003;50:1199–208.
 13.
Berquist RM, Gledhill KM, Peterson MW, Doan AH, Baxter GT, Yopak KE, Kang N, Walker HJ, Hastings PA, Frank LR. The Digital Fish Library: Using MRI to Digitize, Database, and Document the Morphological Diversity of Fish. PLoS ONE. 2012;7:e34499.
 14.
Klein S, Staring M, Murphy K, Viergever MA, Pluim JPW. elastix: a toolbox for intensitybased medical image registration. IEEE Trans Med Imaging. 2010;29:196–205.
 15.
Shamonin DP, Bron EE, Lelieveldt BPF, Smits M, Klein S, Staring M, Alzheimer’s Disease Neuroimaging Initiative. Fast parallel image registration on CPU and GPU for diagnostic classification of Alzheimer’s disease. Front Neuroinform. 2013;7:50.
 16.
Soares JM, Marques P, Alves V, Sousa N. A hitchhiker’s guide to diffusion tensor imaging. Front Neurosci. 2013;7(March):1–14.
 17.
Klooster R Van, Staring M, Klein S, Kwee RM, Kooi ME, Reiber JHC, Lelieveldt BPF. Automated registration of multispectral MR vessel wall images of the carotid artery Automated registration of multispectral MR vessel wall images. Med Phys. 2013;40(12):121904.
 18.
R Core Team. R: A Language and Environment for Statistical Computing. 2013.
 19.
Kingsley PB. Introduction to diffusion tensor imaging mathematics: Part I. Tensors, rotations, and eigenvectors. Concepts Magn Reson Part A. 2006;28A:101–22.
 20.
Kingsley PB. Introduction to diffusion tensor imaging mathematics: Part II. Anisotropy, diffusionweighting factors, and gradient encoding schemes. Concepts Magn Reson Part A. 2006;28A:123–54.
 21.
Kingsley PB. Introduction to diffusion tensor imaging mathematics: Part III. Tensor calculation, noise, simulations, and optimization. Concepts Magn Reson Part A. 2006;28A:155–79.
 22.
Finlay HM, McCullough L, Canham PB. Threedimensional collagen organization of human brain arteries at different transmural pressures. J Vasc Res. 1995;32:301–12.
 23.
Jackson A. Quantitative MRI of the brain: measuring changes caused by disease. By P Tofts. Br J Radiol. 2005;78:87.
 24.
Miller KL, Pauly JM. Nonlinear phase correction for navigated diffusion imaging. Magn Reson Med. 2003;50:343–53.
 25.
Cohenadad J. HighResolution DWI in Brain and Spinal Cord with syngo RESOLVE 1. Siemens Magnetom  Clin Neurol. 2012;510:16–23.
 26.
Hayes LL, Porter DA, Jones RA, Palasis S, GrattanSmith JD. RESOLVE: A Powerful Tool for Imaging the Pediatric Spine. MAGNETOM Flash. 2014;2:38–45.
 27.
Griswold MA, Jakob PM, Heidemann RM, Nittka M, Jellus V, Wang J, Kiefer B, Haase A. Generalized Autocalibrating Partially Parallel Acquisitions (GRAPPA). Magn Reson Med. 2002;47:1202–10.
Acknowledgments
This work was supported, in part by an EU grant (NanoAthero, grant agreement no 30982).
Funding
Not applicable.
Availability of data and materials
Data on which conclusions rely are presented in the paper.
Authors’ contributions
PO and GR were involved in the design of the 2D diffusion gradient direction sampling scheme, conception, study design, data acquisition and postprocessing, drafting the manuscript and figure preparation. HM, RS and HD were involved in conception and revising the manuscript. All authors read and approved the manuscript.
Competing interests
The authors declare that they have no competing interests.
Consent for publication
All authors read and approved the final manuscript. Written informed consent was obtained from the participant for publication of their individual details and accompanying images in this manuscript. The consent form is held by the authors’ institution and is available for review by the EditorinChief.
Ethics approval and consent to participate
The study was approved by the institutional review board (IRB; Medical University of Graz) as well as followed Declaration of Helsinki recommendations and informed consent was signed by each participant prior to the MR measurements.
Author information
Affiliations
Corresponding author
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Opriessnig, P., Mangge, H., Stollberger, R. et al. In vivo cardiovascular magnetic resonance of 2D vessel wall diffusion anisotropy in carotid arteries. J Cardiovasc Magn Reson 18, 81 (2017). https://doi.org/10.1186/s1296801603048
Received:
Accepted:
Published:
Keywords
 Cardiovascular magnetic resonance
 Diffusion tensor imaging
 Cardiovascular diseases
 Atherosclerosis