The current state-of-the-art of spinal cord imaging: applications

Imaging protocols for non-traumatic spinal cord injury: current state of the art and future directions

INTRODUCTION

Non-traumatic spinal cord injury (NTSCI) is a term used to describe damage to the spinal cord from sources other than trauma. Neuroimaging techniques such as computerized tomography (CT) and magnetic resonance imaging (MRI) have improved our ability to diagnose and manage NTSCIs. Several practice guidelines utilize MRI in the diagnostic evaluation of traumatic and non-traumatic SCI to direct surgical intervention.

AREAS COVERED

The authors review practices surrounding the imaging of various causes of NTSCI as well as recent advances and future directions for the use of novel imaging modalities in this realm. The authors also present discussions around the use of simple radiographs and advanced MRI modalities in clinical settings, and briefly highlight areas of active research that seek to advance our understanding and improve patient care.

EXPERT OPINION

Although several obstacles must be overcome, it appears highly likely that novel quantitative imaging features and advancements in artificial intelligence (AI) as well as machine learning (ML) will revolutionize degenerative cervical myelopathy (DCM) care by providing earlier diagnosis, accurate localization, monitoring for deterioration and neurological recovery, outcome prediction, and standardized practice. Some intriguing findings in these areas have been published, including the identification of possible serum and cerebrospinal fluid biomarkers, which are currently in the early phases of translation.

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

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

INTRODUCTION

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

METHODS

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

RESULTS

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

CONCLUSION

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

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

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

Degenerative cervical myelopathy: Where have we been? Where are we now? Where are we going?

Degenerative cervical myelopathy (DCM), a recently coined term, encompasses a group of age-related and genetically associated pathologies that affect the cervical spine, including cervical spondylotic myelopathy and ossification of the posterior longitudinal ligament (OPLL). Given the significant contribution of DCM to global disease and disability, there are worldwide efforts to promote research and innovation in this area. An AO Spine effort termed 'RECODE-DCM' was initiated to create an international multistakeholder consensus group, involving patients, caregivers, physicians and researchers, to focus on launching actionable discourse on DCM. In order to improve the management, treatment and results for DCM, the RECODE-DCM consensus group recently identified ten priority areas for translational research. The current article summarizes recent advancements in the field of DCM. We first discuss the comprehensive definition recently refined by the RECODE-DCM group, including steps taken to arrive at this definition and the supporting rationale. We then provide an overview of the recent advancements in our understanding of the pathophysiology of DCM and modalities to clinically assess and diagnose DCM. A focus will be set on advanced imaging techniques that may offer the opportunity to improve characterization and diagnosis of DCM. A summary of treatment modalities, including surgical and nonoperative options, is then provided along with future neuroprotective and neuroregenerative strategies. This review concludes with final remarks pertaining to the genetics involved in DCM and the opportunity to leverage this knowledge toward a personalized medicine approach.

Spinal neurovascular coupling is preserved despite time-dependent alterations of spinal cord blood flow responses in a rat model of chronic back pain: implications for functional spinal cord imaging

ABSTRACT

Functional magnetic resonance imaging has been used to investigate nociceptive processes in patients with chronic pain. However, the results may be confounded with changes in neurovascular coupling induced by chronic pain. The objective of this study was to examine spinal neurovascular coupling in a rat model of chronic back pain induced by muscle inflammation. Rats received 150 µL intramuscular injections of either complete Freund adjuvant (CFA: n = 18) or saline (control [CTL]: n = 18) in L5-L6 paravertebral muscles. Under 1.2% isoflurane anesthesia, spinal cord blood flow (SCBF) and local field potentials evoked by electrical stimulation of the sciatic nerve were recorded simultaneously in the lumbar enlargement of the spinal cord, 14 or 28 days after the injections. Mechanical hypersensitivity was observed in CFA rats compared with CTL rats for the back ( P < 0.001) and hind paws ( P < 0.01). Spinal cord blood flow response amplitude and local field potential amplitude were not significantly different between groups (day 14: P > 0.5; day 28: P > 0.6). However, the time course of SCBF responses was different between groups on day 14 ( P < 0.001) and day 28 ( P < 0.001). Nevertheless, neurovascular coupling was comparable between groups on days 14 and 28, whether neurovascular coupling was calculated with the amplitude or the area under the curve of SCBF responses (all P > 0.2). These results indicate that spinal hemodynamic changes reflect neuronal activity in this animal model, although the time course of SCBF responses is affected by chronic inflammatory back pain. This warrants a careful use of spinal functional magnetic resonance imaging in animal models and patients with chronic back pain.

Advances in The Application and Research of Magnetic Resonance Diffusion Kurtosis Imaging in The Musculoskeletal System

Magnetic resonance diffusion kurtosis imaging (DKI) is an emerging magnetic resonance imaging (MRI) technique that can reflect microstructural changes in tissue through non-Gaussian diffusion of water molecules. Compared to traditional diffusion weighted imaging (DWI), the DKI model has shown greater sensitivity for diagnosis of musculoskeletal diseases and can help formulate more reasonable treatment plans. Moreover, DKI is an important auxiliary examination for evaluation of the motor function of the musculoskeletal system. This article briefly introduces the basic principles of DKI and reviews the application and research of DKI in the evaluation of disorders of the musculoskeletal system (including bone tumors, soft tissue tumors, spinal lesions, chronic musculoskeletal diseases, musculoskeletal trauma, and developmental disorders) as well as the normal musculoskeletal tissues. EVIDENCE LEVEL: 5 TECHNICAL EFFICACY: 1.

Fully Automatic Method for Reliable Spinal Cord Compartment Segmentation in Multiple Sclerosis

