SCT: Spinal Cord Toolbox, an open-source software for processing spinal cord MRI data

A Follow-up Study of Postoperative DCM Patients Using Diffusion MRI with DTI and NODDI

STUDY DESIGN

A retrospective cohort study.

OBJECTIVE

The aim of this study was to investigate the relationship between spinal cord microstructures and spinal cord dysfunction in degenerative cervical myelopathy (DCM) patients; a follow-up study was carried out using diffusion tensor imaging (DTI) and neurite orientation dispersion and density imaging (NODDI), an advanced diffusion model.

SUMMARY OF BACKGROUND DATA

DTI has been used for diagnosis of DCM. Although DTI has advantages of high sensitivity and quantitative features, it is limited by its low specificity for measuring microstructures.

METHODS

Fifty-eight postoperative DCM patients with 12 to 14 months of surgical recovery were recruited, as well as 14 healthy volunteers for control group. All DTI and NODDI metrics were measured at the most stenotic levels of patients and at all levels of control group. A t test was used to compare the metrics between patient and control groups, and Spearman correlation was used to test the relationship between the metrics and clinical assessment, modified Japanese Orthopedic Association (mJOA) scores.

RESULTS

t test showed that DTI metrics, neurite density Vic, and free water fraction Viso had significant differences between control group and the most stenotic levels in patients. For DTI metrics, fractional anisotropy (FA), radial diffusivity (RD), and mean diffusivity (MD) in patients were significantly correlated with mJOA. For NODDI metrics, only Vic is positively correlated with mJOA.

CONCLUSION

The results of t test and correlation with mJOA suggest that DTI and NODDI are valuable for evaluating spinal cord function. Results of NODDI indicate that the reason for FA reduction in DCM may be decreased neurite density, not increased orientation dispersion.

LEVEL OF EVIDENCE

3.

Filter-probe diffusion imaging improves spinal cord injury outcome prediction

OBJECTIVE

Diffusion-weighted imaging (DWI) is a powerful tool for investigating spinal cord injury (SCI), but has limited specificity for axonal damage, which is the most predictive feature of long-term functional outcome. In this study, a technique designed to detect acute axonal injury, filter-probe double diffusion encoding (FP-DDE), is compared with standard DWI for predicting long-term functional and cellular outcomes.

METHODS

This study extends FP-DDE to predict long-term functional and histological outcomes in a rat SCI model of varying severities (n = 58). Using a 9.4T magnetic resonance imaging (MRI) system, a whole-cord FP-DDE spectroscopic voxel was acquired in 3 minutes at the lesion site and compared to DWI at 48 hours postinjury. Relationships with chronic (30-day) locomotor and histological outcomes were evaluated with linear regression.

RESULTS

The FP-DDE measure of parallel diffusivity (ADC|| ) was significantly related to chronic hind limb locomotor functional outcome (R2  = 0.63, p < 0.0001), and combining this measurement with acute functional scores demonstrated prognostic benefit versus functional testing alone (p = 0.0007). Acute ADC|| measurements were also more closely related to the number of injured axons measured 30 days after the injury than standard DWI. Furthermore, acute FP-DDE images showed a clear and easily interpretable pattern of injury that closely corresponded with chronic MRI and histology observations.

INTERPRETATION

Collectively, these results demonstrate FP-DDE benefits from greater specificity for acute axonal damage in predicting functional and histological outcomes with rapid acquisition and fully automated analysis, improving over standard DWI. FP-DDE is a promising technique compatible with clinical settings, with potential research and clinical applications for evaluation of spinal cord pathology. Ann Neurol 2018;83:37-50.

Focal and diffuse cervical spinal cord damage in patients with early relapsing-remitting MS: A multicentre magnetisation transfer ratio study

BACKGROUND

Studies including patients with well-established multiple sclerosis (MS) have shown a significant and disability-related reduction in the cervical spinal cord (SC) magnetisation transfer ratio (MTR).

OBJECTIVES

The objectives are to (1) assess whether MTR reduction is already measurable in the SC of patients with early relapsing-remitting multiple sclerosis (RRMS) and (2) describe its spatial distribution.

METHODS

We included 60 patients with RRMS <12  months and 34 age-matched controls at five centres. Axial T2*w, sagittal T2w, sagittal phase-sensitive inversion recovery (PSIR), 3DT1w, and axial magnetisation transfer (MT) images were acquired from C1 to C7. Lesions were manually labelled and mean MTR values computed both for the whole SC and for normal-appearing SC in different regions of interest.

RESULTS

Mean whole SC MTR was significantly lower in patients than controls (33.7 vs 34.9  pu, p  =  0.00005), even after excluding lesions (33.9  pu, p  =  0.0003). We observed a greater mean reduction in MTR for vertebral levels displaying the highest lesion loads (C2-C4). In the axial plane, we observed a greater mean MTR reduction at the SC periphery and barycentre.

CONCLUSION

Cervical SC tissue damage measured using MTR is not restricted to macroscopic lesions in patients with early RRMS and is not homogeneously distributed.

Quantitative Magnetization Transfer MRI Measurements of the Anterior Spinal Cord Region are Associated With Clinical Outcomes in Cervical Spondylotic Myelopathy

STUDY DESIGN

A case-control study.

OBJECTIVE

The aim of this study was to understand the role of magnetization transfer ratio (MTR) in identifying patients with clinically significant myelopathy and disability.

SUMMARY OF BACKGROUND DATA

MTR is a quantitative measure that correlates with myelin loss and neural tissue destruction in a variety of neurological diseases. However, the usefulness of MTR in patients with cervical spondylotic myelopathy (CSM) has not been examined.

METHODS

We prospectively enrolled seven CSM patients and seven age-matched controls to undergo magnetic resonance imaging (MRI) of the cervical spine. Nurick, Neck Disability Index (NDI), and modified Japanese Orthopedic Association (mJOA) scores were collected for all patients. Clinical hyperreflexia was tested at the MCP joint, using a six-axis load cell. Reflex was simulated by quickly moving the joint from maximum flexion to maximum extension (300°/second). Anterior, lateral, and posterior cord MTR measurements were compared with clinical outcomes.

RESULTS

Compared with controls, CSM patients had lower anterior cord MTR (38.29 vs. 29.97, Δ = -8.314, P = 0.0022), and equivalent posterior cord (P = 0.2896) and lateral cord (P = 0.3062) MTR. Higher Nurick scores were associated with lower anterior cord MTR (P = 0.0205), but not lateral cord (P = 0.5446) or posterior cord MTR (P = 0.1222). Lower mJOA was associated with lower anterior cord MTR (P = 0.0090), but not lateral cord (P = 0.4864) or posterior cord MTR (P = 0.4819). There was no association between NDI and MTR of the anterior (P = 0.4351), lateral (P = 0.7557), or posterior cord (P = 0.9171). There was a linear relationship between hyperreflexia and anterior cord MTR (slope = -117.3, R = 0.6598, P = 0.0379), but not lateral cord (P = 0.1906, R = 0.4511) or posterior cord (P = 0.2577, R = 0.3957) MTR.

