Next Article in Journal
Eosinophilic Esophagitis: Cytokines Expression and Fibrotic Markers in Comparison to Celiac Disease
Next Article in Special Issue
Dynamic Brain Connectivity in Resting-State FMRI Using Spectral ICA and Graph Approach: Application to Healthy Controls and Multiple Sclerosis
Previous Article in Journal
Chest CT Severity Score and Systemic Inflammatory Biomarkers as Predictors of the Need for Invasive Mechanical Ventilation and of COVID-19 Patients’ Mortality
Previous Article in Special Issue
Custom 3D fMRI Registration Template Construction Method Based on Time-Series Fusion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Utility of Diffusion and Magnetization Transfer MRI in Cervical Spondylotic Myelopathy: A Pilot Study

1
Department of Rehabilitation Medicine, VHS Medical Center, Seoul 05368, Korea
2
Department of Radiology, VHS Medical Center, Seoul 05368, Korea
3
Department and Research Institute of Rehabilitation Medicine, Yonsei University College of Medicine, Seoul 03722, Korea
4
Department of Orthopaedic Surgery, Hallym University Sacred Heart Hospital, Anyang 14068, Korea
5
Department of Physical Medicine and Rehabilitation, Hallym University Sacred Heart Hospital, Anyang 14068, Korea
*
Author to whom correspondence should be addressed.
Diagnostics 2022, 12(9), 2090; https://doi.org/10.3390/diagnostics12092090
Submission received: 28 July 2022 / Revised: 26 August 2022 / Accepted: 28 August 2022 / Published: 29 August 2022
(This article belongs to the Special Issue Advancements in Neuroimaging)

Abstract

:
Diffusion tensor imaging (DTI) and magnetization transfer (MT) magnetic resonance imaging (MRI) can help detect spinal cord pathology, and tract-specific analysis of their parameters, such as fractional anisotropy (FA), mean diffusivity, axial diffusivity (AD), radial diffusivity (RD) and MT ratio (MTR), can give microstructural information. We performed the tract-based acquisition of MR parameters of three major motor tracts: the lateral corticospinal (CS), rubrospinal (RuS) tract, and lateral reticulospinal (RS) tract as well as two major sensory tracts, i.e., the fasciculus cuneatus (FC) and spinal lemniscus, to detect pathologic change and find correlations with clinical items. MR parameters were extracted for each tract at three levels: the most compressed lesion level and above and below the lesion. We compared the MR parameters of eight cervical spondylotic myelopathy patients and 12 normal controls and analyzed the correlation between clinical evaluation items and MR parameters in patients. RuS and lateral RS showed worse DTI parameters at the lesion level in patients compared to the controls. Worse DTI parameters in those tracts were correlated with weaker power grasp at the lesion level. FC and lateral CS showed a correlation between higher RD and lower FA and MTR with a weaker lateral pinch below the lesion level.

1. Introduction

The incidence of cervical spondylotic myelopathy increases with age. The clinical symptoms are diverse, and researchers have obtained various imaging findings associated with the disease [1]; however, the symptoms and imaging findings do not necessarily match [2,3]. Researchers recommend diffusion tensor imaging (DTI) and magnetization transfer (MT) magnetic resonance (MR) imaging (MRI) for the detection of spinal cord pathology that cannot be identified with conventional MRI.
DTI and MT MRI provide variable parameters that enable quantitative studies. Fractional anisotropy (FA) and mean diffusivity (MD) show good consistency and are the most frequently used DTI parameters. FA reveals both demyelination and axonal injury; however, they are not specific to the underlying cause. MD measures membrane density and is sensitive to cellularity, edema, and necrosis. FA values decrease and MD values increase in injured or diseased spinal cords [4,5,6]. Axial diffusivity (AD) and radial diffusivity (RD) are also used; however, they are less frequently used than FA and MD because of their deceptive change in complex tissue architecture and mutual effect on each other. AD specifically reveals axonal content. AD decreases with acute traumatic spinal cord injury because of axon fragmentation, while it increases over time in chronic neurodegenerative disease because axon fragments are cleared. RD reveals the myelin status relatively specifically and increases with demyelination [7,8,9,10].
MT MRI provides an indirect surrogate for myelin content and is specific to demyelination and degeneration. The MT ratio (MTR) is calculated as the relative change in intensity when the images are acquired with and without an off-resonance MT pulse. MTR in the white matter (WM) is reduced in demyelinating diseases. Reduced MTR values increase as remyelination and gliosis occur but do not reach the values at baseline before demyelination [11,12,13,14,15].
Numerous studies using DTI and MT have been conducted to confirm the pathology of the cervical spinal cord. In one study, FA and apparent diffusion coefficient (ADC) were analyzed at the maximal compression level and above the compression level (C2/C3) in patients with cervical spondylotic myelopathy. FA at the compression level and ADC above and at compression levels demonstrated significant differences between the patients and the controls. Among the patients, the FA value differed based on the presence of symptoms; however, the ADC value was not related to the symptoms [4]. Another study on compressive cervical myelopathy reported lower FA values in patients than in controls at the proximal spinal cord level (C2/C3). The FA values at the compression level were not analyzed in this study [5]. In a study on MT MRI in patients with traumatic cervical spinal cord injury (SCI), MT saturation of the cervical cord at the C2/C3 level was associated with the pin-prick score [12]. In patients with cervical spondylotic myelopathy, the MTR was lower in the patients than in the controls. In addition, MTR values were inversely correlated with the Nurick score, indicating that a lower MTR value is correlated with more severe the symptoms [16]. However, the level and columns or tracts of the spinal cord were not differentiated in the aforementioned studies. Previous studies have principally analyzed the values at the relatively wide proximal cord or mean values by combining several levels, and the areas of the cord were not divided and analyzed as the entire cord. This is because performing and accurate analysis in a narrow spinal cord area is difficult, and it is even more challenging if the cord is compressed. However, the pathology in the upper and lower spinal cord of the compression lesion is different owing to their longitudinal structure. In addition, findings obtained within the cord vary based on the location, e.g., whether it is in the ascending or descending tract [17,18]. Successful spinal cord tractography depends on many factors, such as image acquisition settings, region of interest design, and distortion correction [19,20]. With the development of imaging techniques and the availability of excellent software, it is now possible to discriminate between detailed structures, such as columns within the spinal cord, which enables tract-specific analysis [15,20,21,22].
We aimed to examine the pathologic changes in cervical spondylotic myelopathy using DTI and MT imaging techniques and confirmed the correlation between the MR parameters (FA, MD, AD, RD, and MTR) and clinical measures. Individual tracts were separately analyzed at lesion levels at the compression as well as above and below the lesion. We hypothesized that tract-specific analysis using DTI and MT imaging can reveal subtle pathologic changes in cervical spondylotic myelopathy, especially in patients with mild cervical spondylotic myelopathy without obvious motor or sensory impairment affecting activities of daily living.

2. Materials and Methods

We recruited patients with cervical spondylotic myelopathy. Patients who complained of neck pain, upper extremity paresthesia or radicular pain, motor weakness, and spasticity and had confirmed radiographic evidence of cord compression with cervical spondylosis were included. Eight patients were enrolled after excluding patients with previous spinal surgery, cervical cord compression of other than degenerative origin, and peripheral nerve lesions that could cause sensory symptoms. The Institutional Review Board of VHS Medical Center approved all the procedures of the study, and written informed consent was obtained from each participant. The Douleur Neuropathique 4 (DN4) questionnaire was employed to assess neuropathic pain. DN4 consists of seven self-report items and three clinical examination questions and is a valid screening tool for patients with central and peripheral neuropathic pain [23]. Scores range from 0 to 10 points, and a score of 4 points is considered positive [24]. The Japanese Orthopaedic Association (JOA) score, the most frequently used scale to evaluate functional status in patients with cervical myelopathy, was also administered. This scale assesses motor dysfunction in the upper and lower extremities; sensory function in the upper extremities, trunk, and lower extremities; and bladder function. Scores range from 0 to 17 points and a score of 13 points is considered as mild dysfunction [25]. The presence or absence of gait abnormality, spinal ataxia, pyramidal signs, and hyperreflexia were evaluated. None of the patient presented gait abnormality, spinal ataxia, or pyramidal signs, such as Hoffman sign, Babinski sign, and ankle clonus. Increased bicep or brachioradialis reflexes were observed in some patients. We evaluated the American Spinal Cord Injury Association sensory and motor scores, lateral pinch, and power grasp. The maximum compression level was determined by calculating the compression ratios at each intervertebral level using conventional MR images. T2 signal changes occurred in some participants. The clinical data are summarized in Table 1.
The control group consisted of 12 healthy individuals; however, they were not age-matched. However, as spondylosis is attributed to degeneration, age-matched healthy individuals may also have mild spondylosis [1,26]. Therefore, we selected young individuals as controls. Table 2 outlines the differences between the patients and the controls.