BACKGROUND AND PURPOSE

Fully automatic quantification methods of spinal cord compartments are needed to study pathologic changes of the spinal cord GM and WM in MS in vivo. We propose a novel method for automatic spinal cord compartment segmentation (SCORE) in patients with MS.

MATERIALS AND METHODS

The cervical spinal cords of 24 patients with MS and 24 sex- and age-matched healthy controls were scanned on a 3T MR imaging system, including an averaged magnetization inversion recovery acquisition sequence. Three experienced raters manually segmented the spinal cord GM and WM, anterior and posterior horns, gray commissure, and MS lesions. Subsequently, manual segmentations were used to train neural segmentation networks of spinal cord compartments with multidimensional gated recurrent units in a 3-fold cross-validation fashion. Total intracranial volumes were quantified using FreeSurfer.

RESULTS

The intra- and intersession reproducibility of SCORE was high in all spinal cord compartments (eg, mean relative SD of GM and WM: ≤ 3.50% and ≤1.47%, respectively) and was better than manual segmentations (all P < .001). The accuracy of SCORE compared with manual segmentations was excellent, both in healthy controls and in patients with MS (Dice similarity coefficients of GM and WM: ≥ 0.84 and ≥0.92, respectively). Patients with MS had lower total WM areas (P < .05), and total anterior horn areas (P < .01 respectively), as measured with SCORE.

CONCLUSIONS

We demonstrate a novel, reliable quantification method for spinal cord tissue segmentation in healthy controls and patients with MS and other neurologic disorders affecting the spinal cord. Patients with MS have reduced areas in specific spinal cord tissue compartments, which may be used as MS biomarkers.

Effects of variability in manually contoured spinal cord masks on fMRI co-registration and interpretation

Functional magnetic resonance imaging (fMRI) of the human spinal cord (SC) is a unique non-invasive method for characterizing neurovascular responses to stimuli. Group-analysis of SC fMRI data involves co-registration of subject-level data to standard space, which requires manual masking of the cord and may result in bias of group-level SC fMRI results. To test this, we examined variability in SC masks drawn in fMRI data from 21 healthy participants from a completed study mapping responses to sensory stimuli of the C7 dermatome. Masks were drawn on temporal mean functional image by eight raters with varying levels of neuroimaging experience, and the rater from the original study acted as a reference. Spatial agreement between rater and reference masks was measured using the Dice Similarity Coefficient, and the influence of rater and dataset was examined using ANOVA. Each rater's masks were used to register functional data to the PAM50 template. Gray matter-white matter signal contrast of registered functional data was used to evaluate the spatial normalization accuracy across raters. Subject- and group-level analyses of activation during left- and right-sided sensory stimuli were performed for each rater's co-registered data. Agreement with the reference SC mask was associated with both rater (F(7, 140) = 32.12, P < 2 × 10-16, η2 = 0.29) and dataset (F(20, 140) = 20.58, P < 2 × 10-16, η2 = 0.53). Dataset variations may reflect image quality metrics: the ratio between the signal intensity of spinal cord voxels and surrounding cerebrospinal fluid was correlated with DSC results (p < 0.001). As predicted, variability in the manually-drawn masks influenced spatial normalization, and GM:WM contrast in the registered data showed significant effects of rater and dataset (rater: F(8, 160) = 23.57, P < 2 × 10-16, η2 = 0.24; dataset: F(20, 160) = 22.00, P < 2 × 10-16, η2 = 0.56). Registration differences propagated into subject-level activation maps which showed rater-dependent agreement with the reference. Although group-level activation maps differed between raters, no systematic bias was identified. Increasing consistency in manual contouring of spinal cord fMRI data improved co-registration and inter-rater agreement in activation mapping, however our results suggest that improvements in image acquisition and post-processing are also critical to address.

Automated slice-specific z-shimming for functional magnetic resonance imaging of the human spinal cord

Functional magnetic resonance imaging (fMRI) of the human spinal cord faces many challenges, such as signal loss due to local magnetic field inhomogeneities. This issue can be addressed with slice-specific z-shimming, which compensates for the dephasing effect of the inhomogeneities using a slice-specific gradient pulse. Here, we aim to address outstanding issues regarding this technique by evaluating its effects on several aspects that are directly relevant for spinal fMRI and by developing two automated procedures in order to improve upon the time-consuming and subjective nature of manual selection of z-shims: one procedure finds the z-shim that maximizes signal intensity in each slice of an EPI reference-scan and the other finds the through-slice field inhomogeneity for each EPI-slice in field map data and calculates the required compensation gradient moment. We demonstrate that the beneficial effects of z-shimming are apparent across different echo times, hold true for both the dorsal and ventral horn, and are also apparent in the temporal signal-to-noise ratio (tSNR) of EPI time-series data. Both of our automated approaches were faster than the manual approach, lead to significant improvements in gray matter tSNR compared to no z-shimming and resulted in beneficial effects that were stable across time. While the field-map-based approach performed slightly worse than the manual approach, the EPI-based approach performed as well as the manual one and was furthermore validated on an external corticospinal data-set (N > 100). Together, automated z-shimming may improve the data quality of future spinal fMRI studies and lead to increased reproducibility in longitudinal studies.

Imaging of Common Spinal Cord Diseases

