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

Multiparametric in vivo analyses of the brain and spine identify structural and metabolic biomarkers in men with adrenomyeloneuropathy

OBJECTIVE

Progressive myelopathy causes severe handicap in men with adrenomyeloneuropathy (AMN), an X-linked disorder due to ABCD1 pathogenic variants. At present, treatments are symptomatic but disease-modifying therapies are under evaluation. Given the small effect size of clinical scales in AMN, biomarkers with higher effect size are needed. Here we used high-resolution magnetic resonance techniques to identify non-invasive in vivo biomarkers of the brain and spine with high effect sizes.

METHODS

We performed a multiparametric imaging and spectroscopy study in 23 male patients with AMN (age: 44 ± 11) and 23 male controls (age: 43 ± 11) of similar age and body-mass index. We combined (i) macrostructural analyses of the spine, using cross-sectional area (CSA) and magnetization transfer ratio (MTR), (ii) microstructural analyses of the spine and the brain, using diffusion tensor and the newly developed fixel-based analysis, and (iii) advanced metabolic analyses of the spine using metabolite cycling coupled to a semi-LASER sequences.

RESULTS

Macrostructural alterations (decrease in CSA and MTR) were observed in patients at all spinal cord levels studied (C1-T2 for CSA and C1-C5 for MTR) (p < 0.001). Microstructural alterations were observed in the spine and brain on diffusion tensor and fixel-based metrics though the latter showed higher effect sizes. Metabolic alterations were observed in patients as a decreased total N-acetylaspartate/myo-inositol ratio (p < 0.001). Overall, MTR showed the highest effect size.

CONCLUSION

This cross-sectional study supports the use of multiparametric techniques that elucidate the structural, microstructural and metabolic alterations in AMN. These outcome measures should be tested longitudinally and in clinical trials.

Mapping the rest of the human connectome: Atlasing the spinal cord and peripheral nervous system

The emergence of diffusion, structural, and functional neuroimaging methods has enabled major multi-site efforts to map the human connectome, which has heretofore been defined as containing all neural connections in the central nervous system (CNS). However, these efforts are not structured to examine the richness and complexity of the peripheral nervous system (PNS), which arguably forms the (neglected) rest of the connectome. Despite increasing interest in an atlas of the spinal cord (SC) and PNS which is simultaneously stereotactic, interactive, electronically dissectible, scalable, population-based and deformable, little attention has thus far been devoted to this task of critical importance. Nevertheless, the atlasing of these complete neural structures is essential for neurosurgical planning, neurological localization, and for mapping those components of the human connectome located outside of the CNS. Here we recommend a modification to the definition of the human connectome to include the SC and PNS, and argue for the creation of an inclusive atlas to complement current efforts to map the brain's human connectome, to enhance clinical education, and to assist progress in neuroscience research. In addition to providing a critical overview of existing neuroimaging techniques, image processing methodologies and algorithmic advances which can be combined for the creation of a full connectome atlas, we outline a blueprint for ultimately mapping the entire human nervous system and, thereby, for filling a critical gap in our scientific knowledge of neural connectivity.

Spatial correspondence of spinal cord white matter tracts using diffusion tensor imaging, fibre tractography, and atlas-based segmentation

PURPOSE

Neuroimaging provides great utility in complex spinal surgeries, particularly when anatomical geometry is distorted by pathology (tumour, degeneration, etc.). Spinal cord MRI diffusion tractography can be used to generate streamlines; however, it is unclear how well they correspond with white matter tract locations along the cord microstructure. The goal of this work was to evaluate the spatial correspondence of DTI tractography with anatomical MRI in healthy anatomy (where anatomical locations can be well defined in T1-weighted images).

METHODS

Ten healthy volunteers were scanned on a 3T system. T1-weighted (1 × 1 × 1 mm) and diffusion-weighted images (EPI readout, 2 × 2 × 2 mm, 30 gradient directions) were acquired and subsequently registered (Spinal Cord Toolbox (SCT)). Atlas-based (SCT) anatomic label maps of the left and right lateral corticospinal tracts were identified for each vertebral region (C2-C6) from T1 images. Tractography streamlines were generated with a customized approach, enabling seeding of specific spinal tract regions corresponding to individual vertebral levels. Spatial correspondence of generated fibre streamlines with anatomic tract segmentations was compared in unseeded regions of interest (ROIs).

RESULTS

