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

Calibration-free parallel transmission of the cervical, thoracic, and lumbar spinal cord at 7T

PURPOSE

To address the limitations of spinal cord imaging at ultra-high field (UHF) due to time-consuming parallel transmit (pTx) adjustments. This study introduces calibration-free offline computed universal shim modes that can be applied seamlessly for different pTx RF coils and spinal cord target regions, substantially enhancing spinal cord imaging efficiency at UHF.

METHODS

A library of channel-wise relativeB 1 + $$ {B}_1^{+} $$ maps for the cervical spinal cord (six datasets) and thoracic and lumbar spinal cord (nine datasets) was constructed to optimize transmit homogeneity and efficiency for these regions. A tailored B0 shim was optimized for the cervical spine to enhance spatial magnetic field homogeneity further. The performance of the universal shims was validated using absolute saturation basedB 1 + $$ {B}_1^{+} $$ mapping and high-resolution 2D and 3D multi-echo gradient-recalled echo (GRE) data to assess the image quality.

RESULTS

The proposed universal shims demonstrated a 50% improvement inB 1 + $$ {B}_1^{+} $$ efficiency compared to the default (zero phase) shim mode.B 1 + $$ {B}_1^{+} $$ homogeneity was also improved by 20%. The optimized universal shims achieved performance comparable to subject-specific pTx adjustments, while eliminating the need for lengthy pTx calibration times, saving about 10 min per experiment.

CONCLUSION

The development of universal shims represents a significant advance by eliminating time-consuming subject-specific pTx adjustments. This approach is expected to make UHF spinal cord imaging more accessible and user-friendly, particularly for non-pTx experts.

Evaluation of cervical spinal cord atrophy using a modified SIENA approach

Spinal cord (SC) atrophy obtained from structural magnetic resonance imaging has gained relevance as an indicator of neurodegeneration in various neurological disorders. The common method to assess SC atrophy is by comparing numerical differences of the cross-sectional spinal cord area (CSA) between time points. However, this indirect approach leads to considerable variability in the obtained results. Studies showed that this limitation can be overcome by using a registration-based technique. The present study introduces the Structural Image Evaluation using Normalization of Atrophy on the Spinal Cord (SIENA-SC), which is an adapted version of the original SIENA method, designed to directly calculate the percentage of SC volume change over time from clinical brain MRI acquired with an extended field of view to cover the superior part of the cervical SC. In this work, we compared SIENA-SC with the Generalized Boundary Shift Integral (GBSI) and the CSA change. On a scan-rescan dataset, SIENA-SC was shown to have the lowest measurement error than the other two methods. When comparing a group of 190 Healthy Controls with a group of 65 Multiple Sclerosis patients, SIENA-SC provided significantly higher yearly rates of atrophy in patients than in controls and a lower sample size when measured for treatment effect sizes of 50%, 30% and 10%. Our findings indicate that SIENA-SC is a robust, reproducible, and sensitive approach for assessing longitudinal changes in spinal cord volume, providing neuroscientists with an accessible and automated tool able to reduce the need for manual intervention and minimize variability in measurements.

Sensitivity of Advanced Magnetic Resonance Imaging to Progression over Six Months in Early Spinocerebellar Ataxia

BACKGROUND

Clinical trials for upcoming disease-modifying therapies of spinocerebellar ataxias (SCA), a group of rare movement disorders, lack endpoints sensitive to early disease progression, when therapeutics will be most effective. In addition, regulatory agencies emphasize the importance of biological outcomes.

OBJECTIVES

READISCA, a transatlantic clinical trial readiness consortium, investigated whether advanced multimodal magnetic resonance imaging (MRI) detects pathology progression over 6 months in preataxic and early ataxic carriers of SCA mutations.

METHODS

A total of 44 participants (10 SCA1, 25 SCA3, and 9 controls) prospectively underwent 3-T MR scanning at baseline and a median [interquartile range] follow-up of 6.2 [5.9-6.7] months; 44% of SCA participants were preataxic. Blinded analyses of annual changes in structural, diffusion MRI, MR spectroscopy, and the Scale for Assessment and Rating of Ataxia (SARA) were compared between groups using nonparametric testing. Sample sizes were estimated for 6-month interventional trials with 50% to 100% treatment effect size, leveraging existing large cohort data (186 SCA1, 272 SCA3) for the SARA estimate.

RESULTS

Rate of change in microstructural integrity (decrease in fractional anisotropy, increase in diffusivities) in the middle cerebellar peduncle, corona radiata, and superior longitudinal fasciculus significantly differed in SCAs from controls (P < 0.005), with high effect sizes (Cohen's d = 1-2) and moderate-to-high responsiveness (|standardized response mean| = 0.6-0.9) in SCAs. SARA scores did not change, and their rate of change did not differ between groups.

CONCLUSIONS

Diffusion MRI is sensitive to disease progression at very early-stage SCA1 and SCA3 and may provide a >5-fold reduction in sample sizes relative to SARA as endpoint for 6-month-long trials. © 2024 The Author(s). Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Pre-operative spinal cord perfusion quantified by DSC MRI as a predictor of post-operative prognosis in patients with cervical spondylotic myelopathy

OBJECTIVE

This study aims to investigate the potential of preoperative blood supply condition measured by dynamic susceptibility contract (DSC) MRI in prediction of postoperative outcomes for patients with cervical spondylotic myelopathy (CSM).

MATERIALS AND METHOD

Thirty-nine patients (Age: 61 ± 7, male: 23, female: 16) with CSM who underwent laminoplasty were enrolled. All patients received DSC MRI before the operation. Five parameters include Enhance, rEnhance, full width at half maxima (FWHM), Slope1 and Slope2 in DSC MRI, were calculated at all the compressed spinal cord segments. Clinical outcomes were evaluated by modified Japanese Orthopaedic Association (mJOA) scores. Patients were divided into two groups based on mJOA recovery rate of 5 years: good recovery (> 50%) or poor recovery (≤ 50%). The difference between two groups were compared. The value of DSC MRI to CSM was evaluated by logistic and receiver operating characteristic (ROC) curve analysis.

RESULTS

There were 26 patients in good recovery group and 13 patients in poor recovery group. The baseline characteristics, including age, gender, preoperative mJOA score, and smoking status showed no significant difference between the two groups (all p > 0.05). The FWHM was significantly higher in the poor recovery group (9.77 ± 2.78) compared to the good recovery group (6.64 ± 1.65) (p = 0.002). Logistic regression analysis indicated that an increased FWHM was a significant risk factor for poor prognosis recovery (p = 0.013, OR = 0.392, 95%CI: 0.187-0.822). The AUC of FWHM for ROC was 0.843 (95% CI: 0.710-0.975) with a p value of 0.001. In addition, an FWHM greater than 5.87, with a sensitivity of 92.3% and specificity of 69.2%, was found to be an independent risk factor for poor postoperative recovery in patients with CSM.

CONCLUSION

In this study, we successfully quantified the spinal cord blood supply condition by DSC MRI technique. We found that an increase in FWHM was an independent risk factor for poor postoperative recovery in CSM patients. Specifically, patients with FWHM > 5.87 have a poor postoperative recovery.

Cervical spinal cord morphometrics in degenerative cervical myelopathy: quantification using semi-automated normalized technique and correlation with neurological dysfunctions

BACKGROUND CONTEXT

Degenerative cervical myelopathy (DCM) is characterized by spinal cord atrophy. Accurate estimation of spinal cord atrophy is key to the understanding of neurological diseases, including DCM. However, its clinical application is hampered by difficulties in its precise and consistent estimation due to significant variability in spinal cord morphometry along the cervical spine, both within and between individuals.

PURPOSE

To characterize morphometrics of the compressed spinal cord in DCM patients. We employed our semiautomated analysis framework that incorporates the Spinal Cord Toolbox (SCT) and a normalization approach to effectively address the challenges posed by cord compression in these patients. Additionally, we examined the clinical relevance of these morphometric measures to enhance our understanding of DCM pathophysiology.

STUDY DESIGN

Prospective study.

PATIENT SAMPLE

This study investigated 36 DCM patients and 31 healthy controls (HCs).

OUTCOME MEASURES

Clinical scores including 9-hole peg test for hand dexterity, hand grip strength, balance, gait speed, modified Japanese Orthopaedic Association (mJOA) score, and imaging-based spinal cord morphometrics.

METHOD

Using the generic spine acquisition protocol and our semiautomated analysis pipeline, spinal cord morphometrics, including cross-sectional area (CSA), anterior-posterior (AP) and transverse (RL) diameters, eccentricity, and solidity, were estimated from sagittal T2w magnetic resonance imaging (MRI) images using the Spinal Cord Toolbox (SCT). Normalized metrics were extracted from the C1 to C7 vertebral levels and compared between DCM patients and HC. Morphometric data at regions of maximum spinal cord compression (MSCC) were correlated with the clinical scores. A subset of participants underwent follow-up scans at 6 months to monitor longitudinal changes in spinal cord atrophy.

RESULTS

Spinal cord morphometric data were normalized against the healthy population morphometry (PAM50 database) and extracted for all participants. DCM patients showed a notable reduction in CSA, AP, and RL diameter across all vertebral levels compared to HC. MSCC metrics correlated significantly with clinical scores like dexterity, grip strength, and mJOA scores. Longitudinal analysis indicated a decrease in CSA and worsening clinical scores in DCM patients.

CONCLUSION

Our processing pipeline offers a reliable method for assessing spinal cord compression in DCM patients. Normalized spinal cord morphometrics, particularly the CSA could have potential for monitoring DCM disease severity and progression, guiding treatment decisions. Furthermore, to our knowledge our study is the first to apply the generic spinal cord acquisition protocol, ensuring consistent imaging across different MRI scanners and settings. Coupled with our semiautomated analysis pipeline, this protocol is key for the detailed morphometric characterization of compressed spinal cords in patients with DCM, a disease that is both complex and heterogenous. This study was funded by the National Institute of Neurological Disorders and Stroke (NINDS) (K23:NS091430) and (R01: NS129852-01A1).

Advances and Challenges in Spinal Cord Injury Treatments

Spinal cord injury (SCI) is a debilitating condition that is associated with long-term physical and functional disability. Our understanding of the pathogenesis of SCI has evolved significantly over the past three decades. In parallel, significant advances have been made in optimizing the management of patients with SCI. Early surgical decompression, adequate bony decompression and expansile duraplasty are surgical strategies that may improve neurological and functional outcomes in patients with SCI. Furthermore, advances in the non-surgical management of SCI have been made, including optimization of hemodynamic management in the critical care setting. Several promising therapies have also been investigated in pre-clinical studies, with some being translated into clinical trials. Given the recent interest in advancing precision medicine, several investigations have been performed to delineate the role of imaging, cerebral spinal fluid (CSF) and serum biomarkers in predicting outcomes and curating individualized treatment plans for SCI patients. Finally, technological advancements in biomechanics and bioengineering have also found a role in SCI management in the form of neuromodulation and brain-computer interfaces.

Utility of Diffusion Tensor Imaging for Prognosis and Management of Cervical Spondylotic Myelopathy: A PRISMA Review

OBJECTIVE

As advances are made in quantitative magnetic resonance imaging, specifically diffusion tensor imaging, researchers have investigated its potential to serve as a biomarker of disease or prognosticator for postoperative recovery in the management of cervical spondylotic myelopathy. Here, we narratively review the current state of the emerging literature, describing areas of consensus and disagreement.

METHODS

In accordance with the Preferred Reporting Items for Systematic Reviews and Meta-Analyses, we queried 2 large databases for original manuscripts published in English and systematically produced a narrative review of the use of diffusion tensor imaging in the management of cervical spondylotic myelopathy.

RESULTS

Of the 437 manuscripts initially returned in our query, 29 met the final inclusion criteria, and data were extracted regarding diffusion tensor imaging indices and their relationships with clinical outcomes following surgery. Preoperative fractional anisotropy was most commonly found to correlate closely with postsurgical clinical outcomes, though results were mixed.

CONCLUSIONS

Preoperative fractional anisotropy most frequently and best correlates with functional outcomes following surgery for cervical spondylotic myelopathy, according to a review of the current literature. The findings were not universal and at times contradictory, highlighting the need for high-quality future investigations to better define the utility of diffusion tensor imaging in spinal disease.

Spinal cord lesion MRI and behavioral outcomes in a miniature pig model of spinal cord injury: exploring preclinical potential through an ad hoc comparison with human SCI

STUDY DESIGN

prospective case series of Yucatan miniature pig spinal cord contusion injury model with comparison to human cases of spinal cord injury (SCI).

OBJECTIVES

to describe magnetic resonance imaging (MRI) measures of spinal cord lesion severity along with estimates of lateral corticospinal tracts spared neural tissue in both a less severe and more severe contusion SCI model, as well as to describe their corresponding behavioral outcome changes.

SETTING

University laboratory setting.

METHODS

Following a more severe and less severe SCI, each pig underwent spinal cord MRI to measure lesion characteristics, along with locomotor and urodynamics outcomes testing.

RESULTS

In the pig with more severe SCI, locomotor and urodynamic outcomes were poor, and both the spinal cord lesion volume and damage estimates to the lateral corticospinal tracts were large. Conversely, in the pig with less severe SCI, locomotor and urodynamic outcomes were favorable, with the spinal cord lesion volume and damage estimates to the lateral corticospinal tracts being less pronounced. For two human cases matched on estimates of damage to the lateral corticospinal tract regions, the clinical presentations were similar to the pig outcomes, with more limited mobility and more limited bladder functional independence in the more severe case.

CONCLUSIONS

Our initial findings contribute valuable insights to the emergent field of MRI-based evaluation of spinal cord lesions in pig models, offering a promising avenue for understanding and potentially improving outcomes in spinal cord injuries.

Automatic segmentation of dura for quantitative analysis of lumbar stenosis: A deep learning study with 518 CT myelograms

BACKGROUND

The diagnosis of lumbar spinal stenosis (LSS) can be challenging because radicular pain is not often present in the culprit-level localization. Accurate segmentation and quantitative analysis of the lumbar dura on radiographic images are key to the accurate differential diagnosis of LSS. The aim of this study is to develop an automatic dura-contouring tool for radiographic quantification on computed tomography myelogram (CTM) for patients with LSS.

METHODS

A total of 518 CTM cases with or without lumbar stenosis were included in this study. A deep learning (DL) segmentation algorithm 3-dimensional (3D) U-Net was deployed. A total of 210 labeled cases were used to develop the dura-contouring tool, with the ratio of the training, independent testing, and external validation datasets being 150:30:30. The Dice score (DCS) was the primary measure to evaluate the segmentation performance of the 3D U-Net, which was subsequently developed as the dura-contouring tool to segment another unlabeled 308 CTM cases with LSS. Automatic masks of 446 slices on the stenotic levels were then meticulously reviewed and revised by human experts, and the cross-sectional area (CSA) of the dura was compared.

RESULTS

The mean DCS of the 3D U-Net were 0.905 ± 0.080, 0.933 ± 0.018, and 0.928 ± 0.034 in the five-fold cross-validation, the independent testing, and the external validation datasets, respectively. The segmentation performance of the dura-contouring tool was also comparable to that of the second observer (the human expert). With the dura-contouring tool, only 59.0% (263/446) of the automatic masks of the stenotic slices needed to be revised. In the revised cases, there were no significant differences in the dura CSA between automatic masks and corresponding revised masks (p = 0.652). Additionally, a strong correlation of dura CSA was found between the automatic masks and corresponding revised masks (r = 0.805).

CONCLUSIONS

A dura-contouring tool was developed that could automatically segment the dural sac on CTM, and it demonstrated high accuracy and generalization ability. Additionally, the dura-contouring tool has the potential to be applied in patients with LSS because it facilitates the quantification of the dural CSA on stenotic slices.

Biological variation in cervical spinal cord MRI morphometry in healthy individuals and people with multiple sclerosis

BACKGROUND AND PURPOSE

Conclusions from prior literature regarding the impact of sex, age, and height on spinal cord (SC) MRI morphometrics are conflicting, while the effect of body weight on SC morphometrics has been found to be nonsignificant. The purpose of this case-control study is to assess the associations between cervical SC MRI morphometric parameters and age, sex, height, and weight to establish their potential role as confounding variables in a clinical study of people with multiple sclerosis (MS) compared to a cohort of healthy volunteers.

METHODS

Sixty-nine healthy volunteers and 31 people with MS underwent cervical SC MRI at 3 Tesla field strength. Images were centered at the C3/C4 intervertebral disc and processed using Spinal Cord Toolbox v.4.0.2. Mixed-effects linear regression models were used to evaluate the effects of biological variables and disease status on morphometric parameters.

RESULTS

Sex, age, and height had significant effects on cord and gray matter (GM) cross-sectional area (CSA) as well as the GM:cord CSA ratio. There were no significant effects of body weight on morphometric parameters. The effect of MS disease duration on cord CSA in the C4 level was significant when controlling for all other variables.

CONCLUSIONS

Studies of disease-related changes in SC morphometry should control for sex, age, and height to account for physiological variation.

Reproducible Spinal Cord Quantitative MRI Analysis with the Spinal Cord Toolbox

The spinal cord plays a pivotal role in the central nervous system, providing communication between the brain and the body and containing critical motor and sensory networks. Recent advancements in spinal cord MRI data acquisition and image analysis have shown a potential to improve the diagnostics, prognosis, and management of a variety of pathological conditions. In this review, we first discuss the significance of standardized spinal cord MRI acquisition protocol in multi-center and multi-manufacturer studies. Then, we cover open-access spinal cord MRI datasets, which are important for reproducible science and validation of new methods. Finally, we elaborate on the recent advances in spinal cord MRI data analysis techniques implemented in the open-source software package Spinal Cord Toolbox (SCT).

Within-participants reliability and measurement error of magnetization transfer imaging determinations within the healthy cervical spinal cord

PURPOSE

To assess the within-participant reliability and measurement error in the determination of MTR in the healthy human cervical spinal cord.

METHODS AND MATERIALS

A total of twenty healthy controls (10 male, mean ± sd age: 33.9 ± 3.5 years, 10 females, mean ± sd age: 47.5 ± 14.4 years), with no family history of any neurological disorders or a contraindication to MRI scanning were recruited over a period of two months. Each participant was scanned twice with a 3T MRI scanner using standard MTI sequences. Spinal Cord Toolbox (v5.4) was used for image post-processing. Data were first segmented and then registered to a template and then MTR was computed. The within-participant coefficients of variation (CV%), single and average within-participants intraclass correlation coefficients (ICC) and Bland-Altman plots were determined for MT values over the volume between the 2nd and 5th cervical vertebrae for the total WM and for specific WM regions: dorsal column (DC), ventral column (VC) and lateral column (LC).

RESULTS

MTR showed poor to excellent within-participant reliability for the total WM, DC, VC and LC with single/average ICC values of 0.03/0.06, 0.10/0.18, 0.39/0.75, and 0.001/0.002, respectively, and the CV% reported an acceptable variation with values less than 10%. The Bland-Altman plots showed good within-participant agreement between the scan-rescan values.

CONCLUSION

This study demonstrates that clinical trials using MTI technique are feasible and shows that quantitative MTI can monitor tissue changes in degenerative WM patients.

IMPLICATIONS FOR PRACTICE

MTI with its MTR index provide broad assessment of the integrity of white matter tissue and are being studied widely in brain as a diagnostic tool for the assessment of different neurological diseases.

Automated detection, labelling and radiological grading of clinical spinal MRIs

Spinal magnetic resonance (MR) scans are a vital tool for diagnosing the cause of back pain for many diseases and conditions. However, interpreting clinically useful information from these scans can be challenging, time-consuming and hard to reproduce across different radiologists. In this paper, we alleviate these problems by introducing a multi-stage automated pipeline for analysing spinal MR scans. This pipeline first detects and labels vertebral bodies across several commonly used sequences (e.g. T1w, T2w and STIR) and fields of view (e.g. lumbar, cervical, whole spine). Using these detections it then performs automated diagnosis for several spinal disorders, including intervertebral disc degenerative changes in T1w and T2w lumbar scans, and spinal metastases, cord compression and vertebral fractures. To achieve this, we propose a new method of vertebrae detection and labelling, using vector fields to group together detected vertebral landmarks and a language-modelling inspired beam search to determine the corresponding levels of the detections. We also employ a new transformer-based architecture to perform radiological grading which incorporates context from multiple vertebrae and sequences, as a real radiologist would. The performance of each stage of the pipeline is tested in isolation on several clinical datasets, each consisting of 66 to 421 scans. The outputs are compared to manual annotations of expert radiologists, demonstrating accurate vertebrae detection across a range of scan parameters. Similarly, the model's grading predictions for various types of disc degeneration and detection of spinal metastases closely match those of an expert radiologist. To aid future research, our code and trained models are made publicly available.

Quantifying neurodegeneration of the cervical cord and brain in degenerative cervical myelopathy: A multicentre study using quantitative magnetic resonance imaging

BACKGROUND AND PURPOSE

Simultaneous assessment of neurodegeneration in both the cervical cord and brain across multiple centres can enhance the effectiveness of clinical trials. Thus, this study aims to simultaneously assess microstructural changes in the cervical cord and brain above the stenosis in degenerative cervical myelopathy (DCM) using quantitative magnetic resonance imaging (MRI) in a multicentre study.

METHODS

We applied voxelwise analysis with a probabilistic brain/spinal cord template embedded in statistical parametric mappin (SPM-BSC) to process multi parametric mapping (MPM) including effective transverse relaxation rate (R2*), longitudinal relaxation rate (R1), and magnetization transfer (MT), which are indirectly sensitive to iron and myelin content. Regression analysis was conducted to establish associations between neurodegeneration and clinical impairment. Thirty-eight DCM patients (mean age ± SD = 58.45 ± 11.47 years) and 38 healthy controls (mean age ± SD = 41.18 ± 12.75 years) were recruited at University Hospital Balgrist, Switzerland and Toronto Western Hospital, Canada.

RESULTS

Remote atrophy was observed in the cervical cord (p = 0.002) and in the left thalamus (0.026) of the DCM group. R1 was decreased in the periaqueductal grey matter (p = 0.014), thalamus (p = 0.001), corpus callosum (p = 0.0001), and cranial corticospinal tract (p = 0.03). R2* was increased in the primary somatosensory cortices (p = 0.008). Sensory impairments were associated with increased iron-sensitive R2* in the thalamus and periaqueductal grey matter in DCM.

CONCLUSIONS

Simultaneous assessment of the spinal cord and brain revealed DCM-induced demyelination, iron deposition, and atrophy. The extent of remote neurodegeneration was associated with sensory impairment, highlighting the intricate and expansive nature of microstructural neurodegeneration in DCM, reaching beyond the stenosis level.

The relationship between the structural changes in the cervical spinal cord and sensorimotor function of children with thoracolumbar spinal cord injury (TLSCI)

STUDY DESIGN

Cross-sectional study.

OBJECTIVES

To study the relationship between the structural changes in the cervical spinal cord (C2/3 level) and the sensorimotor function of children with traumatic thoracolumbar spinal cord injury (TLSCI) and to discover objective imaging biomarkers to evaluate its functional status.

SETTING

Xuanwu Hospital, Capital Medical University, China; Beijing Key Laboratory of Magnetic Resonance Imaging and Brain Informatics, China.

METHODS

30 children (age range 5-13 years) with TLSCI and 11 typically developing (TD) children (age range 6-12 years) were recruited in this study. Based on whether there is preserved motor function below the neurological level of injury (NLI), the children with TLSCI are divided into the AIS A/B group (motor complete) and the AIS C/D group (motor incomplete). A Siemens Verio 3.0 T MR scanner was used to acquire 3D high-resolution anatomic scans covering the head and upper cervical spinal cord. Morphologic parameters of the spinal cord at the C2/3 level, including cross-sectional area (CSA), anterior-posterior width (APW), and left-right width (LRW) were obtained using the spinal cord toolbox (SCT; https://www.nitrc.org/projects/sct ). Correlation analyses were performed to compare the morphologic spinal cord parameters and clinical scores determined by the International Standard for Neurological Classification of Spinal Cord Injuries (ISNCSCI) examination.

RESULTS

CSA and LRW in the AIS A/B group were significantly lower than those in the TD group and the AIS C/D group. LRW was the most sensitive imaging biomarker to differentiate the AIS A/B group from the AIS C/D group. Both CSA and APW were positively correlated with ISNCSCI sensory scores.

CONCLUSIONS

Quantitative measurement of the morphologic spinal cord parameters of the cervical spinal cord can be used as an objective imaging biomarker to evaluate the neurological function of children with TLSCI. Cervical spinal cord atrophy in children after TLSCI was correlated with clinical grading; CSA and APW can reflect sensory function. Meanwhile, LRW has the potential to be an objective imaging biomarker for evaluating motor function preservation.

DeepRetroMoCo: deep neural network-based retrospective motion correction algorithm for spinal cord functional MRI

Background and purpose

There are distinct challenges in the preprocessing of spinal cord fMRI data, particularly concerning the mitigation of voluntary or involuntary movement artifacts during image acquisition. Despite the notable progress in data processing techniques for movement detection and correction, applying motion correction algorithms developed for the brain cortex to the brainstem and spinal cord remains a challenging endeavor.

Methods

In this study, we employed a deep learning-based convolutional neural network (CNN) named DeepRetroMoCo, trained using an unsupervised learning algorithm. Our goal was to detect and rectify motion artifacts in axial T2*-weighted spinal cord data. The training dataset consisted of spinal cord fMRI data from 27 participants, comprising 135 runs for training and 81 runs for testing.

Results

To evaluate the efficacy of DeepRetroMoCo, we compared its performance against the sct_fmri_moco method implemented in the spinal cord toolbox. We assessed the motion-corrected images using two metrics: the average temporal signal-to-noise ratio (tSNR) and Delta Variation Signal (DVARS) for both raw and motion-corrected data. Notably, the average tSNR in the cervical cord was significantly higher when DeepRetroMoCo was utilized for motion correction, compared to the sct_fmri_moco method. Additionally, the average DVARS values were lower in images corrected by DeepRetroMoCo, indicating a superior reduction in motion artifacts. Moreover, DeepRetroMoCo exhibited a significantly shorter processing time compared to sct_fmri_moco.

Conclusion

Our findings strongly support the notion that DeepRetroMoCo represents a substantial improvement in motion correction procedures for fMRI data acquired from the cervical spinal cord. This novel deep learning-based approach showcases enhanced performance, offering a promising solution to address the challenges posed by motion artifacts in spinal cord fMRI data.

Reliability of task-based fMRI in the dorsal horn of the human spinal cord

The application of functional magnetic resonance imaging (fMRI) to the human spinal cord is still a relatively small field of research and faces many challenges. Here we aimed to probe the limitations of task-based spinal fMRI at 3T by investigating the reliability of spinal cord blood oxygen level dependent (BOLD) responses to repeated nociceptive stimulation across two consecutive days in 40 healthy volunteers. We assessed the test-retest reliability of subjective ratings, autonomic responses, and spinal cord BOLD responses to short heat pain stimuli (1s duration) using the intraclass correlation coefficient (ICC). At the group level, we observed robust autonomic responses as well as spatially specific spinal cord BOLD responses at the expected location, but no spatial overlap in BOLD response patterns across days. While autonomic indicators of pain processing showed good-to-excellent reliability, both β-estimates and z-scores of task-related BOLD responses showed poor reliability across days in the target region (gray matter of the ipsilateral dorsal horn). When taking into account the sensitivity of gradient-echo echo planar imaging (GE-EPI) to draining vein signals by including the venous plexus in the analysis, we observed BOLD responses with fair reliability across days. Taken together, these results demonstrate that heat pain stimuli as short as one second are able to evoke a robust and spatially specific BOLD response, which is however strongly variable within participants across time, resulting in low reliability in the dorsal horn gray matter. Further improvements in data acquisition and analysis techniques are thus necessary before event-related spinal cord fMRI as used here can be reliably employed in longitudinal designs or clinical settings.

Genotype-specific spinal cord damage in spinocerebellar ataxias: an ENIGMA-Ataxia study

BACKGROUND

Spinal cord damage is a feature of many spinocerebellar ataxias (SCAs), but well-powered in vivo studies are lacking and links with disease severity and progression remain unclear. Here we characterise cervical spinal cord morphometric abnormalities in SCA1, SCA2, SCA3 and SCA6 using a large multisite MRI dataset.

METHODS

Upper spinal cord (vertebrae C1-C4) cross-sectional area (CSA) and eccentricity (flattening) were assessed using MRI data from nine sites within the ENIGMA-Ataxia consortium, including 364 people with ataxic SCA, 56 individuals with preataxic SCA and 394 nonataxic controls. Correlations and subgroup analyses within the SCA cohorts were undertaken based on disease duration and ataxia severity.

RESULTS

Individuals in the ataxic stage of SCA1, SCA2 and SCA3, relative to non-ataxic controls, had significantly reduced CSA and increased eccentricity at all examined levels. CSA showed large effect sizes (d>2.0) and correlated with ataxia severity (r<-0.43) and disease duration (r<-0.21). Eccentricity correlated only with ataxia severity in SCA2 (r=0.28). No significant spinal cord differences were evident in SCA6. In preataxic individuals, CSA was significantly reduced in SCA2 (d=1.6) and SCA3 (d=1.7), and the SCA2 group also showed increased eccentricity (d=1.1) relative to nonataxic controls. Subgroup analyses confirmed that CSA and eccentricity are abnormal in early disease stages in SCA1, SCA2 and SCA3. CSA declined with disease progression in all, whereas eccentricity progressed only in SCA2.

CONCLUSIONS

Spinal cord abnormalities are an early and progressive feature of SCA1, SCA2 and SCA3, but not SCA6, which can be captured using quantitative MRI.

Diffusion Tensor Imaging Identifies Cervical Spondylosis, Myelitis, and Spinal Cord Tumors

BACKGROUND

Diffusion tensor imaging (DTI) has been increasingly recognized for its capability to study microstructural changes in the neuropathology of brain diseases. However, the optimal DTI metric and its diagnostic utility for a variety of spinal cord diseases are still under investigation.

PURPOSE

To evaluate the diagnostic efficacy of DTI metrics for differentiating between cervical spondylosis, myelitis, and spinal tumors.

METHODS

This retrospective study analyzed DTI scans from 68 patients (22 with cervical spondylosis, 23 with myelitis, and 23 with spinal tumors). DTI indicators, including fractional anisotropy (FA), mean diffusivity (MD), radial diffusivity (RD) and axial diffusivity (AD), were calculated. The Kruskal-Wallis test was used to compare these indicators, followed by Receiver Operating Characteristic (ROC) curve analysis, to evaluate the diagnostic efficacy of each indicator across disease pairs. Additionally, we explored the correlations of DTI indicators with specific clinical measurements.

RESULTS

FA values were significantly lower in tumor patients compared to those with cervical spondylosis (p < 0.0001) and myelitis (p < 0.05). Additionally, tumor patients exhibited significantly elevated MD and RD values relative to the spondylosis and myelitis groups. ROC curve analysis underscored FA's superior discriminative performance, with an area under the curve (AUC) of 0.902 for differentiating tumors from cervical spondylosis, and an AUC of 0.748 for distinguishing cervical myelitis from spondylosis. Furthermore, a significant negative correlation was observed between FA values and Expanded Disability Status Scores (EDSSs) in myelitis patients (r = -0.62, p = 0.002), as well as between FA values and Ki-67 scores in tumor patients (r = -0.71, p = 0.0002).

CONCLUSION

DTI indicators, especially FA, have the potential in distinguishing spondylosis, myelitis, and spinal cord tumors. The significant correlation between FA values and clinical indicators highlights the value of FA in the clinical assessment and prognosis of spinal diseases and may be applied in diagnostic protocols in the future.

Resting-state brain plasticity is associated with the severity in cervical spondylotic myelopathy

OBJECTIVE

To investigate the brain mechanism of non-correspondence between imaging presentations and clinical symptoms in cervical spondylotic myelopathy (CSM) patients and to test the utility of brain imaging biomarkers for predicting prognosis of CSM.

METHODS

Forty patients with CSM (22 mild-moderate CSM, 18 severe CSM) and 25 healthy controls (HCs) were recruited for rs-fMRI and cervical spinal cord diffusion tensor imaging (DTI) scans. DTI at the spinal cord (level C2/3) with fractional anisotropy (FA) and degree centrality (DC) were recorded. Then one-way analysis of covariance (ANCOVA) was conducted to detect the group differences in the DC and FA values across the three groups. Pearson correlation analysis was then separately performed between JOA with FA and DC.

RESULTS

Among them, degree centrality value of left middle temporal gyrus exhibited a progressive increase in CSM groups compared with HCs, the DC value in severe CSM group was higher compared with mild-moderate CSM group. (P < 0.05), and the DC values of the right superior temporal gyrus and precuneus showed a decrease after increase. Among them, DC values in the area of precuneus in severe CSM group were significantly lower than those in mild-moderate CSM and HCs. (P < 0.05). The fractional anisotropy (FA) values of the level C2/3 showed a progressive decrease in different clinical stages, that severe CSM group was the lowest, significantly lower than those in mild-moderate CSM and HCs (P < 0.05). There was negative correlation between DC value of left middle temporal gyrus and JOA scores (P < 0.001), and the FA values of dorsal column in the level C2/3 positively correlated with the JOA scores (P < 0.001).

CONCLUSION

Structural and functional changes have taken place in the cervical spinal cord and brain of CSM patients. The Brain reorganization plays an important role in maintaining the symptoms and signs of CSM, aberrant DC values in the left middle temporal gyrus may be the possible mechanism of inconsistency between imaging findings and clinical symptoms. Degree centrality is a potentially useful prognostic functional biomarker in cervical spondylotic myelopathy.

Diffusion tensor imaging within the healthy cervical spinal cord: Within- participants reliability and measurement error

BACKGROUND

Diffusion tensor imaging (DTI) is a promising technique for the visualization of the cervical spinal cord (CSC) in vivo. It provides information about the tissue structure of axonal white matter, and it is thought to be more sensitive than other MR imaging techniques for the evaluation of damage to tracts in the spinal cord.

AIM

The purpose of this study was to determine the within-participants reliability and error magnitude of measurements of DTI metrics in healthy human CSC.

METHODS

A total of twenty healthy controls (10 male, mean age: 33.9 ± 3.5 years, 10 females, mean age: 47.5 ± 14.4 years), with no family history of any neurological disorders or a contraindication to MRI scanning were recruited over a period of two months. Each participant was scanned twice with an MRI 3 T scanner using standard DTI sequences. Spinal Cord Toolbox (SCT) software was used for image post-processing. Data were first corrected for motion artefact, then segmented, registered to a template, and then the DTI metrics were computed. The within-participants coefficients of variation (CV%), the single and average within-participants intraclass correlation coefficients (ICC) and Bland-Altman plots for WM, VC, DC and LC fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD), and radial diffusivity (RD) were determined for the cervical spinal cord (between the 2nd and 5th cervical vertebrae).

RESULTS

DTI metrics showed poor to excellent within-participants reliability for both single and average ICC and moderate to high reproducibility for CV%, all variation dependent on the location of the ROI. The BA plots showed good within-participants agreement between the scan-rescan values.

CONCLUSION

Results from this reliability study demonstrate that clinical trials using the DTI technique are feasible and that DTI, in particular regions of the cord is suitable for use for the monitoring of degenerative WM changes.

MR Imaging in Ataxias: Consensus Recommendations by the Ataxia Global Initiative Working Group on MRI Biomarkers

With many viable strategies in the therapeutic pipeline, upcoming clinical trials in hereditary and sporadic degenerative ataxias will benefit from non-invasive MRI biomarkers for patient stratification and the evaluation of therapies. The MRI Biomarkers Working Group of the Ataxia Global Initiative therefore devised guidelines to facilitate harmonized MRI data acquisition in clinical research and trials in ataxias. Recommendations are provided for a basic structural MRI protocol that can be used for clinical care and for an advanced multi-modal MRI protocol relevant for research and trial settings. The advanced protocol consists of modalities with demonstrated utility for tracking brain changes in degenerative ataxias and includes structural MRI, magnetic resonance spectroscopy, diffusion MRI, quantitative susceptibility mapping, and resting-state functional MRI. Acceptable ranges of acquisition parameters are provided to accommodate diverse scanner hardware in research and clinical contexts while maintaining a minimum standard of data quality. Important technical considerations in setting up an advanced multi-modal protocol are outlined, including the order of pulse sequences, and example software packages commonly used for data analysis are provided. Outcome measures most relevant for ataxias are highlighted with use cases from recent ataxia literature. Finally, to facilitate access to the recommendations by the ataxia clinical and research community, examples of datasets collected with the recommended parameters are provided and platform-specific protocols are shared via the Open Science Framework.

Cervical spinal cord angiography and vessel-selective perfusion imaging in the rat

Arterial spin labeling (ASL) has been widely used to evaluate arterial blood and perfusion dynamics, particularly in the brain, but its application to the spinal cord has been limited. The purpose of this study was to optimize vessel-selective pseudocontinuous arterial spin labeling (pCASL) for angiographic and perfusion imaging of the rat cervical spinal cord. A pCASL preparation module was combined with a train of gradient echoes for dynamic angiography. The effects of the echo train flip angle, label duration, and a Cartesian or radial readout were compared to examine their effects on visualizing the segmental arteries and anterior spinal artery (ASA) that supply the spinal cord. Lastly, vessel-selective encoding with either vessel-encoded pCASL (VE-pCASL) or super-selective pCASL (SS-pCASL) were compared. Vascular territory maps were obtained with VE-pCASL perfusion imaging of the spinal cord, and the interanimal variability was evaluated. The results demonstrated that longer label durations (200 ms) resulted in greater signal-to-noise ratio in the vertebral arteries, improved the conspicuity of the ASA, and produced better quality maps of blood arrival times. Cartesian and radial readouts demonstrated similar image quality. Both VE-pCASL and SS-pCASL adequately labeled the right or left vertebral arteries, which revealed the interanimal variability in the segmental artery with variations in their location, number, and laterality. VE-pCASL also demonstrated unique interanimal variations in spinal cord perfusion with a right-sided dominance across the six animals. Vessel-selective pCASL successfully achieved visualization of the arterial inflow dynamics and corresponding perfusion territories of the spinal cord. These methodological developments provide unique insights into the interanimal variations in the arterial anatomy and dynamics of spinal cord perfusion.

Early neurological changes in aging cervical spine: insights from PROMIS mobility assessment

Many studies have shown that the prevalence of degenerative spinal cord compression increases with age. However, most cases at early stages are asymptomatic, and their diagnosis remains challenging. Asymptomatic cervical spinal cord compression (ASCC) patients are more likely to experience annular tears, herniated disks, and later develop symptomatic compression. Asymptomatic individuals do not typically undergo spinal cord imaging; therefore, an assessment test that is both sensitive and specific in diagnosing ASCC may be helpful. It has been demonstrated that the Patient Reported Outcome Measure Information System (PROMIS) mobility test is sensitive in detecting degenerative cervical myelopathy (DCM) symptoms. We investigated the use of the PROMIS mobility test in assessing clinical dysfunction in ASCC. In this study, 51 DCM patients and 42 age-matched healthy control (HC) were enrolled. The degree of cervical spinal cord compression was assessed using the high-resolution cervical spinal cord T2 Weighted (T2w) MRIs, which were available for 14 DCM patients. Measurements of the spinal cords anterior-posterior (AP) diameter at the region(s) that were visibly compressed as well as at different cervical spine levels were used to determine the degree of compression. The age-matched HC cohort had a similar MRI to establish the normal range for AP diameter. Twelve (12) participants in the HC cohort had MRI evidence of cervical spinal cord compression; these individuals were designated as the ASCC cohort. All participants completed the PROMIS mobility, PROMIS pain interference (PI), PROMIS upper extremity (UE), modified Japanese orthopedic association (mJOA), and neck disability index (NDI) scoring scales. We examined the correlation between the AP diameter measurements and the clinical assessment scores to determine their usefulness in the diagnosis of ASCC. Furthermore, we examine the sensitivity and specificity of PROMIS mobility test and mJOA. Compared to the HC group, the participants in the ASCC and DCM cohorts were significantly older (p = 0.006 and p < 0.0001, respectively). Age differences were not observed between ASCC and DCM (p > 0.999). Clinical scores between the ASCC and the HC group were not significantly different using the mJOA (p > 0.99), NDI (p > 0.99), PROMIS UE (p = 0.23), and PROMIS PI (p = 0.82). However, there were significant differences between the ASCC and HC in the PROMIS mobility score (p = 0.01). The spinal cord AP diameter and the PROMIS mobility score showed a significant correlation (r = 0.44, p = 0.002). Decreasing PROMIS mobility was significantly associated with a decrease in cervical spinal cord AP diameter independent of other assessment measures. PROMIS mobility score had a sensitivity of 77.3% and specificity of 79.4% compared to 59.1% and 88.2%, respectively, for mJOA in detecting cervical spinal cord compression. Certain elements of ASCC are not adequately captured with the traditional mJOA and NDI scales used in DCM evaluation. In contrast to other evaluation scales utilized in this investigation, PROMIS mobility score shows a significant association with the AP diameter of the cervical spinal cord, suggesting that it is a sensitive tool for identifying early disability associated with degenerative change in the aging spine. In a comparative analysis of PROMIS mobility test against the standard mJOA, the PROMIS mobility demonstrated higher sensitivity for detecting cervical spinal cord compression. These findings underscore the potential use of PROMIS mobility score in clinical evaluation of the aging spine.

Structural changes in spinal cord following optic neuritis: Insights from quantitative spinal MRI

OBJECTIVES

Previous studies have demonstrated that optic neuritis (ON) affects brain plasticity. However, whether ON affects the spinal cord remains unclear. We aimed to investigate the spinal cord changes in ON and their associations with disability.

METHODS

A total of 101 ON patients, and 41 healthy controls (HC) were retrospectively recruited. High-resolution imaging was conducted using a Magnetization Prepared Rapid Acquisition Gradient-Echo (MP-RAGE) sequence for T1-weighted images and an echo planar imaging (EPI) sequence for Diffusion Tensor Imaging (DTI) data collection. Additionally, patients' disability and cognitive impairment were evaluated using the Expanded Disability Status Scale (EDSS) and the Paced Auditory Serial Addition Test (PASAT), respectively. The quantitative spinal MRI was employed to examine the cross-sectional area (CSA) and diffusion indicators, with a specific focus on calculating the average values across the C2-C7 cervical spinal cord segments. CSA, fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD), and radial diffusivity (RD) were compared between groups. Correlation analyses were performed between CSA, diffusion indicators, and clinical variables.

RESULTS

No significant differences were found in CSA between ON patients and HCs. MD (p = 0.007) and RD (p = 0.018) were increased in ON patients compared with HCs, and AD was decreased in ON (p = 0.013). The AD values of the ON patients were significantly positively correlated with PASAT scores (r = 0.37, p < 0.001).

CONCLUSIONS

This study provided imaging evidence for DTI abnormalities in patients with ON. Spinal cord DTI can improve our knowledge of the path physiology of ON, and clinical progression.

Evidence-based commentary on the diagnosis, management, and further research of degenerative cervical spinal cord compression in the absence of clinical symptoms of myelopathy

Degenerative cervical myelopathy (DCM) represents the final consequence of a series of degenerative changes in the cervical spine, resulting in cervical spinal canal stenosis and mechanical stress on the cervical spinal cord. This process leads to subsequent pathophysiological processes in the spinal cord tissues. The primary mechanism of injury is degenerative compression of the cervical spinal cord, detectable by magnetic resonance imaging (MRI), serving as a hallmark for diagnosing DCM. However, the relative resilience of the cervical spinal cord to mechanical compression leads to clinical-radiological discordance, i.e., some individuals may exhibit MRI findings of DCC without the clinical signs and symptoms of myelopathy. This degenerative compression of the cervical spinal cord without clinical signs of myelopathy, potentially serving as a precursor to the development of DCM, remains a somewhat controversial topic. In this review article, we elaborate on and provide commentary on the terminology, epidemiology, natural course, diagnosis, predictive value, risks, and practical management of this condition-all of which are subjects of ongoing debate.

Neurodegeneration within the rostral spinal cord is associated with brain gray matter volume atrophy in the early stage of cervical spondylotic myelopathy

STUDY DESIGN

Case-control study.

OBJECTIVES

Investigating the association between neurodegeneration within rostral spinal cord and brain gray matter volume (GMV) and assessing the relationship between remote neurodegenerative changes and clinical outcomes at the early phase of Cervical Spondylotic Myelopathy (CSM).

SETTING

University/hospital.

METHODS

Using Spinal Cord Toolbox, spinal cord morphometrics (cross-sectional area [CSA], gray matter area [GMA], white matter area [WMA]) of 40 patients with CSM and 28 healthy controls (HCs) were computed and compared using two-sample t test. Brain GMV of the two groups was analyzed using voxel-based morphometry approach. Pearson's correlation between spinal cord morphometrics and altered brain GMV and Spearman's relationship between remote neurodegenerations and clinical outcomes were conducted in CSM group.

RESULTS

Compared to HCs, CSA and WMA at C2/3 and GMV in right postcentral gyrus (PoCG.R) and left supplementary motor area (SMA.L) were significantly decreased in patients with CSM. CSA and WMA at C2/3 were associated with GMV in SMA.L and MCG.R in patients with CSM. CSA at C2/3 and GMV in PoCG.R were related to modified Japanese Orthopedic Association score in patients with CSM.

CONCLUSIONS

The associations between CSA and WMA at C2/3 and GMV in SMA.L and MCG.R suggest a concordant change pattern and adaptive mechanisms for neuronal plasticity underlying remote neurodegeneration in early CSM. The atrophy of CSA at C2/3 and GMV loss in PoCG.R can serve as potential neuroimaging biomarkers of early structural changes within spinal cord and brain preceding marked clinical disabilities in patients with CSM.

Influence of preprocessing, distortion correction and cardiac triggering on the quality of diffusion MR images of spinal cord

Diffusion MRI of the spinal cord (SC) is susceptible to geometric distortion caused by field inhomogeneities, and prone to misalignment across time series and signal dropout caused by biological motion. Several modifications of image acquisition and image processing techniques have been introduced to overcome these artifacts, but their specific benefits are largely unproven and warrant further investigations. We aim to evaluate two specific aspects of image acquisition and processing that address image quality in diffusion studies of the spinal cord: susceptibility corrections to reduce geometric distortions, and cardiac triggering to minimize motion artifacts. First, we evaluate 4 distortion preprocessing strategies on 7 datasets of the cervical and lumbar SC and find that while distortion correction techniques increase geometric similarity to structural images, they are largely driven by the high-contrast cerebrospinal fluid, and do not consistently improve the geometry within the cord nor improve white-to-gray matter contrast. We recommend at a minimum to perform bulk-motion correction in preprocessing and posit that improvements/adaptations are needed for spinal cord distortion preprocessing algorithms, which are currently optimized and designed for brain imaging. Second, we design experiments to evaluate the impact of removing cardiac triggering. We show that when triggering is foregone, images are qualitatively similar to triggered sequences, do not have increased prevalence of artifacts, and result in similar diffusion tensor indices with similar reproducibility to triggered acquisitions. When triggering is removed, much shorter acquisitions are possible, which are also qualitatively and quantitatively similar to triggered sequences. We suggest that removing cardiac triggering for cervical SC diffusion can be a reasonable option to save time with minimal sacrifice to image quality.

Reproducibility of 3D pH-weighted chemical exchange saturation transfer contrast in the healthy cervical spinal cord

Spinal cord ischemia and hypoxia can be caused by compression, injury, and vascular alterations. Measuring ischemia and hypoxia directly in the spinal cord noninvasively remains challenging. Ischemia and hypoxia alter tissue pH, providing a physiologic parameter that may be more directly related to tissue viability. Chemical exchange saturation transfer (CEST) is an MRI contrast mechanism that can be made sensitive to pH. More specifically, amine/amide concentration independent detection (AACID) is a recently developed endogenous CEST contrast that has demonstrated sensitivity to intracellular pH at 9.4 T. The goal of this study was to evaluate the reproducibility of AACID CEST measurements at different levels of the healthy cervical spinal cord at 3.0 T incorporating B1 correction. Using a 3.0 T MRI scanner, two 3D CEST scans (saturation pulse train followed by a 3D snapshot gradient-echo readout) were performed on 12 healthy subjects approximately 10 days apart, with the CEST volume centered at the C4 level for all subjects. Scan-rescan reproducibility was evaluated by examining between and within-subject coefficients of variation (CVs) and absolute AACID value differences. The C4 level of the spinal cord demonstrated the lowest within-subject CVs (4.1%-4.3%), between-subject CVs (5.6%-6.3%), and absolute AACID percent difference (5.8-6.1%). The B1 correction scheme significantly improved reproducibility (adjusted p-value = 0.002) compared with the noncorrected data, suggesting that implementing B1 corrections in the spinal cord is beneficial. It was concluded that pH-weighted AACID measurements, incorporating B1-inhomogeneity correction, were reproducible within subjects along the healthy cervical spinal cord and that optimal image quality was achieved at the center of the 3D CEST volume.

Wallerian Degeneration Assessed by Multi-Modal Magnetic Resonance Imaging of Cervical Spinal Cord Is Associated With Neurological Impairment After Spinal Cord Injury

While Wallerian degeneration (WD) is a crucial pathological process induced with spinal cord injury (SCI), its underlying mechanisms is still understudied. In this study, we aim to assess structural alterations and clinical significance of WD in the cervical cord following SCI using multi-modal magnetic resonance imaging (MRI), which combines T2*-weighted imaging and diffusion tensor imaging (DTI). T2*-weighted images allow segmentation of anatomical structures and the detection of WD on macrostructural level. DTI, on the other hand, can identify the reduction in neuroaxonal integrity by measuring the diffusion of water molecules on the microstructural level. In this prospective study, 35 SCI patients (19 paraplegic and 16 tetraplegic patients) and 12 healthy controls were recruited between July 2020 and May 2022. The hyperintensity voxels in the dorsal column was manually labeled as WD on T2*-weighted images. The mean cross-sectional area (CSA) and mean DTI indexes of WD at the C2 level were calculated and compared between groups. Correlation analysis was used to determine the associations of the magnitude of WD with lesion characteristics and clinical outcomes. Compared with controls, SCI patients showed evident hyperintensity (35/35) and decreased neuroaxonal integrity (p < 0.05) within the dorsal column at the C2 level. A higher neurological level of injury was associated with a larger mean CSA and reduction in neuroaxonal integrity within WD (p < 0.05). Smaller total and dorsal tissue bridges were related to greater mean CSA and lower fractional anisotropy values in WD (p < 0.05), respectively. Moreover, SCI participants with significantly larger CSAs and significantly lower microstructural integrity had worse sensory outcomes (p < 0.05). This comprehensive evaluation of WD can help us better understand the mechanisms of WD, monitor progression, and assess the effectiveness of therapeutic interventions after SCI.

Body size interacts with the structure of the central nervous system: A multi-center in vivo neuroimaging study

Clinical research emphasizes the implementation of rigorous and reproducible study designs that rely on between-group matching or controlling for sources of biological variation such as subject's sex and age. However, corrections for body size (i.e. height and weight) are mostly lacking in clinical neuroimaging designs. This study investigates the importance of body size parameters in their relationship with spinal cord (SC) and brain magnetic resonance imaging (MRI) metrics. Data were derived from a cosmopolitan population of 267 healthy human adults (age 30.1±6.6 years old, 125 females). We show that body height correlated strongly or moderately with brain gray matter (GM) volume, cortical GM volume, total cerebellar volume, brainstem volume, and cross-sectional area (CSA) of cervical SC white matter (CSA-WM; 0.44≤r≤0.62). In comparison, age correlated weakly with cortical GM volume, precentral GM volume, and cortical thickness (-0.21≥r≥-0.27). Body weight correlated weakly with magnetization transfer ratio in the SC WM, dorsal columns, and lateral corticospinal tracts (-0.20≥r≥-0.23). Body weight further correlated weakly with the mean diffusivity derived from diffusion tensor imaging (DTI) in SC WM (r=-0.20) and dorsal columns (-0.21), but only in males. CSA-WM correlated strongly or moderately with brain volumes (0.39≤r≤0.64), and weakly with precentral gyrus thickness and DTI-based fractional anisotropy in SC dorsal columns and SC lateral corticospinal tracts (-0.22≥r≥-0.25). Linear mixture of sex and age explained 26±10% of data variance in brain volumetry and SC CSA. The amount of explained variance increased at 33±11% when body height was added into the mixture model. Age itself explained only 2±2% of such variance. In conclusion, body size is a significant biological variable. Along with sex and age, body size should therefore be included as a mandatory variable in the design of clinical neuroimaging studies examining SC and brain structure.

Neurodesk: an accessible, flexible and portable data analysis environment for reproducible neuroimaging

Neuroimaging research requires purpose-built analysis software, which is challenging to install and may produce different results across computing environments. The community-oriented, open-source Neurodesk platform ( https://www.neurodesk.org/ ) harnesses a comprehensive and growing suite of neuroimaging software containers. Neurodesk includes a browser-accessible virtual desktop, command-line interface and computational notebook compatibility, allowing for accessible, flexible, portable and fully reproducible neuroimaging analysis on personal workstations, high-performance computers and the cloud.

Advanced MRI metrics improve the prediction of baseline disease severity for individuals with degenerative cervical myelopathy

BACKGROUND CONTEXT

Degenerative cervical myelopathy (DCM) is the most common form of atraumatic spinal cord injury globally. Degeneration of spinal discs, bony osteophyte growth and ligament pathology results in physical compression of the spinal cord contributing to damage of white matter tracts and grey matter cellular populations. This results in an insidious neurological and functional decline in patients which can lead to paralysis. Magnetic resonance imaging (MRI) confirms the diagnosis of DCM and is a prerequisite to surgical intervention, the only known treatment for this disorder. Unfortunately, there is a weak correlation between features of current commonly acquired MRI scans ("community MRI, cMRI") and the degree of disability experienced by a patient.

PURPOSE

This study examines the predictive ability of current MRI sequences relative to "advanced MRI" (aMRI) metrics designed to detect evidence of spinal cord injury secondary to degenerative myelopathy. We hypothesize that the utilization of higher fidelity aMRI scans will increase the effectiveness of machine learning models predicting DCM severity and may ultimately lead to a more efficient protocol for identifying patients in need of surgical intervention.

STUDY DESIGN/SETTING

Single institution analysis of imaging registry of patients with DCM.

PATIENT SAMPLE

A total of 296 patients in the cMRI group and 228 patients in the aMRI group.

OUTCOME MEASURES

Physiologic measures: accuracy of machine learning algorithms to detect severity of DCM assessed clinically based on the modified Japanese Orthopedic Association (mJOA) scale.

METHODS

Patients enrolled in the Canadian Spine Outcomes Research Network registry with DCM were screened and 296 cervical spine MRIs acquired in cMRI were compared with 228 aMRI acquisitions. aMRI acquisitions consisted of diffusion tensor imaging, magnetization transfer, T2-weighted, and T2*-weighted images. The cMRI group consisted of only T2-weighted MRI scans. Various machine learning models were applied to both MRI groups to assess accuracy of prediction of baseline disease severity assessed clinically using the mJOA scale for cervical myelopathy.

RESULTS

Through the utilization of Random Forest Classifiers, disease severity was predicted with 41.8% accuracy in cMRI scans and 73.3% in the aMRI scans. Across different predictive model variations tested, the aMRI scans consistently produced higher prediction accuracies compared to the cMRI counterparts.

CONCLUSIONS

aMRI metrics perform better in machine learning models at predicting disease severity of patients with DCM. Continued work is needed to refine these models and address DCM severity class imbalance concerns, ultimately improving model confidence for clinical implementation.

SCIseg: Automatic Segmentation of T2-weighted Intramedullary Lesions in Spinal Cord Injury

Purpose

To develop a deep learning tool for the automatic segmentation of T2-weighted intramedullary lesions in spinal cord injury (SCI).

Material and Methods

This retrospective study included a cohort of SCI patients from three sites enrolled between July 2002 and February 2023. A deep learning model, SCIseg, was trained in a three-phase process involving active learning for the automatic segmentation of intramedullary SCI lesions and the spinal cord. The data consisted of T2-weighted MRI acquired using different scanner manufacturers with heterogeneous image resolutions (isotropic/anisotropic), orientations (axial/sagittal), lesion etiologies (traumatic/ischemic/hemorrhagic) and lesions spread across the cervical, thoracic and lumbar spine. The segmentations from the proposed model were visually and quantitatively compared with other open-source baselines. Wilcoxon signed-rank test was used to compare quantitative MRI biomarkers (lesion volume, lesion length, and maximal axial damage ratio) computed from manual lesion masks and those obtained automatically with SCIseg predictions.

Results

MRI data from 191 SCI patients (mean age, 48.1 years ± 17.9 [SD]; 142 males) were used for model training and evaluation. SCIseg achieved the best segmentation performance for both the cord and lesions. There was no statistically significant difference between lesion length and maximal axial damage ratio computed from manually annotated lesions and those obtained using SCIseg.

Conclusion

Automatic segmentation of intramedullary lesions commonly seen in SCI replaces the tedious manual annotation process and enables the extraction of relevant lesion morphometrics in large cohorts. The proposed model segments lesions across different etiologies, scanner manufacturers, and heterogeneous image resolutions. SCIseg is open-source and accessible through the Spinal Cord Toolbox.

Systemic inflammation associates with and precedes cord atrophy in progressive multiple sclerosis

In preclinical models of multiple sclerosis, systemic inflammation has an impact on the compartmentalized inflammatory process within the central nervous system and results in axonal loss. It remains to be shown whether this is the case in humans, specifically whether systemic inflammation contributes to spinal cord or brain atrophy in multiple sclerosis. Hence, an observational longitudinal study was conducted to delineate the relationship between systemic inflammation and atrophy using magnetic resonance imaging: the SIMS (Systemic Inflammation in Multiple Sclerosis) study. Systemic inflammation and progression were assessed in people with progressive multiple sclerosis (n = 50) over two and a half years. Eligibility criteria included: (i) primary or secondary progressive multiple sclerosis; (ii) age ≤ 70; and (iii) Expanded Disability Status Scale ≤ 6.5. First morning urine was collected weekly to quantify systemic inflammation by measuring the urinary neopterin-to-creatinine ratio using a validated ultra-performance liquid chromatography mass spectrometry technique. The urinary neopterin-to-creatinine ratio temporal profile was characterized by short-term responses overlaid on a background level of inflammation, so these two distinct processes were considered as separate variables: background inflammation and inflammatory response. Participants underwent MRI at the start and end of the study, to measure cervical spinal cord and brain atrophy. Brain and cervical cord atrophy occurred on the study, but the most striking change was seen in the cervical spinal cord, in keeping with the corticospinal tract involvement that is typical of progressive disease. Systemic inflammation predicted cervical cord atrophy. An association with brain atrophy was not observed in this cohort. A time lag between systemic inflammation and cord atrophy was evident, suggesting but not proving causation. The association of the inflammatory response with cord atrophy depended on the level of background inflammation, in keeping with experimental data in preclinical models where the effects of a systemic inflammatory challenge on tissue injury depended on prior exposure to inflammation. A higher inflammatory response was associated with accelerated cord atrophy in the presence of background systemic inflammation below the median for the study population. Higher background inflammation, while associated with cervical cord atrophy itself, subdued the association of the inflammatory response with cord atrophy. Findings were robust to sensitivity analyses adjusting for potential confounders and excluding cases with new lesion formation. In conclusion, systemic inflammation associates with, and precedes, multiple sclerosis progression. Further work is needed to prove causation since targeting systemic inflammation may offer novel treatment strategies for slowing neurodegeneration in multiple sclerosis.

Improved inter-subject alignment of the lumbosacral cord for group-level in vivo gray and white matter assessments: A scan-rescan MRI study at 3T

INTRODUCTION

Magnetic resonance imaging (MRI) enables the investigation of pathological changes in gray and white matter at the lumbosacral enlargement (LSE) and conus medullaris (CM). However, conducting group-level analyses of MRI metrics in the lumbosacral spinal cord is challenging due to variability in CM length, lack of established image-based landmarks, and unknown scan-rescan reliability. This study aimed to improve inter-subject alignment of the lumbosacral cord to facilitate group-level analyses of MRI metrics. Additionally, we evaluated the scan-rescan reliability of MRI-based cross-sectional area (CSA) measurements and diffusion tensor imaging (DTI) metrics.

METHODS

Fifteen participants (10 healthy volunteers and 5 patients with spinal cord injury) underwent axial T2*-weighted and diffusion MRI at 3T. We assessed the reliability of spinal cord and gray matter-based landmarks for inter-subject alignment of the lumbosacral cord, the inter-subject variability of MRI metrics before and after adjusting for the CM length, the intra- and inter-rater reliability of CSA measurements, and the scan-rescan reliability of CSA measurements and DTI metrics.

RESULTS

The slice with the largest gray matter CSA as an LSE landmark exhibited the highest reliability, both within and across raters. Adjusting for the CM length greatly reduced the inter-subject variability of MRI metrics. The intra-rater, inter-rater, and scan-rescan reliability of MRI metrics were the highest at and around the LSE (scan-rescan coefficient of variation <3% for CSA measurements and <7% for DTI metrics within the white matter) and decreased considerably caudal to it.

CONCLUSIONS

To facilitate group-level analyses, we recommend using the slice with the largest gray matter CSA as a reliable LSE landmark, along with an adjustment for the CM length. We also stress the significance of the anatomical location within the lumbosacral cord in relation to the reliability of MRI metrics. The scan-rescan reliability values serve as valuable guides for power and sample size calculations in future longitudinal studies.

Diffusion Weighted Magnetic Resonance Imaging of Spinal Cord Injuries After Instrumented Fusion Stabilization

Diffusion-weighted magnetic resonance imaging (DW-MRI) is a promising technique for assessing spinal cord injury (SCI) that has historically been challenged by the presence of metallic stabilization hardware. This study leverages recent advances in metal-artifact resistant multi-spectral DW-MRI to enable diffusion quantification throughout the spinal cord even after fusion stabilization. Twelve participants with cervical spinal cord injuries treated with fusion stabilization and 49 asymptomatic able-bodied control participants underwent multi-spectral DW-MRI evaluation. Apparent diffusion coefficient (ADC) values were calculated in axial cord sections. Statistical modeling assessed ADC differences across cohorts and within distinct cord regions of the SCI participants (at, above, or below injured level). Computed models accounted for subject demographics and injury characteristics. ADC was found to be elevated at injured levels compared with non-injured levels (z = 3.2, p = 0.001), with ADC at injured levels decreasing over time since injury (z = -9.2, p < 0.001). Below the injury level, ADC was reduced relative to controls (z = -4.4, p < 0.001), with greater reductions after more severe injuries that correlated with lower extremity motor scores (z = 2.56, p = 0.012). No statistically significant differences in ADC above the level of injury were identified. By enabling diffusion analysis near fusion hardware, the multi-spectral DW-MRI technique allowed intuitive quantification of cord diffusion changes after SCI both at and away from injured levels. This demonstrates the approach's potential for assessing post-surgical spinal cord integrity throughout stabilized regions.

What contributes to disability in progressive MS? A brain and cervical cord-matched quantitative MRI study

BACKGROUND

We assessed the ability of a brain-and-cord-matched quantitative magnetic resonance imaging (qMRI) protocol to differentiate patients with progressive multiple sclerosis (PMS) from controls, in terms of normal-appearing (NA) tissue abnormalities, and explain disability.

METHODS

A total of 27 patients and 16 controls were assessed on the Expanded Disability Status Scale (EDSS), 25-foot timed walk (TWT), 9-hole peg (9HPT) and symbol digit modalities (SDMT) tests. All underwent 3T brain and (C2-C3) cord structural imaging and qMRI (relaxometry, quantitative magnetisation transfer, multi-shell diffusion-weighted imaging), using a fast brain-and-cord-matched protocol with brain-and-cord-unified imaging readouts. Lesion and NA-tissue volumes and qMRI metrics reflecting demyelination and axonal loss were obtained. Random forest analyses identified the most relevant volumetric/qMRI measures to clinical outcomes. Confounder-adjusted linear regression estimated the actual MRI-clinical associations.

RESULTS

Several qMRI/volumetric differences between patients and controls were observed (p < 0.01). Higher NA-deep grey matter quantitative-T1 (EDSS: beta = 7.96, p = 0.006; 9HPT: beta = -0.09, p = 0.004), higher NA-white matter orientation dispersion index (TWT: beta = -3.21, p = 0.005; SDMT: beta = -847.10, p < 0.001), lower whole-cord bound pool fraction (9HPT: beta = 0.79, p = 0.001) and higher NA-cortical grey matter quantitative-T1 (SDMT = -94.31, p < 0.001) emerged as particularly relevant predictors of greater disability.

CONCLUSION

Fast brain-and-cord-matched qMRI protocols are feasible and identify demyelination - combined with other mechanisms - as key for disability accumulation in PMS.

Association of Mean Upper Cervical Spinal Cord Cross-Sectional Area With Cerebral Small Vessel Disease: A Community-Based Cohort Study

BACKGROUND

The purpose of this study was to investigate the association between the mean upper cervical spinal cord cross-sectional area (MUCCA) and the risk and severity of cerebral small vessel disease (CSVD).

METHODS

Community-dwelling residents in Lishui City, China, from the cross-sectional survey in the PRECISE cohort study (Polyvascular Evaluation for Cognitive Impairment and Vascular Events) conducted from 2017 to 2019. We included 1644 of 3067 community-dwelling adults in the PRECISE study after excluding those with incorrect, incomplete, insufficient, or missing clinical or imaging data. Total and modified total CSVD scores, as well as magnetic resonance imaging features, including white matter hyperintensity, lacunes, cerebral microbleeds, enlarged perivascular spaces, and brain atrophy, were assessed at the baseline. The Spinal Cord Toolbox was used to measure the upper cervical spinal cord cross-sectional area of the C1 to C3 segments of the spinal cord and its average value was taken as MUCCA. Participants were divided into 4 groups according to quartiles of MUCCA. Associations were analyzed using linear regression models adjusted for age, sex, current smoking and drinking, medical history, intracranial volume, and total cortical volume.

RESULTS

The means±SD age of the participants was 61.4±6.5 years, and 635 of 1644 participants (38.6%) were men. The MUCCA was smaller in patients with CSVD than those without CSVD. Using the total CSVD score as a criterion, the MUCCA was 61.78±6.12 cm2 in 504 of 1644 participants with CSVD and 62.74±5.94 cm2 in 1140 of 1644 participants without CSVD. Using the modified total CSVD score, the MUCCA was 61.81±6.04 cm2 in 699 of 1644 participants with CSVD and 62.91±5.94 cm2 in 945 of 1644 without CSVD. There were statistical differences between the 2 groups after adjusting for covariates in 3 models. The MUCCA was negatively associated with the total and modified total CSVD scores (adjusted β value, -0.009 [95% CI, -0.01 to -0.003] and -0.007 [95% CI, -0.01 to -0.0006]) after adjustment for covariates. Furthermore, the MUCCA was negatively associated with the white matter hyperintensity burden (adjusted β value, -0.01 [95% CI, -0.02 to -0.003]), enlarged perivascular spaces in the basal ganglia (adjusted β value, -0.005 [95% CI, -0.009 to -0.001]), lacunes (adjusted β value, -0.004 [95% CI, -0.007 to -0.0007]), and brain atrophy (adjusted β value, -0.009 [95% CI, -0.01 to -0.004]).

CONCLUSIONS

The MUCCA and CSVD were correlated. Spinal cord atrophy may serve as an imaging marker for CSVD; thus, small vessel disease may involve the spinal cord in addition to being intracranial.

Quantitative brainstem and spinal MRI in amyotrophic lateral sclerosis: implications for predicting noninvasive ventilation needs

BACKGROUND

Respiratory complications resulting from motor neurons degeneration are the primary cause of death in amyotrophic lateral sclerosis (ALS). Predicting the need for non-invasive ventilation (NIV) in ALS is important for advance care planning and clinical trial design. The aim of this study was to assess the potential of quantitative MRI at the brainstem and spinal cord levels to predict the need for NIV during the first six months after diagnosis.

METHODS

Forty-one ALS patients underwent MRI and spirometry shortly after diagnosis. The need for NIV was monitored according to French health guidelines for 6 months. The performance of four regression models based on: clinical variables, brainstem structures volumes, cervical spinal measurements, and combined variables were compared to predict the need for NIV within this period.

RESULTS

Both the clinical model (R2 = 0.28, AUC = 0.85, AICc = 42.67, BIC = 49.8) and the brainstem structures' volumes model (R2 = 0.30, AUC = 0.85, AICc = 40.13, BIC = 46.99) demonstrated good predictive performance. In addition, cervical spinal cord measurements model similar performance (R2 = 0.338, AUC = 0.87, AICc = 37.99, BIC = 44.49). Notably, the combined model incorporating predictors from all three models yielded the best performance (R2 = 0.60, AUC = 0.959, AICc = 36.38, BIC = 44.8). These findings are supported by observed positive correlations between brainstem volumes, cervical (C4/C7) cross-sectional area, and spirometry-measured lung volumes.

CONCLUSIONS

Our study shows that brainstem volumes and spinal cord area are promising measures to predict respiratory intervention needs in ALS.

Diagnostic efficacy of tract-specific diffusion tensor imaging in cervical spondylotic myelopathy with electrophysiological examination validation

PURPOSE

This study aimed to investigate the effectiveness of tract-specific diffusion tensor imaging (DTI) metrics in identifying the responsible segments for neurological dysfunction in cervical spondylotic myelopathy (CSM).

METHODS

The study encompassed nineteen participants diagnosed with CSM, including 10 males and 9 females. Additionally, a control group consisting of ten healthy caregivers (5 males and 5 females) were recruited with no symptoms and no compressions on magnetic resonance imaging (MRI). All participants underwent a comprehensive physical examination, MRI assessment, and DTI examination conducted by a senior chief physician. Several parameters were collected from the MR images, including the aspect ratio (defined as the anteroposterior diameter / the transverse diameter of the corresponding segment's spinal cord), transverse ratio (defined as the transverse diameter of the corresponding segment's spinal cord / the transverse diameter of the spinal cord at C2/3), and T2 high signal of the spinal cord. Furthermore, quantitative DTI metrics, such as axial diffusivity (AD), mean diffusivity (MD), radial diffusivity (RD), and fractional anisotropy (FA), were calculated using automatic region-of-interest (ROI) analysis for both whole spinal cord column and dorsal column. Receiver operating characteristic (ROC) curves were constructed to evaluate the diagnostic efficacy of the aspect ratio, transverse ratio, and DTI parameters. The area under the curve (AUC), sensitivity, and specificity were calculated. Intraoperative spinal cord electrophysiological examination was performed as the objective measure of spinal cord function during surgery.

RESULTS

As determined by electrophysiological examination, neurological dysfunction was found in 2 patients due to C3/4 compression, in 10 patients due to C4/5 compression, in 6 patients due to C5/6 compression, and in 1 patient due to C6/7 compression. The modified Japanese Orthopedic Association scale (mJOA) was 12.71 ± 1.55 in the CSM group, with 4.87 ± 0.72 for sensory nerve function and 5.05 ± 1.35 for motor nerve function. For the control group, none of the volunteers had neurological dysfunction. T2 high signal was found at the most stenotic segment in 13 patients of the CSM group. Considering all the cervical segments, the aspect ratio (AUC = 0.823, P = 0.001, Sensitivity = 68.42%, Specificity = 82.47%) was more capable of determining the responsible segment than transverse ratio (AUC = 0.661, P = 0.027, Sensitivity = 68.42%, Specificity = 67.01%). AD, MD, and RD were significantly higher while FA was significantly lower in the responsible segment than in the irresponsible segment (P < 0.05). The AUC of DTI-Dorsal column parameters (AD, MD, RD, FA) was larger than the corresponding parameters of the DTI (Whole spinal cord). AD of DTI-Dorsal Column possessed the greatest efficacy (AUC = 0.823, sensitivity = 84.21%, specificity = 77.32%) to determine the responsible segment, larger than AD of DTI-Whole spinal cord (AUC = 0.822, P = 0.001, Sensitivity = 89.47%, Specificity = 77.32%), aspect ratio (AUC = 0.823, P = 0.001, Sensitivity = 68.42%, Specificity = 82.47%) and transverse ratio (AUC = 0.661, P = 0.027, Sensitivity = 68.42%, Specificity = 67.01%). Subgroup analysis revealed that the diagnostic efficacy of DTI and MRI parameters was influenced by cervical spine segment.

CONCLUSIONS

When considering all cervical segments, AD from the DTI-Dorsal Column exhibited the most significant potential in identifying responsible segments. This potential was found to be superior to that of DTI-Whole spinal cord, aspect ratio, the most stenotic segment, T2 high signals, transverse ratio, motor nerve dysfunction, and sensory nerve dysfunction. The diagnostic effectiveness of both DTI and MRI parameters was notably influenced by the specific cervical spine segment.

Practical Applications of Artificial Intelligence in Spine Imaging: A Review

Artificial intelligence (AI), a transformative technology with unprecedented potential in medical imaging, can be applied to various spinal pathologies. AI-based approaches may improve imaging efficiency, diagnostic accuracy, and interpretation, which is essential for positive patient outcomes. This review explores AI algorithms, techniques, and applications in spine imaging, highlighting diagnostic impact and challenges with future directions for integrating AI into spine imaging workflow.

MRI mapping of hemodynamics in the human spinal cord

Impaired spinal cord vascular function contributes to numerous neurological pathologies, making it important to be able to noninvasively characterize these changes. Here, we propose a functional magnetic resonance imaging (fMRI)-based method to map spinal cord vascular reactivity (SCVR). We used a hypercapnic breath-holding task, monitored with end-tidal CO2 (PETCO2), to evoke a systemic vasodilatory response during concurrent blood oxygenation level-dependent (BOLD) fMRI. SCVR amplitude and hemodynamic delay were mapped at the group level in 27 healthy participants as proof-of-concept of the approach, and then in two highly-sampled participants to probe feasibility/stability of individual SCVR mapping. Across the group and the highly-sampled individuals, a strong ventral SCVR amplitude was initially observed without accounting for local regional variation in the timing of the vasodilatory response. Shifted breathing traces (PETCO2) were used to account for temporal differences in the vasodilatory response across the spinal cord, producing maps of SCVR delay. These delay maps reveal an earlier ventral and later dorsal response and demonstrate distinct gray matter regions concordant with territories of arterial supply. The SCVR fMRI methods described here enable robust mapping of spatiotemporal hemodynamic properties of the human spinal cord. This noninvasive approach has exciting potential to provide early insight into pathology-driven vascular changes in the cord, which may precede and predict future irreversible tissue damage and guide the treatment of several neurological pathologies involving the spine.

Thermal stimulus task fMRI in the cervical spinal cord at 7 Tesla

Although functional magnetic resonance imaging (fMRI) 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 blood oxygenation level-dependent (BOLD) signal, but amplifies the effects of B0 inhomogeneity. In this work, we present the first task fMRI in the spinal cord at 7 T. Further, we compare the performance of single-shot and multi-shot 2D echo-planar imaging (EPI) protocols, which differ in sensitivity to spatial and temporal B0 inhomogeneity. 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. All protocols used 3 mm slice thickness. For each protocol, the BOLD response to 13 10-s noxious thermal stimuli applied to the right thumb was acquired in a 10-min fMRI run. Image quality, temporal signal to noise ratio (SNR), and BOLD activation (percent signal change and z-stat) at both individual- and group-level were evaluated between the protocols. 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. 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. PRACTITIONER POINTS: First stimulus task fMRI results in the spinal cord at 7 T. Single-shot 0.75 mm 2D EPI produced the highest mean z-statistic. Multi-shot 0.60 mm 2D EPI provided the best-localized activation and least distortion.

Moving toward elucidating alternative motor pathway structures post-stroke: the value of spinal cord neuroimaging

Stroke results in varying levels of motor and sensory disability that have been linked to the neurodegeneration and neuroinflammation that occur in the infarct and peri-infarct regions within the brain. Specifically, previous research has identified a key role of the corticospinal tract in motor dysfunction and motor recovery post-stroke. Of note, neuroimaging studies have utilized magnetic resonance imaging (MRI) of the brain to describe the timeline of neurodegeneration of the corticospinal tract in tandem with motor function following a stroke. However, research has suggested that alternate motor pathways may also underlie disease progression and the degree of functional recovery post-stroke. Here, we assert that expanding neuroimaging techniques beyond the brain could expand our knowledge of alternate motor pathway structure post-stroke. In the present work, we will highlight findings that suggest that alternate motor pathways contribute to post-stroke motor dysfunction and recovery, such as the reticulospinal and rubrospinal tract. Then we review imaging and electrophysiological techniques that evaluate alternate motor pathways in populations of stroke and other neurodegenerative disorders. We will then outline and describe spinal cord neuroimaging techniques being used in other neurodegenerative disorders that may provide insight into alternate motor pathways post-stroke.

Spinal fMRI demonstrates segmental organisation of functionally connected networks in the cervical spinal cord: A test-retest reliability study

Resting functional magnetic resonance imaging (fMRI) studies have identified intrinsic spinal cord activity, which forms organised motor (ventral) and sensory (dorsal) resting-state networks. However, to facilitate the use of spinal fMRI in, for example, clinical studies, it is crucial to first assess the reliability of the method, particularly given the unique anatomical, physiological, and methodological challenges associated with acquiring the data. Here, we characterise functional connectivity relationships in the cervical cord and assess their between-session test-retest reliability in 23 young healthy volunteers. Resting-state networks were estimated in two ways (1) by estimating seed-to-voxel connectivity maps and (2) by calculating seed-to-seed correlations. Seed regions corresponded to the four grey matter horns (ventral/dorsal and left/right) of C5-C8 segmental levels. Test-retest reliability was assessed using the intraclass correlation coefficient. Spatial overlap of clusters derived from seed-to-voxel analysis between sessions was examined using Dice coefficients. Following seed-to-voxel analysis, we observed distinct unilateral dorsal and ventral organisation of cervical spinal resting-state networks that was largely confined in the rostro-caudal extent to each spinal segmental level, with more sparse connections observed between segments. Additionally, strongest correlations were observed between within-segment ipsilateral dorsal-ventral connections, followed by within-segment dorso-dorsal and ventro-ventral connections. Test-retest reliability of these networks was mixed. Reliability was poor when assessed on a voxelwise level, with more promising indications of reliability when examining the average signal within clusters. Reliability of correlation strength between seeds was highly variable, with the highest reliability achieved in ipsilateral dorsal-ventral and dorso-dorsal/ventro-ventral connectivity. However, the spatial overlap of networks between sessions was excellent. We demonstrate that while test-retest reliability of cervical spinal resting-state networks is mixed, their spatial extent is similar across sessions, suggesting that these networks are characterised by a consistent spatial representation over time.

Spinal cord perfusion is associated with microstructural damage in cervical spondylotic myelopathy patients who underwent cervical laminoplasty

OBJECTIVES

To investigate the association between spinal cord perfusion and microstructural damage in CSM patients who underwent cervical laminoplasty using MR dynamic susceptibility contrast (DSC), diffusion tensor imaging (DTI), and neurite orientation dispersion and density imaging (NODDI) techniques.

METHODS

A follow-up cohort study was conducted with 53 consecutively recruited CSM patients who had undergone cervical laminoplasty 12-14 months after the surgery from April 2016 to December 2016. Twenty-one aged-matched healthy volunteers were recruited as controls. For each patient, decompressed spinal cord levels were imaged on a 3.0-T MRI scanner by diffusion and DSC sequences to quantify the degrees of microstructural damage and perfusion conditions, respectively. The diffusion data were analyzed by DTI and NODDI models to produce diffusion metrics. Classic indicator dilution model was used to quantify the DSC metrics. Mann-Whitney U test was performed for comparison of diffusion metrics between patients and healthy controls. Pearson correlation was used to explore the associations between the metrics of spinal cord perfusion and microstructural damage.

RESULTS

DTI metrics, neurite density, and isotropic volume fraction had significant differences between postoperative patients and healthy controls. Pearson correlation test showed that SCBV was significantly positively correlated with RD, MD, and ODI, and negatively correlated with FA and NDI. SCBF was found to be significantly positively correlated with RD and MD, and negatively correlated with FA.

CONCLUSIONS

Increased spinal cord perfusion quantified by DSC is associated with microstructural damage assessed by diffusion MRI in CSM patients who underwent cervical laminoplasty.

CLINICAL RELEVANCE STATEMENT

This study found that the spinal cord perfusion is associated with microstructural damage in postoperative cervical spondylotic myelopathy patients, indicating that high perfusion may play a role in the pathophysiological process of cervical spondylotic myelopathy and deserves more attention.

KEY POINTS

• Spinal cord microstructural damage can be persistent despite the compression had been relieved 12-14 months after the cervical laminoplasty in cervical spondylotic myelopathy (CSM) patients. • Spinal cord perfusion is associated with microstructural damage in CSM patients after the cervical laminoplasty. • Inflammation in the decompressed spinal cord may be a cause of increased perfusion and is associated with microstructural damage during the recovery period of CSM.

MSEF-Net: Multi-scale edge fusion network for lumbosacral plexus segmentation with MR image

Nerve damage of spine areas is a common cause of disability and paralysis. The lumbosacral plexus segmentation from magnetic resonance imaging (MRI) scans plays an important role in many computer-aided diagnoses and surgery of spinal nerve lesions. Due to the complex structure and low contrast of the lumbosacral plexus, it is difficult to delineate the regions of edges accurately. To address this issue, we propose a Multi-Scale Edge Fusion Network (MSEF-Net) to fully enhance the edge feature in the encoder and adaptively fuse multi-scale features in the decoder. Specifically, to highlight the edge structure feature, we propose an edge feature fusion module (EFFM) by combining the Sobel operator edge detection and the edge-guided attention module (EAM), respectively. To adaptively fuse the multi-scale feature map in the decoder, we introduce an adaptive multi-scale fusion module (AMSF). Our proposed MSEF-Net method was evaluated on the collected spinal MRI dataset with 89 patients (a total of 2848 MR images). Experimental results demonstrate that our MSEF-Net is effective for lumbosacral plexus segmentation with MR images, when compared with several state-of-the-art segmentation methods.

FASB: an integrated processing pipeline for Functional Analysis of simultaneous Spinal cord-Brain fMRI

Simultaneous functional magnetic resonance imaging (fMRI) of the spinal cord and brain represents a powerful method for examining both ascending sensory and descending motor pathways in humans in vivo . However, its image acquisition protocols, and processing pipeline are less well established. This limitation is mainly due to technical difficulties related to spinal cord fMRI, and problems with the logistics stemming from a large field of view covering both brain and cervical cord. Here, we propose an acquisition protocol optimized for both anatomical and functional images, as well as an optimized integrated image processing pipeline, which consists of a novel approach for automatic modeling and mitigating the negative impact of spinal voxels with low temporal signal to noise ratio (tSNR). We validate our integrated pipeline, named FASB, using simultaneous fMRI data acquired during the performance of a motor task, as well as during resting-state conditions. We demonstrate that FASB outperforms the current spinal fMRI processing methods in three domains, including motion correction, registration to the spinal cord template, and improved detection power of the group-level analysis by removing the effects of participant-specific low tSNR voxels, typically observed at the disk level. Using FASB, we identify significant task-based activations in the expected sensorimotor network associated with a unilateral handgrip force production task across the entire central nervous system, including the contralateral sensorimotor cortex, thalamus, striatum, cerebellum, brainstem, as well as ipsilateral ventral horn at C5-C8 cervical levels. Additionally, our results show significant task-based functional connectivity between the key sensory and motor brain areas and the dorsal and ventral horns of the cervical cord. Overall, our proposed acquisition protocol and processing pipeline provide a robust method for characterizing the activation and functional connectivity of distinct cortical, subcortical, brainstem and spinal cord regions in humans.