2.1. MRI Acquisition

2.1.1. Conventional Imaging

Acquisitions were conducted using a 3-Tesla MRI system (Skyra, Siemens, Germany). T1-weighted images were obtained according to the following parameters: axial orientation, one slab of 192 slices; slice thickness, 1.0 mm; field of view (FOV), 320 mm; time repetition (TR), 2000 ms; echo time (TE), 3.72 ms; voxel size, 1.0 × 1.0 × 1.0 mm3; flip angle, 9°; phase oversampling, 0%; slice oversampling, 0%; bandwidth, 150 Hz/pixel; and turbo factor, 192. T2-weighed images were obtained according to the following parameters: axial orientation, one slab of 64 slices; slice thickness, 0.8 mm; FOV, 256 mm; TR, 1500 ms; TE, 119 ms; voxel size, 0.8 × 0.8 × 0.8 mm3; flip angle, 120°; phase oversampling, 80%; slice oversampling, 12.5%; bandwidth, 625 Hz/Pixel; and turbo factor, 88.

2.1.2. DTI

DTI was obtained using single-shot echo-planar imaging with the following parameters: b value, 800 s/mm2; 64 diffusion gradient directions; slice thickness, 5 mm; FOV, 86 mm; TR, 3300 ms; TE, 93 ms; voxel size, 0.9 × 0.9 × 5.0 mm; and bandwidth, 882 Hz/Pixel.

2.1.3. MT Imaging

We acquired T1-weighted three-dimensional gradient echo images with slab-selective excitation, with and without MT saturation pulses using the following parameters: axial orientation, 17 slices; FOV, 160 mm; TR, 600 ms; TE, 4.83 ms; voxel size, 0.6 × 0.6 × 5.0 mm3; flip angle, 25°; phase oversampling, 80%; and bandwidth, 380 Hz/Pix.
Images were processed using open-source software packages, Spinal Cord Toolbox (SCT), and the PAM50 template [27]. The program is available on the website; https://spinalcordtoolbox.com/user_section/installation.html (accessed on 5 August 2021). SCT has been used in cervical cord analysis, even in the presence of pathologies such as traumatic spinal cord injury and myelopathic spinal cord compression, and has proven variability and validity [11,15,17]. DTI parameters, including FA, MD, AD, and RD values, were obtained. We computed the voxel-wise MTR using the following equation: ((S0−SMT)/S0) × 100, where S0 and SMT are T1-weighted images with and without the MT pre-saturation pulse, respectively [15].
MR parameters from three levels were selected as follows: above the lesion level (C2/C3), at the most compressed lesion level, and below the lesion level (C7/T1). The level of the most compressed lesion was determined by calculating the compression ratio [28]. In healthy controls, above the lesion and below the lesion level were analyzed in C2/3 and C7/T1 as in the patient group. As there was no lesion level in the control group, the C5/6 level was selected for comparison with the lesion level of the patients [4,20]. Figure 1 shows the sagittal slices and axial slices of two patients at above the lesion level, the maximum compression level, and below the lesion level. The average of parameters measured from all axial slices contained in each segment were analyzed, and three to four slices were included per segment.
Based on the WM atlas of the human cervical spinal cord, we performed tract-based acquisition of the MR parameters [29]. Among the 15 different WM tracts on each side, three major motor tracts and two major sensory tracts were selected for the analysis. The analyzed motor tracts were as follows: the lateral corticospinal (CS) tract, the rubrospinal (RuS) tract, and the lateral reticulospinal (RS) tract; the sensory tracts were as follows: fasciculus cuneatus (FC) and spinal lemniscus (spinothalamic and spinoreticular tracts) (Figure 2). Alteration of diffusion metrics is significant in the lateral and dorsal columns in the presence of degenerative spinal cord compression [30]. The compressive stress increases in the lateral and dorsal columns when there is acute compression on chronic compressed spinal cord owing to degeneration [31]. The lateral CS tract, RuS tract, lateral RS, spinal lemniscus, and FC tract are located in the lateral and dorsal column, and these are expected to be the most affected.

2.2. Data Processing

Template-based analysis was performed using the SCT, a comprehensive and open-source software for analyzing multiparametric MRI data of the spinal cord [27]. The SCT also offers motion correction of MRI data. The analysis consists of six main steps following motion correction: (1) template registration to an anatomic image, (2) template registration to a WM/gray matter (GM) contrasted image, (3) auto segmentation of the WM and GM, (4) template registration of probabilistic WM/GM atlases to the WM/GM segmentations using a multi-label registration method, (5) template registration to a multiparametric MRI parameter image (FA, MD, AD, RD, and MTR), and (6) atlas-based analysis to quantify the parameters in specific sub-regions of the spinal cord [32]. The registration framework is described in Figure 3.
Although automated designation of the spinal cord area was well performed in the non-compressed normal-appearing spinal cord, an error occurred at the deformed spinal cord owing to compression. The spinal cord area was manually masked if there was an error following a thorough review of the automatically selected cord area.

2.3. Statistical Analysis

Statistical analysis was performed using the SPSS 19.0 Statistical Package for the Social Sciences (IBM Corp., Armonk, NY, USA) by combining the left and right sides as pooled data [5]. We conducted between-group comparisons using the Mann–Whitney U test, chi–square test, and Fisher’s exact test. Statistical significance was set at p < 0.05. Spearman correlation analysis was performed to assess the correlation of each of the three categories of variables (clinical scales, clinical measures, and MR parameters) in patients with cervical spondylotic myelopathy. The false discovery rate was used to perform corrections for multiple comparisons and to better interpret the results.

3. Results

3.1. Group Comparison

FA was lower and MD, AD, and RD were higher in the patient group in the lateral RS tract at each level. Similar changes in the DTI parameters were noted in the RuS tract and spinal lemniscus at the lesion level. Lower FA and higher MD and RD were observed in the RuS tract, and lower FA was observed in the spinal lemniscus. FC also showed group differences in the DTI parametric values; however, the differences were contrary to expectations. FC showed lower MD, AD, and RD and so-called better DTI parametric values in the patient group above the lesion. The MTR was higher in the patient group, contrary to expectations, in the lateral CS tract, RuS tract, lateral RS tract, and spinal lemniscus below the lesion level. Table 3 summarizes the parameters that were different in each tract at each level.

3.2. Correlations between the MR Parameters and Clinical Evaluation Items

All the parameters showing correlations with the clinical evaluation items are summarized in Table 4. In one case, an inverse correlation was observed: the higher the MTR, the worse the symptoms. As the distribution of the correlation seems sporadic, they were grouped as follows and are described in detail: (1) tracts in which the DTI parameters or MTR showed differences between the groups, and correlations with clinical items were also confirmed. (2) DTI parameters and MTR correlated with the clinical items simultaneously.
The RuS tract, lateral RS tract, and spinal lemniscus, which showed group differences in the DTI parameters at the lesion level, also revealed a correlation with the power grasp at the lesion level. A group difference is expressed as a bar graph, and the correlation is presented as a scatter plot in Figure 4. The scatter plots contain 16 data points because the data of the left and right sides were pooled together, as aforementioned. The parameters showing the group differences and the parameters showing the correlation did not necessarily match.
Below the lesion level, the lateral CS tract showed a correlation of the lateral pinch with the DTI parameters and MTR simultaneously. The correlation is presented as a scatter plot in Figure 5. Among the DTI parameters, FA and RD were associated with the MTR.

4. Discussion

4.1. Group Comparison