CONCLUSION

Anterior cord MTR correlates with clinical outcomes as measured by mJOA index, Nurick score, and quantitative hyperreflexia, and could play a role in the preoperative assessment of CSM.

LEVEL OF EVIDENCE

2.

Scan-rescan of axcaliber, macromolecular tissue volume, and g-ratio in the spinal cord

PURPOSE

Recent MRI techniques have been introduced that can extract microstructural information in the white matter, such as the density or macromolecular content. Translating quantitative MRI to the clinic raises many challenges in terms of acquisition strategy, modeling of the MRI signal, artifact corrections, and metric extraction (template registration and partial volume effects). In this work, we investigated the scan-rescan repeatability of several quantitative MRI techniques in the human spinal cord.

METHODS

AxCaliber metrics, macromolecular tissue volume, and the fiber g-ratio were estimated in the spinal cord of eight healthy subjects, scanned and rescanned the same day in two different sessions.

RESULTS

Scan-rescan repeatability deviation was 3% for all metrics, in average in the white matter of all subjects. Intraclass correlation coefficient was up to 0.9. A three-way analysis of variance showed significant effects of white matter pathway, laterality, and subject.

CONCLUSION

The present study suggests that quantitative MRI gives stable measurements of white matter microstructure in the spinal cord of healthy subjects. Our findings remain to be evaluated in diseased populations. Magn Reson Med 79:2759-2765, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Neuroinflammation of the spinal cord and nerve roots in chronic radicular pain patients

Numerous preclinical studies support the role of spinal neuroimmune activation in the pathogenesis of chronic pain, and targeting glia (eg, microglia/astrocyte)- or macrophage-mediated neuroinflammatory responses effectively prevents or reverses the establishment of persistent nocifensive behaviors in laboratory animals. However, thus far, the translation of those findings into novel treatments for clinical use has been hindered by the scarcity of data supporting the role of neuroinflammation in human pain. Here, we show that patients suffering from a common chronic pain disorder (lumbar radiculopathy), compared with healthy volunteers, exhibit elevated levels of the neuroinflammation marker 18 kDa translocator protein, in both the neuroforamina (containing dorsal root ganglion and nerve roots) and spinal cord. These elevations demonstrated a pattern of spatial specificity correlating with the patients' clinical presentation, as they were observed in the neuroforamen ipsilateral to the symptomatic leg (compared with both contralateral neuroforamen in the same patients as well as to healthy controls) and in the most caudal spinal cord segments, which are known to process sensory information from the lumbosacral nerve roots affected in these patients (compared with more superior segments). Furthermore, the neuroforaminal translocator protein signal was associated with responses to fluoroscopy-guided epidural steroid injections, supporting its role as an imaging marker of neuroinflammation, and highlighting the clinical significance of these observations. These results implicate immunoactivation at multiple levels of the nervous system as a potentially important and clinically relevant mechanism in human radicular pain, and suggest that therapies targeting immune cell activation may be beneficial for chronic pain patients.

Future Brain and Spinal Cord Volumetric Imaging in the Clinic for Monitoring Treatment Response in MS

PURPOSE OF REVIEW

Volumetric analysis of brain imaging has emerged as a standard approach used in clinical research, e.g., in the field of multiple sclerosis (MS), but its application in individual disease course monitoring is still hampered by biological and technical limitations. This review summarizes novel developments in volumetric imaging on the road towards clinical application to eventually monitor treatment response in patients with MS.

RECENT FINDINGS

In addition to the assessment of whole-brain volume changes, recent work was focused on the volumetry of specific compartments and substructures of the central nervous system (CNS) in MS. This included volumetric imaging of the deep brain structures and of the spinal cord white and gray matter. Volume changes of the latter indeed independently correlate with clinical outcome measures especially in progressive MS. Ultrahigh field MRI and quantitative MRI added to this trend by providing a better visualization of small compartments on highly resolving MR images as well as microstructural information. New developments in volumetric imaging have the potential to improve sensitivity as well as specificity in detecting and hence monitoring disease-related CNS volume changes in MS.

Monitoring for myelopathic progression with multiparametric quantitative MRI

Background

Patients with mild degenerative cervical myelopathy (DCM) are often managed non-operatively, and surgery is recommended if neurological progression occurs. However, detection of progression is often subjective. Quantitative MRI (qMRI) directly measures spinal cord (SC) tissue changes, detecting axonal injury, demyelination, and atrophy. This longitudinal study compared multiparametric qMRI with clinical measures of progression in non-operative DCM patients.

Methods

26 DCM patients were followed. Clinical data included modified Japanese Orthopedic Association (mJOA) and additional assessments. 3T qMRI data included cross sectional area, diffusion fractional anisotropy, magnetization transfer ratio, and T2*-weighted white/grey matter signal ratio, extracted from the compressed SC and above/below. Progression was defined as 1) patients’ subjective impression, 2) 2-point mJOA decrease, 3) ≥3 clinical measures worsening ≥5%, 4) increased compression on MRI, or 5) ≥1 of 10 qMRI measures or composite score worsening (p < 0.004, corrected).

Results

Follow-up (13.5 ± 4.9 months) included mJOA in all 26 patients, MRI in 25, and clinical/qMRI in 22. 42.3% reported subjective worsening, compared with mJOA (11.5%), MRI (20%), comprehensive assessments (54.6%), and qMRI (68.2%). Relative to subjective worsening, qMRI showed 100% sensitivity and 53.3% specificity compared with comprehensive assessments (75%, 60%), mJOA (27.3%, 100%), and MRI (18.2%, 81.3%). A decision-making algorithm incorporating qMRI identified progression and recommended surgery for 11 subjects (42.3%).

Conclusions

Quantitative MRI shows high sensitivity to detect myelopathic progression. Our results suggest that neuroplasticity and behavioural adaptation may mask progressive SC tissue injury. qMRI appears to be a useful method to confirm subtle myelopathic progression in individual patients, representing an advance toward clinical translation of qMRI.

Can microstructural MRI detect subclinical tissue injury in subjects with asymptomatic cervical spinal cord compression? A prospective cohort study

OBJECTIVES

Degenerative cervical myelopathy (DCM) involves extrinsic spinal cord compression causing tissue injury and neurological dysfunction. Asymptomatic spinal cord compression (ASCC) is more common, but its significance is poorly defined. This study investigates if: (1) ASCC can be automatically diagnosed using spinal cord shape analysis; (2) multiparametric quantitative MRI can detect similar spinal cord tissue injury as previously observed in DCM.

DESIGN

Prospective observational longitudinal cohort study.

SETTING

Single centre, tertiary care and research institution.

PARTICIPANTS

40 neurologically intact subjects (19 female, 21 male) divided into groups with and without ASCC.

INTERVENTIONS

None.

OUTCOME MEASURES

