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

Recent developments and future avenues for human corticospinal neuroimaging

Non-invasive neuroimaging serves as a valuable tool for investigating the mechanisms within the central nervous system (CNS) related to somatosensory and motor processing, emotions, memory, cognition, and other functions. Despite the extensive use of brain imaging, spinal cord imaging has received relatively less attention, regardless of its potential to study peripheral communications with the brain and the descending corticospinal systems. To comprehensively understand the neural mechanisms underlying human sensory and motor functions, particularly in pathological conditions, simultaneous examination of neuronal activity in both the brain and spinal cord becomes imperative. Although technically demanding in terms of data acquisition and analysis, a growing but limited number of studies have successfully utilized specialized acquisition protocols for corticospinal imaging. These studies have effectively assessed sensorimotor, autonomic, and interneuronal signaling within the spinal cord, revealing interactions with cortical processes in the brain. In this mini-review, we aim to examine the expanding body of literature that employs cutting-edge corticospinal imaging to investigate the flow of sensorimotor information between the brain and spinal cord. Additionally, we will provide a concise overview of recent advancements in functional magnetic resonance imaging (fMRI) techniques. Furthermore, we will discuss potential future perspectives aimed at enhancing our comprehension of large-scale neuronal networks in the CNS and their disruptions in clinical disorders. This collective knowledge will aid in refining combined corticospinal fMRI methodologies, leading to the development of clinically relevant biomarkers for conditions affecting sensorimotor processing in the CNS.

Evaluating tissue injury in cervical spondylotic myelopathy with spinal cord MRI: a systematic review

PURPOSE

Cervical Spondylotic Myelopathy (CSM) is a degenerative condition that leads to loss of cervical spinal cord (CSC) integrity. Various spinal cord Magnetic Resonance Imaging (MRI) methods can identify and characterize the extent of this damage. This systematic review aimed to evaluate the diagnostic, biomarker, and predictive utilities of different spinal cord MRI methods in clinical research studies of CSM. The aim was to provide a comprehensive understanding of the progress in this direction for future studies and effective diagnosis and management of CSM.

METHODS

A comprehensive literature search was conducted on PubMed and EMBASE from 2010 to 2022 according to PRISMA guidelines. Studies with non-human subjects, less than 3T magnetic field strength, non-clinical design, or not quantitatively focusing on the structural integrity of CSC were excluded. The extracted data from each study included demographics, disease severity, MRI machine characteristics, quantitative metrics, and key findings in terms of diagnostic, biomarker, and predictive utilities of each MRI method. The risk of bias was performed using the guide from AHRQ. The quality of evidence was assessed separately for each type of utility for different MRI methods using GRADE.

RESULTS

Forty-seven studies met the inclusion criteria, utilizing diffusion-weighted imaging (DTI) (n = 39), magnetization transfer (MT) (n = 6), MR spectroscopy (n = 3), and myelin water imaging (n = 1), as well as a combination of MRI methods (n = 12). The metric fractional anisotropy (FA) showed the highest potential in all facets of utilities, followed by mean diffusivity. Other promising metrics included MT ratio and intracellular volume fraction, especially in multimodal studies. However, the level of evidence for these promising metrics was low due to a small number of studies. Some studies, mainly DTI, also reported the usefulness of spinal cord MRI in mild CSM.

CONCLUSIONS

Spinal cord MRI methods can potentially facilitate the diagnosis and management of CSM by quantitatively interrogating the structural integrity of CSC. DTI is the most promising MRI method, and other techniques have also shown promise, especially in multimodal configurations. However, this field is in its early stages, and more studies are needed to establish the usefulness of spinal cord MRI in CSM.

Application of the T1w/T2w mapping technique for spinal cord assessment in patients with degenerative cervical myelopathy

STUDY DESIGN

Retrospective case-control study.

OBJECTIVES

To investigate signal changes on T1w/T2w signal intensity ratio maps within cervical cord in patients with degenerative cervical myelopathy (DCM).

SETTING

Novosibirsk Neurosurgery Center, Russia.

METHODS

A total of 261 patients with DCM and 42 age- and sex-matched healthy controls were evaluated using the T1w/T2w mapping method and spinal cord automatic morphometry. The T1w/T2w signal intensity ratio, which reflects white matter integrity, and the spinal cord cross-sectional area (CSA) were calculated and compared between the patients and the controls. In patients with DCM, the correlations between these parameters and neurological scores were also evaluated.

RESULTS

The regional T1w/T2w ratio values from the cervical spinal cord at the level of maximal compression in patients with DCM were significantly lower than those in healthy controls (p < 0.001), as were the regional CSA values (p < 0.001). There was a positive correlation between the regional values of the T1w/T2w ratio and the values of the CSA at the level of maximal spinal cord compression.

CONCLUSIONS

T1w/T2w mapping revealed that spinal cord tissue damage exists at the level of maximal compression in patients with DCM in association with spinal cord atrophy according to automatic morphometry. These changes were correlated with each other.

Investigation of perfusion impairment in degenerative cervical myelopathy beyond the site of cord compression

The aim of this study was to determine tissue-specific blood perfusion impairment of the cervical cord above the compression site in patients with degenerative cervical myelopathy (DCM) using intravoxel incoherent motion (IVIM) imaging. A quantitative MRI protocol, including structural and IVIM imaging, was conducted in healthy controls and patients. In patients, T2-weighted scans were acquired to quantify intramedullary signal changes, the maximal canal compromise, and the maximal cord compression. T2*-weighted MRI and IVIM were applied in all participants in the cervical cord (covering C1-C3 levels) to determine white matter (WM) and grey matter (GM) cross-sectional areas (as a marker of atrophy), and tissue-specific perfusion indices, respectively. IVIM imaging resulted in microvascular volume fraction ([Formula: see text]), blood velocity ([Formula: see text]), and blood flow ([Formula: see text]) indices. DCM patients additionally underwent a standard neurological clinical assessment. Regression analysis assessed associations between perfusion parameters, clinical outcome measures, and remote spinal cord atrophy. Twenty-nine DCM patients and 30 healthy controls were enrolled in the study. At the level of stenosis, 11 patients showed focal radiological evidence of cervical myelopathy. Above the stenosis level, cord atrophy was observed in the WM (- 9.3%; p = 0.005) and GM (- 6.3%; p = 0.008) in patients compared to healthy controls. Blood velocity (BV) and blood flow (BF) indices were decreased in the ventral horns of the GM (BV: - 20.1%, p = 0.0009; BF: - 28.2%, p = 0.0008), in the ventral funiculi (BV: - 18.2%, p = 0.01; BF: - 21.5%, p = 0.04) and lateral funiculi (BV: - 8.5%, p = 0.03; BF: - 16.5%, p = 0.03) of the WM, across C1-C3 levels. A decrease in microvascular volume fraction was associated with GM atrophy (R = 0.46, p = 0.02). This study demonstrates tissue-specific cervical perfusion impairment rostral to the compression site in DCM patients. IVIM indices are sensitive to remote perfusion changes in the cervical cord in DCM and may serve as neuroimaging biomarkers of hemodynamic impairment in future studies. The association between perfusion impairment and cervical cord atrophy indicates that changes in hemodynamics caused by compression may contribute to the neurodegenerative processes in DCM.

Association between activity in the ventral premotor cortex and spinal cord activation during force generation-A combined cortico-spinal fMRI study

Force generation is a crucial element of dexterity and a highly relevant skill of the human motor system. How cerebral and spinal components interact and how spinal activation is associated with the activity in the cerebral primary motor and premotor areas is poorly understood. Here, we conducted combined cortico-spinal functional magnetic resonance imaging during a simple visually guided isometric force generation task in 20 healthy young subjects. Activation was localized in the right cervical spinal cord and left primary motor and premotor areas. The main finding is that spinal activation was negatively correlated with ventral premotor cortex activation. Spinal activation was furthermore significantly correlated with primary motor cortex activation, while increasing target forces led to an increase in the amount of activation. These data indicate that human premotor areas such as the ventral premotor cortex might be functionally connected to the lower cervical spinal cord contributing to distal upper limb functions, a finding that extends our understanding of human motor function beyond the animal literature.

Spatial distribution of hand-grasp motor task activity in spinal cord functional magnetic resonance imaging

Upper extremity motor paradigms during spinal cord functional magnetic resonance imaging (fMRI) can provide insight into the functional organization of the cord. Hand-grasping is an important daily function with clinical significance, but previous studies of similar squeezing movements have not reported consistent areas of activity and are limited by sample size and simplistic analysis methods. Here, we study spinal cord fMRI activation using a unimanual isometric hand-grasping task that is calibrated to participant maximum voluntary contraction (MVC). Two task modeling methods were considered: (1) a task regressor derived from an idealized block design (Ideal) and (2) a task regressor based on the recorded force trace normalized to individual MVC (%MVC). Across these two methods, group motor activity was highly lateralized to the hemicord ipsilateral to the side of the task. Activation spanned C5-C8 and was primarily localized to the C7 spinal cord segment. Specific differences in spatial distribution are also observed, such as an increase in C8 and dorsal cord activity when using the %MVC regressor. Furthermore, we explored the impact of data quantity and spatial smoothing on sensitivity to hand-grasp motor task activation. This analysis shows a large increase in number of active voxels associated with the number of fMRI runs, sample size, and spatial smoothing, demonstrating the impact of experimental design choices on motor activation.