Spinal cord evaluation is an integral part of spine assessment, and its reliable imaging work-up is mandatory because even localized lesions may produce serious effects with potentially irreversible sequelae. Spinal cord alterations are found both incidentally during spine evaluation in otherwise neurologically asymptomatic patients or during neurologic/neuroradiologic assessment in myelopathic patients. Myelopathy (an umbrella term for any neurologic deficit that refers to spinal cord impairment) can be caused by intrinsic lesions or extrinsic mechanical compression, and its etiology may be both traumatic and/or nontraumatic. The symptoms largely depend on the size/extension of lesions, ranging from incontinence to ataxia, from spasticity to hyperreflexia, from numbness to weakness. Magnetic resonance imaging is the reference imaging modality in spinal cord evaluation, ensuring the best signal and spatial resolution. We provide an overview of the most common spinal cord disorders encountered by radiologists and describe the technical measures that offer optimal spinal cord visualization.

Imaging and Electrophysiology for Degenerative Cervical Myelopathy [AO Spine RECODE-DCM Research Priority Number 9]

STUDY DESIGN

Narrative review.

OBJECTIVE

The current review aimed to describe the role of existing techniques and emerging methods of imaging and electrophysiology for the management of degenerative cervical myelopathy (DCM), a common and often progressive condition that causes spinal cord dysfunction and significant morbidity globally.

METHODS

A narrative review was conducted to summarize the existing literature and highlight future directions.

RESULTS

Anatomical magnetic resonance imaging (MRI) is well established in the literature as the key imaging tool to identify spinal cord compression, disc herniation/bulging, and inbuckling of the ligamentum flavum, thus facilitating surgical planning, while radiographs and computed tomography (CT) provide complimentary information. Electrophysiology techniques are primarily used to rule out competing diagnoses. However, signal change and measures of cord compression on conventional MRI have limited utility to characterize the degree of tissue injury, which may be helpful for diagnosis, prognostication, and repeated assessments to identify deterioration. Early translational studies of quantitative imaging and electrophysiology techniques show potential of these methods to more accurately reflect changes in spinal cord microstructure and function.

CONCLUSION

Currently, clinical management of DCM relies heavily on anatomical MRI, with additional contributions from radiographs, CT, and electrophysiology. Novel quantitative assessments of microstructure, perfusion, and function have the potential to transform clinical practice, but require robust validation, automation, and standardization prior to uptake.

Extent of cord pathology in the lumbosacral enlargement in non-traumatic versus traumatic spinal cord injury

This study compares remote neurodegenerative changes caudal to a cervical injury in degenerative cervical myelopathy (DCM) (i.e., non-traumatic) and incomplete traumatic spinal cord injury (tSCI) patients, using MRI-based tissue area measurements and diffusion tensor imaging (DTI). Eighteen mild to moderate DCM patients with sensory impairments (mJOA score: 16.2±1.9), 14 incomplete tetraplegic tSCI patients (AIS C&D), and 20 healthy controls were recruited. All participants received DTI and T2*-weighted scans in the lumbosacral enlargement (caudal to injury) and at C2/C3 (rostral to injury). MRI readouts included DTI metrics in the white matter (WM) columns and cross-sectional WM and gray matter area. One-way ANOVA with Tukey's post-hoc comparison (p<0.05) was used to assess group differences. In the lumbosacral enlargement, compared to DCM, tSCI patients exhibited decreased fractional anisotropy in the lateral (tSCI vs. DCM, -11.9%, p=0.007) and ventral WM column (-8.0%, p=0.021), and showed trend toward lower values in the dorsal column (-8.9%, p=0.068). At C2/C3, compared to controls, fractional anisotropy was lower in both groups in the dorsal (DCM vs. controls, -7.9%, p=0.024; tSCI vs. controls, -10.0%, p=0.007) and in the lateral column (DCM: -6.2%, p=0.039; tSCI: -13.3%, p<0.001), while tSCI patients had lower fractional anisotropy than DCM patients in the lateral column (-7.6%, p=0.029). WM areas were not different between patient groups but were lower compared to controls in the lumbosacral enlargement (DCM: -16.9%, p<0.001; tSCI, -10.5%, p=0.043) and at C2/C3 (DCM: -16.0%, p<0.001; tSCI: -18.1%, p<0.001). In conclusion, mild to moderate DCM and incomplete tSCI lead to similar degree of degeneration of the dorsal and lateral columns at C2/C3, but tSCI results in more widespread white matter damage in the lumbosacral enlargement. These remote changes are likely to contribute to the patients' impairment and recovery. DTI is a sensitive tool to assess remote pathological changes in DCM and tSCI patients.

Role of MRI evaluation in acute secondary inability to walk in children

Background

Acute bilateral lower limb weakness is a common problem in children which necessitates a rapid method for diagnosis. MRI is a non-invasive imaging technique that produces high-quality images of the internal structure of the brain and spinal cord.

Results

MRI was very helpful in reaching rapid and prompt diagnosis in children with acute inability to walk. Acute disseminated encephalomyelitis (ADEM), Guillain–Barré syndrome (GBS), and acute transverse myelitis (ATM) were the most common causes in our study. MRI proved to be of high sensitivity in detecting the lesions and reaching the diagnosis in ADEM and GBS; however, there was no significant relation between the lesions’ size, enhancement pattern, and severity of the disease or prognosis, yet in ATM the site of the lesion and number of cord segment affection were significantly related to the severity of the disease and prognosis.

Conclusion

MRI is a quick tool to reach the diagnosis of children with acute secondary inability to walk, and to eliminate other differential diagnosis which is essential for proper treatment and rapid full recovery. It is highly sensitive in detecting the lesions, their site and size.

Generic acquisition protocol for quantitative MRI of the spinal cord

Quantitative spinal cord (SC) magnetic resonance imaging (MRI) presents many challenges, including a lack of standardized imaging protocols. Here we present a prospectively harmonized quantitative MRI protocol, which we refer to as the spine generic protocol, for users of 3T MRI systems from the three main manufacturers: GE, Philips and Siemens. The protocol provides guidance for assessing SC macrostructural and microstructural integrity: T1-weighted and T2-weighted imaging for SC cross-sectional area computation, multi-echo gradient echo for gray matter cross-sectional area, and magnetization transfer and diffusion weighted imaging for assessing white matter microstructure. In a companion paper from the same authors, the spine generic protocol was used to acquire data across 42 centers in 260 healthy subjects. The key details of the spine generic protocol are also available in an open-access document that can be found at https://github.com/spine-generic/protocols . The protocol will serve as a starting point for researchers and clinicians implementing new SC imaging initiatives so that, in the future, inclusion of the SC in neuroimaging protocols will be more common. The protocol could be implemented by any trained MR technician or by a researcher/clinician familiar with MRI acquisition.

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.

A pipeline to quantify spinal cord atrophy with deep learning: Application to differentiation of MS and NMOSD patients

PURPOSE

Quantitative measurement of various anatomical regions of the brain and spinal cord (SC) in MRI images are used as unique biomarkers to consider progress and effects of demyelinating diseases of the central nervous system. This paper presents a fully-automated image processing pipeline which quantifies the SC volume of MRI images.

METHODS

In the proposed pipeline, after conducting some pre-processing tasks, a deep convolutional network is utilized to segment the spinal cord cross-sectional area (SCCSA) of each slice. After full segmentation, certain extra slices interpolate between each two adjacent slices using the shape-based interpolation method. Then, a 3D model of the SC is reconstructed, and, by counting the voxels of it, the SC volume is calculated. The performance of the proposed method for the SCCSA segmentation is evaluated on 140 MRI images. Subsequently, to demonstrate the application of the proposed pipeline, we study the differentiations of SC atrophy between 38 Multiple Sclerosis (MS) and 25 Neuromyelitis Optica Spectrum Disorder (NMOSD) patients.

RESULTS

The experimental results of the SCCSA segmentation indicate that the proposed method, adapted by Mask R-CNN, presented the most satisfactory result with the average Dice coefficient of 0.96. For this method, statistical metrics including sensitivity, specificity, accuracy, and precision are 97.51%, 99.98%, 99.92%, and 98.04% respectively. Moreover, the t-test result (p-value = 0.00089) verified a significant difference between the SC atrophy of MS and NMOSD patients.

CONCLUSION

The pipeline efficiently quantifies the SC volume of MRI images and can be utilized as an affordable computer-aided tool for diagnostic purposes.

Contribution of astrocytes to neurovascular coupling in the spinal cord of the rat

Functional magnetic resonance imaging (fMRI) of the spinal cord relies on the integrity of neurovascular coupling (NVC) to infer neuronal activity from hemodynamic changes. Astrocytes are a key component of cerebral NVC, but their role in spinal NVC is unclear. The objective of this study was to examine whether inhibition of astrocyte metabolism by fluorocitrate alters spinal NVC. In 14 rats, local field potential (LFP) and spinal cord blood flow (SCBF) were recorded simultaneously in the lumbosacral enlargement during noxious stimulation of the sciatic nerve before and after a local administration of fluorocitrate (N  = 7) or saline (N  = 7). Fluorocitrate significantly reduced SCBF responses (p  < 0.001) but not LFP amplitude (p  = 0.22) compared with saline. Accordingly, NVC was altered by fluorocitrate compared with saline (p  < 0.01). These results support the role of astrocytes in spinal NVC and have implications for spinal cord imaging with fMRI for conditions in which astrocyte metabolism may be altered.

Imaging Considerations in Spinal Cord Evaluation

Multiple diverse pathologies result in the clinical presentation of myelopathy. The preferred way to image the spinal cord depends on clinical history, anatomic site of interest, and patient issues limiting certain imaging modalities. This radiology-focused article discusses pertinent physiological considerations, reviews basic and newer imaging techniques, and examines several distinct disease entities in order to highlight the key role of imaging in the work-up of myelopathy.

Spinal cord fMRI with MB-SWIFT for assessing epidural spinal cord stimulation in rats

Purpose

Electrical epidural spinal cord stimulation (SCS) is used as a treatment for chronic pain as well as to partially restore motor function after a spinal cord injury. Monitoring the spinal cord activity during SCS with fMRI could provide important and objective measures of integrative responses to treatment. Unfortunately, spinal cord fMRI is severely challenged by motion and susceptibility artifacts induced by the implanted electrode and bones. This pilot study introduces multi‐band sweep imaging with Fourier transformation (MB‐SWIFT) technique for spinal cord fMRI during SCS in rats. Given the close to zero acquisition delay and high bandwidth in 3 dimensions, MB‐SWIFT is demonstrated to be highly tolerant to motion and susceptibility‐induced artifacts and thus holds promise for fMRI during SCS.

Methods

MB‐SWIFT with 0.78 × 0.78 × 1.50 mm³ spatial resolution and 3‐s temporal resolution was used at 9.4 Tesla in rats undergoing epidural SCS at different frequencies. Its performance was compared with spin echo EPI. The origin of the functional contrast was also explored using suppression bands.

Results

MB‐SWIFT was tolerant to electrode‐induced artifacts and respiratory motion, leading to substantially higher fMRI sensitivity than spin echo fMRI. Clear stimulation frequency‐dependent responses to SCS were detected in the rat spinal cord close to the stimulation site. The origin of MB‐SWIFT fMRI signals was consistent with dominant inflow effects.

Conclusion