Spatial correspondence of the lateral corticospinal tract streamlines was good over a single vertebral ROI (Dice's similarity coefficient (DSC) = 0.75 ± 0.08, Hausdorff distance = 1.08 ± 0.17 mm). Over larger ROI, fair agreement between tractography and anatomical labels was achieved (two levels: DSC = 0.67 ± 0.13, three levels: DSC = 0.52 ± 0.19).

CONCLUSION

DTI tractography produced good spatial correspondence with anatomic white matter tracts, superior to the agreement between multiple manual tract segmentations (DSC ~ 0.5). This supports further development of spinal cord tractography for computer-assisted neurosurgery.

Comparison of multi echo T2 relaxation and steady state approaches for myelin imaging in the central nervous system

The traditional approach for measuring myelin-associated water with quantitative magnetic resonance imaging (MRI) uses multi-echo T₂ relaxation data to calculate the myelin water fraction (MWF). A fundamentally different approach, abbreviated “mcDESPOT”, uses a more efficient steady-state acquisition to generate an equivalent metric (f_M). Although previous studies have demonstrated inherent instability and bias in the complex mcDESPOT analysis procedure, f_M has often been used as a surrogate for MWF. We produced and compared multivariate atlases of MWF and f_M in healthy human brain and cervical spinal cord (available online) and compared their ability to detect multiple sclerosis pathology. A significant bias was found in all regions (p < 10^–5), albeit reversed for spinal cord (f_M-MWF =  − 3.4%) compared to brain (+ 6.2%). MWF and f_M followed an approximately linear relationship for regions with MWF <  ~ 10%. For MWF >  ~ 10%, the relationship broke down and f_M no longer increased in tandem with MWF. For multiple sclerosis patients, MWF and f_M Z score maps showed overlapping areas of low Z score and similar trends between patients and brain regions, although those of f_M generally had greater spatial extent and magnitude of severity. These results will guide future choice of myelin-sensitive quantitative MRI and improve interpretation of studies using either myelin imaging approach.

Development of new outcome measures for adult SMA type III and IV: a multimodal longitudinal study

OBJECTIVE

The aim of this study was the comprehensive characterisation of longitudinal clinical, electrophysiological and neuroimaging measures in type III and IV adult spinal muscular atrophy (SMA) with a view to propose objective monitoring markers for future clinical trials.

METHODS

Fourteen type III or IV SMA patients underwent standardised assessments including muscle strength testing, functional evaluation (SMAFRS and MFM), MUNIX (abductor pollicis brevis, APB; abductor digiti minimi, ADM; deltoid; tibialis anterior, TA; trapezius) and quantitative cervical spinal cord MRI to appraise segmental grey and white matter atrophy. Patients underwent a follow-up assessment with the same protocol 24 months later. Longitudinal comparisons were conducted using the Wilcoxon-test for matched data. Responsiveness was estimated using standardized response means (SRM) and a composite score was generated based on the three most significant variables.

RESULTS

Significant functional decline was observed based on SMAFRS (p = 0.019), pinch and knee flexion strength (p = 0.030 and 0.027), MUNIX and MUSIX value in the ADM (p = 0.0006 and 0.043) and in TA muscle (p = 0.025). No significant differences were observed based on cervical MRI measures. A significant reduction was detected in the composite score (p = 0.0005, SRM = -1.52), which was the most responsive variable and required a smaller number of patients than single variables in the estimation of sample size for clinical trials.

CONCLUSIONS

Quantitative strength testing, SMAFRS and MUNIX readily capture disease progression in adult SMA patients. Composite multimodal scores increase predictive value and may reduce sample size requirements in clinical trials.

Minimum detectable spinal cord atrophy with automatic segmentation: Investigations using an open-access dataset of healthy participants

•

Evaluate the robustness of an automated analysis pipeline for detecting SC atrophy.

•

Simulate spinal cord atrophy and scan-rescan variability.

•

Fully automated analysis method available on an open access database.

•

Evaluation of sample size and inter/intra-subject variability for T1w and T2w images.

Spinal cord atrophy is a well-known biomarker in multiple sclerosis (MS) and other diseases. It is measured by segmenting the spinal cord on an MRI image and computing the average cross-sectional area (CSA) over a few slices. Introduced about 25 years ago, this procedure is highly sensitive to the quality of the segmentation and is prone to rater-bias. Recently, fully-automated spinal cord segmentation methods, which remove the rater-bias and enable the automated analysis of large populations, have been introduced. A lingering question related to these automated methods is: How reliable are they at detecting atrophy? In this study, we evaluated the precision and accuracy of automated atrophy measurements by simulating scan-rescan experiments.

Spinal cord MRI data from the open-access spine-generic project were used. The dataset aggregates 42 sites worldwide and consists of 260 healthy subjects and includes T1w and T2w contrasts. To simulate atrophy, each volume was globally rescaled at various scaling factors. Moreover, to simulate patient repositioning, random rigid transformations were applied. Using the DeepSeg algorithm from the Spinal Cord Toolbox, the spinal cord was segmented and vertebral levels were identified. Then, the average CSA between C3-C5 vertebral levels was computed for each Monte Carlo sample, allowing us to derive measures of atrophy, intra/inter-subject variability, and sample-size calculations.

The minimum sample size required to detect an atrophy of 2% between unpaired study arms, commonly seen in MS studies, was 467 +/− 13.9 using T1w and 467 +/− 3.2 using T2w images. The minimum sample size to detect a longitudinal atrophy (between paired study arms) of 0.8% was 60 +/− 25.1 using T1w and 10 +/− 1.2 using T2w images. At the intra-subject level, the estimated CSA, observed in this study, showed good precision compared to other studies with COVs (across Monte Carlo transformations) of 0.8% for T1w and 0.6% for T2w images.

While these sample sizes seem small, we would like to stress that these results correspond to a “best case” scenario, in that the dataset used here was of particularly good quality and the model for simulating atrophy does not encompass all the variability met in real-life datasets. The simulated atrophy and scan-rescan variability may over-simplify the biological reality. The proposed framework is open-source and available at https://csa-atrophy.readthedocs.io/.

Multiple sclerosis lesions in motor tracts from brain to cervical cord: spatial distribution and correlation with disability

Despite important efforts to solve the clinico-radiological paradox, correlation between lesion load and physical disability in patients with multiple sclerosis remains modest. One hypothesis could be that lesion location in corticospinal tracts plays a key role in explaining motor impairment. In this study, we describe the distribution of lesions along the corticospinal tracts from the cortex to the cervical spinal cord in patients with various disease phenotypes and disability status. We also assess the link between lesion load and location within corticospinal tracts, and disability at baseline and 2-year follow-up. We retrospectively included 290 patients (22 clinically isolated syndrome, 198 relapsing remitting, 39 secondary progressive, 31 primary progressive multiple sclerosis) from eight sites. Lesions were segmented on both brain (T2-FLAIR or T2-weighted) and cervical (axial T2- or T2*-weighted) MRI scans. Data were processed using an automated and publicly available pipeline. Brain, brainstem and spinal cord portions of the corticospinal tracts were identified using probabilistic atlases to measure the lesion volume fraction. Lesion frequency maps were produced for each phenotype and disability scores assessed with Expanded Disability Status Scale score and pyramidal functional system score. Results show that lesions were not homogeneously distributed along the corticospinal tracts, with the highest lesion frequency in the corona radiata and between C2 and C4 vertebral levels. The lesion volume fraction in the corticospinal tracts was higher in secondary and primary progressive patients (mean = 3.6 ± 2.7% and 2.9 ± 2.4%), compared to relapsing-remitting patients (1.6 ± 2.1%, both P < 0.0001). Voxel-wise analyses confirmed that lesion frequency was higher in progressive compared to relapsing-remitting patients, with significant bilateral clusters in the spinal cord corticospinal tracts (P < 0.01). The baseline Expanded Disability Status Scale score was associated with lesion volume fraction within the brain (r = 0.31, P < 0.0001), brainstem (r = 0.45, P < 0.0001) and spinal cord (r = 0.57, P < 0.0001) corticospinal tracts. The spinal cord corticospinal tracts lesion volume fraction remained the strongest factor in the multiple linear regression model, independently from cord atrophy. Baseline spinal cord corticospinal tracts lesion volume fraction was also associated with disability progression at 2-year follow-up (P = 0.003). Our results suggest a cumulative effect of lesions within the corticospinal tracts along the brain, brainstem and spinal cord portions to explain physical disability in multiple sclerosis patients, with a predominant impact of intramedullary lesions.

In vivo detection of microstructural spinal cord lesions in dogs with degenerative myelopathy using diffusion tensor imaging

Background

Degenerative myelopathy (DM) in dogs is a progressive neurodegenerative condition that causes white matter spinal cord lesions. These lesions are undetectable on standard magnetic resonance imaging (MRI), limiting diagnosis and monitoring of the disease. Spinal cord lesions cause disruption to the structural integrity of the axons causing water diffusion to become more random and less anisotropic. These changes are detectable by the technique of diffusion tensor imaging (DTI) which is highly sensitive to diffusion alterations secondary to white matter lesion development.

Objective

Perform spinal DTI on cohorts of dogs with and without DM to identify if lesions caused by DM will cause a detectable alteration in spinal cord diffusivity that correlates with neurological status.

Animals

Thirteen dogs with DM and 13 aged‐matched controls.

Methods

All animals underwent MRI with DTI of the entire spine. Diffusivity parameters fractional anisotropy (FA) and mean diffusivity (MD) were measured at each vertebral level and statistically compared between groups.

Results

Dogs with DM had significant decreases in FA within the regions of the spinal cord that had high expected lesion load. Decreases in FA were most significant in dogs with severe forms of the disease and correlated with neurological grade.

Conclusions and Clinical Importance

Findings suggest that FA has the potential to be a biomarker for spinal cord lesion development in DM and could play an important role in improving diagnosis and monitoring of this condition.

Macromolecular changes in spinal cord white matter characterize whiplash outcome at 1-year post motor vehicle collision

Each year, whiplash injuries from motor vehicle collisions (MVC) affect millions worldwide, with no strong evidence of pathology. While the majority recover soon after the injury, the same is not true for roughly 20% reporting higher levels of pain and distress, without diagnostic options. This study used magnetization transfer (MT) imaging to quantify white matter integrity in 78 subjects with varying levels of pain, 1 year after MVC. MT images of the cervical spinal cord were collected parallel to the intervertebral disks. MT ratios (MTR) were calculated in select white matter tracts along with MTR homogeneity (MTRh) at each level. Significant differences were observed between clinical outcome groups in the left and right spinothalamic tracts (p = 0.003 and 0.020) and MTRh (p = 0.009). MTRh was elevated in females with poor recovery versus females reporting recovery (p < 0.001) or milder symptoms (p < 0.001), and in males reporting recovery (p = 0.007) or no recovery (p < 0.001). There was a significant interaction between recovery status and sex (p = 0.015). MT imaging identified tract specific and regional changes in white matter integrity suggesting potential insults to the cord. Additionally, significant MTRh differences between sexes were observed, characterizing the heterogeneity of whiplash recovery and worse outcomes in females.

A comparison of the effectiveness of functional MRI analysis methods for pain research: The new normal

Studies of the neural basis of human pain processing present many challenges because of the subjective and variable nature of pain, and the inaccessibility of the central nervous system. Neuroimaging methods, such as functional magnetic resonance imaging (fMRI), have provided the ability to investigate these neural processes, and yet commonly used analysis methods may not be optimally adapted for studies of pain. Here we present a comparison of model-driven and data-driven analysis methods, specifically for the study of human pain processing. Methods are tested using data from healthy control participants in two previous studies, with separate data sets spanning the brain, and the brainstem and spinal cord. Data are analyzed by fitting time-series responses to predicted BOLD responses in order to identify significantly responding regions (model-driven), as well as with connectivity analyses (data-driven) based on temporal correlations between responses in spatially separated regions, and with connectivity analyses based on structural equation modeling, allowing for multiple source regions to explain the signal variations in each target region. The results are assessed in terms of the amount of signal variance that can be explained in each region, and in terms of the regions and connections that are identified as having BOLD responses of interest. The characteristics of BOLD responses in identified regions are also investigated. The results demonstrate that data-driven approaches are more effective than model-driven approaches for fMRI studies of pain.

What are the gray and white matter volumes of the human spinal cord?

The gray matter of the spinal cord is the seat of somata of various types of neurons devoted to the sensory and motor activities of the limbs and trunk as well as a part of the autonomic nervous system. The volume of the spinal gray matter is an indicator of the local neuronal processing, and this can decrease due to atrophy associated with degenerative diseases and injury. Nevertheless, the absolute volume of the human spinal cord has rarely been reported, if ever. Here, we use high-resolution magnetic resonance imaging, with a cross-sectional resolution of 50 × 50 μm and a voxel size of 0.0005 mm³ to estimate the total gray and white matter volume of a post mortem human female spinal cord. Segregation of gray and white matter was accomplished using deep learning image segmentation. Furthermore, we include data from a male spinal cord of a previously published study. The gray and white matter volumes were found to be 2.87 and 11.33 mL, respectively, for the female and 3.55 and 19.33 mL, respectively, for the male. The gray and white matter profiles along the vertebral axis were found to be strikingly similar, and the volumes of the cervical, thoracic, and lumbosacral sections were almost equal.NEW & NOTEWORTHY Here, we combine high-field MRI (9.4 T) and deep learning for a post mortem reconstruction of the gray and white matter in human spinal cords. We report a minuscule total gray matter volume of 2.87 mL for a female and 3.55 mL for a male. For comparison, these volumes correspond approximately to the distal digit of the little finger.

Injury volume extracted from MRI predicts neurologic outcome in acute spinal cord injury: A prospective TRACK-SCI pilot study

Conventional MRI measures of traumatic spinal cord injury severity largely rely on 2-dimensional injury characteristics such as intramedullary lesion length and cord compression. Recent advances in spinal cord (SC) analysis have led to the development of a robust anatomic atlas incorporated into an open-source platform called the Spinal Cord Toolbox (SCT) that allows for quantitative volumetric injury analysis. In the current study, we evaluate the prognostic value of volumetric measures of spinal cord injury on MRI following registration of T2-weighted (T2w) images and segmented lesions from acute SCI patients with a standardized atlas. This IRB-approved prospective cohort study involved the image analysis of 60 blunt cervical SCI patients enrolled in the TRACK-SCI clinical research protocol. Axial T2w MRI data obtained within 24 h of injury were processed using the SCT. Briefly, SC MRIs were automatically segmented using the sct_deepseg_sc tool in the SCT and segmentations were manually corrected by a neuro-radiologist. Lesion volume data were used as predictor variables for correlation with lower extremity motor scores at discharge. Volumetric MRI measures of T2w signal abnormality comprising the SCI lesion accurately predict lower extremity motor scores at time of patient discharge. Similarly, MRI measures of injury volume significantly correlated with motor scores to a greater degree than conventional 2-D metrics of lesion size. The volume of total injury and of injured spinal cord motor regions on T2w MRI is significantly and independently associated with neurologic outcome at discharge after injury.

Sensitivity of the Inhomogeneous Magnetization Transfer Imaging Technique to Spinal Cord Damage in Multiple Sclerosis

BACKGROUND AND PURPOSE

The inhomogeneous magnetization transfer technique has demonstrated high specificity for myelin, and has shown sensitivity to multiple sclerosis-related impairment in brain tissue. Our aim was to investigate its sensitivity to spinal cord impairment in MS relative to more established MR imaging techniques (volumetry, magnetization transfer, DTI).

MATERIALS AND METHODS

Anatomic images covering the cervical spinal cord from the C1 to C6 levels and DTI, magnetization transfer/inhomogeneous magnetization transfer images at the C2/C5 levels were acquired in 19 patients with MS and 19 paired healthy controls. Anatomic images were segmented in spinal cord GM and WM, both manually and using the AMU₄₀ atlases. MS lesions were manually delineated. MR metrics were analyzed within normal-appearing and lesion regions in anterolateral and posterolateral WM and compared using Wilcoxon rank tests and z scores. Correlations between MR metrics and clinical scores in patients with MS were evaluated using the Spearman rank correlation.

RESULTS

AMU₄₀-based C1-to-C6 GM/WM automatic segmentations in patients with MS were evaluated relative to manual delineation. Mean Dice coefficients were 0.75/0.89, respectively. All MR metrics (WM/GM cross-sectional areas, normal-appearing and lesion diffusivities, and magnetization transfer/inhomogeneous magnetization transfer ratios) were observed altered in patients compared with controls (P < .05). Additionally, the absolute inhomogeneous magnetization transfer ratio z scores were significantly higher than those of the other MR metrics (P < .0001), suggesting a higher inhomogeneous magnetization transfer sensitivity toward spinal cord impairment in MS. Significant correlations with the Expanded Disability Status Scale (ρ = -0.73/P = .02, ρ = -0.81/P = .004) and the total Medical Research Council scale (ρ = 0.80/P = .009, ρ = -0.74/P = .02) were observed for inhomogeneous magnetization transfer and magnetization transfer ratio z scores, respectively, in normal-appearing WM regions, while weaker and nonsignificant correlations were obtained for DTI metrics.

CONCLUSIONS

With inhomogeneous magnetization transfer being highly sensitive to spinal cord damage in MS compared with conventional magnetization transfer and DTI, it could generate great clinical interest for longitudinal follow-up and potential remyelinating clinical trials. In line with other advanced myelin techniques with which it could be compared, it opens perspectives for multicentric investigations.

Effect of Physiological Noise on Thoracolumbar Spinal Cord Functional Magnetic Resonance Imaging in 3T Magnetic Field

Introduction:

Functional Magnetic Resonance Imaging (fMRI) methods have been used to study sensorimotor processing in the spinal cord. However, these techniques confront unwanted noises to the measured signal from the physiological fluctuations. In the spinal cord imaging, most of the challenges are consequences of cardiac and respiratory movement artifacts that are considered as significant sources of noise, especially in the thoracolumbar region. In this study, we investigated the effect of each source of physiological noise and their contribution to the outcome of the analysis of the blood-oxygen-level-dependent signal in the human thoracolumbar spinal cord.

Methods:

Fifteen young healthy male volunteers participated in the study, and pain stimuli were delivered on the L5 dermatome between the two malleoli. Respiratory and cardiac signals were recorded during the imaging session, and the generated respiration and cardiac regressors were included in the general linear model for quantification of the effect of each of them on the task-analysis results. The sum of active voxels of the clusters was calculated in the spinal cord in three correction states (respiration correction only, cardiac correction only, and respiration and cardiac noise corrections) and analyzed with analysis of variance statistical test and receiver operating characteristic curve.

Results:

The results illustrated that cardiac noise correction had an effective role in increasing the active voxels (Mean±SD = 23.46±9.46) compared to other noise correction methods. Cardiac effects were higher than other physiological noise sources

Conclusion:

In summary, our results indicate great respiration effects on the lumbar and thoracolumbar spinal cord fMRI, and its contribution to the heartbeat effect can be a significant variable in the individual fMRI data analysis. Displacement of the spinal cord and the effects of this noise in the thoracolumbar and lumbar spinal cord fMRI results are significant and cannot be ignored.

Spinal Cord Imaging Markers and Recovery of Volitional Leg Movement With Spinal Cord Epidural Stimulation in Individuals With Clinically Motor Complete Spinal Cord Injury

Previous studies have shown that epidural stimulation of the lumbosacral spinal cord (scES) can re-enable lower limb volitional motor control in individuals with chronic, clinically motor complete spinal cord injury (SCI). This observation entails that residual supraspinal connectivity to the lumbosacral spinal circuitry still persisted after SCI, although it was non-detectable when scES was not provided. In the present study, we aimed at exploring further the mechanisms underlying scES-promoted recovery of volitional lower limb motor control by investigating neuroimaging markers at the spinal cord lesion site via magnetic resonance imaging (MRI). Spinal cord MRI was collected prior to epidural stimulator implantation in 13 individuals with chronic, clinically motor complete SCI, and the spared tissue of specific regions of the spinal cord (anterior, posterior, right, left, and total cord) was assessed. After epidural stimulator implantation, and prior to any training, volitional motor control was evaluated during left and right lower limb flexion and ankle dorsiflexion attempts. The ability to generate force exertion and movement was not correlated to any neuroimaging marker. On the other hand, spared tissue of specific cord regions significantly and importantly correlated with some aspects of motor control that include activation amplitude of antagonist (negative correlation) muscles during left ankle dorsiflexion, and electromyographic coordination patterns during right lower limb flexion. The fact that amount and location of spared spinal cord tissue at the lesion site were not related to the ability to generate volitional lower limb movements may suggest that supraspinal inputs through spared spinal cord regions that differ across individuals can result in the generation of lower limb volitional motor output prior to any training when epidural stimulation is provided.

HARDI-ZOOMit protocol improves specificity to microstructural changes in presymptomatic myelopathy

Diffusion magnetic resonance imaging (dMRI) proved promising in patients with non-myelopathic degenerative cervical cord compression (NMDCCC), i.e., without clinically manifested myelopathy. Aim of the study is to present a fast multi-shell HARDI-ZOOMit dMRI protocol and validate its usability to detect microstructural myelopathy in NMDCCC patients. In 7 young healthy volunteers, 13 age-comparable healthy controls, 18 patients with mild NMDCCC and 15 patients with severe NMDCCC, the protocol provided higher signal-to-noise ratio, enhanced visualization of white/gray matter structures in microstructural maps, improved dMRI metric reproducibility, preserved sensitivity (SE = 87.88%) and increased specificity (SP = 92.31%) of control-patient group differences when compared to DTI-RESOLVE protocol (SE = 87.88%, SP = 76.92%). Of the 56 tested microstructural parameters, HARDI-ZOOMit yielded significant patient-control differences in 19 parameters, whereas in DTI-RESOLVE data, differences were observed in 10 parameters, with mostly lower robustness. Novel marker the white-gray matter diffusivity gradient demonstrated the highest separation. HARDI-ZOOMit protocol detected larger number of crossing fibers (5–15% of voxels) with physiologically plausible orientations than DTI-RESOLVE protocol (0–8% of voxels). Crossings were detected in areas of dorsal horns and anterior white commissure. HARDI-ZOOMit protocol proved to be a sensitive and practical tool for clinical quantitative spinal cord imaging.

Cortico-spinal imaging to study pain

Functional magnetic resonance imaging of the brain has helped to reveal mechanisms of pain perception in health and disease. Recently, imaging approaches have been developed that allow recording neural activity simultaneously in the brain and in the spinal cord. These approaches offer the possibility to examine pain perception in the entire central pain system and in addition, to investigate cortico-spinal interactions during pain processing. Although cortico-spinal imaging is a promising technique, it bears challenges concerning data acquisition and data analysis strategies. In this review, we discuss studies that applied simultaneous imaging of the brain and spinal cord to explore central pain processing. Furthermore, we describe different MR-related acquisition techniques, summarize advantages and disadvantages of approaches that have been implemented so far and present software that has been specifically developed for the analysis of spinal fMRI data to address challenges of spinal data analysis.

Axial MRI biomarkers of spinal cord damage to predict future walking and motor function: a retrospective study

Study design

Retrospective.

Objectives

Primary: to assess if axial damage ratios are predictors of future walking after spinal cord injury (SCI), and if they add any predictive value if initial neurological impairment grades are available. Secondary: to determine if lateral spinal cord regions are predictors of future lower extremity motor scores (LEMS).

Setting

University/hospital.

Methods

Axial T₂-weighted MRIs were used. Axial damage ratios and non-damaged lateral cord volumes were calculated. Each participant answered at 1 year after SCI, “Are you able to walk for 150 feet? (45.72 meters)” For the secondary aim, right and left LEMS were used.

Results

In total, 145 participants were selected. Individuals that could walk had smaller ratios than those that were unable. Walking and axial damage ratios were negatively correlated. A 0.374 ratio cut-off showed optimal sensitivity/specificity. When initial neurological grades were used, axial damage ratios did not add predictive value. Forty-two participants had LEMS available and were included for the secondary aim. Right cord regions and right LEMS were positively correlated and left regions and left LEMS, but these variables were also correlated with each other.

Conclusions

Axial damage ratios were significant predictors of walking ability 1 year after SCI. However, this measure did not add predictive value over initial neurological grades. Lateral cord regions correlated with same-side LEMS, but the opposite was also found, calling this biomarker’s specificity into question. Axial damage ratios may be useful in predicting walking after SCI if initial neurological grades are unavailable.

Sponsorship

This research was funded by a National Institutes of Health award, National Institute of Child Health and Development—NIH R03HD094577.

Magnetization transfer and diffusion tensor imaging in dogs with intervertebral disk herniation

Background

Quantitative magnetic resonance imaging (QMRI) techniques of magnetization transfer ratio (MTR) and diffusion tensor imaging (DTI) provide microstructural information about the spinal cord.

Objective

Compare neurologic grades using the modified Frankel scale with MTR and DTI measurements in dogs with thoracolumbar intervertebral disk herniation (IVDH).

Animals

Fifty‐one dogs with thoracolumbar IVDH.

Methods

Prospective cohort study. Quantitative MRI measurements of the spinal cord were obtained at the region of compression. A linear regression generalized estimating equations model was used to compare QMRI measurements between different neurological grades after adjusting for age, weight, duration of clinical signs, and lesion location.

Results

Grade 5 (.79  ×  10⁻³ mm²/s [median], .43−.91 [range]) and axial (1.47 × 10⁻³ mm²/s, .58−1.8) diffusivity were lower compared to grades 2 (1.003, .68−1.36; P = .02 and 1.81 × 10⁻³ mm²/s, 1.36−2.12; P < .001, respectively) and 3 (1.07 × 10⁻³ mm²/s, .77−1.5; P = .04 and 1.92 × 10⁻³ mm²/s, 1.83−2.37;P < .001, respectively). Compared to dogs with acute myelopathy, chronic myelopathy was associated with higher mean (1.02 × 10⁻³ mm²/s, .77−1.36 vs. .83 × 10⁻³ mm²/s, .64−1.5; P = .03) and radial diffusivity (.75 × 10⁻³ mm²/s, .38−1.04 vs. .44 × 10⁻³ mm²/s, .22−1.01; P = .008) and lower MTR (46.76, 31.8−56.43 vs. 54.4, 45.2−62.27; P = .004) and fractional anisotropy (.58, .4−0.75 vs. .7, .46−.85; P = .02). Fractional anisotropy was lower in dogs with a T2‐weighted intramedullary hyperintensity compared to those without (.7, .45−.85 vs. .54, .4−.8; P = .01).

Conclusion and Clinical Relevance

Mean diffusivity and AD could serve as surrogates of severity of spinal cord injury and are complementary to the clinical exam in dogs with thoracolumbar IVDH.

Magnetization Transfer Ratio and Morphometrics of the Spinal Cord Associates with Surgical Recovery in Patients with Degenerative Cervical Myelopathy

OBJECTIVES

We assessed the prognostic value of the preoperative magnetization transfer ratio (MTR) and morphometrics of the spinal cord in patients with degenerative cervical myelopathy (DCM) in a longitudinal cohort study.

METHODS

Thirteen subjects with DCM underwent 3T magnetization transfer imaging. The MTR was calculated for the spinal cord regions and specific white matter tracts. Morphometric measures were extracted. Clinical (modified Japanese Orthopaedics Association [mJOA] and Nurick scale scores) and health-related quality of life scores were assessed before and after cervical decompression surgery. The association between the magnetic resonance imaging (MRI) metrics and postoperative recovery was assessed (Spearman's correlation). Receiver operating characteristics were used to assess the accuracy of MRI metrics in identifying ≥50% recovery in function.

RESULTS

Preoperative anterior cord MTRs were associated with recovery in mJOA scores (ρ = 0.608; P = 0.036; area under the curve [AUC], 0.66). Preoperative lateral cord MTR correlated with the neck disability index (ρ = 0.699; P = 0.011) and pain interference scale (ρ = 0.732; P = 0.007). Preoperative rubrospinal tract MTR was associated with mJOA score recovery (ρ = 0.573; P = 0.041; AUC, 0.86). Preoperative corticospinal tract and reticulospinal MTRs were related to recovery in pain interference scores (ρ = 0.591; P = 0.033; and ρ = 0.583; P = 0.035, respectively). Eccentricity of the cord was associated with Nurick scores (ρ = 0.606; P = 0.028) and mJOA scores (ρ = 0.651; P = 0.025; AUC, 0.92).

CONCLUSIONS

Preoperative MTR and eccentricity measurements of the spinal cord have prognostic value in assessing the response to surgery and recovery in patients with DCM. Advanced MRI and atlas-based postprocessing techniques can inform interventions and advance the healthcare received by patients with DCM.

Simultaneous voxel-wise analysis of brain and spinal cord morphometry and microstructure within the SPM framework

To validate a simultaneous analysis tool for the brain and cervical cord embedded in the statistical parametric mapping (SPM) framework, we compared trauma-induced macro- and microstructural changes in spinal cord injury (SCI) patients to controls. The findings were compared with results obtained from existing processing tools that assess the brain and spinal cord separately. A probabilistic brain-spinal cord template (BSC) was generated using a generative semi-supervised modelling approach. The template was incorporated into the pre-processing pipeline of voxel-based morphometry and voxel-based quantification analyses in SPM. This approach was validated on T1-weighted scans and multiparameter maps, by assessing trauma-induced changes in SCI patients relative to controls and comparing the findings with the outcome from existing analytical tools. Consistency of the MRI measures was assessed using intraclass correlation coefficients (ICC). The SPM approach using the BSC template revealed trauma-induced changes across the sensorimotor system in the cord and brain in SCI patients. These changes were confirmed with established approaches covering brain or cord, separately. The ICC in the brain was high within regions of interest, such as the sensorimotor cortices, corticospinal tracts and thalamus. The simultaneous voxel-wise analysis of brain and cervical spinal cord was performed in a unique SPM-based framework incorporating pre-processing and statistical analysis in the same environment. Validation based on a SCI cohort demonstrated that the new processing approach based on the brain and cord is comparable to available processing tools, while offering the advantage of performing the analysis simultaneously across the neuraxis.

Spinal cord atrophy in a primary progressive multiple sclerosis trial: Improved sample size using GBSI

•

The GBSI provided clinically meaningful measurements of spinal cord atrophy, with low sample size.

•

Deriving spinal cord atrophy from brain MRI using the GBSI is easier than spinal cord MRI.

•

Spinal cord atrophy on GBSI could be used as a secondary outcome measure.

Background

We aimed to evaluate the implications for clinical trial design of the generalised boundary-shift integral (GBSI) for spinal cord atrophy measurement.

Methods

We included 220 primary-progressive multiple sclerosis patients from a phase 2 clinical trial, with baseline and week-48 3DT1-weighted MRI of the brain and spinal cord (1 × 1 × 1 mm³), acquired separately. We obtained segmentation-based cross-sectional spinal cord area (CSA) at C1-2 (from both brain and spinal cord MRI) and C2-5 levels (from spinal cord MRI) using DeepSeg, and, then, we computed corresponding GBSI.

Results

Depending on the spinal cord segment, we included 67.4–98.1% patients for CSA measurements, and 66.9–84.2% for GBSI. Spinal cord atrophy measurements obtained with GBSI had lower measurement variability, than corresponding CSA. Looking at the image noise floor, the lowest median standard deviation of the MRI signal within the cerebrospinal fluid surrounding the spinal cord was found on brain MRI at the C1-2 level. Spinal cord atrophy derived from brain MRI was related to the corresponding measures from dedicated spinal cord MRI, more strongly for GBSI than CSA. Spinal cord atrophy measurements using GBSI, but not CSA, were associated with upper and lower limb motor progression.

Discussion

Notwithstanding the reduced measurement variability, the clinical correlates, and the possibility of using brain acquisitions, spinal cord atrophy using GBSI should remain a secondary outcome measure in MS studies, until further advancements increase the quality of acquisition and reliability of processing.

Dynamic Functional Connectivity of Resting-State Spinal Cord fMRI Reveals Fine-Grained Intrinsic Architecture

The neuroimaging community has shown tremendous interest in exploring the brain's spontaneous activity using functional magnetic resonance imaging (fMRI). On the contrary, the spinal cord has been largely overlooked despite its pivotal role in processing sensorimotor signals. Only a handful of studies have probed the organization of spinal resting-state fluctuations, always using static measures of connectivity. Many innovative approaches have emerged for analyzing dynamics of brain fMRI, but they have not yet been applied to the spinal cord, although they could help disentangle its functional architecture. Here, we leverage a dynamic connectivity method based on the clustering of hemodynamic-informed transients to unravel the rich dynamic organization of spinal resting-state signals. We test this approach in 19 healthy subjects, uncovering fine-grained spinal components and highlighting their neuroanatomical and physiological nature. We provide a versatile tool, the spinal innovation-driven co-activation patterns (SpiCiCAP) framework, to characterize spinal circuits during rest and task, as well as their disruption in neurological disorders.

Osprey: Open-source processing, reconstruction & estimation of magnetic resonance spectroscopy data

BACKGROUND

Processing and quantitative analysis of magnetic resonance spectroscopy (MRS) data are far from standardized and require interfacing with third-party software. Here, we present Osprey, a fully integrated open-source data analysis pipeline for MRS data, with seamless integration of pre-processing, linear-combination modelling, quantification, and data visualization.

NEW METHOD

Osprey loads multiple common MRS data formats, performs phased-array coil combination, frequency-and phase-correction of individual transients, signal averaging and Fourier transformation. Linear combination modelling of the processed spectrum is carried out using simulated basis sets and a spline baseline. The MRS voxel is coregistered to an anatomical image, which is segmented for tissue correction and quantification is performed based upon modelling parameters and tissue segmentation. The results of each analysis step are visualized in the Osprey GUI. The analysis pipeline is demonstrated in 12 PRESS, 11 MEGA-PRESS, and 8 HERMES datasets acquired in healthy subjects.

RESULTS

Osprey successfully loads, processes, models, and quantifies MRS data acquired with a variety of conventional and spectral editing techniques.

COMPARISON WITH EXISTING METHOD(S)

Osprey is the first MRS software to combine uniform pre-processing, linear-combination modelling, tissue correction and quantification into a coherent ecosystem. Compared to existing compiled, often closed-source modelling software, Osprey's open-source code philosophy allows researchers to integrate state-of-the-art data processing and modelling routines, and potentially converge towards standardization of analysis.

CONCLUSIONS

Osprey combines robust, peer-reviewed data processing methods into a modular workflow that is easily augmented by community developers, allowing the rapid implementation of new methods.

The Evaluation and Prediction of Laminoplasty Surgery Outcome in Patients with Degenerative Cervical Myelopathy Using Diffusion Tensor MRI

BACKGROUND AND PURPOSE

DTI has been proved valuable for the diagnosis of degenerative cervical myelopathy, whereas its capacity for predicting the outcome of surgery is still under debate. Here we conduct a prospective cohort study to analyze the capacity of DTI for evaluating and predicting laminoplasty surgery outcome for degenerative cervical myelopathy.

MATERIALS AND METHODS

We recruited 55 patients with degenerative cervical myelopathy who underwent DTI before surgery and at 3- and 6-month follow-up stages, and 20 healthy subjects. For clinical assessment, the modified Japanese Orthopedic Association scale was recorded for each patient at different stages. DTI metrics were compared between patients before surgery and healthy subjects. Spearman correlation and receiver operating characteristic were used to analyze the evaluation and prediction capacity of DTI for the modified Japanese Orthopedic Association scale, respectively. We analyzed different vertebral levels: maximal compression level, average of all compression levels, and C2 level.

RESULTS

DTI metrics were significantly different between patients before surgery and healthy subjects. Before surgery, DTI for the maximal compression level or DTI for the average of all compression levels had no significant correlation with the modified Japanese Orthopedic Association scale. For all stages, DTI at the C2 level was correlated with the modified Japanese Orthopedic Association scale. DTI metrics at the C2 level before surgery were significantly correlated with the postoperative modified Japanese Orthopedic Association scale recovery rate. Receiver operating characteristic analysis demonstrated that fractional anisotropy at C2 was capable of predicting the postoperative modified Japanese Orthopedic Association scale recovery rate (P = .04).

CONCLUSIONS

The DTI metrics before laminoplasty surgery, especially fractional anisotropy at the C2 level, have the potential for evaluating and predicting the degenerative cervical myelopathy surgery outcome.

Multi-parametric quantitative in vivo spinal cord MRI with unified signal readout and image denoising

Multi-parametric quantitative MRI (qMRI) of the spinal cord is a promising non-invasive tool to probe early microstructural damage in neurological disorders. It is usually performed in vivo by combining acquisitions with multiple signal readouts, which exhibit different thermal noise levels, geometrical distortions and susceptibility to physiological noise. This ultimately hinders joint multi-contrast modelling and makes the geometric correspondence of parametric maps challenging. We propose an approach to overcome these limitations, by implementing state-of-the-art microstructural MRI of the spinal cord with a unified signal readout in vivo (i.e. with matched spatial encoding parameters across a range of imaging contrasts). We base our acquisition on single-shot echo planar imaging with reduced field-of-view, and obtain data from two different vendors (vendor 1: Philips Achieva; vendor 2: Siemens Prisma). Importantly, the unified acquisition allows us to compare signal and noise across contrasts, thus enabling overall quality enhancement via multi-contrast image denoising methods. As a proof-of-concept, here we provide a demonstration with one such method, known as Marchenko-Pastur (MP) Principal Component Analysis (PCA) denoising. MP-PCA is a singular value (SV) decomposition truncation approach that relies on redundant acquisitions, i.e. such that the number of measurements is large compared to the number of components that are maintained in the truncated SV decomposition. Here we used in vivo and synthetic data to test whether a unified readout enables more efficient MP-PCA denoising of less redundant acquisitions, since these can be denoised jointly with more redundant ones. We demonstrate that a unified readout provides robust multi-parametric maps, including diffusion and kurtosis tensors from diffusion MRI, myelin metrics from two-pool magnetisation transfer, and T1 and T2 from relaxometry. Moreover, we show that MP-PCA improves the quality of our multi-contrast acquisitions, since it reduces the coefficient of variation (i.e. variability) by up to 17% for mean kurtosis, 8% for bound pool fraction (myelin-sensitive), and 13% for T1, while enabling more efficient denoising of modalities limited in redundancy (e.g. relaxometry). In conclusion, multi-parametric spinal cord qMRI with unified readout is feasible and provides robust microstructural metrics with matched resolution and distortions, whose quality benefits from multi-contrast denoising methods such as MP-PCA.

Cervical Canal Morphology: Effects of Neck Flexion in Normal Condition: New Elements for Biomechanical Simulations and Surgical Management

STUDY DESIGN

Continuous measurements and computation of absolute metrics of cervical subarachnoid space (CSS) and spinal cord (SC) geometries proposed are based on in vivo magnetic resonance imaging and 3D reconstruction.

OBJECTIVE

The aim of the study is to offer a new methodology to continuously characterize and to quantify the detailed morphology of the CSS and the cervical SC in 3D for healthy subjects in both neutral supine and flexion.

SUMMARY OF BACKGROUND DATA

To the best of our knowledge, no study provides a morphological quantification by absolute indices based on the 3D reconstruction of SC and CSS thanks to in vivo magnetic resonance imaging. Moreover, no study provides a continuous description of the geometries.

METHODS

Absolute indices of SC (cross-sectional area, compression ratio, position in the canal, length) and of CSS (cross-sectional area, occupational ratio, lengths) were computed by measures from 3D semi-automatic reconstructions of high resolution in vivo magnetic resonance images (3D T2-SPACE sequence) on healthy subjects (N = 11) for two postures: supine neutral and flexion neck positions. The variability induced by the semi-automatic reconstruction and by the landmarks positioning were investigated by preliminary sensitivity analyses. Inter and intra-variability were also quantified on a randomly chosen part of our population (N = 5).

RESULTS

The length and cross-sectional area of SC are significantly different (P < 0.05) in flexion compared with neutral neck position. Spinal cord stays centered in the canal for both postures. However, the cross-sectional area of CSS is submitted to low variation after C3 vertebra for both postures. Occupational ratio (OR) and compression ratio (CR) after C3 are significantly lower in flexion.

CONCLUSION

This study presented interpretations of morphological measures: (1) left-right stability (described by the Left-Right eccentricity index) ensured by the denticulate ligaments and the nerve roots attached to the dural sheaths, (2) a Poisson effect of the SC was partially notified through its axial (antero-posterior [AP] diameter, OR, CR) and its longitudinal geometrical descriptions (length of spinal cord [LSC]). Such morphological data can be useful for geometrical finite element modeling and could now be used to compare with injured or symptomatic subjects.

LEVEL OF EVIDENCE

3.

Ultrahigh-resolution quantitative spinal cord MRI at 9.4T

PURPOSE

To present the results of the first human spinal cord in vivo MRI scans at 9.4T.

METHODS

A human brain coil was used to image the human spinal cord at 9.4T. All anatomical images were acquired with a T₂ *-weighted gradient-echo sequence. A comparison of the influence of four different B₀ shimming routines on the image quality was performed. Intrinsic signal-to-noise-ratio maps were determined using a pseudo-multiple replica approach. Measurements with different echo times were compared and processed to one multiecho data image combination image. Based on the multiecho acquisitions, T₂ *-relaxation time maps were calculated. Algorithmic spinal cord detection and gray matter/white matter segmentation were tested.

RESULTS

An echo time between 9 and 13.8 ms compromised best between gray matter/white matter contrast and image quality. A maximum in-plane resolution of 0.15 × 0.15 mm² was achieved for anatomical images. These images offered excellent image quality and made small structures of the spinal cord visible. The scanner vendor implemented B₀ shimming routine performed best during this work. Intrinsic signal-to-noise-ratio values of between 6600 and 8060 at the upper cervical spinal cord were achieved. Detection and segmentation worked reliably. An average T₂ *-time of 24.88 ms ± 6.68 ms for gray matter and 19.37 ms ± 8.66 ms for white matter was calculated.

CONCLUSION

The proposed human brain coil can be used to image the spinal cord. The maximum in-plane resolution in this work was higher compared with the 7T results from the literature. The 9.4T acquisitions made the small structures of the spinal cord clearly visible.

Spinal Cord Atrophy in Neuromyelitis Optica Spectrum Disorders Is Spatially Related to Cord Lesions and Disability

Background The spinal cord is commonly involved in patients with neuromyelitis optica spectrum disorders (NMOSDs). However, the relationship between inflammation and atrophy remains unclear. Purpose To characterize the spatial distribution of T1-hypointense lesions in the spinal cord at MRI, its association with cord atrophy, and its correlation with disability in participants with NMOSDs. Materials and Methods This prospective study evaluated three-dimensional T1-weighted spinal cord MRI scans in seropositive participants with NMOSDs and in age-matched healthy control participants acquired between February 2010 and July 2018. Binary masks of T1-hypointense lesions and lesion probability maps were produced. Cross-sectional area of the cervical and upper thoracic cord (down to T3 level) was calculated with the active-surface method. Full factorial models were used to assess cord atrophy in participants with NMOSDs. Correlations between cord atrophy and clinical and brain MRI measures were investigated with multiple regression models. Results A total of 52 participants with NMOSDs (mean age ± standard deviation, 44 years ± 15; 45 women) and 28 age-matched healthy control participants (mean age, 44 years ± 13; 16 women) were evaluated. Thirty-eight of 52 (73%) participants with NMOSDs had T1-hypointense cord lesions. No cord lesions were detected in the healthy control participants. Lesion probability maps showed a predominant involvement of the upper cervical (C2-C4) and upper thoracic (T1-T3 level) cord. The greater involvement of C1-C4 survived Bonferroni correction (P value range, .007-.04), with a higher percentage lesion extent in the gray matter (P < .001). Atrophy colocalized with focal cord lesions and correlated with pyramidal subscore (r ranging from -0.53 to -0.40; P < .001) and sensitive subscore (r ranging from -0.48 to -0.46; P = .001) of the Expanded Disability Status Scale. Participants without cord lesions had no cord atrophy. Conclusion In participants with neuromyelitis optica spectrum disorders, focal areas of spinal cord atrophy at MRI were topographically associated with lesions and correlated to motor and sensory disability. Participants without visible cord lesions had no atrophy. © RSNA, 2020 Online supplemental material is available for this article.

Voxelwise analysis of diffusion MRI of cervical spinal cord using tract-based spatial statistics

Robust voxelwise analysis using tract-based spatial statistics (TBSS) together with permutation statistical method is standardly used in analyzing diffusion tensor imaging (DTI) of brain. A similar analytical method could be useful when studying DTI of cervical spinal cord. Based on anatomical data of sixty-four healthy volunteers, white (WM) and gray matter (GM) masks were created and subsequently registered into DTI space. Using TBSS, two skeleton types were created (single line and dilated for WM as well as GM). From anatomical data, percentage rates of overlap were calculated for all skeletons in relation to WM and GM masks. Voxelwise analysis of fractional anisotropy values depending on age and sex was conducted. Correlation of fraction anisotropy values with age of subjects was also evaluated. The two WM skeleton types showed a high overlap rate with WM masks (~94%); GM skeletons showed lower rates (56% and 42%, respectively, for single line and dilated). WM and GM areas where fraction anisotropy values differ between sexes were identified (p < .05). Furthermore, using voxelwise analysis such WM voxels were identified where fraction anisotropy values differ depending on age (p < .05) and in these voxels linear dependence of fraction anisotropy and age (r = -0.57, p < .001) was confirmed by regression analysis. This dependence was not proven when using WM anatomical masks (r = -0.21, p = .10). The analytical approach presented shown to be useful for group analysis of DTI data for cervical spinal cord.

Tract-specific analysis improves sensitivity of spinal cord diffusion MRI to cross-sectional and longitudinal changes in amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is a late-onset fatal neurodegenerative disease that causes progressive degeneration of motor neurons in the brain and the spinal cord. Corticospinal tract degeneration is a defining feature of ALS. However, there have been very few longitudinal, controlled studies assessing diffusion MRI (dMRI) metrics in different fiber tracts along the spinal cord in general or the corticospinal tract in particular. Here we demonstrate that a tract-specific analysis, with segmentation of ascending and descending tracts in the spinal cord white matter, substantially increases the sensitivity of dMRI to disease-related changes in ALS. Our work also identifies the tracts and spinal levels affected in ALS, supporting electrophysiologic and pathologic evidence of involvement of sensory pathways in ALS. We note changes in diffusion metrics and cord cross-sectional area, with enhanced sensitivity to disease effects through a multimodal analysis, and with strong correlations between these metrics and spinal components of ALSFRS-R.

Resting-state brain and spinal cord networks in humans are functionally integrated

In the absence of any task, both the brain and spinal cord exhibit spontaneous intrinsic activity organised in a set of functionally relevant neural networks. However, whether such resting-state networks (RSNs) are interconnected across the brain and spinal cord is unclear. Here, we used a unique scanning protocol to acquire functional images of both brain and cervical spinal cord (CSC) simultaneously and examined their spatiotemporal correspondence in humans. We show that the brain and spinal cord activities are strongly correlated during rest periods, and specific spinal cord regions are functionally linked to consistently reported brain sensorimotor RSNs. The functional organisation of these networks follows well-established anatomical principles, including the contralateral correspondence between the spinal hemicords and brain hemispheres as well as sensory versus motor segregation of neural pathways along the brain–spinal cord axis. Thus, our findings reveal a unified functional organisation of sensorimotor networks in the entire central nervous system (CNS) at rest.

This neuroimaging study reveals novel insights into the functional organization of resting-state networks in the brain and spinal cord, such as the contralateral correspondence between the two halves of the brain and spinal cord, and segregation of sensory versus motor neural pathways along this axis.

Spinal cord atrophy as a measure of severity of myelopathy in adrenoleukodystrophy

All men and most women with X-linked adrenoleukodystrophy (ALD) develop myelopathy in adulthood. As clinical trials with new potential disease-modifying therapies are emerging, sensitive outcome measures for quantifying myelopathy are needed. This prospective cohort study evaluated spinal cord size (cross-sectional area - CSA) and shape (eccentricity) as potential new quantitative outcome measures for myelopathy in ALD. Seventy-four baseline magnetic resonance imaging (MRI) scans, acquired in 42 male ALD patients and 32 age-matched healthy controls, and 26 follow-up scans of ALD patients were included in the study. We used routine T1 -weighted MRI sequences to measure mean CSA, eccentricity, right-left and anteroposterior diameters in the cervical spinal cord. We compared MRI measurements between groups and correlated CSA with clinical outcome measures of disease severity. Longitudinally, we compared MRI measurements between baseline and 1-year follow-up. CSA was significantly smaller in patients compared to controls on all measured spinal cord levels (P < .001). The difference was completely explained by the effect of the symptomatic subgroup. Furthermore, the spinal cord showed flattening (higher eccentricity and smaller anteroposterior diameters) in patients. CSA correlated strongly with all clinical measures of severity of myelopathy. There was no detectable change in CSA after 1-year follow-up. The cervical spinal cord in symptomatic ALD patients is smaller and flattened compared to controls, possibly due to atrophy of the dorsal columns. CSA is a reliable marker of disease severity and can be a valuable outcome measure in long-term follow-up studies in ALD. SYNOPSIS: A prospective cohort study in 42 adrenoleukodystrophy (ALD) patients and 32 controls demonstrated that the spinal cord cross-sectional area of patients is smaller compared to healthy controls and correlates with severity of myelopathy in patients, hence it could be valuable as a much needed surrogate outcome measure.

Prospective Follow-up of Intramedullary Slitlike Cavities: A Consecutive Series of 48 Patients

Object: Predicting whether intramedullary slitlike cavity (SC) will worsen over time or remain stable is an outstanding clinical challenge. The aim of this study was to identify early features of SC (clinical and magnetic resonance imaging [MRI] findings).

Methods: We prospectively included all patients referred to our institution following the discovery of a SC and divided them in two groups: typical SC (defined as a cavity spanning fewer than three vertebrae, not enlarging the spinal cord, and located at the midline between the anterior third and posterior two-thirds of the spinal cord) or atypical SC (all others). Clinical evolution and changes in MRI features were evaluated during follow-up. In some patients, diffusion tensor imaging was performed and cervical cord cross-sectional area was analyzed.

Results: A total of 48 consecutive patients were included in the study. The mean follow-up was 58 months. Of the seven patients presenting with deficits at first consultation, two worsened and five remained stable. Of the 41 patients without deficits, seven worsened and 34 remained stable. None of the patients developed severe motor deficits or experienced enlargement of the cavity; 7% of patients who presented with typical SC worsened compared with 35% with atypical SC. The negative predictive value was 0.93 (P = 0.02).

Conclusion: Most patients remained stable and a subset of patients developed minor motor deficits. For clinical management, we propose surveillance of patients with a typical SC and close follow-up of those with an atypical SC and/or presenting with deficits.

Assessing the spatial distribution of cervical spinal cord activity during tactile stimulation of the upper extremity in humans with functional magnetic resonance imaging

Dermatomal maps are a mainstay of clinical practice and provide information on the spatial distribution of the cutaneous innervation of spinal nerves. Dermatomal deficits can help isolate the level of spinal nerve root involvement in spinal conditions and guide clinicians in diagnosis and treatment. Dermatomal maps, however, have limitations, and the spatial distribution of spinal cord sensory activity in humans remains to be quantitatively assessed. Here we used spinal cord functional MRI to map and quantitatively compare the spatial distribution of sensory spinal cord activity during tactile stimulation of the left and right lateral shoulders (i.e. C5 dermatome) and dorsal third digits of the hands (i.e., C7 dermatome) in healthy humans (n = 24, age = 36.8 ± 11.8 years). Based on the central sites for processing of innocuous tactile sensory information, we hypothesized that the activity would be localized more to the ipsilateral dorsal spinal cord with the lateral shoulder stimulation activity being localized more superiorly than the dorsal third digit. The findings demonstrate lateralization of the activity with the left- and right-sided stimuli having more activation in the ipsilateral hemicord. Contradictory to our hypotheses, the activity for both stimulation sites was spread across the dorsal and ventral hemicords and did not demonstrate a clear superior-inferior localization. Instead, the activity for both stimuli had a broader than expected distribution, extending across the C5, C6, and C7 spinal cord segments. We highlight the complexity of the human spinal cord neuroanatomy and several sources of variability that may explain the observed patterns of activity. While the findings were not completely consistent with our a priori hypotheses, this study provides a foundation for continued work and is an important step towards developing normative quantitative spinal cord measures of sensory function, which may become useful objective MRI-based biomarkers of neurological injury and improve the management of spinal disorders.

Combined brain and spinal FDG PET allows differentiation between ALS and ALS mimics

PURPOSE

Amyotrophic lateral sclerosis (ALS) is a progressive neurodegenerative disorder with on average a 1-year delay between symptom onset and diagnosis. Studies have demonstrated the value of [¹⁸F]-FDG PET as a sensitive diagnostic biomarker, but the discriminatory potential to differentiate ALS from patients with symptoms mimicking ALS has not been investigated. We investigated the combination of brain and spine [¹⁸F]-FDG PET-CT for differential diagnosis between ALS and ALS mimics in a real-life clinical diagnostic setting.

METHODS

Patients with a suspected diagnosis of ALS (n = 98; 64.8 ± 11 years; 61 M) underwent brain and spine [¹⁸F]-FDG PET-CT scans. In 62 patients, ALS diagnosis was confirmed (67.8 ± 10 years; 35 M) after longitudinal follow-up (average 18.1 ± 8.4 months). In 23 patients, another disease was diagnosed (ALS mimics, 60.9 ± 12.9 years; 17 M) and 13 had a variant motor neuron disease, primary lateral sclerosis (PLS; n = 4; 53.6 ± 2.5 years; 2 M) and progressive muscular atrophy (PMA; n = 9; 58.4 ± 7.3 years; 7 M). Spine metabolism was determined after manual and automated segmentation. VOI- and voxel-based comparisons were performed. Moreover, a support vector machine (SVM) approach was applied to investigate the discriminative power of regional brain metabolism, spine metabolism and the combination of both.

RESULTS

Brain metabolism was very similar between ALS mimics and ALS, whereas cervical and thoracic spine metabolism was significantly different (in standardised uptake values; cervical: ALS 2.1 ± 0.5, ALS mimics 1.9 ± 0.4; thoracic: ALS 1.8 ± 0.3, ALS mimics 1.5 ± 0.3). As both brain and spine metabolisms were very similar between ALS mimics and PLS/PMA, groups were pooled for accuracy analyses. Mean discrimination accuracy was 65.4%, 80.0% and 81.5%, using only brain metabolism, using spine metabolism and using both, respectively.

CONCLUSION

The combination of brain and spine FDG PET-CT with SVM classification is useful as discriminative biomarker between ALS and ALS mimics in a real-life clinical setting.

Heritability of cervical spinal cord structure

OBJECTIVE

Measures of spinal cord structure can be a useful phenotype to track disease severity and development; this observational study measures the hereditability of cervical spinal cord anatomy and its correlates in healthy human beings.

METHODS

Twin data from the Human Connectome Project were analyzed with semiautomated spinal cord segmentation, evaluating test-retest reliability and broad-sense heritability with an AE model. Relationships between spinal cord metrics, general physical measures, regional brain structural measures, and motor function were assessed.

RESULTS

We found that the spinal cord C2 cross-sectional area (CSA), left-right width (LRW), and anterior-posterior width (APW) are highly heritable (85%-91%). All measures were highly correlated with the brain volume, and CSA only was positively correlated with thalamic volumes (p = 0.005) but negatively correlated with the occipital cortex area (p = 0.001). LRW was correlated with the participant's height (p = 0.00027). The subjects' sex significantly influenced these metrics. Analyses of a test-retest data set confirmed validity of the approach.

CONCLUSIONS

This study provides the evidence of genetic influence on spinal cord structure. MRI metrics of cervical spinal cord anatomy are robust and not easily influenced by nonpathological environmental factors, providing a useful metric for monitoring normal development and progression of neurodegenerative disorders affecting the spinal cord, including-but not limited to-spinal cord injury and MS.

MAGNIMS consensus recommendations on the use of brain and spinal cord atrophy measures in clinical practice

Early evaluation of treatment response and prediction of disease evolution are key issues in the management of people with multiple sclerosis (MS). In the past 20 years, MRI has become the most useful paraclinical tool in both situations and is used clinically to assess the inflammatory component of the disease, particularly the presence and evolution of focal lesions - the pathological hallmark of MS. However, diffuse neurodegenerative processes that are at least partly independent of inflammatory mechanisms can develop early in people with MS and are closely related to disability. The effects of these neurodegenerative processes at a macroscopic level can be quantified by estimation of brain and spinal cord atrophy with MRI. MRI measurements of atrophy in MS have also been proposed as a complementary approach to lesion assessment to facilitate the prediction of clinical outcomes and to assess treatment responses. In this Consensus statement, the Magnetic Resonance Imaging in MS (MAGNIMS) study group critically review the application of brain and spinal cord atrophy in clinical practice in the management of MS, considering the role of atrophy measures in prognosis and treatment monitoring and the barriers to clinical use of these measures. On the basis of this review, the group makes consensus statements and recommendations for future research.

Tracking the neurodegenerative gradient after spinal cord injury

OBJECTIVE

To quantify neurodegenerative changes along the cervical spinal cord rostral to a spinal cord injury (SCI) by means of quantitative MRI (qMRI) and to determine its relationship with clinical impairment.

METHODS

Thirty chronic SCI patients (15 tetraplegics and 15 paraplegics) and 23 healthy controls underwent a high-resolution T1-weighted and myelin-sensitive magnetization transfer (MT) MRI. We assessed macro- and microstructural changes along the cervical cord from levels C1 to C4, calculating cross-sectional spinal cord area, its anterior-posterior and left-right width and myelin content (i.e. MT). Regression analysis determined associations between qMRI parameters and clinical impairment.

RESULTS

In SCI patients, cord area decreased by 2.67 mm² (p = 0.004) and left-right width decreased by 0.35 mm (p = 0.002) per cervical cord level in the caudal direction when compared to the healthy controls. This gradient of neurodegeneration was greater in tetraplegic than paraplegics in the cross-sectional cervical cord area (by 3.28 mm², p = 0.011), left-right width (by 0.36 mm, p = 0.03), and mean cord MT (by 0.13%, p = 0.04), but independant of lesion severity (p > 0.05). Higher lesion level was associated with greater magnitudes of neurodegeneration. Greater loss in myelin content in the dorsal columns and spinothalamic tract was associated with worse light touch (p = 0.016) and pin prick score (p = 0.024), respectively.

CONCLUSIONS

A gradient of neurodegeneration is evident in the cervical cord remote from a SCI. Tract-specific associations with appropriate clinical outcomes highlight that remote neurodegenerative changes are clinically eloquent. Monitoring the neurodegenerative gradient could be used to track treatment effects of regenerative and neuroprotective agents, both in trials targeting cervical and thoracic SCI patients.

Multiscale Imaging Approach for Studying the Central Nervous System: Methodology and Perspective

Non-invasive imaging methods have become essential tools for understanding the central nervous system (CNS) in health and disease. In particular, magnetic resonance imaging (MRI) techniques provide information about the anatomy, microstructure, and function of the brain and spinal cord in vivo non-invasively. However, MRI is limited by its spatial resolution and signal specificity. In order to mitigate these shortcomings, it is crucial to validate MRI with an array of ancillary ex vivo imaging techniques. These techniques include histological methods, such as light and electron microscopy (EM), which can provide specific information on the tissue structure in healthy and diseased brain and spinal cord, at cellular and subcellular level. However, these conventional histological techniques are intrinsically two-dimensional (2D) and, as a result of sectioning, lack volumetric information of the tissue. This limitation can be overcome with genuine three-dimensional (3D) imaging approaches of the tissue. 3D highly resolved information of the CNS achievable by means of other imaging techniques can complement and improve the interpretation of MRI measurements. In this article, we provide an overview of different 3D imaging techniques that can be used to validate MRI. As an example, we introduce an approach of how to combine diffusion MRI and synchrotron X-ray phase contrast tomography (SXRPCT) data. Our approach paves the way for a new multiscale assessment of the CNS allowing to validate and to improve our understanding of in vivo imaging (such as MRI).

Evolution of Magnetic Resonance Imaging as Predictors and Correlates of Functional Outcome after Spinal Cord Contusion Injury in the Rat

Clinical methods for determining the severity of traumatic spinal cord injury (SCI) and long-term functional outcome in the acute setting are limited in their prognostic accuracy because of the heterogeneity of injury and dynamic injury progression. The aim of this study was to evaluate the time course and sensitivity of advanced magnetic resonance imaging (MRI) methods to neurological function after SCI in a rat contusion model. Rats received a graded contusion injury at T10 using a weight-drop apparatus. MRI consisted of morphological measures from T₂-weighted imaging, quantitative T₂ imaging, and diffusion-weighted imaging (DWI) at 1, 30, and 90 days post-injury (dpi). The derived metrics were compared with neurological function assessed using weekly Basso, Beattie, and Bresnahan (BBB) locomotor scoring and return of reflexive micturition function. At the acute time point (1 dpi), diffusion metrics sensitive to axonal injury at the injury epicenter had the strongest correlation with time-matched BBB scores and best predicted 90-dpi BBB scores. At 30 dpi, axonal water fraction derived from DWI and T₂ values were both correlated with time-matched locomotor scores. At the chronic time point (90 dpi), cross-sectional area was most closely correlated to BBB. Overall, the results demonstrate differential sensitivity of MRI metrics at different time points after injury, but the metrics follow the expected pathology of acute axonal injury followed by continued degeneration and finally a terminal level of atrophy. Specificity of DWI in the acute setting may make it impactful as a prognostic tool while T₂ imaging provided the most information about injury severity in chronic injury.

Unsuspected Involvement of Spinal Cord in Alzheimer Disease

Objective: Brain atrophy is an established biomarker for dementia, yet spinal cord involvement has not been investigated to date. As the spinal cord is relaying sensorimotor control signals from the cortex to the peripheral nervous system and vice-versa, it is indeed a very interesting question to assess whether it is affected by atrophy due to a disease that is known for its involvement of cognitive domains first and foremost, with motor symptoms being clinically assessed too. We, therefore, hypothesize that in Alzheimer's disease (AD), severe atrophy can affect the spinal cord too and that spinal cord atrophy is indeed an important in vivo imaging biomarker contributing to understanding neurodegeneration associated with dementia. Methods: 3DT1 images of 31 AD and 35 healthy control (HC) subjects were processed to calculate volume of brain structures and cross-sectional area (CSA) and volume (CSV) of the cervical cord [per vertebra as well as the C2-C3 pair (CSA23 and CSV23)]. Correlated features (ρ > 0.7) were removed, and the best subset identified for patients' classification with the Random Forest algorithm. General linear model regression was used to find significant differences between groups (p ≤ 0.05). Linear regression was implemented to assess the explained variance of the Mini-Mental State Examination (MMSE) score as a dependent variable with the best features as predictors. Results: Spinal cord features were significantly reduced in AD, independently of brain volumes. Patients classification reached 76% accuracy when including CSA23 together with volumes of hippocampi, left amygdala, white and gray matter, with 74% sensitivity and 78% specificity. CSA23 alone explained 13% of MMSE variance. Discussion: Our findings reveal that C2-C3 spinal cord atrophy contributes to discriminate AD from HC, together with more established features. The results show that CSA23, calculated from the same 3DT1 scan as all other brain volumes (including right and left hippocampi), has a considerable weight in classification tasks warranting further investigations. Together with recent studies revealing that AD atrophy is spread beyond the temporal lobes, our result adds the spinal cord to a number of unsuspected regions involved in the disease. Interestingly, spinal cord atrophy explains also cognitive scores, which could significantly impact how we model sensorimotor control in degenerative diseases with a primary cognitive domain involvement. Prospective studies should be purposely designed to understand the mechanisms of atrophy and the role of the spinal cord in AD.

Myelin water imaging data analysis in less than one minute

PURPOSE

Based on a deep learning neural network (NN) algorithm, a super fast and easy to implement data analysis method was proposed for myelin water imaging (MWI) to calculate the myelin water fraction (MWF).

METHODS

A NN was constructed and trained on MWI data acquired by a 32-echo 3D gradient and spin echo (GRASE) sequence. Ground truth labels were created by regularized non-negative least squares (NNLS) with stimulated echo corrections. Voxel-wise GRASE data from 5 brains (4 healthy, 1 multiple sclerosis (MS)) were used for NN training. The trained NN was tested on 2 healthy brains, 1 MS brain with segmented lesions, 1 healthy spinal cord, and 1 healthy brain acquired from a different scanner.

RESULTS

Production of whole brain MWF maps in approximately 33 ​s can be achieved by a trained NN without graphics card acceleration. For all testing regions, no visual differences between NN and NNLS MWF maps were observed, and no obvious regional biases were found. Quantitatively, all voxels exhibited excellent agreement between NN and NNLS (all R²>0.98, p ​< ​0.001, mean absolute error <0.01).

CONCLUSION

The time for accurate MWF calculation can be dramatically reduced to less than 1 ​min by the proposed NN, addressing one of the barriers facing future clinical feasibility of MWI.

Spatial distribution of multiple sclerosis lesions in the cervical spinal cord

Spinal cord lesions detected on MRI hold important diagnostic and prognostic value for multiple sclerosis. Previous attempts to correlate lesion burden with clinical status have had limited success, however, suggesting that lesion location may be a contributor. Our aim was to explore the spatial distribution of multiple sclerosis lesions in the cervical spinal cord, with respect to clinical status. We included 642 suspected or confirmed multiple sclerosis patients (31 clinically isolated syndrome, and 416 relapsing-remitting, 84 secondary progressive, and 73 primary progressive multiple sclerosis) from 13 clinical sites. Cervical spine lesions were manually delineated on T2- and T2*-weighted axial and sagittal MRI scans acquired at 3 or 7 T. With an automatic publicly-available analysis pipeline we produced voxelwise lesion frequency maps to identify predilection sites in various patient groups characterized by clinical subtype, Expanded Disability Status Scale score and disease duration. We also measured absolute and normalized lesion volumes in several regions of interest using an atlas-based approach, and evaluated differences within and between groups. The lateral funiculi were more frequently affected by lesions in progressive subtypes than in relapsing in voxelwise analysis (P < 0.001), which was further confirmed by absolute and normalized lesion volumes (P < 0.01). The central cord area was more often affected by lesions in primary progressive than relapse-remitting patients (P < 0.001). Between white and grey matter, the absolute lesion volume in the white matter was greater than in the grey matter in all phenotypes (P < 0.001); however when normalizing by each region, normalized lesion volumes were comparable between white and grey matter in primary progressive patients. Lesions appearing in the lateral funiculi and central cord area were significantly correlated with Expanded Disability Status Scale score (P < 0.001). High lesion frequencies were observed in patients with a more aggressive disease course, rather than long disease duration. Lesions located in the lateral funiculi and central cord area of the cervical spine may influence clinical status in multiple sclerosis. This work shows the added value of cervical spine lesions, and provides an avenue for evaluating the distribution of spinal cord lesions in various patient groups.

Generalised boundary shift integral for longitudinal assessment of spinal cord atrophy

Spinal cord atrophy measurements obtained from structural magnetic resonance imaging (MRI) are associated with disability in many neurological diseases and serve as in vivo biomarkers of neurodegeneration. Longitudinal spinal cord atrophy rate is commonly determined from the numerical difference between two volumes (based on 3D surface fitting) or two cross-sectional areas (CSA, based on 2D edge detection) obtained at different time-points. Being an indirect measure, atrophy rates are susceptible to variable segmentation errors at the edge of the spinal cord. To overcome those limitations, we developed a new registration-based pipeline that measures atrophy rates directly. We based our approach on the generalised boundary shift integral (GBSI) method, which registers 2 scans and uses a probabilistic XOR mask over the edge of the spinal cord, thereby measuring atrophy more accurately than segmentation-based techniques. Using a large cohort of longitudinal spinal cord images (610 subjects with multiple sclerosis from a multi-centre trial and 52 healthy controls), we demonstrated that GBSI is a sensitive, quantitative and objective measure of longitudinal spinal cord volume change. The GBSI pipeline is repeatable, reproducible, and provides more precise measurements of longitudinal spinal cord atrophy than segmentation-based methods in longitudinal spinal cord atrophy studies.