Abstract
BACKGROUND AND PURPOSE: Diffusion tensor and diffusion-weighted spinal cord imaging remain relatively unexplored techniques despite demonstrations that such images can be obtained and may yield clinically relevant findings. In this study, we examined the temporal dynamics of spinal cord motion and their impact on diffusion tensor image quality.
METHODS: Four healthy volunteers underwent phase contrast-based velocity mapping and segmented echo-planar diffusion tensor scans of the cervical spinal cord. Regions of interest in the cord were used to identify the temporal patterns of motion. The delay of data acquisition after the cardiac trigger was varied to correspond to either quiescence or motion of the cord.
RESULTS: The cervical spinal cord consistently displayed maximal velocities in the range of 0.5 cm/s and accelerations of up to 25 cm/s2. In both these respects, the cervical cord values were greater than those of the medulla. Despite this pronounced motion, approximately 40% of the cardiac cycle can be described as relatively calm, with absolute velocities and accelerations less than 20% of the maximum values. Confining image acquisition to this window reduced ghosting artifacts and increased the consistency with which the dominant direction of diffusion was along the rostral-caudal axis in both gray and white matter of the spine. Preliminary clinical application and fiber tracking in pathologic cases was feasible, and alterations of the diffusion properties by multiple sclerosis lesions, tumor, and syringomyelia were seen.
CONCLUSIONS: Acquiring DTI data during the quiescent phase of spinal cord motion can reduce ghosting artifacts and improve fiber tracking.
Diffusion tensor imaging (DTI), and diffusion-weighted imaging (DWI) noninvasively provide information about structural integrity at a cellular level that is of significant clinical and research interest.1–8 Although demonstrated in application to the spinal cord by several groups, the clinical roles of both DTI and DWI remain poorly established.9–16 To date, demonstrations of fiber tracking in the spinal cord have been especially limited.17 It is our interest to expand the use of DWI and DTI in the spinal cord, where alterations in tract morphology may be relevant to injury diagnosis and in light of the arrival of new therapeutic measures in animal models, to future monitoring of recovery.
Three factors have contributed to the slow advancement of diffusion imaging in the investigation of patients with spinal cord pathology.18–20 First, magnetic field inhomogeneities as a result of nearby bone, lung, and body contour at the neck are exacerbated by the use of long echo times to allow diffusion encoding. Second, spinal imaging coils, because of their primarily planar geometry, are markedly less sensitive at the depth of interest (ie, the spinal cord) than the volume coils that are used for the brain. Third, the spinal cord and surrounding CSF undergo significant pulsations, whereas surrounding tissues are also prone to respiratory and swallowing motions, all of which can corrupt the motion-sensitive diffusion encoding and imaging processes. Because of these factors, single-shot echo-planar images (EPI) of the spine on clinical and high-field imaging systems are generally severely distorted. Multishot, echo-planar, and fast spin-echo sequences have yielded better results for diffusion imaging of the human spinal cord.21–24 On the other hand, strategies to restrict the imaging volume with25 or without the use of phase encoding13,26–29 have been used but are not yet routinely available for clinical use. All these techniques reduce the number of echoes acquired and thus the duration of the sampling window, leading to a reduction in susceptibility artifacts. Spin-echo based approaches have the advantage of better inherent refocusing of T2*-related signal intensity losses.17 Perhaps as a result of shorter acquisition train duration, EPI has been found to outperform fast spin-echo when a multishot strategy is used for DWI of the spinal cord.22
A problematic aspect of multishot acquisitions is the need to ensure that different acquisition “shots” for a given image undergo similar degrees of corruption by motion and that images sensitive to different directions of motion are not differentially affected. The dynamics of the spinal cord consist typically of caudal displacement starting soon after arterial systole in the cranium, and a return of the cord to its original position during diastole.30,31 This dependence on the cardiac cycle has been exploited through the use of cardiac gating.32 Further improvements can be achieved by the incorporation of a flow-compensated gradient pattern33 and correcting residual phase shift differences between shots by navigator echoes.12,34,35
Cardiac triggering is recognized to reduce signal intensity drop-out in single-shot diffusion-weighted acquisitions,36 and the addition of navigator correction is an essential step in segmented EPI data acquisition.32,37 Bammer et al38 and Holder et al10 combined cardiac-gated and navigator correction to multishot DTI of the spinal cord. However, they did not restrict the acquisition window to ensure that all data were acquired during comparable conditions of spinal cord motion. Thus, although cardiac cycle modulation of motion was avoided, they did not account for the possibility that different periods of the cardiac cycle may be more prone to motion than others and hence less desirable for data acquisition. Within the brain, the timing of data acquisition relative to the cardiac cycle has recently been shown to influence the quality of DTI studies.39 Herein, we consider the temporal dynamics of the CSF and spinal cord motion as a factor in optimizing cardiac-triggered, navigator-corrected, multishot, spin-echo echo-planar DTI of the human spinal cord. In particular, we aimed to identify an optimal period for data acquisition in the upper cervical region of the spinal cord by examining cardiac-related movement of the spinal cord and seeking to establish whether different trigger window timings affected DTI image quality or the results of fiber tracking.
Methods
Subjects
Spinal cord dynamics were examined in 4 healthy adult volunteers free of neurologic symptoms (3 men, 1 woman; mean age, 38 years [SD = 6.4 years]). Details of the motion and diffusion-sensitive imaging sequences and postprocessing are given below. A conventional T2-weighted, sagittal fast spin-echo acquisition (echo time [TE]/repetition time [TR]/echo-train length, 80 ms/3000 ms/13) was performed to rule out any pathologic conditions. DTI was then performed using the optimized acquisition strategy targeted to the quiescent period and again with the trigger delay set to begin acquisition during the motion of late spinal systole/early diastole.
In addition, the timing of the optimized triggering strategy obtained from the healthy subjects was applied in 5 patients with cervical spinal cord pathology undergoing spinal MR imaging investigation. In this subgroup, measurement of the spinal cord dynamics was performed in consideration of the time required for the routine clinical imaging. The patient referrals included 2 with multiple sclerosis, 1 with cervical syringomyelia, 1 with cervical intramedullary glioma, and 1 with dolichoectatic vertebrobasilar system compressing the medulla oblongata and the cervicomedullary junction. Each subject provided informed consent in keeping with the practices of the local institution.
All imaging was performed on a 3T MR scanner with high-performance gradient system (Gmax: 80 mT/m, Gslew: 200 T/m/s). A phased-array spine coil was used for reception, using the 4 cervical coil elements. Parallel image acquisition and reconstruction (eg, SENSitivity Encoding40) were not supported with this coil array at the time of the study.
Spinal Cord Movement
Acquisition.
A midline sagittal velocity mapping study of the superior-inferior component of velocity was performed using a retrospectively triggered phase-contrast sequence. The velocity encoding range was set to 8 cm/s, and 36 time points were reconstructed per cardiac cycle. The imaging parameters were: TE/TR/flip angle/number of excitations (NEX), 5.8–6.1 ms/20 ms/60°/3, field of view [FoV]/acquired resolution [AcqRes]/reconstructed resolution [ReconRes]/section thickness −200 mm/100 × 256/332 × 512/5 mm.
Image Analysis.
Using the scanner’s built-in analysis software, rectangular regions of interest (ROIs) containing at least 40 voxels were drawn in the spinal cord at the levels of the medulla and of the C2 and C6 vertebral bodies on the first frame of the velocity mapping study of each subject. The regions were copied to the following frames, and spatial mean velocity within each ROI was determined for each frame. The time course of the velocity was assessed to identify the onset and end of the longest plateau of velocity. As the criteria for identifying this range, we adopted a velocity change of 20% (typically ± 0.06 cm/s) of the maximal values. That any concomitant gradient and eddy current-induced errors would not alter the identified plateau was confirmed by the drawing of a background ROI, but no background correction was made. The time-course data were exported as a text file for preparation of figures.
Diffusion Tensor Imaging
Acquisition.
The diffusion tensor acquisition closely followed the protocols of Bammer et al38 and Holder et al.10 The sequence was cardiac-gated using the finger pulse-oximeter signal intensity. Spatial presaturation bands were applied outside all 4 edges of the FOV to reduce contamination from structures outside the spinal column. The sequence included a navigator echo for correction of patient motion over the course of the scan. The imaging parameters were: TE/TR/flip angle/NEX, 66–72 ms/2 or 3 RR intervals/90°/3; FOV/AcqRes/ReconRes/section thickness, 11 cm/82 × 96/128 × 128/6 mm. Diffusion encoding was performed in 6 directions with a maximal b-factor of 650 seconds/mm2. Previous studies with b-values of ≤750 seconds/mm2 have shown the strong anisotropy difference between spinal gray and white matter to be well resolved41 and have the advantage of allowing a shorter echo time than with the higher b-values typically used in DTI of the brain. This gradient scheme, standard on our system, always used at least 2 gradient coils simultaneously. This rotates the diffusion-encoding axes away from the imaging axes, making interpretation of the source images more difficult but reducing the echo time as much as possible. In addition, for each section, a baseline image with minimal diffusion weighting (b <20 seconds/mm2) was acquired.
Timing and duration of the DTI studies were dictated by the phase-contrast results. In each subject, the above multisection, axial diffusion tensor study was performed once with the gating window set to acquire the sections during quiescence, and once with acquisition including the motion of late spinal systole/early diastole. These corresponded to trigger delays of 290 to 320 ms for quiescence, and the system’s shortest trigger delay (90 ms) respectively. In both cases, the data acquisition window was limited to 350 ms (consistent with the minimum quiescent interval observed); allowing 2 sections was acquired in each cardiac cycle. The order of the scans alternated between subjects. To minimize cross-talk between sections, images were acquired in an interleaved order, with 2 passes through the sections, allowing the scanner software to acquire the odd-numbered sections first and then the even-numbered sections. Because 8 or 12 sections were planned, the TR was correspondingly set to either 2 or 3 RR intervals. Total imaging time was typically 9–12 minutes depending on the heart rate.
The DTI data were first reconstructed as the individual diffusion-weighted images, the reference T2-weighted (“b = 0”) EPI. After the reconstruction process described below, fractional anisotropy (FA) maps and mean diffusion-weighted images were generated. In the volunteers, for whom DTI acquisition was performed with cardiac triggering set in one instance to the scanner default (approximately 60 ms) and in one to the quiescent period of the cardiac cycle, 2 sets of the above DTI data were obtained. Pulsation artifacts were less obvious with gating to the quiescent window in all subjects. Because the ghosts tended to overlap other anatomic structures, however, we could not define an unambiguous quantitative metric of the change in the level of ghosting. An illustrative example of the difference in gating strategies is therefore provided. With the patients, only the optimized gating approach was performed. In assessing these, we focus therefore on the correspondence of the cross-sectional DTI with the conventional MR images.
Fiber Tracking.
Before performing fiber tracking, the data were transferred to an off-line workstation. Dewarping of eddy current-induced image deformations was performed with a correlation-based 2D-affine registration algorithm.42 The independent elements of the diffusion tensor were obtained on a pixel-by-pixel basis by singular value decomposition, using the diffusion-weighted images and the control images with b <20. After diagonalization, the eigenvalues and eigenvectors were determined, and FA maps and color-coded orientation field maps43 were created.
Tracking and subsequent visualization were performed using custom-made software. A line propagation algorithm based on FACT (Fiber Assignment by Continuous Tracking7) was the basis for reconstructing fiber pathways. Stopping criteria based on FA values were set to be either intermediate between gray and white matter (to selectively track in white matter), or at the lower range of values measured for gray matter (to allow tracking through both white and gray matter). Two 2D ROIs were manually defined, one encompassing the spinal cord on the most superior section and the other on the most inferior section. For each ROI, a user-specified number of seed points per voxel was generated; from each of these seed points was the line propagation algorithm launched. For all studies, 20 seed points per voxel were specified, and these were randomly distributed over the voxels of the ROIs. The resulting fibers were generated and saved separately. The 2 ROIs were then used in conjunction for a 2 ROI analysis of fibers propagating from end to end within the imaged volume. This last approach produced a subset of the fibers identified in the single ROI seeding but with reduced visibility of emerging nerve roots. The results reported therefore focus on the single and combined single ROI results.
Results
Spinal Cord Movement
Velocity-sensitive imaging of the spinal cord showed very good consistency between subjects for the overall pattern of spinal cord motion. The velocity time-course curves for each cervical level are shown in Fig 1 along with the derived acceleration time-course curves. Although there was no prolonged period of spinal cord stasis, between 36% and 56% of the cardiac cycle could be described as quiescent, showing velocities less than 0.06 cm/s. Over these time windows, the accelerations were also lowest, typically less than 5 cm/s2. Normalizing the time courses showed the onset and duration of the quiescent period as a fraction of the cardiac cycle to be strongly stereotypical across subjects (Fig. 1G, -H).
The timing and velocity characteristics of the 3 cervical spinal levels examined are summarized in the Table. Progressing from superior to inferior, the onset of quiescence is seen to be progressively delayed, whereas its duration became progressively shorter. At the same time, the maximum velocities achieved over the cardiac cycle progressively increased from less than 2 mm/s to more than 5 mm/s. The subject heart rates were in a narrow range (67–81 beats/min).
Timing of DTI Data Collection
DTI data acquired during a mix of highly motional and quiescent periods of the cardiac cycle and strictly during the quiescent period are exemplified by the results shown in Figs 2 and 3, respectively. Three features are apparent from the acquired images. First, gating greatly reduces CSF-related ghosting in the T2-weighted, low b-value image used as a reference. Second, although ghosting of the spinal cord is not a prominent feature, the average intensity and the consistency of the cord appearance across sections and b-values is much higher in the images acquired entirely during the quiescent period. Finally, there is a notable ghosting of residual extraspinal tissues, particularly subcutaneous fat. To reduce the impact of ghosting from fat, subsequent DTI acquisitions were performed with fat suppression.
Diffusion Tensor Imaging
Examples of DTI and tracking in the spinal cord are given in Figs 4–8, illustrating some preliminary observations based on the small number of subjects studied to date. Acquisition during quiescence results in greater consistency of the dominant direction of diffusion tensor laying in the rostral-cranial throughout the spinal cord. This is seen in the color-coded FA maps and first eigenvector displays from a healthy control subject (Fig 4). In all control subjects and patients, the spinal cord was clearly distinguished from adjacent fluid in the DTI data and maps. The expected butterfly configuration of the central gray matter however, was visible on most sections for each volunteer only with acquisition timed to the quiescent period.
Fiber Tracking.
Fiber tracking results representative of a healthy adult obtained using the optimized trigger timing detailed above are shown in Fig 5. When the ROI contains the entire cross-section of the spinal cord in the most superior (or inferior) section, it is notable that relatively few fibers extend from end to end of the comparatively short imaging volume (typically 2-3 vertebral bodies in the current implementation). At least part of this attrition appears to be due to fiber convergence and crossing at the emergence of nerve roots (Fig 5B). In light of the residual errors in the eigenvectors seen in Fig 4, however, one cannot rule out the presence of artifactual attrition as well. By limiting the initial ROI to the lateral white matter in the most superior section, a subset of the tracks can be more clearly seen. These tracks remain lateralized, and with a suitable stopping threshold, the depicted fibers do not appear to enter the central gray matter of the spinal cord. Blue in Fig 5B, -C, depicts the cranial-caudal orientation of the local first eigenvector, so the dominant direction of diffusion is readily apparent. Intermittent zones of green and red correspond to the origins and progressions of nerve roots.
In Fig 6, DTI and tracking results from a patient with hydromyelia are shown. We note that the diffusion-weighted images show signal intensity reduction in the dilated central canal that are not as pronounced as for the CSF, suggesting a lower diffusion coefficient in this region. The timing of the acquisition window (320 ms after peripheral pulse) was based on the results in controls. In this patient, ghosting of the spinal cord was evident on some sections, possibly as a result of altered spinal cord dynamics associated with its pathologic condition. No such ghosts were seen in the volunteers. Over the cross-section of the cord tissue, the direction of the first eigenvector was rostral-caudal, suggesting that the data were adequate for fiber tracking. Making use of a surface-shaded display of the central canal and high-intensity CSF based on the T2-weighted EPI images (low diffusion encoding), it is apparent that few fiber tracks extend into either of these regions (Fig 6D, -E). We attribute these incursions to differences in the effective boundaries generated in the processes of surface rendering and evaluating the stopping criteria of low FA.
In the patients with MS, the DTI data were suggestive of localized differences in diffusion properties as illustrated in Fig 7. The location and extent of the lesions, however, was not readily visible from fiber tracking. Conversely, in a patient with an intramedullary glioma (Fig 8), the diffusion tracking results are clearly disrupted in proximity with the lesion. As with the syrinx, altered spinal cord dynamics may also play a role here, in that the image quality was exceptionally good, clearly depicting the expected pattern of gray and white matter in the mean diffusion-weighted, apparent diffusion coefficient, and FA maps at levels below the tumor. Immobilization of the thickened cord within the spinal canal was probably the important factor in achieving the good quality images in this patient. The presence of edema cranially and caudally to the tumor may also have accentuated the FA differences between gray and white matter.
Discussion
We have demonstrated a strategy for optimizing cardiac synchronization of the acquisition of multishot echo-planar, DTI of the spinal cord. The strategy takes advantage of a period of relative quiescence in the motion of the spinal cord to reduce artifactual signal intensity losses in the diffusion-weighted images and to reduce pulsation ghosting in the low b-value images acquired. This leads to visible improvement in orientation of the resulting diffusion eigenvectors, and the depiction of cross-sectional spinal cord diffusion properties. In addition, we present some initial results of fiber tracking in the spinal cord.
The motion of the cervical spinal cord we have observed is in very good agreement with the previous MR studies.44–46 Of particular interest here is the intraindividual consistency, across our small group of volunteers, of a period of relative quiescence lasting slightly less than half of the cardiac cycle and beginning roughly 320 ms after the trigger signal intensity. The use of peripheral triggering may be responsible for some of the variability between subjects in the onset of the quiescent window (Fig 1) in that the waveform does not have a sharp triggering feature, will have different delays relative to spinal systole depending on subject size, and is prone to variations due, for example, to cooling of the subject’s arms. Although ECG triggering may reduce this variability, an element of intersubject variation will arise from subject age, height, details of their skull-base anatomy, and pathologic conditions47 and will differ according to the spinal level being examined.30,31 For our purposes, a key advantage of peripheral triggering was that the quiescent window was consistently framed between the peripheral triggers, so that the split in the data acquisition required for prospective triggering coincided with spinal motion (ie, time during which we did not wish to acquire data). With ECG gating, part of the quiescent window would have been sacrificed to allow the scanner to unambiguously detect subsequent trigger signals, with a further loss in scan time efficiency.
We used the average onset and duration of the quiescent window observed in the volunteer studies as a basis for examining patients with different pathologic conditions of the spinal cord. This resulted in good quality DTI images, with consistent cranio-caudal orientation of the primary eigenvectors within the white matter of the cord. The presence of ghosting, specifically of the CSF and cord, in our DTI studies of patients with space-occupying lesions (eg, Fig 6), suggests that alterations in the spinal cord dynamics are present in these cases. Focal increases in cord velocity have previously been reported by Tanaka et al48 in cases of cord compression associated with lower extremity motor function impairment, whereas Levy49 has documented decreases in CSF and cord motion in symptomatic cord compression and increases during secondary atrophy. Ideally, therefore, timing of DTI data acquisition should be based on velocity measurements performed on an individual basis.
In a clinical setting, the need for broad anatomic coverage along the spinal cord within a scanning time acceptable to patients raises a conflict between spatial resolution and attaining adequate signal-to-noise ratio (SNR). With the combination of cardiac triggering, data acquisition limited to part of the cardiac cycle, and multishot acquisition, the resultant scan time was roughly 12 minutes for DTI. Sagittal imaging planes may be used to extend the range of the spinal cord examined, but with currently demonstrated section thicknesses, this exacerbates partial-volume mixing both between gray and white matter and between the cord and CSF in the through-plane direction. As a result of cord curvature, this partial volume effect may be difficult to distinguish from pathologic processes. Equally, such partial volume effects are likely to have a detrimental effect on fiber tracking results. We have therefore focused on axial imaging planes for the present work.
Our preliminary experience shows that limiting data acquisition to spinal cord quiescence increased the fraction of spinal cord voxels for which the first eigenvector of the diffusion tensor had the expected rostral—caudal orientation throughout both the gray and white matter of the cervical spine. However, because axons may extend the entire length of the spinal cord, with shorter axons extending from the level of their point of entry or exit through spinal foramina to the brain,50 it is notable that relatively few fibers propagated from end to end in the comparatively short imaging volume (typically less than 3 vertebral bodies). It is not clear to what extent this is due to noise in the data and to the staircase effect from the large section thickness coupled with the curvature of the spinal cord. In some cases (eg, Fig 5), however, the emergence of nerve roots and the associated partial volume effects of crossing and bending fiber tracts seem to be likely to play a role. These issues are the focus of considerable attention in connection with intracranial fiber-tracking.51–56 In the presence of the overwhelming diffusion anisotropy, a second eigenvector representation may be better suited to demonstrating the sub-structures of the spinal cord. Lee and Alexander57 have recently shown that diffusion encoding patterns may be optimized to improve the estimates of the minor eigenvectors
The good correspondence between the DTI results and the pathologic appearances in conventional MR imaging in the limited pathologic cases studied so far hints at the possible future clinical application of spinal cord DTI. DTI seems to provide a useful window in the depiction of spinal cord anatomy, through the distinction of gray and white matter (ie, the butterfly-shaped gray matter configuration) within the cervical spinal cord. Indeed, several authors have already presented results from individual patients and small patient groups investigating the nature of diffusion and diffusion tensor changes in specific pathologic condtions, including multiple sclerosis and cervical spondylosis.16,23,58–61 In the presence of pathologic conditions, the ability to examine both diffusivity and anisotropy changes within the both gray and white matter compartments is of interest. A shortcoming of the present technique, therefore, is that despite the measures taken and the improvement relative to simple cardiac gating, clear differentiation of gray and white matter was not achieved in all sections in all subjects. Imperfect fat suppression in the present study resulting in the overlay of chemical shift of unsuppressed paraspinal fat signals onto the spinal cord or the inconsistent suppression of lipid signals within the spinal cord signals is a possible contributor to the lack of contrast. Subsequent investigations of quadratic phase radio-frequency pulses have shown improved saturation behavior, allowing the acquisition FOV to be more closely targeted to the spinal cord.62 Higher order shimming63,64 and the placement of materials around the neck to improve the magnetic field homogeneity65 may also prove helpful.
We consider it more likely however, that motion caused by respiration over the course of the scan, or due to the remaining slight movements of the spinal cord during the quiescent period, are responsible for the inconsistent depiction of gray-white contrast. Spuentrup et al66 have demonstrated the advantages of respiratory gating via navigator echoes for diffusion imaging of the spinal vertebrae.31 Residual ghosting artifacts attributed to respiration suggest that navigator echoes did not correct for all residual motion in the present study. A possible future direction therefore would be to account for both the respiratory and cardiac cycles during data acquisition, though in general this would be expected to prolong scan times. Alternatively, the navigator process may be modified to better suit correcting the motion of the small, low-contrast spinal cord that is not rigidly coupled to the surrounding, larger tissue mass, as has been done for coronary imaging studies. Auto-focus algorithms,34,67 inherently less motion sensitive sequences,68 and parallel data acquisition strategies40,69,70 may also be avenues for reducing motion artifacts in spinal cord DTI.
We have not investigated the comparative scan time or SNR efficiencies between our gated approach and other acquisition strategies. Limiting data acquisition to the portion of the cardiac cycle during which the spinal cord is quiescent reduces the scan efficiency to achieve the benefit of reduced artifact levels. Two and a half times as many averages of ungated data could potentially be performed in a comparable scan time, but this is unlikely to be sufficient to average out the spinal cord motion. By rejecting roughly a third of ungated averages on the grounds of motion-related reduction of their navigator echo intensities, Wang et al71 demonstrated ghost reduction and improved image homogeneity in brain DTI studies. If the same level of performance was achieved in the spine, this would allow a scan time reduction of approximately one third compared with the present study. Based on the cord motion time-courses, however, we expect higher levels of motion contamination in the spine, which would reduce this temporal advantage. Improving the SNR of the acquired data, by so-called PROPELLER acquisitions, which have the benefit of reinforcing the SNR by repeated sampling of the center of k-space,72,73 may allow the number of averages to be reduced. Alternatively, targeting image acquisition to the spinal canal by volume-selective or FOV limitation can help with the time reduction but tend to compromise SNR.13,26–29,74 A compromise on repetition time and imaging volume size with a 3D acquisition may also bring signal intensity averaging benefits.75 With the short readout train of just 9 echoes, susceptibility artifacts were not significant on axial acquisitions in the present study. Accepting the intrusion of susceptibility artifacts would allow more echoes to be acquired per shot, thus reducing the number of shots required and overall imaging time. Finally, because the present work was performed with a conventional quadrature array spine coil, the growing availability of coils and receiver chains suitable for parallel data acquisition (eg, SENSE, GRAPPA, ASSET, etc) opens a further avenue to reducing scan times.
Conclusion
In this preliminary study, we have attempted to demonstrate the specific utility of optimizing the choice of cardiac triggering delay with the aim of reducing motion-related artifacts in DTI studies of the spine. Considerable challenges remain if a greater range of coverage along the spinal cord is to be achieved while maintaining or improving the high spatial resolution and contrast properties of the present study. The use of high performance gradient systems, more sensitive coil geometries, and refinements to the data acquisition strategy may contribute to this progress. Data acquisition during the quiescent phase of spinal cord motion should be used to minimize ghosting artifacts and improve fiber tracking.
Acknowledgments
We express our gratitude to Dr. Gerard Crelier for the use of software for diffusion tensor visualization and analysis.
Footnotes
This study was funded by the NCCR on Neural Plasticity, and Repair.
Previously presented in a preliminary form at the Annual Meeting of the American Society of Neuradiology, June 5–11, 2004; Seattle, Wash.
References
- Received September 29, 2005.
- Accepted after revision January 24, 2006.
- Copyright © American Society of Neuroradiology