It was hypothesized that symptoms of cervical spondylotic myelopathy result from pathologic changes in the specific tract rather than from the compression itself [33]. We intended to detect the changes in the spinal cord tracts in patients exposed to chronic compression due to spondylosis compared to healthy controls. To adjust for the age difference between the groups, a regression analysis with age as a covariate was conducted but was not appropriate because of multicollinearity. To complement this, regression analysis was conducted to control for within-group age differences. Tracts showing expected results were the lateral RS tract at all levels, the RuS tract, and spinal lemniscus at the lesion level. A slow neuroinflammatory response of the spinal cord under chronic compression has been reported in previous animal and human studies [34,35,36]. This slow neuroinflammatory process results in chronic demyelination and causes activation of the microglia and the release of toxic stimuli from cells, such as astrocytes and neurons [33]. Our study demonstrated consistent results through changes in the DTI values of the RuS tract, lateral RS tract, and spinal lemniscus. In particular, the lateral RS tract showed consistent results at all levels, with the most pronounced group difference. Unlike the lateral RS tract, the RuS tract and spinal lemniscus showed differences only at the lesion level, which could occur because the lesion level was most directly affected by compression. The majority of the literature on the function of specific tracts was performed in primates rather than in humans. The RS tract plays an important role in the coordination of activity primarily performed by the corticospinal tract, such as reaching movements in the upper extremities and possibly finger movements, and has been shown to be associated with locomotion and postural adjustment [12,37,38]. The RuS tract is directly connected with motor neurons of the distal upper limb muscles and is associated with dynamic movement control in primates [39]. Both the RS and RuS tracts are also known to play an important role in recovery following corticospinal cord lesions in primates and humans, and increased activity of these tracts is a favorable prognostic factor for restoration of hand function [40,41,42]. In our study, lateral CS did not show any difference in DTI parameters between groups; therefore, the RuS and RS tracts did not show a compensatory increase in activity either. This difference is thought to be attributed to the fact that the patients in our study had relatively milder pathology and symptoms compared to previous studies performed with primates with induced SCI and human SCI showing clearly impaired hand function. None of the patients in our study had significant motor weakness or limitations in independent daily activities. The spinal lemniscus refers to the anterior and lateral spinothalamictracts. The spinothalamic tract is a major pathway of pain transmission, and neuropathic pain is worse if the function of the spinothalamic tract decreases in SCI [43]. Degeneration of the spinothalamic tract was confirmed by MRI, and the degree of degeneration correlated with the degree of impairment in traumatic SCI [44]. DTI parameters reflect the pathology of the spinothalamic tract. Interestingly, Hatem et al. reported that the FA value of the anterior cord area, where the spinothalamic tract is located, correlated better with somatosensory dysfunction than the FA value of the whole cord or posterior column [45]. Our study also showed lower FA in the spinal lemniscus in patients, which is consistent with the results of previous studies. Considering the DTI parameters, the RuS and lateral RS tracts, which are motor tracts, showed differences in most DTI parameters, whereas spinal lemniscus, which is the sensory tract, only showed differences in the FA values. Our finding confirms that FA is the most sensitive parameter among the DTI parameters. Judging from their proximity to each other, another possible explanation is that the motor tract is more vulnerable to pressure than the sensory tract.
FC in the dorsal column showed group differences in DTI parameters above the level of the lesion. Animal studies have demonstrated that the damaged tract differs based on the position, which is above and below the compression lesion. A previous study reported changes above and below the compression level with DTI and histopathology by artificially inducing SCI with the application of compression to the spinal cord in dogs. The study pathologically confirmed spinal cord degeneration in the dorsal funiculus above the lesion and in the ventral funiculus below the lesion [18]. In a study on cats, changes in FA, generalized fractional anisotropy (GFA), and RD were predominant; however, they were also detected rostral to the injury on the dorsal aspect and caudal to the injury on the ventrolateral aspect [46]. Degeneration does not necessarily occur in the anterograde direction only, and retrograde degeneration also occurs; however, antegrade Wallerian degeneration could explain the location of the pathologic change based on the positional relation with the lesion [15]. FC is an ascending tract that transmits proprioceptive and cutaneous feedback signals to the upper extremities; therefore, the group difference above the level of the lesion, which was confirmed in our study, was consistent with previous studies [47,48]. However, we were unable to explain why the DTI values of some tracts were better in patients than in the controls because this phenomenon has not been previously reported. The mild symptoms of the patients enrolled in this study compared to other studies could have contributed to this result. The possibility of repair or regeneration of the spinal cord tracts was assessed; however, it seems unlikely. A literature review on the mechanisms of central nervous system axonal regeneration and remyelination by Uyeda et al. identified that astrocytes, oligodendrocytes, and glial cells play critical roles in the repair process [49]. With degeneration, myelin content decreases and astrocytes become hypertrophied, and these cellular changes are reflected in the DTI parameters as a decrease in FA and an increase in MD, AD, and RD. With regeneration, the myelin content and hypertrophied astrocytes recover, and these changes are reflected in the DTI parameters as an increase in FA, as well as a decrease in MD, AD, and RD. However, recovery of the DTI parameter values did not reach or exceed the baseline level [50,51]. Better DTI values in the patient group in our study cannot be explained as regeneration because they cannot exceed the baseline value, even when regeneration occurs. A compensation mechanism may exist in FC; however, it is difficult to explain this unexpected finding because the number of patients was small and the patients had only mild symptoms.
The MTR, lateral CS tract, RuS tract, lateral RS tract, and spinal lemniscus tracts showed group differences in the MTR below the lesion level, and a higher MTR was observed in the patient group, contrary to expectations. Our findings are not consistent with those of Cloney et al., who reported that the MTR of the anterior cord region was associated with the clinical outcomes [52]. A post-mortem study of humans demonstrated that the MTR value has a strong correlation with the myelin content [53]. MTR indicates demyelination and reduction sensitively, even before demyelination becomes visible on T2-weighted images and increases again with remyelination; however, the value does not reach the baseline, similar to the DTI parameters [54]. One possible explanation for this is that, although myelin layers were formed, the lipid and/or protein content of the newly repaired myelin is biochemically different from mature, compact myelin, affecting the magnetization transfer ratio [55]. Therefore, it is difficult to determine the extent of remyelination using only the MTR values [56]. Moreover, despite its sensitivity to the tissue content, the MTR is not specific to the underlying pathology that is affected by free water MR parameters, which are modulated by other factors, such as inflammation [57]. The level or the location of the tracts with higher MTR in our study is inconsistent, and considering the characteristics of MTR, it is difficult to ascertain the significance of this result.
As the control group was not age-matched, it was necessary to clarify the effect of aging. Previous studies have demonstrated that FA and MTR increase with aging; in contrast, MD and RD increases with aging. The effect of age on AD value is not consistent across reports [58,59,60]. Age itself may have influenced the changes in MR metrics in the patient group. However, as the differences in MR in our study were correlated with the lesion level, it is clear that this change was attributed to spondylotic compression.

4.2. Correlations between the MR Parameters and Clinical Items