fMRI of the rat spinal cord during SCS can be consistently achieved with MB‐SWIFT, thus providing a valuable experimental framework for assessing the effects of SCS on the central nervous system.

Sensitivity and specificity of CEST and NOE MRI in injured spinal cord in monkeys

•

Compare sensitivity and specificity of CEST and NOE measures from 6-pool fitting.

•

Differentiate regional molecular signatures at and around spinal cord injury.

•

Provide parameters that improve the diagnostic accuracy of molecular alteration.

•

Down-sampled data acquisition can capture the characteristic molecular profile.

•

High translational potential for clinical assessment of spinal cord injury.

Purpose

The sensitivity and accuracy of chemical exchange saturation transfer (CEST) and nuclear Overhauser enhancement (NOE) effects for assessing injury-associated changes in cervical spinal cords were evaluated in squirrel monkeys. Multiple interacting pools of protons, including one identified by an NOE at −1.6 ppm relative to water (NOE(-1.6)), were derived and quantified from fitting proton Z-spectra. The effects of down-sampled data acquisitions and corrections for non-specific factors including T₁, semi-solid magnetization transfer, and direct saturation of free water (DS), were investigated. The overall goal is to develop a protocol for rapid data acquisition for assessing the molecular signatures of the injured spinal cord and its surrounding regions.

Methods

MRI scans were recorded of anesthetized squirrel monkeys at 9.4 T, before and after a unilateral dorsal column sectioning of the cervical spinal cord. Z-spectral images at 51 different RF offsets were acquired. The amplitudes of CEST and NOE effects from multiple proton pools were quantified using a six-pool Lorenzian fitting of each Z-spectrum (MTR_mfit). In addition, down-sampled data using reduced selections of RF offsets were analyzed and compared. An apparent exchange-dependent relaxation (AREX_mfit) method was also used to correct for non-specific factors in quantifying regional spectra around lesion sites.

Results

The parametric maps from multi-pool fitting using the complete sampling data (P51e) detected unilateral changes at and around the injury. The maps derived from selected twofold down-sampled data with appropriate interpolation (P26sI51) revealed quite similar spatial distributions of different pools as those obtained using P51e at each resonance shift. Across 10 subjects, both data acquisition schemes detected significant decreases in NOE(-3.5) and NOE(-1.6) and increases in DS(0.0) and CEST(3.5) at the lesion site relative to measures of the normal tissues before injury. AREX_mfit of cysts and other abnormal tissues at and around the lesion site also exhibited significant changes, especially at 3.5, −1.6 and −3.5 ppm RF offsets.

Conclusion

These results confirm that a reduced set of RF offsets and down sampling are adequate for CEST imaging of injured spinal cord and allow shorter imaging times and/or permit additional signal averaging. AREX_mfit correction improved the accuracy of CEST and NOE measures. The results provide a rapid (~13 mins), sensitive, and accurate protocol for deriving multiple NOE and CEST effects simultaneously in spinal cord imaging at high field.

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.

The neuroanatomical-functional paradox in spinal cord injury

Although lesion size is widely considered to be the most reliable predictor of outcome after CNS injury, lesions of comparable size can produce vastly different magnitudes of functional impairment and subsequent recovery. This neuroanatomical-functional paradox is likely to contribute to the many failed attempts to independently replicate findings from animal models of neurotrauma. In humans, the analogous clinical-radiological paradox could explain why individuals with similar injuries can respond differently to rehabilitation. We describe the neuroanatomical-functional paradox in the context of traumatic spinal cord injury (SCI) and discuss the underlying mechanisms of the paradox, including the concepts of lesion-affected and recovery-related networks. We also consider the various secondary complications that further limit the accuracy of outcome prediction in SCI and provide suggestions for how to increase the predictive, translational value of preclinical SCI models.

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

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

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

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

Correlation of a new hydrodynamic index with other effective indexes in Chiari I malformation patients with different associations

This study aimed to find a new CSF hydrodynamic index to assess Chiari type I malformation (CM-I) patients’ conditions and examine the relationship of this new index with morphometric and volumetric changes in these patients and their clinical symptoms. To this end, 58 CM-I patients in four groups and 20 healthy subjects underwent PC-MRI. Ten morphometric and three volumetric parameters were calculated. The CSF hydrodynamic parameters were also analyzed through computational fluid dynamic (CFD) simulation. The maximum CSF pressure was identified as a new hydrodynamic parameter to assess the CM-I patients’ conditions. This parameter was similar in patients with the same symptoms regardless of the group to which they belonged. The result showed a weak correlation between the maximum CSF pressure and the morphometric parameters in the patients. Among the volumetric parameters, PCF volume had the highest correlation with the maximum CSF pressure, which its value being higher in patients with CM-I/SM/scoliosis (R² = 65.6%, P = 0.0022) than in the other patients. PCF volume was the more relevant volumetric parameter to assess the patients’ symptoms. The values of PCF volume were greater in patients that headache symptom was more obvious than other symptoms, as compared to the other patients.

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

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

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

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

Imaging of the thoracic spinal cord using radially sampled averaged magnetization inversion recovery acquisitions

BACKGROUND

Spinal cord (SC) gray and white matter atrophy quantification by advanced morphometric MRI can help to better characterize the course of neurodegenerative diseases in vivo, such as e.g. lower motor neuron disorders. Imaging the lower thoracic cord - containing those motor neurons that control leg function - could be particularly informative, however, is challenging due to tissue composition, physiological motion and large field of views.

NEW METHOD

An "averaged magnetization inversion recovery acquisitions" (AMIRA) approach with a radial k-space acquisition scheme was developed. The method is designed for morphometric SC imaging with a focus on the thoracic SC.

RESULTS

