Diffusion MR tractography of the heart
Journal of Cardiovascular Magnetic Resonance volume 11, Article number: 47 (2009)
Histological studies have shown that the myocardium consists of an array of crossing helical fiber tracts. Changes in myocardial fiber architecture occur in ischemic heart disease and heart failure, and can be imaged non-destructively with diffusion-encoded MR. Several diffusion-encoding schemes have been developed, ranging from scalar measurements of mean diffusivity to a 6-dimensional imaging technique known as diffusion spectrum imaging or DSI. The properties of DSI make it particularly suited to the generation of 3-dimensional tractograms of myofiber architecture. In this article we review the physical basis of diffusion-tractography in the myocardium and the attributes of the available techniques, placing particular emphasis on DSI. The application of DSI in ischemic heart disease is reviewed, and the requisites for widespread clinical translation of diffusion MR tractography in the heart are discussed.
The myocardium can be studied at several spatial scales. New techniques, such as molecular imaging, are providing important insights into cardiac disease at the cellular and subcellular levels [1–3]. At the other end of the spectrum, parameters of regional and whole organ function such as ejection fraction, perfusion, viability and strain are now routinely used in clinical practice [4, 5]. The microstructural organization of the myocardium, however, has been less extensively studied, although changes at this scale could provide important biological insights and a mechanism linking cellular and whole-organ pathology [6–8]. Here we describe our initial ex-vivo experience with a relatively new magnetic resonance (MR) technique, diffusion spectrum MR tractography, capable of imaging myocardial fiber architecture at the microstructural level.
In a series of breakthrough histological studies, Streeter and colleagues demonstrated that cardiomyocytes form tracts with a crossing helical architecture [9, 10]. Myofiber tracts in the subendocardium have a positive or right-handed helix angle, those in the mid-myocardium are circumferential and those in the subepicardium have a negative or left-handed helix angle [9, 10]. These fiber tracts form laminar sheets [11–14], and it is the shear, extension, thickening and radial reorientation of these sheets that allows the myocardium to thicken in systole [12, 15–17]. Changes in scalar indices of diffusion and myofiber anatomy have been documented in a variety of small and large animal models of cardiac disease [18–23], as well as in humans [24–26]. In the majority of these studies, however, fiber anatomy was visualized only at discrete points in the myocardium. In the current article we focus on the use of diffusion-encoded MR to create continuous 3-dimensional tractograms of myocardial fiber architecture. We place particular emphasis on our recent experience in the heart with diffusion spectrum MR tractography . We review the rationale and theoretical basis of MR tractography, its application in animal models of ischemic heart disease, the properties of other diffusion-encoding schemes such as diffusion tensor and q-ball imaging, and the pathway towards clinical translation of MR tractography in the heart.
Diffusion Spectrum MR
Diffusion imaging can be performed at several levels of complexity, ranging from the simple acquisition of a single diffusion-weighted image to the complex but robust acquisition scheme used in diffusion spectrum imaging (DSI) [28, 29]. Diffusion tensor imaging (DTI) and q-ball imaging can be thought of as formalisms that sample diffusion or q-space with an intermediate level of complexity [28, 29]. While complex, DSI is the only technique derived directly from first principles , is hypothesis free [28, 30], broadly generalizable, and is regarded by many (including the authors) as the gold standard diffusion imaging technique [28, 31–33]. In the current implementation of DSI, q-space is sampled with 515 diffusion-encoding vectors or q-vectors, although simulations containing up to 925 q-vectors have been performed . The angular resolution and accuracy of DSI result in large part from the number and distribution of the samples acquired in q-space, analogous to the manner in which the region of support and sample density of k-space influence the spatial resolution and field-of-view of an image.
The q-vectors used to sample q-space in a DSI experiment vary in both their strength (b-value) and spatial orientation, sampling q-space in a dense 3D lattice. The b-value of a q-vector is proportional to the product of the square of the gradient strength and the diffusion time interval. (b ~ q2Δ, where q = γδG and γ is the gyromagnetic ratio, δ is the gradient duration, G is the gradient strength and Δ is the diffusion time interval). The b-value in a diffusion-encoded acquisition is in some ways analogous to the degree of velocity encoding (venc) used during phase contrast imaging. Higher b-values increase the resolution of the diffusion spectrum, and are thus desirable. However, an excessively high b-value can reduce image signal-to-noise ratio (SNR) severely, and an optimal balance between these two competing factors must thus be struck [31, 33]. B-values greater than 10,000 s/mm2 have been used both in-vivo and ex-vivo in the brain and myocardium [27, 30, 33].
In a DSI acquisition, each voxel in the spatial (x, y, z) domain of an image has its own 3D q-space associated with it. DSI images are thus 6-dimensional in which the 3 dimensions of q-space are superimposed on the 3 dimensions of image space [28, 30]. Q-space is sampled in a DSI experiment by applying a B0 gradient field (q-vector) along a specific spatial orientation. The duration and intensity of the applied field (b-value) determines the length of the q-vector relative to the origin of q-space, while the vector direction is determined by the spatial orientation of the applied gradients. The resulting signal intensity determines the coefficient for that value of Q in each voxel. A 3D image with 96 × 96 × 96 voxels will thus have greater than 8 × 105 q-space datasets, with each dataset containing 515 coefficients.
The q-space formalism, for a single voxel and its associated q-space dataset, is depicted graphically in figure 1. The value of each coefficient in the q-space dataset of the voxel is determined by the signal intensity in the voxel during the application of the q-vector defined by that q-space coefficient. The impact of diffusion-encoding on the intensity and phase of the MR signal will be easily understood by those familiar with phase-contrast imaging in cardiovascular magnetic resonance (CMR). Velocity encoding during phase contrast CMR involves the application of two gradient lobes of equal magnitude but opposite polarity, producing a phase differential in the presence of flow. Likewise, as described by Stejskal and Tanner , if a pair of diffusion-encoding gradients with equal magnitude but opposite polarity is applied in the direction of diffusion, the moving spin undergoes a phase shift resulting in incomplete rephasing and attenuation of the MR signal [28, 29]. Water diffuses primarily along the long axis of myofibers , as shown in the schematic in figure 1. The intensity of a q-space coefficient will thus be greatest when a small q-vector is directed orthogonal to the myofibers in the voxel, and lowest when a large q-vector is applied in the direction of myofiber orientation (figure 1) [28, 35].
Q-space and k-space are both discretely sampled Cartesian 3-spaces that are Fourier transformed to yield a useful result. While inverse Fourier transformation of k-space produces an anatomical image, the inverse Fourier transform of q-space produces a probability density function (PDF) of water diffusion (and indirectly fiber orientation) in the voxel (figure 1) [28, 30]. The PDF describes the probability that a water molecule will diffuse a certain distance in a particular direction during the diffusion MR acquisition. Since only the directional information is required for DSI tractography, the PDF is usually reduced by radial integration to an orientation distribution function (ODF) [28, 30]. The PDF and ODF in each voxel can have multiple local maxima, each resolving an individual fiber population in the voxel. DSI is thus able to resolve multiple fiber populations in a voxel with both high angular and spatial resolution, including crossing and converging fibers [30, 33, 36]. The reader interested in a more mathematical description of DSI tractography is referred to figure 2 of this article and to prior work in the field .
The merits and limitations of other diffusion-encoding schemes (q-ball imaging, spherical deconvolution, diffusion tensor imaging) will be discussed in detail later in the article. These techniques sample q-space less fully than DSI, enforce a diffusivity model (Funk transform, spherical harmonics, tensor) on the sampled data and degeneralize DSI to some degree. The simplest of these models that still includes directional information is DTI, but this differs however from DSI in several important ways: With DSI each distinct fiber population in a voxel is represented by a unique local maximum in the PDF/ODF [30, 37, 38]. As described below, however, the principal eigenvector of a diffusion tensor reflects the average direction of myofiber orientation in the voxel . In addition, while the spatial resolution of DTI is determined by/identical to the resolution of the image, the 6-dimensional nature of DSI provides subvoxel resolution, determined by the resolution of the ODF . The attributes of DSI are thus inherently suited to the generation of 3D tractograms of myocardial fiber architecture.
Tractograms can be generated using deterministic or probabilistic algorithms [29, 39]. Probabilistic models address the uncertainty associated with the pathways of the reconstructed tractograms, which can be significant when DTI is used [29, 39]. The properties of DSI, however, provide a robust platform for the creation of deterministic myofiber tractograms since the high angular and spatial resolution of DSI significantly reduces uncertainty. Local maxima in the ODFs of each voxel form a vector field that can be integrated into streamlines or tractograms, representing pathways of maximum diffusion coherence [30, 32, 36–38]. The streamline will follow the path of minimal angular difference between adjacent ODFs (voxels), with a threshold angle of 18-35° used to halt further propagation of the streamline. Fiber tracts in the myocardium can be depicted in terms of their direction along Cartesian axes , but should optimally be depicted in terms of the helix or spiral angle they make with the long axis of the left ventricle, in accordance with the histological pattern described by Streeter [27, 41, 42]. DSI tractography of a normal rat heart, color coded by helix angle, is shown in figure 3.
DSI Tractography of the Myocardium
The largest experience with DSI tractography in the myocardium to date has been in the imaging of normal and infarcted rat hearts ex-vivo . Preliminary experience with the technique in excised hearts from mice, lambs and sheep, however, has been highly encouraging . DSI tractography was able to robustly resolve the anisotropy of myocardial fiber architecture in excised fixed rat hearts . As shown in figures 3, 4, 5 and 6 of normal rat hearts, the arrangement of fibers in the myocardium into an array of crossing helical structures was resolved in exquisite detail. Myofibers in the subendocardium can be seen to form a positive or right-handed helix, those in the subepicardium to form a negative or left-handed helix, and myofibers in the mid-myocardium to be circumferential (zero helix angle). The 6-dimensional nature of DSI and its subvoxel resolution resulted in dense tractographic datasets without gaps or discontinuities between adjacent myofibers . The high angular and spatial resolution of DSI tractography can be fully appreciated in the magnified images shown in figures 4 and 5.
DSI tractograms can be visualized as projection images or as tomographic reconstructions of the 3D dataset. Tomographic representation involves the selection of a plane in the 3D field, the thickness of which is defined by the user. Only those myofibers that intersect the plane are shown in the reconstructed image (figures 4 and 5). The density of fibers in the myocardium, however, can make projection images and even tomographic reconstructions complex to interpret. Spherical or discoid regions-of-interest (ROIs) can thus be defined to visualize only those fiber tracts intersecting the ROI (figures 3 and 6). As shown in figure 6, myofiber tracts in the subendocardium and subepicardium form half-turns of a spiral but have orthogonal helix angles. Fibers in the subendocardium of the lateral wall track from the posterior-base to the anterior apex, while those in the subepicardium track from the anterior-base to the posterior apex .
Fiber tracts in a given transmural plane in the septum have the identical helix angle to those in the lateral wall but an opposite alignment (subendocardium: anterior-base to the posterior apex and subepicardium: posterior-base to the anterior apex) and thus complete a turn of their respective helices. Throughout the myocardium a smooth evolution in myofiber helix angle was seen, as shown in figures 3 and 6. Adjacent myofiber tracts with similar helix angles could be seen to form a sheet-like structure both with DSI and histologically (figures 4 and 7) . Little dispersion in myofiber helix angle was seen in a given transmural plane and in normal myocardium fibers with orthogonal helix angles made no contact and were separated from each other by myofibers with intermediate helix angles (figures 6 and 7) .
The impact of ischemic injury on myofiber architecture in infarcted rat hearts was also examined in this study . The infarcted hearts were perfused-fixed with paraformaldehyde three weeks after permanent left coronary artery ligation, and imaged at 4.7 Tesla with a maximum b-value of approximately 10,000. Myocardial fiber architecture was extremely perturbed in the infarcted hearts (figure 8). However, in all cases a large number of residual myofibers were seen within the infarcts, particularly in the most basal portions of the infarct . We hypothesize that these residual myofibers persist within the infarct due to pre-existing collateral networks, but this will require further study in a range of animal models.
The residual myofibers within the infarcts showed several interesting microstructural properties. Many of the residual myofibers had helix angles consistent with subendocardial fibers (figures 8 and 9). Moreover, the residual subendocardial myofibers frequently made contact with residual mid and subepicardial myofibers in the infarct (figures 8 and 9) . A wide dispersion in helix angles was seen in the residual fibers within the infarct. Moreover, presumably due to thinning and expansion of the infarct, residual myofibers with orthogonal helix angles frequently came into contact with one another (figures 8 and 9) . Nodes of orthogonal myofiber intersection were also seen in the infarct, presumably in areas with severe local deformation in myofiber anatomy (figure 9). The residual myofibers within the infarct thus frequently formed a mesh-like structure, producing nodes of orthogonal myofiber intersection or contact (NOMIC) within the infarct . The presence of similar networks of residual myofibers will need to be confirmed in other animal species as well as in humans. Further study will also be needed to determine the mechanical and electrophysiological implications of these networks of residual myofibers and the nodes of orthogonal myofiber contact within them.
Diffusion Tensor MR (DTI)
DTI is a formalism that samples q-space more rapidly than DSI but is based on the assumption that the diffusion in a voxel is Gaussian [28, 29]. This assumption holds true when a voxel contains only a single myofiber population but introduces bias and uncertainty when this is not the case. Several useful scalar parameters of diffusion in the myocardium can be derived from DTI [18–23, 28, 29], but the technique is less suited to the generation of tractograms than DSI [28, 33, 38], as discussed below. The diffusion tensor is a symmetric second order (3 × 3) tensor (figure 10). Values along the diagonal (D xx , D yy and D zz ) describe the degree of diffusion along the principal laboratory axes (x, y, z). Off-diagonal elements of the tensor describe the degree of correlation between diffusion in two directions. The nature of diffusion results in the tensor being symmetric and the transpose of the tensor is thus identical to the original tensor (alternatively stated the off-diagonal elements above and below the diagonal are equal to each other). The 3 × 3 diffusion tensor thus contains 6 independent values, and requires diffusion encoding to be performed in a minimum of 6 independent (non-linear) directions. DTI can thus be viewed conceptually as a linear approximation of DSI, in which q-space is sampled with 6 rather than 515 q-vectors .
Diagonalization of the diffusion tensor rotates its axes out of the laboratory frame (x, y, z) and along the eigenvectors of the tensor, which reflect the diffusion properties of the tissue in the voxel (figure 10) [28, 29]. The principal or largest eigenvector, designated λ1, describes the direction of diffusion along the long axis of the myofibers in the voxel, the second largest eigenvector the direction of the myofiber sheets and the third eigenvector the direction normal to the myofibers [14, 16, 35]. Several important scalar parameters, such as mean diffusivity and fractional anisotropy, can be derived from the eigenvalues associated with these eigenvectors [28, 29]. Studies in animals and humans have shown that both mean diffusivity and fractional anisotropy change significantly in infarcted and healing myocardium [18–26].
The principal eigenvector of the tensor can be used to estimate the average direction and helix angle of the myofibers in a voxel [14, 35]. A reduction in right-handed (subendocardial) fibers and an increase in left-handed (subepicardial) fibers has been noted in animals and patients with myocardial infarction [22, 25]. While of substantial value, several significant limitations of DTI impact cardiac tractography and merit discussion. DTI tractography can only resolve one fiber population, described by the principal eigenvector, in a given voxel. Moreover, the principal eigenvector is a composite measure of mean diffusion in the voxel. If a voxel contains more than one fiber population the principal eigenvector will represent the average direction of diffusion in the voxel, which may actually not be an accurate representation of the fibers in the voxel [28, 39, 44]. In the context of tractography, DTI thus reduces the information contained in the PDF into a single average value, the principal eigenvector. DTI can thus be viewed conceptually as a linear approximation of the displacement spectrum.
The principal eigenvectors in adjacent voxels can be connected using several algorithms to form streamlines of fiber tracts (figure 10) [39, 44]. DTI tractograms, however, are limited in both angular and spatial resolution and are susceptible to a degree of bias and uncertainty introduced by the undersampling of q-space [39, 44]. Because only one average fiber population per voxel can be resolved with DTI, it is unable to detect complex and converging fiber anatomy [33, 38]. The lower number of seed points produced with DTI also leads to smaller and less complete tractography datasets than DSI. In theory, depending on the number of fiber populations per voxel, dramatic improvements in the spatial resolution of DTI could produce a tractographic dataset similar to DSI. In practice, however, this cannot be done because DTI acquisitions are already highly SNR constrained due to the phase dispersion induced by the diffusion-encoding gradients. A simple isotropic doubling of the resolution of a DTI acquisition, for instance, would require 64 signal averages to maintain SNR. DTI acquisitions with near microscopic resolution would thus require prohibitively long scan times even ex-vivo to maintain adequate SNR. The longer readout duration of these ultra-high resolution scans would also increase the TE and thus reduce SNR even further. The parameters of a diffusion-encoded acquisition, particularly in-vivo, are thus frequently dominated by the need to achieve adequate SNR, as discussed below.
In-Vivo Diffusion Imaging
Diffusion, strain and velocity encoded MR are all displacement encoding techniques that rely on the presence of a residual phase following the application of equal and opposite gradients. The spatial scale of the displacement due to diffusion, however, is far lower and requires significant modifications to be made in the acquisition scheme and/or the strength of the applied gradients. Diffusion encoded acquisitions are usually performed with single-shot readouts such as EPI (echoplanar MRI) or HASTE (half Fourier acquired single shot turbo spin echo) , ensuring all lines of k-space have the same phase. Multishot diffusion imaging has been performed but requires a scheme to detect and correct random shot-to-shot phase changes produced by motion and the diffusion-encoding gradients . The incorporation of diffusion encoding in to a HASTE sequence violates the CPMG condition since δ, and hence the time between the 90° and the first 180° refocusing pulse, is significantly greater than half the time between successive 180° refocusing pulses in the readout. This causes phase cancellation between echoes and poor image quality. Experimental techniques have been developed to address phase incoherence in diffusion-encoded HASTE, but many of these eliminate 50% of the signal (either even or odd echoes) and thus suffer from low SNR . Despite the potential for susceptibility, chemical shift and ghosting artifacts to occur, the vast majority of diffusion-encoded acquisitions are thus performed with single-shot EPI.
The potential resolution of a single-shot EPI readout in a stationary tissue such as the brain is determined by the T2 of the tissue (which limits the TE that can be used) and the overall time available for image acquisition. EPI of the heart, however, must also contend with susceptibility artifacts from the lungs, the high fat content in the chest wall and a finite acquisition window in which image quality is not degraded by motion. Despite the use of parallel imaging and multi-element arrays, the spatial resolution of single shot EPI in the heart is limited. This places a fundamental limit on the resolution of DTI, which has the identical resolution to the EPI image. Diffusion-encoding is extremely sensitive to bulk motion , which in the heart is significantly greater than the motion of water due to diffusion. Diffusion-encoded sequences in the heart must thus be designed to be insensitive to bulk cardiac motion as well as myocardial strain in order to be performed in-vivo [47–49]. Several approaches have been described to accomplish this (figure 11), depending in large part on the gradient strength of the clinical scanner.
Tractography of stationary structures such as the brain and the myocardium ex-vivo can be performed with a Stejskal-Tanner diffusion encoded sequence [28, 29, 34]. In practice, however, diffusion-encoded imaging of the brain is frequently performed with a twice refocused spin echo sequence to limit the effects of eddy currents on image quality . The diffusion-encoding gradients in the Stejskal-Tanner sequence are placed on either side of a 180° refocusing pulse [28, 29, 34], which refocuses the diffusion-encoded signal at the echo time of the spin echo EPI readout (figure 11A). The use of this sequence in the heart in-vivo, however, is precluded by motion. A stimulated echo approach has thus been developed to overcome this (figure 11B) : Three 90° radiofrequency (RF) excitation pulses are applied within two successive RR intervals. The first and third excitation pulses have identical trigger delays from the onset of the successive R-waves and are both followed by unipolar diffusion encoding gradients. The second 90 degree RF pulse is placed a duration of TE/2 from the first, and effectively flips the transverse magnetization back into the longitudinal plane, where it is subject to R1 decay but not to R2* decay and motion induced phase change. The diffusion sensitivity of the sequence is determined by the physical displacement of water between the onset of the first and second unipolar diffusion encoding gradients. Diffusion sensitization thus continues to occur while the magnetization produced by the first excitation pulse is stored and protected in the longitudinal axis during a period known as the mixing time (TM). The application of the third 90-degree RF pulse and the second unipolar gradient produce a diffusion-encoded stimulated echo that has adequate diffusion sensitization and is largely free of motion related artifacts .
Several caveats of this approach, however, need to be considered. Diffusion imaging is frequently signal-to-noise (SNR) constrained and a stimulated echo EPI sequence has half the SNR of a spin-echo EPI readout. In addition, the delay between the R wave and the first and third RF pulses needs to be carefully selected to eliminate the effects of myocardial strain on the diffusion measurements . Diffusion measurements (particularly the second and third eigenvectors) with this technique may thus be optimally made only at certain "sweet spots" in the cardiac cycle, for instance at midsystole, where the strain effects become negligible . To overcome this limitation myocardial strain data can be acquired with the diffusion MR data and be used to retrospectively correct the diffusion measurements [48, 49]. A strain-insensitive stimulated echo sequence has also been developed and uses a pair of bipolar diffusion-encoding gradients, rather than unipolar diffusion-encoding gradients (figure 11C) . Diffusion encoding with this sequence, however, occurs only during the duration of the 2 pairs of diffusion encoding gradients, which are insensitive to first order motion terms (flow compensated) . While this sequence is strain-insensitive and can be performed at any stage of the cardiac cycle, the gradient duration required for adequate diffusion encoding (b-values) significantly lengthens the TE and reduces SNR.
The stimulated echo approaches described above have been successfully used to perform DTI in patients in-vivo [24–26], but on a limited scale in a few centers of expertise. More widespread clinical performance of DSI tractography in the heart will require several technical and scientific advances to be made. Gradient technology on whole body MR scanners will need to be improved by at least a factor of 2. The potential of improved gradient performance was demonstrated recently by Gamper and colleagues on an 3 Tesla clinical system equipped with an 87 mT/m gradient . The strength of this gradient allowed a spin-echo EPI readout to be used, avoiding the SNR penalty of a stimulated echo approach. A pair of bipolar diffusion encoding gradients was placed on either side of a 180° refocusing pulse in a flow-compensated modification of the Stejskal-Tanner approach (figure 11D) . While the use of this approach does impose certain constraints, it overcomes several factors impeding the performance of diffusion tractography in-vivo. Clinical translation will also require improved techniques for whole-heart imaging to be developed. The development of multi-element arrays, including a recently developed 128-element cardiac array , has the potential to facilitate the acquisition of volumetric whole-heart datasets in a single breathold. In addition, improvements in radiofrequency and navigator technology will also facilitate the acquisition of volumetric diffusion encoded data of the heart in-vivo .
Q-Space Sampling In-Vivo
Several q-space sampling schemes have been proposed for tractography of the brain in-vivo [30, 31, 55, 56]. The optimal q-space sampling scheme for in-vivo tractography of the heart, however, needs to be considered in the context of the SNR and motion-imposed constraints specific to the heart. Nevertheless, the experience in the brain with different q-space sampling schemes remains highly relevant to in-vivo tractography of the myocardium and worthy of discussion. Diffusion tractography in the brain is SNR constrained [39, 44]. Thus, even if only 6 independent diffusion encoding vectors are applied in order to perform DTI tractography, several signal averages need to be performed to achieve adequate SNR. This, however, is being less frequently done . Rather than averaging data produced by the same diffusion-encoding vector several times, the scan time is used to acquire data in more than 6 directions [44, 57–59]. The improvement in SNR, which is dependent on the total number of acquisitions, is similar with the two approaches. However, the accuracy of the data derived from the tensor is improved when a greater number of directions are sampled [57–59].
A voxel containing n individual myofiber populations will have 3n degrees of freedom, and likely require up to 6n independent q-vectors to fully resolve diffusion in the voxel. Whether the application of greater than 6n q-vectors would be desirable would depend in large part on the SNR of the data. The application of additional q-vectors would consume time but on the other hand increase the SNR of the image and also potentially de-alias the PDF/ODF. The more densely q-space is sampled the larger the PDF/ODF becomes, reducing the potential for aliasing. (This is analogous to increasing the FOV of an image to reduce aliasing in the spatial domain.) The approach favored by many, including ourselves, is thus to sample q-space as densely as possible within the limitations imposed by acquisition time, patient tolerance and the pulse sequence considerations discussed above.
Several high angular resolution diffusion imaging (HARDI) techniques have been developed that sample q-space more densely than DTI, but less so than DSI [55, 56, 60]. All of these techniques are based on certain hypotheses and assumptions, but are easier to implement in-vivo than DSI and may be well suited to imaging of the myocardium under certain scenarios. Q-ball imaging for instance involves the use of q-vectors that all have an identical and fairly large b-value [28, 56, 60]. Rather than sampling a 3D lattice in q-space, the technique samples the surface of a sphere with a given radius in q-space. The technique is simpler and more rapid than DSI but assumes that the selected b-value is optimal for the detection of all fiber or nerve tracts in the tissue, which may not always be the case. In normal myocardium the myofiber tracts have reasonably similar lengths and morphology (figure 6) and q-ball imaging may perform very well in this scenario. Infarcted myocardium, however, much like the brain has a highly heterogeneous population of myofibers, which may not all be detected optimally at the selected b-value. (A useful analogy to consider is the use of a single preset velocity encoding gradient to image several hemodynamic jets, despite large potential variations in the velocities of these jets). Nevertheless, q-ball imaging maintains many of the attributes of DSI, is easier to implement clinically and has the potential to be of significant value in the myocardium.
The accuracy of q-ball imaging and other heuristic methods will need to be validated against DSI tractography, which should be considered the reference gold-standard approach. Several approaches have now been used to validate the accuracy of DSI tractography ex-vivo (Table 1) and provide a solid basis for the translation of the technique. DSI tractography of the heart, tongue and brain ex-vivo has consistently resolved the fiber/tract patterns known to occur in normal organs [27, 33, 36]. In infarcted rat hearts a high degree of correlation has been seen between fiber architecture by DSI tractography and by histology . In the brain autoradiography and manganese enhanced MR have confirmed the accuracy of DSI tractography [37, 61], and the technique has been correlated with two-photon microscopy in the tongue . A technique to validate DSI tractography in-vivo by detecting excess noise and low confidence in the dataset has also recently been developed . The technique involves the random reshuffling of the voxel ODFs and the generation of tractograms from the scrambled vector field. The comparison of parameters in the tractographic datasets generated from the original and the reshuffled ODF vector fields provides an index of image noise and fidelity . Edge weight (link between two nodes in a neuro-imaging dataset) has been used as a parameter of image confidence in the brain, and analogous parameters will need to be developed to assess the impact of noise on the accuracy of in-vivo tractography datasets in the heart. As with all new technologies, however, the ultimate validation of diffusion-tractography in the heart will be determined by its ability to influence hard clinical endpoints in cardiovascular medicine.
The potential of diffusion-tractography, and in particular DSI tractography, to resolve microstructural fiber anatomy in the heart has been demonstrated in several studies ex-vivo. At present, hardware limitations on most clinical scanners constitute the principal impediment to the clinical translation of diffusion tractography in the heart. Progress in the field will thus be rapidly accelerated by the widespread introduction of clinical scanners with gradients greater than 80 mT/m. Advances in RF technology, multi-element arrays, navigators and parallel acquisition schemes will also facilitate the clinical translation of more advanced q-space acquisition schemes. Clinical translation of this promising technology, however, will be a major challenge, requiring an excellent level of collaboration between engineers, industry and physicians. Diffusion tractography does not involve ionizing radiation or exogenous contrast media, and poses no risk to the patient. The technique images the myocardium at the microstructural scale and provides information that is highly complementary to that provided by other modalities and imaging techniques. Diffusion tractography of the myocardium expands the breadth and scope of cardiovascular magnetic resonance and has the potential to become an extremely powerful tool in both the research and clinical settings.
Sosnovik DE, Nahrendorf M, Weissleder R: Molecular magnetic resonance imaging in cardiovascular medicine. Circulation. 2007, 115: 2076-2086. 10.1161/CIRCULATIONAHA.106.658930.
Nahrendorf M, Sosnovik DE, French BA, Swirski FK, Bengel F, Sadeghi MM, Lindner JR, Wu JC, Kraitchman DL, Fayad ZA, Sinusas AJ: Multimodality Cardiovascular Molecular Imaging, Part II. Circ Cardiovasc Imaging. 2009, 2: 56-70. 10.1161/CIRCIMAGING.108.839092.
Sosnovik DE, Nahrendorf M, Weissleder R: Targeted imaging of myocardial damage. Nat Clin Pract Cardiovasc Med. 2008, 5 (Suppl 2): S63-70. 10.1038/ncpcardio1115.
Assomull RG, Pennell DJ, Prasad SK: Cardiovascular magnetic resonance in the evaluation of heart failure. Heart. 2007, 93: 985-992. 10.1136/hrt.2003.025304.
Lockie T, Nagel E, Redwood S, Plein S: Use of cardiovascular magnetic resonance imaging in acute coronary syndromes. Circulation. 2009, 119: 1671-1681. 10.1161/CIRCULATIONAHA.108.816512.
Buckberg GD, Weisfeldt ML, Ballester M, Beyar R, Burkhoff D, Coghlan HC, Doyle M, Epstein ND, Gharib M, Ideker RE, Ingels NB, LeWinter MM, McCulloch AD, Pohost GM, Reinlib LJ, Sahn DJ, Sopko G, Spinale FG, Spotnitz HM, Torrent-Guasp F, Shapiro EP: Left ventricular form and function: scientific priorities and strategic planning for development of new views of disease. Circulation. 2004, 110: e333-336. 10.1161/01.CIR.0000143625.56882.5C.
Buckberg G, Hoffman JI, Mahajan A, Saleh S, Coghlan C: Cardiac mechanics revisited: the relationship of cardiac architecture to ventricular function. Circulation. 2008, 118: 2571-2587. 10.1161/CIRCULATIONAHA.107.754424.
Kramer CM: Insights into Myocardial Microstructure During Infarct Healing and Remodeling: Pathologists Need Not Apply. Circ Cardiovasc Imaging. 2009, 2: 4-5. 10.1161/CIRCIMAGING.108.841619.
Streeter DD, Spotnitz HM, Patel DP, Ross J, Sonnenblick EH: Fiber orientation in the canine left ventricle during diastole and systole. Circ Res. 1969, 24: 339-347.
Streeter DD, Hanna WT: Engineering mechanics for successive states in canine left ventricular myocardium. II. Fiber angle and sarcomere length. Circ Res. 1973, 33: 656-664.
LeGrice IJ, Smaill BH, Chai LZ, Edgar SG, Gavin JB, Hunter PJ: Laminar structure of the heart: ventricular myocyte arrangement and connective tissue architecture in the dog. Am J Physiol. 1995, 269: H571-582.
Dou J, Reese TG, Tseng WY, Wedeen VJ: Cardiac diffusion MRI without motion effects. Magn Reson Med. 2002, 48: 105-114. 10.1002/mrm.10188.
Helm PA, Tseng HJ, Younes L, McVeigh ER, Winslow RL: Ex vivo 3D diffusion tensor imaging and quantification of cardiac laminar structure. Magn Reson Med. 2005, 54: 850-859. 10.1002/mrm.20622.
Tseng WY, Wedeen VJ, Reese TG, Smith RN, Halpern EF: Diffusion tensor MRI of myocardial fibers and sheets: correspondence with visible cut-face texture. J Magn Reson Imaging. 2003, 17: 31-42. 10.1002/jmri.10223.
LeGrice IJ, Takayama Y, Covell JW: Transverse shear along myocardial cleavage planes provides a mechanism for normal systolic wall thickening. Circ Res. 1995, 77: 182-193.
Dou J, Tseng WY, Reese TG, Wedeen VJ: Combined diffusion and strain MRI reveals structure and function of human myocardial laminar sheets in vivo. Magn Reson Med. 2003, 50: 107-113. 10.1002/mrm.10482.
Chen J, Liu W, Zhang H, Lacy L, Yang X, Song SK, Wickline SA, Yu X: Regional ventricular wall thickening reflects changes in cardiac fiber and sheet structure during contraction: quantification with diffusion tensor MRI. Am J Physiol Heart Circ Physiol. 2005, 289: H1898-1907. 10.1152/ajpheart.00041.2005.
Strijkers GJ, Bouts A, Blankesteijn WM, Peeters TH, Vilanova A, van Prooijen MC, Sanders HM, Heijman E, Nicolay K: Diffusion tensor imaging of left ventricular remodeling in response to myocardial infarction in the mouse. NMR Biomed. 2009, 22: 182-190. 10.1002/nbm.1299.
Chen J, Song SK, Liu W, McLean M, Allen JS, Tan J, Wickline SA, Yu X: Remodeling of cardiac fiber structure after infarction in rats quantified with diffusion tensor MRI. Am J Physiol Heart Circ Physiol. 2003, 285: H946-954.
Ripplinger CM, Li W, Hadley J, Chen J, Rothenberg F, Lombardi R, Wickline SA, Marian AJ, Efimov IR: Enhanced transmural fiber rotation and connexin 43 heterogeneity are associated with an increased upper limit of vulnerability in a transgenic rabbit model of human hypertrophic cardiomyopathy. Circ Res. 2007, 101: 1049-1057. 10.1161/CIRCRESAHA.107.161240.
Helm PA, Younes L, Beg MF, Ennis DB, Leclercq C, Faris OP, McVeigh E, Kass D, Miller MI, Winslow RL: Evidence of structural remodeling in the dyssynchronous failing heart. Circ Res. 2006, 98: 125-132. 10.1161/01.RES.0000199396.30688.eb.
Wu EX, Wu Y, Nicholls JM, Wang J, Liao S, Zhu S, Lau CP, Tse HF: MR diffusion tensor imaging study of postinfarct myocardium structural remodeling in a porcine model. Magn Reson Med. 2007, 58: 687-695. 10.1002/mrm.21350.
Wu Y, Chan CW, Nicholls JM, Liao S, Tse HF, Wu EX: MR study of the effect of infarct size and location on left ventricular functional and microstructural alterations in porcine models. J Magn Reson Imaging. 2009, 29: 305-312. 10.1002/jmri.21598.
Tseng WY, Dou J, Reese TG, Wedeen VJ: Imaging myocardial fiber disarray and intramural strain hypokinesis in hypertrophic cardiomyopathy with MRI. J Magn Reson Imaging. 2006, 23: 1-8. 10.1002/jmri.20473.
Wu MT, Tseng WY, Su MY, Liu CP, Chiou KR, Wedeen VJ, Reese TG, Yang CF: Diffusion tensor magnetic resonance imaging mapping the fiber architecture remodeling in human myocardium after infarction: correlation with viability and wall motion. Circulation. 2006, 114: 1036-1045. 10.1161/CIRCULATIONAHA.105.545863.
Wu MT, Su MY, Huang YL, Chiou KR, Yang P, Pan HB, Reese TG, Wedeen VJ, Tseng WY: Sequential Changes of Myocardial Microstructure in Patients Postmyocardial Infarction by Diffusion-Tensor Cardiac MR: Correlation With Left Ventricular Structure and Function. Circ Cardiovasc Imaging. 2009, 2: 32-40. 10.1161/CIRCIMAGING.108.778902.
Sosnovik DE, Wang R, Dai G, Wang T, Aikawa E, Novikov M, Rosenzweig A, Gilbert RJ, Wedeen VJ: Diffusion Spectrum MRI Tractography Reveals the Presence of a Complex Network of Residual Myofibers in Infarcted Myocardium. Circ Cardiovasc Imaging. 2009, 2: 206-212. 10.1161/CIRCIMAGING.108.815050.
Hagmann P, Jonasson L, Maeder P, Thiran JP, Wedeen VJ, Meuli R: Understanding diffusion MR imaging techniques: from scalar diffusion-weighted imaging to diffusion tensor imaging and beyond. Radiographics. 2006, 26 (Suppl 1): S205-223. 10.1148/rg.26si065510.
Mukherjee P, Berman JI, Chung SW, Hess CP, Henry RG: Diffusion tensor MR imaging and fiber tractography: theoretic underpinnings. AJNR Am J Neuroradiol. 2008, 29: 632-641. 10.3174/ajnr.A1051.
Wedeen VJ, Hagmann P, Tseng WY, Reese TG, Weisskoff RM: Mapping complex tissue architecture with diffusion spectrum magnetic resonance imaging. Magn Reson Med. 2005, 54: 1377-1386. 10.1002/mrm.20642.
Kuo LW, Chen JH, Wedeen VJ, Tseng WY: Optimization of diffusion spectrum imaging and q-ball imaging on clinical MRI system. Neuroimage. 2008, 41: 7-18. 10.1016/j.neuroimage.2008.02.016.
Hagmann P, Kurant M, Gigandet X, Thiran P, Wedeen VJ, Meuli R, Thiran JP: Mapping human whole-brain structural networks with diffusion MRI. PLoS ONE. 2007, 2: e597-10.1371/journal.pone.0000597.
Wedeen VJ, Wang RP, Schmahmann JD, Benner T, Tseng WY, Dai G, Pandya DN, Hagmann P, D'Arceuil H, de Crespigny AJ: Diffusion spectrum magnetic resonance imaging (DSI) tractography of crossing fibers. Neuroimage. 2008, 41: 1267-1277. 10.1016/j.neuroimage.2008.03.036.
Stejskal EO, Tanner JE: Spin diffusion measurements: Spin echoes in the presence of a time- dependent field gradient. The Journal of Chemical Physics. 1965, 42: 288-292. 10.1063/1.1695690.
Garrido L, Wedeen VJ, Kwong KK, Spencer UM, Kantor HL: Anisotropy of water diffusion in the myocardium of the rat. Circ Res. 1994, 74: 789-793.
Gilbert RJ, Wedeen VJ, Magnusson LH, Benner T, Wang R, Dai G, Napadow VJ, Roche KK: Three-dimensional myoarchitecture of the bovine tongue demonstrated by diffusion spectrum magnetic resonance imaging with tractography. Anat Rec A Discov Mol Cell Evol Biol. 2006, 288: 1173-1182.
Lin CP, Wedeen VJ, Chen JH, Yao C, Tseng WY: Validation of diffusion spectrum magnetic resonance imaging with manganese-enhanced rat optic tracts and ex vivo phantoms. Neuroimage. 2003, 19: 482-495. 10.1016/S1053-8119(03)00154-X.
Gilbert RJ, Magnusson LH, Napadow VJ, Benner T, Wang R, Wedeen VJ: Mapping complex myoarchitecture in the bovine tongue with diffusion-spectrum magnetic resonance imaging. Biophys J. 2006, 91: 1014-1022. 10.1529/biophysj.105.068015.
Mori S, van Zijl PC: Fiber tracking: principles and strategies - a technical review. NMR Biomed. 2002, 15: 468-480. 10.1002/nbm.781.
Schmid P, Jaermann T, Boesiger P, Niederer PF, Lunkenheimer PP, Cryer CW, Anderson RH: Ventricular myocardial architecture as visualised in postmortem swine hearts using magnetic resonance diffusion tensor imaging. Eur J Cardiothorac Surg. 2005, 27: 468-472. 10.1016/j.ejcts.2004.11.036.
Helm P, Beg MF, Miller MI, Winslow RL: Measuring and mapping cardiac fiber and laminar architecture using diffusion tensor MR imaging. Ann N Y Acad Sci. 2005, 1047: 296-307. 10.1196/annals.1341.026.
Rohmer D, Sitek A, Gullberg GT: Reconstruction and visualization of fiber and laminar structure in the normal human heart from ex vivo diffusion tensor magnetic resonance imaging (DTMRI) data. Invest Radiol. 2007, 42: 777-789. 10.1097/RLI.0b013e3181238330.
Hung J, Solis J, Guerrero JL, Braithwaite GJ, Muratoglu OK, Chaput M, Fernandez-Friera L, Handschumacher MD, Wedeen VJ, Houser S, Vlahakes GJ, Levine RA: A novel approach for reducing ischemic mitral regurgitation by injection of a polymer to reverse remodel and reposition displaced papillary muscles. Circulation. 2008, 118: S263-269. 10.1161/CIRCULATIONAHA.107.756502.
Mukherjee P, Chung SW, Berman JI, Hess CP, Henry RG: Diffusion tensor MR imaging and fiber tractography: technical considerations. AJNR Am J Neuroradiol. 2008, 29: 843-852. 10.3174/ajnr.A1052.
de Crespigny AJ, Marks MP, Enzmann DR, Moseley ME: Navigated diffusion imaging of normal and ischemic human brain. Magn Reson Med. 1995, 33: 720-728. 10.1002/mrm.1910330518.
Alsop DC: Phase insensitive preparation of single-shot RARE: application to diffusion imaging in humans. Magn Reson Med. 1997, 38: 527-533. 10.1002/mrm.1910380404.
Edelman RR, Gaa J, Wedeen VJ, Loh E, Hare JM, Prasad P, Li W: In vivo measurement of water diffusion in the human heart. Magn Reson Med. 1994, 32: 423-428. 10.1002/mrm.1910320320.
Reese TG, Wedeen VJ, Weisskoff RM: Measuring Diffusion in the Presence of Material Strain. J Magn Reson B. 1996, 112: 253-258. 10.1006/jmrb.1996.0139.
Reese TG, Weisskoff RM, Smith RN, Rosen BR, Dinsmore RE, Wedeen VJ: Imaging myocardial fiber architecture in vivo with magnetic resonance. Magn Reson Med. 1995, 34: 786-791. 10.1002/mrm.1910340603.
Reese TG, Heid O, Weisskoff RM, Wedeen VJ: Reduction of eddy-current-induced distortion in diffusion MRI using a twice-refocused spin echo. Magn Reson Med. 2003, 49: 177-182. 10.1002/mrm.10308.
Tseng WY, Reese TG, Weisskoff RM, Wedeen VJ: Cardiac diffusion tensor MRI in vivo without strain correction. Magn Reson Med. 1999, 42: 393-403. 10.1002/(SICI)1522-2594(199908)42:2<393::AID-MRM22>3.0.CO;2-F.
Gamper U, Boesiger P, Kozerke S: Diffusion imaging of the in vivo heart using spin echoes--considerations on bulk motion sensitivity. Magn Reson Med. 2007, 57: 331-337. 10.1002/mrm.21127.
Schmitt M, Potthast A, Sosnovik DE, Polimeni JR, Wiggins GC, Triantafyllou C, Wald LL: A 128-channel receive-only cardiac coil for highly accelerated cardiac MRI at 3 Tesla. Magn Reson Med. 2008, 59: 1431-1439. 10.1002/mrm.21598.
Reese TG, Benner T, Wang R, Feinberg DA, Wedeen VJ: Halving imaging time of whole brain diffusion spectrum imaging and diffusion tractography using simultaneous image refocusing in EPI. J Magn Reson Imaging. 2009, 29: 517-522.
Tournier JD, Calamante F, Gadian DG, Connelly A: Direct estimation of the fiber orientation density function from diffusion-weighted MRI data using spherical deconvolution. Neuroimage. 2004, 23: 1176-1185. 10.1016/j.neuroimage.2004.07.037.
Tuch DS: Q-ball imaging. Magn Reson Med. 2004, 52: 1358-1372. 10.1002/mrm.20279.
Papadakis NG, Xing D, Houston GC, Smith JM, Smith MI, James MF, Parsons AA, Huang CL, Hall LD, Carpenter TA: A study of rotationally invariant and symmetric indices of diffusion anisotropy. Magn Reson Imaging. 1999, 17: 881-892. 10.1016/S0730-725X(99)00029-6.
Skare S, Hedehus M, Moseley ME, Li TQ: Condition number as a measure of noise performance of diffusion tensor data acquisition schemes with MRI. J Magn Reson. 2000, 147: 340-352. 10.1006/jmre.2000.2209.
Jones DK: The effect of gradient sampling schemes on measures derived from diffusion tensor MRI: a Monte Carlo study. Magn Reson Med. 2004, 51: 807-815. 10.1002/mrm.20033.
Tuch DS, Reese TG, Wiegell MR, Wedeen VJ: Diffusion MRI of complex neural architecture. Neuron. 2003, 40: 885-895. 10.1016/S0896-6273(03)00758-X.
Schmahmann JD, Pandya DN, Wang R, Dai G, D'Arceuil HE, de Crespigny AJ, Wedeen VJ: Association fibre pathways of the brain: parallel observations from diffusion spectrum imaging and autoradiography. Brain. 2007, 130: 630-653. 10.1093/brain/awl359.
Gaige TA, Kwon HS, Dai G, Cabral VC, Wang R, Nam YS, Engelward BP, Wedeen VJ, So PT, Gilbert RJ: Multiscale structural analysis of mouse lingual myoarchitecture employing diffusion spectrum magnetic resonance imaging and multiphoton microscopy. J Biomed Opt. 2008, 13: 064005-10.1117/1.3046724.
Gigandet X, Hagmann P, Kurant M, Cammoun L, Meuli R, Thiran JP: Estimating the confidence level of white matter connections obtained with MRI tractography. PLoS One. 2008, 3: e4006-10.1371/journal.pone.0004006.
This study was supported in part by the following grant from the National Institutes of Health: DES (R01 HL093038 and K08 HL079984), VJW (RO1 MH64044), and (NCRR P41RR14075) to the Martinos Center for Biomedical Imaging.
The authors declare that they have no competing interests.
All authors were involved in the previous acquisition of data used in this review. All the authors wrote and/or edited the current article.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Sosnovik, D.E., Wang, R., Dai, G. et al. Diffusion MR tractography of the heart. J Cardiovasc Magn Reson 11, 47 (2009). https://doi.org/10.1186/1532-429X-11-47