Although the analysis was performed using advanced analysis technology, it is difficult to ascertain the significance of every single result demonstrating the correlation with the MR parameters. Moreover, previous studies have claimed that caution is needed in the interpretation of the MR parameters. Several physical parameters can influence the DTI parameters, including myelination, axonal density, axonal diameter, and orientation of the fiber bundles [61,62]. RD and AD can cause a deceptive change in the voxels characterized by crossing fibers [9]. The direction of diffusivity is not always preserved in pathological tissue and is not always aligned with the underlying expected tissue architecture [63]. It is necessary to pay attention to the interpretation of the DTI parameters; however, it can be considered a significant finding when there is meaningful consistency. Therefore, in our study, the following cases were judged to be significant, and the results were examined in detail: (1) the tracts in which group differences were confirmed and correlations with clinical items were also confirmed within the patient groups; (2) the tracts in which the DTI parameter and the MTR showed a correlation in the same direction as the clinical items, although there was no group difference. Evidence exists that combining DTI and MT is more specific and reproducible. Reich et al. reported that combining measures of MT and diffusion-weighted MRI is a means to becoming more specific to WM pathology [64]. Smith et al. also reported that MT and diffusion-weighted MRI together provide high reproducibility in the human cervical cord at 3T as well as a robust assessment of the WM pathology [65]. Cohen-Adad et al. identified that combining high angular diffusion imaging with MT is a promising approach to gaining specificity in characterizing spinal cord pathways in traumatic injury [15]. Figure 6 presents the schematic interpretation process of this study concerning the MR parameters.
The lateral RS tract, RuS tract, and spinal lemniscus, which showed worse DTI parameters in patients than in controls at the lesion level, showed a correlation of the DTI parameters with the power grasp at the lesion level. It is known that the RS and RuS tracts are involved in hand motor function. The RS tract contributes to hand motor tasks involving gross finger manipulations, but to a lesser extent, to more dexterous finger movements in humans with cervical SCI [32]. Recovery of gross hand movement is impaired when RuS tract involvement is present [34]. Our results corroborate previous findings that the lateral RS and RuS tracts contribute to gross hand function. The DTI parameters of the lateral RS and RuS tract were well correlated with the power grasp, explaining the patient’s subtle symptoms well. It is difficult to specify tract-specific pathology via conventional MRI; therefore, tract-specific analysis provides a good chance to identify the pathology in detail and understand the associated symptoms. Spinal lemniscus was associated with lower FA values in patients than controls and FA values were correlated with power grasp at the lesion level. There are also some differences in the results of the spinal lemniscus compared with the lateral RS and RuS tracts. In the group comparison, there were differences in three or more DTI parameters in the lateral RS and RuS tracts; however, only one DTI parameter, FA, was different in the spinal lemniscus. While the parameters showing group differences and those showing correlation overlapped in the lateral RS and RuS tracts, parameters showing group differences and parameters showing correlation did not match in the spinal lemniscus. Even if correlation is confirmed in the tracts showing differences between the groups, it is difficult to determine the unity between the DTI parameters.
The DTI parameters and the MTR showed a correlation in the same direction with clinical features in the lateral CS tract. DTI parameters and MTR showed a positive correlation with lateral pinch below the lesion level in the lateral CS tract. Among the DTI parameters, FA showed correlations; thus, we can see that FA is closely associated with MTR. FA is the most commonly used DTI parameter because it sensitively measures microstructural integrity, although it is not specific to the underlying pathology. In our study, FA showed the most significant DTI results, confirming that FA is a sensitive parameter. MTR increases when there are more proteins and lipids in the myelinated sheath [57]. This is the only significant finding that was demonstrated in the lateral CS tract, the major motor pathway, in our study. There was no difference between the groups in the lateral CS tract, and we assumed that this was attributed to the fact that the patients included in our study only had mild symptoms. However, as the MTR values of the lateral CS tract below the lesion level correlated with the lateral pinch, a change in the myelin state of the lateral CS tract was observed in the patients.
One of the greatest challenge when acquiring MR images in the spinal cord is the inhomogeneous magnetic field in the region [65,66,67], and there are certain possible errors in DTI-based tractography due to gradient inhomogeneity in the spatial context-B matrix spatial distribution. The SCT we used in our study provides few correction techniques. Thus, we applied the optimal MRI acquisition protocol for the processing and analysis using SCT. Correction was performed using slice-wise translation in the axial plane, which in the SCT can be performed using regularization along the z direction. In addition to this feature, correction in SCT includes the methods of Xuet et al. for grouping successive volumes, thereby improving the robustness of registration in high b-value diffusion MRI data. This approach is particularly successful for correcting slow drifts in long acquisitions [67]. Nevertheless, the possibility that diffusion tensor values were affected by using a single device cannot be completely excluded. Even if the average eigenvalue shifted, the difference between the patient and the control groups can be confirmed without looking at the absolute value [68].
The small sample size and lack of age-matched controls are certain limitations of this study. To increase the reliability of our results, future studies with a larger number of patients with age-matched controls should be conducted.

5. Conclusions

Through tract-specific analysis, the difference between the control group and patients, and the correlation between individual tracts and clinical symptoms, could be identified in cervical spondylotic myelopathy. Tract-specific analysis and interpretation of MR parameters should be considered from various perspectives. If the DTI parameter and MTR show unity, it would be reliable, and among the DTI parameters, FA is the most related to MTR. If only the DTI parameters are reviewed without considering the MTR, uniformity should exist among the DTI parameters. The positional relationship between the tracts within an axial plane must be considered when determining the MR parameters of each tract. Considering all the aforementioned points, we obtained the following: RuS and lateral RS showed worse DTI parameters at the lesion level in patients than in the controls, and worse DTI parameters in those tracts were correlated with weaker power grasp at the lesion level. Lateral CS showed a correlation between higher RD and lower FA and MTR, with a weaker lateral pinch below the lesion level.

Author Contributions

Conceptualization, W.-K.Y.; Data curation, D.-H.K.; Methodology, W.-T.K.; Software, D.-H.K.; Supervision, S.-W.K.; Writing–original draft, H.-E.Y. All authors have read and agreed to the published version of the manuscript.

Funding

The study was supported by a VHS Medical Center Research Grant, Republic of Korea (Grant No: VSHMC 18017).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board of VHS Medical Center (IRB File No. 2018-02-010).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Data available on request to restrictions. The data presented in this study are available on request from the corresponding author. The data are not pubulicaly available due to participant’s confidentiality.

Acknowledgments