Spinal Cord fMRI: A New Window into the Central Nervous System

With the brain, the spinal cord forms the central nervous system. Initially considered a passive relay between the brain and the periphery, the spinal cord is now recognized as being active and plastic. Yet, it remains largely overlooked by the human neuroscience community, in stark contrast with the wealth of research investigating the brain. In this review, we argue that fMRI, traditionally used to image cerebral function, can be extended beyond the brain to help unravel spinal mechanisms involved in human behaviors. To this end, we first outline strategies that have been proposed to tackle the challenges inherent to spinal cord fMRI. Then, we discuss how they have been utilized to provide insights into the functional organization of spinal sensorimotor circuits, highlighting their potential to address fundamental and clinical questions. By summarizing guidelines and applications of spinal cord fMRI, we hope to stimulate and support further research into this promising yet underexplored field.

Structural and Physiological Modeling (SAPM) for the Analysis of Functional MRI Data Applied to a Study of Human Nociceptive Processing

A novel method has been developed for analyzing connectivity between regions based on functional magnetic resonance imaging (fMRI) data. This method, termed structural and physiological modeling (SAPM), combines information about blood oxygenation-level dependent (BOLD) responses, anatomy, and physiology to model coordinated signaling across networks of regions, including input and output signaling from each region and whether signaling is predominantly inhibitory or excitatory. The present study builds on a prior proof-of-concept demonstration of the SAPM method by providing evidence for the choice of network model and anatomical sub-regions, demonstrating the reproducibility of the results and identifying statistical thresholds needed to infer significance. The method is further validated by applying it to investigate human nociceptive processing in the brainstem and spinal cord and comparing the results to the known neuroanatomy, including anatomical regions and inhibitory and excitatory signaling. The results of this analysis demonstrate that it is possible to obtain reliable information about input and output signaling from anatomical regions and to identify whether this signaling has predominantly inhibitory or excitatory effects. SAPM provides much more detailed information about neuroanatomy than was previously possible based on fMRI data.

Thermal Stimulus Task fMRI in the Cervical Spinal Cord at 7 Tesla

PURPOSE

Although functional MRI is widely applied in the brain, fMRI of the spinal cord is more technically demanding. Proximity to the vertebral column and lungs results in strong spatial inhomogeneity and temporal fluctuations in B0. Increasing field strength enables higher spatial resolution and improved sensitivity to BOLD signal, but amplifies the effects of B0 inhomogeneity. In this work, we present the first stimulus task fMRI in the spinal cord at 7 T. Further, we compare the performance of single-shot and multi-shot 2D EPI protocols, as they differ in sensitivity to spatial and temporal B0 inhomogeneity.

METHODS

The cervical spinal cords of 11 healthy volunteers were scanned at 7 T using single-shot 2D EPI at 0.75 mm in-plane resolution and multi-shot 2D EPI at 0.75 and 0.6 mm in-plane resolutions. For each protocol, the BOLD response to thirteen 10-second noxious thermal stimuli applied to the right thumb was acquired in a 10-minute fMRI run. Image quality, temporal SNR, and BOLD activation (percent signal change and z-stat) at both individual- and group-level were evaluated between the protocols.

RESULTS

Temporal SNR was highest in single-shot and multi-shot 0.75 mm protocols. In group-level analyses, activation clusters appeared in all protocols in the ipsilateral dorsal quadrant at the expected C6 neurological level. In individual-level analyses, activation clusters at the expected level were detected in some, but not all subjects and protocols. Single-shot 0.75 mm generally produced the highest mean z-statistic, while multi-shot 0.60 mm produced the best-localized activation clusters and the least geometric distortion. Larger than expected within-subject segmental variation of BOLD activation along the cord was observed.

CONCLUSION

Group-level sensory task fMRI of the cervical spinal cord is feasible at 7 T with single-shot or multi-shot EPI. The best choice of protocol will likely depend on the relative importance of sensitivity to activation versus spatial localization of activation for a given experiment.

The CALIPR framework for highly accelerated myelin water imaging with improved precision and sensitivity

Quantitative magnetic resonance imaging (MRI) techniques are powerful tools for the study of human tissue, but, in practice, their utility has been limited by lengthy acquisition times. Here, we introduce the Constrained, Adaptive, Low-dimensional, Intrinsically Precise Reconstruction (CALIPR) framework in the context of myelin water imaging (MWI); a quantitative MRI technique generally regarded as the most rigorous approach for noninvasive, in vivo measurement of myelin content. The CALIPR framework exploits data redundancy to recover high-quality images from a small fraction of an imaging dataset, which allowed MWI to be acquired with a previously unattainable sequence (fully sampled acquisition 2 hours:57 min:20 s) in 7 min:26 s (4.2% of the dataset, acceleration factor 23.9). CALIPR quantitative metrics had excellent precision (myelin water fraction mean coefficient of variation 3.2% for the brain and 3.0% for the spinal cord) and markedly increased sensitivity to demyelinating disease pathology compared to a current, widely used technique. The CALIPR framework facilitates drastically improved MWI and could be similarly transformative for other quantitative MRI applications.

The current state of spinal cord functional magnetic resonance imaging and its application in clinical research

Since its development, spinal cord functional magnetic resonance imaging (fMRI) has utilized various methodologies and stimulation protocols to develop a deeper understanding of a healthy human spinal cord that lays a foundation for its use in clinical research and practice. In this review, we conducted a comprehensive literature search on spinal cord fMRI studies and summarized the recent advancements and resulting scientific achievements of spinal cord fMRI in the following three aspects: the current state of spinal cord fMRI methodologies and stimulation protocols, knowledge about the healthy spinal cord's functions obtained via spinal cord fMRI, and fMRI's exemplary usage in spinal cord diseases and injuries. We conclude with a discussion that, while technical challenges exist, novel fMRI technologies for and new knowledge about the healthy human spinal cord have been established. Empowered by these developments, investigations of pathological and injury states within the spinal cord have become the next important direction of spinal cord fMRI. Recent clinical investigations into spinal cord pathologies, for example, fibromyalgia, multiple sclerosis, spinal cord injury, and cervical spondylotic myelopathy, have already provided deep insights into spinal cord impairments and the time course of impairment-caused changes. We expect that future spinal cord fMRI advancement and research development will further enhance our understanding of various spinal cord diseases and provide the foundation for evaluating existing and developing new treatment plans.

An 8-channel Tx dipole and 20-channel Rx loop coil array for MRI of the cervical spinal cord at 7 Tesla

The quality of cervical spinal cord images can be improved by the use of tailored radiofrequency (RF) coil solutions for ultrahigh field imaging; however, very few commercial and research 7-T RF coils currently exist for the spinal cord, and in particular, those with parallel transmission (pTx) capabilities. This work presents the design, testing, and validation of a pTx/Rx coil for the human neck and cervical/upper thoracic spinal cord. The pTx portion is composed of eight dipoles to ensure high homogeneity over this large region of the spinal cord. The Rx portion is made up of twenty semiadaptable overlapping loops to produce high signal-to-noise ratio (SNR) across the patient population. The coil housing is designed to facilitate patient positioning and comfort, while also being tight fitting to ensure high sensitivity. We demonstrate RF shimming capabilities to optimize B1 + uniformity, power efficiency, and/or specific absorption rate efficiency. B1 + homogeneity, SNR, and g-factor were evaluated in adult volunteers and demonstrated excellent performance from the occipital lobe down to the T4-T5 level. We compared the proposed coil with two state-of-the-art head and head/neck coils, confirming its superiority in the cervical and upper thoracic regions of the spinal cord. This coil solution therefore provides a convincing platform for producing the high image quality necessary for clinical and research scanning of the upper spinal cord.

A spinal cord neuroprosthesis for locomotor deficits due to Parkinson's disease

People with late-stage Parkinson's disease (PD) often suffer from debilitating locomotor deficits that are resistant to currently available therapies. To alleviate these deficits, we developed a neuroprosthesis operating in closed loop that targets the dorsal root entry zones innervating lumbosacral segments to reproduce the natural spatiotemporal activation of the lumbosacral spinal cord during walking. We first developed this neuroprosthesis in a non-human primate model that replicates locomotor deficits due to PD. This neuroprosthesis not only alleviated locomotor deficits but also restored skilled walking in this model. We then implanted the neuroprosthesis in a 62-year-old male with a 30-year history of PD who presented with severe gait impairments and frequent falls that were medically refractory to currently available therapies. We found that the neuroprosthesis interacted synergistically with deep brain stimulation of the subthalamic nucleus and dopaminergic replacement therapies to alleviate asymmetry and promote longer steps, improve balance and reduce freezing of gait. This neuroprosthesis opens new perspectives to reduce the severity of locomotor deficits in people with PD.

Effectiveness of regional diffusion MRI measures in distinguishing multiple sclerosis abnormalities within the cervical spinal cord

INTRODUCTION

Multiple sclerosis (MS) is an inflammatory disorder of the central nervous system. Although conventional magnetic resonance imaging (MRI) is widely used for MS diagnosis and clinical follow-up, quantitative MRI has the potential to provide valuable intrinsic values of tissue properties that can enhance accuracy. In this study, we investigate the efficacy of diffusion MRI in distinguishing MS lesions within the cervical spinal cord, using a combination of metrics extracted from diffusion tensor imaging and Ball-and-Stick models.

METHODS