In a typical setting, radial AMIRA acquires transverse slices with a high 0.50 × 0.50mm² in-plane resolution and a pronounced positive contrast between thoracic gray and white matter, within typically 2:39 min. Additional proof-of-concept measurements in patients demonstrate that such contrast and resolving capability is indeed necessary to assess potential atrophy of the anterior horns.

COMPARISON WITH EXISTING METHOD(S)

Radial AMIRA utilizes two benefits of radial MRI techniques: being generally less prone to motion effects and that fold over artifacts can manifest less intrusively. These benefits are united with the original AMIRA approach which allows the contrast to be 'tuned' and improved based on the combination of five simultaneously acquired images of different tissue contrast.

CONCLUSIONS

Radial AMIRA is a promising approach for in vivo SC gray and white matter atrophy visualization and quantification in lower motor neuron diseases and other autoimmune or genetic diseases involving the entire (not only cervical) spinal cord.

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

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

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

What role should spinal cord MRI take in the future of multiple sclerosis surveillance?

INTRODUCTION

In multiple sclerosis (MS), inflammatory, demyelinating, and neurodegenerative phenomena affect the spinal cord, with detrimental effects on patients' clinical disability. Although spinal cord imaging may be challenging, improvements in MRI technologies have contributed to better evaluate spinal cord involvement in MS.

AREAS COVERED

This review summarizes the current state-of-art of the application of conventional and advanced MRI techniques to evaluate spinal cord damage in MS. Typical features of spinal cord lesions, their role in the diagnostic work-up of suspected MS, their predictive role for subsequent disease course and clinical worsening, and their utility to define treatment response are discussed. The role of spinal cord atrophy and of other advanced MRI techniques to better evaluate the associations between spinal cord abnormalities and the accumulation of clinical disability are also evaluated. Finally, how spinal cord assessment could evolve in the future to improve monitoring of disease progression and treatment effects is examined.

EXPERT OPINION

Spinal cord MRI provides relevant additional information to brain MRI in understanding MS pathophysiology, in allowing an earlier and more accurate diagnosis of MS, and in identifying MS patients at higher risk to develop more severe disability. A future role in monitoring the effects of treatments is also foreseen.

Generalised boundary shift integral for longitudinal assessment of spinal cord atrophy

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

Diffusion Kurtosis Imaging maps neural damage in the EAE model of multiple sclerosis

Diffusion kurtosis imaging (DKI) is an imaging modality that yields novel disease biomarkers and in combination with nervous tissue modeling, provides access to microstructural parameters. Recently, DKI and subsequent estimation of microstructural model parameters has been used for assessment of tissue changes in neurodegenerative diseases and associated animal models. In this study, mouse spinal cords from the experimental autoimmune encephalomyelitis (EAE) model of multiple sclerosis (MS) were investigated for the first time using DKI in combination with biophysical modeling to study the relationship between microstructural metrics and degree of animal dysfunction. Thirteen spinal cords were extracted from animals with varied grades of disability and scanned in a high-field MRI scanner along with five control specimen. Diffusion weighted data were acquired together with high resolution T₂* images. Diffusion data were fit to estimate diffusion and kurtosis tensors and white matter modeling parameters, which were all used for subsequent statistical analysis using a linear mixed effects model. T₂* images were used to delineate focal demyelination/inflammation. Our results reveal a strong relationship between disability and measured microstructural parameters in normal appearing white matter and gray matter. Relationships between disability and mean of the kurtosis tensor, radial kurtosis, radial diffusivity were similar to what has been found in other hypomyelinating MS models, and in patients. However, the changes in biophysical modeling parameters and in particular in extra-axonal axial diffusivity were clearly different from previous studies employing other animal models of MS. In conclusion, our data suggest that DKI and microstructural modeling can provide a unique contrast capable of detecting EAE-specific changes correlating with clinical disability.

Traumatic and nontraumatic spinal cord injury: pathological insights from neuroimaging

Pathophysiological changes in the spinal cord white and grey matter resulting from injury can be observed with MRI techniques. These techniques provide sensitive markers of macrostructural and microstructural tissue integrity, which correlate with histological findings. Spinal cord MRI findings in traumatic spinal cord injury (tSCI) and nontraumatic spinal cord injury - the most common form of which is degenerative cervical myelopathy (DCM) - have provided important insights into the pathophysiological processes taking place not just at the focal injury site but also rostral and caudal to the spinal injury. Although tSCI and DCM have different aetiologies, they show similar degrees of spinal cord pathology remote from the injury site, suggesting the involvement of similar secondary degenerative mechanisms. Advanced quantitative MRI protocols that are sensitive to spinal cord pathology have the potential to improve diagnosis and, more importantly, predict outcomes in patients with tSCI or nontraumatic spinal cord injury. This Review describes the insights into tSCI and DCM that have been revealed by neuroimaging and outlines current activities and future directions for the field.

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.

Automatic Spinal Cord Gray Matter Quantification: A Novel Approach

BACKGROUND AND PURPOSE

Currently, accurate and reproducible spinal cord GM segmentation remains challenging and a noninvasive broadly accepted reference standard for spinal cord GM measurements is still a matter of ongoing discussion. Our aim was to assess the reproducibility and accuracy of cervical spinal cord GM and WM cross-sectional area measurements using averaged magnetization inversion recovery acquisitions images and a fully-automatic postprocessing segmentation algorithm.

MATERIALS AND METHODS

The cervical spinal cord of 24 healthy subjects (14 women; mean age, 40 ± 11 years) was scanned in a test-retest fashion on a 3T MR imaging system. Twelve axial averaged magnetization inversion recovery acquisitions slices were acquired over a 48-mm cord segment. GM and WM were both manually segmented by 2 experienced readers and compared with an automatic variational segmentation algorithm with a shape prior modified for 3D data with a slice similarity prior. Precision and accuracy of the automatic method were evaluated using coefficients of variation and Dice similarity coefficients.

RESULTS

The mean GM area was 17.20 ± 2.28 mm² and the mean WM area was 72.71 ± 7.55 mm² using the automatic method. Reproducibility was high for both methods, while being better for the automatic approach (all mean automatic coefficients of variation, ≤4.77%; all differences, P < .001). The accuracy of the automatic method compared with the manual reference standard was excellent (mean Dice similarity coefficients: 0.86 ± 0.04 for GM and 0.90 ± 0.03 for WM). The automatic approach demonstrated similar coefficients of variation between intra- and intersession reproducibility as well as among all acquired spinal cord slices.

CONCLUSIONS

Our novel approach including the averaged magnetization inversion recovery acquisitions sequence and a fully-automated postprocessing segmentation algorithm demonstrated an accurate and reproducible spinal cord GM and WM segmentation. This pipeline is promising for both the exploration of longitudinal structural GM changes and application in clinical settings in disorders affecting the 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.

Spatiotemporal trajectories of quantitative magnetization transfer measurements in injured spinal cord using simplified acquisitions

Purpose

This study aims to systematically evaluate the accuracy and precision of pool size ratio (PSR) measurements from quantitative magnetization transfer (qMT) acquisitions using simplified models in the context of assessing injury-associated spatiotemporal changes in spinal cords of non-human primates. This study also aims to characterize changes in the spinal tissue pathology in individual subjects, both regionally and longitudinally, in order to demonstrate the relationship between regional tissue compositional changes and sensorimotor behavioral recovery after cervical spinal cord injury (SCI).

Methods

MRI scans were recorded on anesthetized monkeys at 9.4 T, before and serially after a unilateral section of the dorsal column tract. Images were acquired following saturating RF pulses at different offset frequencies. Models incorporating two pools of protons but with differing numbers of variable parameters were used to fit the data to derive qMT parameters. The results using different amounts of measured data and assuming different numbers of variable model parameters were compared. Behavioral impairments and recovery were assessed by a food grasping-retrieving task. Histological sections were obtained post mortem for validation of the injury.

Results

QMT fitting provided maps of pool size ratio (PSR), the relative amounts of immobilized protons exchanging magnetization compared to the “free” water. All the selected modeling approaches detected a lesion/cyst at the site of injury as significant reductions in PSR values. The regional contrasts in the PSR maps obtained using the different fittings varied, but the 2-parameter fitting results showed strong positive correlations with results from 5-parameter modeling. 2-parameter fitting results with modest (>3) RF offsets showed comparable sensitivity for detecting demyelination in white matter and loss of macromolecules in gray matter around lesion sites compared to 5-parameter fitting with fully-sampled data acquisitions. Histology confirmed that decreases of PSR corresponded to regional demyelination around lesion sites, especially when demyelination occurred along the dorsal column on the injury side. Longitudinally, PSR values of injured dorsal column tract and gray matter horns exhibited remarkable recovery that associated with behavioral improvement.

Conclusion

Simplified qMT modeling approaches provide efficient and sensitive means to detect and characterize injury-associated demyelination in white matter tracts and loss of macromolecules in gray matter and to monitor its recovery over time.

•

Simplified 2-parameter and fully sampled 5-parameter qMT modeling achieved comparable accuracy and precision of PSR values.

•

Successfully tracked and differentiated myelination states of specific WM tracts and macromolecular changes in GM horns.

•

Recovery of WM and GM pathology assessed by qMT correlated with improvements in hand uses after injury.

•

High translational potential for clinical studies of human patients with spinal cord injury.

Advances in spinal cord imaging in multiple sclerosis

The spinal cord is frequently affected in multiple sclerosis (MS), causing motor, sensory and autonomic dysfunction. A number of pathological abnormalities, including demyelination and neuroaxonal loss, occur in the MS spinal cord and are studied in vivo with magnetic resonance imaging (MRI). The aim of this review is to summarise and discuss recent advances in spinal cord MRI. Advances in conventional spinal cord MRI include improved identification of MS lesions, recommended spinal cord MRI protocols, enhanced recognition of MRI lesion characteristics that allow MS to be distinguished from other myelopathies, evidence for the role of spinal cord lesions in predicting prognosis and monitoring disease course, and novel post-processing methods to obtain lesion probability maps. The rate of spinal cord atrophy is greater than that of brain atrophy (-1.78% versus -0.5% per year), and reflects neuroaxonal loss in an eloquent site of the central nervous system, suggesting that it can become an important outcome measure in clinical trials, especially in progressive MS. Recent developments allow the calculation of spinal cord atrophy from brain volumetric scans and evaluation of its progression over time with registration-based techniques. Fully automated analysis methods, including segmentation of grey matter and intramedullary lesions, will facilitate the use of spinal cord atrophy in trial designs and observational studies. Advances in quantitative imaging techniques to evaluate neuroaxonal integrity, myelin content, metabolic changes, and functional connectivity, have provided new insights into the mechanisms of damage in MS. Future directions of research and the possible impact of 7T scanners on spinal cord imaging will be discussed.

Intrinsic functional architecture of the non-human primate spinal cord derived from fMRI and electrophysiology

Resting-state functional MRI (rsfMRI) has recently revealed correlated signals in the spinal cord horns of monkeys and humans. However, the interpretation of these rsfMRI correlations as indicators of functional connectivity in the spinal cord remains unclear. Here, we recorded stimulus-evoked and spontaneous spiking activity and local field potentials (LFPs) from monkey spinal cord in order to validate fMRI measures. We found that both BOLD and electrophysiological signals elicited by tactile stimulation co-localized to the ipsilateral dorsal horn. Temporal profiles of stimulus-evoked BOLD signals covaried with LFP and multiunit spiking in a similar way to those observed in the brain. Functional connectivity of dorsal horns exhibited a U-shaped profile along the dorsal-intermediate-ventral axis. Overall, these results suggest that there is an intrinsic functional architecture within the gray matter of a single spinal segment, and that rsfMRI signals at high field directly reflect this underlying spontaneous neuronal activity.

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.

Isoflurane anesthesia does not affect spinal cord neurovascular coupling: evidence from decerebrated rats

Neurological examination remains the primary clinical investigation in patients with spinal cord injury. However, neuroimaging methods such as functional magnetic resonance imaging (fMRI) are promising tools for following functional changes in the course of injury, disease and rehabilitation. However, the relationship between neuronal activity and blood flow in the spinal cord on which fMRI relies has been largely overlooked. The objective of this study was to examine neurovascular coupling in the spinal cord of decerebrated rats during electrical stimulation of the sciatic nerve with and without isoflurane anesthesia (1.2%). Local field potentials (LFP) and spinal cord blood flow (SCBF) were recorded simultaneously in the lumbosacral enlargement. Isoflurane did not significantly alter LFP (p = 0.53) and SCBF (p = 0.57) amplitude. Accordingly, neurovascular coupling remained comparable with or without isoflurane anesthesia (p = 0.39). These results support the use of isoflurane in rodents to investigate nociceptive functions of the spinal cord using fMRI.

Dynamic Computational Model of the Human Spinal Cord Connectome

Connectomes abound, but few for the human spinal cord. Using anatomical data in the literature, we constructed a draft connectivity map of the human spinal cord connectome, providing a template for the many calibrations of specialized behavior to be overlaid on it and the basis for an initial computational model. A thorough literature review gleaned cell types, connectivity, and connection strength indications. Where human data were not available, we selected species that have been studied. Cadaveric spinal cord measurements, cross-sectional histology images, and cytoarchitectural data regarding cell size and density served as the starting point for estimating numbers of neurons. Simulations were run using neural circuitry simulation software. The model contains the neural circuitry in all ten Rexed laminae with intralaminar, interlaminar, and intersegmental connections, as well as ascending and descending brain connections and estimated neuron counts for various cell types in every lamina of all 31 segments. We noted the presence of highly interconnected complex networks exhibiting several orders of recurrence. The model was used to perform a detailed study of spinal cord stimulation for analgesia. This model is a starting point for workers to develop and test hypotheses across an array of biomedical applications focused on the spinal cord. Each such model requires additional calibrations to constrain its output to verifiable predictions. Future work will include simulating additional segments and expanding the research uses of the model.

Gray vs. White Matter Segmentation of the Conus Medullaris: Reliability and Variability in Healthy Volunteers

BACKGROUND AND PURPOSE

Magnetic resonance imaging (MRI)-derived spinal cord (SC) gray and white matter (GM/WM) volume are useful indirect measures of atrophy and neurodegeneration over time, typically obtained in the upper SC. Neuropathological evidence suggests that in certain neurological conditions, early degeneration may occur as low as the sacral SC. In this study, the feasibility of GM/WM segmentation of the conus medullaris (CM) was assessed in vivo.

METHODS

Twenty-three healthy volunteers (11 female, mean age 47 years) underwent high-resolution 3T MRI of the CM using a 3-dimensional fast field echo sequence. Reproducibility of the volume measurements was assessed in 5 subjects (2 female, 25-37 years) by one rater who repeated the analysis 3 times and also with 2 additional raters working independently in order to calculate the intra- and interrater coefficient of variation (COV), respectively. Furthermore, the influence of age, gender, spine and SC metrics on tissue-specific measures of the CM was investigated.

RESULTS

Volumetric CM analyses (N = 23) for the SC, GM, and WM revealed a mean (SD) total volume of CM-TV = 1746.9 (296.7) mm³ , CM-GM-TV = 731.2 (106.0) mm³ , and CM-WM-TV = 1014.6 (211.3) mm³ , respectively. The intra-rater COV for measuring the CM-TV and CM-GM-TV was 3.38% and 7.42%, respectively; the interrater COV was 3.43% and 10.80%, respectively. Using age, gender, spine and SC metrics in regression models substantially reduced group variability for CM-TV, CM-WM-TV, and CM-GM-TV by up to 39.2%, 42.7%, and 21.2%, respectively.

CONCLUSIONS

The results from this study demonstrate the feasibility of obtaining tissue-specific volume measurements in the CM by means of MRI with good reproducibility and provide normative data for future applications in neurological diseases affecting the lower SC.

Monitoring Progressive Multiple Sclerosis with Novel Imaging Techniques

Imaging markers for monitoring disease progression in progressive multiple sclerosis (PMS) are scarce, thereby limiting the possibility to monitor disease evolution and to test effective treatments in clinical trials. Advanced imaging techniques that have the advantage of metrics with increased sensitivity to short-term tissue changes and increased specificity to the structural abnormalities characteristic of PMS have recently been applied in clinical trials of PMS. In this review, we (1) provide an overview of the pathological features of PMS, (2) summarize the findings of research and clinical trials conducted in PMS which have applied conventional and advanced magnetic resonance imaging techniques and (3) discuss recent advancements and future perspectives in monitoring PMS with imaging techniques.