Clinical assessments: modified Japanese Orthopaedic Association score and physical examination. 3T MRI assessments: automated morphometric analysis compared with consensus ratings of spinal cord compression, and measures of tissue injury: cross-sectional area, diffusion fractional anisotropy, magnetisation transfer ratio and T2*-weighted imaging white to grey matter signal intensity ratio (T2*WI WM/GM) extracted from rostral (C1-3), caudal (C6-7) and maximally compressed levels.

RESULTS

ASCC was present in 20/40 subjects. Diagnosis with automated shape analysis showed area under the curve >97%. Five MRI metrics showed differences suggestive of tissue injury in ASCC compared with uncompressed subjects (p<0.05), while a composite of all 10 measures (average of z scores) showed highly significant differences (p=0.002). At follow-up (median 21 months), two ASCC subjects developed DCM.

CONCLUSIONS

ASCC appears to be common and can be accurately and objectively diagnosed with automated morphometric analysis. Quantitative MRI appears to detect subclinical tissue injury in ASCC prior to the onset of neurological symptoms and signs. These findings require further validation, but offer the intriguing possibility of presymptomatic diagnosis and treatment of DCM and other spinal pathologies.

Spinal cord gray matter segmentation using deep dilated convolutions

Gray matter (GM) tissue changes have been associated with a wide range of neurological disorders and were recently found relevant as a biomarker for disability in amyotrophic lateral sclerosis. The ability to automatically segment the GM is, therefore, an important task for modern studies of the spinal cord. In this work, we devise a modern, simple and end-to-end fully-automated human spinal cord gray matter segmentation method using Deep Learning, that works both on in vivo and ex vivo MRI acquisitions. We evaluate our method against six independently developed methods on a GM segmentation challenge. We report state-of-the-art results in 8 out of 10 evaluation metrics as well as major network parameter reduction when compared to the traditional medical imaging architectures such as U-Nets.

Microstructural imaging in the spinal cord and validation strategies

In vivo histology using magnetic resonance imaging (MRI) is a newly emerging research field that aims to non-invasively characterize tissue microstructure. The implications of in vivo histology are many, from discovering novel biomarkers to studying human development, to providing tools for disease diagnosis and monitoring the effects of novel treatments on tissue. This review focuses on quantitative MRI (qMRI) techniques that are used to map spinal cord microstructure. Opening with a rationale for non-invasive imaging of the spinal cord, this article continues with a brief overview of the existing MRI techniques for axon and myelin imaging, followed by the specific challenges and potential solutions for acquiring and processing such data. The final part of this review focuses on histological validation, with suggested tissue preparation, acquisition and processing protocols for large-scale microscopy.

Dorsal and ventral horn atrophy is associated with clinical outcome after spinal cord injury

OBJECTIVE

To investigate whether gray matter pathology above the level of injury, alongside white matter changes, also contributes to sensorimotor impairments after spinal cord injury.

METHODS

A 3T MRI protocol was acquired in 17 tetraplegic patients and 21 controls. A sagittal T2-weighted sequence was used to characterize lesion severity. At the C2-3 level, a high-resolution T2*-weighted sequence was used to assess cross-sectional areas of gray and white matter, including their subcompartments; a diffusion-weighted sequence was used to compute voxel-based diffusion indices. Regression models determined associations between lesion severity and tissue-specific neurodegeneration and associations between the latter with neurophysiologic and clinical outcome.

RESULTS

Neurodegeneration was evident within the dorsal and ventral horns and white matter above the level of injury. Tract-specific neurodegeneration was associated with prolonged conduction of appropriate electrophysiologic recordings. Dorsal horn atrophy was associated with sensory outcome, while ventral horn atrophy was associated with motor outcome. White matter integrity of dorsal columns and corticospinal tracts was associated with daily-life independence.

CONCLUSION

Our results suggest that, next to anterograde and retrograde degeneration of white matter tracts, neuronal circuits within the spinal cord far above the level of injury undergo transsynaptic neurodegeneration, resulting in specific gray matter changes. Such improved understanding of tissue-specific cord pathology offers potential biomarkers with more efficient targeting and monitoring of neuroregenerative (i.e., white matter) and neuroprotective (i.e., gray matter) agents.

Application of Quantitative Microstructural MR Imaging with Atlas-based Analysis for the Spinal Cord in Cervical Spondylotic Myelopathy

Mapping of MR fiber g-ratio, which is the ratio of the diameter of the axon to the diameter of the neuronal fiber, is introduced in this article. We investigated the MR fiber g-ratio, the axon volume fraction (AVF) and the myelin volume fraction (MVF) to evaluate microstructural changes in the spinal cord in patients with cervical spondylotic myelopathy (CSM) in vivo, using atlas-based analysis. We used diffusion MRI data acquired with a new simultaneous multi-slice accelerated readout-segmented echo planar imaging sequence for diffusion analysis for AVF calculation and magnetization transfer saturation imaging for MVF calculation. The AVFs of fasciculus gracilis in the affected side spinal cord, fasciculus cuneatus and lateral corticospinal tracts (LSCT) in the affected and unaffected side spinal cord were significantly lower (P = 0.019, 0.001, 0019, 0.000, and 0.002, respectively) than those of normal controls. No difference was found in the MVFs. The fiber g-ratio of LSCT was significantly lower (P = 0.040) in the affected side spinal cords than in the normal controls. The pathological microstructural changes in the spinal cord in patients with CSM, presumably partial axonal degenerations with preserved myelin. This technique has the potential to be a clinical biomarker in patients with CSM in vivo.

Spinal cord MRI at 7T

Magnetic resonance imaging (MRI) of the human spinal cord at 7T has been demonstrated by a handful of research sites worldwide, and the spinal cord remains one of the areas in which higher fields and resolution could have high impact. The small diameter of the cord (∼1 cm) necessitates high spatial resolution to minimize partial volume effects between gray and white matter, and so MRI of the cord can greatly benefit from increased signal-to-noise ratio and contrasts at ultra-high field (UHF). Herein we review the current state of UHF spinal cord imaging. Technical challenges to successful UHF spinal cord MRI include radiofrequency (B1) nonuniformities and a general lack of optimized radiofrequency coils, amplified physiological noise, and an absence of methods for robust B0 shimming along the cord to mitigate image distortions and signal losses. Numerous solutions to address these challenges have been and are continuing to be explored, and include novel approaches for signal excitation and acquisition, dynamic shimming and specialized shim coils, and acquisitions with increased coverage or optimal slice angulations.

Analysis strategies for high-resolution UHF-fMRI data

