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

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Recent developments and future avenues for human corticospinal neuroimaging

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

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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

STUDY DESIGN

Retrospective case-control study.

OBJECTIVES

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

SETTING

Novosibirsk Neurosurgery Center, Russia.

METHODS

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

RESULTS

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

CONCLUSIONS

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

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 good 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.

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

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

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

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

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

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

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

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

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

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

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

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

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

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

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