Diffusion Tensor Imaging of the Spondylotic Cervical Spinal Cord: A Preliminary Study of Quantifiable Markers in the Evaluation for Surgical Decompression
M Kang, E Anderer, R Elliott, S Kalhorn, P Cooper, A Frempong-Boadu, S Hesseltine, Y Ge, J Marcus, M Law
Citation
M Kang, E Anderer, R Elliott, S Kalhorn, P Cooper, A Frempong-Boadu, S Hesseltine, Y Ge, J Marcus, M Law. Diffusion Tensor Imaging of the Spondylotic Cervical Spinal Cord: A Preliminary Study of Quantifiable Markers in the Evaluation for Surgical Decompression. The Internet Journal of Head and Neck Surgery. 2010 Volume 5 Number 1.
Abstract
Objective: Conventional MRI for the evaluation of cervical spondylotic myelopathy (CSM) may be poorly correlated with a patient’s symptoms and prognosis. Interpretation of canal stenosis, contour deformities, and the presence or absence of T2 signal change, can make it difficult to decide who needs decompressive surgery. Diffusion tensor imaging (DTI) provides quantitative measurements that could help clarify the degree and chronicity of spinal cord disease as a result of compression from degenerative spondylosis. DTI can also detect diffusion abnormalities in areas of acute spondylotic compression occurring without T2 signal change. The purpose of this study is to compare these quantitative DTI measures (i.e. metrics) in patients with severe clinical and radiographic evidence of CSM with controls. Methods: DTI of the cervical spine was performed in 11 patients with severe radiographic and clinical multilevel spondylosis who were planned for surgical decompression versus 10 healthy volunteers (as determined by 2 neurosurgeons A.F.B. and P.R.C., and a neuroradiologist M.L.), using pulsed gradient, double spin echo, echo planar imaging. At the C2-3, C3-4 and C4-5 levels, average FA, MD, E1 (longitudinal diffusion), E2 and E3 (transverse diffusion) were calculated within regions of interest at bilateral anterior, lateral, and posterior regions of the cord. Levels caudal to C4-5 were not analyzed due to artifact on DTI. The average age of the spondylosis patients was 67.2±9.8 years vs. 33.4±15.2 years in the control group (p<.001).
Supported by: Grants RO1CA093992 and RO1111996 from the National Cancer Institute/National Institute of Health.
Introduction
Cervical spondylosis is a common degenerative condition of the spine. While it is present in 75% of people over the age of 65, the vast majority of whom are asymptomatic.6 Conventional T2-weighted MRI is frequently used to analyze cervical spondylosis. However, prior studies report sensitivities of only 15-65% in predicting myelopathy.6 Furthermore, T2 signal change, often indicative of severe compression, appears in the late stages of cord compression and cervical myelopathy limiting the benefit of decompressive surgery. It would be of great benefit to the patient if the anatomic and physiologic changes due to cord compression could be identified earlier and treated before these changes became irreversible.
Diffusion tensor imaging (DTI) is helpful in assessing white matter abnormalities in various brain pathologies including multiple sclerosis, cancer, infection, and ischemic brain injury.9,10,13,15,20,21,24 The use of DTI in spinal cord pathology has been limited due to susceptibility effects of CSF motion, cardio-pulmonary cycles, and surrounding bony structures. Several investigators have explored the use of DTI in the evaluation of spinal cord lesions, including trauma, multiple sclerosis, and spondylosis.6,7,8,12,16 To our knowledge, there is limited published data on using the major and minor eigenvalues from DTI in characterizing spinal cord damage due to spondylosis.
Using individual eigenvalues for cord diffusion, it may be possible to detect very early changes in cord edema before changes are detectable on conventional T2-weighted imaging. This will be of benefit in triaging patients for earlier decompression. We hypothesize that these DTI metrics will be significantly altered in patients with severe spondylosis versus controls. The major eigenvalue (E1) represents longitudinal white matter diffusion and minor eigenvalues (E2 and E3) represent transverse collateral white matter diffusion perpendicular to the long axis of the cervical spine. The ultimate purpose of our study is to compare these eigenvalues, fractional anisotropy, and mean diffusivity in patients with severe spondylosis versus controls.
Terminology
DTI is a relatively new application in spinal cord imaging, not commonly used in the neurosurgical evaluation of cervical spondylotic myelopathy (CSM). In simple terms, DTI provides the ability to detect the direction and ease of water movement in neural tissue. The direction and degree of water movement is not equal in all parts of neural tissues and can also change with age and in pathologic states.11
Methods
Study Subjects
Approval for this study was obtained from the Institutional Board of Research Associates. Informed consent was obtained from all study subjects. DTI of the cervical spine was performed in 10 healthy volunteers and in 11 patients with symptomatic myelopathy with severe multilevel spondylosis on conventional MRI imaging, as determined by two neurosurgeons (A.F.B. and P.R.C.) and a neuroradiologist (M.L.) at our institution. The average age of the spondylosis patients was 67.2 ± 9.8 years (average ± standard deviation) vs. 33.4 ± 15.2 years in the control group. There were 4 males and 6 females in the control group and 6 males and 5 females in the spondylosis group. Control subjects had no history of prior central nervous system disease
MR Imaging and Data Processing
DTI of the cervical spine was performed at 1.5 Tesla using pulsed gradient, double spin echo, echo planar imaging (repetition time [TR]/echo time [TE], 2000/74; 128x128 matrix; 140x140 mm field of view; 10 contiguous 4 mm slices; b= 1000 s/mm2; 2 minutes, 20 seconds acquisition time, using a parallel imaging with an IPAT factor of 2). Leonardo VD10B software on a Syngo VX49B imaging software platform (Siemens Medical Solutions, Erlangen, Germany) was used to process the images. Fractional anisotropy (FA) and mean diffusivity (MD) were calculated using the following formulas:
Six regions of interest (ROIs) were placed at each of the C2-3, C3-4, C4-5 levels, corresponding to the white matter’s left and right anterior, left and right lateral, and left and right posterior regions (Figure 1).
The C5-6 level and below were not able to be analyzed due to artifact. A small distance between the ROIs and the edge of the spinal cord was maintained to avoid CSF partial volume effects. A seventh ROI was placed centrally to cover both the gray and white matter.
Calculations
At the C2-3, C3-4 and C4-5 levels, average FA, MD, and the three principle eigenvalues (E1, E2 and E3) were calculated within paired ROIs at the anterior, lateral, and posterior regions of the spinal cord in the axial plane. These same values were also calculated for the central ROI, which included both gray and white matter. E2 and E3 values were averaged to represent diffusion perpendicular to the axis of the spinal cord. The statistical significances of the computed DTI metrics of gray versus white matter as well as that of the different white matter regions were assessed using a paired, double-tailed Student's
Results
DTI analysis of the spinal cord was done in the axial plane (E2 and E3). There were no significant differences in laterality (left vs. right). Cranial versus caudal diffusion values (E1) were averaged and were not significantly different. Analysis of longitudinal diffusion (E1) showed no significant differences between spondylosis patients versus controls (Table 1).
Figure 3
The minor eigenvalues (E2 and E3) showed strikingly significant differences between spondylosis patients and controls in all regions of interest at C4-5 (Table 2).
Figure 4
In addition, significant differences were also seen in the posterior, lateral, and central ROIs at C2-3. Significant differences in fractional anisotropy were seen only at the C4-5 level where, like the minor eigenvalues, the ROI’s had the greatest magnitude of difference in the posterior and lateral regions (Table 3).
Figure 5
Mean diffusivity measurements showed significant differences between control and spondylosis patients at C2-3 in the lateral and central ROIs and at C4-5 in the anterior and central ROIs (Table 4).
Figure 6
At C4-5, correlations of age as the independent variable and FA or MD as the dependent variable were performed. FA in both groups showed a decrease, which was greater in the spondylosis patients as a function of age (Figures 2 & 3).
Mean diffusivity showed a gradual increase in the normal patients (Figure 4) and a trace decrease as a function of age in the spondylosis patients (Figure 5).
Discussion
Several authors have speculated on the meaning of signal change on conventional MRI stating that it could represent edema, inflammation, gliosis, or myelomalacia depending on the nature and timing of spinal cord disease.19 DTI of the cervical spinal cord provides quantitative information, which can clarify the nature and chronicity of this signal change. This information in combination with a clinician’s evaluation of stenosis, contour deformity, and signal change (if present) on MRI could improve the sensitivity of using image findings to predict the clinical manifestations of CSM. For instance, DTI can detect abnormal diffusion in areas where T2 signal change is absent potentially allowing clinicians to identify surgical candidates earlier.5 In a study by Facon et al., 8 DTI provided information regarding the timing of compressive spondylotic spinal cord disease. They reported that FA and ADC (apparent diffusion coefficient) values could differentiate acute versus chronic compression of neural tissue. For example, FA, which measures the degree of anisotropic flow within white matter tracts, was found to be increased in early spinal cord compression and decreased in later stages of compression. In early compression, extracellular water is displaced into cells leading to decreased ADC and higher FA values as a result of cellular edema and tighter packing of axon bundles. These changes model the diffusion changes of the brain in acute versus chronic infarct. With prolonged compression and neuronal death, there is cell lysis, demyelination, gliosis, cystic degeneration, edema, and Wallerian degeneration leading to an outflow of water into the expanded extracellular space leading to increased ADC and lower FA values.8,14,27,32,33,34 Thus, CSM patients with imaging evidence of compression who present with higher FA and lower ADC values have acute compression may be more likely to recover with surgery than patients with similar imaging with low FA and high ADC values likely representative of gliosis or myelomalacia.
In our study, FA did not differ from controls at C2-3 or C3-4, but showed significant decreases in value versus controls at the C4-5 level, which most consistently demonstrated radiographic spinal cord compression in the levels analyzed in our disease group. This data is encouraging as it implies that DTI metrics correlate strongly with known location of disease. Mamata et al.16 made a similar finding and noted that the average ADC increases and FA decreases significantly versus controls in high signal lesions, due to compression, on T2-weighted images. It is important to note that one of the differences between our study and that of others cited in this paper is that we used axial ROIs while others used sagittal images.8,16 Using sagittal images, ROIs tend to include more grey matter in measurements of intended white matter diffusion.16 Hence, depending on the imaging technique (axial versus sagittal), DTI values which significantly differentiate normal from abnormal may vary from study to study.
We found that mean diffusivity significantly differed at the anterior and central ROIs of the C4-5 level. MD also showed significant changes at the lateral and central regions of C2-3. This finding is unusual in that spondylotic compression generally does not occur at this level. These MD values are likely a product of CSF distortion as the naturally widened canal in the upper cervical spine allows for more CSF pulsation artifact to affect volume averaging when making quantitative measurements in the cord.8 As the canal narrows at C3-4, this effect is minimized. It is also possible that there is an element of retrograde Wallerian degeneration, where spondylosis and myelopathy at a lower level within the cervical spine may result in changes along the spinal tracts within the spinal at more rostral levels. Thus, in using DTI of the spine, the clinician must analyze all metrics (i.e. FA, ADC, eigenvalues, etc.) in conjunction with T2 sequences and the clinical exam to avoid being misled by data that may have been distorted by artifact.
The preservation of E1, or longitudinal diffusion, reflects the maintenance of axonal integrity.26 The often slow microvascular and microstructural changes associated with cervical spondylosis may be insufficient to significantly affect longitudinal diffusion along axons. This makes sense as CSM usually manifests insidiously over several years to decades with progressive motor and sensory disturbances and not acute paraplegia. Animal evidence has shown that it takes a powerful neurotoxin like Cuprizone (a copper chelator which destroys myelin) in order to decrease E1 in the mouse corpus callosum.28 Unlike this toxin, degenerative cord compression causes a lesser and more gradual degree of damage to the white matter tracts and is reflected by significant changes in the minor eigenvalues, E2 and E3, which represent transverse water diffusion within the spinal cord.11,26 The magnitude of significant differences in minor eigenvalues at C4-5 was greater than that of any other DTI metric used in this study—again consistent with the presence of maximal spinal cord compression at the C4-5 level in the levels analyzed in our study group.
DTI metrics may provide novel biomarkers in the radiographic evaluation of CSM. The goal of this paper is to examine the relationship between these markers and a group of patients with clinical and conventional MRI findings of CSM. With that in mind, we need to discuss the shortcomings of this study. We were not able to analyze the cervical spinal cord below C4-5 due to coil design and pulsation artifact. Several studies have shown, however, that C4-5 is commonly affected in multilevel cervical spondylotic myelopathy and all 11 patients in our study had compression at that level.4,18,29
There is also a statistical difference between the ages of the patients in the spondylosis group versus the control group (p < 0.001). The mean age for the spondylotic patients is 67 years versus 33 years for the controls. When reviewing our database it was not possible to find a 67 year-old patient with a normal spine. There is usually, if not always, an element of spondylosis in the aging patient with varying degrees of stenosis. To ensure that we were comparing patients with stenosis with those having nonstenotic canals, our controls had to be younger patients. Investigators have shown a slight age-related decline in FA.16 However, the degree of change between the two groups in our study is far greater than that seen in prior studies in which the effect of age was examined. As demonstrated in Figures 2 and 3, in both the control and spondylosis groups, FA showed an inverse relationship with age. MD showed an age-related trace increase in the control group, but marginally decreased in the spondylosis group. This implies that although the average ages between our study groups significantly differ, age does not seem to independently affect the DTI values within the framework of our analysis.
A second limitation is the relatively small number of patients and controls studied. We are in the process of combining multi-center datasets to increase the number of patients and this will be published in a separate paper. Finally, we were not able to accurately analyze the cervical spinal cord below C4-5 (where there is usually maximal cord compression) due to coil design and pulsation artifact. Improvements in coil design as well as faster-navigated pulse sequences are being developed which will allow examination of the lower cervical spinal cord.
Conclusion
The uniform increase in the minor eigenvalues in cervical spondylosis patients represents a novel biomarker of this disease. These biomarkers can help elucidate the anatomic and physiologic changes in the spinal cord parenchyma due to spondylotic compression earlier and describe the nature of spinal cord signal change better than routine T2-weighted imaging. Currently, there are efforts to overcome the technical limitations of DTI imaging of the spine at our institution (e.g. coils to image the entire cervical spine in the axial plane and sequence modification to minimize hardware artifact on post operative imaging). As these limitations are overcome with future innovations, DTI of the spine could become a powerful adjunct in the diagnosis, triage, and surgical treatment of CSM.