We analyzed spinal cord data acquired from multiple hospitals and extracted average diffusion MRI metrics per vertebral level using a collection of image processing methods and an atlas-based approach. We then performed a statistical analysis to evaluate the feasibility of these metrics for detecting lesions, exploring the usefulness of combining different metrics to improve accuracy.

RESULTS

Our study demonstrates the sensitivity of each metric to underlying microstructure changes in MS patients. We show that selecting a specific subset of metrics, which provide complementary information, significantly improves the prediction score of lesion presence in the cervical spinal cord. Furthermore, the Ball-and-Stick model has the potential to provide novel information about the microstructure of damaged tissue.

CONCLUSION

Our results suggest that diffusion measures, particularly combined measures, are sensitive in discriminating abnormal from healthy cervical vertebral levels in patients. This information could aid in improving MS diagnosis and clinical follow-up. Our study highlights the potential of the Ball-and-Stick model in providing additional insights into the microstructure of the damaged tissue.

A systematic review of repeatability and reproducibility studies of diffusion tensor imaging of cervical spinal cord

OBJECTIVES

Diffusion tensor imaging (DTI) techniques are being studied as a possible diagnostic and predictive tool for the evaluation of cervical spinal cord disease. This systematic review aims to evaluate the previous DTI studies that specifically investigated the repeatability and reproducibility of DTI in the cervical spinal cord.

METHODS AND MATERIALS

A search in the PubMed, Scopus, Web of Science and Ovid electronic databases was conducted for articles published between January 1990 and February 2022 that related to the repeatability and reproducibility of DTI in evaluating the cervical spinal cord using one of the following measurements: the intraclass correlation coefficient (ICC) and/or the coefficient of variation (CV), and/or Bland-Altman (BA) differences analysis methods. DTI studies that presented full statistical analysis of repeatability and/or reproducibility tests of the cervical spinal cord in peer-reviewed full-text publications published in journals were included. Articles that included at least one of the keywords within the titles or abstracts were identified. Additional full-text papers were found by searching the citations and reference lists of related articles. This review has followed the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidance. Risk of bias was evaluated with 13 criteria weighted toward methodological quality of reported studies using the QuADS assessment criteria. This assessment only included full-text articles written in English.

RESULTS

A total of 11 studies were included and assessed for different characteristics, including sample size,(3-34) re-test time interval (<1 h to >3 months), test-retest reproducibility scores and acquisition method. Six studies used ICC which ranged from poor (ICC<0.37) to excellent reproducibility (ICC 0.91-0.99). Four studies reported an overall CV lower than 40% for all DTI metrics. Three studies reported the Bland-Altman (BA) differences and reported a minimum percentage showing no strong differences between repeated measurements. Quantitative analysis was not undertaken due to heterogeneity of methods. Repeatability and reproducibility measures were generally found to be good.

CONCLUSION

This study revealed that the application of DTI and its related measures in a clinical setting in the assessment of cervical spinal cord changes is feasible and reproducible. However, cervical spinal cord DTI suffers from some existing limitations that prevent it from being routinely used in research and clinical settings.

ADVANCES IN KNOWLEDGE

DTI with its parametric maps provide broad evaluation of the tissue structure of axonal white matter and are being studied as a possible diagnostic and predictive tool for the assessment of cervical spinal cord (CSC) disease.

Detection of resting-state functional connectivity in the lumbar spinal cord with 3T MRI

Functional MRI (fMRI) of the spinal cord is an expanding area of research with potential to investigate neuronal activity in the central nervous system. We aimed to characterize the functional connectivity features of the human lumbar spinal cord using resting-state fMRI (rs-fMRI) at 3T, using region-based and data-driven analysis approaches. A 3D multi-shot gradient echo resting-state blood oxygenation level dependent-sensitive rs-fMRI protocol was implemented in 26 healthy participants. Average temporal signal-to-noise ratio in the gray matter was 16.35 ± 4.79 after denoising. Evidence of synchronous signal fluctuations in the ventral and dorsal horns with their contralateral counterparts was observed in representative participants using interactive, exploratory seed-based correlations. Group-wise average in-slice Pearson's correlations were 0.43 ± 0.17 between ventral horns, and 0.48 ± 0.16 between dorsal horns. Group spatial independent component analysis (ICA) was used to identify areas of coherent activity¸ and revealed components within the gray matter corresponding to anatomical regions. Lower-dimensionality ICA revealed bilateral components corresponding to ventral and dorsal networks. Additional separate ICAs were run on two subsets of the participant group, yielding two sets of components that showed visual consistency and moderate spatial overlap. This work shows feasibility of rs-fMRI to probe the functional features and organization of the lumbar spinal cord.

Detecting macromolecular differences of the CSF in low disability multiple sclerosis using quantitative MT MRI at 3T

Background

Imaging investigation of cerebrospinal fluid (CSF) in multiple sclerosis (MS) is understudied. Development of noninvasive methods to detect pathological CSF changes would have a profound effect on MS diagnosis and would offer insight into MS pathophysiology and mechanisms of neurological impairment.

Objective

We propose magnetization transfer (MT) MRI as a tool to detect macromolecular changes in spinal CSF.

Methods

MT and quantitative MT (qMT) data were acquired in the cervical region in 27 people with relapsing-remitting multiple sclerosis (pwRRMS) and 38 age and sex-matched healthy controls (HCs). MT ratio (MTR), the B1, B0, and R1 corrected qMT-derived pool size ratio (PSR) were quantified in the spinal cord and CSF of each group.

Results

Both CSF MTR and CSF qMT-derived PSR were significantly increased in pwRRMS compared to HC (p = 0.027 and p = 0.020, respectively). CSF PSR of pwRRMS was correlated to Expanded Disability Status Scale Scores (p = 0.045, R = 0.352).

Conclusion

Our findings demonstrate increased CSF macromolecular content in pwRRMS and link CSF macromolecular content with clinical impairment. This highlights the potential role of CSF in processing products of demyelination.

Advanced Structural Magnetic Resonance Imaging of the Spinal Cord: Technical Aspects and Clinical Use

For a long time, technical obstacles have hampered the acquisition of high-resolution images and the development of reliable processing protocols for spinal cord (SC) MRI. Fortunately, this scenario has changed in the past 5-10 years, due to hardware and software improvements. Nowadays, with advanced protocols, SC MRI is considered a useful tool for several inherited and acquired neurologic diseases, not only for diagnosis approach but also for pathophysiological unraveling and as a biomarker for disease monitoring and clinical trials. In this review, we address advanced SC MRI sequences for macrostructural and microstructural evaluation, useful semiautomatic and automatic processing tools and clinical applications on several neurologic conditions such as hereditary cerebellar ataxia, hereditary spastic paraplegia, motor neuron diseases and multiple sclerosis.

Attention-gated U-Net networks for simultaneous axial/sagittal planes segmentation of injured spinal cords

Magnetic resonance imaging is currently the gold standard for the evaluation of spinal cord injuries. Automatic analysis of these injuries is however challenging, as MRI resolutions vary for different planes of analysis and physiological features are often distorted around these injuries. This study proposes a new CNN-based segmentation method in which information is exchanged between two networks analyzing the scans from different planes. Our aim was to develop a robust method for automatic segmentation of the spinal cord in patients having suffered traumatic injuries. The database consisted of 106 sagittal MRI scans from 94 patients with traumatic spinal cord injuries. Our method used an innovative approach where the scans were analyzed in series under the axial and sagittal plane by two different convolutional networks. The results were compared with those of Deepseg 2D from the Spinal Cord Toolbox (SCT), which was taken as state-of-the-art. Comparisons were evaluated using K-Fold cross-validation combined with statistical t-test results on separate test data. Our method achieved significantly better results than Deepseg 2D, with an average Dice coefficient of 0.95 against 0.88 for Deepseg 2D (p <0.001). Other metrics were also used to compare the segmentations, all of which showed significantly better results for our approach. In this study, we introduce a robust method for spinal cord segmentation which is capable of adequately segmenting spinal cords affected by traumatic injuries, improving upon the methods contained in SCT.

Brainstem and subcortical regions volume loss in patients with degenerative cervical myelopathy and its association with spinal cord compression severity

BACKGROUND

In recent years, structural and functional reorganization of the brain and changes in brainstem structural connectivity have been shown in patients with degenerative cervical myelopathy (DCM). We hypothesized that volume loss in the basal ganglia, thalami, and brainstem structures exists and is associated with spinal cord compression severity in patients with DCM.

METHODS

Forty-seven patients with DCM and 25 patients with cervical radiculopathy were evaluated using cervical spinal cord and brain magnetic resonance imaging (MRI). Brainstem structures, basal ganglia, and thalami volumes were evaluated with FreeSurfer and compared between groups with correction for individual intracranial volume, as well as patient age and sex. Additionally, spinal cord MRI data were analysed with the Spinal Cord Toolbox, and cross-sectional area (CSA) and fractional anisotropy (FA) values were calculated. Correlations between MR-morphometry data and spinal cord structural changes, as well as disease duration, were also evaluated in patients with DCM.

RESULTS

A statistically significant reduction in the volume of the whole brainstem was revealed in the DCM group compared to the radiculopathy group (p < 0.01, FDR-corrected). Additionally, reductions in medulla oblongata, pons and midbrain volumes were found in patients with DCM (p < 0.01, p < 0.01 and p < 0.05, respectively, FDR-corrected). Additionally, a trend in the loss of volume of the left putamen was found (p = 0.087, FDR-corrected). Furthermore, medulla oblongata volume was correlated with spinal cord compression severity (R = 0.54, adjusted p < 0.001) and white matter damage (R = 0.46, adjusted p < 0.05) in patients with DCM. Negative correlations between the duration of the disease and the severity of spinal cord compression (R = -0.42, adjusted p < 0.05) and white matter damage (R = -0.49, adjusted p < 0.05) were also revealed, as well as a trend toward a negative association between the duration of the disease and the volume of the medulla oblongata (R = -0.35; adjusted p < 0.1).

CONCLUSIONS

We revealed a reduction in the volume of brainstem structures in patients with DCM compared to patients with radiculopathy. Moreover, we found that these changes are associated with cord compression severity.

Longitudinal multiparametric MRI of traumatic spinal cord injury in animal models

Multi-parametric MRI (mpMRI) technology enables non-invasive and quantitative assessments of the structural, molecular, and functional characteristics of various neurological diseases. Despite the recognized importance of studying spinal cord pathology, mpMRI applications in spinal cord research have been somewhat limited, partly due to technical challenges associated with spine imaging. However, advances in imaging techniques and improved image quality now allow longitudinal investigations of a comprehensive range of spinal cord pathological features by exploiting different endogenous MRI contrasts. This review summarizes the use of mpMRI techniques including blood oxygenation level-dependent (BOLD) functional MRI (fMRI), diffusion tensor imaging (DTI), quantitative magnetization transfer (qMT), and chemical exchange saturation transfer (CEST) MRI in monitoring different aspects of spinal cord pathology. These aspects include cyst formation and axonal disruption, demyelination and remyelination, changes in the excitability of spinal grey matter and the integrity of intrinsic functional circuits, and non-specific molecular changes associated with secondary injury and neuroinflammation. These approaches are illustrated with reference to a nonhuman primate (NHP) model of traumatic cervical spinal cord injuries (SCI). We highlight the benefits of using NHP SCI models to guide future studies of human spinal cord pathology, and demonstrate how mpMRI can capture distinctive features of spinal cord pathology that were previously inaccessible. Furthermore, the development of mechanism-based MRI biomarkers from mpMRI studies can provide clinically useful imaging indices for understanding the mechanisms by which injured spinal cords progress and repair. These biomarkers can assist in the diagnosis, prognosis, and evaluation of therapies for SCI patients, potentially leading to improved outcomes.

Multiparametric magnetic resonance imaging for detection of pathological changes in the central nervous system of a mouse model of multiple sclerosis in vivo

Multiple sclerosis (MS) is an autoimmune disease involving demyelination and axonal damage in the central nervous system (CNS). In this study, we investigated pathological changes in the lumbar spinal cord of C57BL/6 mice induced with progressive experimental autoimmune encephalomyelitis (EAE) disease using 9.4-T magnetic resonance imaging (MRI). Multiparametric MRI measurements including MR spectroscopy, diffusion tensor imaging (DTI) and volumetric analyses were applied to detect metabolic changes in the CNS of EAE mice. Compared with healthy mice, EAE mice showed a significant reduction in N-acetyl aspartate and increases in choline, glycine, taurine and lactate. DTI revealed a significant reduction in fractional anisotropy and axial diffusivity and an increase in radial diffusivity in the lumbar spinal cord white matter (WM), while in the grey matter (GM), fractional anisotropy increased. High-resolution structural imaging also revealed lumbar spinal cord WM hypertrophy and GM atrophy. Importantly, these MRI changes were strongly correlated with EAE disease scoring and pathological changes in the lumbar (L2-L6), particularly WM demyelination lesions and aggregation of immune cells (microglia/macrophages and astrocytes) in this region. This study identified changes in MRI biomarker signatures that can be useful for evaluating the efficacy of novel drugs using EAE models in vivo.

The use of deep learning in medical imaging to improve spine care: A scoping review of current literature and clinical applications

Background

Artificial intelligence is a revolutionary technology that promises to assist clinicians in improving patient care. In radiology, deep learning (DL) is widely used in clinical decision aids due to its ability to analyze complex patterns and images. It allows for rapid, enhanced data, and imaging analysis, from diagnosis to outcome prediction. The purpose of this study was to evaluate the current literature and clinical utilization of DL in spine imaging.

Methods

This study is a scoping review and utilized the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) methodology to review the scientific literature from 2012 to 2021. A search in PubMed, Web of Science, Embased, and IEEE Xplore databases with syntax specific for DL and medical imaging in spine care applications was conducted to collect all original publications on the subject. Specific data was extracted from the available literature, including algorithm application, algorithms tested, database type and size, algorithm training method, and outcome of interest.

Results

A total of 365 studies (total sample of 232,394 patients) were included and grouped into 4 general applications: diagnostic tools, clinical decision support tools, automated clinical/instrumentation assessment, and clinical outcome prediction. Notable disparities exist in the selected algorithms and the training across multiple disparate databases. The most frequently used algorithms were U-Net and ResNet. A DL model was developed and validated in 92% of included studies, while a pre-existing DL model was investigated in 8%. Of all developed models, only 15% of them have been externally validated.

Conclusions

Based on this scoping review, DL in spine imaging is used in a broad range of clinical applications, particularly for diagnosing spinal conditions. There is a wide variety of DL algorithms, database characteristics, and training methods. Future studies should focus on external validation of existing models before bringing them into clinical use.

The longitudinal progression of MRI changes in pre-ataxic carriers of SCA3/MJD

BACKGROUND

The natural history of magnetic resonance imaging (MRI) in pre-ataxic stages of spinocerebellar ataxia type 3/Machado-Joseph disease (SCA3/MJD) is not well known. We report cross-sectional and longitudinal data obtained at this stage.

METHODS

Baseline (follow-up) observations included 32 (17) pre-ataxic carriers (SARA < 3) and 20 (12) related controls. The mutation length was used to estimate the time to onset (TimeTo) of gait ataxia. Clinical scales and MRIs were performed at baseline and after a median (IQR) of 30 (7) months. Cerebellar volumetries (ACAPULCO), deep gray-matter (T1-Multiatlas), cortical thickness (FreeSurfer), cervical spinal cord area (SCT) and white matter (DTI-Multiatlas) were assessed. Baseline differences between groups were described; variables that presented a p < 0.1 after Bonferroni correction were assessed longitudinally, using TimeTo and study time. For TimeTo strategy, corrections for age, sex and intracranial volume were done with Z-score progression. A significance level of 5% was adopted.

RESULTS

SCT at C1 level distinguished pre-ataxic carriers from controls. DTI measures of the right inferior cerebellar peduncle (ICP), bilateral middle cerebellar peduncles (MCP) and bilateral medial lemniscus (ML), also distinguished pre-ataxic carriers from controls, and progressed over TimeTo, with effect sizes varying from 0.11 to 0.20, larger than those of the clinical scales. No MRI variable showed progression over study time.

DISCUSSION

DTI parameters of the right ICP, left MCP and right ML were the best biomarkers for the pre-ataxic stage of SCA3/MJD. TimeTo is an interesting timescale, since it captured the longitudinal worsening of these structures.

Pontomedullary junction as a reference for spinal cord cross-sectional area: validation across neck positions

Spinal cord cross-sectional area (CSA) is an important MRI biomarker to assess spinal cord atrophy in various neurodegenerative and traumatic spinal cord diseases. However, the conventional method of computing CSA based on vertebral levels is inherently flawed, as the prediction of spinal levels from vertebral levels lacks reliability, leading to considerable variability in CSA measurements. Computing CSA from an intrinsic neuroanatomical reference, the pontomedullary junction (PMJ), has been proposed in previous work to overcome limitations associated with using a vertebral reference. However, the validation of this alternative approach, along with its variability across and within participants under variable neck extensions, remains unexplored. The goal of this study was to determine if the variability of CSA across neck flexions/extensions is reduced when using the PMJ, compared to vertebral levels. Ten participants underwent a 3T MRI T2w isotropic scan at 0.6 mm3 for 3 neck positions: extension, neutral and flexion. Spinal cord segmentation, vertebral labeling, PMJ labeling, and CSA were computed automatically while spinal segments were labeled manually. Mean coefficient of variation for CSA across neck positions was 3.99 ± 2.96% for the PMJ method vs. 4.02 ± 3.01% for manual spinal segment method vs. 4.46 ± 3.10% for the disc method. These differences were not statistically significant. The PMJ method was slightly more reliable than the disc-based method to compute CSA at specific spinal segments, although the difference was not statistically significant. This suggests that the PMJ can serve as a valuable alternative and reliable method for estimating CSA when a disc-based approach is challenging or not feasible, such as in cases involving fused discs in individuals with spinal cord injuries.

Multiplex core of the human brain using structural, functional and metabolic connectivity derived from hybrid PET-MR imaging

With the increasing success of mapping brain networks and availability of multiple MR- and PET-based connectivity measures, the need for novel methodologies to unravel the structure and function of the brain at multiple spatial and temporal scales is emerging. Therefore, in this work, we used hybrid PET-MR data of healthy volunteers (n = 67) to identify multiplex core nodes in the human brain. First, monoplex networks of structural, functional and metabolic connectivity were constructed, and consequently combined into a multiplex SC-FC-MC network by linking the same nodes categorically across layers. Taking into account the multiplex nature using a tensorial approach, we identified a set of core nodes in this multiplex network based on a combination of eigentensor centrality and overlapping degree. We introduced a coreness coefficient, which mitigates the effect of modeling parameters to obtain robust results. The proposed methodology was applied onto young and elderly healthy volunteers, where differences observed in the monoplex networks persisted in the multiplex as well. The multiplex core showed a decreased contribution to the default mode and salience network, while an increased contribution to the dorsal attention and somatosensory network was observed in the elderly population. Moreover, a clear distinction in eigentensor centrality was found between young and elderly healthy volunteers.

Spatial distribution of hand-grasp motor task activity in spinal cord functional magnetic resonance imaging

Upper extremity motor paradigms during spinal cord functional magnetic resonance imaging (fMRI) can provide insight into the functional organization of the cord. Hand-grasping is an important daily function with clinical significance, but previous studies of similar squeezing movements have not reported consistent areas of activity and are limited by sample size and simplistic analysis methods. Here, we study spinal cord fMRI activation using a unimanual isometric hand-grasping task that is calibrated to participant maximum voluntary contraction (MVC). Two task modeling methods were considered: (1) a task regressor derived from an idealized block design (Ideal) and (2) a task regressor based on the recorded force trace normalized to individual MVC (%MVC). Across these two methods, group motor activity was highly lateralized to the hemicord ipsilateral to the side of the task. Activation spanned C5-C8 and was primarily localized to the C7 spinal cord segment. Specific differences in spatial distribution are also observed, such as an increase in C8 and dorsal cord activity when using the %MVC regressor. Furthermore, we explored the impact of data quantity and spatial smoothing on sensitivity to hand-grasp motor task activation. This analysis shows a large increase in number of active voxels associated with the number of fMRI runs, sample size, and spatial smoothing, demonstrating the impact of experimental design choices on motor activation.

Validation of a New Semiautomated Segmentation Pipeline Based on the Spinal Cord Toolbox DeepSeg Algorithm to Estimate the Cervical Canal Area

BACKGROUND AND PURPOSE

As in the brain reserve concept, a larger cervical canal area may also protect against disability. In this context, a semiautomated pipeline has been developed to obtain quantitative estimations of the cervical canal area. The aim of the study was to validate the pipeline, to evaluate the consistency of the cervical canal area measurements during a 1-year period, and to compare cervical canal area estimations obtained from brain and cervical MRI acquisitions.

MATERIALS AND METHODS

Eight healthy controls and 18 patients with MS underwent baseline and follow-up 3T brain and cervical spine sagittal 3D MPRAGE. The cervical canal area was measured in all acquisitions, and estimations obtained with the proposed pipeline were compared with manual segmentations performed by 1 evaluator using the Dice similarity coefficient. The cervical canal area estimations obtained on baseline and follow-up T1WI were compared; brain and cervical cord acquisitions were also compared using the individual and average intraclass correlation coefficients.

RESULTS

The agreement between the manual cervical canal area masks and the masks provided by the proposed pipeline was excellent, with a mean Dice similarity coefficient mean of 0.90 (range, 0.73-0.97). The cervical canal area estimations obtained from baseline and follow-up scans showed a good level of concordance (intraclass correlation coefficient = 0.76; 95% CI, 0.44-0.88); estimations obtained from brain and cervical MRIs also had good agreement (intraclass correlation coefficient = 0.77; 95% CI, 0.45-0.90).

CONCLUSIONS

The proposed pipeline is a reliable tool to estimate the cervical canal area. The cervical canal area is a stable measure across time; moreover, when cervical sequences are not available, the cervical canal area could be estimated using brain T1WI.

Effect of Non-invasive Spinal Stimulation on Self-sustained Firing Motoneuron Model: In-Silico Study Using Human Body Model

Transcutaneous spinal electrical stimulation (tSCS) is a non-invasive neuromodulation approach using a low intensity direct current. Recent developments in the technique have opened the possibility that tSCS can help restore motor function after spinal cord injury (SCI). However, the exact mechanism of action tSCS has on the spinal circuits is still unknown. Due to the complexity of experimental synthesis in a human model to delineate the mechanisms, models that link the stimulation paradigm and circuit behaviors are advantageous. Thus, this study aims to simulate the underlying changes in motor circuit firing rates in response to external stimuli induced by tSCS. Serial stimulations combining a high-fidelity finite element model with the human torso and spinal cord with a lumped motor neuron model is constructed. The parameters for both components of the model were derived from previous studies. We focused our analysis on a lumped motor neuron model that describes sustained firing behavior of the motor neuron driven primarily by persistent inward current (PIC), a signature behavior of the motor neuron after SCI. Modulation of the PIC behaviors was achieved by stimulating voltage-dependent calcium and sodium channels in the dendrite using a tSCS-induced electric field (E-field) expressed at different a spatial locations of the motor neuron in the gray matter. The PIC behaviors of spinal motor neurons in the left ventral horn were suppressed, while for the most part invariant in the right ventral horn. These initial simulations will provide a steppingstone for future examinations that incorporate additional neuronal models of inhibitory and excitatory interneurons to access the circuit-level effect of spinal stimulation.

A second-order and slice-specific linear shimming technique to improve spinal cord fMRI

PURPOSE

To develop a second-order and slice-specific linear shimming technique and investigate its efficiency in the mitigation of signal loss and distortions, and the increase of temporal signal-to-noise ratio (tSNR) within the spinal cord during functional Magnetic Resonance Imaging (fMRI) of the human cervical spinal cord.

METHODS

All scans were performed on a General Electric Discovery MR750 3 T scanner, using a head, neck and spine coil and a neurovascular array. To improve B₀ homogeneity, a field map was acquired, and second-order shims (SOS) were optimized over manually defined regions of interest (ROIs). Signal loss from dephasing by susceptibility-induced gradients was reduced by optimizing slice-specific x-, y- and z-shims to maximize signal within the spinal cord. Spectral-spatial excitation pulses were used in both the slice-specific linear shimming calibration scan and fMRI acquisitions. The shimming technique's efficiency was initially tested on eight healthy volunteers by comparing tSNR between images acquired with the manufacturer's standard linear shimming and with our SOS and xyz-shimming technique. Subsequently, using an increased spatial resolution as needed for fMRI of the spinal cord, tSNR measurements were performed on resting-state fMRI images from 14 healthy participants.

RESULTS

Spinal fMRI images acquired with only the standard linear shimming suffered from severe signal loss below the C5 vertebral level. The developed shimming technique compensated for this loss especially at levels C6 and C7, while tSNR was significantly higher at all vertebral levels with SOS and xyz-shimming than without it.

CONCLUSION

A comprehensive shimming approach which includes the use of spectral-spatial excitation pulses along with both second-order and slice-specific linear shim optimization reduces regional signal loss and increases tSNR along the c-spine (C3-C7), improving the ability to record functional signals from the human spinal cord.

Feasibility of spinal cord imaging at 7 T using rosette trajectory with magnetization transfer preparation and compressed sensing

MRI is a valuable diagnostic tool to investigate spinal cord (SC) pathology. SC MRI can benefit from the increased signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) at ultra-high fields such as 7 T. However, SC MRI acquisitions with routine Cartesian readouts are prone to image artifacts caused by physiological motion. MRI acquisition techniques with non-Cartesian readouts such as rosette can help reduce motion artifacts. The purpose of this study was to demonstrate the feasibility of high-resolution SC imaging using rosette trajectory with magnetization transfer preparation (MT-prep) and compressed sensing (CS) at 7 T. Five healthy volunteers participated in the study. Images acquired with rosette readouts demonstrated reduced motion artifacts compared to the standard Cartesian readouts. The combination of multi-echo rosette-readout images improved the CNR by approximately 50% between the gray matter (GM) and white matter (WM) compared to single-echo images. MT-prep images showed excellent contrast between the GM and WM with magnetization transfer ratio (MTR) and cerebrospinal fluid normalized MT signal (MTCSF) = 0.12 ± 0.017 and 0.74 ± 0.013, respectively, for the GM; and 0.18 ± 0.011 and 0.58 ± 0.009, respectively, for the WM. Under-sampled acquisition using rosette readout with CS reconstruction demonstrated up to 6 times faster scans with comparable image quality as the fully-sampled acquisition.

Diffusion-weighted MRI of the spinal cord in cervical spondylotic myelopathy after instrumented fusion

Introduction

This study investigated tissue diffusion properties within the spinal cord of individuals treated for cervical spondylotic myelopathy (CSM) using post-decompression stabilization hardware. While previous research has indicated the potential of diffusion-weighted MRI (DW-MRI) markers of CSM, the metallic implants often used to stabilize the decompressed spine hamper conventional DW-MRI.

Methods

Utilizing recent developments in DW-MRI metal-artifact suppression technologies, imaging data was acquired from 38 CSM study participants who had undergone instrumented fusion, as well as asymptomatic (non-instrumented) control participants. Apparent diffusion coefficients were determined in axial slice sections and split into four categories: a) instrumented levels, b) non-instrumented CSM levels, c) adjacent-segment (to instrumentation) CSM levels, and d) non-instrumented control levels. Multi-linear regression models accounting for age, sex, and body mass index were used to investigate ADC measures within each category. Furthermore, the cord diffusivity within CSM subjects was correlated with symptom scores and the duration since fusion procedures.

Results

ADC measures of the spinal cord in CSM subjects were globally reduced relative to control subjects (p = 0.005). In addition, instrumented levels within the CSM subjects showed reduced diffusivity relative to controls (p = 0.003), while ADC within non-instrumented CSM levels did not statistically deviate from control levels (p = 0.107).

Discussion

Multi-spectral DW-MRI technology can be effectively employed to evaluate cord diffusivity near fusion hardware in subjects who have undergone surgery for CSM. Leveraging this advanced technology, this study had identified significant reductions in cord diffusivity, relative to control subjects, in CSM patients treated with conventional metallic fusion instrumentation.

Decoding cerebro-spinal signatures of human behavior: application to motor sequence learning

Mapping the neural patterns that drive human behavior is a key challenge in neuroscience. Even the simplest of our everyday actions stem from the dynamic and complex interplay of multiple neural structures across the central nervous system (CNS). Yet, most neuroimaging research has focused on investigating cerebral mechanisms, while the way the spinal cord accompanies the brain in shaping human behavior has been largely overlooked. Although the recent advent of functional magnetic resonance imaging (fMRI) sequences that can simultaneously target the brain and spinal cord has opened up new avenues for studying these mechanisms at multiple levels of the CNS, research to date has been limited to inferential univariate techniques that cannot fully unveil the intricacies of the underlying neural states. To address this, we propose to go beyond traditional analyses and instead use a data-driven multivariate approach leveraging the dynamic content of cerebro-spinal signals using innovation-driven coactivation patterns (iCAPs). We demonstrate the relevance of this approach in a simultaneous brain-spinal cord fMRI dataset acquired during motor sequence learning (MSL), to highlight how large-scale CNS plasticity underpins rapid improvements in early skill acquisition and slower consolidation after extended practice. Specifically, we uncovered cortical, subcortical and spinal functional networks, which were used to decode the different stages of learning with a high accuracy and, thus, delineate meaningful cerebro-spinal signatures of learning progression. Our results provide compelling evidence that the dynamics of neural signals, paired with a data-driven approach, can be used to disentangle the modular organization of the CNS. While we outline the potential of this framework to probe the neural correlates of motor learning, its versatility makes it broadly applicable to explore the functioning of cerebro-spinal networks in other experimental or pathological conditions.

A Comparison of Functional Connectivity in the Human Brainstem and Spinal Cord Associated with Noxious and Innocuous Thermal Stimulation Identified by Means of Functional MRI

The somatosensory system is multidimensional and processes important information for survival, including the experience of pain. The brainstem and spinal cord serve pivotal roles in both transmitting and modulating pain signals from the periphery; although, they are studied less frequently with neuroimaging when compared to the brain. In addition, imaging studies of pain often lack a sensory control condition, failing to differentiate the neural processes associated with pain versus innocuous sensations. The purpose of this study was to investigate neural connectivity between key regions involved in descending modulation of pain in response to a hot, noxious stimulus as compared to a warm, innocuous stimulus. This was achieved with functional magnetic resonance imaging (fMRI) of the brainstem and spinal cord in 20 healthy men and women. Functional connectivity was observed to vary between specific regions across painful and innocuous conditions. However, the same variations were not observed in the period of anticipation prior to the onset of stimulation. Specific connections varied with individual pain scores only during the noxious stimulation condition, indicating a significant role of individual differences in the experience of pain which are distinct from that of innocuous sensation. The results also illustrate significant differences in descending modulation before and during stimulation in both conditions. These findings contribute to a deeper understanding of the mechanisms underlying pain processing at the level of the brainstem and spinal cord, and how pain is modulated.

Utility of MRI in Quantifying Tissue Injury in Cervical Spondylotic Myelopathy

Cervical spondylotic myelopathy (CSM) is a progressive disease that worsens over time if untreated. However, the rate of progression can vary among individuals and may be influenced by various factors, such as the age of the patients, underlying conditions, and the severity and location of the spinal cord compression. Early diagnosis and prompt treatment can help slow the progression of CSM and improve symptoms. There has been an increased use of magnetic resonance imaging (MRI) methods in diagnosing and managing CSM. MRI methods provide detailed images and quantitative structural and functional data of the cervical spinal cord and brain, allowing for an accurate evaluation of the extent and location of tissue injury. This review aims to provide an understanding of the use of MRI methods in interrogating functional and structural changes in the central nervous system in CSM. Further, we identified several challenges hindering the clinical utility of these neuroimaging methods.

Reliability of resting-state functional connectivity in the human spinal cord: assessing the impact of distinct noise sources

The investigation of spontaneous fluctuations of the blood-oxygen-level-dependent (BOLD) signal has recently been extended from the brain to the spinal cord, where it has stimulated interest from a clinical perspective. A number of resting-state functional magnetic resonance imaging (fMRI) studies have demonstrated robust functional connectivity between the time series of BOLD fluctuations in bilateral dorsal horns and between those in bilateral ventral horns, in line with the functional neuroanatomy of the spinal cord. A necessary step prior to extension to clinical studies is assessing the reliability of such resting-state signals, which we aimed to do here in a group of 45 healthy young adults at the clinically prevalent field strength of 3T. When investigating connectivity in the entire cervical spinal cord, we observed fair to good reliability for dorsal-dorsal and ventral-ventral connectivity, whereas reliability was poor for within- and between-hemicord dorsal-ventral connectivity. Considering how prone spinal cord fMRI is to noise, we extensively investigated the impact of distinct noise sources and made two crucial observations: removal of physiological noise led to a reduction in functional connectivity strength and reliability - due to the removal of stable and participant-specific noise patterns - whereas removal of thermal noise considerably increased the detectability of functional connectivity without a clear influence on reliability. Finally, we also assessed connectivity within spinal cord segments and observed that while the pattern of connectivity was similar to that of whole cervical cord, reliability at the level of single segments was consistently poor. Taken together, our results demonstrate the presence of reliable resting-state functional connectivity in the human spinal cord even after thoroughly accounting for physiological and thermal noise, but at the same time urge caution if focal changes in connectivity (e.g. due to segmental lesions) are to be studied, especially in a longitudinal manner.

Reduced Spinal Cord Gray Matter in Patients with Fibromyalgia Using Opioids Long-term

Objective

Chronic pain involves alterations in brain gray matter volume (GMV). Moreover, opioid medications are known to reduce GMV in numerous brain regions involved in pain processing. However, no research has evaluated (1) chronic pain-related GMV alterations in the spinal cord or (2) the effect of opioids on spinal cord GMV. Accordingly, this study evaluated spinal cord GMV in health controls and patients with fibromyalgia who were using and not using opioids long-term.

Methods

We analyzed average C5 - C7 GMV of the spinal cord dorsal and ventral horns in separate female cohorts of healthy controls (HC, n = 30), fibromyalgia patients not using opioids (FMN, n = 31), and fibromyalgia patients using opioids long-term (FMO, n = 27). To assess the effect of group on average dorsal and ventral horn GMV, we conducted a one-way multivariate analysis of covariance.

Results

After controlling for age, we observed a significant effect of group on ventral horn GMV (p = 0.03, η2 = 0.09), and on dorsal horn GMV (p = 0.05, η2 = 0.08). Tukey's posthoc comparisons showed that, compared to HC participants, FMOs had significantly lower ventral (p = 0.01) and dorsal (p = 0.02) GMVs. Among FMOs only, ventral horn GMV was significantly positively associated with pain severity and interference, and both dorsal and ventral GMVs were significantly positively associated with cold pain tolerance.

Conclusion

Long-term opioid use may impact sensory processing in fibromyalgia via gray matter changes within the cervical spinal cord.

Assessing the Equivalence of Brain-Derived Measures from Two 3D T1-Weighted Acquisitions: One Covering the Brain and One Covering the Brain and Spinal Cord

BACKGROUND AND PURPOSE

In MS, it is common to acquire brain and spinal cord MR imaging sequences separately to assess the extent of the disease. The goal of this study was to see how replacing the traditional brain T1-weighted images (brain-T1) with an acquisition that included both the brain and the cervical spinal cord (cns-T1) affected brain- and spinal cord-derived measures.

MATERIALS AND METHODS

Thirty-six healthy controls (HC) and 42 patients with MS were included. Of those, 18 HC and 35 patients with MS had baseline and follow-up at 1 year acquired on a 3T magnet. Two 3D T1-weighted images (brain-T1 and cns-T1) were acquired at each time point. Regional cortical thickness and volumes were determined with FastSurfer, and the percentage brain volume change per year was obtained with SIENA. The spinal cord area was estimated with the Spinal Cord Toolbox. Intraclass correlation coefficients (ICC) were calculated to check for consistency of measures obtained from brain-T1 and cns-T1.

RESULTS

Cortical thickness measures showed an ICC >0.75 in 94% of regions in healthy controls and 80% in patients with MS. Estimated regional volumes had an ICC >0.88, and the percentage brain volume change had an ICC >0.79 for both groups. The spinal cord area measures had an ICC of 0.68 in healthy controls and 0.92 in patients with MS.

CONCLUSIONS

Brain measurements obtained from 3D cns-T1 are highly equivalent to those obtained from a brain-T1, suggesting that it could be feasible to replace the brain-T1 with cns-T1.

Differential Trajectory of Diffusion and Perfusion Magnetic Resonance Imaging of Rat Spinal Cord Injury

Traumatic spinal cord injury causes rapid neuronal and vascular injury, and predictive biomarkers are needed to facilitate acute patient management. This study examined the progression of magnetic resonance imaging (MRI) biomarkers after spinal cord injury and their ability to predict long-term neurological outcomes in a rodent model, with an emphasis on diffusion-weighted imaging (DWI) markers of axonal injury and perfusion-weighted imaging of spinal cord blood flow (SCBF). Adult Sprague-Dawley rats received a cervical contusion injury of varying severity (injured = 30, sham = 9). MRI at 4 h, 48-h, and 12-weeks post-injury included T1, T2, perfusion, and DWI. Locomotor outcome was assessed up to 12 weeks post-injury. At 4 h, the deficit in SCBF was larger than the DWI lesion, and although SCBF partially recovered by 48 h, the DWI lesion expanded. At 4 h, the volume of the SCBF deficit (R2 = 0.56, padj < 0.01) was significantly correlated with 12-week locomotor outcome, whereas DWI (R2 = 0.30, padj < 0.01) was less predictive of outcome. At 48 h, SCBF (R2 = 0.41, padj < 0.01) became less associated with outcome, and DWI (R2 = 0.38, padj < 0.01) lesion volume became more closely related to outcome. Spinal cord perfusion has unique spatiotemporal dynamics compared with diffusion measures of axonal damage and highlights the importance of acute perfusion abnormalities. Perfusion and diffusion offer complementary and clinically relevant insight into physiological and structural abnormalities following spinal cord injury beyond those afforded by T1 or T2 contrasts.

Could spinal cord oscillation contribute to spinal cord injury in degenerative cervical myelopathy?

Introduction

Degenerative Cervical Myelopathy [DCM] is a slow-motion spinal cord injury. Compression and dynamic compression have been considered disease hallmarks. However, this is likely an oversimplification, as compression is more commonly incidental and has only modest correlation to disease severity. MRI studies have recently suggested spinal cord oscillation could play a role.

Research question

To determine if spinal cord oscillation could contribute to spinal cord injury in degenerative cervical myelopathy.

Material and methods

A computational model of an oscillating spinal cord was developed from imaging of a healthy volunteer. Using finite element analysis, the observed implications of stress and strain, were measured in the context of a simulated disc herniation. The significance was bench marked by comparison to a more recognised dynamic injury mechanism; a flexion extension model of dynamic compression.

Results

Spinal cord oscillation altered both compressive and shear strain on the spinal cord. Following initial compression, compressive strain moves from within the spinal cord to the spinal cord surface, whilst shear strain is magnified by 0.1-0.2, depending on the amplitude of oscillation. These orders of magnitude are equivalent to a dynamic compression model.

Discussion and conclusion

Spinal cord oscillation could significantly contribute to spinal cord damage across DCM. Its repeated occurrence with every heartbeat, draws parallels to the concept of fatigue damage, which could reconcile differing theories on the origins of DCM. This remains hypothetical at this stage, and further investigations are required.

Shimming toolbox: An open-source software toolbox for B0 and B1 shimming in MRI

PURPOSE