We would like to acknowledge and thank VHS Medical Center Research Grant.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Theodore, N. Degenerative cervical spondylosis. N. Engl. J. Med. 2020, 383, 159–168. [Google Scholar] [CrossRef] [PubMed]
  2. McCormick, J.R.; Sama, A.J.; Schiller, N.C.; Butler, A.J.; Donnally, C.J., 3rd. Cervical spondylotic myelopathy: A guideguide to diagnosis and management. J. Am. Board Fam. Med. 2020, 33, 303–313. [Google Scholar] [CrossRef] [PubMed]
  3. Hashizume, Y.; Iijima, S.; Kishimoto, H.; Yanagi, T. Pathology of spinal cord lesions caused by ossification of the posterior longitudinal ligament. Acta Neuropathol. 1984, 63, 123–130. [Google Scholar] [CrossRef] [PubMed]
  4. Kerkovský, M.; Bednarík, J.; Dušek, L.; Sprláková-Puková, A.; Urbánek, I.; Mechl, M.; Válek, V.; Kadanka, Z. Magnetic resonance diffusion tensor imaging in patients with cervical spondylotic spinal cord compression: Correlations between clinical and electrophysiological findings. Spine 2012, 37, 48–56. [Google Scholar] [CrossRef]
  5. Gohmann, R.F.; Blume, C.; Zvyagintsev, M.; Mainz, V.; Clusmann, H.; Wiesmann, M.; Brockmann, M.A.; Mueller, C.A. Cervical spondylotic myelopathy: Changes of fractional anisotropy in the spinal cord and magnetic resonance spectroscopy of the primary motor cortex in relation to clinical symptoms and their duration. Eur. J. Radiol. 2019, 116, 55–60. [Google Scholar] [CrossRef]
  6. Costanzo, R.; Brunasso, L.; Paolini, F.; Benigno, U.E.; Porzio, M.; Giammalva, G.R.; Gerardi, R.M.; Umana, G.E.; di Bonaventura, R.; Sturiale, C.L.; et al. Spinal tractography as a potential prognostic tool in spinal cord injury. A systematic review. World Neurosurg. 2022, 29, 25–32. [Google Scholar] [CrossRef]
  7. Della Della Nave, R.D.; Ginestroni, A.; Diciotti, S.; Salvatore, E.; Soricelli, A.; Mascalchi, M. Axial diffusivity is increased in the degenerating superior cerebellar peduncles of Friedreich’s ataxia. Neuroradiology 2011, 53, 367–372. [Google Scholar] [CrossRef]
  8. Aung, W.Y.; Mar, S.; Benzinger, T.L. Diffusion tensor MRI as a biomarker in axonal and myelin damage. Imaging Med. 2013, 5, 427–440. [Google Scholar] [CrossRef]
  9. Cercignani, M.; Wheeler-Kingshott, C.A.M. About “axial” and “radial” diffusivities. Magn. Reason. Med. 2009, 61, 1255–1260. [Google Scholar] [CrossRef]
  10. Murgoci, A.N.; Baciak, L.; Cubinkova, V.; Smolek, T.; Tvrdik, T.; Juranek, I.; Kafka, J.; Cizkova, D. Diffusion tensor imaging: Tool for tracking injured spinal cord fibres in rat. Neurochem. Res. 2020, 45, 180–187. [Google Scholar] [CrossRef]
  11. Seif, M.; GandiniGandini Wheeler-Kingshott, C.A.G.; Cohen-Adad, J.; Flanders, A.E.; Freund, P. Guidelines for the conduct of clinical trials in spinal cord injury: Neuroimaging biomarkers. Spinal Cord 2019, 57, 717–728. [Google Scholar] [CrossRef] [PubMed]
  12. Grabher, P.; Callaghan, M.F.; Ashburner, J.; Weiskopf, N.; Thompson, A.J.; Curt, A.; Freund, P. Tracking sensory system atrophy and outcome prediction in spinal cord injury. Ann. Neurol. 2015, 78, 751–761. [Google Scholar] [CrossRef] [PubMed]
  13. Nouri, A.; Martin, A.R.; Mikulis, D.; Fehlings, M.G. Magnetic resonance imaging assessment of degenerative cervical myelopathy: A review of structural changes and measurement techniques. Neurosurg. Focus 2016, 40, E5. [Google Scholar] [CrossRef]
  14. Tomiak, M.M.; Rosenblum, J.D.; Prager, J.M.; Metz, C.E. Magnetization transfer: A potential method to determine the age of multiple sclerosis lesions. AJNR Am. J. Neuroradiol. 1994, 15, 1569–1574. [Google Scholar]
  15. Cohen-Adad, J.; El Mendili, M.M.; Lehéricy, S.; Pradat, P.F.; Blancho, S.; Rossignol, S.; Benali, H. Demyelination and degeneration in the injured human spinal cord detected with diffusion and magnetization transfer MRI. Neuroimage 2011, 55, 1024–1033. [Google Scholar] [CrossRef]
  16. Suleiman, L.I.; Weber, K.A.; Rosenthal, B.D.; Bhatt, S.A.; Savage, J.W.; Hsu, W.K.; Patel, A.A.; Parrish, T.B. High-resolution magnetization transfer MRI in patients with cervical spondylotic myelopathy. J. Clin. Neurosci. 2018, 51, 57–61. [Google Scholar] [CrossRef]
  17. Cohen-Adad, J.; Leblond, H.; Delivet-Mongrain, H.; Martinez, M.; Benali, H.; Rossignol, S. Wallerian degeneration after spinal cord lesions in cats detected with diffusion tensor imaging. Neuroimage 2011, 57, 1068–1076. [Google Scholar] [CrossRef]
  18. Yoon, H.; Kim, J.; Moon, W.J.; Nahm, S.S.; Zhao, J.; Kim, H.M.; Eom, K. Characterization of chronic axonal degeneration using diffusion tensor imaging in canine spinal cord injury: A quantitative analysis of Diffusion Tensor Imaging parameters according to histopathological differences. J. Neurotrauma 2017, 34, 3041–3050. [Google Scholar] [CrossRef]
  19. Dauleac, C.; Frindel, C.; Mertens, P.; Jacquesson, T.; Cotton, F. Overcoming challenges of the human spinal cord tractography for routine clinical use: A review. Neuroradiology 2020, 62, 1079–1094. [Google Scholar] [CrossRef]
  20. Dauleac, C.; Bannier, E.; Cotton, F.; Frindel, C. Effect of distortion corrections on the tractography quality in spinal cord diffusion-weighted imaging. Magn. Reason. Med. 2021, 85, 3241–3255. [Google Scholar] [CrossRef]
  21. Cohen-Adad, J. Microstructural imaging in the spinal cord and validation strategies. Neuroimage 2018, 182, 169–183. [Google Scholar] [CrossRef] [PubMed]
  22. Yang, Y.M.; Yoo, W.K.; Bashir, S.; Oh, J.K.; Kwak, Y.H.; Kim, S.W. Spinal cord changes after laminoplasty in cervical compressive myelopathy: A diffusion tensor imaging study. Front. Neurol. 2018, 9, 696. [Google Scholar] [CrossRef] [PubMed]
  23. VanDenKerkhof, E.G.; Stitt, L.; Clark, A.J.; Gordon, A.; Lynch, M.; Morley-Forster, P.K.; Nathan, H.J.; Smyth, C.; Toth, C.; Ware, M.A.; et al. Sensitivity of the DN4 in screening for neuropathic pain syndromes. Clin. J. Pain 2018, 34, 30–36. [Google Scholar] [CrossRef]
  24. Bouhassira, D.; Attal, N.; Alchaar, H.; Boureau, F.; Brochet, B.; Bruxelle, J.; Cunin, G.; Fermanian, J.; Ginies, P.; Grun-Overdyking, A.; et al. Comparison of pain syndromes associated with nervous or somatic lesions and development of a new neuropathic pain diagnostic questionnaire (DN4). Pain 2005, 114, 29–36. [Google Scholar] [CrossRef] [PubMed]
  25. Kato, S.; Oshima, Y.; Oka, H.; Chikuda, H.; Takeshita, Y.; Miyoshi, K.; Kawamura, N.; Masuda, K.; Kunogi, J.; Okazaki, R.; et al. Comparison of the Japanese orthopaedic association (JOA) score and modified JOA (mJOA) score for the assessment of cervical myelopathy: A multicenter observational study. PLoS ONE 2015, 10, e0123022. [Google Scholar] [CrossRef]
  26. Fehlings, M.G.; Skaf, G. A review of the pathophysiology of cervical spondylotic myelopathy with insights for potential novel mechanisms drawn from traumatic spinal cord injury. Spine 1998, 23, 2730–2737. [Google Scholar] [CrossRef]
  27. De De Leener, B.D.; Fonov, V.S.; Collins, D.L.; Callot, V.; Stikov, N.; Cohen-Adad, J. PAM50: Unbiased multimodal template of the brainstem and spinal cord aligned with the ICBM152 space. Neuroimage 2018, 165, 170–179. [Google Scholar] [CrossRef]
  28. Fujiwara, K.; Yonenobu, K.; Ebara, S.; Yamashita, K.; Ono, K. The prognosis of surgery for cervical compression myelopathy. An analysis of the factors involved. J. Bone Jt. Surg. Br. 1989, 71-B, 393–398. [Google Scholar] [CrossRef]
  29. Lévy, S.; Benhamou, M.; Naaman, C.; Rainville, P.; Callot, V.; Cohen-Adad, J. White matter atlas of the human spinal cord with estimation of partial volume effect. Neuroimage 2015, 119, 262–271. [Google Scholar] [CrossRef]
  30. Valošek, J.; Labounek, R.; Horák, T.; Horáková, M.; Bednařík, P.; Keřkovský, M.; Kočica, J.; Rohan, T.; Lenglet, C.; Cohen-Adad, J.; et al. Diffusion magnetic resonance imaging reveals tract-specific microstructural correlates of electrophysiological impairments in non-myelopathic and myelopathic spinal cord compression. Eur. J. Neurol. 2021, 28, 3784–3797. [Google Scholar] [CrossRef]
  31. Ichihara, K.; Taguchi, T.; Sakuramoto, I.; Kawano, S.; Kawai, S. Mechanism of the spinal cord injury and the cervical spondylotic myelopathy: New approach based on the mechanical features of the spinal cord white and gray matter. J. Neurosurg. 2003, 99, 278–285. [Google Scholar] [CrossRef]
  32. Dupont, S.M.; De Leener, B.D.; Taso, M.; Le Troter, A.; Nadeau, S.; Stikov, N.; Callot, V.; Cohen-Adad, J. Fully-integrated framework for the segmentation and registration of the spinal cord white and gray matter. Neuroimage 2017, 150, 358–372. [Google Scholar] [CrossRef] [PubMed]
  33. Craciunas, S.C.; Gorgan, M.R.; Ianosi, B.; Lee, P.; Burris, J.; Cirstea, C.M. Remote motor system metabolic profile and surgery outcome in cervical spondylotic myelopathy. J. Neurosurg. Spine 2017, 26, 668–678. [Google Scholar] [CrossRef] [PubMed]
  34. Beattie, M.S.; Harrington, A.W.; Lee, R.; Kim, J.Y.; Boyce, S.L.; Longo, F.M.; Bresnahan, J.C.; Hempstead, B.L.; Yoon, S.O. ProNGF induced p75-mediated death of oligodendrocytes following spinal cord injury. Neuron 2002, 36, 375–386. [Google Scholar] [CrossRef]
  35. Beattie, M.S.; Manley, G.T. Tight squeeze, slow burn: Inflammation and the aetiology of cervical myelopathy. Brain 2011, 134, 1259–1261. [Google Scholar] [CrossRef] [PubMed]
  36. Yu, W.R.; Liu, T.; Kiehl, T.R.; Fehlings, M.G. Human neuropathological and animal model evidence supporting a role for Fas-mediated apoptosis and inflammation in cervical spondylotic myelopathy. Brain 2011, 134, 1277–1292. [Google Scholar] [CrossRef]
  37. Erschbamer, M.; Öberg, J.; Westman, E.; Sitnikov, R.; Olson, L.; Spenger, C. 1H-MRS in spinal cord injury: Acute and chronic metabolite alterations in rat brain and lumbar spinal cord. Eur. J. Neurosci. 2011, 33, 678–688. [Google Scholar] [CrossRef]
  38. Fehlings, M.G.; Wilson, J.R.; Kopjar, B.; Yoon, S.T.; Arnold, P.M.; Massicotte, E.M.; Vaccaro, A.R.; Brodke, D.S.; Shaffrey, C.I.; Smith, J.S.; et al. Efficacy and safety of surgical decompression in patients with cervical spondylotic myelopathy: Results of the AOSpine North America prospective multi-center study. J. Bone Jt. Surg. Am. 2013, 95, 1651–1658. [Google Scholar] [CrossRef]
  39. Isa, T.; Kinoshita, M.; Nishimura, Y. Role of Direct vs. indirect Pathways from the Motor Cortex to Spinal Motoneurons in the Control of Hand Dexterity. Front. Neurol. 2013, 4, 191. [Google Scholar] [CrossRef]
  40. Isa, T.; Mitsuhashi, M.; Yamaguchi, R. Alternative routes for recovery of hand functions after corticospinal tract injury in primates and rodents. Curr. Opin. Neurol. 2019, 32, 836–843. [Google Scholar] [CrossRef]
  41. Baker, S.N.; Perez, M.A. Reticulospinal contributions to gross hand function after human spinal cord injury. J. Neurosci. 2017, 37, 9778–9784. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Cowley, K.C.; Lane, M.A.; Meehan, C.F.; Rank, M.M.; Stecina, K. Editorial: Propriospinal neurons: Essential elements in locomotion, autonomic function and plasticity after spinal cord injury and disease. Front. Cell Neurosci. 2021, 15, 695424. [Google Scholar] [CrossRef] [PubMed]
  43. Cruz-Almeida, Y.; Felix, E.R.; Martinez-Arizala, A.; Widerström-Noga, E.G. Decreased spinothalamic and dorsal column medial lemniscus-mediated function is associated with neuropathic pain after spinal cord injury. J. Neurotrauma 2012, 29, 2706–2715. [Google Scholar] [CrossRef]
  44. Fischer, T.; Stern, C.; Freund, P.; Schubert, M.; Sutter, R. Wallerian degeneration in cervical spinal cord tracts is commonly seen in routine T2-weighted MRI after traumatic spinal cord injury and is associated with impairment in a retrospective study. Eur. Radiol. 2021, 31, 2923–2932. [Google Scholar] [CrossRef] [PubMed]
  45. Hatem, S.M.; Attal, N.; Ducreux, D.; Gautron, M.; Parker, F.; Plaghki, L.; Bouhassira, D. Assessment of spinal somatosensory systems with diffusion tensor imaging in syringomyelia. J. Neurol. Neurosurg. Psychiatry 2009, 80, 1350–1356. [Google Scholar] [CrossRef] [PubMed]
  46. Perlson, E.; Maday, S.; Fu, M.M.; Moughamian, A.J.; Holzbaur, E.L. Retrograde axonal transport: Pathways to cell death? Trends Neurosci. 2010, 33, 335–344. [Google Scholar] [CrossRef]
  47. Akutagawa, T.; Tani, T.; Kida, K.; Tadokoro, N.; Enoki, H.; Nagano, Y.; Ikeuchi, M. A new method for characterizing hand dysfunction in cervical spondylotic myelopathy: A preliminary study. Spinal Cord 2016, 54, 221–225. [Google Scholar] [CrossRef]
  48. Lee, J.K.; Liu, D.; Jiang, D.; Kulikowicz, E.; Tekes, A.; Liu, P.; Qin, Q.; Koehler, R.C.; Aggarwal, M.; Zhang, J.; et al. Fractional anisotropy from diffusion tensor imaging correlates with acute astrocyte and myelin swelling in neonatal swine models of excitotoxic and hypoxic-ischemic brain injury. J. Comp. Neurol. 2021, 529, 2750–2770. [Google Scholar] [CrossRef]
  49. Yano, R.; Hata, J.; Abe, Y.; Seki, F.; Yoshida, K.; Komaki, Y.; Okano, H.; Tanaka, K.F. Quantitative temporal changes in DTI values coupled with histological properties in cuprizone-induced demyelination and remyelination. Neurochem. Int. 2018, 119, 151–158. [Google Scholar] [CrossRef]
  50. Cloney, M.B.; Smith, Z.A.; Weber, K.A.; Parrish, T.B. Quantitative magnetization transfer MRI measurements of the anterior spinal cord region are associated with clinical outcomes in cervical spondylotic myelopathy. Spine 2018, 43, 675–680. [Google Scholar] [CrossRef]
  51. Mottershead, J.P.; Schmierer, K.; Clemence, M.; Thornton, J.S.; Scaravilli, F.; Barker, G.J.; Tofts, P.S.; Newcombe, J.; Cuzner, M.L.; Ordidge, R.J.; et al. High field MRI correlates of myelin content and axonal density in multiple sclerosis—A post-mortem study of the spinal cord. J. Neurol. 2003, 250, 1293–1301. [Google Scholar] [CrossRef] [PubMed]
  52. Chen, J.T.; Collins, D.L.; Atkins, H.L.; Freedman, M.S.; Arnold, D.L.; Canadian MS/BMT Study Group. Magnetization transfer ratio evolution with demyelination and remyelination in multiple sclerosis lesions. Ann. Neurol. 2008, 63, 254–262. [Google Scholar] [CrossRef]
  53. McCreary, C.R.; Bjarnason, T.A.; Skihar, V.; Mitchell, J.R.; Yong, V.W.; Dunn, J.F. Multiexponential T2 and magnetization transfer MRI of demyelination and remyelination in murine spinal cord. Neuroimage 2009, 45, 1173–1182. [Google Scholar] [CrossRef] [PubMed]
  54. Barkhof, F.; Brück, W.; De Groot, C.J.; Bergers, E.; Hulshof, S.; Geurts, J.; Polman, C.H.; van der Valk, P. Remyelinated lesions in multiple sclerosis: MagneticMagnetic resonance image appearance. Arch. Neurol. 2003, 60, 1073–1081. [Google Scholar] [CrossRef] [PubMed]
  55. Alexander, A.L.; Hurley, S.A.; Samsonov, A.A.; Adluru, N.; Hosseinbor, A.P.; Mossahebi, P.; Tromp Do, P.M.; Zakszewski, E.; Field, A.S. Characterization of cerebral white matter properties using quantitative magnetic resonance imaging stains. Brain Connect. 2011, 1, 423–446. [Google Scholar] [CrossRef]
  56. Taso, M.; Girard, O.M.; Duhamel, G.; Le Troter, A.; Feiweier, T.; Guye, M.; Ranjeva, J.P.; Callot, V. Tract-specific and age-related variations of the spinal cord microstructure: A multi-parametric MRI study using diffusion tensor imaging (DTI) and inhomogeneous magnetization transfer (ihMT). NMR Biomed. 2016, 29, 817–832. [Google Scholar] [CrossRef]
  57. Sullivan, E.V.; Rohlfing, T.; Pfefferbaum, A. Quantitative fiber tracking of lateral and interhemispheric white matter systems in normal aging: Relations to timed performance. Neurobiol. Aging 2010, 31, 464–481. [Google Scholar] [CrossRef]
  58. Van Hecke, W.V.; Leemans, A.; Sijbers, J.; Vandervliet, E.; Goethem, J.V.; Parizel, P.M. A tracking-based diffusion tensor imaging segmentation method for the detection of diffusion-related changes of the cervical spinal cord with aging. J. Magn. Reason. Imaging 2008, 27, 978–991. [Google Scholar] [CrossRef]
  59. Beaulieu, C. The basis of anisotropic water diffusion in the nervous system—A technical review. NMR Biomed. 2002, 15, 435–455. [Google Scholar] [CrossRef]
  60. Sen, P.N.; Basser, P.J. A model for diffusion in white matter in the brain. Biophys. J. 2005, 89, 2927–2938. [Google Scholar] [CrossRef]
  61. Field, A.S.; Alexander, A.L.; Wu, Y.C.; Hasan, K.M.; Witwer, B.; Badie, B. Diffusion tensor eigenvector directional color imaging patterns in the evaluation of cerebral white matter tracts altered by tumor. J. Magn. Reason. Imaging 2004, 20, 555–562. [Google Scholar] [CrossRef] [PubMed]
  62. Reich, D.S.; Smith, S.A.; Zackowski, K.M.; Gordon-Lipkin, E.M.; Jones, C.K.; Farrell, J.A.; Mori, S.; van Zijl, P.C.; Calabresi, P.A. Multiparametric magnetic resonance imaging analysis of the corticospinal tract in multiple sclerosis. Neuroimage 2007, 38, 271–279. [Google Scholar] [CrossRef] [PubMed]
  63. Smith, S.A.; Jones, C.K.; Gifford, A.; Belegu, V.; Chodkowski, B.; Farrell, J.A.; Landman, B.A.; Reich, D.S.; Calabresi, P.A.; McDonald, J.W.; et al. Reproducibility of tract-specific magnetization transfer and diffusion tensor imaging in the cervical spinal cord at 3 tesla. NMR Biomed. 2010, 23, 207–217. [Google Scholar] [CrossRef] [PubMed]
  64. Song, S.K.; Yoshino, J.; Le, T.Q.; Lin, S.J.; Sun, S.W.; Cross, A.H.; Armstrong, R.C. Demyelination increases radial diffusivity in corpus callosum of mouse brain. Neuroimage 2005, 26, 132–140. [Google Scholar] [CrossRef] [PubMed]
  65. Borkowski, K.; Krzyżak, A.T. The generalized Stejskal-Tanner equation for non-uniform magnetic field gradients. J. Magn. Reson. Reson. 2018, 296, 23–28. [Google Scholar] [CrossRef]
  66. Stroman, P.W.; Wheeler-Kingshott, C.; Bacon, M.; Schwab, J.M.; Bosma, R.; Brooks, J.; Cadotte, D.; Carlstedt, T.; Ciccarelli, O.; Cohen-Adad, J.; et al. The current state-of-the-art of spinal cord imaging: Methods. Neuroimage 2014, 84, 1070–1081. [Google Scholar] [CrossRef]
  67. Leener, B.D.; Lévy, S.; Dupont, S.M.; Fonov, V.S.; Stikov, N.; Collins, D.L.; Callot, V.; Cohen-Adad, J. SCT: Spinal Cord Toolbox, an open-source software for processing spinal cord MRI data. Neuroimage 2017, 145, 24–43. [Google Scholar] [CrossRef]
  68. Borkowski, K.; Krzyżak, A.T. Assessment of the systematic errors caused by diffusion gradient inhomogeneity in DTI-computer simulations. NMR Biomed. 2019, 32, e4130. [Google Scholar] [CrossRef]