Functional MRI (fMRI) benefits from both increased sensitivity and specificity with increasing magnetic field strength, making it a key application for Ultra-High Field (UHF) MRI scanners. Most UHF-fMRI studies utilize the dramatic increases in sensitivity and specificity to acquire high-resolution data reaching sub-millimeter scales, which enable new classes of experiments to probe the functional organization of the human brain. This review article surveys advanced data analysis strategies developed for high-resolution fMRI at UHF. These include strategies designed to mitigate distortion and artifacts associated with higher fields in ways that attempt to preserve spatial resolution of the fMRI data, as well as recently introduced analysis techniques that are enabled by these extremely high-resolution data. Particular focus is placed on anatomically-informed analyses, including cortical surface-based analysis, which are powerful techniques that can guide each step of the analysis from preprocessing to statistical analysis to interpretation and visualization. New intracortical analysis techniques for laminar and columnar fMRI are also reviewed and discussed. Prospects for single-subject individualized analyses are also presented and discussed. Altogether, there are both specific challenges and opportunities presented by UHF-fMRI, and the use of proper analysis strategies can help these valuable data reach their full potential.

Quantifying multiple sclerosis pathology in post mortem spinal cord using MRI

Multiple sclerosis (MS) is a common inflammatory, demyelinating and degenerative disease of the central nervous system. The majority of people with MS present with symptoms due to spinal cord damage, and in more advanced MS a clinical syndrome resembling that of progressive myelopathy is not uncommon. Significant efforts have been undertaken to predict MS-related disability based on short-term observations, for example, the spinal cord cross-sectional area measured using MRI. The histo-pathological correlates of spinal cord MRI changes in MS are incompletely understood, however a surge of interest in tissue microstructure has recently led to new approaches to improve the precision with which MRI indices relate to underlying tissue features, such as myelin content, neurite density and orientation, among others. Quantitative MRI techniques including T₁ and T₂, magnetisation transfer (MT) and a number of diffusion-derived indices have all been successfully applied to post mortem MS spinal cord. Combining advanced quantification of histological features with quantitative - particularly diffusion-based - MRI techniques provide a new platform for high-quality MR/pathology data generation. To more accurately quantify grey matter pathology in the MS spinal cord, a key driver of physical disability in advanced MS, remains an important challenge of microstructural imaging.

Test-retest reliability of myelin imaging in the human spinal cord: Measurement errors versus region- and aging-induced variations

PURPOSE

To implement a statistical framework for assessing the precision of several quantitative MRI metrics sensitive to myelin in the human spinal cord: T1, Magnetization Transfer Ratio (MTR), saturation imposed by an off-resonance pulse (MTsat) and Macromolecular Tissue Volume (MTV).

METHODS

Thirty-three healthy subjects within two age groups (young, elderly) were scanned at 3T. Among them, 16 underwent the protocol twice to assess repeatability. Statistical reliability indexes such as the Minimal Detectable Change (MDC) were compared across metrics quantified within different cervical levels and white matter (WM) sub-regions. The differences between pathways and age groups were quantified and interpreted in context of the test-retest repeatability of the measurements.

RESULTS

The MDC was respectively 105.7ms, 2.77%, 0.37% and 4.08% for T1, MTR, MTsat and MTV when quantified over all WM, while the standard-deviation across subjects was 70.5ms, 1.34%, 0.20% and 2.44%. Even though particular WM regions did exhibit significant differences, these differences were on the same order as test-retest errors. No significant difference was found between age groups for all metrics.

CONCLUSION

While T1-based metrics (T1 and MTV) exhibited better reliability than MT-based measurements (MTR and MTsat), the observed differences between subjects or WM regions were comparable to (and often smaller than) the MDC. This makes it difficult to determine if observed changes are due to variations in myelin content, or simply due to measurement error. Measurement error remains a challenge in spinal cord myelin imaging, but this study provides statistical guidelines to standardize the field and make it possible to conduct large-scale multi-center studies.

Optimizing Filter-Probe Diffusion Weighting in the Rat Spinal Cord for Human Translation

Diffusion tensor imaging (DTI) is a promising biomarker of spinal cord injury (SCI). In the acute aftermath, DTI in SCI animal models consistently demonstrates high sensitivity and prognostic performance, yet translation of DTI to acute human SCI has been limited. In addition to technical challenges, interpretation of the resulting metrics is ambiguous, with contributions in the acute setting from both axonal injury and edema. Novel diffusion MRI acquisition strategies such as double diffusion encoding (DDE) have recently enabled detection of features not available with DTI or similar methods. In this work, we perform a systematic optimization of DDE using simulations and an in vivo rat model of SCI and subsequently implement the protocol to the healthy human spinal cord. First, two complementary DDE approaches were evaluated using an orientationally invariant or a filter-probe diffusion encoding approach. While the two methods were similar in their ability to detect acute SCI, the filter-probe DDE approach had greater predictive power for functional outcomes. Next, the filter-probe DDE was compared to an analogous single diffusion encoding (SDE) approach, with the results indicating that in the spinal cord, SDE provides similar contrast with improved signal to noise. In the SCI rat model, the filter-probe SDE scheme was coupled with a reduced field of view (rFOV) excitation, and the results demonstrate high quality maps of the spinal cord without contamination from edema and cerebrospinal fluid, thereby providing high sensitivity to injury severity. The optimized protocol was demonstrated in the healthy human spinal cord using the commercially-available diffusion MRI sequence with modifications only to the diffusion encoding directions. Maps of axial diffusivity devoid of CSF partial volume effects were obtained in a clinically feasible imaging time with a straightforward analysis and variability comparable to axial diffusivity derived from DTI. Overall, the results and optimizations describe a protocol that mitigates several difficulties with DTI of the spinal cord. Detection of acute axonal damage in the injured or diseased spinal cord will benefit the optimized filter-probe diffusion MRI protocol outlined here.

Automatic spinal cord localization, robust to MRI contrasts using global curve optimization

During the last two decades, MRI has been increasingly used for providing valuable quantitative information about spinal cord morphometry, such as quantification of the spinal cord atrophy in various diseases. However, despite the significant improvement of MR sequences adapted to the spinal cord, automatic image processing tools for spinal cord MRI data are not yet as developed as for the brain. There is nonetheless great interest in fully automatic and fast processing methods to be able to propose quantitative analysis pipelines on large datasets without user bias. The first step of most of these analysis pipelines is to detect the spinal cord, which is challenging to achieve automatically across the broad range of MRI contrasts, field of view, resolutions and pathologies. In this paper, a fully automated, robust and fast method for detecting the spinal cord centerline on MRI volumes is introduced. The algorithm uses a global optimization scheme that attempts to strike a balance between a probabilistic localization map of the spinal cord center point and the overall spatial consistency of the spinal cord centerline (i.e. the rostro-caudal continuity of the spinal cord). Additionally, a new post-processing feature, which aims to automatically split brain and spine regions is introduced, to be able to detect a consistent spinal cord centerline, independently from the field of view. We present data on the validation of the proposed algorithm, known as "OptiC", from a large dataset involving 20 centers, 4 contrasts (T₂-weighted n = 287, T₁-weighted n = 120, T₂^∗-weighted n = 307, diffusion-weighted n = 90), 501 subjects including 173 patients with a variety of neurologic diseases. Validation involved the gold-standard centerline coverage, the mean square error between the true and predicted centerlines and the ability to accurately separate brain and spine regions. Overall, OptiC was able to cover 98.77% of the gold-standard centerline, with a mean square error of 1.02 mm. OptiC achieved superior results compared to a state-of-the-art spinal cord localization technique based on the Hough transform, especially on pathological cases with an averaged mean square error of 1.08 mm vs. 13.16 mm (Wilcoxon signed-rank test p-value < .01). Images containing brain regions were identified with a 99% precision, on which brain and spine regions were separated with a distance error of 9.37 mm compared to ground-truth. Validation results on a challenging dataset suggest that OptiC could reliably be used for subsequent quantitative analyses tasks, opening the door to more robust analysis on pathological cases.

Region-specific impairment of the cervical spinal cord (SC) in amyotrophic lateral sclerosis: A preliminary study using SC templates and quantitative MRI (diffusion tensor imaging/inhomogeneous magnetization transfer)

In this preliminary study, our objective was to investigate the potential of high-resolution anatomical imaging, diffusion tensor imaging (DTI) and conventional/inhomogeneous magnetization transfer imaging [magnetization transfer (MT)/inhomogeneous magnetization transfer (ihMT)] at 3 T, analyzed with template-extracted regions of interest, to measure the atrophy and structural changes of white (WM) and gray (GM) matter spinal cord (SC) occurring in patients with amyotrophic lateral sclerosis (ALS). Ten patients with ALS and 20 age-matched healthy controls were recruited. SC GM and WM areas were automatically segmented using dedicated templates. Atrophy indices were evaluated from T₂ *-weighted images at each vertebral level from cervical C1 to C6. DTI and ihMT metrics were quantified within the corticospinal tract (CST), posterior sensory tract (PST) and anterior GM (aGM) horns at the C2 and C5 levels. Clinical disabilities of patients with ALS were evaluated using the Revised ALS Functional Rating Scale, upper motor neuron (UMN) and Medical Research Council scorings, and correlated with MR metrics. Compared with healthy controls, GM and WM atrophy was observed in patients with ALS, especially at lower cervical levels, where a strong correlation was also observed between GM atrophy and the UMN score (R = -0.75, p = 0.05 at C6). Interestingly, a significant decrease in ihMT ratio was found in all regions of interest (p < 0.0008), fractional anisotropy (FA) and MT ratios decreased significantly in CST, especially at C5 (p < 0.005), and λ_// (axial diffusivity) decreased significantly in CST (p = 0.0004) and PST (p = 0.003) at C2. Strong correlations between MRI metrics and clinical scores were also found (0.47 < |R| < 0.87, p < 0.05). Altogether, these preliminary results suggest that high-resolution anatomical imaging and ihMT imaging, in addition to DTI, are valuable for the characterization of SC tissue impairment in ALS. In this study, in addition to an important SC WM demyelination, we also observed, for the first time in ALS, impairments of cervical aGM.

Spinal Cord Gray Matter Atrophy in Amyotrophic Lateral Sclerosis

BACKGROUND AND PURPOSE

There is an emerging need for biomarkers to better categorize clinical phenotypes and predict progression in amyotrophic lateral sclerosis. This study aimed to quantify cervical spinal gray matter atrophy in amyotrophic lateral sclerosis and investigate its association with clinical disability at baseline and after 1 year.

MATERIALS AND METHODS

Twenty-nine patients with amyotrophic lateral sclerosis and 22 healthy controls were scanned with 3T MR imaging. Standard functional scale was recorded at the time of MR imaging and after 1 year. MR imaging data were processed automatically to measure the spinal cord, gray matter, and white matter cross-sectional areas. A statistical analysis assessed the difference in cross-sectional areas between patients with amyotrophic lateral sclerosis and controls, correlations between spinal cord and gray matter atrophy to clinical disability at baseline and at 1 year, and prediction of clinical disability at 1 year.

RESULTS

Gray matter atrophy was more sensitive to discriminate patients with amyotrophic lateral sclerosis from controls (P = .004) compared with spinal cord atrophy (P = .02). Gray matter and spinal cord cross-sectional areas showed good correlations with clinical scores at baseline (R = 0.56 for gray matter and R = 0.55 for spinal cord; P < .01). Prediction at 1 year with clinical scores (R² = 0.54) was improved when including a combination of gray matter and white matter cross-sectional areas (R² = 0.74).

CONCLUSIONS

Although improvements over spinal cord cross-sectional areas were modest, this study suggests the potential use of gray matter cross-sectional areas as an MR imaging structural biomarker to monitor the evolution of amyotrophic lateral sclerosis.

Imaging spinal cord atrophy in progressive myelopathies: HTLV-I-associated neurological disease (HAM/TSP) and multiple sclerosis (MS)

OBJECTIVE

Previous work measures spinal cord thinning in chronic progressive myelopathies, including human T-lymphotropic virus 1 (HTLV-1)-associated myelopathy/tropical spastic paraparesis (HAM/TSP) and multiple sclerosis (MS). Quantitative measurements of spinal cord atrophy are important in fully characterizing these and other spinal cord diseases. We aimed to investigate patterns of spinal cord atrophy and correlations with clinical markers.

METHODS

Spinal cord cross-sectional area was measured in individuals (24 healthy controls [HCs], 17 asymptomatic carriers of HTLV-1 (AC), 47 HAM/TSP, 74 relapsing-remitting MS [RRMS], 17 secondary progressive MS [SPMS], and 40 primary progressive MS [PPMS]) from C1 to T10. Clinical disability scores, viral markers, and immunological parameters were obtained for patients and correlated with representative spinal cord cross-sectional area regions at the C2 to C3, C4 to C5, and T4 to T9 levels. In 2 HAM/TSP patients, spinal cord cross-sectional area was measured over 3 years.

RESULTS

All spinal cord regions are thinner in HAM/TSP (56 mm² [standard deviation, 10], 59 [10], 23 [5]) than in HC (76 [7], 83 [8], 38 [4]) and AC (71 [7], 78 [9], 36 [7]). SPMS (62 [9], 66 [9], 32 [6]) and PPMS (65 [11], 68 [10], 35 [7]) have thinner cervical cords than HC and RRMS (73 [9], 77 [10], 37 [6]). Clinical disability scores (Expanded Disability Status Scale [p = 0.009] and Instituto de Pesquisas de Cananeia [p = 0.03]) and CD8⁺ T-cell frequency (p = 0.04) correlate with T4 to T9 spinal cord cross-sectional area in HAM/TSP. Higher cerebrospinal fluid HTLV-1 proviral load (p = 0.01) was associated with thinner spinal cord cross-sectional area. Both HAM/TSP patients followed longitudinally showed thoracic thinning followed by cervical thinning.

INTERPRETATION

Group average spinal cord cross-sectional area in HAM/TSP and progressive MS show spinal cord atrophy. We further hypothesize in HAM/TSP that is possible that neuroglial loss from a thoracic inflammatory process results in anterograde and retrograde degeneration of axons, leading to the temporal progression of thoracic to cervical atrophy described here. Ann Neurol 2017;82:719-728.

Thermal Stimulation Alters Cervical Spinal Cord Functional Connectivity in Humans

The spinal cord has an active role in the modulation and transmission of the neural signals traveling between the body and the brain. Recent advancements in functional magnetic resonance imaging (fMRI) have made the in vivo examination of spinal cord function in humans now possible. This technology has been recently extended to the investigation of resting state functional networks in the spinal cord, leading to the identification of distinct patterns of spinal cord functional connectivity. In this study, we expand on the previous work and further investigate resting state cervical spinal cord functional connectivity in healthy participants (n = 15) using high resolution imaging coupled with both seed-based functional connectivity analyses and graph theory-based metrics. Within spinal cord segment functional connectivity was present between the left and right ventral horns (bilateral motor network), left and right dorsal horns (bilateral sensory network), and the ipsilateral ventral and dorsal horns (unilateral sensory-motor network). Functional connectivity between the spinal cord segments was less apparent with the connectivity centered at the region of interest and spanning <20 mm along the superior-inferior axis. In a subset of participants (n = 10), the cervical spinal cord functional network was demonstrated to be state-dependent as thermal stimulation of the right ventrolateral forearm resulted in significant disruption of the bilateral sensory network, increased network global efficiency, and decreased network modularity.

Lateral Corticospinal Tract Damage Correlates With Motor Output in Incomplete Spinal Cord Injury

OBJECTIVE

To investigate the relationship between spinal cord damage and specific motor function in participants with incomplete spinal cord injury (iSCI).

DESIGN

Single-blinded, cross-sectional study design.

SETTING

University setting research laboratory.

PARTICIPANTS

Individuals with chronic cervical iSCI (N=14; 1 woman, 13 men; average age ± SD, 43±12y).

INTERVENTIONS

Not applicable.

MAIN OUTCOME MEASURES

Axial T2-weighted magnetic resonance imaging (MRI) of spinal cord damage was performed in 14 participants with iSCI. Each participant's damage was processed for total damage quantification, lateral corticospinal tract (LCST) and gracile fasciculus (GF) analysis. Plantarflexion and knee extension were quantified using an isokinetic dynamometer. Walking ability was assessed using a 6-minute walk test.

RESULTS

Total damage was correlated with plantarflexion, knee extension, and distance walked in 6 minutes. Right LCST damage was correlated with right plantarflexion and right knee extension, while left LCST damage was correlated with left-sided measures. Right and left GF damage was not correlated with the motor output measures.

CONCLUSIONS

MRI measures of spinal cord damage were correlated to motor function, and this measure appears to have spatial specificity to descending tracts, which may offer prognostic value after SCI.

Test-Retest and Interreader Reproducibility of Semiautomated Atlas-Based Analysis of Diffusion Tensor Imaging Data in Acute Cervical Spine Trauma in Adult Patients

BACKGROUND AND PURPOSE

DTI is a tool for microstructural spinal cord injury evaluation. This study evaluated the reproducibility of a semiautomated segmentation algorithm of spinal cord DTI.

MATERIALS AND METHODS

Forty-two consecutive patients undergoing acute trauma cervical spine MR imaging underwent 2 axial DTI scans in addition to their clinical scan. The datasets were put through a semiautomated probabilistic segmentation algorithm that selected white matter, gray matter, and 24 individual white matter tracts. Regional and white matter tract volume, fractional anisotropy, and mean diffusivity values were calculated. Two readers performed the nonautomated steps to evaluate interreader reproducibility. The coefficient of variation and intraclass correlation coefficient were used to assess test-retest and interreader reproducibility.

RESULTS

Of 42 patients, 30 had useable data. Test-retest reproducibility of fractional anisotropy was high for white matter as a whole (coefficient of variation, 3.8%; intraclass correlation coefficient, 0.93). Test-retest coefficient-of-variation ranged from 8.0%-18.2% and intraclass correlation coefficients from 0.47-0.80 across individual white matter tracts. Mean diffusivity metrics also had high test-retest reproducibility (white matter: coefficient-of-variation, 5.6%; intraclass correlation coefficient, 0.86) with coefficients of variation from 11.6%-18.3% and intraclass correlation coefficients from 0.57-0.74 across individual tracts, with better agreement for larger tracts. The coefficients of variation of fractional anisotropy and mean diffusivity both had significant negative relationships with white matter volume (26%-27% decrease for each doubling of white matter volume, P < .01).

CONCLUSIONS

DTI spinal cord segmentation is reproducible in the setting of acute spine trauma, specifically for larger white matter tracts and total white or gray matter.

Assessing structure and function of myelin in cervical spondylotic myelopathy: Evidence of demyelination

PURPOSE

To assess the extent of demyelination in cervical spondylotic myelopathy (CSM) using myelin water imaging (MWI) and electrophysiologic techniques.

METHODS

Somatosensory evoked potentials (SSEPs) and MWI were acquired in 14 patients with CSM and 18 age-matched healthy controls. MWI was performed on a 3.0T whole body magnetic resonance scanner. Myelin water fraction (MWF) was extracted for the dorsal columns and whole cord. SSEPs and MWF were also compared with conventional MRI outcomes, including T2 signal intensity, compression ratio, maximum spinal cord compression (MSCC), and maximum canal compromise (MCC).

RESULTS

Group analysis showed marked differences in T2 signal intensity, compression ratio, MSCC, and MCC between healthy controls and patients with CSM. There were no group differences in MWF and SSEP latencies. However, patients with CSM with pathologic SSEPs exhibited reduction in MWF (p < 0.05). MWF was also correlated with SSEP latencies.

CONCLUSION

Our findings provide evidence of decreased myelin content in the spinal cord associated with impaired spinal cord conduction in patients with CSM. While conventional MRI are of great value to define the extent of cord compression, they show a limited correlation with functional deficits (i.e., delayed SSEPs). MWI provides independent and complementary readouts to spinal cord compression, with a high specificity to detect impaired conduction.

Diffusion Tensor Imaging of the Spinal Cord: Clinical Value, Investigational Applications, and Technical Limitations

Although diffusion-weighted imaging (DWI) has become a mainstay in modern brain imaging, it remains less utilized in the evaluation of the spinal cord. Many studies have shown promise in using DWI and diffusion-tensor imaging (DTI) for evaluation of the spinal cord; however, application has been stalled by technical obstacles and artifacts, and questions remain regarding its clinical utility on an individual examination level. This review discusses the background, concepts, and technical aspects of DWI and DTI, specifically for imaging of the spinal cord. The clinical and investigational applications of spinal cord DTI, as well as the practical difficulties and limitations of DWI and DTI for the evaluation of the spinal cord are examined.

Neurodegeneration in the Spinal Ventral Horn Prior to Motor Impairment in Cervical Spondylotic Myelopathy

Remote gray matter pathology has been suggested rostral to the compression site in cervical spondylotic myelopathy (CSM). We therefore assessed neurodegeneration in the gray matter ventral and dorsal horns. Twenty patients with CSM and 18 healthy subjects underwent a high-resolution structural and diffusion magnetic resonance imaging protocol at vertebra C2/C3. Patients received comprehensive clinical assessments. T2*-weighted data provided cross-sectional area measurements of gray matter ventral and dorsal horns to identify atrophy. At the identical location, mean diffusivity (MD) and fractional anisotropy (FA) determined the microstructural integrity. Finally, the relationships between neurodegeneration occurring in the gray and white matter and clinical impairment were investigated. Patients suffered from mild-to-moderate CSM with mainly sensory impairment. In the ventral horns, cross-sectional area was not reduced (p = 0.863) but MD was increased (p = 0.045). The magnitude of MD changes within the ventral horn was associated with white matter diffusivity changes (MD: p = 0.013; FA: p = 0.028) within the lateral corticospinal tract. In contrast, dorsal horn cross-sectional area was reduced by 16.0% (p < 0.001) without alterations in diffusivity indices, compared with controls. No associations between the magnitude of ventral and dorsal horn neurodegeneration and clinical impairment were evident. Focal cord gray matter pathology is evident remote to the compression site in vivo in CSM patients. Microstructural changes in the ventral horns (i.e., motoneurons) related to corticospinal tract integrity in the absence of atrophy and marked motor impairment. Dorsal horn atrophy corresponded to main clinical representation of sensory impairment. Thus, neuroimaging biomarkers of cord gray matter integrity reveal focal neurodegeneration prior to marked clinical impairment and thus could serve as predictors of ensuing impairment in CSM patients.

Cervical cord myelin water imaging shows degenerative changes over one year in multiple sclerosis but not neuromyelitis optica spectrum disorder

Spinal cord pathology is a feature of both neuromyelitis optica spectrum disorder (NMOSD) and relapsing-remitting multiple sclerosis (MS). While subclinical disease activity has been described in MS using quantitative magnetic resonance imaging measures, current evidence suggests that neurodegeneration is absent between relapses in NMOSD, although most evidence comes from brain studies. We aimed to assess cross-sectional differences and longitudinal changes in myelin integrity in relapse-free MS and NMOSD subjects over one year. 15 NMOSD, 15 MS subjects, and 17 healthy controls were scanned at 3 T using a cervical cord mcDESPOT protocol. A subset of 8 NMOSD, 11 MS subjects and 14 controls completed follow-up. Measures of the myelin water fraction (f_M) within lesioned and non-lesioned cord segments were collected. At baseline, f_M in lesioned and non-lesioned segments was significantly reduced in MS (lesioned: p = 0.002; non-lesioned: p = 0.03) and NMOSD (lesioned: p = 0.0007; non-lesioned: p = 0.002) compared to controls. Longitudinally, f_M decreased within non-lesioned cord segments in the MS group (− 7.3%, p = 0.02), but not in NMOSD (+ 5.8%, p = 0.1), while change in lesioned segments f_M did not differ from controls' in either patient group. These results suggest that degenerative changes outside of lesioned areas can be observed over a short time frame in MS, but not NMOSD, and support the use of longitudinal myelin water imaging for the assessment of pathological changes in the cervical cord in demyelinating diseases.

•

MS and NMOSD subjects underwent longitudinal cervical cord myelin water imaging.

•

Reduced myelin water fraction in MS and NMOSD normal-appearing and lesioned areas

•

Decrease in myelin in normal-appearing tissue over 1 year in MS, but not NMOSD

•

Further evidence that disease progression is absent between relapses in NMOSD.

Clinically Feasible Microstructural MRI to Quantify Cervical Spinal Cord Tissue Injury Using DTI, MT, and T2*-Weighted Imaging: Assessment of Normative Data and Reliability

BACKGROUND AND PURPOSE

DTI, magnetization transfer, T2*-weighted imaging, and cross-sectional area can quantify aspects of spinal cord microstructure. However, clinical adoption remains elusive due to complex acquisitions, cumbersome analysis, limited reliability, and wide ranges of normal values. We propose a simple multiparametric protocol with automated analysis and report normative data, analysis of confounding variables, and reliability.

MATERIALS AND METHODS

Forty healthy subjects underwent T2WI, DTI, magnetization transfer, and T2*WI at 3T in <35 minutes using standard hardware and pulse sequences. Cross-sectional area, fractional anisotropy, magnetization transfer ratio, and T2*WI WM/GM signal intensity ratio were calculated. Relationships between MR imaging metrics and age, sex, height, weight, cervical cord length, and rostrocaudal level were analyzed. Test-retest coefficient of variation measured reliability in 24 DTI, 17 magnetization transfer, and 16 T2*WI datasets. DTI with and without cardiac triggering was compared in 10 subjects.

RESULTS

T2*WI WM/GM showed lower intersubject coefficient of variation (3.5%) compared with magnetization transfer ratio (5.8%), fractional anisotropy (6.0%), and cross-sectional area (12.2%). Linear correction of cross-sectional area with cervical cord length, fractional anisotropy with age, and magnetization transfer ratio with age and height led to decreased coefficients of variation (4.8%, 5.4%, and 10.2%, respectively). Acceptable reliability was achieved for all metrics/levels (test-retest coefficient of variation < 5%), with T2*WI WM/GM comparing favorably with fractional anisotropy and magnetization transfer ratio. DTI with and without cardiac triggering showed no significant differences for fractional anisotropy and test-retest coefficient of variation.

CONCLUSIONS

Reliable multiparametric assessment of spinal cord microstructure is possible by using clinically suitable methods. These results establish normalization procedures and pave the way for clinical studies, with the potential for improving diagnostics, objectively monitoring disease progression, and predicting outcomes in spinal pathologies.

Spinal cord grey matter segmentation challenge

An important image processing step in spinal cord magnetic resonance imaging is the ability to reliably and accurately segment grey and white matter for tissue specific analysis. There are several semi- or fully-automated segmentation methods for cervical cord cross-sectional area measurement with an excellent performance close or equal to the manual segmentation. However, grey matter segmentation is still challenging due to small cross-sectional size and shape, and active research is being conducted by several groups around the world in this field. Therefore a grey matter spinal cord segmentation challenge was organised to test different capabilities of various methods using the same multi-centre and multi-vendor dataset acquired with distinct 3D gradient-echo sequences. This challenge aimed to characterize the state-of-the-art in the field as well as identifying new opportunities for future improvements. Six different spinal cord grey matter segmentation methods developed independently by various research groups across the world and their performance were compared to manual segmentation outcomes, the present gold-standard. All algorithms provided good overall results for detecting the grey matter butterfly, albeit with variable performance in certain quality-of-segmentation metrics. The data have been made publicly available and the challenge web site remains open to new submissions. No modifications were introduced to any of the presented methods as a result of this challenge for the purposes of this publication.

A Novel MRI Biomarker of Spinal Cord White Matter Injury: T2*-Weighted White Matter to Gray Matter Signal Intensity Ratio

BACKGROUND AND PURPOSE

T2*-weighted imaging provides sharp contrast between spinal cord GM and WM, allowing their segmentation and cross-sectional area measurement. Injured WM demonstrates T2*WI hyperintensity but requires normalization for quantitative use. We introduce T2*WI WM/GM signal-intensity ratio and compare it against cross-sectional area, the DTI metric fractional anisotropy, and magnetization transfer ratio in degenerative cervical myelopathy.

MATERIALS AND METHODS

Fifty-eight patients with degenerative cervical myelopathy and 40 healthy subjects underwent 3T MR imaging, covering C1-C7. Metrics were automatically extracted at maximally compressed and uncompressed rostral/caudal levels. Normalized metrics were compared with t tests, area under the curve, and logistic regression. Relationships with clinical measures were analyzed by using Pearson correlation and multiple linear regression.

RESULTS

The maximally compressed level cross-sectional area demonstrated superior differences (P = 1 × 10^-13), diagnostic accuracy (area under the curve = 0.890), and univariate correlation with the modified Japanese Orthopedic Association score (0.66). T2*WI WM/GM showed strong differences (rostral: P = 8 × 10^-7; maximally compressed level: P = 1 × 10^-11; caudal: P = 1 × 10^-4), correlations (modified Japanese Orthopedic Association score; rostral: -0.52; maximally compressed level: -0.59; caudal: -0.36), and diagnostic accuracy (rostral: 0.775; maximally compressed level: 0.860; caudal: 0.721), outperforming fractional anisotropy and magnetization transfer ratio in most comparisons and cross-sectional area at rostral/caudal levels. Rostral T2*WI WM/GM showed the strongest correlations with focal motor (-0.45) and sensory (-0.49) deficits and was the strongest independent predictor of the modified Japanese Orthopedic Association score (P = .01) and diagnosis (P = .02) in multivariate models (R² = 0.59, P = 8 × 10^-13; area under the curve = 0.954, respectively).

CONCLUSIONS

T2*WI WM/GM shows promise as a novel biomarker of WM injury. It detects damage in compressed and uncompressed regions and contributes substantially to multivariate models for diagnosis and correlation with impairment. Our multiparametric approach overcomes limitations of individual measures, having the potential to improve diagnostics, monitor progression, and predict outcomes.

Topologically preserving straightening of spinal cord MRI

PURPOSE

To propose a robust and accurate method for straightening magnetic resonance (MR) images of the spinal cord, based on spinal cord segmentation, that preserves spinal cord topology and that works for any MRI contrast, in a context of spinal cord template-based analysis.

MATERIALS AND METHODS

The spinal cord curvature was computed using an iterative Non-Uniform Rational B-Spline (NURBS) approximation. Forward and inverse deformation fields for straightening were computed by solving analytically the straightening equations for each image voxel. Computational speed-up was accomplished by solving all voxel equation systems as one single system. Straightening accuracy (mean and maximum distance from straight line), computational time, and robustness to spinal cord length was evaluated using the proposed and the standard straightening method (label-based spline deformation) on 3T T₂ - and T₁ -weighted images from 57 healthy subjects and 33 patients with spinal cord compression due to degenerative cervical myelopathy (DCM).

RESULTS

The proposed algorithm was more accurate, more robust, and faster than the standard method (mean distance = 0.80 vs. 0.83 mm, maximum distance = 1.49 vs. 1.78 mm, time = 71 vs. 174 sec for the healthy population and mean distance = 0.65 vs. 0.68 mm, maximum distance = 1.28 vs. 1.55 mm, time = 32 vs. 60 sec for the DCM population).

CONCLUSION

A novel image straightening method that enables template-based analysis of quantitative spinal cord MRI data is introduced. This algorithm works for any MRI contrast and was validated on healthy and patient populations. The presented method is implemented in the Spinal Cord Toolbox, an open-source software for processing spinal cord MRI data.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2017;46:1209-1219.

g-Ratio weighted imaging of the human spinal cord in vivo

The fiber g-ratio is defined as the ratio of the inner to the outer diameter of the myelin sheath. This ratio provides a measure of the myelin thickness that complements axon morphology (diameter and density) for assessment of demyelination in diseases such as multiple sclerosis. Previous work has shown that an aggregate g-ratio map can be computed using a formula that combines axon and myelin density measured with quantitative MRI. In this work, we computed g-ratio weighted maps in the cervical spinal cord of nine healthy subjects. We utilized the 300mT/m gradients from the CONNECTOM scanner to estimate the fraction of restricted water (fr) with high accuracy, using the CHARMED model. Myelin density was estimated using the lipid and macromolecular tissue volume (MTV) method, derived from normalized proton density (PD) mapping. The variability across spinal level, laterality and subject were assessed using a three-way ANOVA. The average g-ratio value obtained in the white matter was 0.76+/-0.03, consistent with previous histology work. Coefficients of variation of fr and MTV were respectively 4.3% and 13.7%. fr and myelin density were significantly different across spinal tracts (p=3×10-7 and 0.004 respectively) and were positively correlated in the white matter (r=0.42), suggesting shared microstructural information. The aggregate g-ratio did not show significant differences across tracts (p=0.6). This study suggests that fr and myelin density can be measured in vivo with high precision and that they can be combined to produce a g-ratio-weighted map robust to free water pool contamination from cerebrospinal fluid or veins. Potential applications include the study of early demyelination in multiple sclerosis, and the quantitative assessment of remyelination drugs.

Investigating resting-state functional connectivity in the cervical spinal cord at 3T

The study of spontaneous fluctuations in the blood-oxygen-level-dependent (BOLD) signal has recently been extended from the brain to the spinal cord. Two ultra-high field functional magnetic resonance imaging (fMRI) studies in humans have provided evidence for reproducible resting-state connectivity between the dorsal horns as well as between the ventral horns, and a study in non-human primates has shown that these resting-state signals are impacted by spinal cord injury. As these studies were carried out at ultra-high field strengths using region-of-interest (ROI) based analyses, we investigated whether such resting-state signals could also be observed at the clinically more prevalent field strength of 3 T. In a reanalysis of a sample of 20 healthy human participants who underwent a resting-state fMRI acquisition of the cervical spinal cord, we were able to observe significant dorsal horn connectivity as well as ventral horn connectivity, but no consistent effects for connectivity between dorsal and ventral horns, thus replicating the human 7 T results. These effects were not only observable when averaging along the acquired length of the spinal cord, but also when we examined each of the acquired spinal segments separately, which showed similar patterns of connectivity. Finally, we investigated the robustness of these resting-state signals against variations in the analysis pipeline by varying the type of ROI creation, temporal filtering, nuisance regression and connectivity metric. We observed that – apart from the effects of band-pass filtering – ventral horn connectivity showed excellent robustness, whereas dorsal horn connectivity showed moderate robustness. Together, our results provide evidence that spinal cord resting-state connectivity is a robust and spatially consistent phenomenon that could be a valuable tool for investigating the effects of pathology, disease progression, and treatment response in neurological conditions with a spinal component, such as spinal cord injury.