Introduce Shimming Toolbox ( https://shimming-toolbox.org), an open-source software package for prototyping new methods and performing static, dynamic, and real-time B0 shimming as well as B1 shimming experiments.

METHODS

Shimming Toolbox features various field mapping techniques, manual and automatic masking for the brain and spinal cord, B0 and B1 shimming capabilities accessible through a user-friendly graphical user interface. Validation of Shimming Toolbox was demonstrated in three scenarios: (i) B0 dynamic shimming in the brain at 7T using custom AC/DC coils, (ii) B0 real-time shimming in the spinal cord at 3T, and (iii) B1 static shimming in the spinal cord at 7T.

RESULTS

The B0 dynamic shimming of the brain at 7T took about 10 min to perform. It showed a 47% reduction in the standard deviation of the B0 field, associated with noticeable improvements in geometric distortions in EPI images. Real-time dynamic xyz-shimming in the spinal cord took about 5 min and showed a 30% reduction in the standard deviation of the signal distribution. B1 static shimming experiments in the spinal cord took about 10 min to perform and showed a 40% reduction in the coefficient of variation of the B1 field.

CONCLUSION

Shimming Toolbox provides an open-source platform where researchers can collaborate, prototype and conveniently test B0 and B1 shimming experiments. Future versions will include additional field map preprocessing techniques, optimization algorithms, and compatibility across multiple MRI manufacturers.

Reliability of spinal cord measures based on synthetic T1-weighted MRI derived from multiparametric mapping (MPM)

Short MRI acquisition time, high signal-to-noise ratio, and high reliability are crucial for image quality when scanning healthy volunteers and patients. Cross-sectional cervical cord area (CSA) has been suggested as a marker of neurodegeneration and potential outcome measure in clinical trials and is conventionally measured on T₁-weigthed 3D Magnetization Prepared Rapid Acquisition Gradient-Echo (MPRAGE) images. This study aims to reduce the acquisition time for the comprehensive assessment of the spinal cord, which is typically based on MPRAGE for morphometry and multi-parameter mapping (MPM) for microstructure. The MPRAGE is replaced by a synthetic T₁-w MRI (synT₁-w) estimated from the MPM, in order to measure CSA. SynT₁-w images were reconstructed using the MPRAGE signal equation based on quantitative maps of proton density (PD), longitudinal (R₁) and effective transverse (R₂*) relaxation rates. The reliability of CSA measurements from synT₁-w images was determined within a multi-center test-retest study format and validated against acquired MPRAGE scans by assessing the agreement between both methods. The response to pathological changes was tested by longitudinally measuring spinal cord atrophy following spinal cord injury (SCI) for synT₁-w and MPRAGE using linear mixed effect models. CSA measurements based on the synT₁-w MRI showed high intra-site (Coefficient of variation [CoV]: 1.43% to 2.71%) and inter-site repeatability (CoV: 2.90% to 5.76%), and only a minor deviation of -1.65 mm² compared to MPRAGE. Crucially, by assessing atrophy rates and by comparing SCI patients with healthy controls longitudinally, differences between synT₁-w and MPRAGE were negligible. These results demonstrate that reliable estimates of CSA can be obtained from synT₁-w images, thereby reducing scan time significantly.

Multiple sclerosis iron rim lesions are linked to impaired cervical spinal cord integrity using the T1/T2-weighted ratio

BACKGROUND AND PURPOSE

In multiple sclerosis (MS), iron rim lesions (IRLs) are characterized by pronounced tissue matrix damage. The T1/T2-weighted (T1/T2w) ratio represents a postprocessing MRI approach to investigate tissue integrity, but studies investigating spinal cord pathology are missing until now. The aim of this study was to characterize tissue integrity using the T1/T2w ratio in lesions and the normal-appearing white and gray matter (NAWM, NAGM) in the spinal cord and brain in MS patients with and without brain IRLs.

METHODS

Forty MS patients (20 patients with at least one brain IRL and 20 age- and sex-matched patients without IRLs) were included. Normalized cross-sectional area (nCSA) of the upper cervical cord was calculated in addition to T1/T2w values and standard brain and spinal cord MRI parameters.

RESULTS

Patients with IRLs had higher disability scores, a smaller nCSA, and a higher cervical T2 lesion volume. T1/T2w values of brain IRLs were significantly lower compared to non-IRLs (p < .001). Furthermore, T1/T2w values of lesions were significantly lower compared to the NAGM and NAWM, both in the brain and the spinal cord (p < .05 for all comparisons). T1/T2w values of the NAGM and NAWM in the brain and spinal cord did not statistically differ between the IRL group and the non-IRL group.

CONCLUSION

IRLs constitute an imaging marker of disease severity. T1/T2w ratio maps represent an interesting technique to capture diffuse tissue properties. Calculation of T1/T2w ratio maps of the spinal cord might provide additional insights into the pathophysiological processes of MS.

Improving T2*-weighted human cortico-spinal acquisitions with a dedicated algorithm for region-wise shimming

Cortico-spinal fMRI acquisitions aim to investigate direct interactions between brain and spinal cord, e.g. during motor output or pain processing, by covering both regions in a single measurement. Due to their large distance and location in the body, a dynamic shim update of constant and linear shim terms is required when using echo-planar imaging (EPI) to achieve reasonable image quality in both target regions. A previously presented approach with region-wise shim settings is based on a standard single-region shim algorithm and suffers from (i) non-optimal shim settings because it combines linear and second-order shim terms optimized for different volumes, and (ii) significant user interactions making it rather cumbersome, time consuming, and error-prone. Here, a dedicated ("CoSpi") shim algorithm for cortico-spinal fMRI is presented that performs joint optimization of static second-order shim terms and one set of linear and constant shim terms for each region in a single run and with minimal user interaction. Field map and T2*-weighted EPI measurements were performed on a clinical 3 T whole-body MR system in water phantoms and five healthy volunteers using the conventional region-wise and CoSpi shim settings as well as "gold standard" shim settings optimized for one of the target regions only. With CoSpi shim settings, (i) overall field inhomogeneity was reduced by about 65% / 75% (brain / spinal cord volume) compared to the conventional region-wise approach and in vivo was within 5% of the values obtained with the single-volume shim settings, (ii) geometric distortions derived from voxel displacement maps were reduced on average by about 35% / 70%, (iii) the temporal SNR determined from an EPI time series that may reflect the impact of through-slice dephasing, was increased by about 17% / 10%, and (iv) the variation of the mean field between slices, a measure targeting the predisposition to insufficient fat saturation and GRAPPA-related ghosting artifacts, was reduced by about 90% / 45%. Thus, the presented algorithm not only speeds up and simplifies the shim procedure considerably, but also provides a better field homogeneity and image quality, which both could help to significantly improve the applicability of cortico-spinal fMRI.

Deep Learning-Based Auto-Segmentation of Spinal Cord Internal Structure of Diffusion Tensor Imaging in Cervical Spondylotic Myelopathy

Cervical spondylotic myelopathy (CSM) is a chronic disorder of the spinal cord. ROI-based features on diffusion tensor imaging (DTI) provide additional information about spinal cord status, which would benefit the diagnosis and prognosis of CSM. However, the manual extraction of the DTI-related features on multiple ROIs is time-consuming and laborious. In total, 1159 slices at cervical levels from 89 CSM patients were analyzed, and corresponding fractional anisotropy (FA) maps were calculated. Eight ROIs were drawn, covering both sides of lateral, dorsal, ventral, and gray matter. The UNet model was trained with the proposed heatmap distance loss for auto-segmentation. Mean Dice coefficients on the test dataset for dorsal, lateral, and ventral column and gray matter were 0.69, 0.67, 0.57, 0.54 on the left side and 0.68, 0.67, 0.59, 0.55 on the right side. The ROI-based mean FA value based on segmentation model strongly correlated with the value based on manual drawing. The percentages of the mean absolute error between the two values of multiple ROIs were 0.07, 0.07, 0.11, and 0.08 on the left side and 0.07, 0.1, 0.1, 0.11, and 0.07 on the right side. The proposed segmentation model has the potential to offer a more detailed spinal cord segmentation and would be beneficial for quantifying a more detailed status of the cervical spinal cord.

Altered Spinal Cord Functional Connectivity Associated with Parkinson's Disease Progression

BACKGROUND

Parkinson's disease (PD) has traditionally been viewed as an α-synucleinopathy brain pathology. Yet evidence based on postmortem human and animal experimental models indicates that the spinal cord may also be affected.

OBJECTIVE

Functional magnetic resonance imaging (fMRI) seems to be a promising candidate to better characterize spinal cord functional organization in PD patients.

METHODS

Resting-state spinal fMRI was performed in 70 PD patients and 24 age-matched healthy controls, the patients being divided into three groups based on their motor symptom severity: PD_low (n = 24), PD_med (n = 22), and PD_adv (n = 24) groups. A combination of independent component analysis (ICA) and a seed-based approach was applied.

RESULTS

When pooling all participants, the ICA revealed distinct ventral and dorsal components distributed along the rostro-caudal axis. This organization was highly reproducible within subgroups of patients and controls. PD severity, assessed by Unified Parkinson's Disease Rating Scale (UPDRS) scores, was associated with a decrease in spinal functional connectivity (FC). Notably, we observed a reduced intersegmental correlation in PD as compared to controls, the latter being negatively associated with patients' upper-limb UPDRS scores (P = 0.0085). This negative association between FC and upper-limb UPDRS scores was significant between adjacent C4-C5 (P = 0.015) and C5-C6 (P = 0.20) cervical segments, levels associated with upper-limb functions.

CONCLUSIONS

The present study provides the first evidence of spinal cord FC changes in PD and opens new avenues for the effective diagnosis and therapeutic strategies in PD. This underscores how spinal cord fMRI can serve as a powerful tool to characterize, in vivo, spinal circuits for a variety of neurological diseases. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Association of volumetric MRI measures and disability in MS patients of the same age: Descriptions from a birth year cohort

BACKGROUND AND OBJECTIVES

Although MRI-based markers of neuroinflammation have proven crucial for the diagnosis of multiple sclerosis (MS), predicting clinical progression with inflammation remains difficult. Neurodegenerative markers such as brain volume loss show stronger clinical (predictive) correlations, but also harbor age-related variation that must be disentangled from disease duration. In this study we investigated how clinical disability is related to volumetric MRI measures in a cohort of MS patients and healthy controls (HC) of the same age: Project Y.

METHODS

This study included 234 MS patients born in 1966 and 112 HC born between 1965 and 1967 in the Netherlands. Disability was quantified using the expanded disability status scale (EDSS), nine hole peg test (9HPT), and timed 25 foot walking test (T25FWT). Volumes were quantified on 3T MRI as normalized whole brain (NBV) and regional gray matter (GM) volumes using the same scanner and MRI protocol: cortical (normalized cortical gray matter volume; NCGMV), deep (NDGMV), thalamic (NThalV), and cerebellar (NCbV) GM volumes. In addition, mean upper cervical cord area (MUCCA), white matter lesion volume (LV), and spinal cord lesions were assessed. These measures were compared between patients and HC, and related to disability measures using linear regression.

RESULTS

Mean age of people with MS (PwMS) was 52.8 years (SD 0.9) and median disease duration 15.8 years (IQR 8.7-24.8). All global and regional brain measures were lower in MS patients compared to HC. Univariate regression models showed that NDGMV (β = -0.20) and MUCCA (β = -0.38) were most strongly related to the EDSS in all PwMS. After subtype stratification, MUCCA was most strongly related to the EDSS (β = -0.60) and 9HPT (β = -0.55) in secondary progressive PwMS. Multivariate regression models demonstrated that in all PwMS, the EDSS was best explained by lower MUCCA, longer disease durations and a progressive disease course (adjusted-R (Sastre-Garriga et al., 2017) = 0.26, p < 0.001). MUCCA was a consistent correlate in separate models of the EDSS for all PwMS, relapsing and progressive onset PwMS. The 9HPT (adjusted-R (Sastre-Garriga et al., 2017) = 0.20, p < 0.001) was best explained by lower MUCCA, higher LV and pack years, while lower limb disability (adjusted-R (Sastre-Garriga et al., 2017) = 0.11, p < 0.001) was best explained by lower MUCCA, progressive onset MS and female sex.

DISCUSSION

Our results indicate that in a cohort unbiased by age differences, spinal cord and deep gray matter volumes best related to physical disability. Our results support the use of these measures in clinical practice and trials.

8-channel Tx dipole and 20-channel Rx loop coil array for MRI of the cervical spinal cord at 7 Tesla

The quality of cervical spinal cord images can be improved by the use of tailored radiofrequency coil solutions for ultra-high field imaging; however, very few commercial and research 7 Tesla radiofrequency coils currently exist for the spinal cord, and in particular those with parallel transmit capabilities. This work presents the design, testing and validation of a pTx/Rx coil for the human neck and cervical/upper-thoracic spinal cord. The pTx portion is composed of 8 dipoles to ensure high homogeneity over this large region of the spinal cord. The Rx portion is made of 20 semi-adaptable overlapping loops to produce high Signal-to-noise ratio (SNR) across the patient population. The coil housing is designed to facilitate patient positioning and comfort, while being tight fitting to ensure high sensitivity. We demonstrate RF shimming capabilities to optimize B ₁ ⁺ uniformity, power efficiency and/or specific absorption rate (SAR) efficiency. B ₁ ⁺ homogeneity, SNR and g-factor was evaluated in adult volunteers and demonstrated excellent performance from the occipital lobe down to the T4-T5 level. We compared the proposed coil with two state-of-the-art head and head/neck coils, confirming its superiority in the cervical and upper-thoracic regions of the spinal cord. This coil solution therefore provides a convincing platform for producing the high image quality necessary for clinical and research scanning of the upper spinal cord.

Remodeling of the brain correlates with gait instability in cervical spondylotic myelopathy

Introduction

Cervical spondylotic myelopathy (CSM) is a common form of non-traumatic spinal cord injury (SCI) and usually leads to remodeling of the brain and spinal cord. In CSM with gait instability, the remodeling of the brain and cervical spinal cord is unclear. We attempted to explore the remodeling of these patients' brains and spinal cords, as well as the relationship between the remodeling of the brain and spinal cord and gait instability.

Methods

According to the CSM patients' gait, we divided patients into two groups: normal gait patients (nPT) and abnormal gait patients (aPT). Voxel-wise z-score transformation amplitude of low-frequency fluctuations (zALFF) and resting-state functional connectivity (rs-FC) were performed for estimating brain changes. Cross-sectional area (CSA) and fractional anisotropy (FA) of the spinal cord were computed by Spinal cord toolbox. Correlations of these measures and the modified Japanese Orthopedic Association (mJOA) score were analyzed.

Results

We found that the zALFF of caudate nucleus in aPT was higher than that in healthy controls (HC) and lower than that in nPT. The zALFF of the right postcentral gyrus and paracentral lobule in HC was higher than those of aPT and nPT. Compared with the nPT, the aPT showed increased functional connectivity between the caudate nucleus and left angular gyrus, bilateral precuneus and bilateral posterior cingulate cortex (PCC), which constitute a vital section of the default mode network (DMN). No significantly different FA values or CSA of spinal tracts at the C2 level were observed between the HC, nPT and aPT groups. In CSM, the right paracentral lobule's zALFF was negatively correlated with the FA value of fasciculus gracilis (FCG), and the right caudate zALFF was positively correlated with the FA value of the fasciculus cuneatus (FCC). The results showed that the functional connectivity between the right caudate nucleus and DMN was negatively correlated with the CSA of the lateral corticospinal tract (CST).

Discussion

The activation of the caudate nucleus and the strengthening functional connectivity between the caudate nucleus and DMN were associated with gait instability in CSM patients. Correlations between spinal cord and brain function might be related to the clinical symptoms in CSM.