Figure 1. Examples of sagittal and axial slices at each level in two patients. The maximum compression level is at the C6/7 (A) and at the C4/5 (B). Panel (C) shows vertebral labeling, cord area at each level, and white matter/gray matter segmentation of patient A. More details concerning data processing are presented below.
Figure 1. Examples of sagittal and axial slices at each level in two patients. The maximum compression level is at the C6/7 (A) and at the C4/5 (B). Panel (C) shows vertebral labeling, cord area at each level, and white matter/gray matter segmentation of patient A. More details concerning data processing are presented below.
Diagnostics 12 02090 g001
Figure 2. Spinal cord white matter tracts. The colored atlas of the spinal cord white matter has been adapted from Lévy et al., 2015 [29].
Figure 2. Spinal cord white matter tracts. The colored atlas of the spinal cord white matter has been adapted from Lévy et al., 2015 [29].
Diagnostics 12 02090 g002
Figure 3. Registration framework for the extraction of metrics using the template. Template-based elements are shown in red. The figure has been adapted from Dupont et al., 2017 [32].
Figure 3. Registration framework for the extraction of metrics using the template. Template-based elements are shown in red. The figure has been adapted from Dupont et al., 2017 [32].
Diagnostics 12 02090 g003
Figure 4. Group difference and correlation between MR parameters and motor items at rubrospinal tract (a), lateral reticulospinal tract (b), and spinal lemniscus (c). * statistically significant.
Figure 4. Group difference and correlation between MR parameters and motor items at rubrospinal tract (a), lateral reticulospinal tract (b), and spinal lemniscus (c). * statistically significant.
Diagnostics 12 02090 g004aDiagnostics 12 02090 g004b
Figure 5. Correlation between MR parameters and lateral pinch at lateral corticospinal tract.
Figure 5. Correlation between MR parameters and lateral pinch at lateral corticospinal tract.
Diagnostics 12 02090 g005
Figure 6. Schematic interpretation process of MR parameters.
Figure 6. Schematic interpretation process of MR parameters.
Diagnostics 12 02090 g006
Table 1. Clinical information of the patients.
Table 1. Clinical information of the patients.
Patient No.12345678
Level of maximal compressionC3/4C5/6C4/5C4/5C6/7C5/6C5/6C6/7
Compression ratio at level of maximal compression0.300.330.390.400.490.440.610.44
T2 high signal intensity+++-----
DN4 score22844354
JOA score1615131516161616
Hyperreflexia-+--++--
Lateral pinch (Rt./Lt.)19.0/19.619.0/19.69.0/7.320.3/11.616.6/16.818.3/17.319.6/17.020.3/19.6
Power grasp (Rt./Lt.)68.3/61.668.3/61.146.0/33.375.0/46.673.3/67.368.3/80.065.0/66.680.0/73.3
ASIA motor (Rt./Lt.)50/5050/5038/3745/4450/5050/5049/4949/49
ASIA sensory (Rt./Lt.)111/112112/112117/117112/112111/110112/112108/108110/110
ASIA light touch (Rt./Lt.)55/5656/5658/5856/5655/5556/5654/5455/55
ASIA pinprick (Rt./Lt.)56/5656/5659/5956/5656/5556/5654/5455/55
DN4, Douleur Neuropathique 4; JOA, Japanese Orthopaedic Association; ASIA, American Spinal Cord Injury Association.
Table 2. Demographic information of each group.
Table 2. Demographic information of each group.
Patient (n = 8)Control (n = 12)p-Value
Sex (male:female)8:08:40.238
Age (age ± SD)73.00 ± 5.7837.08 ± 5.380.000 *
SD, standard deviation; * statistically significant.
Table 3. Group differences of MR parameters of spinal cord tracts at each level.
Table 3. Group differences of MR parameters of spinal cord tracts at each level.
Above the Lesion LevelAt the Lesion LevelBelow the Lesion Level
ControlPatient ControlPatient ControlPatient
Mean ± SDMean ± SDpMean ± SDMean ± SDpMean ± SDMean ± SDp
Lat CSFA0.689 ± 0.083 0.663 ± 0.071 0.4640.648 ± 0.077 0.597 ± 0.114 0.2170.602 ± 0.073 0.631 ± 0.081 0.411
MD1101 × 10−6 ± 0.000 972 × 10−6 ± 0.000 0.1991074 × 10−6 ± 0.000 1125 × 10−6 ± 0.000 0.7461240 × 10−6 ± 0.001 0.001 ± 1061 × 10−60.232
AD1955 × 10−6 ± 0.000 1788 × 10−6 ± 0.000 0.2561849 × 10−6 ± 0.000 1950 × 10−6 ± 0.000 0.5332039 × 10−6 ± 0.001 0.002 ± 1917 × 10−60.434
RD641 × 10−6 ± 0.000 564 × 10−6 ± 0.000 0.367626 × 10−6 ± 0.000 712 × 10−6 ± 0.000 0.382756 × 10−6 ± 0.000 0.001 ± 633 × 10−6 0.146
MTR44.707 ± 3.955 46.966 ± 3.196 0.15643.300 ± 3.89346.393 ± 4.0330.06042.424 ± 4.48147.884 ± 5.5040.006 *
RuSFA0.748 ± 0.089 0.718 ± 0.088 0.4540.745 ± 0.100 0.614 ± 0.109 0.002 *0.655 ± 0.093 0.688 ± 0.088 0.436
MD859 × 10−6 ± 0.000 905 × 10−6 ± 0.000 0.635787 × 10−6 ± 0.000 114 × 10−6 ± 0.000 0.002 *979 × 10−6 ± 0.000 912 × 10−6 ± 0.000 0.492
AD1763 × 10−6 ± 0.000 1796 × 10−6 ± 0.000 0.8661670 × 10−6 ± 0.000 1957 × 10−6 ± 0.000 0.3001945 × 10−6 ± 0.000 1795 × 10−6 ± 0.000 0.419
RD425 × 10−6 ± 0.000 459 × 10−6 ± 0.000 0.638416 × 10−6 ± 0.000 687 × 10−6 ± 0.000 0.005 *653 × 10−6 ± 0.000 467 × 10−6 ± 0.000 0.082
MTR47.318 ± 4.182 48.942 ± 3.389 0.36745.676 ± 5.033 45.877 ± 4.156 0.95944.450 ± 3.60850.966 ± 5.3240.004 *
Lat RSFA0.796 ± 0.137 0.646 ± 0.082 0.001 *0.733 ± 0.124 0.564 ± 0.072 0.001 *0.781 ± 0.086 0.614 ± 0.057 0.001 *
MD646 × 10−6 ± 0.000 928 × 10−6 ± 0.000 0.000 *652 × 10−6 ± 0.000 1138 × 10−6 ± 0.000 0.000 *779 × 10−6 ± 0.000 1030 × 10−6 ± 0.000 0.015 *
AD1327 × 10−6 ± 0.000 1704 × 10−6 ± 0.000 0.002 *1297 × 10−6 ± 0.000 1919 × 10−6 ± 0.000 0.002 *1550 × 10−6 ± 0.000 1844 × 10−6 ± 0.000 0.006 *
RD254 × 10−6 ± 0.000 540 × 10−6 ± 0.000 0.000 *311 × 10−6 ± 0.000 746 × 10−6 ± 0.000 0.000 *295 × 10−6 ± 0.000 623 × 10−6 ± 0.000 0.000 *
MTR46.813 ± 3.449 47.413 ± 3.810 0.72446.309 ± 2.927 46.376 ± 3.712 0.98944.370 ± 4.27448.131 ± 5.0340.048 *
FCFA0.675 ± 0.070 0.676 ± 0.079 0.9890.597 ± 0.068 0.616 ± 0.101 0.6350.603 ± 0.060 0.602 ± 0.069 0.989
MD1166 × 10−6 ± 0.000920 × 10−6 ± 0.0000.006 *1124 × 10−6 ± 0.000 1079 × 10−6 ± 0.000 0.6351164 × 10−6 ± 0.000 1101 × 10−6 ± 0.000 0.532
AD2189 × 10−6 ± 0.0001753 × 10−6 ± 0.0000.014 *2033 × 10−6 ± 0.000 1934 × 10−6 ± 0.000 0.5692022 × 10−6 ± 0.000 1937 × 10−6 ± 0.000 0.618
RD680 × 10−6 ± 0.000503 × 10−6 ± 0.0000.006 *759 × 10−6 ± 0.000 651 × 10−6 ± 0.000 0.111767 × 10−6 ± 0.000 683 × 10−6 ± 0.000 0.256
MTR44.733 ± 2.768 46.912 ± 4.392 0.20743.319 ± 2.85145.640 ± 3.7740.08843.272 ± 4.441 45.287 ± 3.231 0.256
SpLmFA0.692 ± 0.100 0.654 ± 0.076 0.3670.695 ± 0.109 0.595 ± 0.067 0.009 *0.642 ± 0.086 0.658 ± 0.075 0.644
MD995 × 10−6 ± 0.000 959 × 10−6 ± 0.000 0.7461043 × 10−6 ± 0.001 1120 × 10−6 ± 0.000 0.6381070 × 10−6 ± 0.000 1047 × 10−6 ± 0.000 0.927
AD1718 × 10−6 ± 0.000 1722 × 10−6 ± 0.000 0.9891853 × 10−6 ± 0.001 1860 × 10−6 ± 0.000 0.9421862 × 10−6 ± 0.000 1892 × 10−6 ± 0.000 0.910
RD533 × 10−6 ± 0.000 577 × 10−6 ± 0.000 0.638559 × 10−6 ± 0.000 728 × 10−6 ± 0.000 0.232627 × 10−6 ± 0.000 625 × 10−6 ± 0.000 0.989
MTR44.048 ± 5.53247.904 ± 4.5000.07443.300 ± 3.528 44.623 ± 3.545 0.42943.168 ± 4.52749.285 ± 3.9900.001 *
Lat CS, lateral corticospinal tract; RuS, rubrospinal tract; Lat RS, lateral reticulospinal tract; FC, fasciculus cuneatus; SpLm, spinal lemniscus. Underlined values denote results that are contrary to expectation; * statistically significant.
Table 4. MR metrics displaying correlation with clinical items. Underline denotes an inverse correlation.
Table 4. MR metrics displaying correlation with clinical items. Underline denotes an inverse correlation.
Lat CSRuSLat RSFCSpinal Lemni
Power graspabove AD
at FA, MD, AD, RDAD MD, AD
below
Lateral pinchabove MTR
at
belowFA, MTR
Lat CS, lateral corticospinal tract; RuS, rubrospinal tract; Lat RS, lateral reticulospinal tract; FC, fasciculus cuneatus; Spinal Lemni, spinal lemniscus. FA and MTR: positive correlation was considered as normal correlation and negative correlation was considered as inverse correlation; MD, AD, and RD: negative correlation was considered as normal correlation and positive correlation was considered as inverse correlation.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yang, H.-E.; Kim, W.-T.; Kim, D.-H.; Kim, S.-W.; Yoo, W.-K. Utility of Diffusion and Magnetization Transfer MRI in Cervical Spondylotic Myelopathy: A Pilot Study. Diagnostics 2022, 12, 2090. https://doi.org/10.3390/diagnostics12092090

AMA Style

Yang H-E, Kim W-T, Kim D-H, Kim S-W, Yoo W-K. Utility of Diffusion and Magnetization Transfer MRI in Cervical Spondylotic Myelopathy: A Pilot Study. Diagnostics. 2022; 12(9):2090. https://doi.org/10.3390/diagnostics12092090

Chicago/Turabian Style

Yang, Hea-Eun, Wan-Tae Kim, Dae-Hyun Kim, Seok-Woo Kim, and Woo-Kyoung Yoo. 2022. "Utility of Diffusion and Magnetization Transfer MRI in Cervical Spondylotic Myelopathy: A Pilot Study" Diagnostics 12, no. 9: 2090. https://doi.org/10.3390/diagnostics12092090

APA Style

Yang, H. -E., Kim, W. -T., Kim, D. -H., Kim, S. -W., & Yoo, W. -K. (2022). Utility of Diffusion and Magnetization Transfer MRI in Cervical Spondylotic Myelopathy: A Pilot Study. Diagnostics, 12(9), 2090. https://doi.org/10.3390/diagnostics12092090

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop