Framework for integrated MRI average of the spinal cord white and gray matter: the MNI-Poly-AMU template

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.

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.

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.

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

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

Comparison of spinal magnetic resonance imaging and classical clinical factors in predicting motor capacity in amyotrophic lateral sclerosis

BACKGROUND

Motor capacity is crucial in amyotrophic lateral sclerosis (ALS) clinical trial design and patient care. However, few studies have explored the potential of multimodal MRI to predict motor capacity in ALS. This study aims to evaluate the predictive value of cervical spinal cord MRI parameters for motor capacity in ALS compared to clinical prognostic factors.

METHODS

Spinal multimodal MRI was performed shortly after diagnosis in 41 ALS patients and 12 healthy participants as part of a prospective multicenter cohort study, the PULSE study (NCT00002013-A00969-36). Motor capacity was assessed using ALSFRS-R scores. Multiple stepwise linear regression models were constructed to predict motor capacity at 3 and 6 months from diagnosis, based on clinical variables, structural MRI measurements, including spinal cord cross-sectional area (CSA), anterior-posterior, and left-to-right cross-section diameters at vertebral levels from C1 to T4, and diffusion parameters in the lateral corticospinal tracts (LCSTs) and dorsal columns.

RESULTS

Structural MRI measurements were significantly correlated with the ALSFRS-R score and its sub-scores. And as early as 3 months from diagnosis, structural MRI measurements fit the best multiple linear regression model to predict the total ALSFRS-R (R² = 0.70, p value = 0.0001) and arm sub-score (R² = 0.69, p value = 0.0002), and combined with DTI metric in the LCST and clinical factors fit the best multiple linear regression model to predict leg sub-score (R² = 0.73, p value = 0.0002).

CONCLUSIONS

Spinal multimodal MRI could be promising as a tool to enhance prognostic accuracy and serve as a motor function proxy in ALS.

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

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

Quantitative comparison of corticospinal tracts arising from different cortical areas in humans

The corticospinal tract (CST), which plays a major role in the control of voluntary limb movements, arises from multiple motor- and somatosensory-related areas in monkeys. Although the cortical origin and quantitative differences in CSTs among the cortical areas are well-documented in monkeys, they are unclear in humans. We quantitatively investigated the CSTs from the cerebral cortex to the cervical cord in healthy volunteers using fiber tractography of diffusion-weighted magnetic resonance imaging. The corticospinal (CS) streamlines arose from nine cortical areas: primary motor area (mean ± SD = 49.71 ± 1.61%), dorsal (16.33 ± 1.37%) and ventral (11.02 ± 0.90%) premotor cortex, supplementary motor area (5.14 ± 0.36%), pre-supplementary motor area (2.46 ± 0.26%), primary somatosensory cortex (11.06 ± 0.91%), Brodmann area 5 (0.88 ± 0.09%), caudal cingulate zone (1.70 ± 0.30%), and posterior part of the rostral cingulate zone (1.70 ± 0.34%). In all cortical areas, the number of CS streamlines gradually decreased from the rostral to caudal spinal segments, but the proportion was maintained throughout the cervical cord. Over 75% of CS streamlines arose from the lateral surface of the frontal lobe, which may explain the voluntary control of dexterous and flexible limb movements in humans. (197/200 words).

Tract-specific damage at spinal cord level in pure hereditary spastic paraplegia type 4: a diffusion tensor imaging study

BACKGROUND

SPG4 is a subtype of hereditary spastic paraplegia (HSP), an upper motor neuron disorder characterized by axonal degeneration of the corticospinal tracts and the fasciculus gracilis. The few neuroimaging studies that have focused on the spinal cord in HSP are based mainly on the analysis of structural characteristics.

METHODS

We assessed diffusion-related characteristics of the spinal cord using diffusion tensor imaging (DTI), as well as structural and shape-related properties in 12 SPG4 patients and 14 controls. We used linear mixed effects models up to T3 in order to analyze the global effects of 'group' and 'clinical data' on structural and diffusion data. For DTI, we carried out a region of interest (ROI) analysis in native space for the whole spinal cord, the anterior and lateral funiculi, and the dorsal columns. We also performed a voxelwise analysis of the spinal cord to study local diffusion-related changes.

RESULTS

A reduced cross-sectional area was observed in the cervical region of SPG4 patients, with significant anteroposterior flattening. DTI analyses revealed significantly decreased fractional anisotropy (FA) and increased radial diffusivity at all the cervical and thoracic levels, particularly in the lateral funiculi and dorsal columns. The FA changes in SPG4 patients were significantly related to disease severity, measured as the Spastic Paraplegia Rating Scale score.

CONCLUSIONS

Our results in SPG4 indicate tract-specific axonal damage at the level of the cervical and thoracic spinal cord. This finding is correlated with the degree of motor disability.

Spinal cord imaging markers and recovery of standing with epidural stimulation in individuals with clinically motor complete spinal cord injury

Spinal cord epidural stimulation (scES) is an intervention to restore motor function in those with severe spinal cord injury (SCI). Spinal cord lesion characteristics assessed via magnetic resonance imaging (MRI) may contribute to understand motor recovery. This study assessed relationships between standing ability with scES and spared spinal cord tissue characteristics at the lesion site. We hypothesized that the amount of lateral spared cord tissue would be related to independent extension in the ipsilateral lower limb. Eleven individuals with chronic, clinically motor complete SCI underwent spinal cord MRI, and were subsequently implanted with scES. Standing ability and lower limb activation patterns were assessed during an overground standing experiment with scES. This assessment occurred prior to any activity-based intervention with scES. Lesion hyperintensity was segmented from T2 axial images, and template-based analysis was used to estimate spared tissue in anterior, posterior, right, and left spinal cord regions. Regression analysis was used to assess relationships between imaging and standing outcomes. Total volume of spared tissue was related to left (p = 0.007), right (p = 0.005), and bilateral (p = 0.011) lower limb extension. Spared tissue in the left cord region was related to left lower limb extension (p = 0.019). A positive trend (p = 0.138) was also observed between right spared cord tissue and right lower limb extension. In this study, MRI measures of spared spinal cord tissue were significantly related to standing outcomes with scES. These preliminary results warrant future investigation of roles of supraspinal input and MRI-detected spared spinal cord tissue on lower limb motor responsiveness to scES.

Open-access quantitative MRI data of the spinal cord and reproducibility across participants, sites and manufacturers

In a companion paper by Cohen-Adad et al. we introduce the spine generic quantitative MRI protocol that provides valuable metrics for assessing spinal cord macrostructural and microstructural integrity. This protocol was used to acquire a single subject dataset across 19 centers and a multi-subject dataset across 42 centers (for a total of 260 participants), spanning the three main MRI manufacturers: GE, Philips and Siemens. Both datasets are publicly available via git-annex. Data were analysed using the Spinal Cord Toolbox to produce normative values as well as inter/intra-site and inter/intra-manufacturer statistics. Reproducibility for the spine generic protocol was high across sites and manufacturers, with an average inter-site coefficient of variation of less than 5% for all the metrics. Full documentation and results can be found at https://spine-generic.rtfd.io/ . The datasets and analysis pipeline will help pave the way towards accessible and reproducible quantitative MRI in the spinal cord.

Spinal Cord Resting State Activity in Individuals With Fibromyalgia Who Take Opioids

Chronic pain coincides with myriad functional alterations throughout the brain and spinal cord. While spinal cord mechanisms of chronic pain have been extensively characterized in animal models and in vitro, to date, research in patients with chronic pain has focused only very minimally on the spinal cord. Previously, spinal cord functional magnetic resonance imaging (fMRI) identified regional alterations in spinal cord activity in patients (who were not taking opioids) with fibromyalgia, a chronic pain condition. Here, in patients with fibromyalgia who take opioids (N = 15), we compared spinal cord resting-state fMRI data vs. patients with fibromyalgia not taking opioids (N = 15) and healthy controls (N = 14). We hypothesized that the opioid (vs. non-opioid) patient group would show greater regional alterations in spinal cord activity (i.e., the amplitude of low frequency fluctuations or ALFF, a measure of regional spinal cord activity). However, we found that regional spinal cord activity in the opioid group was more similar to healthy controls, while regional spinal cord activity in the non-opioid group showed more pronounced differences (i.e., ventral increases and dorsal decreases in regional ALFF) vs. healthy controls. Across patient groups, self-reported fatigue correlated with regional differences in spinal cord activity. Additionally, spinal cord functional connectivity and graph metrics did not differ among groups. Our findings suggest that, contrary to our main hypothesis, patients with fibromyalgia who take opioids do not have greater alterations in regional spinal cord activity. Thus, regional spinal cord activity may be less imbalanced in patients taking opioids compared to patients not taking opioids.

Structural and resting state functional connectivity beyond the cortex

Mapping the structural and functional connectivity of the central nervous system has become a key area within neuroimaging research. While detailed network structures across the entire brain have been probed using animal models, non-invasive neuroimaging in humans has thus far been dominated by cortical investigations. Beyond the cortex, subcortical nuclei have traditionally been less accessible due to their smaller size and greater distance from radio frequency coils. However, major neuroimaging developments now provide improved signal and the resolution required to study these structures. Here, we present an overview of the connectivity between the amygdala, brainstem, cerebellum, spinal cord and the rest of the brain. While limitations to their imaging and analyses remain, we also provide some recommendations and considerations for mapping brain connectivity beyond the cortex.

Changes in White Matter of the Cervical Spinal Cord after a Single Season of Collegiate Football

The involvement of the central nervous system (CNS), specifically the white matter tracts in the cervical spinal cord, was examined with diffusion tensor imaging (DTI) following exposure to repetitive head acceleration events (HAEs) after a single season of collegiate football. Fifteen National Collegiate Athletic Association (NCAA) Division 1 football players underwent DTI of the cervical spinal cord (vertebral level C1–4) at pre-season (before any contact practices began) and post-season (within 1 week of the last regular season game) intervals. Helmet accelerometer data were also collected in parallel throughout the season. From pre-season to post-season, a significant decrease (p < 0.05) of axial diffusivity was seen within the right spino-olivary tract. In addition, a significant decrease (p < 0.05) in global white matter fractional anisotropy (FA) along with increases (p < 0.05) in global white matter mean diffusivity (MD) and radial diffusivity (RD) were found. These changes in FA from pre-season to post-season were significantly moderated by previous concussion history (p < 0.05) and number of HAEs over 80 g (p < 0.05). Despite the absence of sports-related concussion (SRC), we present measurable changes in the white matter integrity of the cervical spinal cord suggesting injury from repetitive HAEs, or SRC, may include the entirety of the CNS, not just the brain.

T1 Mapping for Microstructural Assessment of the Cervical Spinal Cord in the Evaluation of Patients with Degenerative Cervical Myelopathy

BACKGROUND AND PURPOSE

Although current radiologic evaluation of degenerative cervical myelopathy by conventional MR imaging accurately demonstrates spondylosis or degenerative disc disease causing spinal cord dysfunction, conventional MR imaging still fails to provide satisfactory anatomic and clinical correlations. In this context, we assessed the potential value of quantitative cervical spinal cord T1 mapping regarding the evaluation of patients with degenerative cervical myelopathy.

MATERIALS AND METHODS

Twenty patients diagnosed with mild and moderate-to-severe degenerative cervical myelopathy and 10 healthy subjects were enrolled in a multiparametric MR imaging protocol. Cervical spinal cord T1 mapping was performed with the MP2RAGE sequence procedure. Retrieved data were processed and analyzed regarding the global spinal cord and white and anterior gray matter on the basis of the clinical severity and the spinal canal stenosis grading.

RESULTS

Noncompressed levels in healthy controls demonstrated significantly lower T1 values than noncompressed, mild, moderate, and severe stenotic levels in patients. Concerning the entire spinal cord T1 mapping, patients with moderate-to-severe degenerative cervical myelopathy had higher T1 values compared with healthy controls. Regarding the specific levels, patients with moderate-to-severe degenerative cervical myelopathy demonstrated a T1 value increase at C1, C7, and the level of maximal compression compared with healthy controls. Patients with mild degenerative cervical myelopathy had lower T1 values than those with moderate-to-severe degenerative cervical myelopathy at the level of maximal compression. Analyses of white and anterior gray matter confirmed similar results. Strong negative correlations between individual modified Japanese Orthopaedic Association scores and T1 values were also observed.

CONCLUSIONS

In this preliminary study, 3D-MP2RAGE T1 mapping demonstrated increased T1 values in the pathology tissue samples, with diffuse medullary alterations in all patients with degenerative cervical myelopathy, especially relevant at C1 (nonstenotic level) and at the maximal compression level. Encouraging correlations observed with the modified Japanese Orthopaedic Association score make this novel approach a potential quantitative biomarker related to clinical severity in degenerative cervical myelopathy. Nevertheless, patients with mild degenerative cervical myelopathy demonstrated nonsignificant results compared with healthy controls and should now be studied in multicenter studies with larger patient populations.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

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

BACKGROUND AND PURPOSE

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSIONS

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

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

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

Cortico-spinal imaging to study pain

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

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

Study design

Retrospective.

Objectives

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

Setting

University/hospital.

Methods

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

Results

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

Conclusions

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

Sponsorship

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

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

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

Unsuspected Involvement of Spinal Cord in Alzheimer Disease

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

Cervical Cord Neurodegeneration in Traumatic and Non-Traumatic Spinal Cord Injury

This study aimed to compare macrostructural and microstructural neurodegenerative changes remote from a cervical spinal cord injury in traumatic spinal cord injury (tSCI) and degenerative cervical myelopathy (DCM) patients using quantitative magnetic resonance imaging (MRI). Twenty-nine tSCI patients, 20 mild/moderate DCM patients, and 22 healthy controls underwent a high-resolution MRI protocol at the cervical cord (C2/C3). High-resolution T2*-weighted and diffusion-weighted scans provided data to calculate tissue-specific cross-sectional areas of the spinal cord and tract-specific diffusion indices of cord white matter, respectively. Regression analysis determined associations between neurodegeneration and clinical impairment. tSCI patients showed more impairment in upper limb strength and manual dexterity when compared with DCM patients. While macrostructural MRI measures revealed a similar extent of remote cord atrophy at cervical level, microstructural measures (diffusion indices) were able to distinguish more pronounced tract-specific neurodegeneration in tSCI patients when compared with DCM patients. Tract-specific neurodegeneration was associated with upper limb impairment. Despite clinical differences between severely impaired tSCI compared with mildly affected DCM patient, extensive cord atrophy is present remotely from the focal spinal cord injury. Diffusion indices revealed greater tract-specific alterations in tSCI patients. Therefore, diffusion indices are more sensitive than macrostructural MRI measures as these are able to distinguish between traumatic and non-traumatic spinal cord injury. Neuroimaging biomarkers of cervical cord integrity hold potential as predictors of recovery and might be suitable biomarkers for interventional trials both in traumatic and non-traumatic SCI.

Regional T1 mapping of the whole cervical spinal cord using an optimized MP2RAGE sequence

The recently-proposed MP2RAGE sequence was purposely optimized for cervical spinal cord imaging at 3T. Sequence parameters were chosen to optimize gray/white matter T₁ contrast with sub-millimetric resolution and scan-time < 10 min while preserving reliable T₁ determination with minimal B₁⁺ variation effects within a range of values compatible with pathologies and surrounding structures. Results showed good agreements with IR-based measurements, high MP2RAGE-based T₁ reproducibility and preliminary evidences of age- and tract-related T₁ variations in the healthy spinal cord.

What Has Been Learned from Magnetic Resonance Imaging Examination of the Injured Human Spinal Cord: A Canadian Perspective

Magnetic resonance imaging (MRI) has transformed the way surgeons and researchers study and treat spinal cord injury. In this narrative review, we explore the historical context of imaging the human spinal cord and describe how MRI has evolved from providing the first visualization of the human spinal cord in the 1980s to a remarkable set of imaging tools today. The article focuses in particular on the role of Canadian researchers to this field. We begin by outlining the clinical context of traumatic injury to the human spinal cord and describe why current MRI standards fall short when it comes to treating this disabling condition. Parts 2 and 3 of this work explore an exciting and dramatic shift in the use of MRI technology to aid in our understanding and treatment of traumatic injury to the spinal cord. We explore the use of functional imaging (part 2) and structural imaging (part 3) and explore how these techniques have evolved, how they are used, and the challenges that we face for continued refinement and application to patients who live with the neurological and functional deficits caused by injury to the delicate spinal cord.

Multiple sclerosis lesions affect intrinsic functional connectivity of the spinal cord

Patients with multiple sclerosis present with focal lesions throughout the spinal cord. There is a clinical need for non-invasive measurements of spinal cord activity and functional organization in multiple sclerosis, given the cord's critical role in the disease. Recent reports of spontaneous blood oxygenation level-dependent fluctuations in the spinal cord using functional MRI suggest that, like the brain, cord activity at rest is organized into distinct, synchronized functional networks among grey matter regions, likely related to motor and sensory systems. Previous studies looking at stimulus-evoked activity in the spinal cord of patients with multiple sclerosis have demonstrated increased levels of activation as well as a more bilateral distribution of activity compared to controls. Functional connectivity studies of brain networks in multiple sclerosis have revealed widespread alterations, which may take on a dynamic trajectory over the course of the disease, with compensatory increases in connectivity followed by decreases associated with structural damage. We build upon this literature by examining functional connectivity in the spinal cord of patients with multiple sclerosis. Using ultra-high field 7 T imaging along with processing strategies for robust spinal cord functional MRI and lesion identification, the present study assessed functional connectivity within cervical cord grey matter of patients with relapsing-remitting multiple sclerosis (n = 22) compared to a large sample of healthy controls (n = 56). Patient anatomical images were rated for lesions by three independent raters, with consensus ratings revealing 19 of 22 patients presented with lesions somewhere in the imaged volume. Linear mixed models were used to assess effects of lesion location on functional connectivity. Analysis in control subjects demonstrated a robust pattern of connectivity among ventral grey matter regions as well as a distinct network among dorsal regions. A gender effect was also observed in controls whereby females demonstrated higher ventral network connectivity. Wilcoxon rank-sum tests detected no differences in average connectivity or power of low frequency fluctuations in patients compared to controls. The presence of lesions was, however, associated with local alterations in connectivity with differential effects depending on columnar location. The patient results suggest that spinal cord functional networks are generally intact in relapsing-remitting multiple sclerosis but that lesions are associated with focal abnormalities in intrinsic connectivity. These findings are discussed in light of the current literature on spinal cord functional MRI and the potential neurological underpinnings.

Presymptomatic spinal cord pathology in c9orf72 mutation carriers: A longitudinal neuroimaging study

OBJECTIVE

C9orf72 hexanucleotide repeats expansions account for almost half of familial amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) cases. Recent imaging studies in asymptomatic C9orf72 carriers have demonstrated cerebral white (WM) and gray matter (GM) degeneration before the age of 40 years. The objective of this study was to characterize cervical spinal cord (SC) changes in asymptomatic C9orf72 hexanucleotide carriers.

METHODS

Seventy-two asymptomatic individuals were enrolled in a prospective study of first-degree relatives of ALS and FTD patients carrying the c9orf72 hexanucleotide expansion. Forty of them carried the pathogenic mutation (C9⁺ ). Each subject underwent quantitative cervical cord imaging. Structural GM and WM metrics and diffusivity parameters were evaluated at baseline and 18 months later. Data were analyzed in C9⁺ and C9^- subgroups, and C9⁺ subjects were further stratified by age.

RESULTS

At baseline, significant WM atrophy was detected at each cervical vertebral level in C9⁺ subjects older than 40 years without associated changes in GM and diffusion tensor imaging parameters. At 18-month follow-up, WM atrophy was accompanied by significant corticospinal tract (CST) fractional anisotropy (FA) reductions. Intriguingly, asymptomatic C9⁺ subjects older than 40 years with family history of ALS (as opposed to FTD) also exhibited significant CST FA reduction at baseline.

INTERPRETATION

Cervical SC imaging detects WM atrophy exclusively in C9⁺ subjects older than 40 years, and progressive CST FA reductions can be identified on 18-month follow-up. Cervical SC magnetic resonance imaging readily captures presymptomatic pathological changes and disease propagation in c9orf72-associated conditions. ANN NEUROL 2019;86:158-167.

Convolutional Neural Network-Based Automated Segmentation of the Spinal Cord and Contusion Injury: Deep Learning Biomarker Correlates of Motor Impairment in Acute Spinal Cord Injury

BACKGROUND AND PURPOSE

Our aim was to use 2D convolutional neural networks for automatic segmentation of the spinal cord and traumatic contusion injury from axial T2-weighted MR imaging in a cohort of patients with acute spinal cord injury.

MATERIALS AND METHODS

Forty-seven patients who underwent 3T MR imaging within 24 hours of spinal cord injury were included. We developed an image-analysis pipeline integrating 2D convolutional neural networks for whole spinal cord and intramedullary spinal cord lesion segmentation. Linear mixed modeling was used to compare test segmentation results between our spinal cord injury convolutional neural network (Brain and Spinal Cord Injury Center segmentation) and current state-of-the-art methods. Volumes of segmented lesions were then used in a linear regression analysis to determine associations with motor scores.

RESULTS

Compared with manual labeling, the average test set Dice coefficient for the Brain and Spinal Cord Injury Center segmentation model was 0.93 for spinal cord segmentation versus 0.80 for PropSeg and 0.90 for DeepSeg (both components of the Spinal Cord Toolbox). Linear mixed modeling showed a significant difference between Brain and Spinal Cord Injury Center segmentation compared with PropSeg (P < .001) and DeepSeg (P < .05). Brain and Spinal Cord Injury Center segmentation showed significantly better adaptability to damaged areas compared with PropSeg (P < .001) and DeepSeg (P < .02). The contusion injury volumes based on automated segmentation were significantly associated with motor scores at admission (P = .002) and discharge (P = .009).

CONCLUSIONS

Brain and Spinal Cord Injury Center segmentation of the spinal cord compares favorably with available segmentation tools in a population with acute spinal cord injury. Volumes of injury derived from automated lesion segmentation with Brain and Spinal Cord Injury Center segmentation correlate with measures of motor impairment in the acute phase. Targeted convolutional neural network training in acute spinal cord injury enhances algorithm performance for this patient population and provides clinically relevant metrics of cord injury.

Dynamic per slice shimming for simultaneous brain and spinal cord fMRI

PURPOSE

Simultaneous brain and spinal cord functional MRI is emerging as a new tool to study the central nervous system but is challenging. Poor B0 homogeneity and small size of the spinal cord are principal obstacles to this nascent technology. Here we extend a dynamic shimming approach, first posed by Finsterbusch, by shimming per slice for both the brain and spinal cord.

METHODS

We shim dynamically by a simple and fast optimization of linear field gradients and frequency offset separately for each slice in order to minimize off-resonance for both the brain and spinal cord. Simultaneous acquisition of brain and spinal cord fMRI is achieved with high spatial resolution in the spinal cord by means of an echo-planar RF pulse for reduced FOV. Brain slice acquisition is full FOV.

RESULTS

T2*-weighted images of brain and spinal cord are acquired with high clarity and minimal observable image artifacts. Fist-clenching fMRI experiments reveal task-consistent activation in motor cortices, cerebellum, and C6-T1 spinal segments.

CONCLUSIONS

High quality functional results are obtained for a sensory-motor task. Consistent activation in both the brain and spinal cord is observed at individual levels, not only at group level. Because reduced FOV excitation is applicable to any spinal cord section, future continuation of these methods holds great potential.

Unbiased age-specific structural brain atlases for Chinese pediatric population

In magnetic resonance (MR) imaging studies of child brain development, structural brain atlases usually serve as important references for the pediatric population, in which individual images are spatially normalized into a common or standard stereotactic space. However, the popular existing pediatric brain atlases (e.g., National Institutes of Health pediatric atlases, NIH-PD) are mostly based on MR images obtained from Caucasian populations and thus are not ideal for the characterization of the brains of Chinese children due to neuroanatomical differences related to genetic and environmental factors. Here, we use an unbiased template construction algorithm to create a set of age-specific Chinese pediatric (CHN-PD) atlases based on high-quality T1-and T2-weighted MR images from 328 cognitively normal Chinese children aged 6-12 years. The CHN-PD brain atlases include asymmetric and symmetric templates, sex-specific templates and tissue probability templates, and contain multiple age-specific templates at one-year intervals. A direct comparison of the CHN-PD and NIH-PD atlases reveals dramatic anatomical differences mainly in the bilateral frontal and parietal regions. After applying the CHN-PD and NIH-PD atlases to two independent Chinese pediatric datasets (N = 114 and N = 71), we find that the CHN-PD atlases result in significantly higher accuracy than the NIH-PD atlases in both predicting "brain age" and guiding brain tissue segmentation. These results suggest that the CHN-PD brain atlases are necessary for studies of the typical and atypical development of the Chinese pediatric population. These CHN-PD atlases have been released on the Neuroimaging Informatics Tools and Resources Clearinghouse (NITRC) website (https://www.nitrc.org/projects/chn-pd).

Ten Key Insights into the Use of Spinal Cord fMRI

A comprehensive review of the literature-to-date on functional magnetic resonance imaging (fMRI) of the spinal cord is presented. Spinal fMRI has been shown, over more than two decades of work, to be a reliable tool for detecting neural activity. We discuss 10 key points regarding the history, development, methods, and applications of spinal fMRI. Animal models have served a key purpose for the development of spinal fMRI protocols and for experimental spinal cord injury studies. Applications of spinal fMRI span from animal models across healthy and patient populations in humans using both task-based and resting-state paradigms. The literature also demonstrates clear trends in study design and acquisition methods, as the majority of studies follow a task-based, block design paradigm, and utilize variations of single-shot fast spin-echo imaging methods. We, therefore, discuss the similarities and differences of these to resting-state fMRI and gradient-echo EPI protocols. Although it is newly emerging, complex connectivity and network analysis is not only possible, but has also been shown to be reliable and reproducible in the spinal cord for both task-based and resting-state studies. Despite the technical challenges associated with spinal fMRI, this review identifies reliable solutions that have been developed to overcome these challenges.

Advancing imaging technologies for patients with spinal pain: with a focus on whiplash injury

BACKGROUND CONTEXT

Radiological observations of soft-tissue changes that may relate to clinical symptoms in patients with traumatic and non-traumatic spinal disorders are highly controversial. Studies are often of poor quality and findings are inconsistent. A plethora of evidence suggests some pathoanatomical findings from traditional imaging applications are common in asymptomatic participants across the life span, which further questions the diagnostic, prognostic, and theranostic value of traditional imaging. Although we do not dispute the limited evidence for the clinical importance of most imaging findings, we contend that the disparate findings across studies may in part be due to limitations in the approaches used in assessment and analysis of imaging findings.

PURPOSE

This clinical commentary aimed to (1) briefly detail available imaging guidelines, (2) detail research-based evidence around the clinical use of findings from advanced, but available, imaging applications (eg, fat and water magnetic resonance imaging and magnetization transfer imaging), and (3) introduce how evolving imaging technologies may improve our mechanistic understanding of pain and disability, leading to improved treatments and outcomes.

STUDY DESIGN/SETTING

A non-systematic review of the literature is carried out.

METHODS

A narrative summary (including studies from the authors' own work in whiplash injuries) of the available literature is provided.

RESULTS

An emerging body of evidence suggests that the combination of existing imaging sequences or the use of developing imaging technologies in tandem with a good clinical assessment of modifiable risk factors may provide important diagnostic information toward the exploration and development of more informed and effective treatment options for some patients with traumatic neck pain.

CONCLUSIONS

Advancing imaging technologies may help to explain the seemingly disconnected spectrum of biopsychosocial signs and symptoms of traumatic neck pain.

Thermal Stimulation Changes Diffusivity of the Spinothalamic Tract

STUDY DESIGN

An experimental study.

OBJECTIVE

This study aimed to investigate task-dependent changes in fractional anisotropy (FA) within the spinal cord during painful stimulation.

SUMMARY OF BACKGROUND DATA

Earlier experiments by Mandl et al (2008, 2013) used non-invasive functional diffusion tensor imaging (fDTI) to detect white matter fibers that were active during functional tasks. In two studies, it was observed that FA of involved white matter tracts exhibited repeatable task-related increases. In this study, we attempted to extend the fDTI work in the spinal cord.

METHODS

Twenty-three healthy, right-handed men (mean age 22 yrs, standard deviation [SD] = 4) were invited to participate in this study. Diffusion-weighted images were collected over spinal levels C2 to T4 during a painful thermal stimulus applied to the left thenar eminence. In order to investigate task-related activity, FA values within the contralateral (right) spinothalamic tract were analyzed using a generalized estimating equations (GEE) procedure. As a control, we also examined activity in the ipsilateral and contralateral corticospinal tracts, which are not considered to be involved in nociception.

RESULTS

Significant task-related decreases in FA were observed in the right spinothalamic tract at vertebral levels C2-C5 (Wald X(1) = 17.754, P < 0.001). There was no change in control regions at levels C7-T2 of the same tract, which are located below the level of input from dermatome C6, Wald X(1) = 0.185, P = 0.667. Results in all other regions assessed, that is, the left spinothalamic tract and bilateral corticospinal tract, were also not significant (P > 0.05).

CONCLUSION

The current findings suggest that task-related changes in FA associated with the transmission of pain signals along the spinal cord can be detected using fDTI. We observed decreased FA values in the contralateral (right) spinothalamic tract following painful stimulation, while no such activity was apparent in control regions.

LEVEL OF EVIDENCE

5.

Reliable and fast volumetry of the lumbar spinal cord using cord image analyser (Cordial)

OBJECTIVE

To validate the precision and accuracy of the semi-automated cord image analyser (Cordial) for lumbar spinal cord (SC) volumetry in 3D T1w MRI data of healthy controls (HC).

MATERIALS AND METHODS

40 3D T1w images of 10 HC (w/m: 6/4; age range: 18-41 years) were acquired at one 3T-scanner in two MRI sessions (time interval 14.9±6.1 days). Each subject was scanned twice per session, allowing determination of test-retest reliability both in back-to-back (intra-session) and scan-rescan images (inter-session). Cordial was applied for lumbar cord segmentation twice per image by two raters, allowing for assessment of intra- and inter-rater reliability, and compared to a manual gold standard.

RESULTS

While manually segmented volumes were larger (mean: 2028±245 mm³ vs. Cordial: 1636±300 mm³, p<0.001), accuracy assessments between manually and semi-automatically segmented images showed a mean Dice-coefficient of 0.88±0.05. Calculation of within-subject coefficients of variation (COV) demonstrated high intra-session (1.22-1.86%), inter-session (1.26-1.84%), as well as intra-rater (1.73-1.83%) reproducibility. No significant difference was shown between intra- and inter-session reproducibility or between intra-rater reliabilities. Although inter-rater reproducibility (COV: 2.87%) was slightly lower compared to all other reproducibility measures, between rater consistency was very strong (intraclass correlation coefficient: 0.974).

CONCLUSION

While under-estimating the lumbar SCV, Cordial still provides excellent inter- and intra-session reproducibility showing high potential for application in longitudinal trials.

KEY POINTS

• Lumbar spinal cord segmentation using the semi-automated cord image analyser (Cordial) is feasible. • Lumbar spinal cord is 40-mm cord segment 60 mm above conus medullaris. • Cordial provides excellent inter- and intra-session reproducibility in lumbar spinal cord region. • Cordial shows high potential for application in longitudinal trials.

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Microstructural imaging in the spinal cord and validation strategies

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

High-resolution magnetization transfer MRI in patients with cervical spondylotic myelopathy

Magnetization transfer (MT) contrast has been established as a marker of myelin integrity, and cervical spondylotic myelopathy is known to cause demyelination. Ten patients with clinical and magnetic resonance imaging (MRI) manifestations of cervical spondylotic myelopathy (CSM) were compared to the MRIs of seven historic healthy controls, using the magnetization transfer ratio (MTR) and Nurick scores as the primary metrics. Transverse slices through the intervertebral discs of the cervical spine were acquired using a gradient echo sequence (MEDIC) with and without an MT saturation pulse on a 3 Tesla Siemens Prisma scanner (TR = 300 ms, TE_eff = 17 ms, flip angle = 30°, in-plane resolution = 0.47 × 0.47 mm²). The CSM patients tended to have a lower mean MTR (30.4 ± 6.5) than the controls (34.8 ± 3.8), but the difference was not significant (independent samples t-test, p = 0.110, Cohen's d = 0.80). The mean MTR across all intervertebral disc levels was not significantly correlated to the Nurick score (Spearman's ρ = -0.489, p = 0.151). The intervertebral level with the lowest MTR in each cohort was not significantly different between groups (equal variances not assumed, t = 1.965, dof = 14.8, p = 0.068, Cohen's d = 0.88), but the CSM patients tended to have a lower MTR. The mean MTR at this level was negatively correlated to the Nurick score among CSM patients (Spearman's ρ = -0.725, p = 0.018). CSM patients tended to have decreased MTR indicating myelin degradation compared to our healthy subjects, and MTR was negatively correlated with the severity of CSM.

Spinal cord MRI at 7T

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

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

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

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

Spinal Cord Gray Matter Atrophy in Amyotrophic Lateral Sclerosis

BACKGROUND AND PURPOSE

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSIONS

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

Thermal Stimulation Alters Cervical Spinal Cord Functional Connectivity in Humans

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

Generative diffeomorphic modelling of large MRI data sets for probabilistic template construction

In this paper we present a hierarchical generative model of medical image data, which can capture simultaneously the variability of both signal intensity and anatomical shapes across large populations. Such a model has a direct application for learning average-shaped probabilistic tissue templates in a fully automated manner. While in principle the generality of the proposed Bayesian approach makes it suitable to address a wide range of medical image computing problems, our work focuses primarily on neuroimaging applications. In particular we validate the proposed method on both real and synthetic brain MR scans including the cervical cord and demonstrate that it yields accurate alignment of brain and spinal cord structures, as compared to state-of-the-art tools for medical image registration. At the same time we illustrate how the resulting tissue probability maps can readily be used to segment, bias correct and spatially normalise unseen data, which are all crucial pre-processing steps for MR imaging studies.

•

We present a generative modelling framework to process large MRI data sets.

•

The proposed framework can serve to learn average-shaped tissue probability maps and empirical intensity priors.

•

We explore semi-supervised learning and variational inference schemes.

•

The method is validated against state-of-the-art tools using publicly available data.

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

BACKGROUND AND PURPOSE

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSIONS

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

The efficiency of retrospective artifact correction methods in improving the statistical power of between-group differences in spinal cord DTI

Diffusion tensor imaging (DTI) is a promising approach for investigating the white matter microstructure of the spinal cord. However, it suffers from severe susceptibility, physiological, and instrumental artifacts present in the cord. Retrospective correction techniques are popular approaches to reduce these artifacts, because they are widely applicable and do not increase scan time.

In this paper, we present a novel outlier rejection approach (reliability masking) which is designed to supplement existing correction approaches by excluding irreversibly corrupted and thus unreliable data points from the DTI index maps. Then, we investigate how chains of retrospective correction techniques including (i) registration, (ii) registration and robust fitting, and (iii) registration, robust fitting, and reliability masking affect the statistical power of a previously reported finding of lower fractional anisotropy values in the posterior column and lateral corticospinal tracts in cervical spondylotic myelopathy (CSM) patients.

While established post-processing steps had small effect on the statistical power of the clinical finding (slice-wise registration: −0.5%, robust fitting: +0.6%), adding reliability masking to the post-processing chain increased it by 4.7%. Interestingly, reliability masking and registration affected the t-score metric differently: while the gain in statistical power due to reliability masking was mainly driven by decreased variability in both groups, registration slightly increased variability. In conclusion, reliability masking is particularly attractive for neuroscience and clinical research studies, as it increases statistical power by reducing group variability and thus provides a cost-efficient alternative to increasing the group size.

•

A novel outlier rejection technique (reliability masking) is introduced.

•

Standard artifact correction has little effect on the statistical power of between-group differences.

•

Reliability masking improves the statistical power of between-group differences.

•

This improvement is driven by decreased group-level variability.