A Prospective Longitudinal Study in Degenerative Cervical Myelopathy Using Quantitative Microstructural MRI with Tract-Specific Metrics
Allan R. Martin1, Benjamin De Leener2, Izabela Aleksanderek1, Julien Cohen-Adad2, David W. Cadotte1, Sukhvinder Kalsi-Ryan1, Lindsay Tetreault1, Adrian Crawley3, Howard Ginsberg1, David J. Mikulis3, and Michael G. Fehlings1

1Neurosurgery, University of Toronto, Toronto, ON, Canada, 2Electrical Engineering, Polytechnique Montreal, Montreal, QC, Canada, 3Medical Imaging, University of Toronto, Toronto, ON, Canada

Synopsis

This study investigates if DTI, MT, and T2*-weighted imaging of the rostral cervical cord can 1) detect injury of WM tracts, 2) correlate with global and focal disability, and 3) predict outcomes in degenerative cervical myelopathy (DCM). Data includes detailed clinical assessments, electrophysiology, and MRI, repeated at 1-year. Quantitative MRI in 37 DCM patients and 29 healthy controls provided reliable results and showed decreased CSA, FA, and MTR, and increased T2* WM/GM ratio. FA of individual tracts correlates well with clinical measures. Quantitative multimodal assessment of WM injury with a clinically feasible protocol is possible, with many potential clinical applications.

Purpose

Spinal cord MRI techniques are emerging that can characterize aspects of microstructure: diffusion tensor imaging (DTI), magnetization transfer (MT), and T2*-weighted imaging (T2*-WI). Degenerative cervical myelopathy (DCM) is a common condition involving multilevel disc degeneration, often with hypertrophy and/or ossification of ligaments within the spinal canal, leading to cord compression (Figure 1). In this prospective longitudinal study, we develop a clinically feasible multimodal MRI protocol, establish normative data, and investigate if these techniques can quantify injury to individual tracts rostral to cord compression, correlate with global and focal disability, and predict outcomes in DCM.

Methods

37 DCM patients (age 56.2; 54% male; 19 mild, 12 moderate, 6 severe; 2 with previous surgery) were enrolled consecutively from outpatient neurosurgery clinic. 29 age/gender-matched healthy subjects served as imaging controls. Data collection included detailed clinical assessments (mJOA, Nurick, MDI, QuickDASH, Berg Balance (BB), grip, upper extremity (UE) strength, modified GRASSP, GaitRITE), electrophysiology (ulnar SSEP), and quantitative MRI (3T GE) with 13 axial slices covering C1-C7 (Figure 2): single-shot EPI DTI (TR/TE=4050/75ms, b=800s/mm2, 25 directions, 1.25x1.25x5mm voxels, 3 acquisitions averaged offline, 2m06s each), axial 2D SPGR (TR/TE=36/5.7, FA 6, NEX=3, 1x1x5mm voxels)+/- MT pre-pulse (3m43s each), and MERGE T2*-WI (TR/TE=550/5,10,15, 0.6x0.6x4mm voxels, 4m15s), co-registered to a T2-weighted image (T2-WI) (6m30s) (Figure 3). Data collection will be repeated at 1-year follow-up. Semi-automated analysis was performed using the Spinal Cord Toolbox1 for segmentation, motion correction, registration to WM atlas, and extraction of metrics with automatic correction for partial volume effects (Figure 1). Reliability was analyzed by calculation of intra-class correlation coefficient (ICC) for 7 MT datasets (2 measurements), 66 DTI datasets (3 measurements), 13 DTI datasets using a 3-run average (2 measurements), and 8 T2*-WI datasets (2 measurements). 5 a priori metrics were selected: cord cross sectional area (CSA) at C1, fractional anisotropy (FA), mean diffusivity (MD), MT ratio (MTR), and T2*-WI WM/GM ratio (a novel metric representing grey-white contrast). Regions of interest (ROIs) in the dorsal columns (DCs) and lateral corticospinal tracts (LCSTs) at C1-C3 were also selected a priori, to avoid potential bias at compressed levels. Univariate/multivariate analysis of variance (ANOVA/MANOVA) was performed to compare DCM patients vs. controls. Spearman correlations were calculated between metrics and measures of global disability (mJOA, BB) and focal impairment (mJOA motor scores, UE strength, grip force, sensation).

Results

Reliability of quantitative assessments was excellent for all 4 metrics extracted from DCs and LCSTs (ICCs: 0.86-0.99) (Table 1). Metrics from individual DTI scans showed lower reliability (0.82-0.94), but averaging over 3 acquisitions improved reliability substantially (0.93-0.99). DCM patients demonstrated evidence of WM injury rostral to the compressed spinal cord (multivariate p<0.0001), including decreases in CSA, FA (DCs, LCSTs), and MTR (DCs, trend in LCSTs), and increased T2*-WM/GM ratio (DCs). MD did not show differences vs. controls. Cross-correlations of MTR, FA, and MD metrics were moderate, ranging from 0.2 to 0.5, whereas T2*-WM/GM ratio did not correlate closely with other metrics. FA provided the most robust results overall, showing highly significant differences vs. controls (p<0.0001 in DCs, LCs) and the strongest correlations among MRI metrics with all clinical measures (Table 2). The strongest correlations were also frequently found in the expected anatomical tract, with FA of the bilateral LCSTs correlating well with mJOA motor scores was (r=0.64), FA of the ipsilateral LCST predicting arm power (left: r=0.44, right: r=0.53) and grip strength (left: r=0.57, right: r=0.54), FA of the DCs correlating well with ipsilateral sensation (left: r=0.66, right: r=0.67), and FA in the DCs correlating with BB (r=0.51) (all p<0.05). Furthermore, 2 patients who had undergone previous surgery with metallic implants (lower cervical spine) showed data of acceptable quality at C1-C3.

Discussion

We have developed a reliable multimodal spinal cord MRI protocol that is feasible to implement in a clinical environment, using standard clinical hardware and with an acquisition time of approximately 30 minutes. Bias has been minimized by utilizing automated tools for quantitative assessment and focusing on the rostral portion of the cervical cord, avoiding difficulties with cord segmentation and susceptibility artefact in the compressed cord and enabling assessment of patients with surgically implanted hardware. WM injury of individual tracts can be successfully quantified in the rostral spinal cord, correlating well with focal neurological deficits and global impairment. Each of our imaging techniques appears to provide complimentary information, and multivariate linear regression analysis of a larger cohort is underway.

Conclusions

Our results demonstrate that accurate and reliable quantitative assessment of white matter damage can be achieved with clinically relevant methods, which could be of benefit for all neurological disorders affecting the spinal cord.

Acknowledgements

Protocol development and analysis techniques were supported by a grant from Rick Hansen Institute: “Development of MRI-based Biomarkers in Patients with Acute Spinal Cord Injury”.

Allan R. Martin received salary support that allowed him to lead this study as part of his PhD thesis project, from the Ontario Ministry of Health Clinician Investigator Program.

References

1. De Leener B, Roux A, Touati J, Levy1, Taso M, Fonov V, Collins DL, Callot V, Cohen-Adad J. Template-based analysis of multi-parametric MRI data with the Spinal Cord Toolbox. Proc. ISMRM, Toronto, Canada 2015.

Figures

Table 1: Summary of Results of Multimodal Quantitative Assessment of the Spinal Cord at C1-C3. A multivariate comparison of healthy controls (N=29) and DCM subjects (N=37) using MANOVA. Units for MD are 10-3 mm2/s. Legend: * denotes univariate trend of p<0.10; ** denotes univariate significance of p<0.05 (uncorrected); *** denotes corrected univariate significance of p<0.0056; **** denotes multivariate significance of p<0.05.

Table 2: Correlation of MRI Metrics with Clinical Assessment Tools. Spearman Rank Correlation Coefficients with p<0.10 are displayed; * Denotes significance of p<0.05 (uncorrected).

Figure 1: Pathogenesis of Degenerative Cervical Myelopathy. Left: T2-weighted sagittal MRI image with disc herniations at C5-6, C6-7 indenting the anterior spinal cord. Right: gross pathological section of herniated disc material indenting spinal cord.

Figure 2: Slice Positioning. Left: prescription of 13 axial slices from C1 to C7, alternating between mid-vertebral body and mid-intervertebral disc. Right: T2*-weighted imaged at C5 showing bony osteophyte causing anterior spinal cord compression.

Figure 3: Multimodal MRI Images with Probabilistic White Matter Tracts. Representative images of a healthy control including FA map (A), MTR map (B), and T2*-WI (C), and corresponding images with probabilistic map of lateral corticospinal tracts (LCSTs) overlaid in blue and dorsal columns (DCs) in red-yellow (D-F).



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
0853